A crystallographic study of metal-binding to yeast phenylalanine transfer RNA

RNase P cleavage of pseudoknot substrates reveals differences in active site architecture that depend on residue N-1 in the 5' leader

We show that a small biotin-binding RNA aptamer that folds into a pseudoknot structure acts as a substrate for bacterial RNase P RNA (RPR) with and without the RNase P C5 protein. Cleavage in the single-stranded region in loop 1 was shown to depend on the presence of a RCCA-motif at the 3' end of the substrate. The nucleobase and the 2'hydroxyl at the position immediately 5' of the cleavage site contribute to both cleavage efficiency and site selection, where C at this position induces significant cleavage at an alternative site, one base upstream of the main cleavage site. The frequencies of cleavage at these two sites and Mg2+ binding change upon altering the structural topology in the vicinity of the cleavage site as well as by replacing Mg2+ with other divalent metal ions. Modelling studies of RPR in complex with the pseudoknot substrates suggest alternative structural topologies for cleavage at the main and the alternative site and a shift in positioning of Mg2+ that activates the H2O nucleophile. Together, our data are consistent with a model where the organization of the active site structure and positioning of Mg2+ is influenced by the identities of residues at and in the vicinity of the site of cleavage.

Unraveling the Structure-Spectrum Relationship of Yeast Phenylalanine Transfer RNA: Insights from Theoretical Modeling of Infrared Spectroscopy

Yeast phenylalanine tRNA (tRNAphe) is a paradigmatic model in structural biology. In this work, we combine molecular dynamics simulations and spectroscopy modeling to establish a direct link between its structure, conformational dynamics, and infrared (IR) spectra. Employing recently developed vibrational frequency maps and coupling models, we apply a mixed quantum/classical treatment of the line shape theory to simulate the IR spectra of tRNAphe in the 1600-1800 cm-1 region across its folded and unfolded conformations and under varying concentrations of Mg2+ ions. The predicted IR spectra of folded and unfolded tRNAphe are in good agreement with experimental measurements, validating our theoretical framework. We then elucidate how the characteristic L-shaped tertiary structure of the tRNA and its modulation in response to diverse chemical environments give rise to distinct IR absorption peaks and line shapes. These calculations effectively bridge IR spectroscopy experiments and atomistic molecular simulations, unraveling the molecular origins of the observed IR spectra of tRNAphe. This work presents a robust theoretical protocol for modeling the IR spectroscopy of nucleic acids, which will facilitate its application as a sensitive probe for detecting the fluctuating secondary and tertiary structures of these essential biological macromolecules.

Structural and Genetic Determinants of Convergence in the Drosophila tRNA Structure-Function Map

The evolution of tRNA multigene families remains poorly understood, exhibiting unusual phenomena such as functional conversions of tRNA genes through anticodon shift substitutions. We improved FlyBase tRNA gene annotations from twelve Drosophila species, incorporating previously identified ortholog sets to compare substitution rates across tRNA bodies at single-site and base-pair resolution. All rapidly evolving sites fell within the same metal ion-binding pocket that lies at the interface of the two major stacked helical domains. We applied our tRNA Structure-Function Mapper (tSFM) method independently to each Drosophila species and one outgroup species Musca domestica and found that, although predicted tRNA structure-function maps are generally highly conserved in flies, one tRNA Class-Informative Feature (CIF) within the rapidly evolving ion-binding pocket-Cytosine 17 (C17), ancestrally informative for lysylation identity-independently gained asparaginylation identity and substituted in parallel across tRNAAsn paralogs at least once, possibly multiple times, during evolution of the genus. In D. melanogaster, most tRNALys and tRNAAsn genes are co-arrayed in one large heterologous gene cluster, suggesting that heterologous gene conversion as well as structural similarities of tRNA-binding interfaces in the closely related asparaginyl-tRNA synthetase (AsnRS) and lysyl-tRNA synthetase (LysRS) proteins may have played a role in these changes. A previously identified Asn-to-Lys anticodon shift substitution in D. ananassae may have arisen to compensate for the convergent and parallel gains of C17 in tRNAAsn paralogs in that lineage. Our results underscore the functional and evolutionary relevance of our tRNA structure-function map predictions and illuminate multiple genomic and structural factors contributing to rapid, parallel and compensatory evolution of tRNA multigene families.

Benefits and Detriments of Gadolinium from Medical Advances to Health and Ecological Risks

Gadolinium (Gd)-containing chelates have been established as diagnostics tools. However, extensive use in magnetic resonance imaging has led to increased Gd levels in industrialized parts of the world, adding to natural occurrence and causing environmental and health concerns. A vast amount of data shows that metal may accumulate in the human body and its deposition has been detected in organs such as brain and liver. Moreover, the disease nephrogenic systemic fibrosis has been linked to increased Gd³⁺ levels. Investigation of Gd³⁺ effects at the cellular and molecular levels mostly revolves around calcium-dependent proteins, since Gd³⁺ competes with calcium due to their similar size; other reports focus on interaction of Gd³⁺ with nucleic acids and carbohydrates. However, little is known about Gd³⁺ effects on membranes; yet some results suggest that Gd³⁺ interacts strongly with biologically-relevant lipids (e.g., brain membrane constituents) and causes serious structural changes including enhanced membrane rigidity and propensity for lipid fusion and aggregation at much lower concentrations than other ions, both toxic and essential. This review surveys the impact of the anthropogenic use of Gd emphasizing health risks and discussing debilitating effects of Gd³⁺ on cell membrane organization that may lead to deleterious health consequences.

Nucleobase carbonyl groups are poor Mg2+ inner-sphere binders but excellent monovalent ion binders-a critical PDB survey

