Identification of dynamical hinge points of the L1 ligase molecular switch

The L1 ligase is an in vitro selected ribozyme that uses a noncanonically base-paired ligation site to catalyze regioselectively and regiospecifically the 5′ to 3′ phosphodiester bond ligation, a reaction relevant to origin of life hypotheses that invoke an RNA world scenario. The L1 ligase crystal structure revealed two different conformational states that were proposed to represent the active and inactive forms. It remains an open question as to what degree these two conformers persist as stable conformational intermediates in solution, and along what pathway are they able to interconvert. To explore these questions, we have performed a series of molecular dynamics simulations in explicit solvent of the inactive–active conformational switch in L1 ligase. Four simulations were performed departing from both conformers in both the reactant and product states, in addition to a simulation where local unfolding in the active state was induced. From these simulations, along with crystallographic data, a set of four virtual torsion angles that span two evolutionarily conserved and restricted regions were identified as dynamical hinge points in the conformational switch transition. The ligation site visits three distinct states characterized by hydrogen bond patterns that are correlated with the formation of specific contacts that may promote catalysis. The insights gained from these simulations contribute to a more detailed understanding of the coupled catalytic/conformational switch mechanism of L1 ligase that may facilitate the design and engineering of new catalytic riboswitches.     

Src kinase conformational activation: thermodynamics, pathways, and mechanisms

Tyrosine kinases of the Src-family are large allosteric enzymes that play a key role in cellular signaling. Conversion of the kinase from an inactive to an active state is accompanied by substantial structural changes. Here, we construct a coarse-grained model of the catalytic domain incorporating experimental structures for the two stable states, and simulate the dynamics of conformational transitions in kinase activation. We explore the transition energy landscapes by constructing a structural network among clusters of conformations from the simulations. From the structural network, two major ensembles of pathways for the activation are identified. In the first transition pathway, we find a coordinated switching mechanism of interactions among the αC helix, the activation-loop, and the β strands in the N-lobe of the catalytic domain. In a second pathway, the conformational change is coupled to a partial unfolding of the N-lobe region of the catalytic domain. We also characterize the switching mechanism for the αC helix and the activation-loop in detail. Finally, we test the performance of a Markov model and its ability to account for the structural kinetics in the context of Src conformational changes. Taken together, these results provide a broad framework for understanding the main features of the conformational transition taking place upon Src activation.     

Spin label studies on phosphorylase B.

Spin label studies upon phosphorylase B have revealed that the temperature-dependent equilibria between different conformational states are strongly influenced by the presence of substrate or allosteric modifier. In the absence of substrate or modifier, the equilibrium involves only two states, while at least three are involved in their presence. Removal of pyridoxal-5′-phosphate produces a conformational change reflected by increased mobility of the label. The apo-enzyme undergoes a structural change in the presence of the allosteric modifier AMP.     

Communication over the Network of Binary Switches Regulates the Activation of A2A Adenosine Receptor

Dynamics and functions of G-protein coupled receptors (GPCRs) are accurately regulated by the type of ligands that bind to the orthosteric or allosteric binding sites. To glean the structural and dynamical origin of ligand-dependent modulation of GPCR activity, we performed total ~ 5 μsec molecular dynamics simulations of A_2A adenosine receptor (A_2AAR) in its apo, antagonist-bound, and agonist-bound forms in an explicit water and membrane environment, and examined the corresponding dynamics and correlation between the 10 key structural motifs that serve as the allosteric hotspots in intramolecular signaling network. We dubbed these 10 structural motifs “binary switches” as they display molecular interactions that switch between two distinct states. By projecting the receptor dynamics on these binary switches that yield 2¹⁰ microstates, we show that (i) the receptors in apo, antagonist-bound, and agonist-bound states explore vastly different conformational space; (ii) among the three receptor states the apo state explores the broadest range of microstates; (iii) in the presence of the agonist, the active conformation is maintained through coherent couplings among the binary switches; and (iv) to be most specific, our analysis shows that W246, located deep inside the binding cleft, can serve as both an agonist sensor and actuator of ensuing intramolecular signaling for the receptor activation. Finally, our analysis of multiple trajectories generated by inserting an agonist to the apo state underscores that the transition of the receptor from inactive to active form requires the disruption of ionic-lock in the DRY motif.     

