Dynamically driven protein allostery

Sending signals dynamically

Proteins mediate transmission of signals along intercellular and intracellular pathways and between the exterior and the interior of a cell. The dynamic properties of signaling proteins are crucial to their functions. We discuss emerging paradigms for the role of protein dynamics in signaling. A central tenet is that proteins fluctuate among many states on evolutionarily selected energy landscapes. Upstream signals remodel this landscape, causing signaling proteins to transmit information to downstream partners. New methods provide insight into the dynamic properties of signaling proteins at the atomic scale. The next stages in the signaling hierarchy-how multiple signals are integrated and how cellular signaling pathways are organized in space and time-present exciting challenges for the future, requiring bold multidisciplinary approaches.

Structural basis for cAMP-mediated allosteric control of the catabolite activator protein

The cAMP-mediated allosteric transition in the catabolite activator protein (CAP; also known as the cAMP receptor protein, CRP) is a textbook example of modulation of DNA-binding activity by small-molecule binding. Here we report the structure of CAP in the absence of cAMP, which, together with structures of CAP in the presence of cAMP, defines atomic details of the cAMP-mediated allosteric transition. The structural changes, and their relationship to cAMP binding and DNA binding, are remarkably clear and simple. Binding of cAMP results in a coil-to-helix transition that extends the coiled-coil dimerization interface of CAP by 3 turns of helix and concomitantly causes rotation, by approximately 60 degrees , and translation, by approximately 7 A, of the DNA-binding domains (DBDs) of CAP, positioning the recognition helices in the DBDs in the correct orientation to interact with DNA. The allosteric transition is stabilized further by expulsion of an aromatic residue from the cAMP-binding pocket upon cAMP binding. The results define the structural mechanisms that underlie allosteric control of this prototypic transcriptional regulatory factor and provide an illustrative example of how effector-mediated structural changes can control the activity of regulatory proteins.

The structural basis of allosteric regulation in proteins

Allosteric regulation of protein function occurs when the regulatory trigger, such as the binding of a small-molecule effector or inhibitor, takes place some distance from the protein's, or protein complex's, active site. This distance can be a few A, or tens of A. Many proteins are regulated in this way and exhibit a variety of allosteric mechanisms. Here we review how analyses of experimentally determined models of protein 3D structures, using either X-ray crystallography or NMR spectroscopy, have revealed some of the mechanisms involved.

Mechanical network in titin immunoglobulin from force distribution analysis

The role of mechanical force in cellular processes is increasingly revealed by single molecule experiments and simulations of force-induced transitions in proteins. How the applied force propagates within proteins determines their mechanical behavior yet remains largely unknown. We present a new method based on molecular dynamics simulations to disclose the distribution of strain in protein structures, here for the newly determined high-resolution crystal structure of I27, a titin immunoglobulin (IG) domain. We obtain a sparse, spatially connected, and highly anisotropic mechanical network. This allows us to detect load-bearing motifs composed of interstrand hydrogen bonds and hydrophobic core interactions, including parts distal to the site to which force was applied. The role of the force distribution pattern for mechanical stability is tested by in silico unfolding of I27 mutants. We then compare the observed force pattern to the sparse network of coevolved residues found in this family. We find a remarkable overlap, suggesting the force distribution to reflect constraints for the evolutionary design of mechanical resistance in the IG family. The force distribution analysis provides a molecular interpretation of coevolution and opens the road to the study of the mechanism of signal propagation in proteins in general.

Protein allostery, signal transmission and dynamics: a classification scheme of allosteric mechanisms

Allostery has come of age; the number, breadth and functional roles of documented protein allostery cases are rising quickly. Since all dynamic proteins are potentially allosteric and allostery plays crucial roles in all cellular pathways, sorting and classifying allosteric mechanisms in proteins should be extremely useful in understanding and predicting how the signals are regulated and transmitted through the dynamic multi-molecular cellular organizations. Classification organizes the complex information thereby unraveling relationships and patterns in molecular activation and repression. In signaling, current classification schemes consider classes of molecules according to their functions; for example, epinephrine and norepinephrine secreted by the central nervous system are classified as neurotransmitters. Other schemes would account for epinephrine when secreted by the adrenal medulla to be hormone-like. Yet, such classifications account for the global function of the molecule; not for the molecular mechanism of how the signal transmission initiates and how it is transmitted. Here we provide a unified view of allostery and the first classification framework. We expect that a classification scheme would assist in comprehension of allosteric mechanisms, in prediction of signaling on the molecular level, in better comprehension of pathways and regulation of the complex signals, in translating them to the cascading events, and in allosteric drug design. We further provide a range of examples illustrating mechanisms in protein allostery and their classification from the cellular functional standpoint.

Allosteric communication occurs via networks of tertiary and quaternary motions in proteins

Allosteric proteins bind an effector molecule at one site resulting in a functional change at a second site. We hypothesize that allosteric communication in proteins relies upon networks of quaternary (collective, rigid-body) and tertiary (residue-residue contact) motions. We argue that cyclic topology of these networks is necessary for allosteric communication. An automated algorithm identifies rigid bodies from the displacement between the inactive and the active structures and constructs "quaternary networks" from these rigid bodies and the substrate and effector ligands. We then integrate quaternary networks with a coarse-grained representation of contact rearrangements to form "global communication networks" (GCNs). The GCN reveals allosteric communication among all substrate and effector sites in 15 of 18 multidomain and multimeric proteins, while tertiary and quaternary networks exhibit such communication in only 4 and 3 of these proteins, respectively. Furthermore, in 7 of the 15 proteins connected by the GCN, 50% or more of the substrate-effector paths via the GCN are "interdependent" paths that do not exist via either the tertiary or the quaternary network. Substrate-effector "pathways" typically are not linear but rather consist of polycyclic networks of rigid bodies and clusters of rearranging residue contacts. These results argue for broad applicability of allosteric communication based on structural changes and demonstrate the utility of the GCN. Global communication networks may inform a variety of experiments on allosteric proteins as well as the design of allostery into non-allosteric proteins.

Conformational changes and flexibility in T-cell receptor recognition of peptide-MHC complexes

