Changes in small-angle X-ray scattering parameters observed upon binding of ligand to rabbit muscle pyruvate kinase are not correlated with allosteric transitions

What is allosteric regulation? Exploring the exceptions that prove the rule!

"Allosteric" was first introduced to mean the other site (i.e., a site distinct from the active or orthosteric site), an adjective for "regulation" to imply a regulatory outcome resulting from ligand binding at another site. That original idea outlines a system with two ligand-binding events at two distinct locations on a macromolecule (originally a protein system), which defines a four-state energy cycle. An allosteric energy cycle provides a quantifiable allosteric coupling constant and focuses our attention on the unique properties of the four equilibrated protein complexes that constitute the energy cycle. Because many observed phenomena have been referenced as "allosteric regulation" in the literature, the goal of this work is to use literature examples to explore which systems are and are not consistent with the two-ligand thermodynamic energy cycle-based definition of allosteric regulation. We emphasize the need for consistent language so comparisons can be made among the ever-increasing number of allosteric systems. Building on the mutually exclusive natures of an energy cycle definition of allosteric regulation versus classic two-state models, we conclude our discussion by outlining how the often-proposed Rube-Goldberg-like mechanisms are likely inconsistent with an energy cycle definition of allosteric regulation.

PYK-SubstitutionOME: an integrated database containing allosteric coupling, ligand affinity and mutational, structural, pathological, bioinformatic and computational information about pyruvate kinase isozymes

Interpreting changes in patient genomes, understanding how viruses evolve and engineering novel protein function all depend on accurately predicting the functional outcomes that arise from amino acid substitutions. To that end, the development of first-generation prediction algorithms was guided by historic experimental datasets. However, these datasets were heavily biased toward substitutions at positions that have not changed much throughout evolution (i.e. conserved). Although newer datasets include substitutions at positions that span a range of evolutionary conservation scores, these data are largely derived from assays that agglomerate multiple aspects of function. To facilitate predictions from the foundational chemical properties of proteins, large substitution databases with biochemical characterizations of function are needed. We report here a database derived from mutational, biochemical, bioinformatic, structural, pathological and computational studies of a highly studied protein family-pyruvate kinase (PYK). A centerpiece of this database is the biochemical characterization-including quantitative evaluation of allosteric regulation-of the changes that accompany substitutions at positions that sample the full conservation range observed in the PYK family. We have used these data to facilitate critical advances in the foundational studies of allosteric regulation and protein evolution and as rigorous benchmarks for testing protein predictions. We trust that the collected dataset will be useful for the broader scientific community in the further development of prediction algorithms. Database URL https://github.com/djparente/PYK-DB.

Reversible domain closure modulates GlnBP ligand binding affinity

Glutamine binding protein (GlnBP) is an Escherichia Coli periplasmic binding protein, which binds and carries glutamine to the inner membrane ATP-binding cassette (ABC) transporter. GlnBP binds the ligand with affinity around 0.1μM measured by isothermal titration calorimetry (ITC) and ligand binding stabilizes protein structure shown by its increase in thermodynamic stability. However, the molecular determinant of GlnBP ligand binding is not known. Electrostatic and hydrophobic interaction between GlnBP and glutamine are critical factors. We propose that the freedome of closure movement is also vital for ligand binding. In order to approve this hypothesis, we generate a series of mutants with different linker length that has different magnitude of domain closure. Mutants show different ligand binding affinity, which indicates that the propensity of domain closure determines the ligand binding affinity. Ligand binding triggers gradual ensemble conformational change. Structural changes upon ligand binding are monitored by combination of small angle x-ray scattering (SAXS) and NMR spectroscopy. Detailed structure characterization of GlnBP contributes to a better understanding of ligand binding and provides the structural basis for biosensor design.

H/D Exchange Characterization of Silent Coupling: Entropy-Enthalpy Compensation in Allostery

The allosteric coupling constant in K-type allosteric systems is defined as a ratio of the binding of substrate in the absence of effector to the binding of the substrate in the presence of a saturating concentration of effector. As a result, the coupling constant is itself an equilibrium value comprised of a ΔH and a TΔS component. In the scenario in which TΔS completely compensates ΔH, no allosteric influence of effector binding on substrate affinity is observed. However, in this "silent coupling" scenario, the presence of effector causes a change in the ΔH associated with substrate binding. A suggestion has now been made that "silent modulators" are ideal drug leads because they can be modified to act as either allosteric activators or inhibitors. Any attempt to rationally design the effector to be an allosteric activator or inhibitor is likely to be benefitted by knowledge of the mechanism that gives rise to coupling. Hydrogen/deuterium exchange with mass spectrometry detection has now been used to identify regions of proteins that experience conformational and/or dynamic changes in the allosteric regulation. Here, we demonstrate the expected temperature dependence of the allosteric regulation of rabbit muscle pyruvate kinase by Ala to demonstrate that this effector reduces substrate (phosphoenolpyruvate) affinity at 35°C and at 10°C but is silent at intermediate temperatures. We then explore the use of hydrogen/deuterium exchange with mass spectrometry to evaluate the areas of the protein that are modified in the mechanism that gives rise to the silent coupling between Ala and phosphoenolpyruvate. Many of the peptide regions of the protein identified as changing in this silent system (Ala as the effector) were included in changes previously identified for allosteric inhibition by Phe.

