TorchMD: A Deep Learning Framework for Molecular Simulations

Molecular dynamics simulations provide a mechanistic description of molecules by relying on empirical potentials. The quality and transferability of such potentials can be improved leveraging data-driven models derived with machine learning approaches. Here, we present TorchMD, a framework for molecular simulations with mixed classical and machine learning potentials. All force computations including bond, angle, dihedral, Lennard-Jones, and Coulomb interactions are expressed as PyTorch arrays and operations. Moreover, TorchMD enables learning and simulating neural network potentials. We validate it using standard Amber all-atom simulations, learning an ab initio potential, performing an end-to-end training, and finally learning and simulating a coarse-grained model for protein folding. We believe that TorchMD provides a useful tool set to support molecular simulations of machine learning potentials. Code and data are freely available at github.com/torchmd.

[1]  J. Chodera,et al.  End-to-End Differentiable Molecular Mechanics Force Field Construction , 2020, ArXiv.

[2]  Wujie Wang,et al.  Temperature-transferable coarse-graining of ionic liquids with dual graph convolutional neural networks. , 2020, The Journal of chemical physics.

[3]  Frank Noé,et al.  Coarse graining molecular dynamics with graph neural networks. , 2020, The Journal of chemical physics.

[4]  Xiang Gao,et al.  TorchANI: A Free and Open Source PyTorch-Based Deep Learning Implementation of the ANI Neural Network Potentials , 2020, J. Chem. Inf. Model..

[5]  Wujie Wang,et al.  Active learning and neural network potentials accelerate molecular screening of ether-based solvate ionic liquids. , 2020, Chemical communications.

[6]  Cecilia Clementi,et al.  Ensemble learning of coarse-grained molecular dynamics force fields with a kernel approach. , 2020, The Journal of chemical physics.

[7]  Lars A. Bratholm,et al.  FCHL revisited: Faster and more accurate quantum machine learning. , 2019, The Journal of chemical physics.

[8]  Joel Nothman,et al.  SciPy 1.0-Fundamental Algorithms for Scientific Computing in Python , 2019, ArXiv.

[9]  E. D. Cubuk,et al.  JAX, M.D.: End-to-End Differentiable, Hardware Accelerated, Molecular Dynamics in Pure Python , 2019, 1912.04232.

[10]  Feliks Nüske,et al.  Coarse-graining molecular systems by spectral matching. , 2019, The Journal of chemical physics.

[11]  S. Pantano,et al.  The SIRAH 2.0 Force Field: Altius, Fortius, Citius. , 2019, Journal of chemical theory and computation.

[12]  Wujie Wang,et al.  Coarse-graining auto-encoders for molecular dynamics , 2018, npj Computational Materials.

[13]  Frank Noé,et al.  Machine Learning of Coarse-Grained Molecular Dynamics Force Fields , 2018, ACS central science.

[14]  Hao Wu,et al.  Variational selection of features for molecular kinetics. , 2018, The Journal of chemical physics.

[15]  K-R Müller,et al.  SchNetPack: A Deep Learning Toolbox For Atomistic Systems. , 2018, Journal of chemical theory and computation.

[16]  Patrice Koehl,et al.  Faculty Opinions recommendation of OpenMM 7: Rapid development of high performance algorithms for molecular dynamics. , 2018, Faculty Opinions – Post-Publication Peer Review of the Biomedical Literature.

[17]  David L. Mobley,et al.  Open Force Field Consortium: Escaping atom types using direct chemical perception with SMIRNOFF v0.1 , 2018, bioRxiv.

[18]  Mohammad M. Sultan,et al.  Automated design of collective variables using supervised machine learning. , 2018, The Journal of chemical physics.

[19]  Linfeng Zhang,et al.  DeePCG: Constructing coarse-grained models via deep neural networks. , 2018, The Journal of chemical physics.

[20]  K-R Müller,et al.  SchNet - A deep learning architecture for molecules and materials. , 2017, The Journal of chemical physics.

