Abstract:Many drug molecules contain biaryl fragments, resulting in a torsional barrier corresponding to rotation around the bond linking the aryls. The potential energy surfaces of these torsions vary significantly due to steric and electronic effects, ultimately affecting the relative stability of the molecular conformations in the protein-bound and solution states. Simulations of protein--ligand binding require accurate computational models to represent the intramolecular interactions to provide accurate predictions of the structure and dynamics of binding. In this paper, we compare four force fields (Generalized AMBER Force Field (GAFF), Open Force Field (OpenFF), CHARMM General Force Field (CGenFF), Optimized Potentials for Liquid Simulations (OPLS)) and two neural network potentials (ANI-2x, ANI-1ccx) in their ability to predict the torsional potential energy surfaces of 88 biaryls extracted from drug fragments. The mean of the absolute deviation over the full PES (MADF) and the mean absolute deviation of the torsion rotational barrier height (MADB) relative to high-level ab initio reference data was used as a measure of accuracy. In comparison to high-level ab-initio data, ANI-1ccx was most accurate for predicting the barrier height (MADF: 0.5~kcal/mol, MADB:~0.8~kcal/mol), followed closely by ANI-2x (MADF: 0.5~kcal/mol, MADB:~1.0~kcal/mol), then CGenFF (MADF: 0.8~kcal/mol, MADB: 1.3~kcal/mol), OpenFF (MADF: 1.5~kcal/mol, MADB: 1.4~kcal/mol), GAFF (MADF: 1.2~kcal/mol, MADB: 2.6~kcal/mol), and finally OPLS (MADF: 1.5~kcal/mol, MADB: 2.8~kcal/mol). Significantly, the NNPs are systematically more accurate and more reliable than any of the force fields. As a practical example, the neural network potential/molecular mechanics (NNP/MM) method was used to simulate the isomerization of ozanimod, a drug used for multiple sclerosis. Multi-nanosecond molecular dynamics (MD) simulations in an explicit aqueous solvent were performed, as well as umbrella sampling and adaptive biasing force enhanced sampling techniques. These theories predicted a rate of isomerization of $4.30 \times 10^{-1}$ ns$^{-1}$, which is consistent with direct MD simulations.

Reweighting from Molecular Mechanics Force Fields to the ANI-2x Neural Network Potential

Reweighting from Molecular Mechanics Force Fields to the ANI-2x Neural Network Potential

Neural Network Potential with Multi-Resolution Approach Enables Accurate Prediction of Reaction Free Energies in Solution

Incorporating Neural Networks into the AMOEBA Polarizable Force Field

Insights and Challenges in Correcting Force Field Based Solvation Free Energies Using A Neural Network Potential

LAMMPS Implementation of Rapid Artificial Neural Network Derived Interatomic Potentials

Benchmarking Force Field and the ANI Neural Network Potentials for the Torsional Potential Energy Surface of Biaryl Drug Fragments

Ab initio Accuracy Neural Network Potential for Drug-like Molecules

ANI-1: an extensible neural network potential with DFT accuracy at force field computational cost

Enhancing Molecular Energy Predictions with Physically Constrained Modifications to the Neural Network Potential

ANI/EFP: Modeling Long-Range Interactions in ANI Neural Network with Effective Fragment Potentials

FeNNol: an Efficient and Flexible Library for Building Force-field-enhanced Neural Network Potentials

Scalable Training of Neural Network Potentials for Complex Interfaces Through Data Augmentation

Transfer learning for chemically accurate interatomic neural network potentials

Fast and Sample-Efficient Interatomic Neural Network Potentials for Molecules and Materials Based on Gaussian Moments

Graph Neural Network Interatomic Potential Ensembles with Calibrated Aleatoric and Epistemic Uncertainty on Energy and Forces

Neural Network Potentials for Chemistry: Concepts, Applications and Prospects

Comparing ANI-2x, ANI-1ccx neural networks, force field, and DFT methods for predicting conformational potential energy of organic molecules

NNP/MM: Accelerating Molecular Dynamics Simulations with Machine Learning Potentials and Molecular Mechanics

NNP/MM: Accelerating molecular dynamics simulations with machine learning potentials and molecular mechanic