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.

Study of Entropy-Driven Polymorphic Stability for Aspirin Using Accurate Neural Network Interatomic Potential

Revisiting Aspirin Polymorphic Stability Using a Machine Learning Potential

New Insights into the Crystallographic Disorder in the Polymorphic Forms of Aspirin from Low-Frequency Vibrational Analysis

Anharmonic and Quantum Fluctuations in Molecular Crystals: A First-Principles Study of the Stability of Paracetamol

Revealing the intergrowth phenomenon of aspirin polymorphs through a swift cooling crystallization process

Modern Semiempirical Electronic Structure Methods and Machine Learning Potentials for Drug Discovery: Conformers, Tautomers, and Protonation States

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

Modeling the Structure of Ketoprofen-Poly(vinylpyrrolidone) Amorphous Solid Dispersions with Empirical Potential Structure Refinements of X-ray Scattering Data

Exploring the crystallisation of aspirin in a confined porous material using solid-state nuclear magnetic resonance

Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs

Anharmonic effects in the low-frequency vibrational modes of aspirin and paracetamol crystals

Mechanistic Insights Into A Classic Wonder Drug-Aspirin

A complete description of thermodynamic stabilities of molecular crystals

Nucleation Control and Isolation of Polymorphic Forms of Aspirin through an Efficient Template‐Assisted Swift Cooling Process

Solution Thermodynamics of L-Glutamic Acid Polymorphs from Finite-Sized Molecular Dynamics Simulations

Thermodynamic Investigation of Carbamazepine-Saccharin Co-Crystal Polymorphs

Integrated Structural, Functional, and ADMET Analysis of 2-Methoxy-4,6-Diphenylnicotinonitrile: The Convergence of X-ray Diffraction, Molecular Docking, Dynamic Simulations, and Advanced Computational Insights

Rationalizing polymorphism with local correlation-based methods: A case study of pnictogen molecular crystals

Ab initio Accuracy Neural Network Potential for Drug-like Molecules

Enantiotropy of Simvastatin as a Result of Weakened Interactions in the Crystal Lattice: Entropy-Driven Double Transitions and the Transient Modulated Phase as Seen by Solid-State NMR Spectroscopy

Crystal Polymorph Search in the NPT Ensemble via a Deposition/Sublimation Alchemical Path