The Impact of Second-Shell Nucleotides on Ligand Specificity in Cyclic Dinucleotide Riboswitches.

Ligand specificity is an essential requirement for all riboswitches. Some variant riboswitches utilize a common structural motif, yet through subtle sequence differences, they are able to selectively respond to different small molecule ligands and regulate downstream gene expression. These variants discriminate between structurally and chemically similar ligands. Crystal structures provide insight into how specificity is achieved. However, ligand specificity cannot always be explained solely by nucleotides in direct contact with the ligand. The cyclic dinucleotide variant family contains two classes, cyclic-di-GMP and cyclic-AMP-GMP riboswitches, that were distinguished based on the identity of a single nucleotide in contact with the ligand. Here we report a variant riboswitch with a mutation at a second ligand-contacting position that is promiscuous for both cyclic-di-GMP and cyclic-AMP-GMP despite a predicted preference for cyclic-AMP-GMP. A high-throughput mutational analysis, SMARTT, was used to quantitatively assess thousands of sites in the first- and second-shells of ligand contact for impacts on ligand specificity and promiscuity. In addition to nucleotides in direct ligand contact, nucleotides more distal from the binding site, within the J1/2 linker and the terminator helix, were identified that impact ligand specificity. These findings provide an example of how nucleotides outside the ligand binding pocket influence the riboswitch specificity. Moreover, these distal nucleotides could be used to predict promiscuous sequences.

Translation regulation by a guanidine-II riboswitch is highly tunable in sensitivity, dynamic range, and apparent cooperativity.

Riboswitches function as important translational regulators in bacteria. Comprehensive mutational analysis of transcriptional riboswitches has been used to probe the energetic intricacies of interplay between the aptamer and expression platform, but translational riboswitches have been inaccessible to massively parallel techniques. The guanidine-II (gdm-II) riboswitch is an exclusively translational class. We have integrated RelE cleavage with next-generation sequencing to quantify ligand-dependent changes in translation initiation for all single and double mutations of the Pseudomonas aeruginosa gdm-II riboswitch, a total of more than 23,000 variants. This extensive mutational analysis is consistent with the prominent features of the bioinformatic consensus. These data indicate, unexpectedly, that direct sequestration of the Shine-Dalgarno sequence is dispensable for riboswitch function. Additionally, this comprehensive data set reveals important positions not identified in previous computational and crystallographic studies. Mutations in the variable linker region stabilize alternate conformations. The double mutant data reveal the functional importance of the previously modeled P0b helix formed by the 5' and 3' tails that serves as the basis for translational control. Additional mutations to GU wobble base pairs in both P1 and P2 reveal how the apparent cooperativity of the system involves an intricate network of communication between the two binding sites. This comprehensive examination of a translational riboswitch's expression platform illuminates how the riboswitch is precisely tuned and tunable with regard to ligand sensitivity, the amplitude of expression between ON and OFF states, and the cooperativity of ligand binding.

The Positively Charged Active Site of the Bacterial Toxin RelE Causes a Large Shift in the General Base pKa.

The bacterial toxin RelE cleaves mRNA in the ribosomal A site. Although it shares a global fold with other microbial RNases, the active site contains several positively charged residues instead of histidines and glutamates that are typical of ribonucleases. The pH dependences of wild-type and mutant RelE indicate it uses general acid-base catalysis, but either the general acid (proposed to be R81) or the general base must have a substantially downshifted pKa. However, which group is shifted cannot be determined using available structural and biochemical data. Here, we use a phosphorothiolate at the scissile phosphate to remove the need for a general acid. We show this modification rescues nearly all of the defect of the R81A mutation, supporting R81 as the general acid. We also find that the observed pKa of the general base is dependent on the charge of the side chain at position 81. This indicates that positive charge in the active site contributes to a general base pKa downshifted by more than 5 units. Although this modestly reduces the effectiveness of general acid-base catalysis, it is strongly supplemented by the role of the positive charge in stabilizing the transition state for cleavage. Furthermore, we show that the ribosome is required for cleavage but not binding of mRNA by RelE. Ribosome functional groups do not directly contact the scissile phosphate, indicating that positioning and charge interactions dominate RelE catalysis. The unusual RelE active site catalyzes phosphoryl transfer at a rate comparable to those of similar enzymes, but in a ribosome-dependent fashion.

