Solution structure and thermodynamics of a divalent metal ion binding site in an RNA pseudoknot

Applicability of the nearest-neighbour model for pseudoknot RNAs

The validity of the nearest-neighbour (NN) model was verified in an RNA pseudoknot (PK) structure. The thermodynamic parameters of the second hairpin stem (S2) region, which separates the PK from a hairpin structure, were monitored using CD and UV melting. Different PKs with identical NN base pairs in the S2 region exhibited similar thermodynamic parameters, highlighting the validity of the NN model in this RNA tertiary structure motif.

Nearest neighbor rules for RNA helix folding thermodynamics: improved end effects

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5′GGUC3′ paired to 3′CUGG5′ separately from other 5′GU3′/3′UG5′ stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37°C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model.

Thermodynamics of unfolding mechanisms of mouse mammary tumor virus pseudoknot from a coarse-grained loop-entropy model

Pseudoknotted RNA molecules play important biological roles that depend on their folded structure. To understand the underlying principles that determine their thermodynamics and folding/unfolding mechanisms, we carried out a study on a variant of the mouse mammary tumor virus pseudoknotted RNA (VPK), a widely studied model system for RNA pseudoknots. Our method is based on a coarse-grained discrete-state model and the algorithm of PK3D (pseudoknot structure predictor in three-dimensional space), with RNA loops explicitly constructed and their conformational entropic effects incorporated. Our loop entropy calculations are validated by accurately capturing previously measured melting temperatures of RNA hairpins with varying loop lengths. For each of the hairpins that constitutes the VPK, we identified alternative conformations that are more stable than the hairpin structures at low temperatures and predicted their populations at different temperatures. Our predictions were validated by thermodynamic experiments on these hairpins. We further computed the heat capacity profiles of VPK, which are in excellent agreement with available experimental data. Notably, our model provides detailed information on the unfolding mechanisms of pseudoknotted RNA. Analysis of the distribution of base-pairing probability of VPK reveals a cooperative unfolding mechanism instead of a simple sequential unfolding of first one stem and then the other. Specifically, we find a simultaneous "loosening" of both stems as the temperature is raised, whereby both stems become partially melted and co-exist during the unfolding process.

Monovalent ions modulate the flux through multiple folding pathways of an RNA pseudoknot

The assembly mechanism of RNA, vital to describing its functions, depends on both the sequence and the metal ion concentration. How the latter influences the folding trajectories remains an important unsolved problem. Here, we examine the folding pathways of an RNA pseudoknot (PK) with key functional roles in transcription and translation, using a combination of experiments and simulations. We demonstrate that the PK, consisting of two hairpins with differing stabilities, folds by parallel pathways. Surprisingly, the flux between them is modulated by monovalent salt concentration. Our work shows that the order of assembly of PKs is determined by the relative stability of the hairpins, implying that the folding landscape can be controlled by sequence and ion concentration.

The functions of RNA pseudoknots (PKs), which are minimal tertiary structural motifs and an integral part of several ribozymes and ribonucleoprotein complexes, are determined by their structure, stability, and dynamics. Therefore, it is important to elucidate the general principles governing their thermodynamics/folding mechanisms. Here, we combine laser temperature-jump experiments and coarse-grained simulations to determine the folding/unfolding pathways of VPK, a variant of the mouse mammary tumor virus (MMTV) PK involved in ribosomal frameshifting. Fluorescent nucleotide analogs (2-aminopurine and pyrrolocytidine) placed at different stem/loop positions in the PK serve as local probes allowing us to monitor the order of assembly of VPK that has two constituent hairpins with different intrinsic stabilities. We show that at 50 mM KCl, the dominant folding pathway populates only the more stable hairpin intermediate; as the salt concentration is increased, a parallel folding pathway emerges involving the less stable hairpin as an alternate intermediate. Notably, the flux between the pathways is modulated by the ionic strength. Our findings support the principle that the order of PK structure formation is determined by the relative stabilities of the hairpins, which can be altered by sequence variations or salt concentrations. The experimental results of salt effects on the partitioning between the two folding pathways are in remarkable agreement with simulations that were performed with no adjustable parameters. Our study not only unambiguously demonstrates that VPK folds by parallel pathways but also showcases the power of combining experiments and simulations for a more enriched description of RNA self-assembly.

Characterization of G-Quadruplex/Hairpin Transitions of RNAs by 19 F NMR Spectroscopy

2'-O-[(4-Trifluoromethyl-triazol-1-yl)methyl] reporter groups have been incorporated into guanosine-rich RNA models (including a known bistable Qd/Hp RNA and two G-rich regions of mRNA of human prion protein, PrP) and applied for the ¹⁹ F NMR spectroscopic characterization of plausible G-quadruplex/hairpin (Qd/Hp) transitions in these RNA structures. For the synthesis of the CF₃ -labeled RNAs, phosphoramidite building blocks of 2'-O-[(4-CF₃ -triazol-1-yl)methyl] nucleosides (cytidine, adenosine, and guanosine) were prepared and used as an integral part of the standard solid-phase RNA synthesis. The obtained ¹⁹ F NMR spectra supported the usual characterization data (obtained by UV- and CD-melting profiles and by ¹ H NMR spectra of the imino regions) and additionally gave more detailed information on the Qd/Hp transitions. The molar fractions of the secondary structural species (Qd, Hp) upon thermal denaturation and under varying ionic conditions could be determined from the intensities and shifts of the ¹⁹ F NMR signals. For a well-behaved Qd/Hp transition, thermodynamic parameters could be extracted.

Ribosomal frameshifting and transcriptional slippage: From genetic steganography and cryptography to adventitious use

Genetic decoding is not ‘frozen’ as was earlier thought, but dynamic. One facet of this is frameshifting that often results in synthesis of a C-terminal region encoded by a new frame. Ribosomal frameshifting is utilized for the synthesis of additional products, for regulatory purposes and for translational ‘correction’ of problem or ‘savior’ indels. Utilization for synthesis of additional products occurs prominently in the decoding of mobile chromosomal element and viral genomes. One class of regulatory frameshifting of stable chromosomal genes governs cellular polyamine levels from yeasts to humans. In many cases of productively utilized frameshifting, the proportion of ribosomes that frameshift at a shift-prone site is enhanced by specific nascent peptide or mRNA context features. Such mRNA signals, which can be 5′ or 3′ of the shift site or both, can act by pairing with ribosomal RNA or as stem loops or pseudoknots even with one component being 4 kb 3′ from the shift site. Transcriptional realignment at slippage-prone sequences also generates productively utilized products encoded trans-frame with respect to the genomic sequence. This too can be enhanced by nucleic acid structure. Together with dynamic codon redefinition, frameshifting is one of the forms of recoding that enriches gene expression.

