Resolving coupled pH titrations using alchemical free energy calculations

Abstract

In a protein, nearby titratable sites can be coupled: the (de)protonation of one may affect the other. The degree of this interaction depends on several factors and can influence the measured (Formula presented.). Here, we derive a formalism based on double free energy differences ((Formula presented.)) for quantifying the individual site (Formula presented.) values of coupled residues. As (Formula presented.) values can be obtained by means of alchemical free energy calculations, the presented approach allows for a convenient estimation of coupled residue (Formula presented.) s in practice. We demonstrate that our approach and a previously proposed microscopic (Formula presented.) formalism, can be combined with alchemical free energy calculations to resolve pH-dependent protein (Formula presented.) values. Toy models and both, regular and constant-pH molecular dynamics simulations, alongside experimental data, are used to validate this approach. Our results highlight the insights gleaned when coupling and microstate probabilities are analyzed and suggest extensions to more complex enzymatic contexts. Furthermore, we find that naïvely computed (Formula presented.) values that ignore coupling, can be significantly improved when coupling is accounted for, in some cases reducing the error by half. In short, alchemical free energy methods can resolve the (Formula presented.) values of both uncoupled and coupled residues.