Metal ion binding to catalytic RNA molecules

Mg2+ Sensing by an RNA Fragment: Role of Mg2+-Coordinated Water Molecules

RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg²⁺, Ca²⁺, and K⁺) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg²⁺ solution. We show that the rigid network of Mg²⁺-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K⁺ and Ca²⁺ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.

The structural ensemble of a Holliday junction determined by X-ray scattering interference

The DNA four-way (Holliday) junction is the central intermediate of genetic recombination, yet key aspects of its conformational and thermodynamic properties remain unclear. While multiple experimental approaches have been used to characterize the canonical X-shape conformers under specific ionic conditions, the complete conformational ensemble of this motif, especially at low ionic conditions, remains largely undetermined. In line with previous studies, our single-molecule Förster resonance energy transfer (smFRET) measurements of junction dynamics revealed transitions between two states under high salt conditions, but smFRET could not determine whether there are fast and unresolvable transitions between distinct conformations or a broad ensemble of related states under low and intermediate salt conditions. We therefore used an emerging technique, X-ray scattering interferometry (XSI), to directly probe the conformational ensemble of the Holliday junction across a wide range of ionic conditions. Our results demonstrated that the four-way junction adopts an out-of-plane geometry under low ionic conditions and revealed a conformational state at intermediate ionic conditions previously undetected by other methods. Our results provide critical information to build toward a full description of the conformational landscape of the Holliday junction and underscore the utility of XSI for probing conformational ensembles under a wide range of solution conditions.

Molecular simulations of the pistol ribozyme: unifying the interpretation of experimental data and establishing functional links with the hammerhead ribozyme

The pistol ribozyme (Psr) is among the most recently discovered RNA enzymes and has been the subject of experiments aimed at elucidating the mechanism. Recent biochemical studies have revealed exciting clues about catalytic interactions in the active site not apparent from available crystallographic data. The present work unifies the interpretation of the existing body of structural and functional data on Psr by providing a dynamical model for the catalytically active state in solution from molecular simulation. Our results suggest that a catalytic Mg2+ ion makes inner-sphere contact with G33:N7 and outer-sphere coordination to the pro-RP of the scissile phosphate, promoting electrostatic stabilization of the dianionic transition state and neutralization of the developing charge of the leaving group through a metal-coordinated water molecule that is made more acidic by a hydrogen bond donated from the 2'OH of P32. This model is consistent with experimental activity-pH and mutagenesis data, including sensitivity to G33(7cG) and phosphorothioate substitution/metal ion rescue. The model suggests several experimentally testable predictions, including the response of cleavage activity to mutations at G42 and P32 positions in the ribozyme, and thio substitutions of the substrate in the presence of different divalent metal ions. Further, the model identifies striking similarities of Psr to the hammerhead ribozyme (HHr), including similar global fold, organization of secondary structure around an active site three-way junction, catalytic metal ion binding mode, and guanine general base. However, the specific binding mode and role of the Mg2+ ion, as well as a conserved 2'-OH in the active site, are interrelated but subtly different between the ribozymes.

Effects of Preferential Counterion Interactions on the Specificity of RNA Folding

The real-time search for native RNA structure is essential for the operation of regulatory RNAs. We previously reported that a fraction of the Azoarcus ribozyme achieves a compact structure in less than a millisecond. To scrutinize the forces that drive initial folding steps, we used time-resolved SAXS to compare the folding dynamics of this ribozyme in thermodynamically isostable concentrations of different counterions. The results show that the size of the fast-folding population increases with the number of available counterions and correlates with the flexibility of initial RNA structures. Within 1 ms of folding, Mg²⁺ exhibits a smaller preferential interaction coefficient per charge, ΔΓ₊/ Z, than Na⁺ or [Co(NH₃)₆]³⁺. The lower ΔΓ₊/ Z corresponds to a smaller yield of folded RNA, although Mg²⁺ stabilizes native RNA more efficiently than other ions at equilibrium. These results suggest that strong Mg²⁺-RNA interactions impede the search for globally native structure during early folding stages.

Evidence of a General Acid-Base Catalysis Mechanism in the 8-17 DNAzyme

DNAzymes are catalytic DNA molecules that can perform a variety of reactions. Although advances have been made in obtaining DNAzymes via in vitro selection and many of them have been developed into sensors and imaging agents for metal ions, bacteria, and other molecules, the structural features responsible for these enzymatic reactions are still not well understood. Previous studies of the 8-17 DNAzyme have suggested conserved guanines close to the phosphodiester transfer site may play a role in the catalytic reaction. To identify the specific guanine and functional group of the guanine responsible for the reaction, we herein report the effects of replacing G1.1 and G14 (G; p Ka,N1 = 9.4) with analogues with a different p Ka at the N1 position, such as inosine (G14I; p Ka,N1 = 8.7), 2,6-diaminopurine (G14diAP; p Ka,N1 = 5.6), and 2-aminopurine (G14AP; p Ka,N1 = 3.8) on pH-dependent reaction rates. A comparison of the pH dependence of the reaction rates of these DNAzymes demonstrated that G14 in the bulge loop next to the cleavage site, is involved in proton transfer at the catalytic site. In contrast, we did not find any evidence of G1.1 being involved in acid-base catalysis. These results support general acid-base catalysis as a feasible strategy used in DNA catalysis, as in RNA and protein enzymes.

Inner-Sphere Coordination of Divalent Metal Ion with Nucleobase in Catalytic RNA