Mimicking Ribosomal Unfolding of RNA Pseudoknot in a Protein Channel

Pseudoknots are a fundamental RNA tertiary structure with important roles in regulation of mRNA translation. Molecular force spectroscopic approaches such as optical tweezers can track the pseudoknot's unfolding intermediate states by pulling the RNA chain from both ends, but the kinetic unfolding pathway induced by this method may be different from that in vivo, which occurs during translation and proceeds from the 5' to 3' end. Here we developed a ribosome-mimicking, nanopore pulling assay for dissecting the vectorial unfolding mechanism of pseudoknots. The pseudoknot unfolding pathway in the nanopore, either from the 5' to 3' end or in the reverse direction, can be controlled by a DNA leader that is attached to the pseudoknot at the 5' or 3' ends. The different nanopore conductance between DNA and RNA translocation serves as a marker for the position and structure of the unfolding RNA in the pore. With this design, we provided evidence that the pseudoknot unfolding is a two-step, multistate, metal ion-regulated process depending on the pulling direction. Most notably, unfolding in both directions is rate-limited by the unzipping of the first helix domain (first step), which is Helix-1 in the 5' → 3' direction and Helix-2 in the 3' → 5' direction, suggesting that the initial unfolding step in either pulling direction needs to overcome an energy barrier contributed by the noncanonical triplex base-pairs and coaxial stacking interactions for the tertiary structure stabilization. These findings provide new insights into RNA vectorial unfolding mechanisms, which play an important role in biological functions including frameshifting.

Structural determinants for alternative splicing regulation of the MAPT pre-mRNA

Alternative splicing at the MAPT gene exon 10 yields similar levels of the 3R and 4R tau protein isoforms. (1) The presence of mutations, particularly in exon 10 and intron 10-11, changes the quantity of tau isoforms. Domination each of the isoform yields tau protein aggregation and frontotemporal dementia and Parkinsonism linked to chromosome 17 (FTDP-17). Here, we report for the first time the secondary structure of the 194/195 nucleotide region for the wild type (WT) and 10 mutants of the MAPT gene pre-mRNA determined using both chemical and microarray mapping. Thermodynamic analyses indicate that single nucleotide mutations in the splicing regulatory element (SRE) that form a hairpin affect its stability by up to 4 and 7 kcal/mol. Moreover, binding the regulatory hairpin of small molecule ligands (neomycin, kanamycin, tobramycin and mitoxantrone) enhance its stability depending on the nature of the ligands and the RNA mutations. Experiments using the cos-7 cell line indicate that the presence of ligands and modified antisense oligonucleotides affect the quantity of 3R and 4R isoforms. This finding correlates with the thermodynamic stability of the regulatory hairpin. An alternative splicing regulation mechanism for exon 10 is postulated based on our experimental data and on published data.

On the interaction between [Ru(NH3)6]3+ and the G-quadruplex forming thrombin binding aptamer sequence

The interaction between the thrombin binding aptamer (TBA), a G-quadruplex forming DNA sequence, and the electroactive hexaammineruthenium(III) cation has been studied by electrochemical methods and circular dichroism spectroscopy. When TBA is immobilised on a gold surface in a typical aptasensor configuration, the [Ru(NH3)6](3+) cation can be bound to the electrode surface through its interaction with the TBA sequence. This interaction is strong enough to enable the ruthenium complex to remain at the surface when the electrode is immersed in an electrolyte free of [Ru(NH3)6](3+), meaning that the complex does not diffuse back into the solution. A stoichiometry of 2 [Ru(NH3)6](3+) per TBA strand has been determined, indicating that the interaction differs from the conventional, non-specific electrostatic charge compensation, for which a 5 to 1 ratio would be expected between the triply charged cation and the 15 bases sequence. It is shown that this interaction takes place not only at the surface, but also when both TBA and hexaammineruthenium(III) are dissolved in solution. Under such conditions, a similar stoichiometry of 2 [Ru(NH3)6](3+) per TBA strand has been evidenced by two independent methods, namely circular dichroism spectroscopy and differential pulse voltammetry.

The structural stabilization of the κ three-way junction by Mg(II) represents the first step in the folding of a group II intron

Folding of group II introns is characterized by a first slow compaction of domain 1 (D1) followed by the rapid docking of other domains to this scaffold. D1 compaction initiates in a small subregion encompassing the κ and ζ elements. These two tertiary elements are also the major interaction sites with domain 5 to form the catalytic core. Here, we provide the first characterization of the structure adopted at an early folding step and show that the folding control element can be narrowed down to the three-way junction with the κ motif. In our nuclear magnetic resonance studies of this substructure derived from the yeast mitochondrial group II intron Sc.ai5γ, we show that a high affinity Mg(II) ion stabilizes the κ element and enables coaxial stacking between helices d′ and d′′, favoring a rigid duplex across the three-way junction. The κ-element folds into a stable GAAA-tetraloop motif and engages in A-minor interactions with helix d′. The addition of cobalt(III)hexammine reveals three distinct binding sites. The Mg(II)-promoted structural rearrangement and rigidification of the D1 core can be identified as the first micro-step of D1 folding.

An RNA hairpin to G-quadruplex conformational transition

RNA molecules can fold into noncanonical structures such as the four-stranded structures known as G-quadruplexes. G-quadruplexes in the transcriptome have recently emerged as relevant regulatory elements of gene expression. Conformational transitions in RNA molecules offer an important way to regulate their biological functions. Here we report on the competition between a canonical hairpin structure and a G-quadruplex structure within an RNA molecule. We show that the conformational preference strongly depends on the relative amounts of mono- and divalent metal ions present in solution. In our system, the G-quadruplex, whose formation is not predicted by available predictive RNA folding programs, is the major conformer at physiologically relevant K(+) and Mg(2+) concentrations. Furthermore, we show that a synthetic small molecule can displace the structural dynamic equilibrium in favor of the hairpin conformer. This work highlights a new and important level of complexity in RNA folding that could be relevant to the biological functions and targeting of RNAs comprising G-quadruplex motifs.

The 3' splice site of influenza A segment 7 mRNA can exist in two conformations: a pseudoknot and a hairpin

