Adenine protonation in domain B of the hairpin ribozyme

A stable uncompleted tetramolecular G-quadruplex formed by d(AGnA) under acidic condition

In this work, the effects of terminal adenines on the formation and stability of tetramolecular G-quadruplexes (G4s) have been studied by electrospray ionization mass spectrometry (ESI-MS), UV, CD and NMR spectroscopy. Several evidences suggested that the sequences d(AG_nA) (n = 4 or 5) form stable uncompleted tetramolecular G4 at acidic condition which is different from the canonical one in the neutral condition. In addition, hydrolysis of guanine has also been observed in acidic condition that may occur for unpaired strands rather than in complete G4. Thus, a new G4 topology containing incomplete G-quartet is proposed that is very stable and particularly presents in acidic ammonium ions solution. The information presented in this study provides the new insight on the polymorphism of G4s in acidic environment, which may help understand of the special role of adenines on the formation of G4s.

RNA Dynamics by NMR Spectroscopy

An ever-increasing number of functional RNAs require a mechanistic understanding. RNA function relies on changes in its structure, so-called dynamics. To reveal dynamic processes and higher energy structures, new NMR methods have been developed to elucidate these dynamics in RNA with atomic resolution. In this Review, we provide an introduction to dynamics novices and an overview of methods that access most dynamic timescales, from picoseconds to hours. Examples are provided as well as insight into theory, data acquisition and analysis for these different methods. Using this broad spectrum of methodology, unprecedented detail and invisible structures have been obtained and are reviewed here. RNA, though often more complicated and therefore neglected, also provides a great system to study structural changes, as these RNA structural changes are more easily defined-Lego like-than in proteins, hence the numerous revelations of RNA excited states.

The effect of adenine protonation on RNA phosphodiester backbone bond cleavage elucidated by deaza-nucleobase modifications and mass spectrometry

The catalytic strategies of small self-cleaving ribozymes often involve interactions between nucleobases and the ribonucleic acid (RNA) backbone. Here we show that multiply protonated, gaseous RNA has an intrinsic preference for the formation of ionic hydrogen bonds between adenine protonated at N3 and the phosphodiester backbone moiety on its 5'-side that facilitates preferential phosphodiester backbone bond cleavage upon vibrational excitation by low-energy collisionally activated dissociation. Removal of the basic N3 site by deaza-modification of adenine was found to abrogate preferential phosphodiester backbone bond cleavage. No such effects were observed for N1 or N7 of adenine. Importantly, we found that the pH of the solution used for generation of the multiply protonated, gaseous RNA ions by electrospray ionization affects phosphodiester backbone bond cleavage next to adenine, which implies that the protonation patterns in solution are at least in part preserved during and after transfer into the gas phase. Our study suggests that interactions between protonated adenine and phosphodiester moieties of RNA may play a more important mechanistic role in biological processes than considered until now.

Base-pair Opening Dynamics of Nucleic Acids in Relation to Their Biological Function

Base-pair opening is a conformational transition that is required for proper biological function of nucleic acids. Hydrogen exchange, observed by NMR spectroscopic experiments, is a widely used method to study the thermodynamics and kinetics of base-pair opening in nucleic acids. The hydrogen exchange data of imino protons are analyzed based on a two-state (open/closed) model for the base-pair, where hydrogen exchange only occurs from the open state. In this review, we discuss examples of how hydrogen exchange data provide insight into several interesting biological processes involving functional interactions of nucleic acids: i) selective recognition of DNA by proteins; ii) regulation of RNA cleavage by site-specific mutations; iii) intermolecular interaction of proteins with their target DNA or RNA; iv) formation of PNA:DNA hybrid duplexes.

•

This review systematically summarizes hydrogen exchange theory for base-paired imino protons of nucleic acids.

•

Base-pair opening kinetics explain how the DNA can be selectively recognized by its target proteins.

•

Base-pair opening kinetics explain the mechanisms by which site-specific mutations regulate RNA cleavage.

•

Hydrogen exchange studies can elucidate the intermolecular interaction of proteins with their target DNA or RNA.

Coarse-grained dynamic RNA titration simulations

Electrostatic interactions play a pivotal role in many biomolecular processes. The molecular organization and function in biological systems are largely determined by these interactions. Owing to the highly negative charge of RNA, the effect is expected to be more pronounced in this system. Moreover, RNA base pairing is dependent on the charge of the base, giving rise to alternative secondary and tertiary structures. The equilibrium between uncharged and charged bases is regulated by the solution pH, which is therefore a key environmental condition influencing the molecule's structure and behaviour. By means of constant-pH Monte Carlo simulations based on a fast proton titration scheme, coupled with the coarse-grained model HiRE-RNA, molecular dynamic simulations of RNA molecules at constant pH enable us to explore the RNA conformational plasticity at different pH values as well as to compute electrostatic properties as local pK _a values for each nucleotide.

Single-Molecule Mechanical Folding and Unfolding of RNA Hairpins: Effects of Single A-U to A·C Pair Substitutions and Single Proton Binding and Implications for mRNA Structure-Induced -1 Ribosomal Frameshifting

A wobble A·C pair can be protonated at near physiological pH to form a more stable wobble A⁺·C pair. Here, we constructed an RNA hairpin (rHP) and three mutants with one A-U base pair substituted with an A·C mismatch on the top (near the loop, U22C), middle (U25C), and bottom (U29C) positions of the stem, respectively. Our results on single-molecule mechanical (un)folding using optical tweezers reveal the destabilization effect of A-U to A·C pair substitution and protonation-dependent enhancement of mechanical stability facilitated through an increased folding rate, or decreased unfolding rate, or both. Our data show that protonation may occur rapidly upon the formation of an apparent mechanical folding transition state. Furthermore, we measured the bulk -1 ribosomal frameshifting efficiencies of the hairpins by a cell-free translation assay. For the mRNA hairpins studied, -1 frameshifting efficiency correlates with mechanical unfolding force at equilibrium and folding rate at around 15 pN. U29C has a frameshifting efficiency similar to that of rHP (∼2%). Accordingly, the bottom 2-4 base pairs of U29C may not form under a stretching force at pH 7.3, which is consistent with the fact that the bottom base pairs of the hairpins may be disrupted by ribosome at the slippery site. U22C and U25C have a similar frameshifting efficiency (∼1%), indicating that both unfolding and folding rates of an mRNA hairpin in a crowded environment may affect frameshifting. Our data indicate that mechanical (un)folding of RNA hairpins may mimic how mRNAs unfold and fold in the presence of translating ribosomes.

Crystal structures of N6-modified-aminoacid/peptide nucleobase analogs: hybrid adenine–glycine and adenine–glycylglycine molecules

Anion–π interactions in crystal structures of N⁶-modified-aminoacid and dipeptide adenine analogs are investigated using X-ray crystallography and DFT calculations.

The crystal structure of the 5΄ functional domain of the transcription riboregulator 7SK