Identification of the function of metal ions and the RNA moieties, particularly nucleobases, that bind metal ions is important in RNA catalysis. Here we combine single-atom and abasic substitutions to probe functions of conserved nucleobases in ribonuclease P (RNase P). Structural and biophysical studies of bacterial RNase P propose direct coordination of metal ions by the nucleobases of conserved uridine and guanosine in helix P4 of the RNA subunit (P RNA). To biochemically probe the function of metal ion interactions, we substituted the universally conserved bulged uridine (U51) in the P4 helix of circularly permuted Bacillus subtilis P RNA with 4-thiouridine, 4-deoxyuridine, and abasic modifications and G378/379 with 2-aminopurine, N7-deazaguanosine, and 6-thioguanosine. The functional group modifications of U51 decrease RNase P-catalyzed phosphodiester bond cleavage 16- to 23-fold, as measured by the single-turnover cleavage rate constant. The activity of the 4-thiouridine RNase P is partially rescued by addition of Cd(II) or Mn(II) ions. This is the first time a metal-rescue experiment provides evidence for inner-sphere divalent metal ion coordination with a nucleobase. Modifications of G379 modestly decrease the cleavage activity of RNase P, suggesting outer-sphere coordination of O6 on G379 to a metal ion. These data provide biochemical evidence for catalytically important interactions of the P4 helix of P RNA with metal ions, demonstrating that the bulged uridine coordinates at least one catalytic metal ion through an inner-sphere interaction. The combination of single-atom and abasic nucleotide substitutions provides a powerful strategy to probe functions of conserved nucleobases in large RNAs.

Crystal structure of a novel RNA motif that allows for precise positioning of a single metal ion

We have determined a crystal structure of an RNA duplex containing a novel metal-binding motif. The motif is composed of two sheared G○A base pairs, two unpaired A residues and four phosphate groups in close proximity. Four A residues make an A-A-A-A stacking column at the minor groove side and two G bases are highly inclined, thereby forming the pocket-shaped motif at the major groove side. In the present structure, a hydrated Sr²⁺ ion exists in the pocket and binds to the O6 and N7 atoms of the two G bases and four phosphate groups. According to the previously-reported metal-binding properties to RNA molecules, many of divalent cations, such as Mg²⁺, Mn²⁺, Co²⁺, Zn²⁺, Ba²⁺, Pb²⁺ and Cd²⁺, may bind to the motif. This metal-binding motif can be used as a modular building block that allows for precise positioning of a single metal ion in functional nucleic acid molecules.

Structural Basis of Catalysis in the Bacterial Monoterpene Synthases Linalool Synthase and 1,8-Cineole Synthase

Terpenoids form the largest and stereochemically most diverse class of natural products, and there is considerable interest in producing these by biocatalysis with whole cells or purified enzymes, and by metabolic engineering. The monoterpenes are an important class of terpenes and are industrially important as flavors and fragrances. We report here structures for the recently discovered Streptomyces clavuligerus monoterpene synthases linalool synthase (bLinS) and 1,8-cineole synthase (bCinS), and we show that these are active biocatalysts for monoterpene production using biocatalysis and metabolic engineering platforms. In metabolically engineered monoterpene-producing E. coli strains, use of bLinS leads to 300-fold higher linalool production compared with the corresponding plant monoterpene synthase. With bCinS, 1,8-cineole is produced with 96% purity compared to 67% from plant species. Structures of bLinS and bCinS, and their complexes with fluorinated substrate analogues, show that these bacterial monoterpene synthases are similar to previously characterized sesquiterpene synthases. Molecular dynamics simulations suggest that these monoterpene synthases do not undergo large-scale conformational changes during the reaction cycle, making them attractive targets for structured-based protein engineering to expand the catalytic scope of these enzymes toward alternative monoterpene scaffolds. Comparison of the bLinS and bCinS structures indicates how their active sites steer reactive carbocation intermediates to the desired acyclic linalool (bLinS) or bicyclic 1,8-cineole (bCinS) products. The work reported here provides the analysis of structures for this important class of monoterpene synthase. This should now guide exploitation of the bacterial enzymes as gateway biocatalysts for the production of other monoterpenes and monoterpenoids.

Metals induce transient folding and activation of the twister ribozyme

Twister is a small ribozyme present in almost all kingdoms of life that rapidly self-cleaves in variety of divalent metal ions. We used activity assays, bulk FRET and single-molecule FRET (smFRET) to understand how different metal ions promote folding and self-cleavage of the Oryza sativa Twister ribozyme. Although most ribozymes require additional Mg²⁺ for catalysis, Twister inverts this expectation, requiring 20–30 times less Mg²⁺ to self-cleave than to fold. Transition metals such as Co²⁺, Ni²⁺ and Zn²⁺ activate Twister more efficiently than Mg²⁺ ions. Although Twister is fully active in ≤ 0.5 mM MgCl₂, smFRET experiments showed that the ribozyme visits the folded state infrequently under these conditions. Comparison of folding and self-cleavage rates indicates that most folding events lead to catalysis, which correlates with metal bond strength. Thus, the robust activity of Twister reports on transient metal ion binding under physiological conditions.

The GlcN6P cofactor plays multiple catalytic roles in the glmS ribozyme