The role of conserved waters in conformational transitions of Q61H K-ras

To investigate the stability and functional role of long-residence water molecules in the Q61H variant of the signaling protein K-ras, we analyzed all available Ras crystal structures and conformers derived from a series of independent explicit solvent molecular dynamics (MD) simulations totaling 1.76 μs. We show that the protein samples a different region of phase space in the presence and absence of several crystallographically conserved and buried water molecules. The dynamics of these waters is coupled with the local as well as the global motions of the protein, in contrast to less buried waters whose exchange with bulk is only loosely coupled with the motion of loops in their vicinity. Aided by two novel reaction coordinates involving the distance (d) between the C(α) atoms of G60 at switch 2 and G10 at the P-loop and the N-C(α)-C-O dihedral (ξ) of G60, we further show that three water molecules located in lobe1, at the interface between the lobes and at lobe2, are involved in the relative motion of residues at the two lobes of Q61H K-ras. Moreover, a d/ξ plot classifies the available Ras x-ray structures and MD-derived K-ras conformers into active GTP-, intermediate GTP-, inactive GDP-bound, and nucleotide-free conformational states. The population of these states and the transition between them is modulated by water-mediated correlated motions involving the functionally critical switch 2, P-loop and helix 3. These results suggest that water molecules act as allosteric ligands to induce a population shift among distinct switch 2 conformations that differ in effector recognition.     

Allosteric Effects of Sodium Ion Binding on Activation of the M3 Muscarinic G-Protein-Coupled Receptor

G-protein-coupled receptors (GPCRs) are important membrane proteins that mediate cellular signaling and represent primary targets for about one-third of currently marketed drugs. Recent x-ray crystallographic studies identified distinct conformations of GPCRs in the active and inactive states. An allosteric sodium ion was found bound to a highly conserved D2.50 residue in inactive GPCRs, whereas the D2.50 allosteric pocket became collapsed in active GPCR structures. However, the dynamic mechanisms underlying these observations remain elusive. In this study, we aimed to understand the mechanistic effects of sodium ion binding on dynamic activation of the M3 muscarinic GPCR through long-timescale accelerated molecular dynamics (aMD) simulations. Results showed that with the D2.50 residue deprotonated, the M3 receptor is bound by an allosteric sodium ion and confined mostly in the inactive state with remarkably reduced flexibility. In contrast, the D2.50-protonated receptor does not exhibit sodium ion binding to the D2.50 allosteric site and samples a significantly larger conformational space. The receptor activation is captured and characterized by large-scale structural rearrangements of the transmembrane helices via dynamic hydrogen bond and salt bridge interactions. The residue motions are highly correlated during receptor activation. Further network analysis revealed that the allosteric signaling between residue D2.50 and key residues in the intracellular, extracellular, and orthosteric pockets is significantly weakened upon sodium ion binding.     

Mapping the nucleotide and isoform-dependent structural and dynamical features of Ras proteins

Ras GTPases are conformational switches controlling cell proliferation, differentiation, and development. Despite their prominent role in many forms of cancer, the mechanism of conformational transition between inactive GDP-bound and active GTP-bound states remains unclear. Here we describe a detailed analysis of available experimental structures and molecular dynamics simulations to quantitatively assess the structural and dynamical features of active and inactive states and their interconversion. We demonstrate that GTP-bound and nucleotide-free G12V H-ras sample a wide region of conformational space, and show that the inactive-to-active transition is a multiphase process defined by the relative rearrangement of the two switches and the orientation of Tyr32. We also modeled and simulated N- and K-ras proteins and found that K-ras is more flexible than N- and H-ras. We identified a number of isoform-specific, long-range side chain interactions that define unique pathways of communication between the nucleotide binding site and the C terminus.     

Conformational flexibility in the allosteric regulation of human UDP-α-D-glucose 6-dehydrogenase

