Abstract:Our study focused on the implementation and testing of machine learning interatomic potentials (MLIPs) into the AMBER software suite. This implementation enables us to perform a novel type of molecular dynamics simulation utilizing the hybrid machine learning/molecular mechanics (ML/MM) potentials. To underpin the capabilities of ML/MM simulations, we first validated our implementation at a fundamental physical level by confirming energy and momentum conservation laws. The successful validation indicates that our implementation is able to produce adequate and physically interpretable samplings. Building upon this, for the first time to the best of our knowledge, we proposed an ML/MM-compatible thermodynamic integration (TI) protocol to tackle real-world challenges, such as solvation free energy calculation. Our results demonstrate that this computational protocol can predict hydration free energies with an accuracy of less than 1.00 kcal/mol compared to experimental data, paving the way for the use of ML/MM in multiscale simulations to addressing future drug design problems. Moreover, by applying ML/MM in molecular dynamics simulations of protein-ligand complexes, we demonstrated that the adequate samplings enable us to accurately reproduce experimental binding free energies. Thus, our implementation can offer new insights into biomolecular systems using the ML/MM "microscope". Last, we demonstrated that our implementation can achieve nanosecond timescale simulations daily after significant effort being put to improve the code performance. In a conclusion, we have successfully implemented ML/MM potential to AMBER software package after overcoming limitations in current multi-scale simulations including low computational efficiency. We have advanced TI theory allowing us to accurately predict free energies with ML/MM potentials.

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman-Rubin Convergence Diagnostics

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman-Rubin Convergence Diagnostics

Coupling Accelerated Molecular Dynamics Methods with Thermodynamic Integration Simulations

Fast Computation of Solvation Free Energies with Molecular Density Functional Theory: Thermodynamic-Ensemble Partial Molar Volume Corrections

An Approach to Computing Solvent Reorganization Energy.

On delivering polar solvation free energy of proteins from energy minimized structures using a regularized super‐Gaussian Poisson–Boltzmann model

A combined steepest descent and genetic algorithm (SD/GA) approach for the optimization of solvation parameters

Solvation Free Energy Calculations with Quantum Mechanics/Molecular Mechanics and Machine Learning Models.

Ion solvation free energy calculations based on first-principles molecular dynamics thermodynamic integration

Calculations Of Free Energies Based On The Combination Of Molecular Dynamics Simulations And Continuum Solvation Model

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Solvation Free Energies from Machine Learning Molecular Dynamics

Effect of the Integration Method on the Accuracy and Computational Efficiency of Free Energy Calculations Using Thermodynamic Integration

Calculation of solvation force in molecular dynamics simulation by deep-learning method

Enhancing Accuracy and Feature Insights in Hydration Free Energy Predictions for Small Molecules with Machine Learning

Computing spatially resolved rotational hydration entropies from atomistic simulations

Free Energy Calculations Using the Movable Type Method with Molecular Dynamics Driven Protein-Ligand Sampling

Electrostatic solvation free energies of charged hard spheres using molecular dynamics with density functional theory interactions

Accurate Free Energy Calculation via Multiscale Simulations Driven by Hybrid Machine Learning and Molecular Mechanics Potentials

A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals

Modeling and analysis of ensemble average solvation energy and solute-solvent interfacial fluctuations