In vertebrates, the 7SK RNA forms the scaffold of a complex, which regulates transcription pausing of RNA-polymerase II. By binding to the HEXIM protein, the complex comprising proteins LARP7 and MePCE captures the positive transcription elongation factor P-TEFb and prevents phosphorylation of pausing factors. The HEXIM-binding site embedded in the 5΄-hairpin of 7SK (HP1) encompasses a short signature sequence, a GAUC repeat framed by single-stranded uridines. The present crystal structure of HP1 shows a remarkably straight helical stack involving several unexpected triples formed at a central region. Surprisingly, two uridines of the signature sequence make triple interactions in the major groove of the (GAUC)₂. The third uridine is turned outwards or inward, wedging between the other uridines, thus filling the major groove. A molecular dynamics simulation indicates that these two conformations of the signature sequence represent stable alternatives. Analyses of the interaction with the HEXIM protein confirm the importance of the triple interactions at the signature sequence. Altogether, the present structural analysis of 7SK HP1 highlights an original mechanism of swapping bases, which could represent a possible ‘7SK signature’ and provides new insight into the functional importance of the plasticity of RNA.

The structure of an E. coli tRNAfMet A1-U72 variant shows an unusual conformation of the A1-U72 base pair

Translation initiation in eukaryotes and archaea involves a methionylated initiator tRNA delivered to the ribosome in a ternary complex with e/aIF2 and GTP. Eukaryotic and archaeal initiator tRNAs contain a highly conserved A₁-U₇₂ base pair at the top of the acceptor stem. The importance of this base pair to discriminate initiator tRNAs from elongator tRNAs has been established previously using genetics and biochemistry. However, no structural data illustrating how the A₁-U₇₂ base pair participates in the accurate selection of the initiator tRNAs by the translation initiation systems are available. Here, we describe the crystal structure of a mutant E. coli initiator tRNA_f^MetA₁-U₇₂, aminoacylated with methionine, in which the C₁:A₇₂ mismatch at the end of the tRNA acceptor stem has been changed to an A₁-U₇₂ base pair. Sequence alignments show that the mutant E. coli tRNA is a good mimic of archaeal initiator tRNAs. The crystal structure, determined at 2.8 Å resolution, shows that the A₁-U₇₂ pair adopts an unusual arrangement. A₁ is in a syn conformation and forms a single H-bond interaction with U₇₂ This interaction requires protonation of the N1 atom of A₁ Moreover, the 5' phosphoryl group folds back into the major groove of the acceptor stem and interacts with the N7 atom of G₂ A possible role of this unusual geometry of the A₁-U₇₂ pair in the recognition of the initiator tRNA by its partners during eukaryotic and archaeal translation initiation is discussed.

The Novel Chemical Mechanism of the Twister Ribozyme

We describe the multifactorial origins of catalysis by the twister ribozyme. We provide evidence that the adenine immediately 3' to the scissile phosphate (A1) acts as a general acid. Substitution of ring nitrogen atoms indicates that very unusually the N3 of A1 is the proton donor to the oxyanion leaving group. A1 is accommodated in a specific binding pocket that raises its pKa toward neutrality, juxtaposes its N3 with the O5' to be protonated, and helps create the in-line trajectory required for nucleophilic attack. A1 performs general acid catalysis while G33 acts as a general base. A 100-fold stereospecific phosphorothioate effect at the scissile phosphate is consistent with a significant stabilization of the transition state by the ribozyme, and functional group substitution at G33 indicates that its exocyclic N2 interacts directly with the scissile phosphate. A model of the ribozyme active site is proposed that accommodates these catalytic strategies.

A Mini-Twister Variant and Impact of Residues/Cations on the Phosphodiester Cleavage of this Ribozyme Class

Nucleolytic ribozymes catalyze site-specific cleavage of their phosphodiester backbones. A minimal version of the twister ribozyme is reported that lacks the phylogenetically conserved stem P1 while retaining wild-type activity. Atomic mutagenesis revealed that nitrogen atoms N1 and N3 of the adenine-6 at the cleavage site are indispensable for cleavage. By NMR spectroscopy, a pKa value of 5.1 was determined for a (13) C2-labeled adenine at this position in the twister ribozyme, which is significantly shifted compared to the pKa of the same adenine in the substrate alone. This finding pinpoints at a potential role for adenine-6 in the catalytic mechanism besides the previously identified invariant guanine-48 and a Mg(2+) ion, both of which are directly coordinated to the non-bridging oxygen atoms of the scissile phosphate; for the latter, additional evidence stems from the observation that Mn(2+) or Cd(2+) accelerated cleavage of phosphorothioate substrates. The relevance of this metal ion binding site is further emphasized by a new 2.6 Å X-ray structure of a 2'-OCH3 -U5 modified twister ribozyme.

Centromeric Alpha-Satellite DNA Adopts Dimeric i-Motif Structures Capped by AT Hoogsteen Base Pairs

Human centromeric alpha-satellite DNA is composed of tandem arrays of two types of 171 bp monomers; type A and type B. The differences between these types are concentrated in a 17 bp region of the monomer called the A/B box. Here, we have determined the solution structure of the C-rich strand of the two main variants of the human alpha-satellite A box. We show that, under acidic conditions, the C-rich strands of two A boxes self-recognize and form a head-to-tail dimeric i-motif stabilized by four intercalated hemi-protonated C:C(+) base pairs. Interestingly, the stack of C:C(+) base pairs is capped by T:T and Hoogsteen A:T base pairs. The two main variants of the A box adopt a similar three-dimensional structure, although the residues involved in the formation of the i-motif core are different in each case. Together with previous studies showing that the B box (known as the CENP-B box) also forms dimeric i-motif structures, our finding of this non-canonical structure in the A box shows that centromeric alpha satellites in all human chromosomes are able to form i-motifs, which consequently raises the possibility that these structures may play a role in the structural organization of the centromere.

Role of tautomerism in RNA biochemistry

Heterocyclic nucleic acid bases and their analogs can adopt multiple tautomeric forms due to the presence of multiple solvent-exchangeable protons. In DNA, spontaneous formation of minor tautomers has been speculated to contribute to mutagenic mispairings during DNA replication, whereas in RNA, minor tautomeric forms have been proposed to enhance the structural and functional diversity of RNA enzymes and aptamers. This review summarizes the role of tautomerism in RNA biochemistry, specifically focusing on the role of tautomerism in catalysis of small self-cleaving ribozymes and recognition of ligand analogs by riboswitches. Considering that the presence of multiple tautomers of nucleic acid bases is a rare occurrence, and that tautomers typically interconvert on a fast time scale, methods for studying rapid tautomerism in the context of nucleic acids under biologically relevant aqueous conditions are also discussed.

Building a stable RNA U-turn with a protonated cytidine

The U-turn is a classical three-dimensional RNA folding motif first identified in the anticodon and T-loops of tRNAs. It also occurs frequently as a building block in other functional RNA structures in many different sequence and structural contexts. U-turns induce sharp changes in the direction of the RNA backbone and often conform to the 3-nt consensus sequence 5'-UNR-3' (N = any nucleotide, R = purine). The canonical U-turn motif is stabilized by a hydrogen bond between the N3 imino group of the U residue and the 3' phosphate group of the R residue as well as a hydrogen bond between the 2'-hydroxyl group of the uridine and the N7 nitrogen of the R residue. Here, we demonstrate that a protonated cytidine can functionally and structurally replace the uridine at the first position of the canonical U-turn motif in the apical loop of the neomycin riboswitch. Using NMR spectroscopy, we directly show that the N3 imino group of the protonated cytidine forms a hydrogen bond with the backbone phosphate 3' from the third nucleotide of the U-turn analogously to the imino group of the uridine in the canonical motif. In addition, we compare the stability of the hydrogen bonds in the mutant U-turn motif to the wild type and describe the NMR signature of the C+-phosphate interaction. Our results have implications for the prediction of RNA structural motifs and suggest simple approaches for the experimental identification of hydrogen bonds between protonated C-imino groups and the phosphate backbone.