The asymmetry and cooperativity of tandem glycine riboswitch aptamers.

Glycine riboswitches utilize both single- and tandem-aptamer architectures. In the tandem system, the relative contribution of each aptamer toward gene regulation is not well understood. To dissect these contributions, the effects of 684 single mutants of a tandem ON switch from Bacillus subtilis were characterized for the wild-type construct and binding site mutations that selectively restrict ligand binding to either the first or second aptamer. Despite the structural symmetry of tandem aptamers, the response to these mutations was frequently asymmetrical. Mutations in the first aptamer often significantly weakened the K1/2, while several mutations in the second aptamer improved the amplitude. These results demonstrate that this ON switch favors ligand binding to the first aptamer. This is in contrast to the tandem OFF switch variant from Vibrio cholerae, which was previously shown to have preferential binding to its second aptamer. A bioinformatic analysis of tandem glycine riboswitches revealed that the two binding pockets are differentially conserved between ON and OFF switches. Altogether, this indicates that tandem ON switch variants preferentially utilize binding to the first aptamer to promote helical switching, while OFF switch variants favor binding to the second aptamer. The data set also revealed a cooperative glycine response when both binding pockets were maximally stabilized with three GC base pairs. This indicates a cooperative response may sometimes be obfuscated by a difference in the affinities of the two aptamers. This conditional cooperativity provides an additional layer of tunability to tandem glycine riboswitches that adds to their versatility as genetic switches.

Gene regulation by a glycine riboswitch singlet uses a finely tuned energetic landscape for helical switching.

Riboswitches contain structured aptamer domains that, upon ligand binding, facilitate helical switching in their downstream expression platforms to alter gene expression. To fully dissect how riboswitches function requires a better understanding of the energetic landscape for helical switching. Here, we report a sequencing-based high-throughput assay for monitoring in vitro transcription termination and use it to simultaneously characterize the functional effects of all 522 single point mutants of a glycine riboswitch type-1 singlet. Mutations throughout the riboswitch cause ligand-dependent defects, but only mutations within the terminator hairpin alter readthrough efficiencies in the absence of ligand. A comprehensive analysis of the expression platform reveals that ligand binding stabilizes the antiterminator by just 2-3 kcal/mol, indicating that the competing expression platform helices must be extremely close in energy to elicit a significant ligand-dependent response. These results demonstrate that gene regulation by this riboswitch is highly constrained by the energetics of ligand binding and conformational switching. These findings exemplify the energetic parameters of RNA conformational rearrangements driven by binding events.

Mycofumigation through production of the volatile DNA-methylating agent N-methyl-N-nitrosoisobutyramide by fungi in the genus Muscodor.

Antagonistic microorganisms produce antimicrobials to inhibit the growth of competitors. Although water-soluble antimicrobials are limited to proximal interactions via aqueous diffusion, volatile antimicrobials are able to act at a distance and diffuse through heterogeneous environments. Here, we identify the mechanism of action of Muscodor albus, an endophytic fungus known for its volatile antimicrobial activity toward a wide range of human and plant pathogens and its potential use in mycofumigation. Proposed uses of the Muscodor species include protecting crops, produce, and building materials from undesired fungal or bacterial growth. By analyzing a collection of Muscodor isolates with varying toxicity, we demonstrate that the volatile mycotoxin, N-methyl-N-nitrosoisobutyramide, is the dominant factor in Muscodor toxicity and acts primarily through DNA methylation. Additionally, Muscodor isolates exhibit higher resistance to DNA methylation compared with other fungi. This work contributes to the evaluation of Muscodor isolates as potential mycofumigants, provides insight into chemical strategies that organisms use to manipulate their environment, and provokes questions regarding the mechanisms of resistance used to tolerate constitutive, long-term exposure to DNA methylation.

Transition State Charge Stabilization and Acid-Base Catalysis of mRNA Cleavage by the Endoribonuclease RelE.