UDP-α-D-xylose (UDX) acts as a feedback inhibitor of human UDP-α-D-glucose 6-dehydrogenase (hUGDH) by activating an unusual allosteric switch, the Thr131 loop. UDX binding induces the Thr131 loop to translate ~5 ? through the protein core, changing packing interactions and rotating a helix (α6(136-144)) to favor the formation of an inactive hexameric complex. But how does to conformational change occur given the steric packing constraints of the protein core? To answer this question, we deleted Val132 from the Thr131 loop to approximate an intermediate state in the allosteric transition. The 2.3 ? resolution crystal structure of the deletion construct (Δ132) reveals an open conformation that relaxes steric constraints and facilitates repacking of the protein core. Sedimentation velocity studies show that the open conformation stabilizes the Δ132 construct as a hexamer with point group symmetry 32, similar to that of the active complex. In contrast, the UDX-inhibited enzyme forms a lower-symmetry, horseshoe-shaped hexameric complex. We show that the Δ132 and UDX-inhibited structures have similar hexamer-building interfaces, suggesting that the hinge-bending motion represents a path for the allosteric transition between the different hexameric states. On the basis of (i) main chain flexibility and (ii) a model of the conformational change, we propose that hinge bending can occur as a concerted motion between adjacent subunits in the high-symmetry hexamer. We combine these results in a structurally detailed model for allosteric feedback inhibition and substrate--product exchange during the catalytic cycle.     

The projection analysis of NMR chemical shifts reveals extended EPAC autoinhibition determinants

EPAC is a cAMP-dependent guanine nucleotide exchange factor that serves as a prototypical molecular switch for the regulation of essential cellular processes. Although EPAC activation by cAMP has been extensively investigated, the mechanism of EPAC autoinhibition is still not fully understood. The steric clash between the side chains of two conserved residues, L273 and F300 in EPAC1, has been previously shown to oppose the inactive-to-active conformational transition in the absence of cAMP. However, it has also been hypothesized that autoinhibition is assisted by entropic losses caused by quenching of dynamics that occurs if the inactive-to-active transition takes place in the absence of cAMP. Here, we test this hypothesis through the comparative NMR analysis of several EPAC1 mutants that target different allosteric sites of the cAMP-binding domain (CBD). Using what to our knowledge is a novel projection analysis of NMR chemical shifts to probe the effect of the mutations on the autoinhibition equilibrium of the CBD, we find that whenever the apo/active state is stabilized relative to the apo/inactive state, dynamics are consistently quenched in a conserved loop (β2-β3) and helix (α5) of the CBD. Overall, our results point to the presence of conserved and nondegenerate determinants of CBD autoinhibition that extends beyond the originally proposed L273/F300 residue pair, suggesting that complete activation necessitates the simultaneous suppression of multiple autoinhibitory mechanisms, which in turn confers added specificity for the cAMP allosteric effector.     

Allosteric regulation of a ribozyme activity through ligand-induced conformational change.

An allosteric ribozyme has been designed using the hammerhead ribozyme as the active site and aflavin-specific RNA aptamer as a regulatory site. We constructed six variants with a series of base pairs in the linker region (stem II). Under single turnover conditions, kinetic studies were carried out in the absence and presence of flavin mononucleotide (FMN). Interestingly, FMN addition did not influence the cleavage rate of constructs with a 5-6 bp linker but stimulated the catalytic activity of those bearing a shorter linker. In particular, the apparent k cat of Rz3 increases by approximately 10-fold upon addition of saturating amounts of FMN. To determine the rate constants( K m4and k cat), the ribozyme regulated most effectively by FMN was further investigated. FMN mainly affected the k cat value, reflecting the rate limiting conformational change step of the overall cleavage reaction, depending on helix formation in stem II. Probably, FMN influences the orientation of structures necessary for the cleavage reaction through stem II formation. The result of chemical modification revealed that binding of FMN to the aptamer domain induced the helix formation in stem II required for catalytic activity. Therefore, a specific FMN-mediated allosteric interaction seems to promote a conformational alteration from an open to a closed structure in stem II. The concept of conformational modification in the allosteric effect is consistent with other allosteric enzymes, suggesting that such a conformational change is a fundamental feature of allosteric enzymes in biological systems.     

Conformational equilibrium of an enzyme catalytic site in the allosteric transition.

