Biomolecular simulation: a computational microscope for molecular biology.

Molecular dynamics simulations capture the behavior of biological macromolecules in full atomic detail, but their computational demands, combined with the challenge of appropriately modeling the relevant physics, have historically restricted their length and accuracy. Dramatic recent improvements in achievable simulation speed and the underlying physical models have enabled atomic-level simulations on timescales as long as milliseconds that capture key biochemical processes such as protein folding, drug binding, membrane transport, and the conformational changes critical to protein function. Such simulation may serve as a computational microscope, revealing biomolecular mechanisms at spatial and temporal scales that are difficult to observe experimentally. We describe the rapidly evolving state of the art for atomic-level biomolecular simulation, illustrate the types of biological discoveries that can now be made through simulation, and discuss challenges motivating continued innovation in this field.

[1]  J. Kirkwood Statistical Mechanics of Fluid Mixtures , 1935 .

[2]  R. Zwanzig High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[3]  M. Karplus,et al.  Dynamics of folded proteins , 1977, Nature.

[4]  William L. Jorgensen,et al.  Efficient computation of absolute free energies of binding by computer simulations. Application to the methane dimer in water , 1988 .

[5]  F Bezanilla,et al.  Solute inaccessible aqueous volume changes during opening of the potassium channel of the squid giant axon. , 1990, Biophysical journal.

[6]  R Fine,et al.  FASTRUN: A special purpose, hardwired computer for molecular simulation , 1991, Proteins.

[7]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[8]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[9]  A. Fersht,et al.  Synergy between simulation and experiment in describing the energy landscape of protein folding. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[10]  Nobuaki Miyakawa,et al.  Development of MD Engine: High-speed accelerator with parallel processor design for molecular dynamics simulations , 1999, J. Comput. Chem..

[11]  A. V. Duin,et al.  ReaxFF: A Reactive Force Field for Hydrocarbons , 2001 .

[12]  A. Roitberg,et al.  All-atom structure prediction and folding simulations of a stable protein. , 2002, Journal of the American Chemical Society.

[13]  K. Schulten,et al.  Control of the Selectivity of the Aquaporin Water Channel Family by Global Orientational Tuning , 2002, Science.

