Abstract:We present a GPU implementation of the continuous constant pH molecular dynamics (CpHMD) based on the most recent generalized Born implicit-solvent model in the pmemd engine of the Amber molecular dynamics package. To test the accuracy of the tool for rapid pKa predictions, a series of 2 ns single-pH simulations were performed for over 120 titratable residues in 10 benchmark proteins that were previously used to test the various continuous CpHMD methods. The calculated pKa's showed a root-mean-square deviation of 0.80 and correlation coefficient of 0.83 with respect to experiment. Also, 90% of the pKa's were converged with estimated errors below 0.1 pH units. Surprisingly, this level of accuracy is similar to our previous replica-exchange simulations with 2 ns per replica and an exchange attempt frequency of 2 ps-1 (Huang, Harris, and Shen J. Chem. Inf. Model. 2018 , 58 , 1372 - 1383 ). Interestingly, for the linked titration sites in two enzymes, although residue-specific protonation state sampling in the single-pH simulations was not converged within 2 ns, the protonation fraction of the linked residues appeared to be largely converged, and the experimental macroscopic pKa values were reproduced to within 1 pH unit. Comparison with replica-exchange simulations with different exchange attempt frequencies showed that the splitting between the two macroscopic pKa's is underestimated with frequent exchange attempts such as 2 ps-1, while single-pH simulations overestimate the splitting. The same trend is seen for the single-pH vs replica-exchange simulations of a hydrogen-bonded aspartyl dyad in a much larger protein. A 2 ns single-pH simulation of a 400-residue protein takes about 1 h on a single NVIDIA GeForce RTX 2080 graphics card, which is over 1000 times faster than a CpHMD run on a single CPU core of a high-performance computing cluster node. Thus, we envision that GPU-accelerated continuous CpHMD may be used in routine pKa predictions for a variety of applications, from assisting MD simulations with protonation state assignment to offering pH-dependent corrections of binding free energies and identifying reactive hot spots for covalent drug design.

Application of Explicit-Solvent Constant Ph Molecular Dynamics to Proteins

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

Including Explicit Water Molecules as Part of the Protein Structure in MM/PBSA Calculations.

Introducing titratable water to all-atom molecular dynamics at constant pH.

Recent Development and Application of Constant Ph Molecular Dynamics.

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

Constant-pH Molecular Dynamics Simulations for Large Biomolecular Systems.

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

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

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

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

Biomolecular simulations at constant pH

On the Validity of Constant pH Simulations

OPEP6: A New Constant-pH Molecular Dynamics Simulation Scheme with OPEP Coarse-Grained Force Field

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics

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

Force Field Limitations of All-Atom Continuous Constant pH Molecular Dynamics

On Validity of Constant pH Simulations

GPU-Accelerated Implementation of Continuous Constant pH Molecular Dynamics in Amber: p Ka Predictions with Single-pH Simulations

Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model.