The dynamic equilibrium of a catalytic site between active and inactive conformations, the missing link between the structure and function of allosteric enzymes, was identified using protein engineering and NMR techniques. Kinetic analyses of the wild-type and three mutants of Thermus L-lactate dehydrogenase established that the allosteric property of the enzyme is associated with a concerted transition between the high-affinity (R) and low-affinity (T) states. By introducing mutations, we prepared an enzyme in which the R and T states were balanced. The conformation of the enzyme-bound coenzyme, NAD+, which interacts directly with the substrate, was analyzed using NMR spectroscopy. NAD+ bound to the mutant enzyme was in a conformational mixture of the active and inactive forms, while NAD+ took on predominantly one of the two forms when it was bound to the other enzymes we had analyzed. We interpret this to mean that the catalytic site is in equilibrium between the two conformations. The ratio of the conformers of each enzyme agreed with the [T]/[R] ratio as determined by kinetic analyses. Therefore, it is the identified conformational equilibrium of the catalytic site that governs the allosteric regulation of the enzyme activity.     

Design of an optical switch for studying conformational dynamics in individual molecules of GroEL

We describe the design of an optical switch in the chaperonin GroEL that is opened and closed by its ATP- and cochaperonin GroES-driven conformational changes. The switch, based on a fluorophore and a quencher, is engineered into the single-ring variant of the chaperone, and shows dramatic modulation of its fluorescent intensity in response to the transition of the protein between its allosteric states. It, therefore, forms a sensitive probe for the dynamics of the allosteric transitions of this machine, both in the bulk and in single molecules.     

Coarse-grained dynamics of the receiver domain of NtrC: fluctuations, correlations and implications for allosteric cooperativity

Receiver domains are key molecular switches in bacterial signaling. Structural studies have shown that the receiver domain of the nitrogen regulatory protein C (NtrC) exists in a conformational equilibrium encompassing both inactive and active states, with phosphorylation of Asp54 allosterically shifting the equilibrium towards the active state. To analyze dynamical fluctuations and correlations in NtrC as it undergoes activation, we have applied a coarse-grained dynamics algorithm using elastic network models. Normal mode analysis reveals possible dynamical pathways for the transition of NtrC from the inactive state to the active state. The diagonalized correlation between the inactive and the active (phosphorylated) state shows that most correlated motions occur around the active site of Asp54 and in the region Thr82 to Tyr101. This indicates a coupled correlation of dynamics in the "Thr82-Tyr101" motion. With phosphorylation inducing significant flexibility changes around the active site and alpha3 and alpha4 helices, we find that this activation makes the active-site region and the loops of alpha3/beta4 and alpha4/beta5 more stable. This means that phosphorylation entropically favors the receiver domain in its active state, and the induced conformational changes occur in an allosteric manner. Analyses of the local flexibility and long-range correlated motion also suggest a dynamics criterion for determining the allosteric cooperativity of NtrC, and may be applicable to other proteins.     

Neutralizing Antibodies Against Allosteric Proteins: Insights From a Bacterial Adhesin

Allosteric proteins transition between ‘inactive’ and ‘active’ states. In general, such proteins assume distinct conformational states at the level of secondary, tertiary and/or quaternary structure. Different conformers of an allosteric protein can be antigenically dissimilar and induce antibodies with a highly distinctive specificities and neutralizing functional effects. Here we summarize studies on various functional types of monoclonal antibodies obtained against different allosteric conformers of the mannose-specific bacterial adhesin FimH – the most common cell attachment protein of Escherichia coli and other enterobacterial pathogens. Included are types of antibodies that activate the FimH function via interaction with ligand-induced binding sites or by wedging between domains as well as antibodies that inhibit FimH through orthosteric, parasteric, or novel dynasteric mechanisms. Understanding the molecular mechanism of antibody action against allosteric proteins provides insights on how to design antibodies with a desired functional effect, including those with neutralizing activity against bacterial and viral cell attachment proteins.     

Segmental Motions, Not a Two-State Concerted Switch, Underlie Allostery in CheY

The switch between an inactive and active conformation is an important transition for signaling proteins, yet the mechanisms underlying such switches are not clearly understood. Escherichia coli CheY, a response regulator protein from the two-component signal transduction system that regulates bacterial chemotaxis, is an ideal protein for the study of allosteric mechanisms. By using (15)N CPMG relaxation dispersion experiments, we monitored the inherent dynamic switching of unphosphorylated CheY. We show that CheY does not undergo a two-state concerted switch between the inactive and active conformations. Interestingly, partial saturation of Mg(2+) enhances the intrinsic allosteric motions. Taken together with chemical shift perturbations, these data indicate that the μs-ms timescale motions underlying CheY allostery are segmental in nature. We propose an expanded allosteric network of residues, including W58, that undergo asynchronous, local switching between inactive and active-like conformations as the primary basis for the allosteric mechanism.     