RNA enzymes (ribozymes) have remarkably diverse biological roles despite having limited chemical diversity. Protein enzymes enhance their reactivity through recruitment of cofactors; likewise, the naturally occurring glmS ribozyme uses the glucosamine-6-phosphate (GlcN6P) organic cofactor for phosphodiester bond cleavage. Prior structural and biochemical studies have implicated GlcN6P as the general acid. Here we describe new catalytic roles of GlcN6P through experiments and calculations. Large stereospecific normal thio effects and a lack of metal-ion rescue in the holoribozyme indicate that nucleobases and the cofactor play direct chemical roles and align the active site for self-cleavage. Large stereospecific inverse thio effects in the aporibozyme suggest that the GlcN6P cofactor disrupts an inhibitory interaction of the nucleophile. Strong metal-ion rescue in the aporibozyme reveals that this cofactor also provides electrostatic stabilization. Ribozyme organic cofactors thus perform myriad catalytic roles, thereby allowing RNA to compensate for its limited functional diversity.

Mediating physicochemical properties and paclitaxel release of pH-responsive H-type multiblock copolymer self-assembly nanomicelles through epoxidation

We focused on modulation of the physicochemical and biomedical properties of copolymer nanomicellesviaepoxidation, which provided significant improvements.

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.

Structure, dynamics, and mechanism of the lead-dependent ribozyme

Leadzyme is a small catalytic RNA that was identified by in vitro selection for Pb2+-dependent cleavage from a tRNA library. Leadzyme employs a unique two-step Pb2+-specific mechanism to cleave within its active site. NMR and crystal structures of the active site revealed different folding patterns, but neither features the in-line alignment for attack by the 2'-OH nucleophilic group. These experimentally determined structures most likely represent ground states and are catalytically inactive. There are significant dynamics of the active site and the motif samples multiple conformations at the ground states. Various metal ion binding sites have been identified, including one that may be occupied by a catalytic Pb2+. Based on functional group analysis, a computational model of the transition state has been proposed. This model features a unique base triple that is consistent with sequence and functional group requirements for catalysis. This structure is likely only populated transiently, but imposing appropriate conformational constraints may significantly stabilize this state thereby promoting catalysis. Other ions may inhibit the cleavage by competing for the Pb2+ binding site, or by stabilizing the ground state thereby suppressing its transition to the catalytically active conformation. Some rare earth ions can enhance the reaction via an unknown mechanism. Because of its unique chemistry and dynamic behavior, leadzyme can continue to serve as an excellent model system for teaching us RNA biology and chemistry.

Comparison of structural, thermodynamic, kinetic and mass transport properties of Mg(2+) ion models commonly used in biomolecular simulations

The prevalence of Mg(2+) ions in biology and their essential role in nucleic acid structure and function has motivated the development of various Mg(2+) ion models for use in molecular simulations. Currently, the most widely used models in biomolecular simulations represent a nonbonded metal ion as an ion-centered point charge surrounded by a nonelectrostatic pairwise potential that takes into account dispersion interactions and exchange effects that give rise to the ion's excluded volume. One strategy toward developing improved models for biomolecular simulations is to first identify a Mg(2+) model that is consistent with the simulation force fields that closely reproduces a range of properties in aqueous solution, and then, in a second step, balance the ion-water and ion-solute interactions by tuning parameters in a pairwise fashion where necessary. The present work addresses the first step in which we compare 17 different nonbonded single-site Mg(2+) ion models with respect to their ability to simultaneously reproduce structural, thermodynamic, kinetic and mass transport properties in aqueous solution. None of the models based on a 12-6 nonelectrostatic nonbonded potential was able to reproduce the experimental radial distribution function, solvation free energy, exchange barrier and diffusion constant. The models based on a 12-6-4 potential offered improvement, and one model in particular, in conjunction with the SPC/E water model, performed exceptionally well for all properties. The results reported here establish useful benchmark calculations for Mg(2+) ion models that provide insight into the origin of the behavior in aqueous solution, and may aid in the development of next-generation models that target specific binding sites in biomolecules.

Effect of Zn2+ binding and enzyme active site on the transition state for RNA 2'-O-transphosphorylation interpreted through kinetic isotope effects

Divalent metal ions, due to their ability to stabilize high concentrations of negative charge, are important for RNA folding and catalysis. Detailed models derived from the structures and kinetics of enzymes and from computational simulations have been developed. However, in most cases the specific catalytic modes involving metal ions and their mechanistic roles and effects on transition state structures remain controversial. Valuable information about the nature of the transition state is provided by measurement of kinetic isotope effects (KIEs). However, KIEs reflect changes in all bond vibrational modes that differ between the ground state and transition state. QM calculations are therefore essential for developing structural models of the transition state and evaluating mechanistic alternatives. Herein, we present computational models for Zn2+ binding to RNA 2'O-transphosphorylation reaction models that aid in the interpretation of KIE experiments. Different Zn2+ binding modes produce distinct KIE signatures, and one binding mode involving two zinc ions is in close agreement with KIEs measured for non-enzymatic catalysis by Zn2+ aquo ions alone. Interestingly, the KIE signatures in this specific model are also very close to those in RNase A catalysis. These results allow a quantitative connection to be made between experimental KIE measurements and transition state structure and bonding, and provide insight into RNA 2'O-ransphosphorylation reactions catalyzed by metal ions and enzymes. This article is part of a Special Issue entitled: Enzyme Transition States from Theory and Experiment.

Distance measurements between manganese(ii) and nitroxide spin-labels by DEER determine a binding site of Mn2+ in the HP92 loop of ribosomal RNA

W-band (95 GHz) double electron–electron resonance (DEER) distance measurements between Mn²⁺ and nitroxide spin labels were used to determine the location of a Mn²⁺ binding site within an RNA molecule.

