Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2024 Jan 11;128(1):109-116.
doi: 10.1021/acs.jpcb.3c06662. Epub 2023 Dec 28.

OpenMM 8: Molecular Dynamics Simulation with Machine Learning Potentials

Affiliations

OpenMM 8: Molecular Dynamics Simulation with Machine Learning Potentials

Peter Eastman et al. J Phys Chem B. .

Abstract

Machine learning plays an important and growing role in molecular simulation. The newest version of the OpenMM molecular dynamics toolkit introduces new features to support the use of machine learning potentials. Arbitrary PyTorch models can be added to a simulation and used to compute forces and energy. A higher-level interface allows users to easily model their molecules of interest with general purpose, pretrained potential functions. A collection of optimized CUDA kernels and custom PyTorch operations greatly improves the speed of simulations. We demonstrate these features in simulations of cyclin-dependent kinase 8 (CDK8) and the green fluorescent protein chromophore in water. Taken together, these features make it practical to use machine learning to improve the accuracy of simulations with only a modest increase in cost.

PubMed Disclaimer

Figures

Figure 1:
Figure 1:
The 53 atom inhibitor bound to CDK8.
Figure 2:
Figure 2:
Anionic GFP chromophore in water. The dihedral ɸP is defined between the CP atom and the carbons of the methine bridge and the dihedral ɸI is defined between the NI atom and the carbons of the methine bridge.
Figure 3:
Figure 3:
Equivariant transformer training and validation loss.
Figure 4:
Figure 4:
Anionic oxygen-oxygen RDF (left) and carbonyl oxygen-oxygen RDF (right) generated from the MLP and AIMD trajectories. The shaded region represents +/− 1 standard deviation across the 10 AIMD trajectories.
Figure 5:
Figure 5:
Distribution of dihedral angles for ɸI and ɸP. The shaded region represents +/− 1 standard deviation across the 10 AIMD trajectories.

Update of

References

    1. Eastman P; Swails J; Chodera JD; McGibbon RT; Zhao Y; Beauchamp KA; Wang L-P; Simmonett AC; Harrigan MP; Stern CD; et al. OpenMM 7: Rapid Development of High Performance Algorithms for Molecular Dynamics. PLOS Comput. Biol 2017, 13 (7), e1005659. 10.1371/journal.pcbi.1005659. - DOI - PMC - PubMed
    1. Smith JS; Isayev O; Roitberg AE ANI-1: An Extensible Neural Network Potential with DFT Accuracy at Force Field Computational Cost. Chem. Sci 2017, 8 (4), 3192–3203. 10.1039/C6SC05720A. - DOI - PMC - PubMed
    1. Batzner S; Musaelian A; Sun L; Geiger M; Mailoa JP; Kornbluth M; Molinari N; Smidt TE; Kozinsky B E(3)-Equivariant Graph Neural Networks for Data-Efficient and Accurate Interatomic Potentials. Nat. Commun 2022, 13 (1), 2453. 10.1038/s41467-022-29939-5. - DOI - PMC - PubMed
    1. Thölke P; De Fabritiis G TorchMD-NET: Equivariant Transformers for Neural Network Based Molecular Potentials. arXiv April 23, 2022. 10.48550/arXiv.2202.02541. - DOI
    1. Musaelian A; Batzner S; Johansson A; Sun L; Owen CJ; Kornbluth M; Kozinsky B Learning Local Equivariant Representations for Large-Scale Atomistic Dynamics. Nat. Commun 2023, 14 (1), 579. 10.1038/s41467-023-36329-y. - DOI - PMC - PubMed

Publication types