Constant pH molecular dynamics of proteins in explicit solvent with proton tautomerism

pH is a ubiquitous regulator of biological activity, including protein-folding, protein-protein interactions, and enzymatic activity. Existing constant pH molecular dynamics (CPHMD) models that were developed to address questions related to the pH-dependent properties of proteins are largely based on implicit solvent models. However, implicit solvent models are known to underestimate the desolvation energy of buried charged residues, increasing the error associated with predictions that involve internal ionizable residue that are important in processes like hydrogen transport and electron transfer. Furthermore, discrete water and ions cannot be modeled in implicit solvent, which are important in systems like membrane proteins and ion channels. We report on an explicit solvent constant pH molecular dynamics framework based on multi-site λ-dynamics (CPHMD(MSλD)). In the CPHMD(MSλD) framework, we performed seamless alchemical transitions between protonation and tautomeric states using multi-site λ-dynamics, and designed novel biasing potentials to ensure that the physical end-states are predominantly sampled. We show that explicit solvent CPHMD(MSλD) simulations model realistic pH-dependent properties of proteins such as the Hen-Egg White Lysozyme (HEWL), binding domain of 2-oxoglutarate dehydrogenase (BBL) and N-terminal domain of ribosomal protein L9 (NTL9), and the pKa predictions are in excellent agreement with experimental values, with a RMSE ranging from 0.72 to 0.84 pKa units. With the recent development of the explicit solvent CPHMD(MSλD) framework for nucleic acids, accurate modeling of pH-dependent properties of both major class of biomolecules-proteins and nucleic acids is now possible.

Photophysical and calorimetric studies on the binding of 9-O-substituted analogs of the plant alkaloid berberine to double stranded poly(A)

This interaction of four novel 9-O-substituted analogs of the plant alkaloid berberine with double stranded poly(A) was studied using a variety of biophysical techniques. Remarkably higher binding of two 9-O-ω-amino alkyl ether analogs compared to the two 9-O-N-aryl/arylalkyl amino carbonyl methyl berberine analogs was observed. Quantum efficiency values suggested that energy was transferred from the adenine base pairs to the analogs on binding. Ferrocyanide quenching and viscosity studies revealed the binding mode to be intercalative for these analogs. Circular dichroism studies showed that these analogs induced significant conformational changes in the secondary structure of ds poly(A). Energetics of the binding suggested that 9-O-N-aryl/arylalkyl amino carbonyl methyl berberines bound very weakly to ds poly(A). The binding of 9-O-ω-amino alkyl ether analogs was entropy dominated with a smaller but favorable enthalpic contribution to the Gibbs energy. Increasing the temperature resulted in weaker binding; the enthalpic contribution increased and the entropic contribution decreased. A small negative heat capacity change with significant enthalpy-entropy compensation established the involvement of multiple weak noncovalent interactions in the binding process.

Towards Accurate Prediction of Protonation Equilibrium of Nucleic Acids

The role of protonated nucleotides in modulating the pH-dependent properties of nucleic acids is one of the emerging frontiers in the field of nucleic acid biology. The recent development of a constant pH molecular dynamics simulation (CPHMD^MSλD) framework for simulating nucleic acids has provided a tool for realistic simulations of pH-dependent dynamics. We enhanced the CPHMD^MSλD framework with pH-based replica exchange (pH-REX), which significantly improves the sampling of both titration and spatial coordinates. The results from our pK_a calculations for the GAAA tetraloop, which was predicted with lower accuracy previously due to sampling challenges, demonstrates that pH-REX reduces the average unsigned error (AUE) to 0.7 pK_a units, and the error of the most poorly predicted residue A17 was drastically reduced from 2.9 to 1.2 pK_a unit. Lastly, we show that pH-REX CPHMD^MSλD simulations can be used to identify the dominant conformation of nucleic acid structures in alternate pH environments. This work suggests that pH-REX CPHMD^MSλD simulations provide a practical tool for predicting nucleic acid protonation equilibrium from first-principles, and offering structural and mechanistic insight into the study of pH-dependent properties of nucleic acids.

pH-dependent dynamics of complex RNA macromolecules

The role of pH-dependent protonation equilibrium in modulating RNA dynamics and function is one of the key unanswered questions in RNA biology. Molecular dynamics (MD) simulations can provide insight into the mechanistic roles of protonated nucleotides, but it is only capable of modeling fixed protonation states and requires prior knowledge of the key residue's protonation state. Recently, we developed a framework for constant pH molecular dynamics simulations (CPHMDMSλD) of nucleic acids, where the nucleotides' protonation states are modeled as dynamic variables that are coupled to the structural dynamics of the RNA. In the present study, we demonstrate the application of CPHMDMSλD to the lead-dependent ribozyme; establishing the validity of this approach for modeling complex RNA structures. We show that CPHMDMSλD accurately predicts the direction of the pKa shifts and reproduces experimentally-measured microscopic pKa values with an average unsigned error of 1.3 pKa units. The effects of coupled titration states in RNA structures are modeled, and the importance of conformation sampling is highlighted. The general accuracy of CPHMDMSλD simulations in reproducing pH-dependent observables reported in this work demonstrates that constant pH simulations provides a powerful tool to investigate pH-dependent processes in nucleic acids.

Mechanisms of RNA catalysis

Ribozymes are RNA molecules that act as chemical catalysts. In contemporary cells, most known ribozymes carry out phosphoryl transfer reactions. The nucleolytic ribozymes comprise a class of five structurally-distinct species that bring about site-specific cleavage by nucleophilic attack of the 2'-O on the adjacent 3'-P to form a cyclic 2',3'-phosphate. In general, they will also catalyse the reverse reaction. As a class, all these ribozymes appear to use general acid-base catalysis to accelerate these reactions by about a million-fold. In the Varkud satellite ribozyme, we have shown that the cleavage reaction is catalysed by guanine and adenine nucleobases acting as general base and acid, respectively. The hairpin ribozyme most probably uses a closely similar mechanism. Guanine nucleobases appear to be a common choice of general base, but the general acid is more variable. By contrast, the larger ribozymes such as the self-splicing introns and RNase P act as metalloenzymes.

Constant pH Molecular Dynamics Simulations of Nucleic Acids in Explicit Solvent

The nucleosides of adenine and cytosine have pKa values of 3.50 and 4.08, respectively, and are assumed to be unprotonated under physiological conditions. However, evidence from recent NMR and X-Ray crystallography studies has revealed the prevalence of protonated adenine and cytosine in RNA macromolecules. Such nucleotides with elevated pKa values may play a role in stabilizing RNA structure and participate in the mechanism of ribozyme catalysis. With the work presented here, we establish the framework and demonstrate the first constant pH MD simulations (CPHMD) for nucleic acids in explicit solvent in which the protonation state is coupled to the dynamical evolution of the RNA system via λ-dynamics. We adopt the new functional form λ(Nexp) for λ that was recently developed for Multi-Site λ-Dynamics (MSλD) and demonstrate good sampling characteristics in which rapid and frequent transitions between the protonated and unprotonated states at pH = pKa are achieved. Our calculated pKa values of simple nucleotides are in a good agreement with experimentally measured values, with a mean absolute error of 0.24 pKa units. This work demonstrates that CPHMD can be used as a powerful tool to investigate pH-dependent biological properties of RNA macromolecules.

