Millisecond dynamics in the allosteric enzyme imidazole glycerol phosphate synthase (IGPS) from Thermotoga maritima

IGPS is a 51 kDa heterodimeric enzyme comprised of two proteins, HisH and HisF, that catalyze the hydrolysis of glutamine to produce NH₃ in the HisH active site and the cyclization of ammonia with N′-[(5′-phosphoribulosyl)formimino]-5-aminoimidazole-4-carboxamide-ribonucleotide (PRFAR) in HisF to produce imidazole glycerol phosphate (IGP) and 5-aminoimidazole-4-carboxamide ribotide (AICAR). Binding of PRFAR and IGP stimulates glutaminase activity in the HisH enzyme over 5,000 and 100-fold, respectively, despite the active sites being >25 Å apart. The details of this long-range protein communication process were investigated by solution NMR spectroscopy and CPMG relaxation dispersion experiments. Formation of the heterodimer enzyme results in a reduction in millisecond motions in HisF that extend throughout the protein. Binding of lGP results in an increase in protein-wide millisecond dynamics evidenced as severe NMR line broadening and elevated R _ex values. Together, these data demonstrate a grouping of flexible residues that link the HisF active site with the protein interface to which HisH binds and provide a model for the path of communication between the IGPS active sites.     

A network of conserved interactions regulates the allosteric signal in a glutamine amidotransferase

