Abstract
The pH-dependence of enzyme fold stability and catalytic activity is a fundamentally dynamic, structural property which is difficult to study. The challenges and expense of investigating dynamic, atomic scale behavior experimentally means that computational methods, particularly constant pH molecular dynamics (CpHMD), are well situated tools for this. However, these methods often struggle with affordable sampling of sufficiently long time scales while also obtaining accurate pK(a) prediction and verifying the structures they generate. We introduce Titr-DMD, an affordable CpHMD method that combines the quasi-all-atom coarse-grained discrete molecular dynamics (DMD) method for conformational sampling with Propka for pK(a) prediction, to circumvent these issues. The combination enables rapid sampling on limited computational resources, while simulations are still performed on the atomic scale. We benchmark the method on a set of proteins with experimentally attested pK(a) and on the pH triggered conformational change in a staphylococcal nuclease mutant, a rare experimental study of such behavior. Our results show Titr-DMD to be an effective and inexpensive method to study pH-coupled protein dynamics.