Abstract:Recent advances in computational neuroscience have demonstrated the usefulness and importance of stochastic, spatial reaction-diffusion simulations. However, ever increasing model complexity renders traditional serial solvers, as well as naive parallel implementations, inadequate. This paper introduces a new generation of the STochastic Engine for Pathway Simulation (STEPS) project (http://steps.sourceforge.net/), denominated STEPS 4.0, and its core components which have been designed for improved scalability, performance, and memory efficiency. STEPS 4.0 aims to enable novel scientific studies of macroscopic systems such as whole cells while capturing their nanoscale details. This class of models is out of reach for serial solvers due to the vast quantity of computation in such detailed models, and also out of reach for naive parallel solvers due to the large memory footprint. Based on a distributed mesh solution, we introduce a new parallel stochastic reaction-diffusion solver and a deterministic membrane potential solver in STEPS 4.0. The distributed mesh, together with improved data layout and algorithm designs, significantly reduces the memory footprint of parallel simulations in STEPS 4.0. This enables massively parallel simulations on modern HPC clusters and overcomes the limitations of the previous parallel STEPS implementation. Current and future improvements to the solver are not sustainable without following proper software engineering principles. For this reason, we also give an overview of how the STEPS codebase and the development environment have been updated to follow modern software development practices. We benchmark performance improvement and memory footprint on three published models with different complexities, from a simple spatial stochastic reaction-diffusion model, to a more complex one that is coupled to a deterministic membrane potential solver to simulate the calcium burst activity of a Purkinje neuron. Simulation results of these models suggest that the new solution dramatically reduces the per-core memory consumption by more than a factor of 30, while maintaining similar or better performance and scalability.

Steps and bumps: precision extraction of discrete states of molecular machines using physically-based, high-throughput time series analysis

An accurate probabilistic step finder for time-series analysis

Waiting time distributions in hybrid models of motor-bead assays: A concept and tool for inference

Renewal processes and fluctuation analysis of molecular motor stepping

Automated Event Detection and Activity Monitoring in Long Molecular Dynamics Simulations

Single-molecule Automata: Harnessing Kinetic-Thermodynamic Discrepancy for Temporal Pattern Recognition

Dwell Time Symmetry in Random Walks and Molecular Motors

Stepping and crowding of molecular motors: statistical kinetics from an exclusion process perspective

Efficient use of single molecule time traces to resolve kinetic rates, models and uncertainties

Detecting substeps in the rotary motors of FoF1-ATP synthase by Hidden Markov Models

Functional Protein Dynamics in a Crystal

Biomolecular Simulation: A Computational Microscope for Molecular Biology

High-speed Parameter Search of Dynamic Biological Pathways from Time-Course Transcriptomic Profiles Using High-Level Petri Net.

STEPS 4.0: Fast and memory-efficient molecular simulations of neurons at the nanoscale

Diagnosing Heterogeneous Dynamics in Single Molecule/Particle Trajectories with Multiscale Wavelets

The pathway elaboration method for mean first passage time estimation in large continuous-time Markov chains with applications to nucleic acid kinetics

Molecular free energies, rates, and mechanisms from data-efficient path sampling simulations

Markov field models: Scaling molecular kinetics approaches to large molecular machines

Dynamic coarse-graining fills the gap between atomistic simulations and experimental investigations of mechanical unfolding

Molecular Dynamics Simulation for All

In-Silico Analyses of Molecular Force Sensors for Mechanical Characterization of Biological Systems