[14]  T. Narumi,et al.  Protein Explorer: A Petaflops Special-Purpose Computer System for Molecular Dynamics Simulations , 2003, ACM/IEEE SC 2003 Conference (SC'03).

[15]  Marc Snir,et al.  A Note on N-Body Computations with Cutoffs , 2004, Theory of Computing Systems.

[16]  B. Roux,et al.  Control of ion selectivity in potassium channels by electrostatic and dynamic properties of carbonyl ligands , 2004, Nature.

[17]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[18]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[19]  J Andrew McCammon,et al.  Discovery of a novel binding trench in HIV integrase. , 2004, Journal of medicinal chemistry.

[20]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[21]  Martin Gruebele,et al.  Engineering a β-sheet protein toward the folding speed limit , 2005 .

[22]  S. Joseph,et al.  Simulating movement of tRNA into the ribosome during decoding. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[23]  Eric J. Sorin,et al.  How well can simulation predict protein folding kinetics and thermodynamics? , 2005, Annual review of biophysics and biomolecular structure.

[24]  B. Roux,et al.  Calculation of absolute protein-ligand binding free energy from computer simulations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[25]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[26]  M. DePristo,et al.  Simultaneous determination of protein structure and dynamics , 2005, Nature.

[27]  R. Dror,et al.  Gaussian split Ewald: A fast Ewald mesh method for molecular simulation. , 2005, The Journal of chemical physics.

[28]  John Mongan,et al.  Biomolecular simulations at constant pH. , 2005, Current opinion in structural biology.

[29]  David E. Shaw,et al.  A fast, scalable method for the parallel evaluation of distance‐limited pairwise particle interactions , 2005, J. Comput. Chem..

[30]  M. Seibert,et al.  Reproducible polypeptide folding and structure prediction using molecular dynamics simulations. , 2005, Journal of molecular biology.

[31]  Richard A. Friesner,et al.  Integrated Modeling Program, Applied Chemical Theory (IMPACT) , 2005, J. Comput. Chem..

[32]  Alan Grossfield,et al.  Retinal counterion switch mechanism in vision evaluated by molecular simulations. , 2006, Journal of the American Chemical Society.

[33]  Stewart A. Adcock,et al.  Molecular dynamics: survey of methods for simulating the activity of proteins. , 2006, Chemical reviews.

[34]  Peter L. Freddolino,et al.  Molecular dynamics simulations of the complete satellite tobacco mosaic virus. , 2006, Structure.

[35]  Federico D. Sacerdoti,et al.  Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[36]  Robert S. Germain,et al.  Blue Matter: Strong Scaling of Molecular Dynamics on Blue Gene/L , 2006, International Conference on Computational Science.

[37]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[38]  Stewart A. Adcock,et al.  Molecular Dynamics: Survey of Methods for Simulating the Activity of Proteins , 2006 .

[39]  Brent A. Gregersen,et al.  Mechanism of Na+/H+ Antiporting , 2007, Science.

[40]  David E. Shaw,et al.  Zonal methods for the parallel execution of range-limited N-body simulations , 2007, J. Comput. Phys..

[41]  H. Stern Molecular simulation with variable protonation states at constant pH. , 2007, The Journal of chemical physics.

[42]  V. Pande,et al.  Heterogeneity even at the speed limit of folding: large-scale molecular dynamics study of a fast-folding variant of the villin headpiece. , 2007, Journal of molecular biology.

[43]  Benoît Roux,et al.  On the importance of a funneled energy landscape for the assembly and regulation of multidomain Src tyrosine kinases , 2007, Proceedings of the National Academy of Sciences.

[44]  David L Bostick,et al.  Selectivity in K+ channels is due to topological control of the permeant ion's coordinated state , 2007, Proceedings of the National Academy of Sciences.

[45]  J. Šponer,et al.  Refinement of the AMBER Force Field for Nucleic Acids: Improving the Description of α/γ Conformers , 2007 .

[46]  Nils G Walter,et al.  Molecular dynamics simulations of RNA: an in silico single molecule approach. , 2007, Biopolymers.

[47]  Y. Duan,et al.  Folding free-energy landscape of villin headpiece subdomain from molecular dynamics simulations , 2007, Proceedings of the National Academy of Sciences.

[48]  R. Dror,et al.  Dynamic control of slow water transport by aquaporin 0: Implications for hydration and junction stability in the eye lens , 2008, Proceedings of the National Academy of Sciences.

[49]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[50]  R. Larson,et al.  The MARTINI Coarse-Grained Force Field: Extension to Proteins. , 2008, Journal of chemical theory and computation.

[51]  William L. Jorgensen,et al.  Optimization of azoles as anti-human immunodeficiency virus agents guided by free-energy calculations. , 2008, Journal of the American Chemical Society.

[52]  Charles L Brooks,et al.  Recent advances in implicit solvent-based methods for biomolecular simulations. , 2008, Current opinion in structural biology.

[53]  Arun K Shukla,et al.  A crystal clear view of the β2-adrenergic receptor , 2008, Nature Biotechnology.

[54]  Edmond Chow,et al.  Desmond Performance on a Cluster of Multicore Processors , 2008 .

[55]  Peter L. Freddolino,et al.  Ten-microsecond molecular dynamics simulation of a fast-folding WW domain. , 2008, Biophysical journal.

[56]  Klaus Schulten,et al.  Adapting a message-driven parallel application to GPU-accelerated clusters , 2008, 2008 SC - International Conference for High Performance Computing, Networking, Storage and Analysis.

[57]  Joshua A. Anderson,et al.  General purpose molecular dynamics simulations fully implemented on graphics processing units , 2008, J. Comput. Phys..

[58]  R. Dror,et al.  A conserved protonation-dependent switch controls drug binding in the Abl kinase , 2009, Proceedings of the National Academy of Sciences.

[59]  Vijay S. Pande,et al.  Folding@home: Lessons from eight years of volunteer distributed computing , 2009, 2009 IEEE International Symposium on Parallel & Distributed Processing.

[60]  Fang Zheng,et al.  Free-energy perturbation simulation on transition states and redesign of butyrylcholinesterase. , 2009, Biophysical journal.

[61]  Ron O Dror,et al.  Identification of two distinct inactive conformations of the β2-adrenergic receptor reconciles structural and biochemical observations , 2009, Proceedings of the National Academy of Sciences.

[62]  Peter L. Freddolino,et al.  Common structural transitions in explicit-solvent simulations of villin headpiece folding. , 2009, Biophysical journal.

[63]  A. Cavalli,et al.  Protein conformational transitions: the closure mechanism of a kinase explored by atomistic simulations. , 2009, Journal of the American Chemical Society.

[64]  John Kuriyan,et al.  Equally potent inhibition of c-Src and Abl by compounds that recognize inactive kinase conformations. , 2009, Cancer research.

[65]  Benoît Roux,et al.  Mapping the conformational transition in Src activation by cumulating the information from multiple molecular dynamics trajectories , 2009, Proceedings of the National Academy of Sciences.

[66]  J. P. Grossman,et al.  Millisecond-scale molecular dynamics simulations on Anton , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[67]  M J Harvey,et al.  ACEMD: Accelerating Biomolecular Dynamics in the Microsecond Time Scale. , 2009, Journal of chemical theory and computation.

[68]  Stefano Vanni,et al.  Observation of "ionic lock" formation in molecular dynamics simulations of wild-type beta 1 and beta 2 adrenergic receptors. , 2009, Biochemistry.

[69]  G. Hummer,et al.  Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides. , 2009, The journal of physical chemistry. B.

[70]  G. Voth,et al.  A role for a specific cholesterol interaction in stabilizing the Apo configuration of the human A(2A) adenosine receptor. , 2009, Structure.

[71]  Vijay S. Pande,et al.  Accelerating molecular dynamic simulation on graphics processing units , 2009, J. Comput. Chem..

[72]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[73]  Kresten Lindorff-Larsen,et al.  Principles of conduction and hydrophobic gating in K+ channels , 2010, Proceedings of the National Academy of Sciences.

[74]  E. Tajkhorshid,et al.  Simulation of spontaneous substrate binding revealing the binding pathway and mechanism and initial conformational response of GlpT. , 2010, Biochemistry.

[75]  Alexander D. MacKerell,et al.  Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types. , 2010, The journal of physical chemistry. B.

[76]  Margaret E. Johnson,et al.  Current status of the AMOEBA polarizable force field. , 2010, The journal of physical chemistry. B.

[77]  Kyle A. Beauchamp,et al.  Molecular simulation of ab initio protein folding for a millisecond folder NTL9(1-39). , 2010, Journal of the American Chemical Society.

[78]  Tod D Romo,et al.  A Lipid Pathway for Ligand Binding Is Necessary for a Cannabinoid G Protein-coupled Receptor* , 2010, The Journal of Biological Chemistry.

[79]  Ron O. Dror,et al.  Exploring atomic resolution physiology on a femtosecond to millisecond timescale using molecular dynamics simulations , 2010, The Journal of general physiology.

[80]  Tod D Romo,et al.  Concerted interconversion between ionic lock substates of the beta(2) adrenergic receptor revealed by microsecond timescale molecular dynamics. , 2010, Biophysical journal.

[81]  C. Baragli,et al.  Selecting sequences that fold into a defined 3D structure: A new approach for protein design based on molecular dynamics and energetics. , 2010, Biophysical chemistry.

[82]  Edmond Chow,et al.  Exploiting 162-Nanosecond End-to-End Communication Latency on Anton , 2010, 2010 ACM/IEEE International Conference for High Performance Computing, Networking, Storage and Analysis.

[83]  David Baker,et al.  Evaluation and ranking of enzyme designs , 2010, Protein science : a publication of the Protein Society.

[84]  Alexander D. MacKerell,et al.  CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields , 2009, J. Comput. Chem..

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

[86]  Joseph A. Bank,et al.  Supporting Online Material Materials and Methods Figs. S1 to S10 Table S1 References Movies S1 to S3 Atomic-level Characterization of the Structural Dynamics of Proteins , 2022 .

[87]  R Dustin Schaeffer,et al.  Dynameomics: a comprehensive database of protein dynamics. , 2010, Structure.

[88]  D. Turner,et al.  Benchmarking AMBER Force Fields for RNA: Comparisons to NMR Spectra for Single-Stranded r(GACC) Are Improved by Revised χ Torsions , 2011, The journal of physical chemistry. B.

[89]  Eric T. Kim,et al.  How does a drug molecule find its target binding site? , 2011, Journal of the American Chemical Society.

[90]  Albert C. Pan,et al.  Pathway and mechanism of drug binding to G-protein-coupled receptors , 2011, Proceedings of the National Academy of Sciences.

[91]  Nagarajan Vaidehi,et al.  The role of conformational ensembles in ligand recognition in G-protein coupled receptors. , 2011, Journal of the American Chemical Society.

[92]  Gebhard F. X. Schertler,et al.  Two distinct conformations of helix 6 observed in antagonist-bound structures of a β1-adrenergic receptor , 2011, Proceedings of the National Academy of Sciences.

[93]  Cheng Zhang,et al.  Structure and Function of an Irreversible Agonist-β2 Adrenoceptor complex , 2010, Nature.

[94]  Tobin R Sosnick,et al.  The folding of single domain proteins--have we reached a consensus? , 2011, Current opinion in structural biology.

[95]  Klaus Schulten,et al.  High-performance scalable molecular dynamics simulations of a polarizable force field based on classical Drude oscillators in NAMD. , 2011, The journal of physical chemistry letters.

[96]  Paul S. Crozier,et al.  General-purpose molecular dynamics simulations on GPU-based clusters , 2011 .

[97]  G. de Fabritiis,et al.  Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations , 2011, Proceedings of the National Academy of Sciences.

[98]  S. Rasmussen,et al.  Structure of a nanobody-stabilized active state of the β2 adrenoceptor , 2010, Nature.

[99]  Daniel M Zuckerman,et al.  Equilibrium sampling in biomolecular simulations. , 2011, Annual review of biophysics.

[100]  R. Dror,et al.  How Fast-Folding Proteins Fold , 2011, Science.

[101]  Arieh Warshel,et al.  The empirical valence bond model: theory and applications , 2011 .

[102]  Krzysztof Fidelis,et al.  CASP9 results compared to those of previous casp experiments , 2011, Proteins.

[103]  K. Lindorff-Larsen,et al.  How robust are protein folding simulations with respect to force field parameterization? , 2011, Biophysical journal.

[104]  Albert C. Pan,et al.  Activation mechanism of the β2-adrenergic receptor , 2011, Proceedings of the National Academy of Sciences.

[105]  J. Šponer,et al.  Refinement of the Cornell et al. Nucleic Acids Force Field Based on Reference Quantum Chemical Calculations of Glycosidic Torsion Profiles , 2011, Journal of chemical theory and computation.

[106]  William L. Jorgensen,et al.  Efficient discovery of potent anti-HIV agents targeting the Tyr181Cys variant of HIV reverse transcriptase. , 2011, Journal of the American Chemical Society.

[107]  Alexander D. MacKerell,et al.  Impact of 2′‐hydroxyl sampling on the conformational properties of RNA: Update of the CHARMM all‐atom additive force field for RNA , 2011, J. Comput. Chem..

[108]  M. Gruebele,et al.  Computational design and experimental testing of the fastest-folding β-sheet protein. , 2011, Journal of molecular biology.

[109]  R. Dror,et al.  Systematic Validation of Protein Force Fields against Experimental Data , 2012, PloS one.

[110]  J. Changeux,et al.  to the activation mechanism , 2022 .