Solution structure of the HIV-1 exon splicing silencer 3

Alternative splicing of the human immunodeficiency virus type 1 (HIV-1) genomic RNA is necessary to produce the complete viral protein complement, and aberrations in the splicing pattern impair HIV-1 replication. Genome splicing in HIV-1 is tightly regulated by the dynamic assembly/disassembly of trans host factors with cis RNA control elements. The host protein, heterogeneous nuclear ribonucleoprotein (hnRNP) A1, regulates splicing at several highly conserved HIV-1 3' splice sites by binding 5'-UAG-3' elements embedded within regions containing RNA structure. The physical determinants of hnRNP A1 splice site recognition remain poorly defined in HIV-1, thus precluding a detailed understanding of the molecular basis of the splicing pattern. Here, the three-dimensional structure of the exon splicing silencer 3 (ESS3) from HIV-1 has been determined using NMR spectroscopy. ESS3 adopts a 27-nucleotide hairpin with a 10-bp A-form stem that contains a pH-sensitive A(+)C wobble pair. The seven-nucleotide hairpin loop contains the high-affinity hnRNP-A1-responsive 5'-UAGU-3' element and a proximal 5'-GAU-3' motif. The NMR structure shows that the heptaloop adopts a well-organized conformation stabilized primarily by base stacking interactions reminiscent of a U-turn. The apex of the loop is quasi-symmetric with UA dinucleotide steps from the 5'-GAU-3' and 5'-UAGU-3' motifs stacking on opposite sides of the hairpin. As a step towards understanding the binding mechanism, we performed calorimetric and NMR titrations of several hnRNP A1 subdomains into ESS3. The data show that the UP1 domain forms a high-affinity (K(d)=37.8±1.1 nM) complex with ESS3 via site-specific interactions with the loop.

NMR structure of a 4 x 4 nucleotide RNA internal loop from an R2 retrotransposon: identification of a three purine-purine sheared pair motif and comparison to MC-SYM predictions

The NMR solution structure is reported of a duplex, 5'GUGAAGCCCGU/3'UCACAGGAGGC, containing a 4 × 4 nucleotide internal loop from an R2 retrotransposon RNA. The loop contains three sheared purine-purine pairs and reveals a structural element found in other RNAs, which we refer to as the 3RRs motif. Optical melting measurements of the thermodynamics of the duplex indicate that the internal loop is 1.6 kcal/mol more stable at 37°C than predicted. The results identify the 3RRs motif as a common structural element that can facilitate prediction of 3D structure. Known examples include internal loops having the pairings: 5'GAA/3'AGG, 5'GAG/3'AGG, 5'GAA/3'AAG, and 5'AAG/3'AGG. The structural information is compared with predictions made with the MC-Sym program.

glmS Riboswitch binding to the glucosamine-6-phosphate α-anomer shifts the pKa toward neutrality

The glmS riboswitch regulates gene expression through a self-cleavage activity. The reaction is catalyzed with the assistance of the metabolite cofactor glucosamine-6-phosphate (GlcN6P), whose amino group is proposed to serve as the general acid during the reaction. This reaction is pH-dependent with a pK(a) that is lower than the observed pK(a) for the amine of GlcN6P in solution. GlcN6P, like other pyranose sugars, undergoes spontaneous and rapid interconversion between the α and β anomers at the C1 position. Here we demonstrate by NMR that the Bacillus anthracis glmS riboswitch selectively binds the α-anomer of GlcN6P with a maximum binding affinity of 0.36 mM and that binding is pH-dependent. We also report that the anomeric ratio between α and β is pH-dependent and the pK(a)s of the two amines differ by 0.5 pH units, α being the higher of the two (pK(a)=8.3). The pH dependence of binding reveals a pK(a) of 6.7, suggesting that the glmS RNA reduces the pK(a) of the GlcN6P amine by 1.6 units in the ground state. We reevaluated previously obtained kinetic data and found the reaction pK(a) is 6.9, within error of the binding data. The data support a model where the reaction pK(a) corresponds to that of the GlcN6P amine. This observation has broader relevance for considering how the microenvironment of an RNA, despite its anionic character, can reduce the pK(a)s of functional groups for use in catalysis.

The pH dependence of hairpin ribozyme catalysis reflects ionization of an active site adenine

Understanding how self-cleaving ribozymes mediate catalysis is crucial in light of compelling evidence that human and bacterial gene expression can be regulated through RNA self-cleavage. The hairpin ribozyme catalyzes reversible phosphodiester bond cleavage through a mechanism that does not require divalent metal cations. Previous structural and biochemical evidence implicated the amidine group of an active site adenosine, A38, in a pH-dependent step in catalysis. We developed a way to determine microscopic pK(a) values in active ribozymes based on the pH-dependent fluorescence of 8-azaadenosine (8azaA). We compared the microscopic pK(a) for ionization of 8azaA at position 38 with the apparent pK(a) for the self-cleavage reaction in a fully functional hairpin ribozyme with a unique 8azaA at position 38. Microscopic and apparent pK(a) values were virtually the same, evidence that A38 protonation accounts for the decrease in catalytic activity with decreasing pH. These results implicate the neutral unprotonated form of A38 in a transition state that involves formation of the 5'-oxygen-phosphorus bond.

NMR structure of the A730 loop of the Neurospora VS ribozyme: insights into the formation of the active site

The Neurospora VS ribozyme is a small nucleolytic ribozyme with unique primary, secondary and global tertiary structures, which displays mechanistic similarities to the hairpin ribozyme. Here, we determined the high-resolution NMR structure of a stem–loop VI fragment containing the A730 internal loop, which forms part of the active site. In the presence of magnesium ions, the A730 loop adopts a structure that is consistent with existing biochemical data and most likely reflects its conformation in the VS ribozyme prior to docking with the cleavage site internal loop. Interestingly, the A730 loop adopts an S-turn motif that is also present in loop B within the hairpin ribozyme active site. The S-turn appears necessary to expose the Watson–Crick edge of a catalytically important residue (A756) so that it can fulfill its role in catalysis. The A730 loop and the cleavage site loop of the VS ribozyme display structural similarities to internal loops found in the active site of the hairpin ribozyme. These similarities provided a rationale to build a model of the VS ribozyme active site based on the crystal structure of the hairpin ribozyme.

Direct measurement of the ionization state of an essential guanine in the hairpin ribozyme

Active site guanines are critical for self-cleavage reactions of several ribozymes, but their precise functions in catalysis are unclear. To learn whether protonated or deprotonated forms of guanine predominate in the active site, microscopic pKa values were determined for ionization of 8-azaguanosine substituted for G8 in the active site of a fully functional hairpin ribozyme in order to determine microscopic pKa values for 8-azaguanine deprotonation from the pH dependence of fluorescence. Microscopic pKa values above 9 for deprotonation of 8-azaguanine in the active site were about 3 units higher than apparent pKa values determined from the pH dependence of self-cleavage kinetics. Thus, the increase in activity with increasing pH does not correlate with deprotonation of G8, and most of G8 is protonated at neutral pH. These results do not exclude a role in proton transfer, but a simple interpretation is that G8 functions in the protonated form, perhaps by donating hydrogen bonds.

