Abstract
Like temperature and pressure, the solution pH is an important thermodynamic variable that is commonly varied in experiments and is used by cells to influence biochemical function. It is now becoming feasible to carry out practical molecular dynamics simulations that mimic the thermodynamics of such experiments, by allowing proton transfer between the system of interest and a hypothetical bath of protons at a given pH. These are demanding calculations, because the energetics of charge changes upon protonation or deprotonation must be accurately modeled, and because such simulations must sample both molecular configurations and the large number of protonation states that are possible for a molecule with many titrating sites.