The 3′ splice site of influenza A segment 7 is used to produce mRNA for the M2 ion-channel protein, which is critical to the formation of viable influenza virions. Native gel analysis, enzymatic/chemical structure probing, and oligonucleotide binding studies of a 63 nt fragment, containing the 3′ splice site, key residues of an SF2/ASF splicing factor binding site, and a polypyrimidine tract, provide evidence for an equilibrium between pseudoknot and hairpin structures. This equilibrium is sensitive to multivalent cations, and can be forced towards the pseudoknot by addition of 5 mM cobalt hexammine. In the two conformations, the splice site and other functional elements exist in very different structural environments. In particular, the splice site is sequestered in the middle of a double helix in the pseudoknot conformation, while in the hairpin it resides in a two-by-two nucleotide internal loop. The results suggest that segment 7 mRNA splicing can be controlled by a conformational switch that exposes or hides the splice site.

Characterization of metal ion-nucleic acid interactions in solution

Metal ions are inextricably involved with nucleic acids due to their polyanionic nature. In order to understand the structure and function of RNAs and DNAs, one needs to have detailed pictures on the structural, thermodynamic, and kinetic properties of metal ion interactions with these biomacromolecules. In this review we first compile the physicochemical properties of metal ions found and used in combination with nucleic acids in solution. The main part then describes the various methods developed over the past decades to investigate metal ion binding by nucleic acids in solution. This includes for example hydrolytic and radical cleavage experiments, mutational approaches, as well as kinetic isotope effects. In addition, spectroscopic techniques like EPR, lanthanide(III) luminescence, IR and Raman as well as various NMR methods are summarized. Aside from gaining knowledge about the thermodynamic properties on the metal ion-nucleic acid interactions, especially NMR can be used to extract information on the kinetics of ligand exchange rates of the metal ions applied. The final section deals with the influence of anions, buffers, and the solvent permittivity on the binding equilibria between metal ions and nucleic acids. Little is known on some of these aspects, but it is clear that these three factors have a large influence on the interaction between metal ions and nucleic acids.

Fast folding of RNA pseudoknots initiated by laser temperature-jump

RNA pseudoknots are examples of minimal structural motifs in RNA with tertiary interactions that stabilize the structures of many ribozymes. They also play an essential role in a variety of biological functions that are modulated by their structure, stability, and dynamics. Therefore, understanding the global principles that determine the thermodynamics and folding pathways of RNA pseudoknots is an important problem in biology, both for elucidating the folding mechanisms of larger ribozymes as well as addressing issues of possible kinetic control of the biological functions of pseudoknots. We report on the folding/unfolding kinetics of a hairpin-type pseudoknot obtained with microsecond time-resolution in response to a laser temperature-jump perturbation. The kinetics are monitored using UV absorbance as well as fluorescence of extrinsically attached labels as spectroscopic probes of the transiently populated RNA conformations. We measure folding times of 1-6 ms at 37 °C, which are at least 100-fold faster than previous observations of very slow folding pseudoknots that were trapped in misfolded conformations. The measured relaxation times are remarkably similar to predictions of a computational study by Thirumalai and co-workers (Cho, S. S.; Pincus, D.L.; Thirumalai, D. Proc. Natl. Acad. Sci. U. S. A. 2009, 106, 17349-17354). Thus, these studies provide the first observation of a fast-folding pseudoknot and present a benchmark against which computational models can be refined.

Solution structure of an alternate conformation of helix27 from Escherichia coli16S rRNA

Helix (H)27 of 16S ribosomal (r)RNA from Escherichia coli was dubbed the "switch helix" when mutagenesis suggested that two alternative base pair registers may have distinct functional roles in the bacterial ribosome. Although more recent genetic analyses suggest that H27 conformational switching is not required for translation, previous solution studies demonstrated that the isolated E. coli H27 can dynamically convert between the 885 and 888 conformations. Here, we have solved the nuclear magnetic resonance solution structure of a locked 888 conformation. NOE and residual dipolar coupling restraints reveal an architecture that markedly differs from that of the 885 conformation found in crystal structures of the bacterial ribosome. In place of the loop E motif that characterizes the 885 conformer and that the 888 conformer cannot adopt, we find evidence for an asymmetrical A-rich internal loop stabilized by stacking interactions among the unpaired A's. Comparison of the isolated H27 888 solution structure with the 885 crystal structure within the context of the ribosome suggests a difference in overall length of H27 that presents one plausible reason for the absence of H27 conformational switching within the sterically confining ribosome.

Consecutive terminal GU pairs stabilize RNA helices

Consecutive GU pairs at the ends of RNA helices provide significant thermodynamic stability between -1.0 and -3.8 kcal/mol at 37 °C, which is equivalent to approximately 2 orders of magnitude in the value of a binding constant. The thermodynamic stabilities of GU pairs depend on the sequence, stacking orientation, and position in the helix. In contrast to GU pairs in the middle of a helix that may be destabilizing, all consecutive terminal GU pairs contribute favorable thermodynamic stability. This work presents measured thermodynamic stabilities for 30 duplexes containing two, three, or four consecutive GU pairs at the ends of RNA helices and a model to predict the thermodynamic stabilities of terminal GU pairs. Imino proton NMR spectra show that the terminal GU nucleotides form hydrogen-bonded pairs. Different orientations of terminal GU pairs can have different conformations with equivalent thermodynamic stabilities. These new data and prediction model will help improve RNA secondary structure prediction, identification of miRNA target sequences with GU pairs, and efforts to understand the fundamental physical forces directing RNA structure and energetics.

The identity of the nucleophile substitution may influence metal interactions with the cleavage site of the minimal hammerhead ribozyme

