Abstract:The growing interest in machine learning (ML) tools within chemistry and material science stems from their novelty and ability to predict properties almost as accurately as underlying electronic structure calculations or experiments. Transition path sampling (TPS) offers a practical way to explore transition routes between metastable minima such as conformers and isomers on the multidimensional potential energy surface. However, TPS has historically suffered from the computational cost vs. accuracy trade-off between affordable force-field simulations and expensive high-fidelity quantum mechanical calculations. ML interatomic potentials combined with TPS offer a new approach for the exploration of transition pathways at near-quantum mechanical accuracy, while keeping the computational cost comparable to classical force fields. In this study, we employ the HIP-NN-TS and ANI-1x neural network-based ML potentials, both trained on the ANI-1x dataset of 5 million HCNO structures. We first verify the correctness of our approach by applying it to alanine dipeptide and compare the resulting energy surface and transition paths to the literature. Our findings suggest that proposed approach holds promise for conformational searches, as evidenced by the chemical accuracy (errors ≲ 1 kcal/mol) for thermal molecular dynamics trajectories of alanine dipeptide. While we were able to successfully reconstruct alanine dipeptide’s potential landscape using both HIP-NN-TS and ANI-1x frameworks, we observed that ML models with lower accuracy may still recover additional important conformations. We also find that active learning, augmenting the training data by structures taken from TPS trajectories, improved the accuracy by ~30% with small amounts of additional data. Finally, we evaluated a more intricate case, azobenzene, and reiterated that seemingly simple torsions may bear a challenge for ML potentials and limit their applications in TPS. Inability of HIP-NN-TS to correctly describe the energetics of major rotational pathway in azobenzene isomerization highlights deficiencies of the reference method in describing the electronic degrees of freedom. Our study underscores the importance of domain expertise in selecting physically meaningful pathways for benchmarking ML potentials, especially considering the intricacies of electronic structure in chemical dynamics and non-equilibrium processes.

Unveiling interatomic distances influencing the reaction coordinates in alanine dipeptide isomerization: An explainable deep learning approach

Investigating the hyperparameter space of deep neural network models for reaction coordinates: Revisiting the solvent coordinate in alanine dipeptide isomerization

Learning reaction coordinates via cross-entropy minimization: Application to alanine dipeptide

Machine Learning-Guided Adaptive Parametrization for Coupling Terms in a Mixed United-Atom/Coarse-Grained Model for Diphenylalanine Self-Assembly in Aqueous Ionic Liquids

Peptide Isomerization is Suppressed at the Air-Water Interface

Three Stages of Dynamic Assembly Process of Dipeptide-Based Supramolecular Gel Revealed by In Situ Infrared Spectroscopy

Analysis of the adaptive multilevel splitting method on the isomerization of alanine dipeptide

Conformational Dynamics of an Alanine Dipeptide Analog: an Ab Initio Molecular Dynamics Study

Energy Landscapes Associated with the Self-Aggregation of an Alanine-Based Oligopeptide (Aaka)(4)

Machine learning based implicit solvent model for aqueous-solution alanine dipeptide molecular dynamics simulations

Predicting equilibrium distributions for molecular systems with deep learning

Towards Predicting Equilibrium Distributions for Molecular Systems with Deep Learning

Effect of the environment and role of the pi-pi stacking interactions in the stabilization of the 3(10)-helix conformation in dehydroalanine oligopeptides.

Extensive Conformational Searches of 13 Representative Dipeptides and an Efficient Method for Dipeptide Structure Determinations Based on Amino Acid Conformers.

Computational Exploration of Conformations of Glycine-Arginine and A Deduced Model on Global Minimum Configurations of Dipeptides in Gas Phase

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

Interpretable deep-learning pKa prediction for small molecule drugs via atomic sensitivity analysis

Amino Acid Transfer Free Energies Reveal Thermodynamic Driving Forces in Biomolecular Condensate Formation

The competing reaction mechanisms of peptide bond formation in water revealed by deep potential molecular dynamics and path sampling

Learned Mappings for Targeted Free Energy Perturbation between Peptide Conformations

Correlating molecular structure and self-assembly mechanism via Temporal Analysis of Multidimensional Chemical Interactions Space: Understanding the difference between assembly behaviors of isomeric peptides in water