Identification and characterization of metal ion binding sites in RNA

On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules

RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.

Magnesium ions mediate ligand binding and conformational transition of the SAM/SAH riboswitch

The SAM/SAH riboswitch binds S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH) with similar affinities. Mg2+ is generally known to stabilize RNA structures by neutralizing phosphates, but how it contributes to ligand binding and conformational transition is understudied. Here, extensive molecular dynamics simulations (totaling 120 μs) predicted over 10 inner-shell Mg2+ ions in the SAM/SAH riboswitch. Six of them line the two sides of a groove to widen it and thereby pre-organize the riboswitch for ligand entry. They also form outer-shell coordination with the ligands and stabilize an RNA-ligand hydrogen bond, which effectively diminishes the selectivity between SAM and SAH. One Mg2+ ion unique to the apo form maintains the Shine-Dalgarno sequence in an autonomous mode and thereby facilitates its release for ribosome binding. Mg2+ thus plays vital roles in SAM/SAH riboswitch function.

Magnesium ions mediate ligand binding and conformational transition of the SAM/SAH riboswitch

The SAM/SAH riboswitch binds S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH) with similar affinities. Mg 2+ is generally known to stabilize RNA structures by neutralizing phosphates, but how it contributes to ligand binding and conformational transition is understudied. Here, extensive molecular dynamics simulations (totaling 120 μs) identified over 10 inner-shell Mg 2+ ions in the SAM/SAH riboswitch. Six of them line the two sides of a groove to widen it and thereby pre-organize the riboswitch for ligand entry. They also form outer-shell coordination with the ligands and stabilize an RNA-ligand hydrogen bond, which effectively diminish the selectivity between SAM and SAH. One Mg 2+ ion unique to the apo form maintains the Shine-Dalgarno sequence in an autonomous mode and thereby facilitates its release for ribosome binding. Mg 2+ thus plays vital roles in SAM/SAH riboswitch function.

Structural and Functional Annotation of Hypothetical Proteins from the Microsporidia Species Vittaforma corneae ATCC 50505 Using in silico Approaches

Microsporidia are spore-forming eukaryotes that are related to fungi but have unique traits that set them apart. They have compact genomes as a result of evolutionary gene loss associated with their complete dependency on hosts for survival. Despite having a relatively small number of genes, a disproportionately high percentage of the genes in microsporidia genomes code for proteins whose functions remain unknown (hypothetical proteins-HPs). Computational annotation of HPs has become a more efficient and cost-effective alternative to experimental investigation. This research developed a robust bioinformatics annotation pipeline of HPs from Vittaforma corneae, a clinically important microsporidian that causes ocular infections in immunocompromised individuals. Here, we describe various steps to retrieve sequences and homologs and to carry out physicochemical characterization, protein family classification, identification of motifs and domains, protein-protein interaction network analysis, and homology modelling using a variety of online resources. Classification of protein families produced consistent findings across platforms, demonstrating the accuracy of annotation utilizing in silico methods. A total of 162 out of 2034 HPs were fully annotated, with the bulk of them categorized as binding proteins, enzymes, or regulatory proteins. The protein functions of several HPs from Vittaforma corneae were accurately inferred. This improved our understanding of microsporidian HPs despite challenges related to the obligate nature of microsporidia, the absence of fully characterized genes, and the lack of homologous genes in other systems.

Petkov P, Milenov TI, Vayssilov GN. Sodium and Magnesium Ion Location at the Backbone and at the Nucleobase of RNA: Ab Initio Molecular Dynamics in Water Solution