Is allostery an intrinsic property of all dynamic proteins?

Allostery involves coupling of conformational changes between two widely separated binding sites. The common view holds that allosteric proteins are symmetric oligomers, with each subunit existing in "at least" two conformational states with a different affinity for ligands. Recent observations such as the allosteric behavior of myoglobin, a classical example of a nonallosteric protein, call into question the existing allosteric dogma. Here we argue that all (nonfibrous) proteins are potentially allosteric. Allostery is a consequence of re-distributions of protein conformational ensembles. In a nonallosteric protein, the binding site shape may not show a concerted second-site change and enzyme kinetics may not reflect an allosteric transition. Nevertheless, appropriate ligands, point mutations, or external conditions may facilitate a population shift, leading a presumably nonallosteric protein to behave allosterically. In principle, practically any potential drug binding to the protein surface can alter the conformational redistribution. The question is its effectiveness in the redistribution of the ensemble, affecting the protein binding sites and its function. Here, we review experimental observations validating this view of protein allostery.     

Molecular Dynamics Simulations and Classical Multidimensional Scaling Unveil New Metastable States in the Conformational Landscape of CDK2

Protein kinases are key regulatory nodes in cellular networks and their function has been shown to be intimately coupled with their structural flexibility. However, understanding the key structural mechanisms of large conformational transitions remains a difficult task. CDK2 is a crucial regulator of cell cycle. Its activity is finely tuned by Cyclin E/A and the catalytic segment phosphorylation, whereas its deregulation occurs in many types of cancer. ATP competitive inhibitors have failed to be approved for clinical use due to toxicity issues raised by a lack of selectivity. However, in the last few years type III allosteric inhibitors have emerged as an alternative strategy to selectively modulate CDK2 activity. In this study we have investigated the conformational variability of CDK2. A low dimensional conformational landscape of CDK2 was modeled using classical multidimensional scaling on a set of 255 crystal structures. Microsecond-scale plain and accelerated MD simulations were used to populate this landscape by using an out-of-sample extension of multidimensional scaling. CDK2 was simulated in the apo-form and in complex with the allosteric inhibitor 8-anilino-1-napthalenesulfonic acid (ANS). The apo-CDK2 landscape analysis showed a conformational equilibrium between an Src-like inactive conformation and an active-like form. These two states are separated by different metastable states that share hybrid structural features with both forms of the kinase. In contrast, the CDK2/ANS complex landscape is compatible with a conformational selection picture where the binding of ANS in proximity of the αC helix causes a population shift toward the inactive conformation. Interestingly, the new metastable states could enlarge the pool of candidate structures for the development of selective allosteric CDK2 inhibitors. The method here presented should not be limited to the CDK2 case but could be used to systematically unmask similar mechanisms throughout the human kinome.     

A molecular clamp ensures allosteric coordination of peptidyltransfer and ligand binding to the ribosomal A-site

Although the ribosome is mainly comprised of rRNA and many of its critical functions occur through RNA–RNA interactions, distinct domains of ribosomal proteins also participate in switching the ribosome between different conformational/functional states. Prior studies demonstrated that two extended domains of ribosomal protein L3 form an allosteric switch between the pre- and post-translocational states. Missing was an explanation for how the movements of these domains are communicated among the ribosome''s functional centers. Here, a third domain of L3 called the basic thumb, that protrudes roughly perpendicular from the W-finger and is nestled in the center of a cagelike structure formed by elements from three separate domains of the large subunit rRNA is investigated. Mutagenesis of basically charged amino acids of the basic thumb to alanines followed by detailed analyses suggests that it acts as a molecular clamp, playing a role in allosterically communicating the ribosome''s tRNA occupancy status to the elongation factor binding region and the peptidyltransferase center, facilitating coordination of their functions through the elongation cycle. The observation that these mutations affected translational fidelity, virus propagation and cell growth demonstrates how small structural changes at the atomic scale can propagate outward to broadly impact the biology of cell.