Fluorine substituted adenosines as probes of nucleobase protonation in functional RNAs

Ionized nucleobases are required for folding, conformational switching, or catalysis in a number of functional RNAs. A common strategy to study these sites employs nucleoside analogues with perturbed pKa, but the interpretation of these studies is often complicated by the chemical modification introduced, in particular modifications that add, remove, or translocate hydrogen bonding groups in addition to perturbing pKa values. In the present study we present a series of fluorine substituted adenosine analogues that produce large changes in N1 pKa values with minimal structural perturbation. These analogues include fluorine for hydrogen substitutions in the adenine ring of adenosine and 7-deaza-adenosine with resulting N1 pKa values spanning more than 4 pKa units. To demonstrate the utility of these analogues we have conducted a nucleotide analogue interference mapping (NAIM) study on a self-ligating construct of the Varkud Satellite (VS) ribozyme. We find that each of the analogues is readily incorporated by T7 RNA polymerase and produces fully active transcripts when substituted at the majority of sites. Strong interferences are observed for three sites known to be critical for VS ribozyme function, most notably A756. Substitutions at A756 lead to slight enhancements in activity for elevated pKa analogues and dramatic interferences in activity for reduced pKa analogues, supporting the proposed catalytic role for this base. The structural similarity of these analogues, combined with their even incorporation and selective interference, provides an improved method for identifying sites of adenosine protonation in a variety of systems.

A chemo-genetic approach for the study of nucleobase participation in nucleolytic ribozymes

A novel chemo-genetic approach for the analysis of general acid-base catalysis by nucleobases in ribozymes is reviewed. This involves substitution of a C-nucleoside with imidazole in place of a natural nucleobase. The Varkud satellite ribozyme in which the nucleobase at the critical 756 position has been replaced by imidazole is active in both cleavage and ligation reactions. Similarly, a modified hairpin ribozyme with the nucleobase at position 8 substituted by imidazole is active in cleavage and ligation reactions. Although the rates are lower than those of the natural ribozymes, they are significantly greater than other variants at these positions. The dependence of the hairpin ribozyme reaction rates on pH has been studied. Both cleavage and ligation reactions display a bell-shaped pH dependence, consistent with general acid-base catalysis involving the nucleotide at position 8.

A guanine nucleobase important for catalysis by the VS ribozyme

A guanine (G638) within the substrate loop of the VS ribozyme plays a critical role in the cleavage reaction. Replacement by any other nucleotide results in severe impairment of cleavage, yet folding of the substrate is not perturbed, and the variant substrates bind the ribozyme with similar affinity, acting as competitive inhibitors. Functional group substitution shows that the imino proton on the N1 is critical, suggesting a possible role in general acid-base catalysis, and this in accord with the pH dependence of the reaction rate for the natural and modified substrates. We propose a chemical mechanism for the ribozyme that involves general acid-base catalysis by the combination of the nucleobases of guanine 638 and adenine 756. This is closely similar to the probable mechanism of the hairpin ribozyme, and the active site arrangements for the two ribozymes appear topologically equivalent. This has probably arisen by convergent evolution.

8-Azaguanine reporter of purine ionization states in structured RNAs

The fluorescent nucleotide analogue 8-azaguanosine-5'-triphosphate (8azaGTP) is prepared easily by in vitro enzymatic synthesis methods. 8azaGTP is an efficient substrate for T7 RNA polymerase and is incorporated specifically opposite cytosine in the transcription template, as expected for a nucleobase analogue with the same Watson-Crick hydrogen bonding face as guanine. 8-Azaguanine (8azaG) in oligonucleotides also is recognized as guanine during ribonuclease T1 digestion. Moreover, replacement of guanine by 8azaG does not alter the melting temperature of base-paired RNAs significantly, evidence that 8azaG does not disrupt stacking and hydrogen bonding interactions. 8azaGTP displays a high fluorescent quantum yield when the N1 position is deprotonated at high pH, but fluorescence intensity decreases significantly when N1 is protonated at neutral pH. Fluorescence is quenched 10-fold to 100-fold when 8azaG is incorporated into base-paired RNA and remains pH-dependent, although apparent pKa values determined from the pH dependence of fluorescence intensity shift in the basic direction. Thus, 8azaG is a guanine analogue that does not perturb RNA structure and displays pH-dependent fluorescence that can be used to probe the ionization states of nucleobases in structured RNAs. A key application will be in determining the ionization state of active site nucleobases that have been implicated in the catalytic mechanisms of RNA enzymes.

Active Site Labeling of G8 in the Hairpin Ribozyme: Implications for Structure and Mechanism