A necessary feature of the immune system, TCR (T-cell receptor) cross-reactivity has been implicated in numerous autoimmune pathologies and is an underlying cause of transplant rejection. Early studies of the interactions of alphabeta TCRs (T-cell receptors) with their peptide-MHC ligands suggested that conformational plasticity in the TCR CDR (complementarity determining region) loops is a dominant contributor to T-cell cross-reactivity. Since these initial studies, the database of TCRs whose structures have been solved both bound and free is now large enough to permit general conclusions to be drawn about the extent of TCR plasticity and the types and locations of motion that occur. In the present paper, we review the conformational differences between free and bound TCRs, quantifying the structural changes that occur and discussing their possible roles in specificity and cross-reactivity. We show that, rather than undergoing major structural alterations or 'folding' upon binding, the majority of TCR CDR loops shift by relatively small amounts. The structural changes that do occur are dominated by hinge-bending motions, with loop remodelling usually occurring near loop apexes. As predicted from previous studies, the largest changes are in the hypervariable CDR3alpha and CDR3beta loops, although in some cases the germline-encoded CDR1alpha and CDR2alpha loops shift in magnitudes that approximate those of the CDR3 loops. Intriguingly, the smallest shifts are in the germline-encoded loops of the beta-chain, consistent with recent suggestions that the TCR beta domain may drive ligand recognition.

Moving towards high-resolution descriptions of the molecular interactions and structural rearrangements of the human hepatitis B core protein

The human hepatitis B virus core protein (HBc) forms icosahedral capsids and plays central roles in viral replication. The critical interactions that HBc makes prior to capsid formation (potential drug targets) have proved refractory to structural characterisation as HBc aggressively forms capsids. Our current structural understanding of HBc interactions is therefore capsid-centric, and this view has been limited by the resolution of cryo-electron microscopy and the inherent difficulties in getting high-quality crystals of viral capsids. To augment these approaches, we used capsid-dissociating conditions, solution NMR, and biophysical methodologies to directly characterise, at atomic resolution, the structural properties of dimeric HBc, as well as its dynamics and intermolecular interactions. Dimeric HBc recapitulates the structural properties and binding interactions of HBc within the context of capsids. Antiviral peptides induced long-range structural asymmetry in dimeric HBc, providing new insights into how ligand binding can effect communication between different regions of HBc and, therefore, between the capsid interior and the capsid exterior. Our work also paves the way for detailed descriptions of the previously invisible early stages of replication involving soluble HBc.

Ligand selectivity of D2 dopamine receptors is modulated by changes in local dynamics produced by sodium binding

We have uncovered a significant allosteric response of the D(2) dopamine receptor to physiologically relevant concentrations of sodium (140 mM), characterized by a sodium-enhanced binding affinity for a D(4)-selective class of agonists and antagonists. This enhancement is significantly more pronounced in a D(2)-V2.61(91)F mutant and cannot be mimicked by an equivalent concentration of the sodium replacement cation N-methyl-D-glucamine. This phenomenon was explored computationally at the molecular level by analyzing the effect of sodium binding on the dynamic properties of D(2) receptor model constructs. Normal mode analysis identified one mode (M(19)), which is involved in the open/closed motions of the binding cleft as being particularly sensitive to the sodium effect. To examine the consequences for D(2) receptor ligand recognition, one of the ligands, L-745,870 [3-{[4-(4-chlorophenyl) piperazin-1-yl]-methyl}-1H-pyrrolo[2,3-b]pyridine or CPPMA, chlorophenylpiperazinyl methylazaindole], was docked into conformers along the M(19) trajectory. Structurally and pharmacologically well established ligand-receptor interactions, including the ionic interaction with D3.32(114) and interactions between the ligand aryl moieties and V2.61(91)F, were achieved only in "open" phase conformers. The docking of (-)-raclopride [3,5-dichloro-N-(1-ethylpyrrolidin-2-ylmethyl)-2-hydroxy-6-methoxybenzamide] suggests that the same binding cleft changes in response to sodium-binding perturbation account as well for the enhancements in binding affinity for substituted benzamides in the wild-type D(2) receptor. Our findings demonstrate how key interactions can be modulated by occupancy at an allosteric site and are consistent with a mechanism in which sodium binding enhances the affinity of selected ligands through dynamic changes that increase accessibility of substituted benzamides and 1,4-DAP ligands to the orthosteric site and accessibility of 1,4-DAPs to V2.61(91)F.

Variations in clique and community patterns in protein structures during allosteric communication: investigation of dynamically equilibrated structures of methionyl tRNA synthetase complexes

The allosteric concept has played a key role in understanding the biological functions of proteins. The rigidity or plasticity and the conformational population are the two important ideas invoked in explaining the allosteric effect. Although molecular insights have been gained from a large number of structures, a precise assessment of the ligand-induced conformational changes in proteins at different levels, ranging from gross topology to intricate details, remains a challenge. In this study, we have explored the conformational changes in the complexes of methionyl tRNA synthetase (MetRS) through novel network parameters such as cliques and communities, which identify the rigid regions in the protein structure networks (PSNs) constructed from the noncovalent interactions of amino acid side chains. MetRS belongs to the aminoacyl tRNA synthetase (aaRS) family that plays a crucial role in the translation of genetic code. These enzymes are modular with distinct domains from which extensive genetic, kinetic, and structural data are available, highlighting the role of interdomain communication. The network parameters evaluated here on the conformational ensembles of MetRS complexes, generated from molecular dynamics simulations, have enabled us to understand the interdomain communication in detail. Additionally, the characterization of conformational changes in terms of cliques and communities has also become possible, which had eluded conventional analyses. Furthermore, we find that most of the residues participating in cliques and communities are strikingly different from those that take part in long-range communication. The cliques and communities evaluated here for the first time on PSNs have beautifully captured the local geometries in detail within the framework of global topology. Here the allosteric effect is revealed at the residue level via identification of the important residues specific for structural rigidity and functional flexibility in MetRS. This ought to enhance our understanding of the functioning of aaRS in general.

Residual ligand entropy in the binding of p-substituted benzenesulfonamide ligands to bovine carbonic anhydrase II

In studies on the thermodynamics of ligand-protein interactions, it is often assumed that the configurational and conformational entropy of the ligand is zero in the bound state (i.e., the ligand is rigidly fixed in the binding pocket). However, there is little direct experimental evidence for this assumption, and in the case of binding of p-substituted benzenesulfonamide inhibitors to bovine carbonic anhydrase II (BCA II), the observed thermodynamic binding signature derived from isothermal titration calorimetry experiments leads indirectly to the conclusion that a considerable degree of residual entropy remains in the bound ligand. Specifically, the entropy of binding increases with glycine chain length n, and strong evidence exists that this thermodynamic signature is not driven by solvent reorganization. By use of heteronuclear (15)N NMR relaxation measurements in a series (n = 1-6) of (15)N-glycine-enriched ligands, we find that the observed thermodynamic binding signature cannot be explained by residual ligand dynamics in the bound state, but rather results from the indirect influence of ligand chain length on protein dynamics.

