Predicting pKa values with continuous constant pH molecular dynamics.

Area of Science:

• Biophysics
• Computational Biology
• Protein Chemistry

Background:

• Accurate prediction of protein pK(a) values is crucial for understanding protein structure-function relationships.
• Traditional macroscopic models for pK(a) calculations have limitations in describing protein environments.
• Recent advances include microscopic approaches like constant pH molecular dynamics.

Purpose of the Study:

• To introduce and detail the Continuous Constant pH Molecular Dynamics (CPHMD) methodology.
• To demonstrate the accuracy and robustness of CPHMD for predicting protein pK(a)s.
• To explore the link between ionization equilibria and pH-modulated protein dynamics.

Main Methods:

• Implementation of the Continuous Constant pH Molecular Dynamics (CPHMD) method.
• Application of CPHMD to calculate pK(a) values for various protein residues (surface, buried, flexible).
• Comparison of CPHMD results with experimental data and other theoretical methods (PB-based, empirical).

Main Results:

• CPHMD accurately predicts pK(a) values for diverse protein systems, including ribonuclease A, staphylococcal nuclease, and α-lactalbumin.
• The method effectively captures the correlation between ionization equilibria and conformational dynamics.
• CPHMD demonstrates robustness in accounting for local dielectric responses during structural changes.

Conclusions:

• CPHMD is a powerful and accurate tool for predicting protein pK(a)s.
• The methodology provides insights into pH-dependent protein behavior and dynamics.
• CPHMD offers a robust framework for future advancements in computational protein biophysics.