There is mounting evidence that suggests that general acid/base catalysis is operative in the hairpin ribozyme, with analogy to the protein enzyme RNaseA. Nevertheless, the extent of general base catalysis as well as the identity of the specific chemical groups responsible remains the subject of some controversy. An affinity label has previously been used to alkylate histidine 12 (His12), the active general base in RNaseA. To date, no such experiment has been applied to a ribozyme. We have synthesized the analogous affinity label for the hairpin ribozyme with an electrophilic 2'-bromoacetamide group in lieu of the 2'-hydroxyl (2'OH) at the substrate cleavage site and show that guanosine 8 (G8) of the hairpin ribozyme is specifically alkylated, most likely at the N1 position. This evidence strongly implicates N1 of G8 in active site chemistry. By direct analogy to RNase A, these findings could be consistent with the hypothesis that deprotonated G8 residue functions as a general base in the hairpin ribozyme. Other mechanistic possibilities for N1 of G8 such as indirect general base catalysis mediated by a water molecule or transition state stabilization could also be consistent with our findings.

Calculation of pKas in RNA: on the structural origins and functional roles of protonated nucleotides

pK(a) calculations based on the Poisson-Boltzmann equation have been widely used to study proteins and, more recently, DNA. However, much less attention has been paid to the calculation of pK(a) shifts in RNA. There is accumulating evidence that protonated nucleotides can stabilize RNA structure and participate in enzyme catalysis within ribozymes. Here, we calculate the pK(a) shifts of nucleotides in RNA structures using numerical solutions to the Poisson-Boltzmann equation. We find that significant shifts are predicted for several nucleotides in two catalytic RNAs, the hairpin ribozyme and the hepatitis delta virus ribozyme, and that the shifts are likely to be related to their functions. We explore how different structural environments shift the pK(a)s of nucleotides from their solution values. RNA structures appear to use two basic strategies to shift pK(a)s: (a) the formation of compact structural motifs with structurally-conserved, electrostatic interactions; and (b) the arrangement of the phosphodiester backbone to focus negative electrostatic potential in specific regions.

Non-identical electronic characters of the internucleotidic phosphates in RNA modulate the chemical reactivity of the phosphodiester bonds

We here show that the electronic properties and the chemical reactivities of the internucleotidic phosphates in the heptameric ssRNAs are dissimilar in a sequence-specific manner because of their non-identical microenvironments, in contrast with the corresponding isosequential ssDNAs. This has been evidenced by monitoring the delta H8(G) shifts upon pH-dependent ionization (pK(a1)) of the central 9-guaninyl (G) to the 9-guanylate ion (G-), and its electrostatic effect on each of the internucleotidic phosphate anions, as measured from the resultant delta 31P shifts (pKa2) in the isosequential heptameric ssRNAs vis-à-vis ssDNAs: [d/r(5'-Cp1Ap2Q1p3Gp4Q2p5Ap6C-3'): Q1 = Q2 = A (5a/5b) or C (8a/8b), Q(1) = A, Q(2) = C (6a/6b), Q1 = C, Q2 = A (7a/7b)]. These oligos with single ionizable G in the centre are chosen because of the fact that the pseudoaromatic character of G can be easily modulated in a pH-dependent manner by its transformation to G- (the 2'-OH to 2-O- ionization effect is not detectable below pH 11.6 as evident from the N(1-Me)-G analog), thereby modulating/titrating the nature of the electrostatic interactions of G to G- with the phosphates, which therefore constitute simple models to interrogate how the variable pseudoaromatic characters of nucleobases under different sequence context (J. Am. Chem. Soc., 2004, 126, 8674-8681) can actually influence the reactivity of the internucleotide phosphates as a result of modulation of sequence context-specific electrostatic interactions. In order to better understand the impact of the electrostatic effect of the G to G- on the tunability of the electronic character of internucleotidic phosphates in the heptameric ssRNAs 5b, 6b, 7b and 8b, we have also performed their alkaline hydrolysis at pH 12.5 at 20 degrees C, and have identified the preferences of the cleavage sites at various phosphates, which are p2, p3 and p4 (Fig.3). The results of these alkaline hydrolysis studies have been compared with the hydrolysis of analogous N(1-Me)-G heptameric ssRNA sequences 5c, 7c and 8c under identical conditions in order to establish the role of the electrostatic effect of the 9-guanylate ion (and the 2'-OH to 2-O- ionization) on the internucleotidic phosphate. It turned out that the relative alkaline hydrolysis rate at those particular phosphates (p2, p3 and p4) in the N(1-Me)-G heptamers was reduced from 16-78% compared to those in the native counterparts [Fig. 4, and ESI 2 (Fig. S11)]. Thus, these physico-chemical studies have shown that those p2, p3 and p4 phosphates in the native heptameric RNAs, which show pKa2 as well as more deshielding (owing to weaker 31P screening) in the alkaline pH compared to those at the neutral pH, are more prone to the alkaline hydrolysis because of their relatively enhanced electrophilic character resulting from weaker 31P screening. This screening effect originates as a result of the systematic charge repulsion effect between the electron cloud in the outermost orbitals of phosphorus and the central guanylate ion, leading to delocalization of the phosphorus p(pi) charge into its dpi orbitals. It is thus likely that, just as in the non-enzymatic hydrolysis, the enzymatic hydrolysis of a specific phosphate in RNA by general base-catalysis in RNA-cleaving proteins (RNase A, RNA phosphodiesterase or nuclease) can potentially be electrostatically influenced by tuning the transient charge on the nucleobase in the steric proximity or as a result of specific sequence context owing to nearest-neighbor interactions.

Ribozyme catalysis: not different, just worse

Evolution has resoundingly favored protein enzymes over RNA-based catalysts, yet ribozymes occupy important niches in modern cell biology that include the starring role in catalysis of protein synthesis on the ribosome. Recent results from structural and biochemical studies show that natural ribozymes use an impressive range of catalytic mechanisms, beyond metalloenzyme chemistry and analogous to more chemically diverse protein enzymes. These findings make it increasingly possible to compare details of RNA- and protein-based catalysis.

Linkage between proton binding and folding in RNA: a thermodynamic framework and its experimental application for investigating pKa shifting

Perturbation of pKa values can change the favored protonation states of the nucleobases at biological pH and thereby modulate the function of RNA and DNA molecules. In an effort to understand the driving forces for pKa shifting specific to nucleic acids, we developed a thermodynamic framework that relates proton binding to the nucleobases and the helix-coil transition. Key features that emerge from the treatment are a comprehensive description of all the actions of proton binding on RNA folding: acid and alkaline denaturation of the helix and pKa shifting in the folded state. Practical experimental approaches for measuring pKas from thermal denaturation experiments are developed. Microscopic pka values (where ka is the acid dissociation constant) for the unfolded state were determined directly by experiments on unstructured oligonucleotides, which led to a macroscopic pKa for the ensemble of unfolded states shifted toward neutrality. The formalism was then applied to pH-dependent UV melting data for model DNA oligonucleotides. Folded-state pka) values were in good agreement with the outcome of pH titrations, and the acid and alkaline denaturation regions were well described. The formalism developed here is similar to that of Draper and coworkers for Mg2+ binding to RNA, except that the unfolded state is described explicitly owing to the presence of specific proton-binding sites on the bases. A principal conclusion is that it should be possible to attain large pKa shifts by designing RNA molecules that fold cooperatively.

Use of both direct and indirect 13C tags for probing nitrogen interactions in hairpin ribozyme models by 15N NMR

We have used the synthesis and 15N NMR study of separate loop A and loop B domains of the hairpin ribozyme to demonstrate that multiple 15N atoms can be incorporated into an RNA strand and be unambiguously distinguished through a combination of direct and indirect tagging by 13C atoms. Absence of 15N chemical shift changes shows that the G8N1 in loop A does not become deprotonated up to pH 8, and that the G21N7 of loop B does not bind to Mg2+.

Significant pKa perturbation of nucleobases is an intrinsic property of the sequence context in DNA and RNA

The pH titration and NMR studies (pH 6.6-12.5) in the heptameric isosequential ssDNA and ssRNA molecules, [d/r(5'-CAQ1GQ2AC-3', with variable Q1/Q2)], show that the pKa of the central G residue within the heptameric ssDNAs (DeltapKa = 0.67 +/- 0.03) and ssRNAs (DeltapKa = 0.49 +/- 0.02) is sequence-dependent. This variable pKa of the G clearly shows that its pseudoaromatic character, hence, its chemical reactivity, is strongly modulated and tuned by its sequence context. In contradistinction to the ssDNAs, the electrostatic transmission of the pKa of the G moiety to the neighboring A or C residues in the heptameric ssRNAs (as observed by the response of the aromatic marker protons of As or Cs) is found to be uniquely dependent upon the sequence composition. This demonstrates that the neighboring As or Cs in ssRNAs have variable electrostatic efficiency to interact with the central G/G-, which is owing to the variable pseudoaromatic characters (giving variable chemical reactivities) of the flanking As or Cs compared to those of the isosequential ssDNAs. The sequence-dependent variation of pKa of the central G and the modulation of its pKa transmission through the nearest-neighbors by variable electrostatic interaction is owing to the electronically coupled nature of the constituent nucleobases across the single strand, which demonstrates the unique chemical basis of the sequence context specificity of DNA or RNA in dictating the biological interaction, recognition, and function with any specific ligand.

NMR structure of the 101-nucleotide core encapsidation signal of the Moloney murine leukemia virus

The full length, positive-strand genome of the Moloney Murine Leukemia Virus contains a "core encapsidation signal" that is essential for efficient genome packaging during virus assembly. We have determined the structure of a 101-nucleotide RNA that contains this signal (called mPsi) using a novel isotope-edited NMR approach. The method is robust and should be generally applicable to larger RNAs. mPsi folds into three stem loops, two of which (SL-C and SL-D) co-stack to form an extended helix. The third stem loop (SL-B) is connected to SL-C by a flexible, four-nucleotide linker. The structure contains five mismatched base-pairs, an unusual C.CG base-triple platform, and a novel "A-minor K-turn," in which unpaired adenosine bases A340 and A341 of a GGAA bulge pack in the minor groove of a proximal stem, and a bulged distal uridine (U319) forms a hydrogen bond with the phosphodiester of A341. Phylogenetic analyses indicate that these essential structural elements are conserved among the murine C-type retroviruses.