Precise knowledge of Mg²⁺ inner-sphere binding site properties is vital for understanding the structure and function of nucleic acid systems. Unfortunately, the PDB, which represents the main source of Mg²⁺ binding sites, contains a substantial number of assignment issues that blur our understanding of the functions of these ions. Here, following a previous study devoted to Mg²⁺ binding to nucleobase nitrogens, we surveyed nucleic acid X-ray structures from the PDB with resolutions ≤2.9 Å to classify the Mg²⁺ inner-sphere binding patterns to nucleotide carbonyl, ribose hydroxyl, cyclic ether, and phosphodiester oxygen atoms. From this classification, we derived a set of "prior-knowledge" nucleobase Mg²⁺ binding sites. We report that crystallographic examples of trustworthy nucleobase Mg²⁺ binding sites are fewer than expected since many of those are associated with misidentified Na⁺ or K⁺ We also emphasize that binding of Na⁺ and K⁺ to nucleic acids is much more frequent than anticipated. Overall, we provide evidence derived from X-ray structures that nucleobases are poor inner-sphere binders for Mg²⁺ but good binders for monovalent ions. Based on strict stereochemical criteria, we propose an extended set of guidelines designed to help in the assignment and validation of ions directly contacting nucleobase and ribose atoms. These guidelines should help in the interpretation of X-ray and cryo-EM solvent density maps. When borderline Mg²⁺ stereochemistry is observed, alternative placement of Na⁺, K⁺, or Ca²⁺ must be considered. We also critically examine the use of lanthanides (Yb³⁺, Tb³⁺) as Mg²⁺ substitutes in crystallography experiments.

Biochemical and Biophysical Understanding of Metal Ion Selectivity of DNAzymes

This review summarizes research into the metal-binding properties of catalytic DNAzymes, towards the goal of understanding the structural properties leading to metal ion specificity. Progress made and insight gained from a range of biochemical and biophysical techniques are covered, and promising directions for future investigations are discussed.

Understanding nucleic acid-ion interactions

Ions surround nucleic acids in what is referred to as an ion atmosphere. As a result, the folding and dynamics of RNA and DNA and their complexes with proteins and with each other cannot be understood without a reasonably sophisticated appreciation of these ions' electrostatic interactions. However, the underlying behavior of the ion atmosphere follows physical rules that are distinct from the rules of site binding that biochemists are most familiar and comfortable with. The main goal of this review is to familiarize nucleic acid experimentalists with the physical concepts that underlie nucleic acid-ion interactions. Throughout, we provide practical strategies for interpreting and analyzing nucleic acid experiments that avoid pitfalls from oversimplified or incorrect models. We briefly review the status of theories that predict or simulate nucleic acid-ion interactions and experiments that test these theories. Finally, we describe opportunities for going beyond phenomenological fits to a next-generation, truly predictive understanding of nucleic acid-ion interactions.

The crystal structure of unmodified tRNAPhe from Escherichia coli

Post-transcriptional nucleoside modifications fine-tune the biophysical and biochemical properties of transfer RNA (tRNA) so that it is optimized for participation in cellular processes. Here we report the crystal structure of unmodified tRNA^Phe from Escherichia coli at a resolution of 3 Å. We show that in the absence of modifications the overall fold of the tRNA is essentially the same as that of mature tRNA. However, there are a number of significant structural differences, such as rearrangements in a triplet base pair and a widened angle between the acceptor and anticodon stems. Contrary to previous observations, the anticodon adopts the same conformation as seen in mature tRNA.

Nature of guanine oxidation in RNA via the flash-quench technique versus direct oxidation by a metal oxo complex

Oxidation of RNA can be effected by two different techniques: a photochemical, electron-transfer method termed "flash-quench" and direct oxidation by metal oxo complexes. The flash-quench method produces selective oxidation using a metal photosensitizer, tris(bipyridyl)ruthenium(III) trichloride (Ru(bpy)(3)(3+)), and quencher, pentaamminechlorocobalt(III) chloride (Co(NH(3))(5)Cl(2+)). We have optimized the flash-quench technique for the following RNAs: tRNA(Phe), human ferritin iron-responsive element (IRE), and a mutated human ferritin IRE. We have also employed a chemical footprinting technique involving the oxoruthenium(IV) complex (Ru(tpy)(bpy)O(2+) (tpy = 2,2',2''-terpyridine; bpy = 2,2'-bipyridine)) to oxidize guanine. Comparison of the two methods shows that the flash-quench technique provides a visualization of nucleotide accessibility for a static conformation of RNA while the Ru(tpy)(bpy)O(2+) complex selectively oxidizes labile guanines and gives a visualization of a composite of multiple conformations of the RNA structure.

Frequency and effect of the binding of Mg2+, Mn2+, and Co2+ ions on the guanine base in Watson-Crick and reverse Watson-Crick base pairs

We performed MP2 calculations to elucidate the structure and energetics of the Mg(2+), Mn(2+), and Co(2+) hexahydrated aquaions, and the effect of the metal binding to the N7 atom of (i) a single guanine, (ii) a guanine involved in a Watson-Crick pair, and (iii) a guanine involved in a reverse Watson-Crick base pair. Our comparative analysis of the three aquaions indicates a clear inverse correlation between the radius of the cation and the binding energy, that indeed increases in the order Mn(2+) < Co(2+) < Mg(2+). The trend in the binding energies of the pentahydrated cations to the N7 atom of the guanine is instead Mg(2+) < Mn(2+) < Co(2+), suggesting a rather different bonding scheme that, for the two transition metals, involves back-donation from the aromatic ring of the guanine to their empty d orbitals. In the gas phase, the three hydrated metals significantly stabilize both G-C base pair geometries, Watson-Crick and reverse Watson-Crick, we investigated. Inclusion of a continuous solvent model, however, remarkably reduces this additional stabilization, which becomes almost negligible in the case of the Mg(2+) cation coordinated to the guanine in the standard Watson-Crick geometry. Conversely, all three metal ions sensibly stabilize the reverse Watson-Crick geometry, also in water. Our results are supported by a screening of the structures available in the Protein Data Bank, which clearly indicates that the two transition metals we investigated have a tendency greater than Mg(2+) to coordinate to the N7 atom of guanines, and that there is no clear correlation between the number of guanines in experimental structures with a metal bound to N7 atom and their involvement in Watson-Crick base pairs.

