Substrate-modulated thermal fluctuations affect long-range allosteric signaling in protein homodimers: exemplified in CAP

cAMP is an allosteric modulator of DNA-binding specificity in the cAMP receptor protein from Mycobacterium tuberculosis

Allosteric proteins with multiple subunits and ligand-binding sites are central in regulating biological signals. The cAMP receptor protein from Mycobacterium tuberculosis (CRP_MTB) is a global regulator of transcription composed of two identical subunits, each one harboring structurally conserved cAMP- and DNA-binding sites. The mechanisms by which these four binding sites are allosterically coupled in CRP_MTB remain unclear. Here, we investigate the binding mechanism between CRP_MTB and cAMP, and the linkage between cAMP and DNA interactions. Using calorimetric and fluorescence-based assays, we find that cAMP binding is entropically driven and displays negative cooperativity. Fluorescence anisotropy experiments show that apo-CRP_MTB forms high-order CRP_MTB–DNA oligomers through interactions with nonspecific DNA sequences or preformed CRP_MTB–DNA complexes. Moreover, we find that cAMP prevents and reverses the formation of CRP_MTB–DNA oligomers, reduces the affinity of CRP_MTB for nonspecific DNA sequences, and stabilizes a 1-to-1 CRP_MTB–DNA complex, but does not increase the affinity for DNA like in the canonical CRP from Escherichia coli (CRP_Ecoli). DNA-binding assays as a function of cAMP concentration indicate that one cAMP molecule per homodimer dissociates high-order CRP_MTB–DNA oligomers into 1-to-1 complexes. These cAMP-mediated allosteric effects are lost in the double-mutant L47P/E178K found in CRP from Mycobacterium bovis Bacille Calmette-Guérin (CRP_BCG). The functional behavior, thermodynamic stability, and dimerization constant of CRP_BCG are not due to additive effects of L47P and E178K, indicating long-range interactions between these two sites. Altogether, we provide a previously undescribed archetype of cAMP-mediated allosteric regulation that differs from CRP_Ecoli, illustrating that structural homology does not imply allosteric homology.

Allosteric Response of DNA Recognition Helices of Catabolite Activator Protein to cAMP and DNA Binding

The homodimeric catabolite activator protein (CAP) regulates the transcription of several bacterial genes based on the cellular concentration of cyclic adenosine monophosphate (cAMP). The binding of cAMP to CAP triggers allosteric communication between the cAMP binding domains (CBD) and DNA binding domains (DBD) of CAP, which entails repositioning of DNA recognition helices (F-helices) in the DBD to dock favorably to the target DNA. Despite considerable progress, much remains to be understood about the mechanistic details of DNA recognition by CAP and about the map of allosteric pathways involved in CAP-mediated gene transcription. The present study uses molecular dynamics and umbrella sampling simulations to investigate the mechanism of cAMP- and DNA-induced changes in the conformation and energetics of F-helices observed during the allosteric regulation of CAP by cAMP and the subsequent binding to the DNA promoter region. Using novel collective variables, the free energy profiles associated with the orientation and dynamics of F-helices in the unliganded, cAMP-bound, and cAMP-DNA-bound states of CAP are calculated and compared. The binding-induced alterations in the resultant free energy profiles reveal important flexibility constraints imposed on DBD upon cAMP and DNA binding. A comprehensive analysis of residue-wise interaction maps reveals potential allosteric pathways between CBD and DBD that facilitate the allosteric transduction of regulatory signals in CAP. The revelation that the predicted allosteric pathways crisscross the intersubunit interface offers important clues on the microscopic origin of the intersubunit cooperativity and dimer stability of CAP.

Allostery in Its Many Disguises: From Theory to Applications

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

Ligand-regulated oligomerisation of allosterically interacting proteins

The binding of ligands to distinct sites at proteins or at protein clusters is often cooperative or anti-cooperative due to allosteric signalling between those sites. The allostery is usually attributed to a configurational change of the proteins from a relaxed to a configurationally different tense state. Alternatively, as originally proposed by Cooper and Dryden, a tense state may be achieved by merely restricting the thermal vibrations of the protein around its mean configuration. In this work, we provide theoretical tools to investigate fluctuation allostery using cooling and titration experiments in which ligands regulate dimerisation, or ring or chain formation. We discuss in detail how ligands may regulate the supramolecular (co)polymerisation of liganded and unliganded proteins.