Sulphate removal induces a major conformational change in Leishmania mexicana pyruvate kinase in the crystalline state

We report X-ray structures of pyruvate kinase from Leishmania mexicana (LmPYK) that are trapped in different conformations. These, together with the previously reported structure of LmPYK in its inactive (T-state) conformation, allow comparisons of three different conformers of the same species of pyruvate kinase (PYK). Four new site point mutants showing the effects of side-chain alteration at subunit interfaces are also enzymatically characterised. The LmPYK tetramer crystals grown with ammonium sulphate as precipitant adopt an active-like conformation, with sulphate ions at the active and effector sites. The sulphates occupy positions similar to those of the phosphates of ligands bound to active (R-state) and constitutively active (nonallosteric) PYKs from several species, and provide insight into the structural roles of the phosphates of the substrates and effectors. Crystal soaking in sulphate-free buffers was found to induce major conformational changes in the tetramer. In particular, the unwinding of the Aalpha6' helix and the inward hinge movement of the B domain are coupled with a significant widening (4 A) of the tetramer caused by lateral movement of the C domains. The two new LmPYK structures and the activity studies of site point mutations described in this article are consistent with a developing picture of allosteric activity in which localised changes in protein flexibility govern the distribution of conformer families adopted by the tetramer in its active and inactive states.

Molecular basis of halorespiration control by CprK, a CRP-FNR type transcriptional regulator

Certain bacteria are able to conserve energy via the reductive dehalogenation of halo-organic compounds in a respiration-type metabolism. The transcriptional regulator CprK from Desulfitobacterium spp. induces expression of halorespiratory genes upon binding of o-chlorophenol ligands and is reversibly inactivated by oxygen through disulphide bond formation. We report crystal structures of D. hafniense CprK in the ligand-free (both oxidation states), ligand-bound (reduced) and DNA-bound states, making it the first member of the widespread CRP-FNR superfamily for which a complete structural description of both redox-dependent and allosteric molecular rearrangements is available. In conjunction with kinetic and thermodynamic ligand binding studies, we provide a model for the allosteric mechanisms underpinning transcriptional control. Amino acids that play a key role in this mechanism are not conserved in functionally distinct CRP-FNR members. This suggests that, despite significant structural homology, distinct allosteric mechanisms are used, enabling this protein family to control a very wide range of processes.

Allosteric regulation and catalysis emerge via a common route

Allosteric regulation of protein function is a mechanism by which an event in one place of a protein structure causes an effect at another site, much like the behavior of a telecommunications network in which a collection of transmitters, receivers and transceivers communicate with each other across long distances. For example, ligand binding or an amino acid mutation at an allosteric site can alter enzymatic activity or binding affinity in a distal region such as the active site or a second binding site. The mechanism of this site-to-site communication is of great interest, especially since allosteric effects must be considered in drug design and protein engineering. In this review, conformational mobility as the common route between allosteric regulation and catalysis is discussed. We summarize recent experimental data and the resulting insights into allostery within proteins, and we discuss the nature of future studies and the new applications that may result from increased understanding of this regulatory mechanism.

Conformational and dynamics changes induced by bile acids binding to chicken liver bile acid binding protein

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.

Hydrophobic core mutations in CI2 globally perturb fast side-chain dynamics similarly without regard to position

Protein dynamics is currently an area of intense research because of its importance as complementary information to the huge quantity of available data relating protein structure and function. Because it is usually the influence of dynamics on function that is studied, the physical determinants of the distribution of flexibility in proteins have not been explored as thoroughly. In the present NMR study, an expanded suite of five (2)H relaxation experiments was used to characterize the picosecond-to-nanosecond side-chain dynamics of chymotrypsin inhibitor 2 (CI2) and five hydrophobic core mutants, some of which are members of the folding nucleus. Because CI2 is a homologue of the serine protease inhibitor eglin c, which has already been extensively characterized in terms of its dynamics, it was possible to compare not only side-chain dynamics but also the responses of these dynamics to analogous mutations. Remarkably, each of the five core mutations in CI2 led to similar and reproducible increases in side-chain flexibility throughout the entire structure. Although the expanded suite of (2)H relaxation experiments does not affect model selection for the vast majority of residues, it did enable the detection of increasing levels of nanosecond-scale motions in CI2's reactive site binding loop as the L68 side chain was progressively shortened by mutation. Collectively, we observed that the CI2 mutants are more dynamically similar to each other than to the more rigid wild-type CI2, from which we propose that wild-type CI2 has been optimized to a specific level of rigidity which may aid in its function as a serine protease inhibitor. We also observed that the pattern of side-chain dynamics of CI2 is quantitatively similar to eglin c, but that this similarity is lost upon mutating both proteins at an equivalent position. Finally, (15)N relaxation was used to characterize the backbone dynamics of wild-type and mutant CI2. Interestingly, mutation at folding nucleus positions led to widespread increases in backbone flexibility, whereas non-folding-nucleus positions led to increases in flexibility in the C-terminal half of the protein only.

Allostery and cooperativity revisited

Although phenomenlogical models that account for cooperativity in allosteric systems date back to the early and mid-60's (e.g., the KNF and MWC models), there is resurgent interest in the topic due to the recent experimental and computational studies that attempted to reveal, at an atomistic level, how allostery actually works. In this review, using systems for which atomistic simulations have been carried out in our groups as examples, we describe the current understanding of allostery, how the mechanisms go beyond the classical MWC/Pauling-KNF descriptions, and point out that the "new view" of allostery, emphasizing "population shifts," is, in fact, an "old view." The presentation offers not only an up-to-date description of allostery from a theoretical/computational perspective, but also helps to resolve several outstanding issues concerning allostery.

Mechanism of homotropic control to coordinate hydrolysis in a hexameric AAA+ ring ATPase

AAA(+) proteins are ubiquitous mechanochemical ATPases that use energy from ATP hydrolysis to remodel their versatile substrates. The AAA(+) characteristic hexameric ring assemblies raise important questions about if and how six often identical subunits coordinate hydrolysis and associated motions. The PspF AAA(+) domain, PspF(1-275), remodels the bacterial sigma(54)-RNA polymerase to activate transcription. Analysis of ATP substrate inhibition kinetics on ATP hydrolysis in hexameric PspF(1-275) indicates negative homotropic effects between subunits. Functional determinants required for allosteric control identify: (i) an important link between the ATP bound ribose moiety and the SensorII motif that would allow nucleotide-dependent *-helical */beta subdomain dynamics; and (ii) establishes a novel regulatory role for the SensorII helix in PspF, which may apply to other AAA(+) proteins. Consistent with functional data, homotropic control appears to depend on nucleotide state-dependent subdomain angles imposing dynamic symmetry constraints in the AAA(+) ring. Homotropic coordination is functionally important to remodel the sigma(54) promoter. We propose a structural symmetry-based model for homotropic control in the AAA(+) characteristic ring architecture.