Analyzing the flexibility of RNA structures by constraint counting

RNA requires conformational dynamics to undergo its diverse functional roles. Here, a new topological network representation of RNA structures is presented that allows analyzing RNA flexibility/rigidity based on constraint counting. The method extends the FIRST approach, which identifies flexible and rigid regions in atomic detail in a single, static, three-dimensional molecular framework. Initially, the network rigidity of a canonical A-form RNA is analyzed by counting on constraints of network elements of increasing size. These considerations demonstrate that it is the inclusion of hydrophobic contacts into the RNA topological network that is crucial for an accurate flexibility prediction. The counting also explains why a protein-based parameterization results in overly rigid RNA structures. The new network representation is then validated on a tRNA(ASP) structure and all NMR-derived ensembles of RNA structures currently available in the Protein Data Bank (with chain length >/=40). The flexibility predictions demonstrate good agreement with experimental mobility data, and the results are superior compared to predictions based on two previously used network representations. Encouragingly, this holds for flexibility predictions as well as mobility predictions obtained by constrained geometric simulations on these networks. Potential applications of the approach to analyzing the flexibility of DNA and RNA/protein complexes are discussed.

The 1.3 A resolution structure of the RNA tridecamer r(GCGUUUGAAACGC): metal ion binding correlates with base unstacking and groove contraction

Metal ions play a key role in RNA folding and activity. Elucidating the rules that govern the binding of metal ions is therefore an essential step for better understanding the RNA functions. High-resolution data are a prerequisite for a detailed structural analysis of ion binding on RNA and, in particular, the observation of monovalent cations. Here, the high-resolution crystal structures of the tridecamer duplex r(GCGUUUGAAACGC) crystallized under different conditions provides new structural insights on ion binding on GAAA/UUU sequences that exhibit both unusual structural and functional properties in RNA. The present study extends the repertory of RNA ion binding sites in showing that the two first bases of UUU triplets constitute a specific site for sodium ions. A striking asymmetric pattern of metal ion binding in the two equivalent halves of the palindromic sequence demonstrates that sequence and its environment act together to bind metal ions. A highly ionophilic half that binds six metal ions allows, for the first time, the observation of a disodium cluster in RNA. The comparison of the equivalent halves of the duplex provides experimental evidences that ion binding correlates with structural alterations and groove contraction.

Conformational energy and structure in canonical and noncanonical forms of tRNA determined by temperature analysis of the rate of s(4)U8-C13 photocrosslinking

Bacterial tRNAs frequently have 4-thiouridine (s(4)U) modification at position 8, which is adjacent to the C13-G22-m(7)G46 base triple in the elbow region of the tRNA tertiary structure. Irradiation with light in the UVA range induces an efficient photocrosslink between s(4)U8 and C13. The temperature dependence of the rate constants for photocrosslinking between the s(4)U8 and C13 has been used to investigate the tRNA conformational energy and structure in Escherichia coli tRNA(Val), tRNA(Phe), and tRNA(fMet) under different conditions. Corrections have been made in the measured rate constants to compensate for differences in the excited state lifetimes due to tRNA identity, buffer conditions, and temperature. The resulting rate constants are related to the rate at which the s(4)U8 and C13 come into the alignment needed for photoreaction; this depends on an activation energy, attributable to the conformational potential energy that occurs during the photoreaction, and on the extent of the structural change. Different photocrosslinking rate constants and temperature dependencies occur in the three tRNAs, and these differences are due both to modest differences in the activation energies and in the apparent s(4)U8-C13 geometries. Analysis of tRNA(Val) in buffers without Mg(2+) indicate a smaller activation energy (~13 kJ mol(-1)) and a larger apparent s(4)U8-C13 distance (~12 A) compared to values for the same parameters in buffers with Mg(2+) (~26 kJ mol(-1) and 0.36 A, respectively). These measurements are a quantitative indication of the strong constraint that Mg(2+) imposes on the tRNA flexibility and structure.

Mg2+ binding and archaeosine modification stabilize the G15 C48 Levitt base pair in tRNAs

The G15-C48 Levitt base pair, located at a crucial position in the core of canonical tRNAs, assumes a reverse Watson-Crick (RWC) geometry. By means of bioinformatics analysis and quantum mechanics calculations we show here that such a geometry is moderately more stable than an alternative bifurcated trans geometry, involving the guanine Watson-Crick face and the cytosine keto group, which we have also found in known RNA structures. However we also demonstrate that the RWC geometry can take advantage of additional stabilizing effects such as metal binding or post-transcriptional chemical modification. One of the few strong metal binding sites characterized for cytosolic tRNAs is localized on G15, and a domain-specific complex modification known as archaeosine is widespread at position 15 in archaeal tRNAs. We have found that both the bound Mg2+ ion and the archaeosine modification induce an analogous electron density redistribution, which results in an effective stabilization of the RWC geometry. Metal binding and chemical modification thus play an interchangeable role in stabilizing the G15-C48 correct geometry. Interestingly, these different but convergent strategies are selectively adopted in the different life domains.

Structural metals in the group I intron: a ribozyme with a multiple metal ion core