The 'allosteron' model for entropic allostery of self-assembly

Using the simple 'allosteron' model, we show that it is possible, in principle, to elicit pathways by which fluctuation allostery affects self-assembly of protein complexes. We treat the cases of (i) protein fibrils and nucleation, (ii) n-mer protein complexes, and (iii) weakly attractive allosteric interactions in protein-like soft nanoscale objects that can be tuned to define exclusive self-associating families.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Dynamic Transmission of Protein Allostery without Structural Change: Spatial Pathways or Global Modes?

We examine the contrast between mechanisms for allosteric signaling that involve structural change, and those that do not, from the perspective of allosteric pathways. In particular we treat in detail the case of fluctuation-allostery by which amplitude modulation of the thermal fluctuations of the elastic normal modes conveys the allosteric signal, and address the question of what an allosteric pathway means in this case. We find that a perturbation theory of thermal elastic solids and nonperturbative approach (by super-coarse-graining elasticity into internal bending modes) have opposite signatures in their structure of correlated pathways. We illustrate the effect from analysis of previous results from GlxR of Corynebacterium glutamicum, an example of the CRP/FNR transcription family of allosteric homodimers. We find that the visibility of both correlated pathways and disconnected sites of correlated motion in this protein suggests that mechanisms of local elastic stretch and bend are recruited for the purpose of creating and controlling allosteric cooperativity.

The Role of Protein-Ligand Contacts in Allosteric Regulation of the Escherichia coli Catabolite Activator Protein

Allostery is a fundamental process by which ligand binding to a protein alters its activity at a distant site. Both experimental and theoretical evidence demonstrate that allostery can be communicated through altered slow relaxation protein dynamics without conformational change. The catabolite activator protein (CAP) of Escherichia coli is an exemplar for the analysis of such entropically driven allostery. Negative allostery in CAP occurs between identical cAMP binding sites. Changes to the cAMP-binding pocket can therefore impact the allosteric properties of CAP. Here we demonstrate, through a combination of coarse-grained modeling, isothermal calorimetry, and structural analysis, that decreasing the affinity of CAP for cAMP enhances negative cooperativity through an entropic penalty for ligand binding. The use of variant cAMP ligands indicates the data are not explained by structural heterogeneity between protein mutants. We observe computationally that altered interaction strength between CAP and cAMP variously modifies the change in allosteric cooperativity due to second site CAP mutations. As the degree of correlated motion between the cAMP-contacting site and a second site on CAP increases, there is a tendency for computed double mutations at these sites to drive CAP toward noncooperativity. Naturally occurring pairs of covarying residues in CAP do not display this tendency, suggesting a selection pressure to fine tune allostery on changes to the CAP ligand-binding pocket without a drive to a noncooperative state. In general, we hypothesize an evolutionary selection pressure to retain slow relaxation dynamics-induced allostery in proteins in which evolution of the ligand-binding site is occurring.

Terahertz optical measurements of correlated motions with possible allosteric function

A suggested mechanism for allosteric response is the distortion of the energy landscape with agonist binding changing the protein structure's access to functional configurations. Intramolecular vibrations are indicative of the energy landscape and may have trajectories that enable functional conformational change. Here, we discuss the development of an optical method to measure the intramolecular vibrations in proteins, namely, crystal anisotropy terahertz microscopy, and the various approaches which can be used to identify the spectral data with specific structural motions.

Global low-frequency motions in protein allostery: CAP as a model system

Allostery is a fundamental process by which ligand binding to a protein alters its activity at a distant site. There is considerable evidence that allosteric cooperativity can be communicated by the modulation of protein dynamics without conformational change. The Catabolite Activator Protein (CAP) of Escherichia coli is an important experimental exemplar for entropically driven allostery. Here we discuss recent experimentally supported theoretical analysis that highlights the role of global low-frequency dynamics in allostery in CAP and identify how allostery arises as a natural consequence of changes in global low-frequency protein fluctuations on ligand binding.

