Department home
Homepage Institute
Logo  
Constant pH molecular dynamics simulations
    / home / Research Projects / Methods / Const pH MD /

Molecular dynamics simulations at constant pH

lupe
Serena Donnini, Florian Tegeler, Gerrit Groenhof, and Helmut Grubmüller

Typically, in a molecular dynamics simulation, the protonation states of ionizable groups of a protein are set in the beginning of the simulation. This 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 molecular dynamics 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 pH of the solution. The pKa values of the ionizable groups can then be obtained from the distributions of the protonation states.


Simulation at constant pH

The downloadable movie is a 5 ns constant pH simulation (pH=4) of turkey ovomucoid third domain in water. The acidic amino acids (aspartatic acid and glutamatic acid) are drawn with sticks. They are coloured red when they are negatively charged and blue when they are protonated. The proton is highlighted as a ball.


linkPfeil download mpeg-1 movie file (15 MB)
linkPfeil download QuickTime movie file (20 MB)


Average protonation state of the acidic amino acids of turkey ovomucoid third domain

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 in the movie shows the extent of deprotonation as a function of pH. This represents 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 pKas of the amino acids are obtained.

References

  1. Donnini S, Tegeler F, Groenhof G, Grubmüller H. Constant pH Molecular Dynamics in Explicit Solvent with lambda-Dynamics. J. Chem Theory and Comp 7: 1962-1978 (2011) [pdf]


linkPfeil How to use constant pH lambda dynamics
Here usage instructions are given for our "in house" version of Gromacs 3.3.1 for performing constant pH MD as described in the paper. [more]

print Print version topPfeil  Top
© 2011, Max Planck Institute for Biophysical Chemistry, Göttingen