Potential metal interactions with the cleavage site of a minimal hammerhead ribozyme (mHHRz) were probed using (31)P NMR-detected Cd(2+) titration studies of HHRz constructs containing a phosphorothioate (PS) modification at the cleavage site. The mHHRz nucleophile position was replaced by either a 2'-F or a 2'-NH(2) in order to block cleavage activity during the study. The 2'-F/PS cleavage site mHHRz construct, in which the 2'-F should closely imitate the atom size and electronegativity of a 2'-OH, demonstrates low levels of metal ion association (<1 ppm (31)P chemical shift changes). This observation indicates that having an atom size and electrostatic properties that are similar to the 2'-OH are not the governing factors in allowing metal interactions with the scissile phosphate of the mHHRz. With a 2'-NH(2) substitution, a large upfield change in (31)P NMR chemical shift of the phosphorothioate peak (Delta approximately 3 ppm with 6 equiv of added Cd(2+)) indicates observable Cd(2+) interactions with the substituted site. Since a 2'-NH(2), but not a 2'-F, can serve as a metal ligand, these data suggest that a metal ion interaction with the HHRz cleavage site may include both the scissile phosphate and the 2' nucleophile. Control samples in which the 2'-NH(2)/PS unit is placed either next to the mHHRz cleavage site (at U16.1), in a duplex, or in a (am)U(PS)U dinucleotide show much weaker interactions with Cd(2+). Results with these control samples indicate that simply the presence of a 2'-NH(2)/PS unit does not create a strong metal binding site, reinforcing the possibility that the 2'-NH(2)-moderated Cd-PS interaction is specific to the mHHRz cleavage site. Upfield chemical shifts of both (31)P and H-2' (1)H resonances in (am)U(PS)U are observed with addition of Cd(2+), consistent with the predicted metal coordination to both 2'-NH(2) and phosphorothioate ligands. These data suggest that metal ion association with the HHRz cleavage site may include an interaction with the 2'-OH nucleophile.

Assembly mechanisms of RNA pseudoknots are determined by the stabilities of constituent secondary structures

Understanding how RNA molecules navigate their rugged folding landscapes holds the key to describing their roles in a variety of cellular functions. To dissect RNA folding at the molecular level, we performed simulations of three pseudoknots (MMTV and SRV-1 from viral genomes and the hTR pseudoknot from human telomerase) using coarse-grained models. The melting temperatures from the specific heat profiles are in good agreement with the available experimental data for MMTV and hTR. The equilibrium free energy profiles, which predict the structural transitions that occur at each melting temperature, are used to propose that the relative stabilities of the isolated helices control their folding mechanisms. Kinetic simulations, which corroborate the inferences drawn from the free energy profiles, show that MMTV folds by a hierarchical mechanism with parallel paths, i.e., formation of one of the helices nucleates the assembly of the rest of the structure. The SRV-1 pseudoknot, which folds in a highly cooperative manner, assembles in a single step in which the preformed helices coalesce nearly simultaneously to form the tertiary structure. Folding occurs by multiple pathways in the hTR pseudoknot, the isolated structural elements of which have similar stabilities. In one of the paths, tertiary interactions are established before the formation of the secondary structures. Our work shows that there are significant sequence-dependent variations in the folding landscapes of RNA molecules with similar fold. We also establish that assembly mechanisms can be predicted using the stabilities of the isolated secondary structures.

Characterizing riboswitch function: identification of Mg2+ binding site in T box antiterminator RNA

T box bacterial genes utilize a riboswitch mechanism to regulate gene expression at the transcriptional level. Complementary base pairing of the 5'-untranslated mRNA with uncharged cognate tRNA stabilizes formation of an antiterminator element and permits complete transcription. In the absence of tRNA, a mutually exclusive RNA terminator element forms and results in transcription termination. This regulatory mechanism requires divalent metal ions at the antitermination event. The structural effects of Mg(2+) binding to antiterminator model RNA were investigated to ascertain if this requirement is due to the presence of a specific metal ion binding site in the antiterminator. Spectroscopic analysis identified the presence of a hydrated, diffuse Mg(2+) binding site. The results indicate that the mechanistic requirement for divalent metal ions is not due to Mg(2+)-induced pre-formation of a functional antiterminator receptor; rather, Mg(2+) binds in a helical region of high phylogenetic sequence conservation adjacent to the tRNA binding site.

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 Mg2+-ions. Furthermore, a role for Mg2+-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 Mg2+, but becomes partially pre-organized for ligand binding in the presence of Mg2+. 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 Mg2+ is more pronounced.

Distinct reaction pathway promoted by non-divalent-metal cations in a tertiary stabilized hammerhead ribozyme

Divalent ion sensitivity of hammerhead ribozymes is significantly reduced when the RNA structure includes appropriate tertiary stabilization. Therefore, we investigated the activity of the tertiary stabilized "RzB" hammerhead ribozyme in several nondivalent ions. Ribozyme RzB is active in spermidine and Na(+) alone, although the cleavage rates are reduced by more than 1,000-fold relative to the rates observed in Mg(2+) and in transition metal ions. The trivalent cobalt hexammine (CoHex) ion is often used as an exchange-inert analog of hydrated magnesium ion. Trans-cleavage rates exceeded 8 min(-1) in 20 mM CoHex, which promoted cleavage through outersphere interactions. The stimulation of catalysis afforded by the tertiary structural interactions within RzB does not require Mg(2+), unlike other extended hammerhead ribozymes. Site-specific interaction with at least one Mg(2+) ion is suggested by CoHex competition experiments. In the presence of a constant, low concentration of Mg(2+), low concentrations of CoHex decreased the rate by two to three orders of magnitude relative to the rate in Mg(2+) alone. Cleavage rates increased as CoHex concentrations were raised further, but the final fraction cleaved was lower than what was observed in CoHex or Mg(2+) alone. These observations suggest that Mg(2+) and CoHex compete for binding and that they cause misfolded structures when they are together. The results of this study support the existence of an alternate catalytic mechanism used by nondivalent ions (especially CoHex) that is distinct from the one promoted by divalent metal ions, and they imply that divalent metals influence catalysis through a specific nonstructural role.

Probing Na(+)-induced changes in the HIV-1 TAR conformational dynamics using NMR residual dipolar couplings: new insights into the role of counterions and electrostatic interactions in adaptive recognition

Many regulatory RNAs undergo large changes in structure upon recognition of proteins and ligands, but the mechanism by which this occurs remains poorly understood. Using NMR residual dipolar coupling (RDCs), we characterized Na+-induced changes in the structure and dynamics of the bulge-containing HIV-1 transactivation response element (TAR) RNA that mirrors changes induced by small molecules bearing a different number of cationic groups. Increasing the Na+ concentration from 25 to 320 mM led to a continuous reduction in the average inter-helical bend angle (from 46 degrees to 22 degrees ), inter-helical twist angle (from 66 degrees to -18 degrees ), and inter-helix flexibility (as measured by an increase in the internal generalized degree of order from 0.56 to 0.74). Similar conformational changes were observed with Mg2+, indicating that nonspecific electrostatic interactions drive the conformational transition, although results also suggest that Na+ and Mg2+ may associate with TAR in distinct modes. The transition can be rationalized on the basis of a population-weighted average of two ensembles comprising an electrostatically relaxed bent and flexible TAR conformation that is weakly associated with counterions and a globally rigid coaxial conformation that has stronger electrostatic potential and association with counterions. The TAR inter-helical orientations that are stabilized by small molecules fall around the metal-induced conformational pathway, indicating that counterions may help predispose the TAR conformation for target recognition. Our results underscore the intricate sensitivity of RNA conformational dynamics to environmental conditions and demonstrate the ability to detect subtle conformational changes using NMR RDCs.