The crystal structures of apo and cAMP-bound GlxR from Corynebacterium glutamicum reveal structural and dynamic changes upon cAMP binding in CRP/FNR family transcription factors

The cyclic AMP-dependent transcriptional regulator GlxR from Corynebacterium glutamicum is a member of the super-family of CRP/FNR (cyclic AMP receptor protein/fumarate and nitrate reduction regulator) transcriptional regulators that play central roles in bacterial metabolic regulatory networks. In C. glutamicum, which is widely used for the industrial production of amino acids and serves as a non-pathogenic model organism for members of the Corynebacteriales including Mycobacterium tuberculosis, the GlxR homodimer controls the transcription of a large number of genes involved in carbon metabolism. GlxR therefore represents a key target for understanding the regulation and coordination of C. glutamicum metabolism. Here we investigate cylic AMP and DNA binding of GlxR from C. glutamicum and describe the crystal structures of apo GlxR determined at a resolution of 2.5 Å, and two crystal forms of holo GlxR at resolutions of 2.38 and 1.82 Å, respectively. The detailed structural analysis and comparison of GlxR with CRP reveals that the protein undergoes a distinctive conformational change upon cyclic AMP binding leading to a dimer structure more compatible to DNA-binding. As the two binding sites in the GlxR homodimer are structurally identical dynamic changes upon binding of the first ligand are responsible for the allosteric behavior. The results presented here show how dynamic and structural changes in GlxR lead to optimization of orientation and distance of its two DNA-binding helices for optimal DNA recognition.

Allostery without conformation change: modelling protein dynamics at multiple scales

The original ideas of Cooper and Dryden, that allosteric signalling can be induced between distant binding sites on proteins without any change in mean structural conformation, has proved to be a remarkably prescient insight into the rich structure of protein dynamics. It represents an alternative to the celebrated Monod-Wyman-Changeux mechanism and proposes that modulation of the amplitude of thermal fluctuations around a mean structure, rather than shifts in the structure itself, give rise to allostery in ligand binding. In a complementary approach to experiments on real proteins, here we take a theoretical route to identify the necessary structural components of this mechanism. By reviewing and extending an approach that moves from very coarse-grained to more detailed models, we show that, a fundamental requirement for a body supporting fluctuation-induced allostery is a strongly inhomogeneous elastic modulus. This requirement is reflected in many real proteins, where a good approximation of the elastic structure maps strongly coherent domains onto rigid blocks connected by more flexible interface regions.

Modulation of global low-frequency motions underlies allosteric regulation: demonstration in CRP/FNR family transcription factors

Allostery in bacterial transcription factors arises from changes in global low-frequency protein dynamics. Amino acids that regulate low-frequency dynamics are identified and seen to be evolutionarily conserved.

Allostery is a fundamental process by which ligand binding to a protein alters its activity at a distinct site. There is growing evidence that allosteric cooperativity can be communicated by modulation of protein dynamics without conformational change. The mechanisms, however, for communicating dynamic fluctuations between sites are debated. We provide a foundational theory for how allostery can occur as a function of low-frequency dynamics without a change in structure. We have generated coarse-grained models that describe the protein backbone motions of the CRP/FNR family transcription factors, CAP of Escherichia coli and GlxR of Corynebacterium glutamicum. The latter we demonstrate as a new exemplar for allostery without conformation change. We observe that binding the first molecule of cAMP ligand is correlated with modulation of the global normal modes and negative cooperativity for binding the second cAMP ligand without a change in mean structure. The theory makes key experimental predictions that are tested through an analysis of variant proteins by structural biology and isothermal calorimetry. Quantifying allostery as a free energy landscape revealed a protein “design space” that identified the inter- and intramolecular regulatory parameters that frame CRP/FNR family allostery. Furthermore, through analyzing CAP variants from diverse species, we demonstrate an evolutionary selection pressure to conserve residues crucial for allosteric control. This finding provides a link between the position of CRP/FNR transcription factors within the allosteric free energy landscapes and evolutionary selection pressures. Our study therefore reveals significant features of the mechanistic basis for allostery. Changes in low-frequency dynamics correlate with allosteric effects on ligand binding without the requirement for a defined spatial pathway. In addition to evolving suitable three-dimensional structures, CRP/FNR family transcription factors have been selected to occupy a dynamic space that fine-tunes biological activity and thus establishes the means to engineer allosteric mechanisms driven by low-frequency dynamics.