Allosteric signaling in the biotin repressor occurs via local folding coupled to global dampening of protein dynamics

The biotin repressor is an allosterically regulated, site-specific DNA-binding protein. Binding of the small ligand bio-5'-AMP activates repressor dimerization, which is a prerequisite to DNA binding. Multiple disorder-to-order transitions, some of which are known to be important for the functional allosteric response, occur in the vicinity of the ligand-binding site concomitant with effector binding to the repressor monomer. In this work, the extent to which these local changes are coupled to additional changes in the structure/dynamics of the repressor was investigated using hydrogen/deuterium exchange coupled to mass spectrometry. Measurements were performed on the apo-protein and on complexes of the protein bound to four different effectors that elicit a range of thermodynamic responses in the repressor. Global exchange measurements indicate that binding of any effector to the intact protein is accompanied by protection from exchange. Mass spectrometric analysis of pepsin-cleavage products generated from the exchanged complexes reveals that the protection is distributed throughout the protein. Furthermore, the magnitude of the level of protection in each peptide from hydrogen/deuterium exchange correlates with the magnitude of the functional allosteric response elicited by a ligand. These results indicate that local structural changes in the binding site that occur concomitant with effector binding nucleate global dampening of dynamics. Moreover, the magnitude of dampening of repressor dynamics tracks with the magnitude of the functional response to effector binding.

Two-state allosteric modeling suggests protein equilibrium as an integral component for cyclic AMP (cAMP) specificity in the cAMP receptor protein of Escherichia coli

Activation of the cAMP receptor protein (CRP) from Escherichia coli is highly specific to its allosteric ligand, cAMP. Ligands such as adenosine and cGMP, which are structurally similar to cAMP, fail to activate wild-type CRP. However, several cAMP-independent CRP variants (termed CRP*) exist that can be further activated by both adenosine and cGMP, as well as by cAMP. This has remained a puzzle because the substitutions in many of these CRP* variants lie far from the cAMP-binding pocket (>10 A) and therefore should not directly affect that pocket. Here we show a surprising similarity in the altered ligand specificity of four CRP* variants with a single substitution in D53S, G141K, A144T, or L148K, and we propose a common basis for this phenomenon. The increased active protein population caused by an equilibrium shift in these variants is hypothesized to preferentially stabilize ligand binding. This explanation is completely consistent with the cAMP specificity in the activation of wild-type CRP. The model also predicts that wild-type CRP should be activated even by the lower-affinity ligand, adenosine, which we experimentally confirmed. The study demonstrates that protein equilibrium is an integral factor for ligand specificity in an allosteric protein, in addition to the direct effects of ligand pocket residues.

Allostery: absence of a change in shape does not imply that allostery is not at play

Allostery is essential for controlled catalysis, signal transmission, receptor trafficking, turning genes on and off, and apoptosis. It governs the organism's response to environmental and metabolic cues, dictating transient partner interactions in the cellular network. Textbooks taught us that allostery is a change of shape at one site on the protein surface brought about by ligand binding to another. For several years, it has been broadly accepted that the change of shape is not induced; rather, it is observed simply because a larger protein population presents it. Current data indicate that while side chains can reorient and rewire, allostery may not even involve a change of (backbone) shape. Assuming that the enthalpy change does not reverse the free-energy change due to the change in entropy, entropy is mainly responsible for binding.

Entropy-driven cAMP-dependent allosteric control of inhibitory interactions in exchange proteins directly activated by cAMP

Exchange proteins directly activated by cAMP (EPACs) are guanine nucleotide-exchange factors for the small GTPases Rap1 and Rap2 and represent a key receptor for the ubiquitous cAMP second messenger in eukaryotes. The cAMP-dependent activation of apoEPAC is typically rationalized in terms of a preexisting equilibrium between inactive and active states. Structural and mutagenesis analyses have shown that one of the critical determinants of the EPAC activation equilibrium is a cluster of salt bridges formed between the catalytic core and helices alpha1 and alpha2 at the N terminus of the cAMP binding domain and commonly referred to as ionic latch (IL). The IL stabilizes the inactive states in a closed topology in which access to the catalytic domain is sterically occluded by the regulatory moiety. However, it is currently not fully understood how the IL is allosterically controlled by cAMP. Chemical shift mapping studies consistently indicate that cAMP does not significantly perturb the structure of the IL spanning sites within the regulatory region, pointing to cAMP-dependent dynamic modulations as a key allosteric carrier of the cAMP-signal to the IL sites. Here, we have therefore investigated the dynamic profiles of the EPAC1 cAMP binding domain in its apo, cAMP-bound, and Rp-cAMPS phosphorothioate antagonist-bound forms using several 15N relaxation experiments. Based on the comparative analysis of dynamics in these three states, we have proposed a model of EPAC activation that incorporates the dynamic features allosterically modulated by cAMP and shows that cAMP binding weakens the IL by increasing its entropic penalty due to dynamic enhancements.

Monitoring aromatic picosecond to nanosecond dynamics in proteins via 13C relaxation: expanding perturbation mapping of the rigidifying core mutation, V54A, in eglin c

Long-range effects, such as allostery, have evolved in proteins as a means of regulating function via communication between distal sites. An NMR-based perturbation mapping approach was used to more completely probe the dynamic response of the core mutation V54A in the protein eglin c by monitoring changes in picosecond to nanosecond aromatic side-chain dynamics and H/D exchange stabilities. Previous side-chain dynamics studies on this mutant were limited to methyl-bearing residues, most of which were found to rigidify on the picosecond to nanosecond time scale in the form of a contiguous "network". Here, high precision (13)C relaxation data from 13 aromatic side chains were acquired by applying canonical relaxation experiments to a newly developed carbon labeling scheme [Teilum et al. (2006) J. Am. Chem. Soc. 128, 2506-2507]. The fitting of model-free parameters yielded S (2) variability which is intermediate with respect to backbone and methyl-bearing side-chain variability and tau e values that are approximately 1 ns. Inclusion of the aromatic dynamic response results in an expanded network of dynamically coupled residues, with some aromatics showing increases in flexibility, which partially offsets the rigidification in methyl side chains. Using amide hydrogen exchange, dynamic propagation on a slower time scale was probed in response to the V54A perturbation. Surprisingly, regional stabilization (slowed exchange) 10-12 A from the site of mutation was observed despite a global destabilization of 1.5 kcal x mol (-1). Furthermore, this unlikely pocket of stabilized residues colocalizes with increases in aromatic flexibility on the faster time scale. Because the converse is also true (destabilized residues colocalize with rigidification on the fast time scale), a plausible entropy-driven mechanism is discussed for relating colocalization of opposing dynamic trends on vastly different time scales.