What Mutagenesis Can and Cannot Reveal About Allostery

Allosteric regulation of protein function is recognized to be widespread throughout biology; however, knowledge of allosteric mechanisms, the molecular changes within a protein that couple one binding site to another, is limited. Although mutagenesis is often used to probe allosteric mechanisms, we consider herein what the outcome of a mutagenesis study truly reveals about an allosteric mechanism. Arguably, the best way to evaluate the effects of a mutation on allostery is to monitor the allosteric coupling constant (Qax), a ratio of the substrate binding constants in the absence versus presence of an allosteric effector. A range of substitutions at a given residue position in a protein can reveal when a particular substitution causes gain-of-function, which addresses a key challenge in interpreting mutation-dependent changes in the magnitude of Qax. Thus, whole-protein mutagenesis studies offer an acceptable means of identifying residues that contribute to an allosteric mechanism. With this focus on monitoring Qax, and keeping in mind the equilibrium nature of allostery, we consider alternative possibilities for what an allosteric mechanism might be. We conclude that different possible mechanisms (rotation-of-solid-domains, movement of secondary structure, side-chain repacking, changes in dynamics, etc.) will result in different findings in whole-protein mutagenesis studies.

Exploring the limits of the usefulness of mutagenesis in studies of allosteric mechanisms

The outcome of structure-guided mutational analyses is often used in support of postulated mechanisms of protein allostery. However, the limits of how informative mutations can be in understanding allosteric mechanisms are not completely clear. Here, we report an exercise to evaluate whether mutational data can support a simplistic mechanistic model, developed with minimal data inputs. Due to the lack of a mechanism to explain how alanine allosterically modifies the affinity of human liver pyruvate kinase (approved symbol PKLR) for its substrate, phosphoenolpyruvate, we proposed a speculative allosteric mechanism for this system. Within the allosteric amino-acid-binding site (something in the effector site must, of necessity, contribute to the allosteric mechanism), we implemented multiple mutational strategies: (1) site-directed random mutagenesis at positions that contact bound alanine and (2) mutations to probe specific questions. Despite acknowledged inadequacies used to formulate the speculative mechanism, many mutations modified the allosteric coupling constant (Q_ax ) consistent with that mechanism. The observed support for this speculative mechanism leaves us to ponder the best use of mutational data in structure-function studies of allosteric mechanisms. The mutational databank derived from this exercise has an independent value for training and testing algorithms specific to allostery.

Conformational Dynamics and Allostery in Pyruvate Kinase

Pyruvate kinase catalyzes the final step in glycolysis and is allosterically regulated to control flux through the pathway. Two models are proposed to explain how Escherichia coli pyruvate kinase type 1 is allosterically regulated: the "domain rotation model" suggests that both the domains within the monomer and the monomers within the tetramer reorient with respect to one another; the "rigid body reorientation model" proposes only a reorientation of the monomers within the tetramer causing rigidification of the active site. To test these hypotheses and elucidate the conformational and dynamic changes that drive allostery, we performed time-resolved electrospray ionization mass spectrometry coupled to hydrogen-deuterium exchange studies followed by mutagenic analysis to test the activation mechanism. Global exchange experiments, supported by thermostability studies, demonstrate that fructose 1,6-bisphosphate binding to the allosteric domain causes a shift toward a globally more dynamic ensemble of conformations. Mapping deuterium exchange to peptides within the enzyme highlight site-specific regions with altered conformational dynamics, many of which increase in conformational flexibility. Based upon these and mutagenic studies, we propose an allosteric mechanism whereby the binding of fructose 1,6-bisphosphate destabilizes an α-helix that bridges the allosteric and active site domains within the monomeric unit. This destabilizes the β-strands within the (β/α)8-barrel domain and the linked active site loops that are responsible for substrate binding. Our data are consistent with the domain rotation model but inconsistent with the rigid body reorientation model given the increased flexibility at the interdomain interface, and we can for the first time explain how fructose 1,6-bisphosphate affects the active site.

Accurate SAXS profile computation and its assessment by contrast variation experiments

