myPresto/omegagene: a GPU-accelerated molecular dynamics simulator tailored for enhanced conformational sampling methods with a non-Ewald electrostatic scheme

Molecular dynamics (MD) is a promising computational approach to investigate dynamical behavior of molecular systems at the atomic level. Here, we present a new MD simulation engine named “myPresto/omegagene” that is tailored for enhanced conformational sampling methods with a non-Ewald electrostatic potential scheme. Our enhanced conformational sampling methods, e.g., the virtual-system-coupled multi-canonical MD (V-McMD) method, replace a multi-process parallelized run with multiple independent runs to avoid inter-node communication overhead. In addition, adopting the non-Ewald-based zero-multipole summation method (ZMM) makes it possible to eliminate the Fourier space calculations altogether. The combination of these state-of-the-art techniques realizes efficient and accurate calculations of the conformational ensemble at an equilibrium state. By taking these advantages, myPresto/omegagene is specialized for the single process execution with Graphics Processing Unit (GPU). We performed benchmark simulations for the 20-mer peptide, Trp-cage, with explicit solvent. One of the most thermodynamically stable conformations generated by the V-McMD simulation is very similar to an experimentally solved native conformation. Furthermore, the computation speed is four-times faster than that of our previous simulation engine, myPresto/psygene-G. The new simulator, myPresto/omegagene, is freely available at the following URLs: http://www.protein.osaka-u.ac.jp/rcsfp/pi/omegagene/ and http://presto.protein.osaka-u.ac.jp/myPresto4/.

[1]  Haruki Nakamura,et al.  Theory for trivial trajectory parallelization of multicanonical molecular dynamics and application to a polypeptide in water , 2011, J. Comput. Chem..

[2]  Haruki Nakamura,et al.  The zero-multipole summation method for estimating electrostatic interactions in molecular dynamics: analysis of the accuracy and application to liquid systems. , 2014, The Journal of chemical physics.

[3]  Haruki Nakamura,et al.  Improved Estimation of Protein-Ligand Binding Free Energy by Using the Ligand-Entropy and Mobility of Water Molecules , 2013, Pharmaceuticals.

[4]  Haruki Nakamura,et al.  A virtual-system coupled multicanonical molecular dynamics simulation: principles and applications to free-energy landscape of protein-protein interaction with an all-atom model in explicit solvent. , 2013, The Journal of chemical physics.

[5]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born , 2012, Journal of chemical theory and computation.

[6]  Ian J. Grimstead,et al.  GPU‐accelerated molecular mechanics computations , 2013, J. Comput. Chem..

[7]  Yuji Sugita,et al.  Surface-Tension Replica-Exchange Molecular Dynamics Method for Enhanced Sampling of Biological Membrane Systems. , 2013, Journal of chemical theory and computation.

[8]  Haruki Nakamura,et al.  Integration of Ligand-Based Drug Screening with Structure-Based Drug Screening by Combining Maximum Volume Overlapping Score with Ligand Doscking , 2012, Pharmaceuticals.

[9]  Haruki Nakamura,et al.  A Novel Approach of Dynamic Cross Correlation Analysis on Molecular Dynamics Simulations and Its Application to Ets1 Dimer–DNA Complex , 2014, PloS one.

[10]  A. Kidera,et al.  Multicanonical Ensemble Generated by Molecular Dynamics Simulation for Enhanced Conformational Sampling of Peptides , 1997 .

[11]  Ron O. Dror,et al.  Mechanism of Voltage Gating in Potassium Channels , 2012, Science.

[12]  Haruki Nakamura,et al.  Enhanced conformational sampling to visualize a free-energy landscape of protein complex formation , 2016, The Biochemical journal.

[13]  Haruki Nakamura,et al.  Molecular Dynamics Simulations Accelerated by GPU for Biological Macromolecules with a Non-Ewald Scheme for Electrostatic Interactions. , 2013, Journal of chemical theory and computation.

[14]  J. W. Neidigh,et al.  Designing a 20-residue protein , 2002, Nature Structural Biology.

[15]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[16]  Haruki Nakamura,et al.  A free-energy landscape for coupled folding and binding of an intrinsically disordered protein in explicit solvent from detailed all-atom computations. , 2011, Journal of the American Chemical Society.