Adenosylcobinamide Plus Exogenous, Sterically Hindered, Putative Axial Bases: A Reinvestigation into the Cause of Record Levels of Co−C Heterolysis

A reinvestigation of an earlier Ph.D. thesis (Sirovatka, J. M. Ph.D. Thesis, Colorado State University, Fort Collins, CO, 1999) is reported herein. That thesis examined the thermolysis reaction of AdoCbi(+)BF(4)(-) in ethylene glycol solution with exogenous bases, N-methylimidazole (N-Me-Im) and the sterically hindered 1,2-dimethylimidazole, (1,2-Me(2)-Im), 2-methylpyridine (2-Me-py), and 2,6-dimethylpyridine (2,6-Me(2)-py). In the present work, multiple purities of each base have been utilized as a check to see if impurities in the nitrogenous bases are causing the observed homolysis and heterolysis product distributions as others have implied (Trommel, J. S.; Warncke, K.; Marzilli, L. G. J. Am. Chem. Soc. 2001, 123, 3358). The "impurity hypothesis" is disproven by a series of results, including the following: N-Me-Im displays an invariant 52 +/- 10% heterolysis and the 1,2-Me(2)-Im system displays an invariant 83 +/- 7% heterolysis as a function of different base purification methods. Moreover, 2-Me-py and 2,6-Me(2)-py also display an invariant approximately 16 +/- 5% heterolysis as a function of different purification methods. What is responsible for the high levels of Co-C heterolysis in the AdoCbi(+) plus sterically bulky base thermolyses was uncovered via a revisitation of our four, earlier alternative hypotheses for the enhanced Co-C heterolysis (Sirovatka, J. M.; Finke, R. G. Inorg. Chem. 1999, 38, 1697). Our prior number one alternative hypothesis is shown to be correct: the added bases simply deprotonate the ethylene glycol solvent, forming ethylene glycolate anion and base-H(+)() as the key agents behind the previously ill-understood Co-C heterolyses. Also reported are Co(II)Cbi(+) titrations with five bases (1,2-Me(2)-Im, N-Me-Im, pyridine, 2-MePy, and 2,6-Me(2)-py). These experiments confirm Marzilli and co-workers' (op. cit.) results by showing that sterically hindered bases do not bind to Co(II)Cbi(+); therefore, Co(II)Cbi(+) EPR literature showing binding of bulky pyridines is erroneous as is the previously reported binding of bulky pyridine bases to Co(II)Cbi(+) by UV-vis spectroscopy (Sirovatka, J. Ph.D. Thesis, op. cit.). Also reported is our current best synthesis and purification of AdoCbi(+)BF(4)(-), work that builds off our 1987 synthesis of AdoCbi(+)BF(4)(-) (Hay, B. P.; Finke, R. G. J. Am. Chem. Soc. 1987, 109, 8012). Finally, the multiple, compounding errors which have caused problems in this project are listed, notably errors in the protein X-ray crystallography literature, the EXAFS literature, the Co(II)Cbi(+) plus bulky-bases EPR literature, the misleading B(12)-model literature, the erroneous experimental work (Sirovatka, op. cit.) and thus incorrect conclusions in one of our prior papers, as well as the erroneous implications in parts of the Marzilli and co-workers paper (op. cit.). It is hoped that a forthright reporting of these errors will help others avoid similar mistakes in the future when studying complex, bioinorganic systems.

A proton-coupled dynamic conformational switch in the HIV-1 dimerization initiation site kissing complex

In HIV type 1 (HIV-1), the dimerization initiation site (DIS) is the sequence primarily responsible for initiating the noncovalent linkage of two homologous strands of genomic RNA during viral assembly. The DIS loop contains an autocomplementary hexanucleotide sequence and forms a symmetric homodimer through a loop-loop kissing interaction. In a structural rearrangement catalyzed by the HIV-1 nucleocapsid protein (NCp7) and suggested to be associated with maturation of the budded viral particle, the DIS converts from a metastable kissing dimer to an extended duplex. Here, we demonstrate that the DIS kissing dimer displays localized conformational dynamics that result from the specific protonation of the N1 base nitrogen of the DIS loop residue A272 at near-physiological pH. The rate of NCp7-catalyzed maturation of the DIS kissing dimer is also shown to directly correlate with the observed proton-coupled conformational dynamics, where NCp7 is found to convert the dynamic A272 protonated state with a faster rate. Taken together, these results reveal a previously undescribed role for base protonation in modulating local RNA structure and demonstrate a mechanism for promoting the chaperone-mediated structural rearrangement of a kinetically trapped RNA conformational state.

Naturally-occurring modification restricts the anticodon domain conformational space of tRNA(Phe)