The bacterial toxin RelE is a ribosome-dependent endoribonuclease. It is part of a type II toxin-antitoxin system that contributes to antibiotic resistance and biofilm formation. During amino acid starvation, RelE cleaves mRNA in the ribosomal A-site, globally inhibiting protein translation. RelE is structurally similar to microbial RNases that employ general acid-base catalysis to facilitate RNA cleavage. The RelE active site is atypical for acid-base catalysis, in that it is enriched with positively charged residues and lacks the prototypical histidine-glutamate catalytic pair, making the mechanism of mRNA cleavage unclear. In this study, we use a single-turnover kinetic analysis to measure the effect of pH and phosphorothioate substitution on the rate constant for cleavage of mRNA by wild-type RelE and seven active-site mutants. Mutation and thio effects indicate a major role for stabilization of increased negative change in the transition state by arginine 61. The wild-type RelE cleavage rate constant is pH-independent, but the reaction catalyzed by many of the mutants is strongly dependent on pH, suggestive of general acid-base catalysis. pH-rate curves indicate that wild-type RelE operates with the pK(a) of at least one catalytic residue significantly downshifted by the local environment. Mutation of any single active-site residue is sufficient to disrupt this microenvironment and revert the shifted pK(a) back above neutrality. pH-rate curves are consistent with K54 functioning as a general base and R81 as a general acid. The capacity of RelE to effect a large pK(a) shift and facilitate a common catalytic mechanism by uncommon means furthers our understanding of other atypical enzymatic active sites.

Minimal transition state charge stabilization of the oxyanion during peptide bond formation by the ribosome.

Peptide bond formation during ribosomal protein synthesis involves an aminolysis reaction between the aminoacyl α-amino group and the carbonyl ester of the growing peptide via a transition state with a developing negative charge, the oxyanion. Structural and molecular dynamic studies have suggested that the ribosome may stabilize the oxyanion in the transition state of peptide bond formation via a highly ordered water molecule. To biochemically investigate this mechanistic hypothesis, we estimated the energetic contribution to catalytic charge stabilization of the oxyanion using a series of transition state mimics that contain different charge distributions and hydrogen bond potential on the functional group mimicking the oxyanion. Inhibitors containing an oxyanion mimic that carried a neutral charge and a mimic that preserved the negative charge but could not form hydrogen bonds had less than a 3-fold effect on inhibitor binding affinity. These observations argue that the ribosome provides minimal transition state charge stabilization to the oxyanion during peptide bond formation via the water molecule. This is in contrast to the substantial level of oxyanion stabilization provided by serine proteases. This suggests that the oxyanion may be neutralized via a proton shuttle, resulting in an uncharged transition state.

The chemical versatility of RNA.

The ability of RNA to both store genetic information and catalyse chemical reactions has led to the hypothesis that it predates DNA and proteins. While there is no doubt that RNA is capable of storing the genetic information of a primitive organism, only two classes of reactions-phosphoryl transfer and peptide bond formation-have been observed to be catalysed by RNA in nature. However, these naturally occurring ribozymes use a wide range of catalytic strategies that could be applied to other reactions. Furthermore, RNA can bind several cofactors that are used by protein enzymes to facilitate a wide variety of chemical processes. Despite its limited functional groups, these observations indicate RNA is a versatile molecule that could, in principle, catalyse the myriad reactions necessary to sustain life.

A two-step chemical mechanism for ribosome-catalysed peptide bond formation.

The chemical step of natural protein synthesis, peptide bond formation, is catalysed by the large subunit of the ribosome. Crystal structures have shown that the active site for peptide bond formation is composed entirely of RNA. Recent work has focused on how an RNA active site is able to catalyse this fundamental biological reaction at a suitable rate for protein synthesis. On the basis of the absence of important ribosomal functional groups, lack of a dependence on pH, and the dominant contribution of entropy to catalysis, it has been suggested that the role of the ribosome is limited to bringing the substrates into close proximity. Alternatively, the importance of the 2'-hydroxyl of the peptidyl-transfer RNA and a Brønsted coefficient near zero have been taken as evidence that the ribosome coordinates a proton-transfer network. Here we report the transition state of peptide bond formation, based on analysis of the kinetic isotope effect at five positions within the reaction centre of a peptidyl-transfer RNA mimic. Our results indicate that in contrast to the uncatalysed reaction, formation of the tetrahedral intermediate and proton transfer from the nucleophilic nitrogen both occur in the rate-limiting step. Unlike in previous proposals, the reaction is not fully concerted; instead, breakdown of the tetrahedral intermediate occurs in a separate fast step. This suggests that in addition to substrate positioning, the ribosome is contributing to chemical catalysis by changing the rate-limiting transition state.