[17]  Haruki Nakamura,et al.  Molecular dynamics scheme for precise estimation of electrostatic interaction via zero-dipole summation principle. , 2011, The Journal of chemical physics.

[18]  Ikuo Fukuda,et al.  Zero-multipole summation method for efficiently estimating electrostatic interactions in molecular system. , 2013, The Journal of chemical physics.

[19]  Yuji Sugita,et al.  Efficient lookup table using a linear function of inverse distance squared , 2013, J. Comput. Chem..

[20]  Makoto Taiji,et al.  MDGRAPE-4: a special-purpose computer system for molecular dynamics simulations , 2014, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[21]  Haruki Nakamura,et al.  Molecular Dynamics Simulations of Double-Stranded DNA in an Explicit Solvent Model with the Zero-Dipole Summation Method , 2013, PloS one.

[22]  Y. Fukunishi,et al.  Similarities among receptor pockets and among compounds: analysis and application to in silico ligand screening. , 2005, Journal of molecular graphics & modelling.

[23]  Mitsunori Ikeguchi,et al.  Partial rigid‐body dynamics in NPT, NPAT and NPγT ensembles for proteins and membranes , 2004, J. Comput. Chem..

[24]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[25]  Haruki Nakamura,et al.  Enhanced and effective conformational sampling of protein molecular systems for their free energy landscapes , 2012, Biophysical Reviews.

[26]  J. P. Grossman,et al.  Anton 2: Raising the Bar for Performance and Programmability in a Special-Purpose Molecular Dynamics Supercomputer , 2014, SC14: International Conference for High Performance Computing, Networking, Storage and Analysis.

[27]  David E. Shaw,et al.  166 Millisecond-long molecular dynamics simulations of proteins on a special-purpose machine , 2013 .

[28]  Takeshi Kawabata,et al.  LigandBox: A database for 3D structures of chemical compounds , 2013, Biophysics.

[29]  John L. Klepeis,et al.  Anton, a special-purpose machine for molecular dynamics simulation , 2007, ISCA '07.

[30]  Haruki Nakamura,et al.  Simple and accurate scheme to compute electrostatic interaction: zero-dipole summation technique for molecular system and application to bulk water. , 2012, The Journal of chemical physics.

[31]  D. Shaw.  Millisecond-Long Molecular Dynamics Simulations of Proteins on a Special-Purpose Machine , 2013 .

[32]  Haruki Nakamura,et al.  Prediction of ligand‐binding sites of proteins by molecular docking calculation for a random ligand library , 2011, Protein science : a publication of the Protein Society.

[33]  Haruki Nakamura,et al.  Non-Ewald methods: theory and applications to molecular systems , 2012, Biophysical Reviews.

[34]  Berk Hess,et al.  A flexible algorithm for calculating pair interactions on SIMD architectures , 2013, Comput. Phys. Commun..

[35]  Genji Kurisu,et al.  Structure of the entire stalk region of the Dynein motor domain. , 2014, Journal of molecular biology.

[36]  Berk Hess,et al.  GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers , 2015 .

[37]  Haruki Nakamura,et al.  Virtual‐system‐coupled adaptive umbrella sampling to compute free‐energy landscape for flexible molecular docking , 2015, J. Comput. Chem..

[38]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. , 2013, Journal of chemical theory and computation.

[39]  Masao Fukushima,et al.  MODYLAS: A Highly Parallelized General-Purpose Molecular Dynamics Simulation Program for Large-Scale Systems with Long-Range Forces Calculated by Fast Multipole Method (FMM) and Highly Scalable Fine-Grained New Parallel Processing Algorithms. , 2013, Journal of chemical theory and computation.

[40]  Yuji Sugita,et al.  Midpoint cell method for hybrid (MPI+OpenMP) parallelization of molecular dynamics simulations , 2014, J. Comput. Chem..

[41]  R. Dror,et al.  Improved side-chain torsion potentials for the Amber ff99SB protein force field , 2010, Proteins.

[42]  Atsushi Matsumoto,et al.  Free-Energy Landscape of Reverse tRNA Translocation through the Ribosome Analyzed by Electron Microscopy Density Maps and Molecular Dynamics Simulations , 2014, PloS one.

[43]  T. Cheatham,et al.  Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations , 2008, The journal of physical chemistry. B.