Principles of ion recognition in RNA: insights from the group II intron structures

Metal ions promote both RNA folding and catalysis, thus being essential in stabilizing the structure and determining the function of large RNA molecules, including group II introns. The latter are self-splicing metalloribozymes, containing a heteronuclear four-metal-ion center within the active site. In addition to these catalytic ions, group II introns bind many other structural ions, including delocalized ions that bind the RNA diffusively and well-ordered ions that bind the RNA tightly with high occupancy. The latter ions, which can be studied by biophysical methods, have not yet been analyzed systematically. Here, we compare crystal structures of the group IIC intron from Oceanobacillus iheyensis and classify numerous site-bound ions, which are primarily localized in the intron core and near long-range tertiary contacts. Certain ion-binding sites resemble motifs observed in known RNA structures, while others are idiosyncratic to the group II intron. Particularly interesting are (1) ions proximal to the active site, which may participate in splicing together with the catalytic four-metal-ion center, (2) organic ions that bind regions predicted to interact with intron-encoded proteins, and (3) unusual monovalent ions bound to GU wobble pairs, GA mismatches, the S-turn, the tetraloop-receptor, and the T-loop. Our analysis extends the general principles by which ions participate in RNA structural organization and it will aid in the determination and interpretation of future RNA structures.

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.

Quantum mechanical/molecular mechanical free energy simulations of the self-cleavage reaction in the hepatitis delta virus ribozyme

The hepatitis delta virus (HDV) ribozyme catalyzes a self-cleavage reaction using a combination of nucleobase and metal ion catalysis. Both divalent and monovalent ions can catalyze this reaction, although the rate is slower with monovalent ions alone. Herein, we use quantum mechanical/molecular mechanical (QM/MM) free energy simulations to investigate the mechanism of this ribozyme and to elucidate the roles of the catalytic metal ion. With Mg²⁺ at the catalytic site, the self-cleavage mechanism is observed to be concerted with a phosphorane-like transition state and a free energy barrier of ∼13 kcal/mol, consistent with free energy barrier values extrapolated from experimental studies. With Na⁺ at the catalytic site, the mechanism is observed to be sequential, passing through a phosphorane intermediate, with free energy barriers of 2–4 kcal/mol for both steps; moreover, proton transfer from the exocyclic amine of protonated C75 to the nonbridging oxygen of the scissile phosphate occurs to stabilize the phosphorane intermediate in the sequential mechanism. To explain the slower rate observed experimentally with monovalent ions, we hypothesize that the activation of the O2′ nucleophile by deprotonation and orientation is less favorable with Na⁺ ions than with Mg²⁺ ions. To explore this hypothesis, we experimentally measure the pK_a of O2′ by kinetic and NMR methods and find it to be lower in the presence of divalent ions rather than only monovalent ions. The combined theoretical and experimental results indicate that the catalytic Mg²⁺ ion may play three key roles: assisting in the activation of the O2′ nucleophile, acidifying the general acid C75, and stabilizing the nonbridging oxygen to prevent proton transfer to it.

Experimental approaches for measuring pKa's in RNA and DNA

RNA and DNA carry out diverse functions in biology including catalysis, splicing, gene regulation, and storage of genetic information. Interest has grown in understanding how nucleic acids perform such sophisticated functions given their limited molecular repertoire. RNA can fold into diverse shapes that often perturb pKa values and allow it to ionize appreciably under biological conditions, thereby extending its molecular diversity. The goal of this chapter is to enable experimental measurement of pKa's in RNA and DNA. A number of experimental methods for measuring pKa values in RNA and DNA have been developed over the last 10 years, including RNA cleavage kinetics; UV-, fluorescence-, and NMR-detected pH titrations; and Raman crystallography. We begin with general considerations for choosing a pKa assay and then describe experimental conditions, advantages, and disadvantages for these assays. Potential pitfalls in measuring a pKa are provided including the presence of apparent pKa's due to a kinetic pKa or coupled acid- and alkali-promoted RNA unfolding, as well as degradation of RNA, precipitation of metal hydroxides and poor baselines. Use of multiple data fitting procedures and the study of appropriate mutants are described as ways to avoid some of these pitfalls. Application of these experimental methods to RNA and DNA will increase the number of available nucleic acid pKa values in the literature, which should deepen insight into biology and provide benchmarks for pKa calculations. Future directions for measuring pKa's in nucleic acids are discussed.

Enhanced group II intron retrohoming in magnesium-deficient Escherichia coli via selection of mutations in the ribozyme core

Mobile group II introns are bacterial retrotransposons thought to be evolutionary ancestors of spliceosomal introns and retroelements in eukaryotes. They consist of a catalytically active intron RNA ("ribozyme") and an intron-encoded reverse transcriptase, which function together to promote RNA splicing and intron mobility via reverse splicing of the intron RNA into new DNA sites ("retrohoming"). Although group II introns are active in bacteria, their natural hosts, they function inefficiently in eukaryotes, where lower free Mg(2+) concentrations decrease their ribozyme activity and constitute a natural barrier to group II intron proliferation within nuclear genomes. Here, we show that retrohoming of the Ll.LtrB group II intron is strongly inhibited in an Escherichia coli mutant lacking the Mg(2+) transporter MgtA, and we use this system to select mutations in catalytic core domain V (DV) that partially rescue retrohoming at low Mg(2+) concentrations. We thus identified mutations in the distal stem of DV that increase retrohoming efficiency in the MgtA mutant up to 22-fold. Biochemical assays of splicing and reverse splicing indicate that the mutations increase the fraction of intron RNA that folds into an active conformation at low Mg(2+) concentrations, and terbium-cleavage assays suggest that this increase is due to enhanced Mg(2+) binding to the distal stem of DV. Our findings indicate that DV is involved in a critical Mg(2+)-dependent RNA folding step in group II introns and demonstrate the feasibility of selecting intron variants that function more efficiently at low Mg(2+) concentrations, with implications for evolution and potential applications in gene targeting.