The energetic contribution of induced electrostatic asymmetry to DNA bending by a site-specific protein.

DNA bending can be promoted by reducing the net negative electrostatic potential around phosphates on one face of the DNA, such that electrostatic repulsion among phosphates on the opposite face drives bending toward the less negative surface. To provide the first assessment of energetic contribution to DNA bending when electrostatic asymmetry is induced by a site-specific DNA binding protein, we manipulated the electrostatics in the EcoRV endonuclease-DNA complex by mutation of cationic side chains that contact DNA phosphates and/or by replacement of a selected phosphate in each strand with uncharged methylphosphonate. Reducing the net negative charge at two symmetrically located phosphates on the concave DNA face contributes -2.3 kcal mol(-1) to -0.9 kcal mol(-1) (depending on position) to complex formation. In contrast, reducing negative charge on the opposing convex face produces a penalty of +1.3 kcal mol(-1). Förster resonance energy transfer experiments show that the extent of axial DNA bending (about 50°) is little affected in modified complexes, implying that modification affects the energetic cost but not the extent of DNA bending. Kinetic studies show that the favorable effects of induced electrostatic asymmetry on equilibrium binding derive primarily from a reduced rate of complex dissociation, suggesting stabilization of the specific complex between protein and markedly bent DNA. A smaller increase in the association rate may suggest that the DNA in the initial encounter complex is mildly bent. The data imply that protein-induced electrostatic asymmetry makes a significant contribution to DNA bending but is not itself sufficient to drive full bending in the specific EcoRV-DNA complex.

Transition states of uncatalyzed hydrolysis and aminolysis reactions of a ribosomal P-site substrate determined by kinetic isotope effects.

The ester bond of peptidyl-tRNA undergoes nucleophilic attack in solution and when catalyzed by the ribosome. To characterize the uncatalyzed hydrolysis reaction, a model of peptide release, the transition state structure for hydrolysis of a peptidyl-tRNA mimic was determined. Kinetic isotope effects were measured at several atoms that potentially undergo a change in bonding in the transition state. Large kinetic isotope effects of carbonyl (18)O and alpha-deuterium substitutions on uncatalyzed hydrolysis indicate the transition state is nearly tetrahedral. Kinetic isotope effects were also measured for aminolysis by hydroxylamine to study a reaction similar to the formation of a peptide bond. In contrast to hydrolysis, the large leaving group (18)O isotope effect indicates the C-O3' bond has undergone significant scission in the transition state. The smaller carbonyl (18)O and alpha-deuterium effects are consistent with a later transition state. The assay developed here can also be used to measure isotope effects for the ribosome-catalyzed reactions. These uncatalyzed reactions serve as a basis for determining what aspects of the transition states are stabilized by the ribosome to achieve a rate enhancement.

Positively charged C-terminal subdomains of EcoRV endonuclease: contributions to DNA binding, bending, and cleavage.

The carboxy-terminal subdomains of the homodimeric EcoRV restriction endonuclease each bear a net charge of +4 and are positioned on the inner concave surface of the 50 degree DNA bend that is induced by the enzyme. A complete kinetic and structural analysis of a truncated EcoRV mutant lacking these domains was performed to assess the importance of this diffuse charge in facilitating DNA binding, bending, and cleavage. At the level of formation of an enzyme-DNA complex, the association rate for the dimeric mutant enzyme was sharply decreased by 10(3)-fold, while the equilibrium dissociation constant was weakened by nearly 10(6)-fold compared with that of wild-type EcoRV. Thus, the C-terminal subdomains strongly stabilize the enzyme-DNA ground-state complex in which the DNA is known to be bent. Further, the extent of DNA bending as observed by fluorescence resonance energy transfer was also significantly decreased. The crystal structure of the truncated enzyme bound to DNA and calcium ions at 2.4 A resolution reveals that the global fold is preserved and suggests that a divalent metal ion crucial to catalysis is destabilized in the active site. This may explain the 100-fold decrease in the rate of metal-dependent phosphoryl transfer observed for the mutant. These results show that diffuse positive charge associated with the C-terminal subdomains of EcoRV plays a key role in DNA association, bending, and cleavage.