Metal ions play key roles in the folding and function for many structured RNAs, including group I introns. We determined the X-ray crystal structure of the Azoarcus bacterial group I intron in complex with its 5' and 3' exons. In addition to 222 nucleotides of RNA, the model includes 18 Mg(2+) and K(+) ions. Five of the metals bind within 12 A of the scissile phosphate and coordinate the majority of the oxygen atoms biochemically implicated in conserved metal-RNA interactions. The metals are buried deep within the structure and form a multiple metal ion core that is critical to group I intron structure and function. Eight metal ions bind in other conserved regions of the intron structure, and the remaining five interact with peripheral structural elements. Each of the 18 metals mediates tertiary interactions, facilitates local bends in the sugar-phosphate backbone or binds in the major groove of helices. The group I intron has a rich history of biochemical efforts aimed to identify RNA-metal ion interactions. The structural data are correlated to the biochemical results to further understand the role of metal ions in group I intron structure and function.

Use of a novel Förster resonance energy transfer method to identify locations of site-bound metal ions in the U2-U6 snRNA complex

U2 and U6 snRNAs pair to form a phylogenetically conserved complex at the catalytic core of the spliceosome. Interactions with divalent metal ions, particularly Mg(II), at specific sites are essential for its folding and catalytic activity. We used a novel Förster resonance energy transfer (FRET) method between site-bound luminescent lanthanide ions and a covalently attached fluorescent dye, combined with supporting stoichiometric and mutational studies, to determine locations of site-bound Tb(III) within the human U2-U6 complex. At pH 7.2, we detected three metal-ion-binding sites in: (1) the consensus ACACAGA sequence, which forms the internal loop between helices I and III; (2) the four-way junction, which contains the conserved AGC triad; and (3) the internal loop of the U6 intra-molecular stem loop (ISL). Binding at each of these sites is supported by previous phosphorothioate substitution studies and, in the case of the ISL site, by NMR. Binding of Tb(III) at the four-way junction and the ISL sites was found to be pH-dependent, with no ion binding observed below pH 6 and 7, respectively. This pH dependence of metal ion binding suggests that the local environment may play a role in the binding of metal ions, which may impact on splicing activity.

MeRNA: a database of metal ion binding sites in RNA structures

Metal ions are essential for the folding of RNA into stable tertiary structures and for the catalytic activity of some RNA enzymes. To aid in the study of the roles of metal ions in RNA structural biology, we have created MeRNA (Metals in RNA), a comprehensive compilation of all metal binding sites identified in RNA 3D structures available from the PDB and Nucleic Acid Database. Currently, our database contains information relating to binding of 9764 metal ions corresponding to 23 distinct elements, in 256 RNA structures. The metal ion locations were confirmed and ligands characterized using original literature references. MeRNA includes eight manually identified metal-ion binding motifs, which are described in the literature. MeRNA is searchable by PDB identifier, metal ion, method of structure determination, resolution and R-values for X-ray structure and distance from metal to any RNA atom or to water. New structures with their respective binding motifs will be added to the database as they become available. The MeRNA database will further our understanding of the roles of metal ions in RNA folding and catalysis and have applications in structural and functional analysis, RNA design and engineering. The MeRNA database is accessible at .

A nature of conformational changes of yeast tRNAPhe

We analysed conformational changes of yeast tRNA(Phe) induced by high hydrostatic pressure (HHP) measured by Fourier-transform infrared (FTIR) and fluorescence spectroscopies. High pressure influences RNA conformation without other cofactors, such as metal ions and salts. FTIR spectra of yeast tRNA(Phe) recorded at high hydrostatic pressure up to 13 kbar with and without magnesium ions showed a shift of the bands towards higher frequencies. That blue shift is due to an increase a higher energy of bonds as a result of shortening of hydrogen bonds followed by dehydration of tRNA. The fluorescence spectra of Y-base tRNA(Phe) at high pressure up to 3 kbar showed a decrease of the intensity band at 430 nm as a consequence of conformational rearrangement of the anticodon loop leading to exposure of Y-base side chain to the solution. We suggest that structural transition of nucleic acids is driven by the changes of water structure from tetrahedral to a cubic-like geometry induced by high pressure and, in consequence, due to economy of hydration.

The structure-function dilemma of the hammerhead ribozyme

A powerful approach to understanding protein enzyme catalysis is to examine the structural context of essential amino acid side chains whose deletion or modification negatively impacts catalysis. In principle, this approach can be even more powerful for RNA enzymes, given the wide variety and subtlety of functionally modified nucleotides now available. Numerous recent success stories confirm the utility of this approach to understanding ribozyme function. An anomaly, however, is the hammerhead ribozyme, for which the structural and functional data do not agree well, preventing a unifying view of its catalytic mechanism from emerging. To delineate the hammerhead structure-function comparison, we have evaluated and distilled the large body of biochemical data into a consensus set of functional groups unambiguously required for hammerhead catalysis. By examining the context of these functional groups within available structures, we have established a concise set of disagreements between the structural and functional data. The number and relative distribution of these inconsistencies throughout the hammerhead reaffirms that an extensive conformational rearrangement from the fold observed in the crystal structure must be necessary for cleavage to occur. The nature and energetic driving force of this conformational isomerization are discussed.

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.

Characterization of structure and metal ions specificity of Co2+-binding RNA aptamers

