Abstract:SUMMARY Performance improvements in biomolecular simulations based on molecular dynamics (MD) codes are widely desired. Unfortunately, the factors, which allowed past performance improvements, particularly the microprocessor clock frequencies, are no longer increasing. Hence, novel software and hardware solutions are being explored for accelerating performance of widely used MD codes. In this paper, we describe our efforts on porting, optimizing and tuning of Large‐scale Atomic/Molecular Massively Parallel Simulator, a popular MD framework, on heterogeneous architectures: multi‐core processors with graphical processing unit (GPU) accelerators. Our implementation is based on accelerating the most computationally expensive non‐bonded interaction terms on the GPUs and overlapping the computation on the CPU and GPUs. This functionality is built on top of message passing interface that allows multi‐level parallelism to be extracted even at the workstation level with the multi‐core CPUs and allows extension of the implementation on GPU‐enabled clusters. We hypothesize that the optimal benefit of heterogeneous architectures for applications will come by utilizing all possible resources (for example, CPU‐cores and GPU devices on GPU‐enabled clusters). Benchmarks for a range of biomolecular system sizes are provided, and an analysis is performed on four generations of NVIDIA's GPU devices. On GPU‐enabled Linux clusters, by overlapping and pipelining computation and communication, we observe up to 10‐folds application acceleration in multi‐core and multi‐GPU environments illustrating significant performance improvements. Detailed analysis of the implementation is presented that allows identification of bottlenecks in algorithm, indicating that code optimization and improvements on GPUs could allow microsecond scale simulation throughput on workstations and inexpensive GPU clusters, putting widely desired biologically relevant simulation time‐scales within reach of a large user community. In order to systematically optimize simulation throughput and to enable performance prediction, we have developed a parameterized performance model that will allow developers and users to explore the performance potential of future heterogeneous systems for biological simulations. Copyright © 2012 John Wiley & Sons, Ltd.

Process Simulation of Complex Biological Pathways in Physical Reactive Space and Reformulated for Massively Parallel Computing Platforms

Computational Framework for In-Silico Study of Virtual Cell Biology Via Process Simulation and Multiscale Modeling

A Systematic Petri Net Approach for Multiple-Scale Modeling and Simulation of Biochemical Processes

Smoldyn on Graphics Processing Units: Massively Parallel Brownian Dynamics Simulations

A Scalable Moment-Closure Approximation for Large-Scale Biochemical Reaction Networks

Fine-Grained Parallel And Distributed Spatial Stochastic Simulation Of Biological Reactions

Exact Parallelization of the Stochastic Simulation Algorithm for Scalable Simulation of Large Biochemical Networks

A Study of Parallel and Evolutionary Framework for Modelling Biochemical Signalling Pathways

Parallelization and High-Performance Computing Enables Automated Statistical Inference of Multi-scale Models.

3D mesh processing using GAMer 2 to enable reaction-diffusion simulations in realistic cellular geometries

Performance modeling of microsecond scale biological molecular dynamics simulations on heterogeneous architectures

Green's Function Reaction Dynamics: a new approach to simulate biochemical networks at the particle level and in time and space

Spatial-Stochastic Simulation of Reaction-Diffusion Systems

Stochastic Simulation of Biomolecular Networks in Dynamic Environments

Multi-scale Approach for Simulating Time-Delay Biochemical Reaction Systems.

pSpatiocyte: A Parallel Stochastic Method for Particle Reaction-Diffusion Systems

Novel Ground-Up 3D Multicellular Simulators for Synthetic Biology CAD Integrating Stochastic Gillespie Simulations Benchmarked with Topologically Variable SBML Models

Generalizing Gillespie's direct method to enable network-free simulations

MCell4 with BioNetGen: A Monte Carlo simulator of rule-based reaction-diffusion systems with Python interface

A Hybrid Parallel Framework for the Cellular Potts Model Simulations

Simulation and computational analysis of multiscale graph agent-based tumor model