Resolving fast and slow motions in the internal loop containing stem-loop 1 of HIV-1 that are modulated by Mg2+ binding: role in the kissing-duplex structural transition

Stem loop 1 (SL1) is a highly conserved hairpin in the 5'-leader of the human immunodeficiency virus type I that forms a metastable kissing dimer that is converted during viral maturation into a stable duplex with the aid of the nucleocapsid (NC) protein. SL1 contains a highly conserved internal loop that promotes the kissing-duplex transition by a mechanism that remains poorly understood. Using NMR, we characterized internal motions induced by the internal loop in an SL1 monomer that may promote the kissing-duplex transition. This includes micro-to-millisecond secondary structural transitions that cause partial melting of three base-pairs above the internal loop making them key nucleation sites for exchanging strands and nanosecond rigid-body stem motions that can help bring strands into spatial register. We show that while Mg2+ binds to the internal loop and arrests these internal motions, it preserves and/or activates local mobility at internal loop residues G272 and G273 which are implicated in NC binding. By stabilizing SL1 without compromising the accessibility of G272 and G273 for NC binding, Mg2+ may increase the dependence of the kissing-duplex transition on NC binding thus preventing spontaneous transitions from taking place and ensuring that viral RNA and protein maturation occur in concert.

Structural plasticity and Mg2+ binding properties of RNase P P4 from combined analysis of NMR residual dipolar couplings and motionally decoupled spin relaxation

The P4 helix is an essential element of ribonuclease P (RNase P) that is believed to bind catalytically important metals. Here, we applied a combination of NMR residual dipolar couplings (RDCs) and a recently introduced domain-elongation strategy for measuring "motionally decoupled" relaxation data to characterize the structural dynamics of the P4 helix from Bacillus subtilis RNase P. In the absence of divalent ions, the two P4 helical domains undergo small amplitude (approximately 13 degrees) collective motions about an average interhelical angle of 10 degrees. The highly conserved U7 bulge and helical residue C8, which are proposed to be important for substrate recognition and metal binding, are locally mobile at pico- to nanosecond timescales and together form the pivot point for the collective domain motions. Chemical shift mapping reveals significant association of Mg2+ ions at the P4 major groove near the flexible pivot point at residues (A5, G22, G23) previously identified to bind catalytically important metals. The Mg2+ ions do not, however, significantly alter the structure or dynamics of P4. Analysis of results in the context of available X-ray structures of the RNA component of RNase P and structural models that include the pre-tRNA substrate suggest that the internal motions observed in P4 likely facilitate adaptive changes in conformation that take place during folding and substrate recognition, possibly aided by interactions with Mg2+ ions. Our results add to a growing view supporting the existence of functionally important internal motions in RNA occurring at nanosecond timescales.

The conformational landscape of the ribosomal protein S15 and its influence on the protein interaction with 16S RNA

The interaction between the ribosomal protein S15 and its binding sites in the 16S RNA was examined from two points of view. First, the isolated protein S15 was studied by comparing NMR conformer sets, available in the PDB and recalculated using the CNS-ARIA protocol. Molecular dynamics (MD) trajectories were then recorded starting from a conformer of each set. The recalculation of the S15 NMR structure, as well as the recording of MD trajectories, reveals that several orientations of the N-terminal alpha-helix alpha1 with respect to the structure core are populated. MD trajectories of the complex between the ribosomal protein S15 and RNA were also recorded in the presence and absence of Mg(2+) ions. The Mg(2+) ions are hexacoordinated by water and RNA oxygens. The coordination spheres mainly interact with the RNA phosphodiester backbone, reducing the RNA mobility and inducing electrostatic screening. When the Mg(2+) ions are removed, the internal mobility of the RNA and of the protein increases at the interaction interface close to the RNA G-U/G-C motif as a result of a gap between the phosphate groups in the UUCG capping tetraloop and of the disruption of S15-RNA hydrogen bonds in that region. On the other hand, several S15-RNA hydrogen bonds are reinforced, and water bridges appear between the three-way junction region and S15. The network of hydrogen bonds observed in the loop between alpha1 and alpha2 is consequently reorganized. In the absence of Mg(2+), this network has the same pattern as the network observed in the isolated protein, where the helix alpha1 is mobile with respect to the protein core. The presence of Mg(2+) ions may thus play a role in stabilizing the orientation of the helix alpha1 of S15.

Role of metal ions in the tetraloop-receptor complex as analyzed by NMR

Metal ions are critical for the proper folding of RNA, and the GAAA tetraloop-receptor is necessary for the optimal folding and function of many RNAs. We have used NMR to investigate the role of metal ions in the structure of the tetraloop-receptor in solution. The NMR data indicate native tertiary structure is formed under a wide range of ionic conditions. The lack of conformational adaptation in response to very different ionic conditions argues against a structural role for divalent ions. Nuclear Overhauser effects to cobalt hexammine and paramagnetic relaxation enhancement induced by manganese ions were used to determine the NMR structures of the tetraloop receptor in association with metal ions, providing the first atomic-level view of these interactions in the solution state. Five manganese and two cobalt hexammine ions could be localized to the RNA surface. The locations of the associated metal ions are similar, but not identical to, those of previously determined crystal structures. The sites of association are in general agreement with nonlinear Poisson-Boltzmann calculations of the electrostatic surface, emphasizing the general importance of diffusely associated ions in RNA tertiary structure.

Probing RNA hairpins with cobalt(III)hexammine and electrospray ionization mass spectrometry

In this work, electrospray ionization mass spectrometry (ESI MS) was employed to study the interactions of cobalt(III) hexammine, Co(NH3)6(3+), with five RNA hairpins representing the 790 loop of 16S ribosomal RNA and 1920 loop of 23S ribosomal RNA. The RNAs varied in mismatch identity (G.U versus A.C) and level of base modification (pseudouridine versus uridine). Co(NH3)6(3+) binding was observed with the four RNA hairpins that contained a G.U wobble pair in the stem region. ESI MS revealed 1:1 and 1:2 complex formation with all RNAs. Weaker binding was observed with the fifth RNA hairpin that contained an A.C wobble pair in the stem region. The effects of pH on Co(NH3)6(3+) binding were also examined.

