Thermodynamic coupling between neighboring binding sites in homo-oligomeric ligand sensing proteins from mass resolved ligand-dependent population distributions

Homo-oligomeric ligand-activated proteins are ubiquitous in biology. The functions of such molecules are commonly regulated by allosteric coupling between ligand-binding sites. Understanding the basis for this regulation requires both quantifying the free energy Delta G transduced between sites, and the structural basis by which it is transduced. We consider allostery in three variants of the model ring-shaped homo-oligomeric trp RNA-binding attenuation protein (TRAP). First, we developed a nearest-neighbor statistical thermodynamic binding model comprising microscopic free energies for ligand binding to isolated sites Delta G(0), and for coupling between adjacent sites, Delta G(alpha). Using the resulting partition function (PF) we explored the effects of these parameters on simulated population distributions for the 2(N) possible liganded states. We then experimentally monitored ligand-dependent population shifts using conventional spectroscopic and calorimetric methods and using native mass spectrometry (MS). By resolving species with differing numbers of bound ligands by their mass, native MS revealed striking differences in their ligand-dependent population shifts. Fitting the populations to a binding polynomial derived from the PF yielded coupling free energy terms corresponding to orders of magnitude differences in cooperativity. Uniquely, this approach predicts which of the possible 2(N) liganded states are populated at different ligand concentrations, providing necessary insights into regulation. The combination of statistical thermodynamic modeling with native MS may provide the thermodynamic foundation for a meaningful understanding of the structure-thermodynamic linkage that drives cooperativity.