Allostery is a process by which a molecule binding to one site of a protein alters the activity of the protein at another site. Allostery is typically thought to occur through a change in protein structure, but there is now clear evidence that the dynamic properties of a protein can also regulate allostery without a change in overall conformation. Here we examine two members of a large family of bacterial transcription factors and provide a mechanism to describe the allosteric binding of their activating ligands. We demonstrate, in these systems, that allostery arises as a natural consequence of changes in global low-frequency protein fluctuations on ligand binding. We further demonstrate that the higher dimensional parameter space that describes all potential variant transcription factors can be reduced to a two-dimensional free energy landscape that determines the key molecular parameters that predominantly regulate allostery. We additionally show that the amino acids we determine as contributing sensitively to allosteric control tend to be conserved in diverse bacteria; thus we identify a link between residues that contribute to low-frequency fluctuations and evolutionary selection pressures.

Thermodynamics of allostery paves a way to allosteric drugs

We overview here our recent work on the thermodynamic view of allosteric regulation and communication. Starting from the geometry-based prediction of regulatory binding sites in a static structure, we move on to exploring a connection between ligand binding and the intrinsic dynamics of the protein molecule. We describe here two recently introduced measures, binding leverage and leverage coupling, which allow one to analyze the molecular basis of allosteric regulation. We discuss the advantages of these measures and show that they work universally in proteins of different sizes, oligomeric states, and functions. We also point the problems that have to be solved before completing an atomic level description of allostery, and briefly discuss ideas for computational design of allosteric drugs. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

Allosteric inhibition of a zinc-sensing transcriptional repressor: insights into the arsenic repressor (ArsR) family

The molecular basis of allosteric regulation remains a subject of intense interest. Staphylococcus aureus CzrA is a member of the ubiquitous arsenic repressor (ArsR) family of bacterial homodimeric metal-sensing proteins and has emerged as a model system for understanding allosteric regulation of operator DNA binding by transition metal ions. Using unnatural amino acid substitution and a standard linkage analysis, we show that a His97' NH(ε2)...O=C His67 quaternary structural hydrogen bond is an energetically significant contributor to the magnitude of the allosteric coupling free energy, ∆Gc. A "cavity" introduced just beneath this hydrogen bond in V66A/L68V CzrA results in a significant reduction in regulation by Zn(II) despite adopting a wild-type global structure and Zn(II) binding and DNA binding affinities only minimally affected from wild type. The energetics of Zn(II) binding and heterotropic coupling free energies (∆Hc, -T∆Sc) of the double mutant are also radically altered and suggest that increased internal dynamics leads to poorer allosteric negative regulation in V66A/L68V CzrA. A statistical coupling analysis of 3000 ArsR proteins reveals a sector that links the DNA-binding determinants and the α5 Zn(II)-sensing sites through V66/L68 in CzrA. We propose that distinct regulatory sites uniquely characteristic of individual ArsR proteins result from evolution of distinct connectivities to this sector, each capable of driving the same biological outcome, transcriptional derepression.

Allosteric activation via kinetic control: potassium accelerates a conformational change in IMP dehydrogenase

Allosteric activators are generally believed to shift the equilibrium distribution of enzyme conformations to favor a catalytically productive structure; the kinetics of conformational exchange is seldom addressed. Several observations suggested that the usual allosteric mechanism might not apply to the activation of IMP dehydrogenase (IMPDH) by monovalent cations. Therefore, we investigated the mechanism of K(+) activation in IMPDH by delineating the kinetic mechanism in the absence of monovalent cations. Surprisingly, the K(+) dependence of k(cat) derives from the rate of flap closure, which increases by ≥65-fold in the presence of K(+). We performed both alchemical free energy simulations and potential of mean force calculations using the orthogonal space random walk strategy to computationally analyze how K(+) accelerates this conformational change. The simulations recapitulate the preference of IMPDH for K(+), validating the computational models. When K(+) is replaced with a dummy ion, the residues of the K(+) binding site relax into ordered secondary structure, creating a barrier to conformational exchange. K(+) mobilizes these residues by providing alternate interactions for the main chain carbonyls. Potential of mean force calculations indicate that K(+) changes the shape of the energy well, shrinking the reaction coordinate by shifting the closed conformation toward the open state. This work suggests that allosteric regulation can be under kinetic as well as thermodynamic control.