A major challenge in structural biology is to characterize structures of proteins and their assemblies in solution. At low resolution, such a characterization may be achieved by small angle x-ray scattering (SAXS). Because SAXS analyses often require comparing profiles calculated from many atomic models against those determined by experiment, rapid and accurate profile computation from molecular structures is needed. We developed fast open-source x-ray scattering (FoXS) for profile computation. To match the experimental profile within the experimental noise, FoXS explicitly computes all interatomic distances and implicitly models the first hydration layer of the molecule. For assessing the accuracy of the modeled hydration layer, we performed contrast variation experiments for glucose isomerase and lysozyme, and found that FoXS can accurately represent density changes of this layer. The hydration layer model was also compared with a SAXS profile calculated for the explicit water molecules in the high-resolution structures of glucose isomerase and lysozyme. We tested FoXS on eleven protein, one DNA, and two RNA structures, revealing superior accuracy and speed versus CRYSOL, AquaSAXS, the Zernike polynomials-based method, and Fast-SAXS-pro. In addition, we demonstrated a significant correlation of the SAXS score with the accuracy of a structural model. Moreover, FoXS utility for analyzing heterogeneous samples was demonstrated for intrinsically flexible XLF-XRCC4 filaments and Ligase III-DNA complex. FoXS is extensively used as a standalone web server as a component of integrative structure determination by programs IMP, Chimera, and BILBOMD, as well as in other applications that require rapidly and accurately calculated SAXS profiles.

A single amino-acid substitution allows endo-polygalacturonase of Fusarium verticillioides to acquire recognition by PGIP2 from Phaseolus vulgaris

Polygalacturonases (PGs) are secreted by phytopathogenic fungi to degrade the plant cell wall homogalacturonan during plant infection. To counteract Pgs, plants have evolved polygalacturonase-inhibiting proteins (PGIPs) that slow down fungal infection and defend cell wall integrity. PGIPs favour the accumulation of oligogalacturonides, which are homogalacturonan fragments that act as endogenous elicitors of plant defence responses. We have previously shown that PGIP2 from Phaseolus vulgaris (PvPGIP2) forms a complex with PG from Fusarium phyllophilum (FpPG), hindering the enzyme active site cleft from substrate. Here we analyse by small angle X-ray scattering (SAXS) the interaction between PvPGIP2 and a PG from Colletotrichum lupini (CluPG1). We show a different shape of the PG-PGIP complex, which allows substrate entry and provides a structural explanation for the different inhibition kinetics exhibited by PvPGIP2 towards the two isoenzymes. The analysis of SAXS structures allowed us to investigate the basis of the inability of PG from Fusarium verticilloides (FvPG) to be inhibited by PvPGIP2 or by any other known PGIP. FvPG is 92.5% identical to FpPG, and we show here, by both loss- and gain-of-function mutations, that a single amino acid site acts as a switch for FvPG recognition by PvPGIP2.

A single amino-acid substitution allows endo-polygalacturonase of Fusarium verticillioides to acquire recognition by PGIP2 from Phaseolus vulgaris

Polygalacturonases (PGs) are secreted by phytopathogenic fungi to degrade the plant cell wall homogalacturonan during plant infection. To counteract Pgs, plants have evolved polygalacturonase-inhibiting proteins (PGIPs) that slow down fungal infection and defend cell wall integrity. PGIPs favour the accumulation of oligogalacturonides, which are homogalacturonan fragments that act as endogenous elicitors of plant defence responses. We have previously shown that PGIP2 from Phaseolus vulgaris (PvPGIP2) forms a complex with PG from Fusarium phyllophilum (FpPG), hindering the enzyme active site cleft from substrate. Here we analyse by small angle X-ray scattering (SAXS) the interaction between PvPGIP2 and a PG from Colletotrichum lupini (CluPG1). We show a different shape of the PG-PGIP complex, which allows substrate entry and provides a structural explanation for the different inhibition kinetics exhibited by PvPGIP2 towards the two isoenzymes. The analysis of SAXS structures allowed us to investigate the basis of the inability of PG from Fusarium verticilloides (FvPG) to be inhibited by PvPGIP2 or by any other known PGIP. FvPG is 92.5% identical to FpPG, and we show here, by both loss- and gain-of-function mutations, that a single amino acid site acts as a switch for FvPG recognition by PvPGIP2.

Identification of regions of rabbit muscle pyruvate kinase important for allosteric regulation by phenylalanine, detected by H/D exchange mass spectrometry

