Abstract:The structural dynamics of biological macromolecules, such as proteins, DNA/RNA, or their complexes, are strongly influenced by protonation changes of their typically many titratable groups, which explains their pH sensitivity. In turn, conformational and environmental changes in the biomolecule affect the protonation state of these groups. With a few exceptions, conventional force field-based molecular dynamics (MD) simulations do not account for these effects, nor do they allow for coupling to a pH buffer. The $\lambda$-dynamics method implements this coupling and thus allows for MD simulations at constant pH. It uses separate Hamiltonians for the protonated and deprotonated states of each titratable group, with a $\lambda$ variable that continuously interpolates between them. However, rigorous implementations of Hamiltonian Interpolation (HI) $\lambda$-dynamics are prohibitively slow when used with Particle Mesh Ewald (PME). To circumvent this problem, it has been proposed to interpolate the charges instead of the Hamiltonians (QI). Here, we propose a rigorous yet efficient Multipole-Accelerated Hamiltonian Interpolation (MAHI) method to perform $\lambda$-dynamics in GROMACS. Starting from a charge-scaled Hamiltonian, precomputed with the Fast Multipole Method (FMM) or with PME, the correct HI forces are calculated with negligible computational overhead. We compare HI with QI and show that HI leads to more frequent transitions between protonation states, resulting in better sampling and accuracy. Our performance benchmarks show that introducing, e.g., 512 titratable sites to a one million atom MD system increases runtime by less than 20% compared to a regular FMM-based simulation. We have integrated the scheme into our GPU-FMM code for the simulation software GROMACS, allowing an easy and effortless transition from standard force field simulations to constant pH simulations.

Best Practices in Constant pH MD Simulations: Accuracy and Sampling

Scalable Constant pH Molecular Dynamics in GROMACS

All-Atom Continuous Constant Ph Molecular Dynamics with Particle Mesh Ewald and Titratable Water.

Constant pH Simulation with FMM Electrostatics in GROMACS. (A) Design and Applications

Biomolecular simulations at constant pH

Recent Development and Application of Constant Ph Molecular Dynamics.

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics

On the Validity of Constant pH Simulations

Development of constant-pH simulation methods in implicit solvent and applications in biomolecular systems

Constant pH Coarse-Grained Molecular Dynamics with Stochastic Charge Neutralization

Application of Explicit-Solvent Constant Ph Molecular Dynamics to Proteins

On Validity of Constant pH Simulations

An equilibrium constant pH molecular dynamics method for accurate prediction of pH-dependence in protein systems: Theory and application

Constant pH Simulation with FMM Electrostatics in GROMACS. (B) GPU Accelerated Hamiltonian Interpolation

Constant-pH Molecular Dynamics Simulations of Closed and Open States of a Proton-gated Ion Channel

How to Correct Systematic Errors in Constant-pH Ensemble Simulations

Editorial: Advanced Sampling and Modeling in Molecular Simulations for Slow and Large-Scale Biomolecular Dynamics

Mimicking titration experiments with MD simulations: A protocol for the investigation of pH-dependent effects on proteins

Effects Of System Net Charge And Electrostatic Truncation On All-Atom Constant Ph Molecular Dynamics

How to Correct Erroneous Symmetry-Breaking in Coarse-Grained Constant-pH Simulations

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field