Metal-ion binding and metal-ion induced folding of the adenine-sensing riboswitch aptamer domain

Divalent cations are important in the folding and stabilization of complex RNA structures. The adenine-sensing riboswitch controls the expression of mRNAs for proteins involved in purine metabolism by directly sensing intracellular adenine levels. Adenine binds with high affinity and specificity to the ligand binding or aptamer domain of the adenine-sensing riboswitch. The X-ray structure of this domain in complex with adenine revealed an intricate RNA-fold consisting of a three-helix junction stabilized by long-range base-pairing interactions and identified five binding sites for hexahydrated Mg²⁺-ions. Furthermore, a role for Mg²⁺-ions in the ligand-induced folding of this RNA was suggested. Here, we describe the interaction of divalent cations with the RNA–adenine complex in solution as studied by high-resolution NMR spectroscopy. Paramagnetic line broadening, chemical shift mapping and intermolecular nuclear Overhauser effects (NOEs) indicate the presence of at least three binding sites for divalent cations. Two of them are similar to those in the X-ray structure. The third site, which is important for the folding of this RNA, has not been observed previously. The ligand-free state of the RNA is conformationally heterogeneous and contains base-pairing patterns detrimental to ligand binding in the absence of Mg²⁺, but becomes partially pre-organized for ligand binding in the presence of Mg²⁺. Compared to the highly similar guanine-sensing riboswitch, the folding pathway for the adenine-sensing riboswitch aptamer domain is more complex and the influence of Mg²⁺ is more pronounced.     

Mg2+-induced conformational changes in the btuB riboswitch from E. coli

Mg²⁺ has been shown to modulate the function of riboswitches by facilitating the ligand-riboswitch interactions. The btuB riboswitch from Escherichia coli undergoes a conformational change upon binding to its ligand, coenzyme B₁₂ (adenosyl-cobalamine, AdoCbl), and down-regulates the expression of the B₁₂ transporter protein BtuB in order to control the cellular levels of AdoCbl. Here, we discuss the structural folding attained by the btuB riboswitch from E. coli in response to Mg²⁺ and how it affects the ligand binding competent conformation of the RNA. The btuB riboswitch notably adopts different conformational states depending upon the concentration of Mg²⁺. With the help of in-line probing, we show the existence of at least two specific conformations, one being achieved in the complete absence of Mg²⁺ (or low Mg²⁺ concentration) and the other appearing above ∼0.5 mM Mg²⁺. Distinct regions of the riboswitch exhibit different dissociation constants toward Mg²⁺, indicating a stepwise folding of the btuB RNA. Increasing the Mg²⁺ concentration drives the transition from one conformation toward the other. The conformational state existing above 0.5 mM Mg²⁺ defines the binding competent conformation of the btuB riboswitch which can productively interact with the ligand, coenzyme B₁₂, and switch the RNA conformation. Moreover, raising the Mg²⁺ concentration enhances the ratio of switched RNA in the presence of AdoCbl. The lack of a AdoCbl-induced conformational switch experienced by the btuB riboswitch in the absence of Mg²⁺ indicates a crucial role played by Mg²⁺ for defining an active conformation of the riboswitch.     

A recurrent magnesium-binding motif provides a framework for the ribosomal peptidyl transferase center