Binding leverage as a molecular basis for allosteric regulation

Allosteric regulation involves conformational transitions or fluctuations between a few closely related states, caused by the binding of effector molecules. We introduce a quantity called binding leverage that measures the ability of a binding site to couple to the intrinsic motions of a protein. We use Monte Carlo simulations to generate potential binding sites and either normal modes or pairs of crystal structures to describe relevant motions. We analyze single catalytic domains and multimeric allosteric enzymes with complex regulation. For the majority of the analyzed proteins, we find that both catalytic and allosteric sites have high binding leverage. Furthermore, our analysis of the catabolite activator protein, which is allosteric without conformational change, shows that its regulation involves other types of motion than those modulated at sites with high binding leverage. Our results point to the importance of incorporating dynamic information when predicting functional sites. Because it is possible to calculate binding leverage from a single crystal structure it can be used for characterizing proteins of unknown function and predicting latent allosteric sites in any protein, with implications for drug design.

Allosteric protein regulation is the mechanism by which binding of a molecule to one site in a protein affects the activity at another site. Although the two classical phenomenological models, Monod-Wyman-Changeux (MWC) and Koshland-Némethy-Filmer (KNF), span from the case of hemoglobin to membrane receptors, they do not describe the intramolecular interactions involved. The coupling between two allosterically connected sites commonly takes place through coherent collective motion involving the whole protein. We therefore introduce a quantity called binding leverage to measure the strength of the coupling between particular binding sites and such motions. We show that high binding leverage is a characteristic of both allosteric sites and catalytic sites, emphasizing that both enzymatic function and allosteric regulation require a coupling between ligand binding and protein dynamics. We also consider the first known case of purely entropic allostery, where ligand binding only affects the amplitudes of fluctuations. We find that the binding site in this protein does not primarily connect to collective motions – instead the modulation of fluctuations is controlled from a deeply buried and highly connected site. Finally, sites with high binding leverage but no known biological function could be latent allosteric sites, and thus drug targets.

Role of fluctuations in ligand binding cooperativity of membrane receptors

Signal transduction upon binding of a ligand to a membrane protein can occur not only via allosteric conformational changes but also through fluctuations. We report a numerical study on the influence of conformational fluctuations on the cooperativity of a binding reaction in a simple model of an integral membrane receptor consisting of transmembrane helices. We find that small fluctuations lateral as well as perpendicular to the membrane can increase the cooperativity, with the former more dominant. Moreover, too much fluctuation induces negative cooperativity. Proteins with fewer than four helices do not show positive cooperativity under any circumstances. This behavior is rather robust, and independent of the receptor topology or ligand size. Fluctuations measured in all-atom molecular dynamics simulations of a G-protein coupled receptor fall within the predicted region of maximum cooperativity.

Positions 94-98 of the lactose repressor N-subdomain monomer-monomer interface are critical for allosteric communication

The central region of the LacI N-subdomain monomer-monomer interface includes residues K84, V94, V95, V96, S97, and M98. The side chains of these residues line the β-strands at this interface and interact to create a network of hydrophobic, charged, and polar interactions that significantly rearranges in different functional states of LacI. Prior work showed that converting K84 to an apolar residue or converting V96 to an acidic residue impedes the allosteric response to inducer. Thus, we postulated that a disproportionate number of substitutions in this region of the monomer-monomer interface would alter the complex features of the LacI allosteric response. To explore this hypothesis, acidic, basic, polar, and apolar mutations were introduced at positions 94-98. Despite their varied locations along the β-strands that flank the interface, ∼70% of the mutations impact allosteric behavior, with the most significant effects found for charged substitutions. Of note, many of the LacI variants with minor functional impact exhibited altered stability to urea denaturation. The results confirm the critical role of amino acids 94-98 and indicate that this N-subdomain interface forms a primary pathway in LacI allosteric response.