Activation of G protein-coupled receptors: beyond two-state models and tertiary conformational changes

Transformation of G protein-coupled receptors (GPCRs) from a quiescent to an active state initiates signal transduction. All GPCRs share a common architecture comprising seven transmembrane-spanning alpha-helices, which accommodates signal propagation from a diverse repertoire of external stimuli across biological membranes to a heterotrimeric G protein. Signal propagation through the transmembrane helices likely involves mechanistic features common to all GPCRs. The structure of the light receptor rhodopsin may serve as a prototype for the transmembrane architecture of GPCRs. Early biochemical, biophysical, and pharmacological studies led to the conceptualization of receptor activation based on the context of two-state equilibrium models and conformational changes in protein structure. More recent studies indicate a need to move beyond these classical paradigms and to consider additional aspects of the molecular character of GPCRs, such as the oligomerization and dynamics of the receptor.

Reassessing a sparse energetic network within a single protein domain

Understanding the molecular principles that govern allosteric communication is an important goal in protein science. One way allostery could be transmitted is via sparse energetic networks of residues, and one such evolutionary conserved network was identified in the PDZ domain family of proteins by multiple sequence alignment [Lockless SW, Ranganathan R (1999) Science 286:295-299]. We have reassessed the energetic coupling of these residues by double mutant cycles together with ligand binding and stability experiments and found that coupling is not a special property of the coevolved network of residues in PDZ domains. The observed coupling for ligand binding is better explained by a distance relationship, where residues close in space are more likely to couple than distal residues. Our study demonstrates that statistical coupling from sequence analysis is not necessarily a reporter of energetic coupling and allostery.

Probing invisible, low-populated States of protein molecules by relaxation dispersion NMR spectroscopy: an application to protein folding

Biological function depends on molecular dynamics that lead to excursions from highly populated ground states to much less populated excited states. The low populations and the transient formation of such excited states render them invisible to the conventional methods of structural biology. Thus, while detailed pictures of ground-state structures of biomolecules have emerged over the years, largely through X-ray diffraction and solution nuclear magnetic resonance (NMR) spectroscopy studies, much less structural data has been accumulated on the conformational properties of the invisible excited states that are necessary to fully explain function. NMR spectroscopy is a powerful tool for studying conformational dynamics because it is sensitive to dynamics over a wide range of time scales, extending from picoseconds to seconds and because information is, in principle, available at nearly every position in the molecule. Here an NMR method for quantifying millisecond time scale dynamics that involve transitions between different molecular conformations is described. The basic experimental approach, termed relaxation dispersion NMR spectroscopy, is outlined to provide the reader with an intuitive feel for the technology. A variety of different experiments that probe conformational exchange at different sites in proteins are described, including a brief summary of data-fitting procedures to extract both the kinetic and thermodynamic properties of the exchange process and the structural features of the invisible excited states along the exchange pathway. It is shown that the methodology facilitates detection of intermediates and other excited states that are populated at low levels, 0.5% or higher, that cannot be observed directly in spectra, so long as they exchange with the observable ground state of the protein on the millisecond time scale. The power of the methodology is illustrated by a detailed application to the study of protein folding of the small modular SH3 domain. The kinetics and thermodynamics that describe the folding of this domain have been characterized through the effects of temperature, pressure, side-chain deuteration, and mutation, and the structural features of a low-populated folding intermediate have been assessed. Despite the fact that many previous studies have shown that SH3 domains fold via a two-state mechanism, the NMR methods presented unequivocally establish the presence of an on-pathway folding intermediate. The unique capabilities of NMR relaxation dispersion follow from the fact that large numbers of residues can be probed individually in a single experiment. By contrast, many other forms of spectroscopy monitor properties that are averaged over all residues in the molecule or that make use of only one or two reporters. The NMR methodology is not limited to protein folding, and applications to enzymatic catalysis, binding, and molecular recognition are beginning to emerge.

Compensatory and long-range changes in picosecond-nanosecond main-chain dynamics upon complex formation: 15N relaxation analysis of the free and bound states of the ubiquitin-like domain of human plexin-B1 and the small GTPase Rac1

The formation of a complex between Rac1 and the cytoplasmic domain of plexin-B1 is one of the first documented cases of a direct interaction between a small guanosine 5'-triphosphatase (GTPase) and a transmembrane receptor. Structural studies have begun to elucidate the role of this interaction for the signal transduction mechanism of plexins. Mapping of the Rac1 GTPase surface that contacts the Rho GTPase binding domain of plexin-B1 by solution NMR spectroscopy confirms the plexin domain as a GTPase effector protein. Regions neighboring the GTPase switch I and II regions are also involved in the interaction and there is considerable interest to examine the changes in protein dynamics that take place upon complex formation. Here we present main-chain nitrogen-15 relaxation measurements for the unbound proteins as well as for the Rho GTPase binding domain and Rac1 proteins in their complexed state. Derived order parameters, S2, show that considerable motions are maintained in the bound state of plexin. In fact, some of the changes in S2 on binding appear compensatory, exhibiting decreased as well as increased dynamics. Fluctuations in Rac1, already a largely rigid protein on the picosecond-nanosecond timescale, are overall diminished, but isomerization dynamics in the switch I and II regions of the GTPase are retained in the complex and appear to be propagated to the bound plexin domain. Remarkably, fluctuations in the GTPase are attenuated at sites, including helices alpha6 (the Rho-specific insert helix), alpha7 and alpha8, that are spatially distant from the interaction region with plexin. This effect of binding on long-range dynamics appears to be communicated by hinge sites and by subtle conformational changes in the protein. Similar to recent studies on other systems, we suggest that dynamical protein features are affected by allosteric mechanisms. Altered protein fluctuations are likely to prime the Rho GTPase-plexin complex for interactions with additional binding partners.

Efficient coupling of catalysis and dynamics in the E1 component of Escherichia coli pyruvate dehydrogenase multienzyme complex