Metal ion as both a cofactor and a probe of metal-binding sites in a uranyl-specific DNAzyme: a uranyl photocleavage study

DNAzymes are known to bind metal ions specifically to carry out catalytic functions. Despite many studies since DNAzymes were discovered nearly two decades ago, the metal-binding sites in DNAzymes are not fully understood. Herein, we adopt uranyl photocleavage to probe specific uranyl-binding sites in the 39E DNAzyme with catalytically relevant concentrations of uranyl. The results indicate that uranyl binds between T23 and C25 in the bulge loop, G11 and T12 in the stem loop of the enzyme strand, as well as between T2.4 and G3 close to the cleavage site in the substrate strand. Control experiments using two 39E DNAzyme mutants revealed a different cleavage pattern of the mutated region. Another DNAzyme, the 8–17 DNAzyme, which has a similar secondary structure but shows no activity in the presence of uranyl, indicated a different uranyl-dependent photocleavage as well. In addition, a close correlation between the concentration-dependent photocleavage and enzymatic activities is also demonstrated. Together, these experiments suggest that uranyl photocleavage has been successfully used to probe catalytically relevant uranyl-binding sites in the 39E DNAzyme. As uranyl is the cofactor of the 39E DNAzyme as well as the probe, specific uranyl binding has now been identified without disruption of the structure.

Analysis of diffuse K+ and Mg2+ ion binding to a two-base-pair kissing complex by single-molecule mechanical unfolding

The folding and stability of RNA tertiary interactions depend critically on cationic conditions. It is usually difficult, however, to isolate such effects on tertiary interactions from those on the entire RNA. By manipulating conformations of single RNA molecules using optical tweezers, we distinguished individual steps of breaking and forming of a two-base-pair kissing interaction from those of secondary folding. The binding of metal ions to the small tertiary structure appeared to be saturable with an apparent Kd of 160 mM for K(+) and 1.5 mM for Mg(2+). The kissing formation was estimated to be associated with binding of ~2-3 diffuse K(+) or Mg(2+) ions. At their saturated binding, Mg(2+) provided ~3 kcal/mol more stabilizing energy to the structure than K(+). Furthermore, the cations change the unkissing forces significantly more than the kissing ones. For example, the presence of Mg(2+) ions increased the average unkissing force from 21 pN to 44 pN, surprisingly high for breaking merely two base pairs; in contrast, the mean kissing force was changed by only 4.5 pN. Interestingly, the differential salt effects on the transition forces were not caused by different changes in the height of the kinetic barriers but were instead attributed to how different molecular structures respond to the applied force. Our results showed the importance of diffuse cation binding to the stability of tertiary interaction and demonstrated the utility of mechanical unfolding in studying tertiary interactions.

Identification of the catalytic Mg²⁺ ion in the hepatitis delta virus ribozyme

The hepatitis delta virus ribozyme catalyzes an RNA cleavage reaction using a catalytic nucleobase and a divalent metal ion. The catalytic base, C75, serves as a general acid and has a pK(a) shifted toward neutrality. Less is known about the role of metal ions in the mechanism. A recent crystal structure of the precleavage ribozyme identified a Mg²⁺ ion that interacts through its partial hydration sphere with the G25·U20 reverse wobble. In addition, this Mg²⁺ ion is in position to directly coordinate the nucleophile, the 2'-hydroxyl of U(-1), suggesting it can serve as a Lewis acid to facilitate deprotonation of the 2'-hydroxyl. To test the role of the active site Mg²⁺ ion, we replaced the G25·U20 reverse wobble with an isosteric A25·C20 reverse wobble. This change was found to significantly reduce the negative potential at the active site, as supported by electrostatics calculations, suggesting that active site Mg²⁺ binding could be adversely affected by the mutation. The kinetic analysis and molecular dynamics of the A25·C20 double mutant suggest that this variant stably folds into an active structure. However, pH-rate profiles of the double mutant in the presence of Mg²⁺ are inverted relative to the profiles for the wild-type ribozyme, suggesting that the A25·C20 double mutant has lost the active site metal ion. Overall, these studies support a model in which the partially hydrated Mg²⁺ positioned at the G25·U20 reverse wobble is catalytic and could serve as a Lewis acid, a Brønsted base, or both to facilitate deprotonation of the nucleophile.

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.

The importance of peripheral sequences in determining the metal selectivity of an in vitro-selected Co(2+) -dependent DNAzyme