Studies on RNA motifs capable of binding metal ions have largely focused on Mg(2+)-specific motifs, therefore information concerning interactions of other metal ions with RNA is still very limited. Application of the in vitro selection approach allowed us to isolate two RNA aptamers that bind Co(2+) ions. Structural analysis of their secondary structures revealed the presence of two motifs, loop E and "kissing" loop complex, commonly occurring in RNA molecules. The Co(2+)-induced cleavage method was used for identification of Co(2+)-binding sites after the determination of the optimal cleavage conditions. In the aptamers, Co(2+) ions seem to bind to N7 atoms of purines, inducing cleavage of the adjacent phosphodiester bonds, similarly as is the case with yeast tRNA(Phe). Although the in vitro selection experiment was carried out in the presence of Co(2+) ions only, the aptamers displayed broader metal ions specificity. This was shown by inhibition of Co(2+)-induced cleavages in the presence of the following transition metal ions: Zn(2+), Cd(2+), Ni(2+), and Co(NH(3))(6)(3+) complex. On the other hand, alkaline metal ions such as Mg(2+), Ca(2+), Sr(2+), and Ba(2+) affected Co(2+)-induced cleavages only slightly. Multiple metal ions specificity of Co(2+)-binding sites has also been reported for other in vitro selected or natural RNAs. Among many factors that influence metal specificity of the Co(2+)-binding pocket, chemical properties of metal ions, such as their hardness as well as the structure of the coordination site, seem to be particularly important.

Cross talk between the +73/294 interaction and the cleavage site in RNase P RNA mediated cleavage

To monitor functionally important metal ions and possible cross talk in RNase P RNA mediated cleavage we studied cleavage of substrates, where the 2'OH at the RNase P cleavage site (at -1) and/or at position +73 had been replaced with a 2' amino group (or 2'H). Our data showed that the presence of 2' modifications at these positions affected cleavage site recognition, ground state binding of substrate and/or rate of cleavage. Cleavage of 2' amino substituted substrates at different pH showed that substitution of Mg2+ by Mn2+ (or Ca2+), identity of residues at and near the cleavage site, and addition of C5 protein influenced the frequency of miscleavage at -1 (cleavage at the correct site is referred to as +1). From this we infer that these findings point at effects mediated by protonation/deprotonation of the 2' amino group, i.e. an altered charge distribution, at the site of cleavage. Moreover, our data suggested that the structural architecture of the interaction between the 3' end of the substrate and RNase P RNA influence the charge distribution at the cleavage site as well as the rate of cleavage under conditions where the chemistry is suggested to be rate limiting. Thus, these data provide evidence for cross talk between the +73/294 interaction and the cleavage site in RNase P RNA mediated cleavage. We discuss the role metal ions might play in this cross talk and the likelihood that at least one functionally important metal ion is positioned in the vicinity of, and use the 2'OH at the cleavage site as an inner or outer sphere ligand.

The origins of RNA catalysis in ribozymes

The discovery of RNA catalysis provided a paradigm shift in biology, insight into the evolution of life on the planet and a challenge to understand its mechanistic origins. RNA has limited catalytic resources that must be used to maximal effect. Consequently, RNA catalysis tends to be multifactorial, with several processes contributing to an overall significant enhancement of reaction rate. These include general acid-base catalysis, electrostatic effects, and substrate orientation and proximity. The main players are the RNA nucleobases and bound metal ions. Although most ribozymes carry out phosphoryl transfer, the same considerations appear to apply to peptidyl transfer in the ribosome.

A crystallographic study of the binding of 13 metal ions to two related RNA duplexes

Metal ions, and magnesium in particular, are known to be involved in RNA folding by stabilizing secondary and tertiary structures, and, as cofactors, in RNA enzymatic activity. We have conducted a systematic crystallographic analysis of cation binding to the duplex form of the HIV-1 RNA dimerization initiation site for the subtype-A and -B natural sequences. Eleven ions (K+, Pb2+, Mn2+, Ba2+, Ca2+, Cd2+, Sr2+, Zn2+, Co2+, Au3+ and Pt4+) and two hexammines [Co (NH3)6]3+ and [Ru (NH3)6]3+ were found to bind to the DIS duplex structure. Although the two sequences are very similar, strong differences were found in their cation binding properties. Divalent cations bind almost exclusively, as Mg2+, at 'Hoogsteen' sites of guanine residues, with a cation-dependent affinity for each site. Notably, a given cation can have very different affinities for a priori equivalent sites within the same molecule. Surprisingly, none of the two hexammines used were able to efficiently replace hexahydrated magnesium. Instead, [Co (NH3)4]3+ was seen bound by inner-sphere coordination to the RNA. This raises some questions about the practical use of [Co (NH3)6]3+ as a [Mg (H2O)6]2+ mimetic. Also very unexpected was the binding of the small Au3+ cation exactly between the Watson-Crick sites of a G-C base pair after an obligatory deprotonation of N1 of the guanine base. This extensive study of metal ion binding using X-ray crystallography significantly enriches our knowledge on the binding of middleweight or heavy metal ions to RNA, particularly compared with magnesium.

The bulge region of HIV-1 TAR RNA binds metal ions in solution

Binding of Mg2+, Ca2+ and Co(NH3)6(3+) ions to the HIV-1 TAR RNA in solution was analysed by 19F NMR spectroscopy, metal ion-induced RNA cleavages and Brownian dynamics (BD) simulations. Chemically synthesised 29mer oligoribonucleotides of the TAR sequence labelled with 5-fluorouridine (FU) were used for 19F NMR-monitored metal ion titration. The chemical shift changes of fluorine resonances FU-23, FU-25 and FU-40 upon titration with Mg2+ and Ca2+ ions indicated specific, although weak, binding at the bulge region with the dissociation constants (K(d)) of 0.9 +/- 0.6 and 2.7 +/- 1.7 mM, respectively. Argininamide, inducing largest (19)F chemical shifts changes at FU-23, was used as a reference ligand (K(d) = 0.3 +/- 0.1 mM). In the Pb2+-induced TAR RNA cleavage experiment, strong and selective cleavage of the C24-U25 phosphodiester bond was observed, while Mg2+ and Ca2+ induced cuts at all 3-nt residues of the bulge. The inhibition of Pb2+-specific TAR cleavage by di- and trivalent metal ions revealed a binding specificity [in the order Co(NH3)6(3+) > Mg2+ > Ca2+] at the bulge site. A BD simulation search of potential magnesium ion sites within the NMR structure of HIV-1 TAR RNA was conducted on a set of 20 conformers (PDB code 1ANR). For most cases, the bulge region was targeted by magnesium cations.

Kinetic and binding analysis of the catalytic involvement of ribose moieties of a trans-acting delta ribozyme

We have identified ribose 2'-hydroxyl groups (2'-OHs) that are critical for the activity of a trans-cleaving delta ribozyme derived from the antigenomic strand of the hepatitis delta virus. Initially, an RNA-DNA mixed ribozyme composed of 26 deoxyribo- (specifically the nucleotides forming the P2 stem and the P4 stem-loop) and 31 ribonucleotides (those forming the catalytic center) was engineered. This mixed ribozyme catalyzed the cleavage of a small substrate with kinetic parameters virtually identical to those of the all-RNA ribozyme. The further substitution of deoxyribose for ribose residues permitted us to investigate the contribution of all 2'-OHs to catalysis. Determination of the kinetic parameters for the cleavage reaction of the resulting ribozymes revealed (i) 10 2'-OH groups appear to be important in supporting the formation of several hydrogen bonds within the catalytic core, (ii) none of the important 2'-OHs seem to coordinate a magnesium cation, and (iii) 1 of the tested RNA-DNA mixed polymers appeared to stabilize the ribozyme-substrate transition-state complex, resulting in an improvement over the all-RNA counterpart. The contribution of the 2'-OHs to the catalytic mechanism is discussed, and differences with the crystal structure of a genomic delta self-cleaved product are explained. Clearly, the 2'-OHs are essential components of the network of interactions involved in the formation of the catalytic center of the delta ribozyme.

