Analytic gradients for the electrostatic embedding QM/MM in periodic boundary conditions using particle-mesh Ewald sums and electrostatic potential fitted charge operators

Miquel Huix-Rotllant,Simone Bonfrate,Nicolas Ferré

DOI: https://doi.org/10.26434/chemrxiv-2024-5tzsf

2024-02-20

Abstract:Long-range electrostatic effects are fundamental for describing chemical reactivity in the condensed phase. Here, we present the methodology of an efficient quantum mechanical/molecular mechanical (QM/MM) model in periodic boundary conditions (PBC) compatible with QM/MM boundaries at chemical bonds. The method combines electrostatic potential fitted (ESPF) charge operators and electrostatic potentials derived from the smooth particle-mesh Ewald (PME) sum approach. The total energy and its analytic first derivatives with respect to QM, MM and lattice vectors allow QM/MM molecular dynamics (MD) in the most common thermodynamic ensembles. We demonstrate the robustness of the method by performing a QM/MM MD equilibration of methanol in water. We simulate the cis/trans isomerization free energy profiles in water of proline amino acid and a proline-containing oligopeptide, showing a correct description of the reaction barrier. Our PBC-compatible QM/MM model can efficiently be used to study chemical reactivity in condensed phase and enzymatic catalysis.

Chemistry

What problem does this paper attempt to address?

This paper mainly discusses the analysis gradient of the Quantum Mechanics/Molecular Mechanics (QM/MM) embedding model using the Particle Mesh Ewald (PME) method and the Potential Fitting Charge Operator under Periodic Boundary Conditions (PBC). The authors propose an effective method that is compatible with the QM/MM boundary at chemical bonds and can describe chemical reactivity in condensed phases. This method combines the efficiency of the PME method and the advantages of the Electrostatic Potential (ESP) charge operator to approximately represent the behavior of the QM density in replicas. By calculating the total energy and its gradients with respect to QM, MM atomic coordinates, and lattice parameters, QM/MM molecular dynamics simulations can be performed under the most common thermodynamic states. The paper first introduces the importance of the QM/MM embedding method in describing the reactivity of biomolecules, emphasizing the necessity of accurately describing long-range electrostatic interactions for calculating properties such as free energy changes. The authors point out the incomplete long-range interactions in existing methods and propose an improved approach based on Ewald summation methods, particularly using the Smooth Particle Mesh Ewald (sPME) method and the ESP charge operator. The paper provides a detailed description of using the ESP method to describe the electrostatic interactions between QM and MM subsystems, including the expressions for energy and energy gradient, as well as the derivatives with respect to QM, MM atomic coordinates, and lattice parameters. In addition, the paper discusses the treatment of QM-MM boundaries connected by covalent bonds, such as the Link Atom method, to address the problem of overpolarized unsaturated atoms between QM and MM regions. Finally, the paper validates the stability and efficiency of the method through QM/MM MD equilibrium simulations of methanol in water and free energy profiles simulations of the amino acid proline and proline-containing oligopeptides in water. This demonstrates that the method can effectively be used to study condensed phases and enzymatic catalysis.

Analytic gradients for the electrostatic embedding QM/MM in periodic boundary conditions using particle-mesh Ewald sums and electrostatic potential fitted charge operators

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

Assessment of free energies from electrostatic embedding density functional tight binding-based/molecular mechanics in periodic boundary conditions

Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis

Accelerated Quantum Mechanics/Molecular Mechanics Simulations Via Neural Networks Incorporated with Mechanical Embedding Scheme.

An Electrostatic Energy-Based Charge Model For Molecular Dynamics Simulation

Efficient polarizable QM/MM using the direct reaction field Hamiltonian with electrostatic potential fitted multipole operators

Multistate, polarizable QM/MM embedding scheme based on the direct reaction field method: Solvatochromic shifts, analytical gradients and optimization of conical intersections in solution

Assessment of embedding schemes in a hybrid machine learning/classical potentials (ML/MM) scheme

Learning the Quantum Centroid Force Correction in Molecular Systems: A Localized Approach

Analytic Energy Gradient of Excited Electronic State Within TDDFT/MMpol Framework: Benchmark Tests and Parallel Implementation.

Electrostatic Embedding Machine Learning for Ground and Excited State Molecular Dynamics of Solvated Molecules

Efficient Formulation of Polarizable Gaussian Multipole Electrostatics for Biomolecular Simulations

Assessment of Embedding Schemes in a Hybrid Machine Learning/Classical Potentials (ML/MM) Approach

Asymptotic Expansion for Electrostatic Embedding Integrals in QM/MM Calculations.

Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis.

An Electrostatically Embedded QM/MM Scheme for Electrified Interfaces

A Quantum-Mechanics Molecular-Mechanics scheme for extended systems

Development of Effective Quantum Mechanical/molecular Mechanical (QM/MM) Methods for Complex Biological Processes.

Seamless integration of GEM, a density based-force field, for QM/MM simulations via LICHEM, Psi4 and Tinker-HP