DNAzymes are catalytically active DNA molecules that use metal cofactors for their enzymatic functions. While a growing number of DNAzymes with diverse functions and metal selectivities have been reported, the relationships between metal ion selectivity, conserved sequences and structures responsible for selectivity remain to be elucidated. To address this issue, we report biochemical assays of a family of previously reported in vitro selected DNAzymes. This family includes the clone 11 DNAzyme, which was isolated by positive and negative selection, and the clone 18 DNAzyme, which was isolated by positive selection alone. The clone 11 DNAzyme has a higher selectivity for Co(2+) over Pb(2+) compared with clone 18. The reasons for this difference are explored here through phylogenetic comparison, mutational analysis and stepwise truncation. A novel DNAzyme truncation method incorporated a nick in the middle of the DNAzyme to allow for truncation close to the nicked site while preserving peripheral sequences at both ends of the DNAzyme. The results demonstrate that peripheral sequences within the substrate binding arms, most notably the stem loop, loop II, are sufficient to restore its selectivity for Co(2+) over Pb(2+) to levels observed in clone 11. A comparison of these sequences' secondary structures and Co(2+) selectivities suggested that metastable structures affect metal ion selectivity. The Co(2+) selectivity of the clone 11 DNAzyme showed that the metal ion binding and selectivities of small, in vitro selected DNAzymes may be more complex than previously appreciated, and that clone 11 may be more similar to larger ribozymes than to other small DNAzymes in its structural complexity and behavior. These factors should be taken into account when metal-ion selectivity is required in rationally designed DNAzymes and DNAzyme-based biosensors.

Ground-state coordination of a catalytic metal to the scissile phosphate of a tertiary-stabilized Hammerhead ribozyme

Although the Hammerhead ribozyme (HHRz) has long been used as a model system in the field of ribozyme enzymology, several details of its mechanism are still not well understood. In particular, significant questions remain concerning the disposition and role of catalytic metals in the HHRz. Previous metal-rescue experiments using a "minimal" HHRz resulted in prediction of a catalytic metal that is bound in the A9/G10.1 site in the ground state of the reaction and that bridges to the scissile phosphate further along the reaction pathway. "Native" or extended HHRz constructs contain tertiary contacts that stabilize a more compact structure at moderate ionic strength. We performed Cd(2+) rescue experiments on an extended HHRz from Schistosoma mansoni using stereo-pure scissile phosphorothioate-substituted substrates in order to determine whether a metal ion makes contact with the scissile phosphate in the ground state or further along the reaction coordinate. Inhibition in Ca(2+)/Mg(2+) and rescue by thiophilic Cd(2+) was specific for the R(p)-S stereoisomer of the scissile phosphate. The affinity of the rescuing Cd(2+), measured in two different ionic strength backgrounds, increased fourfold to 17-fold when the pro-R(p) oxygen is replaced by sulfur. These data support a model in which the rescuing metal ion makes a ground-state interaction with the scissile phosphate in the native HHRz. The resulting model for Mg(2+) activation in the HHRz places a metal ion in contact with the scissile phosphate, where it may provide ground-state electrostatic activation of the substrate.

Metal-induced specific and nonspecific oligonucleotide folding studied by FRET and related biophysical and bioanalytical implications

Metal induced nucleic acid folding has been extensively studied with ribozymes, DNAzymes, tRNA and riboswitches. These RNA/DNA molecules usually have a high content of double-stranded regions to support a rigid scaffold. On the other hand, such rigid structural features are not available for many in vitro selected or rationally designed DNA aptamers; they adopt flexible random coil structures in the absence of target molecules. Upon target binding, these aptamers adaptively fold into a compact structure with a reduced end-to-end distance, making fluorescence resonance energy transfer (FRET) a popular signaling mechanism. However, nonspecific folding induced by mono- or divalent metal ions can also reduce the end-to-end distance and thus lead to false positive results. In this study we used a FRET pair labeled Hg(II) binding DNA and monitored metal-induced folding in the presence of various cations. While nonspecific electrostatically mediated folding can be very significant, at each tested salt condition, Hg(II) induced folding was still observed with a similar sensitivity. We also studied the biophysical meaning of the acceptor/donor fluorescence ratio that allowed us to explain the experimental observations. Potential solutions for this ionic strength problem have been discussed. For example, probes designed to signal the formation of double-stranded DNA showed a lower dependency on ionic strength.

Metal-ion-dependent folding of a uranyl-specific DNAzyme: insight into function from fluorescence resonance energy transfer studies

Fluorescence resonance energy transfer (FRET) has been used to study the global folding of an uranyl (UO(2)(2+))-specific 39E DNAzyme in the presence of Mg(2+), Zn(2+), Pb(2+), or UO(2)(2+). At pH 5.5 and physiological ionic strength (100 mM Na(+)), two of the three stems in this DNAzyme folded into a compact structure in the presence of Mg(2+) or Zn(2+). However, no folding occurred in the presence of Pb(2+) or UO(2)(2+); this is analogous to the "lock-and-key" catalysis mode first observed in the Pb(2+)-specific 8-17 DNAzyme. However, Mg(2+) and Zn(2+) exert different effects on the 8-17 and 39E DNAzymes. Whereas Mg(2+) or Zn(2+)-dependent folding promoted 8-17 DNAzyme activity, the 39E DNAzyme folding induced by Mg(2+) or Zn(2+) inhibited UO(2)(2+)-specific activity. Group IIA series of metal ions (Mg(2+), Ca(2+), Sr(2+)) also caused global folding of the 39E DNAzyme, for which the apparent binding affinity between these metal ions and the DNAzyme decreases as the ionic radius of the metal ions increases. Because the ionic radius of Sr(2+) (1.12 Å) is comparable to that of Pb(2+) (1.20 Å), but contrary to Pb(2+), Sr(2+) induces the DNAzyme to fold under identical conditions, ionic size alone cannot account for the unique folding behaviors induced by Pb(2+) and UO(2)(2+). Under low ionic strength (30 mM Na(+)), all four metal ions (Mg(2+), Zn(2+), Pb(2+), and UO(2)(2+)), caused 39E DNAzyme folding, suggesting that metal ions can neutralize the negative charge of DNA-backbone phosphates in addition to playing specific catalytic roles. Mg(2+) at low (<2 mM) concentration promoted UO(2)(2+)-specific activity, whereas Mg(2+) at high (>2 mM) concentration inhibited the UO(2)(2+)-specific activity. Therefore, the lock-and-key mode of DNAzymes depends on ionic strength, and the 39E DNAzyme is in the lock-and-key mode only at ionic strengths of 100 mM or greater.

