Abstract:An efficient yet accurate method for producing a large amount of energy data for molecular mechanical force field (MMFF) parameterization is on demand, especially for torsional angle parameters which are typically derived to reproduce ab initio rotational profiles or torsional potential energy surfaces (PESs). Recently, an active learning potential (ANI-1x) for organic molecules which can produce smooth and physically meaningful PESs has been developed. The high efficiency and accuracy make ANI-1x especially attractive for geometry optimization at low cost. To apply the ANI-1x potential in MMFF parameterization, one needs to perform constrained geometry optimization. In this work, we first developed a computational protocol to constrain rotatable torsional angles and other geometric parameters for a molecule whose geometry is described by Cartesian coordinates. The constraint is successfully achieved by force projection for the two conjugated gradient (CG) algorithms. We then conducted large-scale assessments on ANI-1x along with four different optimization algorithms in reproducing DFT energies and geometries for two CG algorithms, CG backtracking line search (CG-BS) and CG Wolfe line search (CG-WS), and two quasi-Newton algorithms, Broyden-Fletcher-Goldfarb-Shanno (BFGS) and low-memory BFGS (L-BFGS). Note that CG-BS is a new algorithm we developed in this work. All four algorithms take the ANI energies and forces to optimize a molecule geometry. Last, we conducted a large-scale assessment of applying ANI-1x in MMFF development in three aspects. First, we performed full optimizations for 100 drug molecules, each consisting of five distinct conformations. The average root-mean-square error (RMSE) between ANI-1x and DFT is about 1.3 kcal/mol, and the root-mean-square displacement (RMSD) of heavy atoms is about 0.35 Å. Second, we generated torsional PESs for 160 organic molecules, and constrained optimizations were performed for up to 18 conformations for each PES. We found that the RMSE of all the conformers is 1.23 kcal/mol. Last, we carried out constrained optimizations for alanine dipeptide with both ϕ and φ angles being frozen. The Ramachandran plots indicate that the two CG algorithms in conjunction with the ANI-1x potential could well reproduce the DFT-optimized geometries and torsional PESs. We concluded that CG-BS and CG-WS are good choices for generating PESs, while CG-WS or BFGS is ideal for performing full geometry optimization. With the continuously increased quality of ANI, it is expected that the computational algorithms and protocols presented in this work will have great applications in improving the quality of an existing small-molecule MMFF.

DeepConf: Leveraging ANI-ML Potentials for Exploring Local Minima with A Focus on Bioactive Conformations

ComABAN: refining molecular representation with the graph attention mechanism to accelerate drug discovery

ConfGen: A Conformational Search Method for Efficient Generation of Bioactive Conformers

Small-Molecule Conformer Generators: Evaluation of Traditional Methods and AI Models on High-Quality Data Sets

Development and Evaluation of Geometry Optimization Algorithms in Conjunction with ANI Potentials.

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

An optimization algorithm for conformer generation based on the bond contribution ranking

Auto3D: Automatic Generation of the Low-Energy 3D Structures with ANI Neural Network Potentials

Do Deep Learning Methods Really Perform Better in Molecular Conformation Generation?

First-Principles Molecular Structure Search with a Genetic Algorithm

DeepMol: An Automated Machine and Deep Learning Framework for Computational Chemistry

Challenges and Opportunities for Machine Learning Potentials in Transition Path Sampling: Alanine Dipeptide and Azobenzene Studies

Less is more: sampling chemical space with active learning

Geometry Optimization Algorithms in Conjunction with the Machine Learning Potential ANI-2x Facilitate the Structure-Based Virtual Screening and Binding Mode Prediction

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

MoleCLUEs: Molecular Conformers Maximally In-Distribution for Predictive Models

Benchmarking ANI potentials as a rescoring function and screening FDA drugs for SARS-CoV-2 Mpro

Bioactive Conformational Generation of Small Molecules: A Comparative Analysis Between Force-Field and Multiple Empirical Criteria Based Methods

AIMD-Chig: Exploring the conformational space of a 166-atom protein Chignolin with ab initio molecular dynamics

Assessing Conformation Validity and Rationality of Deep Learning-Generated 3D Molecules