Solution probing of metal ion binding by helix 27 from Escherichia coli 16S rRNA

Helix (H)27 from Escherichia coli 16S ribosomal (r)RNA is centrally located within the small (30S) ribosomal subunit, immediately adjacent to the decoding center. Bacterial 30S subunit crystal structures depicting Mg(2+) binding sites resolve two magnesium ions within the vicinity of H27: one in the major groove of the G886-U911 wobble pair, and one within the GCAA tetraloop. Binding of such metal cations is generally thought to be crucial for RNA folding and function. To ask how metal ion-RNA interactions in crystals compare with those in solution, we have characterized, using solution NMR spectroscopy, Tb(3+) footprinting and time-resolved fluorescence resonance energy transfer (tr-FRET), location, and modes of metal ion binding in an isolated H27. NMR and Tb(3+) footprinting data indicate that solution secondary structure and Mg(2+) binding are generally consistent with the ribosomal crystal structures. However, our analyses also suggest that H27 is dynamic in solution and that metal ions localize within the narrow major groove formed by the juxtaposition of the loop E motif with the tandem G894-U905 and G895-U904 wobble pairs. In addition, tr-FRET studies provide evidence that Mg(2+) uptake by the H27 construct results in a global lengthening of the helix. We propose that only a subset of H27-metal ion interactions has been captured in the crystal structures of the 30S ribosomal subunit, and that small-scale structural dynamics afforded by solution conditions may contribute to these differences. Our studies thus highlight an example for differences between RNA-metal ion interactions observed in solution and in crystals.

Antiviral 6-amino-quinolones: Molecular basis for potency and selectivity

Structural modifications introduced in 6-amino-quinolones to increase antiviral activity can strongly affect cytotoxicity to host cells. By competition to Tat-TAR complex and binding experiments to viral and cellular DNA and RNA structures, we show that the nature of the substituent at position 7 modifies drug affinity and specificity for the nucleic acid. Interestingly, the basicity of the above substituent modulates chelation of the quinolone template to magnesium ions, which, in turn, critically affects the potency and target selectivity in the antiviral quinolone family.

Ions and RNA Folding

The problem of how ions influence the folding of RNA into specific tertiary structures is being addressed from both thermodynamic (by how much do different salts affect the free energy change of folding) and structural (how are ions arranged on or near an RNA and what kinds of environments do they occupy) points of view. The challenge is to link these different approaches in a theoretical framework that relates the energetics of ion-RNA interactions to the spatial distribution of ions. This review distinguishes three different kinds of ion environments that differ in the extent of direct ion-RNA contacts and the degree to which the ion hydration is perturbed, and summarizes the current understanding of the way each environment relates to the overall energetics of RNA folding.

Toward multiplexing the application of solvent accessibility probes for the investigation of RNA three-dimensional structures by electrospray ionization-Fourier transform mass spectrometry

Multiple solvent accessibility probes can be applied simultaneously to investigate the three-dimensional structure of complex RNA substrates when electrospray ionization-Fourier transform mass spectrometry (ESI-FTMS) is employed in place of polyacrylamide gel electrophoresis (PAGE). We show that classic chemical probes, such as dimethylsulfate, kethoxal, and 1-cyclohexyl-3-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate, can be combined in probing mixtures designed to assess the full spectrum of base pairing and steric protection for the most abundant ribonucleotides included in RNA. After probe-independent hydrolysis of the alkylated substrate, the mixture of oligonucleotide products is mass mapped by ESI-FTMS analysis, which enables the unambiguous identification of probed bases from the unique mass signatures provided by the different chemical modifiers. In this bottom-up approach, any theoretical limit to the size of the possible target RNA will be determined by the effectiveness of the hydrolysis procedure rather than by the performance of the detection technique. Control experiments performed on the stem-loop 4 of human immunodeficiency virus type 1 have shown no adverse interactions between the reagents combined in the probing cocktails. No significant discrepancies between the alkylation patterns offered by the cocktails and the individual reagents could be detected, indicating that multiplexing the probe application does not necessarily lead to structural distortion but provides valid data on base accessibility and protection. To demonstrate the ruggedness of this approach, optimized cocktails were finally employed to assess the stability of the folded structure of mouse mammary tumor virus pseudoknot in the presence of different amounts of Mg2+. Multiplexing the probe application constitutes an essential step toward high-throughput applications, which will take advantage of a strategy that maximizes the information attainable from a single experiment, while minimizing time and sample consumption over PAGE-based methods.

Calculation of the binding free energy for magnesium-RNA interactions

The nature of the interaction between nucleic acids and divalent ions in solution is complex. It includes long-range electrostatic and short-range nonelectrostatic forces. Water molecules can be in an inner coordination shell that intervenes between the ion and its binding site. This work describes a method for calculating the binding free energy and applies it to a specific Mg-RNA system in the presence of monovalent salt. The approach combines high-level ab initio theory with Poisson-Boltzmann calculations and provides an accurate description of all terms of the binding free energy for magnesium ions located at the RNA surface (including nonelectrostatic interactions). Some alternative macroscopic approaches are also discussed.

Structural features and thermodynamics of the J4/5 loop from the Candida albicans and Candida dubliniensis group I introns

The J4/5 loop of group I introns has tertiary interactions with the P1 helix that position the P1 substrate for the self-splicing reaction. The J4/5 loop of Candida albicans and Candida dubliniensis, 5'GAAGG3'/3'UAAUU5', potentially contains two A.A pairs flanked by one G.U pair on one side and two G.U pairs on the other side. Results from optical melting, nuclear magnetic resonance spectroscopy, and functional group substitution experiments with a mimic of the C. albicans and C. dubliniensis J4/5 loop are consistent with the adenosines forming tandem sheared A.A pairs with a cross-strand stack and only the G.U pair not adjacent to an A.A pair forming a static wobble G.U pair. The two G.U pairs adjacent to the tandem A.A pairs are likely in a dynamic equilibrium between multiple conformations. Although Co(NH(3))(6)(3+) stabilizes the loop by several kilocalories per mole at 37 degrees C, addition of Mg(2+) or Co(NH(3))(6)(3+) has no effect on the structure of the loop. The tandem G.U pairs provide a pocket of negative charge for Co(NH(3))(6)(3+) to bind. The results contribute to understanding the structure and dynamics of purine-rich internal loops and potential G.U pairs adjacent to internal loops.