Post-transcriptional modifications contribute chemistry and structure to RNAs. Modifications of tRNA at nucleoside 37, 3'-adjacent to the anticodon, are particularly interesting because they facilitate codon recognition and negate translational frame-shifting. To assess if the functional contribution of a position 37-modified nucleoside defines a specific structure or restricts conformational flexibility, structures of the yeast tRNA(Phe) anticodon stem and loop (ASL(Phe)) with naturally occurring modified nucleosides differing only at position 37, ASL(Phe)-(Cm(32),Gm(34),m(5)C(40)), and ASL(Phe)-(Cm(32),Gm(34),m(1)G(37),m(5)C(40)), were determined by NMR spectroscopy and restrained molecular dynamics. The ASL structures had similarly resolved stems (RMSD approximately 0.6A) of five canonical base-pairs in standard A-form RNA. The "NOE walk" was evident on the 5' and 3' sides of the stems of both RNAs, and extended to the adjacent loop nucleosides. The NOESY cross-peaks involving U(33) H2' and characteristic of tRNA's anticodon domain U-turn were present but weak, whereas those involving the U(33) H1' proton were absent from the spectra of both ASLs. However, ASL(Phe)-(Cm(32),Gm(34),m(1)G(37),m(5)C(40)) exhibited the downfield shifted 31P resonance of U(33)pGm(34) indicative of U-turns; ASL(Phe)-(Cm(32),Gm(34),m(5)C(40)) did not. An unusual "backwards" NOE between Gm(34) and A(35) (Gm(34)/H8 to A(35)/H1') was observed in both molecules. The RNAs exhibited a protonated A(+)(38) resulting in the final structures having C(32).A(+)(38) intra-loop base-pairs, with that of ASL(Phe)-(Cm(32),Gm(34),m(1)G(37),m(5)C(40)) being especially well defined. A single family of low-energy structures of ASL(Phe)-(Cm(32),Gm(34), m(1)G(37),m(5)C(40)) (loop RMSD 0.98A) exhibited a significantly restricted conformational space for the anticodon loop in comparison to that of ASL(Phe)-(Cm(32),Gm(34),m(5)C(40)) (loop RMSD 2.58A). In addition, the ASL(Phe)-(Cm(32),Gm(34),m(1)G(37),m(5)C(40)) average structure had a greater degree of similarity to that of the yeast tRNA(Phe) crystal structure. A comparison of the resulting structures indicates that modification of position 37 affects the accuracy of decoding and the maintenance of the mRNA reading frame by restricting anticodon loop conformational space.

Phylogenetic conservation of RNA secondary and tertiary structure in the trpEDCFBA operon leader transcript in Bacillus

Expression of the trpEDCFBA operon of Bacillus subtilis is regulated by transcription attenuation and translation control mechanisms. We recently determined that the B. subtilis trp leader readthrough transcript can adopt a Mg(2+)-dependent tertiary structure that appears to interfere with TRAP-mediated translation control of trpE. In the present study, sequence comparisons to trp leaders from three other Bacillus sp. were made, suggesting that RNA secondary and tertiary structures are phylogenetically conserved. To test this hypothesis, experiments were carried out with the trp leader transcript from Bacillus stearothermophilus. Structure mapping experiments confirmed the predicted secondary structure. Native gel experiments identified a faster mobility species in the presence of Mg(2+), suggesting that a Mg(2+)-dependent tertiary structure forms. Mg(2+)-dependent protection of residues within the first five triplet repeats of the TRAP binding target and a pyrimidine-rich internal loop were observed, consistent with tertiary structure formation between these regions. Structure mapping in the presence of a competitor DNA oligonucleotide allowed the interacting partners to be identified as a single-stranded portion of the purine-rich TRAP binding target and the large downstream pyrimidine-rich internal loop. Thermal denaturation experiments revealed a Mg(2+)- and pH-dependent unfolding transition that was absent for a transcript missing the first five triplet repeats. The stability of several mutant transcripts allowed a large portion of the base-pairing register for the tertiary interaction to be determined. These data indicate that RNA secondary and tertiary structures involved in TRAP-mediated translation control are conserved in at least four Bacillus species.

Non-Watson-Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E

Explicit solvent and counterion molecular dynamics simulations have been carried out for a total of >80 ns on the bacterial and spinach chloroplast 5S rRNA Loop E motifs. The Loop E sequences form unique duplex architectures composed of seven consecutive non-Watson-Crick basepairs. The starting structure of spinach chloroplast Loop E was modeled using isostericity principles, and the simulations refined the geometries of the three non-Watson-Crick basepairs that differ from the consensus bacterial sequence. The deep groove of Loop E motifs provides unique sites for cation binding. Binding of Mg(2+) rigidifies Loop E and stabilizes its major groove at an intermediate width. In the absence of Mg(2+), the Loop E motifs show an unprecedented degree of inner-shell binding of monovalent cations that, in contrast to Mg(2+), penetrate into the most negative regions inside the deep groove. The spinach chloroplast Loop E shows a marked tendency to compress its deep groove compared with the bacterial consensus. Structures with a narrow deep groove essentially collapse around a string of Na(+) cations with long coordination times. The Loop E non-Watson-Crick basepairing is complemented by highly specific hydration sites ranging from water bridges to hydration pockets hosting 2 to 3 long-residing waters. The ordered hydration is intimately connected with RNA local conformational variations.

Tandem electrostatic effect from the first to the third aglycon in the trimeric RNA owing to the nearest-neighbor interaction

We here show an electrostatic polar-pi interaction from the first to the third aglycon, via the second aglycon, in the ground state in two single stranded trimeric RNAs, 5'-GpA(1)pA(2)-3' (3) and 5'-GpApC-3' (4), as a result of intramolecular nearest neighbor offset-stacking. The experimental evidence in support of this conclusion has been obtained by comparing the pK(a)s of each aglycone in the two trimers with those of guanosine 3'-ethyl phosphate, GpEt (1) and 5'-GpA-3' (2): Thus, the pK(a) of N(1)-H of guanin-9-yl of 5'-GpA(1)pA(2)-3' (3) could be measured by pH titration (pH 7.3-11.6) of its own deltaH8G (pK(a) 9.75 +/- 0.02) as well as from deltaH8A(1) (pK(a) 9.72 +/- 0.02) and deltaH2A(1) (pK(a) 9.83 +/- 0.04) of the neighboring pA(1)p moiety and the deltaH8A(2) (pK(a) 9.83 +/- 0.02) of the terminal pA(2) moiety. Similarly, the pH titration of GpApC (4) shows the pK(a) of N(1)-H of guanin-9-yl from its own deltaH8G (pK(a) 9.88 +/- 0.03) as well as from deltaH8A (pK(a) 9.87 +/- 0.01) of the neighboring pAp moiety, and deltaH5/H6C (pK(a) 9.88 +/- 0.01 and 9.90 +/- 0.01 respectively) of the 3'-terminal cytosin-1-yl. This intramolecular nearest neighbor electrostatic interaction in the single-stranded RNA modulates the pseudoaromaticity of the nearest neighbors by almost total transmission of because they constitute an extended array of offset-stacked coupled aromatic heterocycles within a polyanionic sugar-phosphate backbone at the ground state. The enhanced basicity of Gp residue by ca. 0.6 pK(a) unit in the trimers compared to that of the dimer is a result of the change in the electrostatic microenvironment owing to the nearest neighbors in the former (the nucleobases as well as the phosphates). Thus, the from the 5'-guanylate ion to the 3'-end aglycon via the central adenin-9-yl is 55 to 56 kJ mol(-)(1) in each step through a distance spanning approximately 6.8 A in an unfolded state. As a result, the pK(a) of guanin-9-yl moiety has become 9.25 +/- 0.02 in GpEt (1), 9.17 +/- 0.02 in GpA (2), 9.75 +/- 0.02 in GpApA (3), and 9.88 +/- 0.03 in GpApC (4). This means that guanin-9-yl moiety of trimers 3 and 4 is more basic than in the monomer or the dimer. The net outcome of this electrostatic cross-talk between the two neighboring heterocycles is creation of new hybrid aglycones in an oligo or polynucleotide, whose physicochemical property and the pseudoaromatic character are completely dependent both upon the nearest neighbors, and whether they are stacked or unstacked. Thus, this tunable physicochemical property of an aglycon (an array of the extended genetic code) may have considerable implication in our understanding of the specific ligand binding ability of an aptamer, the pK(a) and the hydrogen bonding ability of nucleic acids in a microenvironment, or in the triplet usage by the anticodon-codon interaction in the protein biosynthesis in the ribosome.

Binding to an RNA aptamer changes the charge distribution and conformation of malachite green

RNA plays a central role in many biological processes and is therefore an important target for drug development. In recent years an increasing wealth of structural and functional information about RNA-ligand complexes has been obtained using in vitro selected RNAs (aptamers). However, all those studies focused on structure and changes of the nucleic acid and mostly considered the ligand as a rigid target. To develop a detailed picture of ligand structure and dynamics in RNA-small molecule complexes, the malachite green binding aptamer was studied. Isotopically labeled ligand in complex with RNA was analyzed by NMR spectroscopy in solution. The surprisingly asymmetric changes in the (13)C chemical shift of the ligand methyl groups indicate that the dye undergoes changes in its conformation and charge distribution upon binding. The role of the RNA electrostatic field in this interaction was explored using ab initio calculations of the ligand structure and charge distribution. The results indicate that the uneven charge distribution in the RNA binding pocket provides a major contribution to the driving force of the ligand structural changes. The observation that not only the RNA adapts to the ligand, in what is called adaptive binding, but that the ligand itself also undergoes conformational changes ("induced fit") is crucial for the rational design of RNA ligands and for understanding the properties of RNA-ligand complexes.