NMR analysis of bovine tRNATrp: conformation dependence of Mg2+ binding

NMR was used to study the solution structure of bovine tRNA(Trp) hyperexpressed in Escherichia coli. With the use of (15)N labeling and site-directed mutagenesis to assign overlapping resonances through the base pair replacement of U(71)A(2) by G(2)C(71), U(27)A(43) by G(27)C(43), and G(12)C(23) by U(12)A(23), the resonances of all 26 observable imino protons in the helical regions and in the tertiary interactions were assigned unambiguously by means of two-dimensional nuclear Overhauser effect spectroscopy and heteronuclear single quantum coherence methods. When the discriminator base A(73) and the G(12)C(23) base pair on the D stem, two identity elements on bovine tRNA(Trp) that are important for effective recognition by tryptophanyl-tRNA synthetase, were mutated to the ineffective forms of G(73) and U(12)A(23), respectively, NMR analysis revealed an important conformational change in the U(12)A(23) mutant but not in the G(73) mutant molecule. Thus A(73) appears to be directly recognized by tryptophanyl-tRNA synthetase, and G(12)C(23) represents an important structural determinant. Mg(2+) effects on the assigned resonances of imino protons allowed the identification of strong, medium, and weak Mg(2+) binding sites in tRNA(Trp). Strong Mg(2+) binding modes were associated with the residues G(7), s(4)U(8) (where s(4)U is 4-thiouridine), G(12), and U(52). The observations that G(42) was associated with strong Mg(2+) binding in only the U(12)A(23) mutant tRNA(Trp) but not the wild type or G(73) mutant tRNA(Trp) and that the G(7), s(4)U(8), G(24), and G(22) imino protons are associated with a two-site Mg(2+) binding mode in wild type and G(73) mutant but only a one-site mode in the U(12)A(23) mutant established the occurrence of conformational change in the U(12)A(23) mutant tRNA(Trp). These observations also established the dependence of Mg(2+) binding on tRNA conformation and the usefulness of Mg(2+) binding sites as conformational probes. The thermal titration of tRNA(Trp) in the presence and absence of 10 mm Mg(2+) indicated that overall tRNA(Trp) structure stability was increased by more than 15 degrees C by the presence of Mg(2+).

The role of metal ions in RNA catalysis

Understanding the catalytic mechanisms of RNA enzymes remains an important and intriguing challenge - one that has grown in importance since the recent demonstration that the ribosome is a ribozyme. At first, it seemed that all RNA enzymes compensate for the limited chemical versatility of ribonucleotide functional groups by recruiting obligatory metal ion cofactors to carry out catalytic chemistry. Mechanistic studies of the large self-splicing and pre-tRNA-processing ribozymes continue to support this idea, yielding increasingly detailed views of RNA active sites as scaffolds for positioning catalytic metal ions. Re-evaluation of the methodologies used to distinguish catalytic and structural roles for metal ions, however, has challenged this notion in the case of the small self-cleaving RNAs. Recent studies of the small ribozymes blur the distinction between catalytic and structural roles for metal ions, and suggest that RNA nucleobases have a previously unrecognized capacity for mediating catalytic chemistry.

59Co solid-state NMR as a new probe for elucidating metal binding in polynucleotides

Although magnesium fulfills several essential biochemical roles, direct studies on this ion are complicated by its unfavorable spectroscopic characteristics. This contribution explores the possibility of monitoring magnesium-nucleic acid binding via a combination of [Co(NH3)6]3+ as surrogate for [Mg(H2O)6]2+, and of high-resolution solid-state 59Co NMR as a spectroscopic probe. Such strategy quenches fast cationic exchanges between bound and free states, while exploiting the superior NMR properties of the 59Co spin. Experiments on relatively small amounts of tRNA can then discern resonances corresponding to different metal binding environments. These characterizations were assisted by studies on model compounds and by multinuclear 31P-59Co recoupling experiments.

NAIM and site-specific functional group modification analysis of RNase P RNA: magnesium dependent structure within the conserved P1-P4 multihelix junction contributes to catalysis