Protein motions are ubiquitous and are intrinsically coupled to catalysis. Their specific roles, however, remain largely elusive. Dynamic loops at the active center of the E1 component of Escherichia coli pyruvate dehydrogenase multienzyme complex are essential for several catalytic functions starting from a predecarboxylation event and culminating in transfer of the acetyl moiety to the E2 component. Monitoring the kinetics of E1 and its loop variants at various solution viscosities, we show that the rate of a chemical step is modulated by loop dynamics. A cysteine-free E1 construct was site-specifically labeled on the inner loop (residues 401-413), and the EPR nitroxide label revealed ligand-induced conformational dynamics of the loop and a slow "open <--> close" conformational equilibrium in the unliganded state. An (19)F NMR label placed at the same residue revealed motion on the millisecond-second time scale and suggested a quantitative correlation of E1 catalysis and loop dynamics for the 200,000-Da protein. Thermodynamic studies revealed that these motions may promote covalent addition of substrate to the enzyme-bound thiamin diphosphate by reducing the free energy of activation. Furthermore, the global dynamics of E1 presumably regulate and streamline the catalytic steps of the overall complex by inducing an entirely entropic (nonmechanical) negative cooperativity with respect to substrate binding at higher temperatures. Our results are consistent with, and reinforce the hypothesis of, coupling of catalysis and regulation with enzyme dynamics and suggest the mechanism by which it is achieved in a key branchpoint enzyme in sugar metabolism.

Allosteric effects in the marginally stable von Hippel-Lindau tumor suppressor protein and allostery-based rescue mutant design

Many multifunctional tumor suppressor proteins have low stability, a property linked to cancer development. The von Hippel-Lindau tumor suppressor protein (pVHL) is one of these proteins. pVHL forms part of the E3 ubiquitin ligase complex that regulates the degradation of the hypoxia-inducible factor (HIF). Under native conditions, free pVHL is a molten globule, but it is stabilized in the E3 complex. By using molecular dynamics simulations, we observed that the interface between the two pVHL domains is the least stable region in unbound pVHL. We designed five stable mutants: one with a mutation at the interdomain interface and the others in the alpha- or beta-domains. Experimentally, type 2B pVHL disease mutant Y98N at the HIF binding site was shown to destabilize pVHL and decrease its binding affinity to HIF. Our simulations showed that the decrease in pVHL stability and binding affinity are allosterically regulated. The mutations designed to stabilize unbound wild-type pVHL, which are away from the elongin C and HIF binding sites, successfully stabilized the Y98N pVHL-elongin C complex and lowered the binding free energy of pVHL with HIF. Our results indicated both the enthalpic and dynamic allosteric components between the elongin C and HIF binding sites in pVHL, in the alpha- and beta-domains, respectively, mediated by the interdomain interface and linker. Drugs mimicking the allosteric effects of these mutants may rescue pVHL function in von Hippel-Lindau disease.

Allosteric cooperativity in protein kinase A

Allosteric signaling in proteins requires long-range communication mediated by highly conserved residues, often triggered by ligand binding. In this article, we map the allosteric network in the catalytic subunit of protein kinase A using NMR spectroscopy. We show that positive allosteric cooperativity is generated by nucleotide and substrate binding during the transitions through the major conformational states: apo, intermediate, and closed. The allosteric network is disrupted by a single site mutation (Y204A), which also decouples the cooperativity of ligand binding. Because protein kinase A is the prototype for the entire kinome, these findings may serve as a paradigm for describing long-range coupling in other protein kinases.

Functional consequences of exchanging domains between LacI and PurR are mediated by the intervening linker sequence

Homologue function can be differentiated by changing residues that affect binding sites or long-range interactions. LacI and PurR are two proteins that represent the LacI/GalR family (>500 members) of bacterial transcription regulators. All members have distinct DNA-binding and regulatory domains linked by approximately 18 amino acids. Each homologue has specificity for different DNA and regulatory effector ligands; LacI and PurR also exhibit differences in allosteric communication between DNA and effector binding sites. A comparative study of LacI and PurR suggested that alterations in the interface between the regulatory domain and linker are important for differentiating their functions. Four residues (equivalent to LacI positions 48, 55, 58, and 61) appear particularly important for creating a unique interface and were predicted to be necessary for allosteric regulation. However, nearby residues in the linker interact with DNA ligand. Thus, differences observed in interactions between linker and regulatory domain may be the cause of altered function or an effect of the two proteins binding different DNA ligands. To separate these possibilities, we created a chimeric protein with the LacI DNA-binding domain/linker and the PurR regulatory domain (LLhP). If the interface requires homologue-specific interactions in order to propagate the signal from effector binding, then LLhP repression should not be allosterically regulated by effector binding. Experiments show that LLhP is capable of repression from lacO1 and, contrary to expectation, allosteric response is intact. Further, restoring the potential for PurR-like interactions via substitutions in the LLhP linker tends to diminish repression. These effects are especially pronounced for residues 58 and 61. Clearly, binding affinity of LLhP for the lacO1 DNA site is sensitive to long-range changes in the linker. This result also raises the possibility that mutations at positions 58 and 61 co-evolved with changes in the DNA-binding site. In addition, repression measured in the absence and presence of effector ligand shows that allosteric response increases for several LLhP variants with substitutions at positions 48 and 55. Thus, while side chain variation at these sites does not generally dictate the presence or absence of allostery, the nature of the amino acid can modulate the response to effector.

Interdomain interaction of cyclic AMP receptor protein in the absence of cyclic AMP

Interdomain interaction of apo-cyclic AMP receptor protein (apo-CRP) was qualified using its isolated domains. The cAMP-binding domain was prepared by a limited proteolysis, while the DNA-binding domain was constructed as a recombinant protein. Three different regions making interdomain contacts in apo-CRP were identified by a sequence-specific comparison of the HSQC spectra. The results indicated that apo-CRP possesses characteristic modules of interdomain interaction that are properly organized to suppress activity and to sense and transfer the cAMP binding signals. Particularly, the inertness of the DNA-binding motif in apo-CRP was attributable to the participation of F-helices in the interdomain contacts.

A computational investigation of allostery in the catabolite activator protein

