[Molecular dynamics method for proteins with ionization-conformation coupling and equilibrium titration]

Area of Science:

• Computational chemistry
• Biophysics
• Molecular modeling

Context:

• Accurate simulation of protein behavior requires accounting for ionization states and environmental factors.
• Constant pH molecular dynamics is crucial for understanding protein function in physiological conditions.
• Existing methods may lack computational efficiency or precise electrostatic energy calculations.

Purpose:

• Introduce a novel constant-pH molecular dynamics simulation method (MD-pH-ET).
• Enhance the accuracy and efficiency of calculating protein free energy and pH-dependent properties.
• Optimize protein structure and conformational analysis in aqueous solutions.

Summary:

• The MD-pH-ET method simulates proteins in their most probable ionization microstates under equilibrium titration.
• It incorporates potential of mean force corrections for accurate ensemble averaging.
• Utilizes a continuous dielectric model with Poisson equation and generalized Born method for precise electrostatic calculations.

Impact:

• Enables optimization of protein structure and total free energy calculations at constant pH.
• Facilitates the computation of pH-dependent protein properties.
• Demonstrates workability using the BPTI protein molecule, advancing molecular dynamics simulations.