The ribosome is an ancient macromolecular machine responsible for the synthesis of all proteins in all living organisms. Here we demonstrate that the ribosomal peptidyl transferase center (PTC) is supported by a framework of magnesium microclusters (Mg²⁺-μc's). Common features of Mg²⁺-μc's include two paired Mg²⁺ ions that are chelated by a common bridging phosphate group in the form Mg_(a)²⁺–(O1P-P-O2P)–Mg_(b)²⁺. This bridging phosphate is part of a 10-membered chelation ring in the form Mg_(a)²⁺–(OP-P-O5′-C5′-C4′-C3′-O3′-P-OP)–Mg_(a)²⁺. The two phosphate groups of this 10-membered ring are contributed by adjacent residues along the RNA backbone. Both Mg²⁺ ions are octahedrally coordinated, but are substantially dehydrated by interactions with additional RNA phosphate groups. The Mg²⁺-μc's in the LSU (large subunit) appear to be highly conserved over evolution, since they are unchanged in bacteria (Thermus thermophilus, PDB entry 2J01) and archaea (Haloarcula marismortui, PDB entry 1JJ2). The 2D elements of the 23S rRNA that are linked by Mg²⁺-μc's are conserved between the rRNAs of bacteria, archaea and eukarya and in mitochondrial rRNA, and in a proposed minimal 23S-rRNA. We observe Mg²⁺-μc's in other rRNAs including the bacterial 16S rRNA, and the P4–P6 domain of the tetrahymena Group I intron ribozyme. It appears that Mg²⁺-μc's are a primeval motif, with pivotal roles in RNA folding, function and evolution.     

In vitro selection and characterization of a highly efficient Zn(II)-dependent RNA-cleaving deoxyribozyme

A group of highly efficient Zn(II)-dependent RNA-cleaving deoxyribozymes has been obtained through in vitro selection. They share a common motif with the ‘8–17’ deoxyribozyme isolated under different conditions, including different design of the random pool and metal ion cofactor. We found that this commonly selected motif can efficiently cleave both RNA and DNA/RNA chimeric substrates. It can cleave any substrate containing rNG (where rN is any ribonucleotide base and G can be either ribo- or deoxyribo-G). The pH profile and reaction products of this deoxyribozyme are similar to those reported for hammerhead ribozyme. This deoxyribozyme has higher activity in the presence of transition metal ions compared to alkaline earth metal ions. At saturating concentrations of Zn²⁺, the cleavage rate is 1.35 min^–1 at pH 6.0; based on pH profile this rate is estimated to be at least ~30 times faster at pH 7.5, where most assays of Mg²⁺-dependent DNA and RNA enzymes are carried out. This work represents a comprehensive characterization of a nucleic acid-based endonuclease that prefers transition metal ions to alkaline earth metal ions. The results demonstrate that nucleic acid enzymes are capable of binding transition metal ions such as Zn²⁺ with high affinity, and the resulting enzymes are more efficient at RNA cleavage than most Mg²⁺-dependent nucleic acid enzymes under similar conditions.     

Reduced Model Captures Mg2+-RNA Interaction Free Energy of Riboswitches

The stability of RNA tertiary structures depends heavily on Mg²⁺. The Mg²⁺-RNA interaction free energy that stabilizes an RNA structure can be computed experimentally through fluorescence-based assays that measure Γ₂₊, the number of excess Mg²⁺ associated with an RNA molecule. Previous explicit-solvent simulations predict that the majority of excess Mg²⁺ ions interact closely and strongly with the RNA, unlike monovalent ions such as K⁺, suggesting that an explicit treatment of Mg²⁺ is important for capturing RNA dynamics. Here we present a reduced model that accurately reproduces the thermodynamics of Mg²⁺-RNA interactions. This model is able to characterize long-timescale RNA dynamics coupled to Mg²⁺ through the explicit representation of Mg²⁺ ions. KCl is described by Debye-Hückel screening and a Manning condensation parameter, which represents condensed K⁺ and models its competition with condensed Mg²⁺. The model contains one fitted parameter, the number of condensed K⁺ ions in the absence of Mg²⁺. Values of Γ₂₊ computed from molecular dynamics simulations using the model show excellent agreement with both experimental data on the adenine riboswitch and previous explicit-solvent simulations of the SAM-I riboswitch. This agreement confirms the thermodynamic accuracy of the model via the direct relation of Γ₂₊ to the Mg²⁺-RNA interaction free energy, and provides further support for the predictions from explicit-solvent calculations. This reduced model will be useful for future studies of the interplay between Mg²⁺ and RNA dynamics.     

Thermodynamic analysis of ligand binding and ligand binding-induced tertiary structure formation by the thiamine pyrophosphate riboswitch

The thi-box riboswitch regulates gene expression in response to the intracellular concentration of thiamine pyrophosphate (TPP) in archaea, bacteria, and eukarya. To complement previous biochemical, genetic, and structural studies of this phylogenetically widespread RNA domain, we have characterized its interaction with TPP by isothermal titration calorimetry. This shows that TPP binding is highly dependent on Mg²⁺ concentration. The dissociation constant decreases from ∼200 nM at 0.5 mM Mg²⁺ concentration to ∼9 nM at 2.5 mM Mg²⁺ concentration. Binding is enthalpically driven, but the unfavorable entropy of binding decreases as Mg²⁺ concentration rises, suggesting that divalent cations serve to pre-organize the RNA. Mutagenesis, biochemical analysis, and a new crystal structure of the riboswitch suggest that a critical element that participates in organizing the riboswitch structure is the tertiary interaction formed between the P3 and L5 regions. This tertiary contact is distant from the TPP binding site, but calorimetric analysis reveals that even subtle mutations in L5 can have readily detectable effects on TPP binding. The thermodynamic signatures of these mutations, namely decreased favorable enthalpy of binding and small effects on entropy of binding, are consistent with the P3–L5 association contributing allosterically to TPP-induced compaction of the RNA.     

Thermodynamic examination of the pyrophosphate sensor helix in the thiamine pyrophosphate riboswitch

Riboswitches are functional mRNA that control gene expression. Thiamine pyrophosphate (TPP) binds to thi-box riboswitch RNA and allosterically inhibits genes that code for proteins involved in the biosynthesis and transport of thiamine. Thiamine binding to the pyrimidine sensor helix and pyrophosphate binding to the pyrophosphate sensor helix cause changes in RNA conformation that regulate gene expression. Here we examine the thermodynamic properties of the internal loop of the pyrophosphate binding domain by comparing the wild-type construct (RNA WT) with six modified 2 × 2 bulged RNA and one 2 × 2 bulged DNA. The wild-type construct retains five conserved bases of the pyrophosphate sensor domain, two of which are in the 2 × 2 bulge (C65 and G66). The RNA WT construct was among the most stable (ΔG°₃₇ = −7.7 kcal/mol) in 1 M KCl at pH 7.5. Breaking the A•G mismatch of the bulge decreases the stability of the construct ∼0.5–1 kcal/mol, but does not affect magnesium binding to the RNA WT. Guanine at position 48 is important for RNA–Mg²⁺ interactions of the TPP-binding riboswitch at pH 7.5. In the presence of 9.5 mM magnesium at pH 5.5, the bulged RNA constructs gained an average of 1.1 kcal/mol relative to 1 M salt. Formation of a single A+•C mismatch base pair contributes about 0.5 kcal/mol at pH 5.5, whereas two tandem A+•C mismatch base pairs together contribute about 2 kcal/mol.     

Reduced Model Captures Mg-RNA Interaction Free Energy of Riboswitches

The stability of RNA tertiary structures depends heavily on Mg²⁺. The Mg²⁺-RNA interaction free energy that stabilizes an RNA structure can be computed experimentally through fluorescence-based assays that measure Γ₂₊, the number of excess Mg²⁺ associated with an RNA molecule. Previous explicit-solvent simulations predict that the majority of excess Mg²⁺ ions interact closely and strongly with the RNA, unlike monovalent ions such as K⁺, suggesting that an explicit treatment of Mg²⁺ is important for capturing RNA dynamics. Here we present a reduced model that accurately reproduces the thermodynamics of Mg²⁺-RNA interactions. This model is able to characterize long-timescale RNA dynamics coupled to Mg²⁺ through the explicit representation of Mg²⁺ ions. KCl is described by Debye-Hückel screening and a Manning condensation parameter, which represents condensed K⁺ and models its competition with condensed Mg²⁺. The model contains one fitted parameter, the number of condensed K⁺ ions in the absence of Mg²⁺. Values of Γ₂₊ computed from molecular dynamics simulations using the model show excellent agreement with both experimental data on the adenine riboswitch and previous explicit-solvent simulations of the SAM-I riboswitch. This agreement confirms the thermodynamic accuracy of the model via the direct relation of Γ₂₊ to the Mg²⁺-RNA interaction free energy, and provides further support for the predictions from explicit-solvent calculations. This reduced model will be useful for future studies of the interplay between Mg²⁺ and RNA dynamics.     

Dissecting the influence of Mg2+ on 3D architecture and ligand-binding of the guanine-sensing riboswitch aptamer domain

Long-range tertiary interactions determine the three-dimensional structure of a number of metabolite-binding riboswitch RNA elements and were found to be important for their regulatory function. For the guanine-sensing riboswitch of the Bacillus subtilis xpt-pbuX operon, our previous NMR-spectroscopic studies indicated pre-formation of long-range tertiary contacts in the ligand-free state of its aptamer domain. Loss of the structural pre-organization in a mutant of this RNA (G37A/C61U) resulted in the requirement of Mg²⁺ for ligand binding. Here, we investigate structural and stability aspects of the wild-type aptamer domain (Gsw) and the G37A/C61U-mutant (Gsw^loop) of the guanine-sensing riboswitch and their Mg²⁺-induced folding characteristics to dissect the role of long-range tertiary interactions, the link between pre-formation of structural elements and ligand-binding properties and the functional stability. Destabilization of the long-range interactions as a result of the introduced mutations for Gsw^loop or the increase in temperature for both Gsw and Gsw^loop involves pronounced alterations of the conformational ensemble characteristics of the ligand-free state of the riboswitch. The increased flexibility of the conformational ensemble can, however, be compensated by Mg²⁺. We propose that reduction of conformational dynamics in remote regions of the riboswitch aptamer domain is the minimal pre-requisite to pre-organize the core region for specific ligand binding.     

Characterization of the metal ion binding site in the anti-terminator protein, HutP, of Bacillus subtilis

HutP is an RNA-binding protein that regulates the expression of the histidine utilization (hut) operon in Bacillus subtilis, by binding to cis-acting regulatory sequences on hut mRNA. It requires L-histidine and an Mg²⁺ ion for binding to the specific sequence within the hut mRNA. In the present study, we show that several divalent cations can mediate the HutP–RNA interactions. The best divalent cations were Mn²⁺, Zn²⁺ and Cd²⁺, followed by Mg²⁺, Co²⁺ and Ni²⁺, while Cu²⁺, Yb²⁺ and Hg²⁺ were ineffective. In the HutP–RNA interactions, divalent cations cannot be replaced by monovalent cations, suggesting that a divalent metal ion is required for mediating the protein–RNA interactions. To clarify their importance, we have crystallized HutP in the presence of three different metal ions (Mg²⁺, Mn²⁺ and Ba²⁺), which revealed the importance of the metal ion binding site. Furthermore, these analyses clearly demonstrated how the metal ions cause the structural rearrangements that are required for the hut mRNA recognition.     

Riboswitch structure: an internal residue mimicking the purine ligand

The adenine and guanine riboswitches regulate gene expression in response to their purine ligand. X-ray structures of the aptamer moiety of these riboswitches are characterized by a compact fold in which the ligand forms a Watson–Crick base pair with residue 65. Phylogenetic analyses revealed a strict restriction at position 39 of the aptamer that prevents the G39–C65 and A39–U65 combinations, and mutational studies indicate that aptamers with these sequence combinations are impaired for ligand binding. In order to investigate the rationale for sequence conservation at residue 39, structural characterization of the U65C mutant from Bacillus subtilis pbuE adenine riboswitch aptamer was undertaken. NMR spectroscopy and X-ray crystallography studies demonstrate that the U65C mutant adopts a compact ligand-free structure, in which G39 occupies the ligand-binding site of purine riboswitch aptamers. These studies present a remarkable example of a mutant RNA aptamer that adopts a native-like fold by means of ligand mimicking and explain why this mutant is impaired for ligand binding. Furthermore, this work provides a specific insight into how the natural sequence has evolved through selection of nucleotide identities that contribute to formation of the ligand-bound state, but ensures that the ligand-free state remains in an active conformation.     

Structure and folding of a rare, natural kink turn in RNA with an A*A pair at the 2b*2n position

The kink turn (k-turn) is a frequently occurring motif, comprising a bulge followed by G•A and A•G pairs that introduces a sharp axial bend in duplex RNA. Natural k-turn sequences exhibit significant departures from the consensus, including the A•G pairs that form critical interactions stabilizing the core of the structure. Kt-23 found in the small ribosomal subunit differs from the consensus in many organisms, particularly in the second A•G pair distal to the bulge (2b•2n). Analysis of many Kt-23 sequences shows that the frequency of occurrence at the 2n position (i.e., on the nonbulged strand, normally G in standard k-turns) is U>C>G>A. Less than 1% of sequences have A at the 2n position, but one such example occurs in Thelohania solenopsae Kt-23. This sequence folds only weakly in the presence of Mg²⁺ ions but is induced to fold normally by the binding of L7Ae protein. Introduction of this sequence into the SAM-I riboswitch resulted in normal binding of SAM ligand, indicating that tertiary RNA contacts have resulted in k-turn folding. X-ray crystallography shows that the T. solenopsae Kt-23 adopts a standard k-turn geometry, making the key, conserved hydrogen bonds in the core and orienting the 1n (of the bulge-proximal A•G pair) and 2b adenine nucleobases in position facing the opposing minor groove. The 2b and 2n adenine nucleobases are not directly hydrogen bonded, but each makes hydrogen bonds to their opposing strands.     

Ion-induced folding of a kink turn that departs from the conventional sequence

Kink turns (k-turns) are important structural motifs that create a sharp axial bend in RNA. Most conform to a consensus in which a three-nucleotide bulge is followed by consecutive G•A and A•G base pairs, and when these G•A pairs are modified in vitro this generally leads to a failure to adopt the k-turn conformation. Kt-23 in the 30S ribosomal subunit of Thermus thermophilus is a rare exception in which the bulge-distal A•G pair is replaced by a non-Watson–Crick A•U pair. In the context of the ribosome, Kt-23 adopts a completely conventional k-turn geometry. We show here that this sequence is induced to fold into a k-turn structure in an isolated RNA duplex by Mg²⁺ or Na⁺ ions. Therefore, the Kt-23 is intrinsically stable despite lacking the key A•G pair; its formation requires neither tertiary interactions nor protein binding. Moreover, the Kt-23 k-turn is stabilized by the same critical hydrogen-bonding interactions within the core of the structure that are found in more conventional sequences such as the near-consensus Kt-7. T. thermophilus Kt-23 has two further non-Watson–Crick base pairs within the non-canonical helix, three and four nucleotides from the bulge, and we find that the nature of these pairs influences the ability of the RNA to adopt k-turn conformation, although the base pair adjacent to the A•U pair is more important than the other.     

Conserved Structural Chemistry for Incision Activity in Structurally Non-homologous Apurinic/Apyrimidinic Endonuclease APE1 and Endonuclease IV DNA Repair Enzymes

Non-coding apurinic/apyrimidinic (AP) sites in DNA form spontaneously and as DNA base excision repair intermediates are the most common toxic and mutagenic in vivo DNA lesion. For repair, AP sites must be processed by 5′ AP endonucleases in initial stages of base repair. Human APE1 and bacterial Nfo represent the two conserved 5′ AP endonuclease families in the biosphere; they both recognize AP sites and incise the phosphodiester backbone 5′ to the lesion, yet they lack similar structures and metal ion requirements. Here, we determined and analyzed crystal structures of a 2.4 Å resolution APE1-DNA product complex with Mg²⁺ and a 0.92 Å Nfo with three metal ions. Structural and biochemical comparisons of these two evolutionarily distinct enzymes characterize key APE1 catalytic residues that are potentially functionally similar to Nfo active site components, as further tested and supported by computational analyses. We observe a magnesium-water cluster in the APE1 active site, with only Glu-96 forming the direct protein coordination to the Mg²⁺. Despite differences in structure and metal requirements of APE1 and Nfo, comparison of their active site structures surprisingly reveals strong geometric conservation of the catalytic reaction, with APE1 catalytic side chains positioned analogously to Nfo metal positions, suggesting surprising functional equivalence between Nfo metal ions and APE1 residues. The finding that APE1 residues are positioned to substitute for Nfo metal ions is supported by the impact of mutations on activity. Collectively, the results illuminate the activities of residues, metal ions, and active site features for abasic site endonucleases.     

Binding of C5 Protein to P RNA Enhances the Rate Constant for Catalysis for P RNA Processing of Pre-tRNAs Lacking a Consensus G(+ 1)/C(+ 72) Pair

The RNA subunit of the ribonucleoprotein enzyme ribonuclease P (RNase P (P RNA) contains the active site, but binding of Escherichia coli RNase P protein (C5) to P RNA increases the rate constant for catalysis for certain pre-tRNA substrates up to 1000-fold. Structure-swapping experiments between a substrate that is cleaved slowly by P RNA alone (pre-tRNA^f-met605) and one that is cleaved quickly (pre-tRNA^met608) pinpoint the characteristic C(+ 1)/A(+ 72) base pair of initiator tRNA^f-met as the sole determinant of slow RNA-alone catalysis. Unlike other substrate modifications that slow RNA-alone catalysis, the presence of a C(+ 1)/A(+ 72) base pair reduces the rate constant for processing at both correct and miscleavage sites, indicating an indirect but nonetheless important role in catalysis. Analysis of the Mg²⁺ dependence of apparent catalytic rate constants for pre-tRNA^met608 and a pre-tRNA^met608 (+ 1)C/(+ 72)A mutant provides evidence that C5 promotes rate enhancement primarily by compensating for the decrease in the affinity of metal ions important for catalysis engendered by the presence of the CA pair. Together, these results support and extend current models for RNase P substrate recognition in which contacts involving the conserved (+ 1)G/C(+ 72) pair of tRNA stabilize functional metal ion binding. Additionally, these observations suggest that C5 protein has evolved to compensate for tRNA variation at positions important for binding to P RNA, allowing for tRNA specialization.     

Biochemical analysis of hatchet self-cleaving ribozymes

Hatchet RNAs are members of a novel self-cleaving ribozyme class that was recently discovered by using a bioinformatics search strategy. The consensus sequence and secondary structure of this class includes 13 highly conserved and numerous other modestly conserved nucleotides interspersed among bulges linking four base-paired substructures. A representative hatchet ribozyme from a metagenomic source requires divalent ions such as Mg²⁺ to promote RNA strand scission with a maximum rate constant of ∼4 min⁻¹. As with all other small self-cleaving ribozymes discovered to date, hatchet ribozymes employ a general mechanism for catalysis involving the nucleophilic attack of a ribose 2′-oxygen atom on an adjacent phosphorus center. Kinetic characteristics of the reaction demonstrate that members of this ribozyme class have an essential requirement for divalent metal ions and that they might have a complex active site that employs multiple catalytic strategies to accelerate RNA cleavage by internal phosphoester transfer.