Two distinct catalytic strategies in the hepatitis δ virus ribozyme cleavage reaction

The hepatitis delta virus (HDV) ribozyme and related RNAs are widely dispersed in nature. This RNA is a small nucleolytic ribozyme that self-cleaves to generate products with a 2',3'-cyclic phosphate and a free 5'-hydroxyl. Although small ribozymes are dependent on divalent metal ions under biologically relevant buffer conditions, they function in the absence of divalent metal ions at high ionic strengths. This characteristic suggests that a functional group within the covalent structure of small ribozymes is facilitating catalysis. Structural and mechanistic analyses have demonstrated that the HDV ribozyme active site contains a cytosine with a perturbed pK(a) that serves as a general acid to protonate the leaving group. The reaction of the HDV ribozyme in monovalent cations alone never approaches the velocity of the Mg(2+)-dependent reaction, and there is significant biochemical evidence that a Mg(2+) ion participates directly in catalysis. A recent crystal structure of the HDV ribozyme revealed that there is a metal binding pocket in the HDV ribozyme active site. Modeling of the cleavage site into the structure suggested that this metal ion can interact directly with the scissile phosphate and the nucleophile. In this manner, the Mg(2+) ion can serve as a Lewis acid, facilitating deprotonation of the nucleophile and stabilizing the conformation of the cleavage site for in-line attack of the nucleophile at the scissile phosphate. This catalytic strategy had previously been observed only in much larger ribozymes. Thus, in contrast to most large and small ribozymes, the HDV ribozyme uses two distinct catalytic strategies in its cleavage reaction.

Brownian dynamics study of the association between the 70S ribosome and elongation factor G

Protein synthesis on the ribosome involves a number of external protein factors that bind at its functional sites. One key factor is the elongation factor G (EF-G) that facilitates the translocation of transfer RNAs between their binding sites, as well as advancement of the messenger RNA by one codon. The details of the EF-G/ribosome diffusional encounter and EF-G association pathway still remain unanswered. Here, we applied Brownian dynamics methodology to study bimolecular association in the bacterial EF-G/70S ribosome system. We estimated the EF-G association rate constants at 150 and 300 mM monovalent ionic strengths and obtained reasonable agreement with kinetic experiments. We have also elucidated the details of EF-G/ribosome association paths and found that positioning of the L11 protein of the large ribosomal subunit is likely crucial for EF-G entry to its binding site.

Mechanistic strategies in the HDV ribozyme: chelated and diffuse metal ion interactions and active site protonation

The crystal structure of the precleaved form of the hepatitis delta virus (HDV) ribozyme reveals two G•U wobbles near the active site: a rare reverse G•U wobble involving a syn G base, and a standard G•U wobble at the cleavage site. The catalytic mechanism for this ribozyme has been proposed to involve a Mg(2+) ion bound to the reverse G•U wobble, as well as a protonated C75 base. We carried out molecular dynamics simulations to analyze metal ion interaction with the reverse and standard G•U wobbles and to investigate the impact of C75 protonation on the structure and motions of the ribozyme. We identified two types of Mg(2+) ions associated with the ribozyme, chelated and diffuse, at the reverse and standard G•U wobbles, respectively, which appear to contribute to catalysis and stability, respectively. These two metal ion sites exhibit relatively independent behavior. Protonation of C75 was observed to locally organize the active site in a manner that facilitates the catalytic mechanism, in which C75(+) acts as a general acid and Mg(2+) as a Lewis acid. The simulations also indicated that the overall structure and thermal motions of the ribozyme are not significantly influenced by the catalytic Mg(2+) interaction or C75 protonation. This analysis suggests that the reaction pathway of the ribozyme is dominated by small local motions at the active site rather than large-scale global conformational changes. These results are consistent with a wealth of experimental data.

Metal binding motif in the active site of the HDV ribozyme binds divalent and monovalent ions

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G·U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg(2+) ion to participate in catalysis. Herein, we used molecular dynamics (MD) and X-ray crystallography to characterize the conformation and metal binding characteristics of this base pair in product and precleaved forms. Beginning with a crystal structure of the product form, we observed formation of the reverse G·U wobble during MD trajectories. We also demonstrated that this base pair is compatible with the diffraction data for the product-bound state. During MD trajectories of the product form, Na(+) ions interacted with the reverse G·U wobble in the RNA active site, and a Mg(2+) ion, introduced in certain trajectories, remained bound at this site. Beginning with a crystal structure of the precleaved form, the reverse G·U wobble with bound Mg(2+) remained intact during MD simulations. When we removed Mg(2+) from the starting precleaved structure, Na(+) ions interacted with the reverse G·U wobble. In support of the computational results, we observed competition between Na(+) and Mg(2+) in the precleaved ribozyme crystallographically. Nonlinear Poisson-Boltzmann calculations revealed a negatively charged patch near the reverse G·U wobble. This anionic pocket likely serves to bind metal ions and to help shift the pK(a) of the catalytic nucleobase, C75. Thus, the reverse G·U wobble motif serves to organize two catalytic elements, a metal ion and catalytic nucleobase, within the active site of the HDV ribozyme.