Untying the FIV frameshifting pseudoknot structure by MS3D

The structure of the putative feline immunodeficiency virus (FIV) ribosomal frameshifting pseudoknot (PK) has been investigated by a mass spectrometric three-dimensional (MS3D) approach, which involves the application of established solvent-accessibility probes and chemical crosslinkers with detection by electrospray ionization (ESI) Fourier transform mass spectrometry (FTMS). Regardless of their size, probed substrates can be treated with ribonucleases and analyzed by ESI-FTMS to obtain the correct position of chemically modified nucleotides. Protection maps and distance information can be utilized to generate 3D models using the constraint satisfaction algorithm provided by MC-SYM and the energy minimization modules included in CNS. Control experiments were performed on a mutant of mouse mammary tumor virus pseudoknot (VPK), for which an NMR structure is available. Comparison between the MS3D model and the high-resolution structure provided a approximately 3A root-mean-square deviation calculated from all the atoms present in double-stranded regions. Applied to FIV-PK, the MS3D approach confirmed that the selected sequence could fold into an actual pseudoknot, supporting the sequence alignment predictions. Characteristic features of H-type pseudoknots were recognized immediately, but a putative A13-U30 pair was not observed at the stem junction, making FIV-PK resemble VPK more closely than the initially suggested simian retrovirus type-1 pseudoknot. In our model, the unpaired U30 protrudes into the medium, while the hinging A13 assumes a stacked conformation that enables the stems to form a approximately 60 degrees bend and relieve the strain caused by a short loop 1. The model provided the basis to explain the different alkylation patterns observed in the absence and presence of Mg(2+), suggesting the possible formation of a specific metal-binding site between loop 1 and stem 2. This instance illustrates how the MS3D model of FIV-PK can be utilized effectively to generate hypotheses and support functional observations in the absence of a high-resolution structure.

Detection of RNA nucleobase metalation by NMR spectroscopy

Nucleic acids possess several metal cation recognition sites, including phosphates, nucleobases and possibly riboses. This article focuses on the detection of nucleobase-metal interactions by NMR spectroscopy.

NMR spectroscopy of RNA

NMR spectroscopy is a powerful tool for studying proteins and nucleic acids in solution. This is illustrated by the fact that nearly half of all current RNA structures were determined by using NMR techniques. Information about the structure, dynamics, and interactions with other RNA molecules, proteins, ions, and small ligands can be obtained for RNA molecules up to 100 nucleotides. This review provides insight into the resonance assignment methods that are the first and crucial step of all NMR studies, into the determination of base-pair geometry, into the examination of local and global RNA conformation, and into the detection of interaction sites of RNA. Examples of NMR investigations of RNA are given by using several different RNA molecules to illustrate the information content obtainable by NMR spectroscopy and the applicability of NMR techniques to a wide range of biologically interesting RNA molecules.

Selection and evolution of NTP-specific aptamers

ATP occupies a central position in biology, for it is both an elementary building block of RNA and the most widely used cofactor in all living organisms. For this reason, it has been a recurrent target for in vitro molecular evolution techniques. The exploration of ATP-binding motifs constitutes both an important step in investigating the plausibility of the 'RNA world' hypothesis and a central starting point for the development of new enzymes. To date, only two RNA motifs that bind ATP have been characterized. The first one is targeted to the adenosine moiety, while the second one recognizes the 'Hoogsteen' face of the base. To isolate aptamers that bind ATP in different orientations, we selected RNAs on an affinity resin that presents ATP in three different orientations. We obtained five new motifs that were characterized and subsequently submitted to a secondary selection protocol designed to isolate aptamers specific for cordycepin. Interestingly, all the ATP-binding motifs selected specifically recognize the sugar-phosphate backbone region of the nucleotides. Three of the aptamers show some selectivity for adenine derivatives, while the remainder recognize any of the four nucleotides with similar efficiency. The characteristics of these aptamers are discussed along with implications for in vitro molecular evolution.

Metal ion stabilization of the U-turn of the A37 N6-dimethylallyl-modified anticodon stem-loop of Escherichia coli tRNAPhe

Nucleoside base modifications can alter the structures, dynamics, and metal ion binding properties of transfer RNA molecules and are important for accurate aminoacylation and for maintaining translational fidelity and efficiency. The unmodified anticodon stem-loop from Escherichia coli tRNA(Phe) forms a trinucleotide loop in solution, but Mg(2+) and dimethylallyl modification of A(37) N6 disrupt the loop conformation and increase the mobility of the loop and loop-proximal nucleotides. We have used NMR spectroscopy to investigate the binding and structural effects of multivalent cations on the unmodified and dimethylallyl-modified anticodon stem-loops from E. coli tRNA(Phe). The divalent cation binding sites were probed using Mn(2+) and Co(NH(3))(6)(3+). These ions bind along the major groove of the stem and associate with the anticodon loop on the major groove side in a nonspecific manner. Co(NH(3))(6)(3+) stabilizes the U-turn conformation of the loop in the dimethylallyl-modified molecule, and the chemical shift changes that accompany Co(NH(3))(6)(3+) binding are similar to those observed with the addition of Mg(2+). The base-phosphate and base-2'-OH hydrogen bonds that characterize the UNR U-turn motif lead to spectral signatures in the form of unusual (15)N and (1)H chemical shifts and reduced solvent exchange of the U(33) 2'-OH and N3H protons. The unmodified molecule also displays spectral features of the U-turn fold in the presence of Co(NH(3))(6)(3+), but the loop has additional conformations and is dynamic. The results indicate that charge neutralization by a polyvalent cation is sufficient to promote formation of the U-turn fold. However, base modification is necessary to destabilize competing alternative conformers even for a purine-rich loop sequence that is predicted to have strongly favorable base stacking energy.

Roles of protein subunits in RNA-protein complexes: lessons from ribonuclease P