The tRNA processing endonuclease ribonuclease P contains an essential and highly conserved RNA molecule (RNase P RNA) that is the catalytic subunit of the enzyme. To identify and characterize functional groups involved in RNase P RNA catalysis, we applied self-cleaving ribozyme-substrate conjugates, on the basis of the RNase P RNA from Escherichia coli, in nucleotide analogue interference mapping (NAIM) and site-specific modification experiments. At high monovalent ion concentrations (3 M) that facilitate protein-independent substrate binding, we find that the ribozyme is largely insensitive to analogue substitution and that concentrations of Mg2+ (1.25 mM) well below that necessary for optimal catalytic rate (>100 mM) are required to produce interference effects because of modification of nucleotide bases. An examination of the pH dependence of the reaction rate at 1.25 mM Mg2+ indicates that the increased sensitivity to analogue interference is not due to a change in the rate-limiting step. The nucleotide positions detected by NAIM under these conditions are located exclusively in the catalytic domain, consistent with the proposed global structure of the ribozyme, and predominantly occur within the highly conserved P1-P4 multihelix junction. Several sensitive positions in J3/4 and J2/4 are proximal to a previously identified site of divalent metal ion binding in the P1-P4 element. Kinetic analysis of ribozymes with site-specific N7-deazaadenosine and deazaguanosine modifications in J3/4 was, in general, consistent with the interference results and also permitted the analysis of sites not accessible by NAIM. These results show that, in this region only, modification of the N7 positions of A62, A65, and A66 resulted in measurable effects on reaction rate and modification at each position displayed distinct sensitivities to Mg2+ concentration. These results reveal a restricted subset of individual functional groups within the catalytic domain that are particularly important for substrate cleavage and demonstrate a close association between catalytic function and metal ion-dependent structure in the highly conserved P1-P4 multihelix junction.

The contribution of 2'-hydroxyls to the cleavage activity of the Neurospora VS ribozyme