The interactions between Na⁺ or Mg²⁺ ions with different parts of single-stranded RNA molecules, namely, the oxygen atoms from the phosphate groups or the guanine base, in water solution have been studied using first-principles molecular dynamics. Sodium ions were found to be much more mobile than Mg²⁺ ions and readily underwent transitions between a state directly bonded to RNA oxygen atoms and a completely solvated state. The inner solvation shell of Na⁺ ions fluctuated stochastically at a femtosecond timescale coordinating on average 5 oxygen atoms for bonded Na⁺ ions and 5.5 oxygen atoms for solvated Na⁺ ions. In contrast, the inner solvation shell of Mg²⁺ ions was stable in both RNA-bonded and completely solvated states. In both cases, Mg²⁺ ions coordinated 6 oxygen atoms from the inner solvation shell. Consistent with their stable solvation shells, Mg²⁺ ions were more effective than Na⁺ ions in stabilizing the RNA backbone conformation. The exclusion zones between the first and second solvation shells, solvation shell widths, and angles for binding to carbonyl oxygen of guanine for solvated Na⁺ or Mg²⁺ ions exhibited a number of quantitative differences when compared with RNA crystallographic data. The presented results support the distinct capacity of Mg²⁺ ions to support the RNA structure not only in the crystal phase but also in the dynamic water environment both on the side of the phosphate moiety and on the side of the nucleobase.

Structural Biology of the Enterovirus Replication-Linked 5'-Cloverleaf RNA and Associated Virus Proteins

Although enteroviruses are associated with a wide variety of diseases and conditions, their mode of replication is well conserved. Their genome is carried as a single, positive-sense RNA strand. At the 5' end of the strand is an approximately 90-nucleotide self-complementary region called the 5' cloverleaf, or the oriL. This noncoding region serves as a platform upon which host and virus proteins, including the 3B, 3C, and 3D virus proteins, assemble in order to initiate replication of a negative-sense RNA strand. The negative strand in turn serves as a template for synthesis of multiple positive-sense RNA strands. Building on structural studies of individual RNA stem-loops, the structure of the intact 5' cloverleaf from rhinovirus has recently been determined via nuclear magnetic resonance/small-angle X-ray scattering (NMR/SAXS)-based methods, while structures have also been determined for enterovirus 3A, 3B, 3C, and 3D proteins. Analysis of these structures, together with structural and modeling studies of interactions between host and virus proteins and RNA, has begun to provide insight into the enterovirus replication mechanism and the potential to inhibit replication by blocking these interactions.

Conformational flexibility in the enterovirus RNA replication platform