We have combined equilibrium and steered molecular dynamics (SMD) simulations with principal component and correlation analyses to probe the mechanism of allosteric regulation in imidazole glycerol phosphate (IGP) synthase. An evolutionary analysis of IGP synthase revealed a conserved network of interactions leading from the effector binding site to the glutaminase active site, forming conserved communication pathways between the remote active sites. SMD simulations of the undocking of the ribonucleotide effector N1-[(5'-phosphoribulosyl)-formino]-5'-aminoimidazole carboxamide ribonucleotide (PRFAR) resulted in a large scale hinge-opening motion at the interface. Principal component analysis and a correlation analysis of the equilibration protein motion indicate that the dynamics involved in the allosteric transition are mediated by coupled motion between sites that are more than 25 A apart. Furthermore, conserved residues at the substrate-binding site, within the barrel, and at the interface were found to exhibit highly correlated motion during the allosteric transition. The coupled motion between PRFAR unbinding and the directed opening of the interface is interpreted in combination with kinetic assays for the wild-type and mutant systems to develop a model of allosteric regulation in IGP synthase that is monitored and investigated with atomic resolution.     

Comparison of conformational features of staphylococcal nuclease in ternary complexes with pdTp, pdGp, and nitrophenyl-pdTp.

The conformations of wild-type staphylococcal nuclease (SNase) in the ternary complexes with thymidine 3',5'-bisphosphate (pdTp), 2'-deoxyguanine 3',5'-bisphosphate (pdGp), and thymidine 3'-phosphate 5'-(p-nitrophenylphosphate) (NpdTp) with Ca2+ were examined by two-dimensional NMR NOESY and ROESY experiments. The results of these experiments indicate that the conformational features of the SNase are quite similar in the three ternary complexes. This suggests that the conformational features of SNase, in these ternary complexes, are not strongly dependent on whether the 5'-phosphate is a mono- or diester. This is in contrast to our prior studies on substitutions of active site charged amino acids which indicated that the conformational features of SNase in the ternary complex are quite sensitive to substitutions for active site charged amino acids (Hibler et al., 1987; Wilde et al., 1988; Pourmotabbed et al., 1990). The similarity of the SNase conformational features in the ternary complexes with pdTp and pdGp indicates that the features of the nucleotide bound at the active site are not strong determinants of the enzyme conformation in the ternary complexes. These conclusions are in general agreement with the results on pdApdT ternary complexes with SNase which suggested that it is the conformational features of the bound nucleic acid which determine the differences in catalysis observed for SNase with different substrates (Weber et al., 1991), more so than the conformational features of the enzyme.     

Evidence for two different mechanisms triggering the change in quaternary structure of the allosteric enzyme, glucosamine-6-phosphate deaminase

The generation and propagation of conformational changes associated with ligand binding in the allosteric enzyme glucosamine-6-phosphate deaminase (GlcN6P deaminase, EC 3.5.99.6) from Escherichia coli were analyzed by fluorescence measurements. Single-tryptophan mutant forms of the enzyme were constructed on the basis of previous structural and functional evidence and used as structural-change probes. The reporter residues were placed in the active-site lid (position 174) and in the allosteric site (254 and 234); in addition, signals from the natural Trp residues (15 and 224) were also studied as structural probes. The structural changes produced by the occupation of either the allosteric or the active site by site-specific ligands were monitored through changes in the spectral center of mass (SCM) of their steady-state emission fluorescence spectra. Binding of the allosteric activator produces only minimal signals in titration experiments. In contrast, measurable spectral signals were found when the active site was occupied by a dead-end inhibitor. The results reveal that the two binary complexes, enzyme-activator (R(A)) and enzyme-inhibitor (R(S)) complexes, have structural differences and that they also differ from the ternary complex (R(AS)). The mobility of the active-site lid motif is shown to be independent of the allosteric transition. The active-site ligand induces cooperative SCM changes even in the enzyme-activator complex, indicating that the propagation pathway of the conformational relaxation triggered from the active site is different from that involved in the heterotropic activation. Analysis of the complete set of mutants shows that the occupation of the active site generates structural perturbations, which are propagated to the whole of the monomer and extend to the other subunits. The accumulative effect of these propagated changes should be responsible for the change in the sign of the DeltaG degrees ' of the T to R transition associated with the progression of the active-site occupation, resulting in the predominance of the R over the T forms in the population of deaminase hexamers.     

The role of slow and fast protein motions in allosteric interactions

Allostery is fundamentally thermodynamic in nature. Long-range communication in proteins may be mediated not only by changes in the mean conformation with enthalpic contribution but also by changes in dynamic fluctuations with entropic contribution. The important role of protein motions in mediating allosteric interactions has been established by NMR spectroscopy. By using CAP as a model system, we have shown how changes in protein structure and internal dynamics can allosterically regulate protein function and activity. The results indicate that changes in conformational entropy can give rise to binding enhancement, binding inhibition, or have no effect in the expected affinity, depending on the magnitude and sign of enthalpy–entropy compensation. Moreover, allosteric interactions can be regulated by the modulation a low-populated conformation states that serve as on-pathway intermediates for ligand binding. Taken together, the interplay between fast internal motions, which are intimately related to conformational entropy, and slow internal motions, which are related to poorly populated conformational states, can regulate protein activity in a way that cannot be predicted on the basis of the protein’s ground-state structure.     

The interaction of ammonia and xenon with the imidazole glycerol phosphate synthase from Thermotoga maritima as detected by NMR spectroscopy

The imidazole glycerol phosphate (ImGP) synthase from the hyperthermophilic bacterium Thermotoga maritima is a 1:1 complex of the glutaminase subunit HisH and the cyclase subunit HisF. It has been proposed that ammonia generated by HisH is transported through a channel to the active site of HisF, which generates intermediates of histidine (ImGP) and de novo biosynthesis of 5‐aminoimidazole‐4‐carboxamideribotide. Solution NMR spectroscopy of ammonium chloride‐titrated samples was used to study the interaction of NH₃ with amino acids inside this channel. Although numerous residues showed ¹⁵N chemical shift changes, most of these changes were caused by nonspecific ionic strength effects. However, several interactions appeared to be specific. Remarkably, the amino acid residue Thr 78—which is located in the central channel—shows a large chemical shift change upon titration with ammonium chloride. This result and the reduced catalytic activity of the Thr78Met mutant indicate a special role of this residue in ammonia channeling. To detect and further characterize internal cavities in HisF, which might for example contribute to ammonia channeling, the interaction of HisF with the noble gas xenon was analyzed by solution NMR spectroscopy using ¹H‐¹⁵N HSQC experiments. The results indicate that HisF contains three distinct internal cavities, which could be identified by xenon‐induced chemical shift changes of the neighboring amino acid residues. Two of these cavities are located at the active site at opposite ends of the substrate N′‐[(5′‐phosphoribulosyl)formimino]‐5‐aminoimidazole‐4‐carboxamide‐ribonucleotide (PRFAR) binding groove. The third cavity is located in the interior of the central β‐barrel of HisF and overlaps with the putative ammonia transport channel.     

Delineation of the allosteric mechanism of a cytidylyltransferase exhibiting negative cooperativity.

The dimeric enzyme CTP:glycerol-3-phosphate cytidylyltransferase (GCT) displays strong negative cooperativity between the first and second binding of its substrate, CTP. Using NMR to study the allosteric mechanism of this enzyme, we observe widespread chemical shift changes for the individual CTP binding steps. Mapping these changes onto the molecular structure allowed the formulation of a detailed model of allosteric conformational change. Upon the second step of ligand binding, NMR experiments indicate an extensive loss of conformational exchange broadening of the backbone resonances of GCT. This suggests that a fraction of the free energy of negative cooperativity is entropic in origin.     

A structural basis for the allosteric regulation of non-hydrolysing UDP-GlcNAc 2-epimerases

The non-hydrolysing bacterial UDP-N-acetylglucosamine 2-epimerase (UDP-GlcNAc 2-epimerase) catalyses the conversion of UDP-GlcNAc into UDP-N-acetylmannosamine, an intermediate in the biosynthesis of several cell-surface polysaccharides. This enzyme is allosterically regulated by its substrate UDP-GlcNAc. The structure of the ternary complex between the Bacillus anthracis UDP-GlcNAc 2-epimerase, its substrate UDP-GlcNAc and the reaction intermediate UDP, showed direct interactions between UDP and its substrate, and between the complex and highly conserved enzyme residues, identifying the allosteric site of the enzyme. The binding of UDP-GlcNAc is associated with conformational changes in the active site of the enzyme. Kinetic data and mutagenesis of the highly conserved UDP-GlcNAc-interacting residues confirm their importance in the substrate binding and catalysis of the enzyme. This constitutes the first example to our knowledge, of an enzymatic allosteric activation by direct interaction between the substrate and the allosteric activator.     

Thumb inhibitor binding eliminates functionally important dynamics in the hepatitis C virus RNA polymerase

Hepatitis C virus (HCV) has infected almost 200 million people worldwide, typically causing chronic liver damage and severe complications such as liver failure. Currently, there are few approved treatments for viral infection. Thus, the HCV RNA‐dependent RNA polymerase (gene product NS5B) has emerged as an important target for small molecule therapeutics. Potential therapeutic agents include allosteric inhibitors that bind distal to the enzyme active site. While their mechanism of action is not conclusively known, it has been suggested that certain inhibitors prevent a conformational change in NS5B that is crucial for RNA replication. To gain insight into the molecular origin of long‐range allosteric inhibition of NS5B, we employed molecular dynamics simulations of the enzyme with and without an inhibitor bound to the thumb domain. These studies indicate that the presence of an inhibitor in the thumb domain alters both the structure and internal motions of NS5B. Principal components analysis identified motions that are severely attenuated by inhibitor binding. These motions may have functional relevance by facilitating interactions between NS5B and RNA template or nascent RNA duplex, with presence of the ligand leading to enzyme conformations with narrower and thus less accessible RNA binding channels. This study provides the first evidence for a mechanistic basis of allosteric inhibition in NS5B. Moreover, we present evidence that allosteric inhibition of NS5B results from intrinsic features of the enzyme free energy landscape, suggesting a common mechanism for the action of diverse allosteric ligands. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Conformational Transitions in Yeast Chorismate Mutase Important for Allosteric Regulation as Identified by Nuclear Magnetic Resonance Spectroscopy

Proteins fluctuate between different conformations in solution, and these conformational fluctuations can be important for protein function and allosteric regulation. The chorismate mutase from Saccharomyces cerevisiae (ScCM), a key enzyme in the biosynthesis of aromatic amino acids, is allosterically activated and inhibited by tryptophan and tyrosine, respectively. It was initially proposed that in the absence of effector, ScCM fluctuates between activated R and inhibited T conformations according to the Monod-Wyman-Changeux (MWC) model, although a more complex regulation pattern was later suggested by mutagenesis and kinetic data. Here we used NMR relaxation dispersion experiments to understand the conformational fluctuations on the microsecond-to-millisecond timescale that occur in ScCM. In the absence of allosteric effectors, ScCM did not exclusively exchange between T and R conformations, suggesting that the two-state MWC model is insufficient to explain conformational dynamics. Addition of tyrosine led to the quenching of much of the motion on this timescale, while new motions were identified in the presence of tryptophan. These new motions are consistent with conformational fluctuations into an alternative conformation that may be important for enzyme activity.     

3-Mercaptopicolinate. A reversible active site inhibitor of avian liver phosphoenolpyruvate carboxykinase

The inhibition of chicken liver phosphoenolpyruvate carboxykinase by 3-mercaptopicolinic acid (3-MP) has been investigated. Kinetic studies show 3-MP to be a noncompetitive inhibitor relative to all substrates and to the activator, Mn2+. EPR studies demonstrate that Mn2+ binding to the enzyme is unaffected by 3-MP. Proton relaxation rate studies demonstrate that 3-MP binds to the binary E X Mn complex with a KD of 0.5 X 10(-6) M and gives a ternary enhancement of 8.0. Additional proton relaxation rate studies detected formation of the quaternary complexes E X Mn X IDP X 3-MP, E X Mn X ITP X 3-MP, and E X Mn X CO2 X 3-MP. High resolution 1H nuclear relaxation rate studies suggest that 3-MP binds in close proximity to the activator cation, Mn2+, but not in the first coordination sphere. Active site models suggest that the 3-MP-binding site may partially overlap the phosphoenolpyruvate-binding site. The NMR studies, which detected formation of the quaternary E X Mn X 3-MP X phosphoenolpyruvate complex, also demonstrated that the binding of one of these ligands affects the interactions of the other ligand with E X Mn. Calorimetric studies of the E X Mn complex demonstrated that 3-MP causes an increase in the transition temperature midpoint without an increase in enthalpy. These results indicate that 3-MP causes a conformational change in the enzyme but does not increase the thermostability of the ternary complex. The experiments reported herein suggest that inhibition by 3-MP is due to specific and reversible binding within the active site of phosphoenolpyruvate carboxykinase.     

Alteration of hydrogen bonding in the vicinity of histidine 48 disrupts millisecond motions in RNase A

The motion of amino acid residues on the millisecond (ms) time scale is involved in the tight regulation of catalytic function in numerous enzyme systems. Using a combination of mutational, enzymological, and relaxation-compensated (15)N Carr-Purcell-Meiboom-Gill (CPMG) methods, we have previously established the conformational significance of the distant His48 residue and the neighboring loop 1 in RNase A function. These studies suggested that RNase A relies on an intricate network of hydrogen bonding interactions involved in propagating functionally relevant, long-range ms motions to the catalytic site of the enzyme. To further investigate the dynamic importance of this H-bonding network, this study focuses on the individual replacement of Thr17 and Thr82 with alanine, effectively altering the key H-bonding interactions that connect loop 1 and His48 to the rest of the protein. (15)N CPMG dispersion studies, nuclear magnetic resonance (NMR) chemical shift analysis, and NMR line shape analysis of point mutants T17A and T82A demonstrate that the evolutionarily conserved single H-bond linking His48 to Thr82 is essential for propagating ms motions from His48 to the active site of RNase A on the time scale of catalytic turnover, whereas the T17A mutation increases the off rate and conformational exchange motions in loop 1. Accumulating evidence from our mutational studies indicates that residues experiencing conformational exchange in RNase A can be grouped into two separate clusters displaying distinct dynamical features, which appear to be independently affected by mutation. Overall, this study illuminates how tightly controlled and finely tuned ms motions are in RNase A, suggesting that designed modulation of protein motions may be possible.     

Toward an understanding of the role of dynamics on enzymatic catalysis in lactate dehydrogenase

The motions of key residues at the substrate binding site of lactate dehydrogenase (LDH) were probed on the 10 ns to 10 ms time scale using laser-induced temperature-jump relaxation spectroscopy employing both UV fluorescence and isotope-edited IR absorption spectroscopy as structural probes. The dynamics of the mobile loop, which closes over the active site and is important for catalysis and binding, were characterized by studies of the inhibitor oxamate binding to the LDH/NADH binary complex monitoring the changes in emission of bound NADH. The bound NAD-pyruvate adduct, whose pyruvate moiety likely interacts with the same residues that interact with pyruvate in its ternary complex with LDH, served as a probe for any relative motions of active site residues against the substrate. The frequencies of its C=O stretch and -COO(-) antisymmetric stretch shift substantially should any relative motion of the polar moieties at the active site (His-195, Asp-168, Arg-109, and Arg-171) occur. The dynamics associated with loop closure are observed to involve several steps with motions from 1 to 300 microms. Apart from the "melting" of a few residues on the protein's surface, no kinetics were observed on any time scale in experiments of the bound NAD-pyr adduct although the measurements were made with a high degree of accuracy, even for final temperatures close to the unfolding transition of the protein. This is contrary to simple physical considerations and models. These results show that, once a productive protein/substrate complex is formed, the binding pocket is very rigid with very little, if any, motion apart from the mobile loop. The results also show that loop opening involves concomitant movement of the substrate out of the binding pocket.     

Dynamic modes of the flipped-out cytosine during HhaI methyltransferase-DNA interactions in solution.

Flipping of a nucleotide out of a B-DNA helix into the active site of an enzyme has been observed for the HhaI and HaeIII cytosine-5 methyltransferases (M.HhaI and M.HaeIII) and for numerous DNA repair enzymes. Here we studied the base flipping motions in the binary M. HhaI-DNA and the ternary M.HhaI-DNA-cofactor systems in solution. Two 5-fluorocytosines were introduced into the DNA in the places of the target cytosine and, as an internal control, a cytosine positioned two nucleotides upstream of the recognition sequence 5'-GCGC-3'. The 19F NMR spectra combined with gel mobility data show that interaction with the enzyme induces partition of the target base among three states, i.e. stacked in the B-DNA, an ensemble of flipped-out forms and the flipped-out form locked in the enzyme active site. Addition of the cofactor analogue S-adenosyl-L-homocysteine greatly enhances the trapping of the target cytosine in the catalytic site. Distinct dynamic modes of the target cytosine have thus been identified along the reaction pathway, which includes novel base-flipping intermediates that were not observed in previous X-ray structures. The new data indicate that flipping of the target base out of the DNA helix is not dependent on binding of the cytosine in the catalytic pocket of M.HhaI, and suggest an active role of the enzyme in the opening of the DNA duplex.     

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.     

ATP binding enables broad antibiotic selectivity of aminoglycoside phosphotransferase(3')-IIIa: an elastic network analysis

The bacterial enzyme aminoglycoside phosphotransferase(3′)-IIIa (APH) confers resistance against a wide range of aminoglycoside antibiotics. In this study, we use the Gaussian network model to investigate how the binding of nucleotides and antibiotics influences the dynamics and thereby the ligand binding properties of APH. Interestingly, in NMR experiments, the dynamics differ significantly in various APH complexes, although crystallographic studies indicate that no larger conformational changes occur upon ligand binding. Isothermal titration calorimetry also shows different thermodynamic contributions to ligand binding. Formation of aminoglycoside-APH complexes is enthalpically driven, while the enthalpic change upon aminoglycoside binding to the nucleotide-APH complex is much smaller. The differential effects of nucleotide binding and antibiotic binding to APH can be explained theoretically by single-residue fluctuations and correlated motions of the enzyme. The surprising destabilization of β-sheet residues upon nucleotide binding, as seen in hydrogen/deuterium exchange experiments, shows that the number of closest neighbors does not fully explain residue flexibility. Additionally, we must consider correlated motions of dynamic protein domains, which show that not only connectivity but also the overall protein architecture is important for protein dynamics.     

Long-range structural and dynamical changes induced by cofactor binding in DNA methyltransferase M.HhaI

The bacterial DNA cytosine methyltransferase M.HhaI sequence-specifically modifies DNA in an S-adenosylmethionine dependent reaction. The enzyme stabilizes the target cytosine (GCGC) into an extrahelical position, with a concomitant large movement of an active site loop involving residues 80-99. We used multidimensional, transverse relaxation-optimized NMR experiments to assign nearly 80% of all residues in the cofactor-bound enzyme form, providing a basis for detailed structural and dynamical characterization. We examined details of the previously unknown effects of the cofactor binding with M.HhaI in solution. Addition of the cofactor results in numerous structural changes throughout the protein, including those decorating the cofactor binding site, and distal residues more than 30 A away. The active site loop is involved in motions both on a picosecond to nanosecond time scale and on a microsecond to millisecond time scale and is not significantly affected by cofactor binding except for a few N-terminal residues. The cofactor also affects residues near the DNA binding cleft, suggesting a role for the cofactor in regulating DNA interactions. The allosteric properties we observed appear to be closely related to the significant amount of dynamics and dynamical changes in response to ligand binding detected in the protein.