Abstract
Allostery pervades macromolecular function and drives cooperative binding of ligands to macromolecules. To decipher the mechanisms of cooperative ligand binding it is necessary to define at a microscopic level the structural and thermodynamic consequences of binding of each ligand to its allosterically coupled site(s). However, dynamic sampling of alternative conformations (microstates) in allosteric molecules complicates interpretation of both structural and thermodynamic data. Isothermal titration calorimetry has the potential to directly quantify the thermodynamics of allosteric interactions, but usually falls short of enabling mechanistic insight. This is because 1) its measurements reflect the sum of overlapping caloric processes involving binding-linked population shifts within and between microstates, and 2) data are generally fit with phenomenological binding polynomials that are underdetermined. Nevertheless, temperature-dependent binding data have the potential to resolve overlapping thermodynamic processes, while mechanistically constrained models enable hypothesis testing and identification of informative parameters. We globally fit temperature-dependent isothermal titration calorimetry data for binding of 11 tryptophan ligands to the homo-undecameric trp RNA-binding Attenuation Protein from Bacillus stearothermophilus using nearest-neighbor statistical thermodynamic models. This approach allowed us to distinguish alternative nearest-neighbor interaction models, and quantifies the thermodynamic contribution of neighboring ligands to individual binding sites. We also perform conventional Hill equation modeling and illustrate how comparatively limited it is in quantitative or mechanistic value. This work illustrates the potential of mechanistically constrained global fitting of binding data to yield the microscopic thermodynamic parameters essential for deciphering mechanisms of cooperativity in a wide range of ligand-regulated homo-oligomeric assemblies.