The catabolite activator protein is a dimer that consists of two cAMP-binding subunits, each containing a C-terminus DNA-binding module and a N-terminus ligand binding domain. The system is well-known to exhibit negative cooperativity, whereby the binding of one cAMP molecule reduces the binding affinity of the other cAMP molecule by 2 orders of magnitude, despite the large separation between the cAMP binding pockets. Here we use extensive explicit-solvent molecular dynamics simulations (135 ns) to investigate the allosteric mechanism of CAP. Six trajectories were carried out for apo, singly liganded, and doubly liganded CAP, both in the presence and absence of DNA. Thorough analyses of the dynamics through the construction of dynamical cross-correlated maps, as well as essential dynamics analyses, indicated that the system experienced a switch in motion as a result of cAMP binding, in accordance with recent NMR experiments carried out on a truncated form of the protein. Analyses of conformer structures collected from the simulations revealed a remarkable event: the DNA-binding module was found to dissociate from the N-terminus ligand binding domain. An interesting aspect of this structural change is that it only occurred in unoccupied subunits, suggesting that the binding of cAMP provides additional stability to the system, consistent with the increase in entropy that was observed in our calculations and from isothermal titration calorimetry. Analysis of the distribution of intrinsic disorder propensities in CAP amino acid sequence using PONDR VLXT and VSL1 predictors revealed that the region connecting ligand-binding and DNA-binding domains of CAP have the potential to exhibit increased flexibility. We complemented these trajectories with free energy calculations following the MM-PBSA approach on more than 2000 snapshots that included 880 normal mode analysis. The resulting free energy differences between the singly liganded and doubly liganded states were in excellent agreement with isothermal titration calorimetry data. When the free energy calculations were carried out in the presence of DNA, we discovered that a switch in cooperativity occurred, so that the binding of the first cAMP promoted the binding of the other cAMP. The components of the free energy reveal that this effect is mainly entropic in nature, whereby the DNA reduces the degree of tightening that is observed in its absence, thereby promoting binding of the second cAMP. This finding prompted us to propose a new mechanism by which CAP triggers the transcription activation that is based on an order to disorder transition mediated by cAMP binding as well as DNA.

Intrinsic dynamics of enzymes in the unbound state and relation to allosteric regulation

In recent years, there has been a surge in the number of studies exploring the relationship between proteins' equilibrium dynamics and structural changes involved in function. An emerging concept, supported by both theory and experiments, is that under native state conditions proteins have an intrinsic ability to sample conformations that meet functional requirements. A typical example is the ability of enzymes to sample open and closed forms, irrespective of substrate, succeeded by the stabilization of one form (usually closed) upon substrate binding. This ability is structure-encoded, and plays a key role in facilitating allosteric regulation, which suggests complementing the sequence-encodes-structure paradigm of protein science by structure-encodes-dynamics-encodes-function. The emerging connection implies an evolutionary role in selecting/conserving structures based on their ability to achieve functional dynamics, and in turn, selecting sequences that fold into such 'apt' structures.

Nonequilibrium dynamic mechanism for allosteric effect

Allosteric regulation is often viewed as thermodynamic in nature. However, protein internal motions during an enzymatic reaction cycle can slow the hoping processes over numerous potential barriers. We propose that regulating molecules may function by modifying the nonequilibrium protein dynamics. The theory predicts that an enzyme under the new mechanism has a different temperature dependence, waiting time distribution of the turnover cycle, and dynamic fluctuation patterns with and without an effector. Experimental tests of the theory are proposed.

All-atom simulations of structures and energetics of c-di-GMP-bound and free PleD

Cyclic diguanosine monophosphate is a bacterial second messenger involved in a lifestyle switch from single cells to biofilm formation. Atomistic simulations are used to characterize inhibited diguanylate cyclase (DGC) PleD with emphasis on the feedback inhibition mechanism. Normal-mode calculations show a rigidification particularly in both the inhibition site and the active site of the protein upon ligand binding. Extensive molecular dynamics simulations in explicit solvent and analysis of the dynamical cross-correlation maps suggest two distinct coupling pathways between the active and the inhibition site: direct information transfer either through the beta-strands beta2 and beta3 of the DGC domain (pathway I) or via the disordered regions connecting domains D2 and DGC (pathway II). In addition, dynamical cross-correlation maps show differences in the correlation between neighboring domains upon ligand binding and upon the point mutation R390A. The correlated motions between domains D1 and D2, which form the dimerization interface, are stronger for free PleD. Complementary to the experimentally observed short-range interactions in ligated PleD, the present work also characterizes the long-range, delocalized interactions between domains that are important for understanding activation and allosteric control of the protein. Based on the results, experimental characterization of the point mutant R353 and of the double mutant N357/H394 is proposed to differentiate between pathways I and II.

Coarse-grained free energy functions for studying protein conformational changes: a double-well network model

In this work, a double-well network model (DWNM) is presented for generating a coarse-grained free energy function that can be used to study the transition between reference conformational states of a protein molecule. Compared to earlier work that uses a single, multidimensional double-well potential to connect two conformational states, the DWNM uses a set of interconnected double-well potentials for this purpose. The DWNM free energy function has multiple intermediate states and saddle points, and is hence a "rough" free energy landscape. In this implementation of the DWNM, the free energy function is reduced to an elastic-network model representation near the two reference states. The effects of free energy function roughness on the reaction pathways of protein conformational change is demonstrated by applying the DWNM to the conformational changes of two protein systems: the coil-to-helix transition of the DB-loop in G-actin and the open-to-closed transition of adenylate kinase. In both systems, the rough free energy function of the DWNM leads to the identification of distinct minimum free energy paths connecting two conformational states. These results indicate that while the elastic-network model captures the low-frequency vibrational motions of a protein, the roughness in the free energy function introduced by the DWNM can be used to characterize the transition mechanism between protein conformations.

Structural dynamics in the gating ring of cyclic nucleotide-gated ion channels

For ligand-gated ion channels, the binding of a ligand to an intracellular or extracellular domain generates changes in transmembrane pore-forming helices, which alters ion flow. The molecular mechanism for this allostery, however, remains unknown. Here we explore the structure and conformational rearrangements of the C-terminal gating ring of the cyclic nucleotide-gated channel CNGA1 during activation by cyclic nucleotides with patch-clamp fluorometry. By monitoring fluorescent resonance energy transfer (FRET) between membrane-resident quenchers and fluorophores attached to the channel, we detected no movement orthogonal to the membrane during channel activation. By monitoring FRET between fluorophores within the C-terminal region, we determined that the C-terminal end of the C-linker and the end of the C-helix move apart when channels open. We conclude that during channel activation, a portion of the gating ring moves parallel to the plasma membrane, hinging toward the central axis of the channel.

From "simple" DNA-protein interactions to the macromolecular machines of gene expression