Non-cognate enzyme-DNA complex: structural and kinetic analysis of EcoRV endonuclease bound to the EcoRI recognition site GAATTC.

The crystal structure of EcoRV endonuclease bound to non-cognate DNA at 2.0 angstroms resolution shows that very small structural adaptations are sufficient to ensure the extreme sequence specificity characteristic of restriction enzymes. EcoRV bends its specific GATATC site sharply by 50 degrees into the major groove at the center TA step, generating unusual base-base interactions along each individual DNA strand. In the symmetric non-cognate complex bound to GAATTC, the center step bend is relaxed to avoid steric hindrance caused by the different placement of the exocyclic thymine methyl groups. The decreased base-pair unstacking in turn leads to small conformational rearrangements in the sugar-phosphate backbone, sufficient to destabilize binding of crucial divalent metal ions in the active site. A second crystal structure of EcoRV bound to the base-analog GAAUTC site shows that the 50 degrees center-step bend of the DNA is restored. However, while divalent metals bind at high occupancy in this structure, one metal ion shifts away from binding at the scissile DNA phosphate to a position near the 3'-adjacent phosphate group. This may explain why the 10(4)-fold attenuated cleavage efficiency toward GAATTC is reconstituted by less than tenfold toward GAAUTC. Examination of DNA binding and bending by equilibrium and stopped-flow florescence quenching and fluorescence resonance energy transfer (FRET) methods demonstrates that the capacity of EcoRV to bend the GAATTC non-cognate site is severely limited, but that full bending of GAAUTC is achieved at only a threefold reduced rate compared with the cognate complex. Together, the structural and biochemical data demonstrate the existence of distinct mechanisms for ensuring specificity at the bending and catalytic steps, respectively. The limited conformational rearrangements observed in the EcoRV non-cognate complex provide a sharp contrast to the extensive structural changes found in a non-cognate BamHI-DNA crystal structure, thus demonstrating a diversity of mechanisms by which restriction enzymes are able to achieve specificity.

Simultaneous DNA binding and bending by EcoRV endonuclease observed by real-time fluorescence.

The complete catalytic cycle of EcoRV endonuclease has been observed by combining fluorescence anisotropy with fluorescence resonance energy transfer (FRET) measurements. Binding, bending, and cleavage of substrate oligonucleotides were monitored in real time by rhodamine-x anisotropy and by FRET between rhodamine and fluorescein dyes attached to opposite ends of a 14-mer DNA duplex. For the cognate GATATC site binding and bending are found to be nearly simultaneous, with association and bending rate constants of (1.45-1.6) x 10(8) M(-1) s(-1). On the basis of the measurement of k(off) by a substrate-trapping approach, the equilibrium dissociation constant of the enzyme-DNA complex in the presence of inhibitory calcium ions was calculated as 3.7 x 10(-12) M from the kinetic constants. Further, the entire DNA cleavage reaction can be observed in the presence of catalytic Mg(2+) ions. These measurements reveal that the binding and bending steps occur at equivalent rates in the presence of either Mg(2+) or Ca(2+), while a slow decrease in fluorescence intensity following bending corresponds to k(cat), which is limited by the cleavage and product dissociation steps. Measurement of k(on) and k(off) in the absence of divalent metals shows that the DNA binding affinity is decreased by 5000-fold to 1.4 x 10(-8) M, and no bending could be detected in this case. Together with crystallographic studies, these data suggest a model for the induced-fit conformational change in which the role of divalent metal ions is to stabilize the sharply bent DNA in an orientation suitable for accessing the catalytic transition state.