[21]  David W Toth,et al.  The TensorMol-0.1 model chemistry: a neural network augmented with long-range physics , 2017, Chemical science.

[22]  V. Pande,et al.  Modeling the mechanism of CLN025 beta-hairpin formation. , 2017, The Journal of chemical physics.

[23]  Toni Giorgino,et al.  Drug Discovery and Molecular Dynamics: Methods, Applications and Perspective Beyond the Second Timescale. , 2017, Current topics in medicinal chemistry.

[24]  Klaus-Robert Müller,et al.  SchNet: A continuous-filter convolutional neural network for modeling quantum interactions , 2017, NIPS.

[25]  V. Pande,et al.  Modeling the Mechanism of CLN025 Beta-Hairpin Formation , 2017, bioRxiv.

[26]  Vijay S. Pande,et al.  OpenMM 7: Rapid development of high performance algorithms for molecular dynamics , 2016, bioRxiv.

[27]  Mohammad M. Sultan,et al.  Optimized parameter selection reveals trends in Markov state models for protein folding. , 2016, The Journal of chemical physics.

[28]  Jian Yin,et al.  Lessons learned from comparing molecular dynamics engines on the SAMPL5 dataset , 2016, bioRxiv.

[29]  Alán Aspuru-Guzik,et al.  Convolutional Networks on Graphs for Learning Molecular Fingerprints , 2015, NIPS.

[30]  C. Simmerling,et al.  ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. , 2015, Journal of chemical theory and computation.

[31]  Jimmy Ba,et al.  Adam: A Method for Stochastic Optimization , 2014, ICLR.

[32]  Pavlo O. Dral,et al.  Quantum chemistry structures and properties of 134 kilo molecules , 2014, Scientific Data.

[33]  S. Doerr,et al.  On-the-Fly Learning and Sampling of Ligand Binding by High-Throughput Molecular Simulations. , 2014, Journal of chemical theory and computation.

[34]  W G Noid,et al.  Perspective: Coarse-grained models for biomolecular systems. , 2013, The Journal of chemical physics.

[35]  D. Tieleman,et al.  Perspective on the Martini model. , 2013, Chemical Society reviews.

[36]  Marissa G. Saunders,et al.  Coarse-graining methods for computational biology. , 2013, Annual review of biophysics.

[37]  Toni Giorgino,et al.  Identification of slow molecular order parameters for Markov model construction. , 2013, The Journal of chemical physics.

[38]  Vijay S Pande,et al.  Simple few-state models reveal hidden complexity in protein folding , 2012, Proceedings of the National Academy of Sciences.

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

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

[41]  David P. Anderson,et al.  High-Throughput All-Atom Molecular Dynamics Simulations Using Distributed Computing , 2010, J. Chem. Inf. Model..

[42]  Klaus Schulten,et al.  Discovery through the computational microscope. , 2009, Structure.

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

[44]  Shinya Honda,et al.  Crystal structure of a ten-amino acid protein. , 2008, Journal of the American Chemical Society.

[45]  Cecilia Clementi,et al.  Coarse-grained models of protein folding: toy models or predictive tools? , 2008, Current opinion in structural biology.

[46]  Gregory A Voth,et al.  A multiscale coarse-graining method for biomolecular systems. , 2005, The journal of physical chemistry. B.

[47]  Brandsdal,et al.  [Advances in Protein Chemistry] Protein Simulations Volume 66 || Free Energy Calculations and Ligand Binding , 2003 .

[48]  Vladimir N. Vapnik,et al.  The Nature of Statistical Learning Theory , 2000, Statistics for Engineering and Information Science.

[49]  Berk Hess,et al.  Improving Efficiency of Large Time-Scale Molecular Dynamics Simulations of Hydrogen-Rich Systems , 1999 .

[50]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

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

[52]  J. A. V. BUTLER,et al.  Advances in Protein Chemistry , 1946, Nature.

[53]  A. Given,et al.  The Time Machine , 1895, Nature.