A presumed RNA cloverleaf (5'CL), located at the 5'-most end of the noncoding region of the enterovirus genome, is the primary established site for initiation of genomic replication. Stem-loop B (SLB) and stem-loop D (SLD), the two largest stem-loops within the 5'CL, serve as recognition sites for protein interactions that are essential for replication. Here we present the solution structure of rhinovirus serotype 14 5'CL using a combination of nuclear magnetic resonance spectroscopy and small-angle X-ray scattering. In the absence of magnesium, the structure adopts an open, somewhat extended conformation. In the presence of magnesium, the structure compacts, bringing SLB and SLD into close contact, a geometry that creates an extensive accessible major groove surface, and permits interaction between the proteins that target each stem-loop.

Dynamics of metal ions around an RNA molecule

Noncoding RNA molecules take part in many biological processes, while metal ions play crucial roles in helping RNAs to perform their functions. However, the statics and dynamics of these metal ions around RNA molecules are still not well understood. In this work, we report a detailed molecular dynamics study of the type-I preQ_{1}-bound riboswitch aptamer domain (PRAD) at different ionic conditions (K^{+}, Na^{+}, and Mg^{2+}). The results show that the structural properties and flexibility of the PRAD molecule greatly influence the distributions and dynamics of metal ions around it. Simultaneously, Na^{+} ions show a stronger competitiveness with Mg^{2+} ions than K^{+} ions, and the three types of metal ions have different modes of interaction with the RNA molecule. Furthermore, we have also investigated specific binding sites of metal ions on the PRAD molecule and found that the dynamics and hydration structures of metal ions located at the ion-binding sites were obviously affected by the RNA structure near these ion-binding sites. These results may be useful to understand the role of the metal ions in noncoding RNA functions.

Influence of Na+ and Mg2+ ions on RNA structures studied with molecular dynamics simulations

The structure of ribonucleic acid (RNA) polymers is strongly dependent on the presence of, in particular Mg2+ cations to stabilize structural features. Only in high-resolution X-ray crystallography structures can ions be identified reliably. Here, we perform molecular dynamics simulations of 24 RNA structures with varying ion concentrations. Twelve of the structures were helical and the others complex folded. The aim of the study is to predict ion positions but also to evaluate the impact of different types of ions (Na+ or Mg2+) and the ionic strength on structural stability and variations of RNA. As a general conclusion Mg2+ is found to conserve the experimental structure better than Na+ and, where experimental ion positions are available, they can be reproduced with reasonable accuracy. If a large surplus of ions is present the added electrostatic screening makes prediction of binding-sites less reproducible. Distinct differences in ion-binding between helical and complex folded structures are found. The strength of binding (ΔG‡ for breaking RNA atom-ion interactions) is found to differ between roughly 10 and 26 kJ/mol for the different RNA atoms. Differences in stability between helical and complex folded structures and of the influence of metal ions on either are discussed.

The X-ray Structures of Six Octameric RNA Duplexes in the Presence of Different Di- and Trivalent Cations

Due to the polyanionic nature of RNA, the principles of charge neutralization and electrostatic condensation require that cations help to overcome the repulsive forces in order for RNA to adopt a three-dimensional structure. A precise structural knowledge of RNA-metal ion interactions is crucial to understand the mechanism of metal ions in the catalytic or regulatory activity of RNA. We solved the crystal structure of an octameric RNA duplex in the presence of the di- and trivalent metal ions Ca(2+), Mn(2+), Co(2+), Cu(2+), Sr(2+), and Tb(3+). The detailed investigation reveals a unique innersphere interaction to uracil and extends the knowledge of the influence of metal ions for conformational changes in RNA structure. Furthermore, we could demonstrate that an accurate localization of the metal ions in the X-ray structures require the consideration of several crystallographic and geometrical parameters as well as the anomalous difference map.

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.

Metal ion-RNA interactions studied via multinuclear NMR

Metal ions are indispensable for ribonucleic acids (RNAs) folding and activity. First they act as charge neutralization agents, allowing the RNA molecule to attain the complex active three dimensional structure. Second, metal ions are eventually directly involved in function. Nuclear magnetic resonance (NMR) spectroscopy offers several ways to study the RNA-metal ion interactions at an atomic level. Here, we first focus on special requirements for NMR sample preparation for this kind of experiments: the practical aspects of in vitro transcription and purification of small (<50 nt) RNA fragments are described, as well as the precautions that must be taken into account when a sample for metal ion titration experiments is prepared. Subsequently, we discuss the NMR techniques to accurately locate and characterize metal ion binding sites in a large RNA. For example, (2) J-[(1)H,(15)N]-HSQC (heteronuclear single quantum coherence) experiments are described to qualitatively distinguish between different modes of interaction. Finally, part of the last section is devoted to data analysis; this is how to calculate intrinsic affinity constants.

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.

Modulation of the stability of the Salmonella fourU-type RNA thermometer

RNA thermometers are translational control elements that regulate the expression of bacterial heat shock and virulence genes. They fold into complex secondary structures that block translation at low temperatures. A temperature increase releases the ribosome binding site and thus permits translation initiation. In fourU-type RNA thermometers, the AGGA sequence of the SD region is paired with four consecutive uridines. We investigated the melting points of the wild-type and mutant sequences. It was decreased by 5°C when a stabilizing GC basepair was exchanged by an AU pair or increased by 11°C when an internal AG mismatch was converted to a GC pair, respectively. Stabilized or destabilized RNA structures are directly correlated with decreased or increased in vivo gene expression, respectively. Mg(2+) also affected the melting point of the fourU thermometer. Variations of the Mg(2+) concentration in the physiological range between 1 and 2 mM translated into a 2.8°C shift of the melting point. Thus, Mg(2+) binding to the hairpin RNA is regulatory relevant. Applying three different NMR techniques, two Mg(2+) binding sites were found in the hairpin structure. One of these binding sites could be identified as outer sphere binding site that is located within the fourU motif. Binding of the two Mg(2+) ions exhibits a positive cooperativity with a Hill coefficient of 1.47. Free energy values ΔG for Mg(2+) binding determined by NMR are in agreement with data determined from CD measurements.

Specific Interactions of Divalent Metal Ions with a DNA Duplex Containing the d(CA)n/(GT)nTandem Repeat

Divalent metal ions are essential for maintaining functional states of the DNA molecule. Their participation in DNA structure is modulated by the base sequence and varies depending on the nature of the ion. The present investigation addresses the interaction of Ca2+ ions with a tandem repeat of two CA dinucleotides, (CA)2/(TG)2. The binding of Ca2+ to the repeat is monitored by nuclear magnetic resonance (NMR) spectroscopy using chemical shift mapping. Parallel experiments monitor binding of Mg2+ ions to the repeat as well as binding of each ion to a DNA duplex in which the (CA)2/(TG)2 repeat is eliminated. The results reveal that the direction and the magnitude of chemical shift changes induced by Ca2+ ions in the NMR spectra of the repeat are different from those induced by Mg2+ ions. The differences between the two cations are significantly diminished by the elimination of the (CA)2/(TG)2 repeat. These findings suggest a specific interaction of Ca2+ ions with the (CA)2/(TG)2 motif. The specificity of the interaction resides in the two A-T base pairs of the repeat, and it involves the major groove of the first A-T base pair and both grooves of the second A-T base pair.

Dynamics and metal ion binding in the U6 RNA intramolecular stem-loop as analyzed by NMR

The U6 RNA intramolecular stem-loop (ISL) is a conserved component of the spliceosome, and contains an essential metal ion binding site centered between a protonated adenine, A79, and U80. Correlated with protonation of A79, U80 undergoes a base-flipping conformational change accompanied by significant helical movement. We have investigated the dynamics of the U6 ISL by analyzing the power dependence of 13C NMR relaxation rates in the rotating frame. The data provide evidence that the conformational transition is centered around an exchange lifetime of 84 micros. The U80 nucleotide displays low internal mobility on the picosecond time-scale at pH 7.0 but high internal mobility at pH 6.0, in agreement with the global transition resulting in the base of U80 adopting a looped-out conformation with increased dynamic disorder. A kinetic analysis suggests that the conformational change, rather than adenine protonation, is the rate-limiting step in the pathway of the conformational transition. Two nucleotides, U70 and U80, were found from chemical shift perturbation mapping to interact with the magnesium ion, with apparent K(d) values in the micromolar to millimolar range. These nucleotides also displayed metal ion-induced elevation of R1 rates, which can be explained by a model that assumes dynamic metal ion coordination concomitant with an induced higher shielding anisotropy for the base 13C nuclei. Addition of Mg2+ shifts the conformational equilibrium toward the high-pH (base-stacked) structure, accompanied by a significant drop in the apparent pK(a) of A79.

Ribosomal RNA in Alzheimer disease is oxidized by bound redox-active iron

Oxidative modification of cytoplasmic RNA in vulnerable neurons is an important, well documented feature of the pathophysiology of Alzheimer disease. Here we report that RNA-bound iron plays a pivotal role for RNA oxidation in vulnerable neurons in Alzheimer disease brain. The cytoplasm of hippocampal neurons showed significantly higher redox activity and iron(II) staining than age-matched controls. Notably, both were susceptible to RNase, suggesting a physical association of iron(II) with RNA. Ultrastructural analysis further suggested an endoplasmic reticulum association. Both rRNA and mRNA showed twice the iron binding as tRNA. rRNA, extremely abundant in neurons, was considered to provide the greatest number of iron binding sites among cytoplasmic RNA species. Interestingly, the difference of iron binding capacity disappeared after denaturation of RNA, suggesting that the higher order structure may contribute to the greater iron binding of rRNA. Reflecting the difference of iron binding capacity, oxidation of rRNA by the Fenton reaction formed 13 times more 8-hydroxyguanosine than tRNA. Consistent with in situ findings, ribosomes purified from Alzheimer hippocampus contained significantly higher levels of RNase-sensitive iron(II) and redox activity than control. Furthermore, only Alzheimer rRNA contains 8-hydroxyguanosine in reverse transcriptase-PCR. Addressing the biological significance of ribosome oxidation by redox-active iron, in vitro translation with oxidized ribosomes from rabbit reticulocyte showed a significant reduction of protein synthesis. In conclusion these results suggest that rRNA provides a binding site for redox-active iron and serves as a redox center within the cytoplasm of vulnerable neurons in Alzheimer disease in advance of the appearance of morphological change indicating neurodegeneration.

Modulation of the distance dependence of paramagnetic relaxation enhancements by CSA x DSA cross-correlation

Paramagnetic metal ions with fast-relaxing electronic spin and anisotropic susceptibility tensor provide a rich source of structural information that can be derived from pseudo-contact shifts, residual dipolar couplings, dipole-dipole Curie spin cross-correlation, and paramagnetic relaxation enhancements. The present study draws attention to a cross-correlation effect between nuclear relaxation due to anisotropic chemical shielding (CSA) and due to the anisotropic dipolar shielding (DSA) caused by the electronic Curie spin. This CSA x DSA cross-correlation contribution seems to have been overlooked in previous interpretations of paramagnetic relaxation enhancements. It is shown to be sufficiently large to compromise the 1/r6 distance dependence usually assumed. The effect cannot experimentally be separated from auto-correlated DSA relaxation. It can increase or decrease the observed paramagnetic relaxation enhancement. Under certain conditions, the effect can dominate the entire paramagnetic relaxation, resulting in nuclear resonances narrower than in the absence of the paramagnetic center. CSAxDSA cross-correlation becomes important when paramagnetic relaxation is predominantly due to the Curie rather than the Solomon mechanism. Therefore the effect is most pronounced for relaxation by metal ions with large magnetic susceptibility and fast-relaxing electron spin. It most strongly affects paramagnetic enhancements of transverse relaxation in macromolecules and of longitudinal relaxation in small molecules.

NMR studies of the structure and Mg2+ binding properties of a conserved RNA motif of EMCV picornavirus IRES element

The structure and Mg(2+) binding properties of a conserved 75mer RNA motif of the internal ribosome entry site (IRES) element of encephalomyocarditis virus picornavirus have been investigated by (1)H-NMR and UV melting experiments. The assignment of the imino proton resonances with characteristic chemical shift dispersion for canonical and non-canonical base pairs confirmed the predicted secondary structure of the 75mer and its fragments. Addition of Mg(2+) resulted in a dramatic increase in apparent melting temperature, with the 75mer RNA registering the biggest increase, from 63 to 80 degrees C, thus providing evidence for enhanced stability arising from Mg(2+) binding. Similarly, addition of Mg(2+) induced selective changes to the chemical shifts of the imino protons of a GCGA tetraloop in the 75mer, that is essential for IRES activity, thereby highlighting a possible structural role for Mg(2+) in the folding of the 75mer. Significantly, the same protons show retarded exchange to water solvent, even at elevated temperature, which suggest that Mg(2+) induces a conformational rearrangement of the 75mer. Thus, we propose that Mg(2+) serves two important roles: (i) enhancing thermodynamic stability of the 75mer RNA (and its submotifs) via non-specific interactions with the phosphate backbone and (ii) promoting the folding of the 75mer RNA by binding to the GCGA tetraloop.

Structure of HCV IRES domain II determined by NMR

Complex RNA structures regulate many biological processes, but are often too large for structure determination by NMR methods. The 5' untranslated region (5' UTR) of the hepatitis C viral (HCV) RNA genome contains an internal ribosome entry site (IRES) that binds to 40S ribosomal subunits with high affinity and specificity to control translation. Domain II of the HCV IRES forms a 25-kDa folded subdomain that may alter ribosome conformation. We report here the structure of domain II as determined using an NMR approach that combines short- and long-range structural data. Domain II adopts a distorted L-shape structure, and its overall shape in the free form is markedly similar to its 40S subunit-bound form; this suggests how domain II may modulate 40S subunit conformation. The results show how NMR can be used for structural analysis of large biological RNAs.

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.