Constant pH Molecular Dynamics Simulations
![Figure 1: This titration curve of the tripeptide Asp-Glu-Asp is obtained by performing constant pH MD simulations in explicit solvent at different pH values.](/627933/original-1640260188.jpg?t=eyJ3aWR0aCI6MjQ2LCJvYmpfaWQiOjYyNzkzM30%3D--fec9756074ee2c4905d3b72e67061e6fea0049ac)
![Figure 2: Scheme of the three protonation states of Glu/Asp and His. The titrating reaction coordinate (λ1) is depicted with black arrows and the switching one (λ2) with red arrows.](/627778/original-1640260188.jpg?t=eyJ3aWR0aCI6MjQ2LCJvYmpfaWQiOjYyNzc3OH0%3D--6e1605f0b4678d80ffd721c6eed4cdfd05ee2ec0)
Typically, in a molecular dynamics (MD) simulation, the protonation states of ionizable groups of a protein are set at the beginning of the simulation. This presetting is not always an obvious task, in particular for histidine, as its pKa is close to the physiological pH. In contrast, in a constant pH MD simulation, the protonation state of an ionizable group of a protein is allowed to change during the simulation according to the local electrostatic environment and the actual pH of the solution.
In this model, the protonation state of an amino acid changes continuously by propagating the motion of a virtual particle (usually referred to as “λ-particle”) that interpolates between different protonation Hamiltonians in explicit solvent MD simulations. The pKa values of the ionizable groups can then be obtained from the distributions of the protonation states.
The Three State Model
Amino acids such as aspartate (Asp), glutamate (Glu), or histidine (His) have two atoms that can carry a proton: two oxygen atoms for Asp and Glu and two nitrogen atom for His, respectively. To ensure correct representation of such residues, we developed a “three state model” that includes two protonated and one deprotonated form for Asp and Glu and one protonated and two deprotonated forms for His, respectively (Dobrev et al. 2016).
Validation
Our approach to calculate the protonation states has been tested in a recent study using challanging control pentapeptide systems. In comparison with results from NMR titration experiments, our method is suitable to not only reproduce the correct pKa values in those systems but also to determine the direction of the shift from their reference values in solution (Dobrev et al. 2020).
![Figure 3: A Illustration of the four pentapetide test systems: Gly-His-Asp-His-Gly, Gly-Glu-Asp-Glu-Gly, Gly-His-Asp-Glu-Gly, and Gly-Glu-Asp-His-Gly. The titrated residues are highlighted in the box. B Comparison of theoretical and experimental titration curves for the Gly-His-Asp-His-Gly peptide obtained in silico (black and red) and in NMR experiments (blue and orange).](/627780/original-1640260188.jpg?t=eyJ3aWR0aCI6MjQ2LCJvYmpfaWQiOjYyNzc4MH0%3D--e20b594037587f8d0b81abe3aa2c7489dfb96e99)
Simulation at constant pH
The downloadable movie is a 5 ns constant-pH MD simulation (pH=4) of the turkey ovomucoid third domain in water. The acidic amino acids Asp and Glu are drawn with sticks. They appear in red colour as negatively charged (deprotonated) and in blue colour as positively charged (protonated) residues, respectively. The proton is illustrated as a ball.
Constant pH MD simulation of turkey ovomucoid third domain in water
Constant pH MD simulation of turkey ovomucoid third domain in water
Average protonation state of the acidic amino acids of turkey ovomucoid third domain
![© MPI-BPC](/627785/original-1593525323.jpg?t=eyJ3aWR0aCI6MjQ2LCJvYmpfaWQiOjYyNzc4NX0%3D--7c58d261da04ab56aaaa3cbc15e7a8d94338aaef)
This animation shows the average protonation state (blue=protonated, red=deprotonated) of the acidic amino acids of turkey ovomucoid third domain at different pH values. The plot shows the extent of deprotonation as a function of pH: the titration curve of the protein. The higher the affinity of the protein for the protons, the more the titration curve is shifted to the right. From the titration curve, the apparent pKa values of the amino acids are obtained.