Mass spectrometry has been used to determine the number of exchangeable backbone amide protons and the associated rate constants that are altered when rabbit muscle pyruvate kinase (rM1-PYK) binds either the allosteric inhibitor (phenylalanine) or a nonallosteric analogue of the inhibitor. Alanine is used as the nonallosteric analogue because it binds competitively with phenylalanine but elicits a negligible allosteric inhibition, i.e., a negligible reduction in the affinity of rM1-PYK for the substrate, phosphoenolpyruvate. This experimental design is expected to distinguish changes in the protein caused by effector binding (i.e., those changes common upon the addition of alanine vs phenylalanine) from changes associated with allosteric regulation (i.e., those elicited by the addition of phenylalanine binding, but not alanine binding). High-quality peptic fragments covering 98% of the protein were identified. Changes in both the number of exchangeable protons per peptide and in the rate constant associated with exchange highlight regions of the protein with allosteric roles. The set of allosterically relevant peptides identified by this technique includes residues previously identified by mutagenesis to have roles in allosteric regulation by phenylalanine.

Integrative structural modeling with small angle X-ray scattering profiles

Recent technological advances enabled high-throughput collection of Small Angle X-ray Scattering (SAXS) profiles of biological macromolecules. Thus, computational methods for integrating SAXS profiles into structural modeling are needed more than ever. Here, we review specifically the use of SAXS profiles for the structural modeling of proteins, nucleic acids, and their complexes. First, the approaches for computing theoretical SAXS profiles from structures are presented. Second, computational methods for predicting protein structures, dynamics of proteins in solution, and assembly structures are covered. Third, we discuss the use of SAXS profiles in integrative structure modeling approaches that depend simultaneously on several data types.

Monitoring allostery in D2O: a necessary control in studies using hydrogen/deuterium exchange to characterize allosteric regulation

There is currently a renewed focus aimed at understanding allosteric mechanisms at atomic resolution. This current interest seeks to understand how both changes in protein conformations and changes in protein dynamics contribute to relaying an allosteric signal between two ligand binding sites on a protein (e.g., active and allosteric sites). Both nuclear magnetic resonance (NMR), by monitoring protein dynamics directly, and hydrogen/deuterium exchange, by monitoring solvent accessibility of backbone amides, offer insights into protein dynamics. Unfortunately, many allosteric proteins exceed the size limitations of standard NMR techniques. Although hydrogen/deuterium exchange as detected by mass spectrometry (H/DX-MS) offers an alternative evaluation method, any application of hydrogen/deuterium exchange requires that the property being measured functions in both H(2)O and D(2)O. Due to the promising future H/DX-MS has in the evaluation of allosteric mechanisms in large proteins, we demonstrate an evaluation of allosteric regulation in D(2)O. Exemplified using phenylalanine inhibition of rabbit muscle pyruvate kinase, we find that binding of the inhibitor is greatly reduced in D(2)O, but the effector continues to elicit an allosteric response.

Structural and functional energetic linkages in allosteric regulation of muscle pyruvate kinase

The understanding of the molecular mechanisms of allostery in rabbit muscle pyruvate kinase (RMPK) is still in its infancy. Although, there is a paucity of knowledge on the ground rules on how its functions are regulated, RMPK is an ideal system to address basic questions regarding the fundamental chemical principles governing the regulatory mechanisms about this enzyme which has a TIM (α/β)(8) barrel structural motif [Copley, R. R., and Bork, P. (2000). Homology among (βα)8 barrels: Implications for the evolution of metabolic pathways. J. Mol. Biol.303, 627-640; Farber, G. K., and Petsko, G. A. (1990). The evolution of α/ß barrel enzymes. Trends Biochem.15, 228-234; Gerlt, J. A., and Babbitt, P. C. (2001). Divergent evolution of enzymatic function: Mechanistically diverse superfamilies and functionally distinct superfamilies. Annu. Rev. Biochem.70, 209-246; Heggi, H., and Gerstein, M. (1999). The relationship between protein structure and function: A comprehensive survey with application to the yeast genome. J. Mol. Biol.288, 147-164; Wierenga, R. K. (2001). The TIM-barrel fold: A versatile framework for efficient enzymes. FEB Lett.492, 193-198]. RMPK is a homotetramer. Each subunit consists of 530 amino acids and multiple domains. The active site resides between the A and B domains. Besides the basic TIM-barrel motif, RMPK also exhibits looped-out regions in the α/β barrel of each monomer forming the B- and C-domains. The two isozymes of PK, namely, the kidney and muscle isozymes, exhibit very different allosteric behaviors under the same experimental condition. The only amino acid sequence differences between the mammalian kidney and muscle PK isozymes are located in the C-domain and are involved in intersubunit interactions. Thus, embedded in these two isozymes of PK are the rules involved in engineering the popular TIM (α/β)(8) motif to modulate its allosteric properties. The PK system exhibits a lot of the properties that will allow mining of the ground rules governing the correlative linkages between sequence-fold-function. In this chapter, we review the approaches to acquire the fundamental functional and structural energetics that establish the linkages among this intricate network of linked multiequilibria. Results from these diverse approaches are integrated to establish a working model to represent the complex network of multiple linked reactions which ultimately leads to the observation of allosteric regulation of PK.