Ribonucleoproteins (RNP) are involved in many essential processes in life. However, the roles of RNA and protein subunits in an RNP complex are often hard to dissect. In many RNP complexes, including the ribosome and the Group II introns, one main function of the protein subunits is to facilitate RNA folding. However, in other systems, the protein subunits may perform additional functions, and can affect the biological activities of the RNP complexes. In this review, we use ribonuclease P (RNase P) as an example to illustrate how the protein subunit of this RNP affects different aspects of catalysis. RNase P plays an essential role in the processing of the precursor to transfer RNA (pre-tRNA) and is found in all three domains of life. While every cell has an RNase P (ribonuclease P) enzyme, only the bacterial and some of the archaeal RNase P RNAs (RNA component of RNase P) are active in vitro in the absence of the RNase P protein. RNase P is a remarkable enzyme in the fact that it has a conserved catalytic core composed of RNA around which a diverse array of protein(s) interact to create the RNase P holoenzyme. This combination of highly conserved RNA and altered protein components is a puzzle that allows the dissection of the functional roles of protein subunits in these RNP complexes.

In vitro characterization of a base pairing interaction between the primer binding site and the minimal packaging signal of avian leukosis virus genomic RNA

The 5' leader region of avian sarcoma-leukosis viruses (ASLVs) folds into a series of RNA secondary structures which are involved in key steps in the viral replication cycle such as reverse transcription, dimerization and packaging of genomic RNA. The O3 stem and three stem-loops (O3SLa, O3SLb and O3SLc) form the minimal packaging signal that is located downstream of the primer binding site (PBS). The U5-PBS region contributes to packaging via a mechanism that remains unknown. In this in vitro study, we have investigated the possibility of interactions between the R-U5-PBS region and the minimal packaging signal using chemical and enzymatic probing, antisense oligonucleotides and site-directed mutagenesis. We have identified a base pairing interaction between the PBS sequence and the terminal loop of O3SLa. It was found that the PBS/O3SLa interaction was intramolecular since it occurred not only in dimeric RNA but also in monomeric RNA. This interaction probably corresponds to a pseudoknot interaction. The PBS/O3SLa interaction may be formed in vivo since the sequences are highly conserved in ASLV strains. The PBS/O3SLa interaction may regulate the processes of primer tRNA annealing, packaging and initiation of Gag translation through its involvement in leader tertiary structure. Interestingly, we found that in other retroviruses the PBS sequence can also base pair with a terminal loop of the stem-loops involved in RNA packaging.

Theoretical Calculation of the NMR Spin−Spin Coupling Constants and the NMR Shifts Allow Distinguishability between the Specific Direct and the Water-Mediated Binding of a Divalent Metal Cation to Guanine

The calculated intermolecular and intramolecular indirect NMR spin-spin coupling constants and NMR shifts were used for the discrimination between the inner-shell and the outer-shell binding motif of hydrated divalent cations Mg(2+) or Zn(2+) with a guanine base. The intermolecular coupling constants (1)J(X,O6) and (1)J(X,N7) (X = Mg(2+), Zn(2+)) can be unambiguously assigned to the specific inner-shell binding motif of the hydrated cation either with oxygen O6 or with nitrogen N7 of guanine. The calculated coupling constants (1)J(Mg,O6) and (1)J(Zn,O6) were 6.2 and -17.5 Hz, respectively, for the inner-shell complex of cation directly interacting with oxygen O6 of guanine. For the inner-shell coordination of the cation at nitrogen N7, the calculated coupling constants (1)J(Mg,N7) and (1)J(Zn,N7) were 5.6 and -36.5 Hz, respectively. When the binding of the cation is water-mediated, the coupling constant is zero. To obtain reliable shifts in NMR parameters, hydrated guanine was utilized as the reference state. The calculated change of NMR spin-spin coupling constants due to the hydration and coordination of the cation with guanine is caused mainly by the variation of Fermi-contact coupling contribution while the variation of diamagnetic spin-orbit, paramagnetic spin-orbit, and spin-dipolar coupling contributions is small. The change of s-character of guanine sigma bonding, sigma antibonding, and lone pair orbitals upon the hydration and cation coordination (calculated using the Natural Bond Orbital analysis) correlates with the variation of the Fermi-contact term. The calculated NMR shifts delta(N7) of -15.3 and -12.2 ppm upon the coordination of Mg(2+) and Zn(2+) ion are similar to the NMR shift of 19.6 ppm toward the high field measured by Tanaka for N7 of guanine upon the coordination of the Cd(2+) cation (Tanaka, Y.; Kojima, C.; Morita, E. H.; Kasai. Y.; Yamasaki, K.; Ono, A.; Kainosho, M.; Taira, K. J. Am. Chem. Soc. 2002, 124, 4595-4601). The present data indicate that measurements of NMR intermolecular coupling constants may be used to discriminate between the specific inner- and outer-shell binding of divalent cations to nucleobases in DNA and RNA.

An eight-stranded helical fragment in RNA crystal structure: implications for tetraplex interaction

Multistranded helical structures in nucleic acids play various functions in biological processes. Here we report the crystal structure of a hexamer, rU(BrdG)r(AGGU),at 1.5 A resolution containing a structural complex of an alternating antiparallel eight-stranded helical fragment that is sandwiched in two tetraplexes. The octaplex is formed by groove binding interaction and base tetrad intercalation between two tetraplexes. Two different forms of octaplexes have been proposed, which display different properties in interaction with proteins and nucleic acids. Adenines form a base tetrad in the novel N6-H em leader N3 conformation and further interact with uridines to form an adenine-uridine octad in the reverse Hoogsteen pairing scheme. The conformational flexibility of adenine tetrad indicates that it can optimize its conformation in different interactions.

NMR structure of the active conformation of the Varkud satellite ribozyme cleavage site

Substrate cleavage by the Neurospora Varkud satellite (VS) ribozyme involves a structural change in the stem-loop I substrate from an inactive to an active conformation. We have determined the NMR solution structure of a mutant stem-loop I that mimics the active conformation of the cleavage site internal loop. This structure shares many similarities, but also significant differences, with the previously determined structures of the inactive internal loop. The active internal loop displays different base-pairing interactions and forms a novel RNA fold composed exclusively of sheared G-A base pairs. From chemical-shift mapping we identified two Mg2+ binding sites in the active internal loop. One of the Mg2+ binding sites forms in the active but not the inactive conformation of the internal loop and is likely important for catalysis. Using the structure comparison program mc-search, we identified the active internal loop fold in other RNA structures. In Thermus thermophilus 16S rRNA, this RNA fold is directly involved in a long-range tertiary interaction. An analogous tertiary interaction may form between the active internal loop of the substrate and the catalytic domain of the VS ribozyme. The combination of NMR and bioinformatic approaches presented here has identified a novel RNA fold and provides insights into the structural basis of catalytic function in the Neurospora VS ribozyme.