The physicochemical concepts that underlie our present ideas on the structure and assembly of the "macromolecular machines of gene expression" are developed, starting with the structure and folding of the individual protein and DNA components, the thermodynamics and kinetics of their conformational rearrangements during complex assembly, and the molecular basis of the sequence specificity and recognition interactions of the final assemblies that include the DNA genome. The role of diffusion in reduced dimensions in the kinetics of the assembly of macromolecular machines from their components is also considered, and diffusion-driven reactions are compared with those fueled by ATP binding and hydrolysis, as well as by the specific covalent chemical modifications involved in rearranging chromatin and modifying signal transduction networks in higher organisms.

Dynamic Allostery in the Ring Protein TRAP

We have discovered distinct, characteristic differences in the thermodynamic signatures of tryptophan binding by trp RNA-binding attenuation protein (TRAP) from two different bacterial species. The TRAP 11mer ring binds 11 molecules of tryptophan at symmetry-related sites. Tryptophan binding to Bacillus stearothermophilus TRAP is not cooperative, but isothermal titration calorimetry shows that filling the first tryptophan binding sites of Bacillus subtilis TRAP has a marked effect on the thermodynamics of subsequent ligand binding. We have identified a single, conservative amino acid replacement (Ile to Leu) in B. subtilis TRAP that abolishes this effect, and suggest the initial ligand binding causes a change throughout the wild-type protein ring.

Nucleation of an allosteric response via ligand-induced loop folding

The Escherichia coli biotin repressor BirA is an allosteric transcription regulatory protein to which binding of the small ligand corepressor biotinyl-5'-AMP promotes homodimerization and subsequent DNA binding. Structural data indicate that the apo or unliganded repressor is characterized by four partially disordered loops that are ordered in the ligand-bound dimer. While three of these loops participate directly in the dimerization, the fourth, consisting of residues 212-234 is distal to the interface. This loop, which is ordered around the adenine ring of the adenylate moiety in the BirA.adenylate structure, is referred to as the adenylate-binding loop (ABL). Although residues in the loop do not interact directly with the ligand, a hydrophobic cluster consisting of a tryptophan and two valine side-chains assembles over the adenine base. Results of previous measurements suggest that folding of the ABL is integral to the allosteric response. This idea and the role of the hydrophobic cluster in the process were investigated by systematic replacement of each side-chain in the cluster with alanine and analysis of the mutant proteins for small ligand binding and dimerization. Isothermal titration calorimetry measurements indicate defects in adenylate binding for all ABL variants. Additionally, sedimentation equilibrium measurements reveal that coupling between adenylate binding and dimerization is compromised in each mutant. Partial proteolysis measurements indicate that the mutants are defective in ligand-linked folding of the ABL. These results indicate that the hydrophobic cluster is critical to the ligand-induced disorder-to-order transition in the ABL and that this transition is integral to the allosteric response in the biotin repressor.

Intrinsic disorder as a mechanism to optimize allosteric coupling in proteins

Transcription factors and other allosteric cell signaling proteins contain a disproportionate number of domains or segments that are intrinsically disordered (ID) under native conditions. In many cases folding of these segments is coupled to binding with one or more of their interaction partners, suggesting that intrinsic disorder plays an important functional role. Despite numerous hypotheses for the role of ID domains in regulation, a mechanistic model has yet to be established that can quantitatively assess the importance of intrinsic disorder for intramolecular site-to-site communication, the hallmark property of allosteric proteins. Here, we present such a model and show that site-to-site allosteric coupling is maximized when intrinsic disorder is present in the domains or segments containing one or both of the coupled binding sites. This result not only explains the prevalence of ID domains in regulatory proteins, it also calls into question the classical mechanical view of energy propagation in proteins, which predicts that site-to-site coupling would be maximized when a well defined pathway of folded structure connects the two sites. Furthermore, in showing that the coupling mechanism conferred by intrinsic disorder is robust and independent of the network of interactions that physically link the coupled sites, unique insights are gained into the energetic ground rules that govern site-to-site communication in all proteins.

A single-quantum methyl 13C-relaxation dispersion experiment with improved sensitivity

A pulse sequence is described for recording single-quantum (13)C-methyl relaxation dispersion profiles of (13)C-selectively labeled methyl groups in proteins that offers significant improvements in sensitivity relative to existing approaches where initial magnetization derives from (13)C polarization. Sensitivity gains in the new experiment are achieved by making use of polarization from (1)H spins and (1)H --> (13)C --> (1)H type magnetization transfers. Its utility has been established by applications involving three different protein systems ranging in molecular weight from 8 to 28 kDa, produced using a number of different selective labeling approaches. In all cases exchange parameters from both (13)C-->(1)H and (1)H --> (13)C --> (1)H classes of experiment are in good agreement, with gains in sensitivity of between 1.7 and 4-fold realized using the new scheme.

Propagation of Dynamic Changes in Barnase Upon Binding of Barstar: An NMR and Computational Study

NMR spectroscopy and computer simulations were used to examine changes in chemical shifts and in dynamics of the ribonuclease barnase that result upon binding to its natural inhibitor barstar. Although the spatial structures of free and bound barnase are very similar, binding results in changes of the dynamics of both fast side-chains, as revealed by (2)H relaxation measurements, and NMR chemical shifts in an extended beta-sheet that is located far from the binding interface. Both side-chain dynamics and chemical shifts are sensitive to variations in the ensemble populations of the inter-converting molecular states, which can escape direct structural observation. Molecular dynamics simulations of free barnase and barnase in complex with barstar, as well as a normal mode analysis of barnase using a Gaussian network model, reveal relatively rigid domains that are separated by the extended beta-sheet mentioned above. The observed changes in NMR parameters upon ligation can thus be rationalized in terms of changes in inter-domain dynamics and in populations of exchanging states, without measurable structural changes. This provides an alternative model for the propagation of a molecular response to ligand binding across a protein that is based exclusively on changes in dynamics.

Proline cis-trans isomerization controls autoinhibition of a signaling protein

Autoinhibition is being widely used in nature to repress otherwise constitutive protein activities and is typically regulated by extrinsic factors. Here we show that autoinhibition can be controlled by an intrinsic intramolecular switch afforded by prolyl cis-trans isomerization. We find that a proline on the linker tethering the two SH3 domains of the Crk adaptor protein interconverts between the cis and trans conformation. In the cis conformation, the two SH3 domains interact intramolecularly, thereby forming the basis of an autoinhibitory mechanism. Conversely, in the trans conformation Crk exists in an extended, uninhibited conformation that is marginally populated but serves to activate the protein upon ligand binding. Interconversion between the cis and trans, and, hence, of the autoinhibited and activated conformations, is accelerated by the action of peptidyl-prolyl isomerases. Proline isomerization appears to make an ideal switch that can regulate the kinetics of activation, thereby modulating the dynamics of signal response.