We have used nucleotide analog interference mapping and site-specific substitution to determine the effect of 2'-deoxynucleotide substitution of each nucleotide in the VS ribozyme on the self-cleavage reaction. A large number of 2'-hydroxyls (2'-OHs) that contribute to cleavage activity of the VS ribozyme were found distributed throughout the core of the ribozyme. The locations of these 2'-OHs in the context of a recently developed helical orientation model of the VS ribozyme suggest roles in multi-stem junction structure, helix packing, internal loop structure and catalysis. The functional importance of three separate 2'-OHs supports the proposal that three uridine turns contribute to local and long-range tertiary structure formation. A cluster of important 2'-OHs near the loop that is the candidate region for the active site and one very important 2'-OH in the loop that contains the cleavage site confirm the functional importance of these two loops. A cluster of important 2'-OHs lining the minor groove of stem-loop I and helix II suggests that these regions of the backbone may play an important role in positioning helices in the active structure of the ribozyme.

Structural analysis of metal ion ligation to nucleotides and nucleic acids using pulsed EPR spectroscopy

Metal ions play key structural and functional roles in many nucleic acid systems, particularly as required cofactors for many catalytic RNA molecules (ribozymes). We apply the pulsed EPR technologies of electron spin-echo envelope modulation and electron spin-echo-electron nuclear double resonance to the structural analysis of the paramagnetic metal ion Mn(II) bound to nucleotides and nucleic acids. We demonstrate that pulsed EPR, supplemented with specific isotope labeling, can characterize ligation to nucleotide base nitrogens, outer-sphere interactions with phosphate groups, distances to sites of specific (2)H atom labels, and the hydration level of the metal ion. These techniques allow a comprehensive structural analysis of the mononucleotide model system MnGMP. Spectra of phenylalanine-specific transfer RNA from budding yeast and of the hammerhead ribozyme demonstrate the applicability of the methods to larger, structured RNA systems. This suite of experiments opens the way to detailed structural characterization of specifically bound metal ions in a variety of ribozymes and other nucleic acids of biological interest.

Structural remodeling of an A + U-rich RNA element by cation or AUF1 binding

Association of AUF1 with A + U-rich elements (AREs) induces rapid cytoplasmic degradation of mRNAs containing these sequences, involving the recruitment or assembly of multisubunit trans-acting complexes on the mRNA. Recently, we reported that Mg(2+)-induced conformational changes in the ARE from tumor necrosis factor alpha mRNA inhibited AUF1 binding and oligomerization activities on this substrate (Wilson, G. M., Sutphen, K., Chuang, K., and Brewer, G. (2001) J. Biol. Chem. 276, 8695-8704). In this study, resonance energy transfer was employed to characterize structural changes in RNA substrates in response to cation- and AUF1-binding events. An RNA substrate containing the tumor necrosis factor alpha ARE displayed a weak conformational transition in the absence of added cations but was cooperatively stabilized by Mg(2+). Additional assays demonstrated a strong preference for small, multivalent cations, suggesting that the folded RNA structure was stabilized by counterion neutralization at discrete regions of high negative charge density. Association of AUF1 with cognate RNA substrates also induced formation of condensed RNA structures, although distinct from the folded structure stabilized by multivalent cations. Taken together, these experiments indicate that association of AUF1 with an ARE may function to remodel local RNA structures, which may be a prerequisite for subsequent recruitment of additional trans-acting factors.

Ribozyme structures and mechanisms

The past few years have seen exciting advances in understanding the structure and function of catalytic RNA. Crystal structures of several ribozymes have provided detailed insight into the folds of RNA molecules. Models of other biologically important RNAs have been constructed based on structural, phylogenetic, and biochemical data. However, many questions regarding the catalytic mechanisms of ribozymes remain. This review compares the structures and possible catalytic mechanisms of four small self-cleaving RNAs: the hammerhead, hairpin, hepatitis delta virus, and in vitro-selected lead-dependent ribozymes. The organization of these small catalysts is contrasted to that of larger ribozymes, such as the group I intron.

Monitoring the structure of Escherichia coli RNase P RNA in the presence of various divalent metal ions

Lead(II)-induced cleavage can be used as a tool to probe conformational changes in RNA. In this report, we have investigated the conformation of M1 RNA, the catalytic subunit of Escherichia coli RNase P, by studying the lead(II)-induced cleavage pattern in the presence of various divalent metal ions. Our data suggest that the overall conformation of M1 RNA is very similar in the presence of Mg(2+), Mn(2+), Ca(2+), Sr(2+) and Ba(2+), while it is changed compared to the Mg(2+)-induced conformation in the presence of other divalent metal ions, Cd(2+) for example. We also observed that correct folding of some M1 RNA domains is promoted by Pb(2+), while folding of other domain(s) requires the additional presence of other divalent metal ions, cobalt(III) hexamine or spermidine. Based on the suppression of Pb(2+) cleavage at increasing concentrations of various divalent metal ions, our findings suggest that different divalent metal ions bind with different affinities to M1 RNA as well as to an RNase P hairpin-loop substrate and yeast tRNA(Phe). We suggest that this approach can be used to obtain information about the relative binding strength for different divalent metal ions to RNA in general, as well as to specific RNA divalent metal ion binding sites. Of those studied in this report, Mn(2+) is generally among the strongest RNA binders.

Altering the intermediate in the equilibrium folding of unmodified yeast tRNAPhe with monovalent and divalent cations

The isothermal equilibrium folding of the unmodified yeast tRNA(Phe) is studied as a function of Na(+), Mg(2+), and urea concentration with hydroxyl radical protection, circular dichroism, and diethyl pyrocarbonate (DEPC) modification. These assays indicate that this tRNA folds in Na(+) alone. Similar to folding in Mg(2+), folding in Na(+) can be described by two transitions, unfolded-to-intermediate-to-native. The I-to-N transition has a Na(+) midpoint of approximately 0.5 M and a Hill constant of approximately 4. Unexpectedly, the urea m-value, the dependence of free energy on urea concentration, for the I-to-N transition is significantly smaller in Na(+) than in Mg(2+), 0.4 versus 1.7 kcal mol(-1) M(-1), indicating that more structure is formed in the Mg(2+)-induced transition. DEPC modification indicates that the I state in Na(+)-induced folding contains all four helices of tRNA and the I-to-N transition primarily corresponds to the formation of the tertiary structure. In contrast, the intermediate in Mg(2+)-induced folding contains only three helices, and the I-to-N transition corresponds to the formation of the acceptor stem plus tertiary structure. The cation dependence of the intermediates arises from the differences in the stability of the acceptor stem and the tertiary structure. The acceptor stem is stable at a lower Na(+) concentration than required for the tertiary structure formation. The relative stability is reversed in Mg(2+) so that the acceptor stem and the tertiary structure form simultaneously in the I-to-N transition. These results demonstrate that formation of the RNA secondary structure can be independent or coupled to the formation of the tertiary structure depending on their relative stability in monovalent and divalent ions.

Electrostatic mechanisms of DNA deformation

The genomes of higher cells consist of double-helical DNA, a densely charged polyelectrolyte of immense length. The intrinsic physical properties of DNA, as well as the properties of its complexes with proteins and ions, are therefore of fundamental interest in understanding the functions of DNA as an informational macromolecule. Because individual DNA molecules often exceed 1 cm in length, it is clear that DNA bending, folding, and interaction with nuclear proteins are necessary for packaging genomes in small volumes and for integrating the nucleotide sequence information that guides genetic readout. This review first focuses on recent experiments exploring how the shape of the densely charged DNA polymer and asymmetries in its surrounding counterion distribution mutually influence one another. Attention is then turned to experiments seeking to discover the degree to which asymmetric phosphate neutralization can lead to DNA bending in protein-DNA complexes. It is argued that electrostatic effects play crucial roles in the intrinsic, sequence-dependent shape of DNA and in DNA shapes induced by protein binding.

Ribozyme structures and mechanisms

The past few years have seen exciting advances in understanding the structure and function of catalytic RNA. Crystal structures of several ribozymes have provided detailed insight into the folds of RNA molecules. Models of other biologically important RNAs have been constructed based on structural, phylogenetic, and biochemical data. However, many questions regarding the catalytic mechanisms of ribozymes remain. This review compares the structures and possible catalytic mechanisms of four small self-cleaving RNAs: the hammerhead, hairpin, hepatitis delta virus, and in vitro-selected lead-dependent ribozymes. The organization of these small catalysts is contrasted to that of larger ribozymes, such as the group I intron.

The crystal structure of yeast phenylalanine tRNA at 2.0 A resolution: cleavage by Mg(2+) in 15-year old crystals

We have re-determined the crystal structure of yeast tRNA(Phe) to 2. 0 A resolution using 15 year old crystals. The accuracy of the new structure, due both to higher resolution data and formerly unavailable refinement methods, consolidates the previous structural information, but also reveals novel details. In particular, the water structure around the tightly bound Mg(2+) is now clearly resolved, and hence provides more accurate information on the geometry of the magnesium-binding sites and the role of water molecules in coordinating the metal ions to the tRNA. We have assigned a total of ten magnesium ions and identified a partly conserved geometry for high-affinity Mg(2+ )binding. In the electron density map there is also clear density for a spermine molecule binding in the major groove of the TPsiC arm and also contacting a symmetry-related tRNA molecule. Interestingly, we have also found that two specific regions of the tRNA in the crystals are partially cleaved. The sites of hydrolysis are within the D and anticodon loops in the vicinity of Mg(2+).