A Guide to the Continuous Constant pH Molecular Dynamics Methods in Amber and CHARMM [Article v1.0]

Area of Science:

• Biomolecular simulations
• Computational chemistry
• Protein dynamics

Background:

• Solution pH is a critical environmental variable influencing protein structure and function.
• Conventional molecular dynamics (MD) simulations often use fixed protonation states for titratable sidechains, which can be inaccurate and limit dynamic studies.
• Accurate representation of protonation states is essential for understanding protein behavior.

Purpose of the Study:

• To provide a tutorial on utilizing continuous constant pH MD methods.
• To demonstrate the application of these methods in predicting pKa values.
• To elucidate proton-coupled conformational dynamics in proteins.

Main Methods:

• Continuous constant pH molecular dynamics (MD) simulations.
• Utilizing implementations within the Amber and CHARMM simulation packages.
• Applying methods to predict pKa values and analyze protein dynamics.

Main Results:

• Demonstrated the utility of constant pH MD for accurate protein simulations.
• Successfully applied methods to predict pKa values for various proteins.
• Elucidated proton-coupled conformational changes influencing protein function.

Conclusions:

• Constant pH MD methods overcome limitations of fixed protonation states in biomolecular simulations.
• These advanced simulation techniques are valuable for studying protein structure-function relationships.
• The tutorial facilitates broader adoption of constant pH MD for protein research.