Rapid steps in the glmS ribozyme catalytic pathway: cation and ligand requirements

The glmS ribozyme is a conserved riboswitch found in numerous Gram-positive bacteria and responds to the cellular concentrations of glucosamine 6-phosphate (GlcN6P). GlcN6P binding promotes site-specific self-cleavage in the 5' UTR of the glmS mRNA, resulting in downregulation of gene expression. The glmS ribozyme has previously been shown to lack strong cation specificity when the rate-limiting folding step of the cleavage reaction pathway is measured. This does not provide data regarding cation and ligand specificities of the glmS ribozyme during the rapid ligand binding chemical catalysis events. Prefolding of the ribozyme in Mg(2+)-containing buffers effectively isolates the rapid ligand binding and catalytic events (k(obs) > 60 min(-1)) from rate-limiting folding (k(obs) < 4 min(-1)). Here we employ this experimental design to assay the cations and ligand requirements for rapid ligand binding and catalysis. We show that molar concentrations of monovalent cations are also capable of inducing the formation of the native GlcN6P binding structure but are unable to promote ligand binding and catalysis rates of >4 min(-1). Our data show that the sole obligatory role for divalent cations, for which there is crystallographic evidence, is coordination of the phosphate moiety of GlcN6P in the ligand-binding pocket. In further support of this hypothesis, our data show that a nonphosphorylated analogue of GlcN6P, glucosamine, is unable to promote rapid ligand binding and catalysis in the presence of divalent cations. Folding of the ribozyme is, therefore, relatively independent of cation identity, but the rapid initiation of catalysis upon the addition of ligand is stricter.

Activity, folding and Z-DNA formation of the 8-17 DNAzyme in the presence of monovalent ions

The effect of monovalent ions on both the reactivity and global folding of the 8-17 DNAzyme is investigated, and the results are compared with those of the hammerhead ribozyme, which has similar size and secondary structure. In contrast to the hammerhead ribozyme, the 8-17 DNAzyme activity is not detectable in the presence of 4 M K(+), Rb(+), or Cs(+) or in the presence of 80 mM, [Co(NH(3))(6)](3+). Only 4 M Li(+), NH(4)(+) and, to a lesser extent, Na(+) conferred detectable activity. The observed rate constants (k(obs) approximately 10(-3) min(-1) for Li(+) and NH(4)(+)) are approximately 1000-fold lower than that in the presence of 10 mM Mg(2+), and approximately 200,000-fold slower than that in the presence of 100 microM Pb(2+). Since the hammerhead ribozyme displays monovalent ion-dependent activity that is often within approximately 10-fold of divalent metal ion-dependent activity, these results suggest that the 8-17 DNAzyme, obtained by in vitro selections, has evolved to have a more stringent divalent metal ion requirement for high activity as compared to the naturally occurring ribozymes, making the 8-17 DNAzyme an excellent choice as a Pb(2+) sensor with high selectivity. In contrast to the activity data, folding was observed in the presence of all the monovalent ions investigated, although those monovalent ions that do not support DNAzyme activity have weaker binding affinity (K(d) approximately 0.35 M for Rb(+) and Cs(+)), while those that confer DNAzyme activity possess stronger affinity (K(d) approximately 0.22 M for Li(+), Na(+) and NH(4)(+)). In addition, a correlation between metal ion charge density, binding affinity and enzyme activity was found among mono- and divalent metal ions except Pb(2+); higher charge density resulted in stronger affinity and higher activity, suggesting that the observed folding and activity is at least partially due to electrostatic interactions between ions and the DNAzyme. Finally, circular dichroism (CD) study has revealed Z-DNA formation with the monovalent metal ions, Zn(2+) and Mg(2+); the K(d) values obtained using CD were in the same range as those obtained from folding studies using FRET. However, Z-DNA formation was not observed with Pb(2+). These results indicate that Pb(2+)-dependent function follows a different mechanism from the monovalent metal ions and other divalent metal ions; in the presence of latter metal ions, metal-ion dependent folding and structural changes, including formation of Z-DNA, play an important role in the catalytic function of the 8-17 DNAzyme.

Computational exploration of mobile ion distributions around RNA duplex

Atomically detailed distributions of ions around an A-form RNA are computed. Different mixtures of monovalent and divalent ions are considered explicitly. Studies of tightly bound and of diffusive (but bound) ions around 25 base pairs RNA are conducted in explicit solvent. Replica exchange simulations provide detailed equilibrium distributions with moderate computing resources (20 ns of simulation using 64 replicas). The simulations show distinct behavior of single and double charged cations. Binding of Mg(2+) ion includes tight binding to specific sites while Na(+) binds only diffusively. The tight binding of Mg(2+) is with a solvation shell while Na(+) can bind directly to RNA. Negative mobile ions can be found near the RNA but must be assisted by proximate and mobile cations. At distances larger than 16 A from the RNA center, a model of RNA as charged rod in a continuum of ionic solution provides quantitative description of the ion density (the same as in atomically detailed simulation). At shorter distances, the structure of RNA (and ions) has a significant impact on the pair correlation functions. Predicted binding sites of Mg(2+) at the RNA surface are in accord with structures from crystallography. Electric field relaxation is investigated. The relaxation due to solution rearrangements is completed in tens of picoseconds, while the contribution of RNA tumbling continues to a few nanoseconds.