Cooperative Interactions in the Hammerhead Ribozyme Drive pKa Shifting of G12 and Its Stacked Base C17

Scaling Catalytic Contributions of Small Self-Cleaving Ribozymes

Nucleolytic ribozymes utilize general acid-base catalysis to perform phosphodiester cleavage. In most ribozyme classes, a conserved active site guanosine is positioned to act as general base, thereby activating the 2'-OH group to attack the scissile phosphate (γ-catalysis). Here, we present an atomic mutagenesis study for the pistol ribozyme class. Strikingly, "general base knockout" by replacement of the guanine N1 atom by carbon results in only 2.7-fold decreased rate. Therefore, the common view that γ-catalysis critically depends on the N1 moiety becomes challenged. For pistol ribozymes we found that γ-catalysis is subordinate in overall catalysis, made up by two other catalytic factors (α and δ). Our approach allows scaling of the different catalytic contributions (α, β, γ, δ) with unprecedented precision and paves the way for a thorough mechanistic understanding of nucleolytic ribozymes with active site guanines.

Scaling Catalytic Contributions of Small Self-Cleaving Ribozymes

Nucleolytic ribozymes utilize general acid-base catalysis to perform phosphodiester cleavage. In most ribozyme classes, a conserved active site guanosine is positioned to act as general base, thereby activating the 2'-OH group to attack the scissile phosphate (γ-catalysis). Here, we present an atomic mutagenesis study for the pistol ribozyme class. Strikingly, "general base knockout" by replacement of the guanine N1 atom by carbon results in only 2.7-fold decreased rate. Therefore, the common view that γ-catalysis critically depends on the N1 moiety becomes challenged. For pistol ribozymes we found that γ-catalysis is subordinate in overall catalysis, made up by two other catalytic factors (α and δ). Our approach allows scaling of the different catalytic contributions (α, β, γ, δ) with unprecedented precision and paves the way for a thorough mechanistic understanding of nucleolytic ribozymes with active site guanines.

Solution Structure of NPSL2, A Regulatory Element in the oncomiR-1 RNA

The miR-17 ∼ 92a polycistron, also known as oncomiR-1, is commonly overexpressed in multiple cancers and has several oncogenic properties. OncomiR-1 encodes six constituent microRNAs (miRs), each enzymatically processed with different efficiencies. However, the structural mechanism that regulates this differential processing remains unclear. Chemical probing of oncomiR-1 revealed that the Drosha cleavage sites of pri-miR-92a are sequestered in a four-way junction. NPSL2, an independent stem loop element, is positioned just upstream of pri-miR-92a and sequesters a crucial part of the sequence that constitutes the basal helix of pri-miR-92a. Disruption of the NPSL2 hairpin structure could promote the formation of a pri-miR-92a structure that is primed for processing by Drosha. Thus, NPSL2 is predicted to function as a structural switch, regulating pri-miR-92a processing. Here, we determined the solution structure of NPSL2 using solution NMR spectroscopy. This is the first high-resolution structure of an oncomiR-1 element. NPSL2 adopts a hairpin structure with a large, but highly structured, apical and internal loops. The 10-bp apical loop contains a pH-sensitive A+·C mismatch. Additionally, several adenosines within the apical and internal loops have elevated pKa values. The protonation of these adenosines can stabilize the NPSL2 structure through electrostatic interactions. Our study provides fundamental insights into the secondary and tertiary structure of an important RNA hairpin proposed to regulate miR biogenesis.

The role of Na+ in catalysis by the 8-17 DNAzyme

The 8-17 DNAzyme is the most studied deoxyribozyme in terms of its molecular mechanism; hence it has become a model system to understand the basis behind DNA catalysis. New functional studies and the recent attainment of high-resolution X-ray structures, in addition to theoretical calculations have offered a great opportunity to gain a broader comprehension of its mechanism; however many aspects are unclear yet, especially regarding the precise role of metal ions in catalysis. Recently, molecular dynamics simulations have suggested for the first time a specific and dynamical participation of Na⁺ in the mechanism through the reaction pathway, besides the roles proposed for divalent metal cofactors. Herein, we present experimental evidence of a cooperative role of the monovalent cation Na⁺ in catalysis that is in line with these theoretical suggestions. Our findings show a clear influence of the concentration of Na⁺ on the activity of the 8-17 DNAzyme when Pb²⁺ is used as the cofactor. Interestingly, this effect is not noticed with Mg²⁺, indicating a particular contribution of the monovalent ion to catalysis that would operate preferentially with Pb²⁺. We have also found that Na⁺ affects the pK_a of the general base and the general acid, indicating its influence on general acid-base catalysis, already identified as part of the mechanism of the 8-17 DNAzyme. Finally, our results emphasize the need to consider Na⁺ carefully in the design and analysis of functional studies of catalytic DNAs and its possible specific role in their mechanisms.

Beyond the Plateau: pL Dependence of Proton Inventories as a Tool for Studying Ribozyme and Ribonuclease Catalysis

Acid/base catalysis is an important catalytic strategy used by ribonucleases and ribozymes; however, understanding the number and identity of functional groups involved in proton transfer remains challenging. The proton inventory (PI) technique analyzes the dependence of the enzyme reaction rate on the ratio of D₂O to H₂O and can provide information about the number of exchangeable sites that produce isotope effects and their magnitude. The Gross-Butler (GB) equation is used to evaluate H/D fractionation factors from PI data typically collected under conditions (i.e., a "plateau" in the pH-rate profile) assuming minimal change in active site residue ionization. However, restricting PI analysis to these conditions is problematic for many ribonucleases, ribozymes, and their variants due to ambiguity in the roles of active site residues, the lack of a plateau within the accessible pL range, or cooperative interactions between active site functional groups undergoing ionization. Here, we extend the integration of species distributions for alternative enzyme states in noncooperative models of acid/base catalysis into the GB equation, first used by Bevilacqua and colleagues for the HDV ribozyme, to develop a general population-weighted GB equation that allows simulation and global fitting of the three-dimensional relationship of the D₂O ratio (n) versus pL versus k_n/k₀. Simulations using the GPW-GB equation of PI results for RNase A, HDVrz, and VSrz illustrate that data obtained at multiple selected pL values across the pL-rate profile can assist in the planning and interpreting of solvent isotope effect experiments to distinguish alternative mechanistic models.

pH Effects and Cooperativity among Key Titratable Residues for Escherichia coli Glycinamide Ribonucleotide Transformylase

Human glycinamide ribonucleotide transformylase (GAR Tfase) is a regulatory enzyme in the de novo purine biosynthesis pathway that has been extensively studied as an anticancer target. To some extent, inhibition of GAR Tfase selectively targets cancer cells over normal cells and inhibits purine formation and DNA replication. In this study, we investigated E. coli GAR Tfase, which shares high sequence similarity with the human GAR Tfase, and most functional residues are conserved. Herein, we aim to predict the pH-activity curve through a computational approach. We carried out pH-replica exchange molecular dynamics (pH-REMD) simulations to investigate pH-dependent functions such as structural changes, ligand binding, and catalytic activity. To compute the pH-activity curve, we identified the catalytic residues in specific protonation states, referred to as the catalytic competent protonation states (CCPS), which maintain the structure, keep ligands bound, and facilitate catalysis. Our computed population of CCPS with respect to pH matches well with the experimental pH-activity curve. To compute the microscopic pK_a values in the catalytically active conformation, we devised a thermodynamic model that considers the coupling between protonation states of CCPS residues and conformational states. These results allow us to correctly identify the general acid and base catalysts and interpret the pH-activity curve at an atomistic level.

Hydrated metal ion as a general acid in the catalytic mechanism of the 8-17 DNAzyme

The RNA-cleaving 8-17 DNAzyme, which is a metalloenzyme that depends on divalent metal ions for its function, is the most studied catalytic DNA in terms of its mechanism. By the end of 2017, a report of the crystal structure of the enzyme-substrate complex in the presence of Pb²⁺ probed some of the previous findings and opened new questions, especially around the participation of the metal ion in the catalytic mechanism and the promiscuity exhibited by the enzyme in terms of the metal cofactor required for catalysis. In this article we explore the role of the divalent metal ion in the mechanism of the 8-17 DNAzyme as a general acid, by measuring the influence of pH over the activity of a slower variant of the enzyme in the presence of Pb²⁺. We replaced G14, which has been identified as a general base in the mechanism of the enzyme, by the unnatural analog 2-aminopurine, with a lower pK_a value of the N1 group. With this approach, we obtained a bell-shaped pH-rate profile with experimental pK_a values of 5.4 and 7.0. Comparing these results with previous pH-rate profiles in the presence of Mg²⁺, our findings suggest the stabilization of the 5'-O leaving group by the hydrated metal ion acting as a general acid, in addition to the activation of the 2'-OH nucleophile by the general base G14.

Prebiotically-relevant low polyion multivalency can improve functionality of membraneless compartments

Multivalent polyions can undergo complex coacervation, producing membraneless compartments that accumulate ribozymes and enhance catalysis, and offering a mechanism for functional prebiotic compartmentalization in the origins of life. Here, we evaluate the impact of lower, more prebiotically-relevant, polyion multivalency on the functional performance of coacervates as compartments. Positively and negatively charged homopeptides with 1–100 residues and adenosine mono-, di-, and triphosphate nucleotides are used as model polyions. Polycation/polyanion pairs are tested for coacervation, and resulting membraneless compartments are analyzed for salt resistance, ability to provide a distinct internal microenvironment (apparent local pH, RNA partitioning), and effect on RNA structure formation. We find that coacervates formed by phase separation of the shorter polyions more effectively generated distinct pH microenvironments, accumulated RNA, and preserved duplexes than those formed by longer polyions. Hence, coacervates formed by reduced multivalency polyions are not only viable as functional compartments for prebiotic chemistries, they can outperform higher molecular weight analogues.

Short cationic peptides and nucleotides can form complex coacervates, but the influence of reduced multivalency on coacervate functionality was not investigated. Here, the authors report that coacervates formed from short polyions generate distinct pH microenvironments, accumulate RNA and preserve nucleic acid duplexes more efficiently than their longer counterparts.

The L-platform/L-scaffold framework: a blueprint for RNA-cleaving nucleic acid enzyme design

We develop an L-platform/L-scaffold framework we hypothesize may serve as a blueprint to facilitate site-specific RNA-cleaving nucleic acid enzyme design. Building on the L-platform motif originally described by Suslov and coworkers, we identify new critical scaffolding elements required to anchor a conserved general base guanine ("L-anchor") and bind functionally important metal ions at the active site ("L-pocket"). Molecular simulations, together with a broad range of experimental structural and functional data, connect the L-platform/L-scaffold elements to necessary and sufficient conditions for catalytic activity. We demonstrate that the L-platform/L-scaffold framework is common to five of the nine currently known naturally occurring ribozyme classes (Twr, HPr, VSr, HHr, Psr), and intriguingly from a design perspective, the framework also appears in an artificially engineered DNAzyme (8-17dz). The flexibility of the L-platform/L-scaffold framework is illustrated on these systems, highlighting modularity and trends in the variety of known general acid moieties that are supported. These trends give rise to two distinct catalytic paradigms, building on the classifications proposed by Wilson and coworkers and named for the implicated general base and acid. The "G + A" paradigm (Twr, HPr, VSr) exclusively utilizes nucleobase residues for chemistry, and the "G + M + " paradigm (HHr, 8-17dz, Psr) involves structuring of the "L-pocket" metal ion binding site for recruitment of a divalent metal ion that plays an active role in the chemical steps of the reaction. Finally, the modularity of the L-platform/L-scaffold framework is illustrated in the VS ribozyme where the "L-pocket" assumes the functional role of the "L-anchor" element, highlighting a distinct mechanism, but one that is functionally linked with the hammerhead ribozyme.

Nucleic Acid Catalysis under Potential Prebiotic Conditions

Catalysis by nucleic acids is indispensable for extant cellular life, and it is widely accepted that nucleic acid enzymes were crucial for the emergence of primitive life 3.5-4 billion years ago. However, geochemical conditions on early Earth must have differed greatly from the constant internal milieus of today's cells. In order to explore plausible scenarios for early molecular evolution, it is therefore essential to understand how different physicochemical parameters, such as temperature, pH, and ionic composition, influence nucleic acid catalysis and to explore to what extent nucleic acid enzymes can adapt to non-physiological conditions. In this article, we give an overview of the research on catalysis of nucleic acids, in particular catalytic RNAs (ribozymes) and DNAs (deoxyribozymes), under extreme and/or unusual conditions that may relate to prebiotic environments.

Exposing Hidden High-Affinity RNA Conformational States

RNA recognition frequently results in conformational changes that optimize intermolecular binding. As a consequence, the overall binding affinity of RNA to its binding partners depends not only on the intermolecular interactions formed in the bound state but also on the energy cost associated with changing the RNA conformational distribution. Measuring these "conformational penalties" is, however, challenging because bound RNA conformations tend to have equilibrium populations in the absence of the binding partner that fall outside detection by conventional biophysical methods. In this study we employ as a model system HIV-1 TAR RNA and its interaction with the ligand argininamide (ARG), a mimic of TAR's cognate protein binding partner, the transactivator Tat. We use NMR chemical shift perturbations and relaxation dispersion in combination with Bayesian inference to develop a detailed thermodynamic model of coupled conformational change and ligand binding. Starting from a comprehensive 12-state model of the equilibrium, we estimate the energies of six distinct detectable thermodynamic states that are not accessible by currently available methods. Our approach identifies a minimum of four RNA intermediates that differ in terms of the TAR conformation and ARG occupancy. The dominant bound TAR conformation features two bound ARG ligands and has an equilibrium population in the absence of ARG that is below detection limit. Consequently, even though ARG binds to TAR with an apparent overall weak affinity (Kdapp ≈ 0.2 mM), it binds the prefolded conformation with a Kd in the nM range. Our results show that conformational penalties can be major determinants of RNA-ligand binding affinity as well as a source of binding cooperativity, with important implications for a predictive understanding of how RNA is recognized and for RNA-targeted drug discovery.

Insights into DNA catalysis from structural and functional studies of the 8-17 DNAzyme

The review examines functional knowledge gathered over two decades of research on the 8-17 DNAzyme, focusing on three aspects: the structural requirements for catalysis, the role of metal ions and the participation of general acid-base catalysis.

Dynamical ensemble of the active state and transition state mimic for the RNA-cleaving 8-17 DNAzyme in solution

We perform molecular dynamics simulations, based on recent crystallographic data, on the 8-17 DNAzyme at four states along the reaction pathway to determine the dynamical ensemble for the active state and transition state mimic in solution. A striking finding is the diverse roles played by Na+ and Pb2+ ions in the electrostatically strained active site that impact all four fundamental catalytic strategies, and share commonality with some features recently inferred for naturally occurring hammerhead and pistol ribozymes. The active site Pb2+ ion helps to stabilize in-line nucleophilic attack, provides direct electrostatic transition state stabilization, and facilitates leaving group departure. A conserved guanine residue is positioned to act as the general base, and is assisted by a bridging Na+ ion that tunes the pKa and facilitates in-line fitness. The present work provides insight into how DNA molecules are able to solve the RNA-cleavage problem, and establishes functional relationships between the mechanism of these engineered DNA enzymes with their naturally evolved RNA counterparts. This adds valuable information to our growing body of knowledge on general mechanisms of phosphoryl transfer reactions catalyzed by RNA, proteins and DNA.

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.

Cellular Concentrations of Nucleotide Diphosphate-Chelated Magnesium Ions Accelerate Catalysis by RNA and DNA Enzymes

RNAs are involved in myriad cellular events. In general, RNA function is affected by cellular conditions. For instance, molecular crowding promotes RNA folding through compaction of the RNA. Metabolites generally destabilize RNA secondary structure, which improves RNA folding cooperativity and increases the fraction of functional RNA. Our recent studies demonstrate that cellular concentrations of amino acid-chelated magnesium (aaCM) stimulate RNA folding and catalysis while protecting RNAs from magnesium ion-induced degradation. However, effects of other cellular magnesium ion chelators on RNA function have not been tested. Herein, we report that nucleotide diphosphate-chelated magnesium, which is of intermediate strength, promotes RNA catalysis much like aaCM. Nucleotides are some of the major metabolites in cells and have one to three phosphates, which have increasingly tight binding of magnesium. On the basis of binding calculations, ∼85% ATP, ∼40% ADP, and only 5% AMP are estimated to possess a magnesium ion under cellular conditions of 0.50 mM Mg²⁺_free. We tested the self-cleaving activity of the hammerhead ribozyme in the presence of these chelated magnesium species. Our results indicate that NTP-chelated magnesium and NMP-chelated magnesium do not appreciably stimulate RNA catalysis, whereas NDP-chelated magnesium promotes RNA catalysis up to 6.5-fold. Inspired by NDP, we observed similar stimulatory effects for several other Mg²⁺ diphosphate-containing metabolites, including UDP-GlcNAC and UDP-Glc; in addition, we found similar effects for a DNAzyme. Thus, rate stimulatory effects are general with respect to the diphosphate and nucleic acid enzyme. These results implicate magnesium-chelated diphosphate metabolites as general facilitators of RNA function inside cells.

Predicting Site-Binding Modes of Ions and Water to Nucleic Acids Using Molecular Solvation Theory

Site binding of ions and water shapes nucleic acids folding, dynamics, and biological function, complementing the more diffuse, nonspecific "territorial" ion binding. Unlike territorial binding, prediction of site-specific binding to nucleic acids remains an unsolved challenge in computational biophysics. This work presents a new toolset based on the 3D-RISM molecular solvation theory and topological analysis that predicts cation and water site binding to nucleic acids. 3D-RISM is shown to accurately capture alkali cations and water binding to the central channel, transversal loops, and grooves of the Oxytricha nova's telomeres' G-quadruplex ( Oxy-GQ), in agreement with high-resolution crystallographic data. To improve the computed cation occupancy along the Oxy-GQ central channel, it was necessary to refine and validate new cation-oxygen parameters using structural and thermodynamic data available for crown ethers and ion channels. This single set of parameters that describes both localized and delocalized binding to various biological systems is used to gain insight into cation occupancy along the Oxy-GQ channel under various salt conditions. The paper concludes with prospects for extending the method to predict divalent cation binding to nucleic acids. This work advances the forefront of theoretical methods able to provide predictive insight into ion atmosphere effects on nucleic acids function.

Bulky cations greatly increase the turnover of a native hammerhead ribozyme

Methods to facilitate the catalytic turnover of ribozymes are required for advancing oligonucleotide-based technologies. This study examined tetraalkylammonium ions for their ability to increase the efficiency of catalytic turnover of a native hammerhead ribozyme. Kinetic analysis showed that large tetraalkylammonium ions significantly increased the turnover rate of the ribozyme and was much more effective than poly(ethylene glycol) (PEG) and urea. The magnitude of the rate increase depended on the concentrations of Mg²⁺ and tetrapentylammonium ions, and the rate was enhanced by more than 180-fold at the optimal concentrations of these salts. The results provide physical insights into interactions of ribozymes with large cationic molecules through electrostatic forces and steric hindrance.

Large tetraalkylammonium ions increase the turnover rate of the ribozyme derived from an intronic ribozyme in the human genome. The rate can be enhanced by more than a hundred-fold at the optimal concentrations of Mg²⁺ and TPeA ions.

Structural Simplicity and Mechanistic Complexity in the Hammerhead Ribozyme

Natural or full-length hammerhead ribozymes are up to 1000-fold more active than their minimal counterparts that lack a complex tertiary interaction that pre-organizes and stabilizes the ribozyme active site, positioning RNA functional groups to facilitate acid-base catalysis. The recent discovery that a single tertiary contact (an AU Hoogsteen pair) between Stems I and II confers essentially all of the enhanced activity greatly simplifies our understanding of the structural requirements for hammerhead ribozyme activity. In contrast, the simplest mechanistic interpretations are challenged with the presentation of more complex alternatives. These alternatives are elucidated and critically analyzed in the context of several of the active hammerhead ribozyme structures now available.

Incorporating uracil and 5-halouracils into short peptide nucleic acids for enhanced recognition of A-U pairs in dsRNAs

Double-stranded RNA (dsRNA) structures form triplexes and RNA-protein complexes through binding to single-stranded RNA (ssRNA) regions and proteins, respectively, for diverse biological functions. Hence, targeting dsRNAs through major-groove triplex formation is a promising strategy for the development of chemical probes and potential therapeutics. Short (e.g., 6-10 mer) chemically-modified Peptide Nucleic Acids (PNAs) have been developed that bind to dsRNAs sequence specifically at physiological conditions. For example, a PNA incorporating a modified base thio-pseudoisocytosine (L) has an enhanced recognition of a G-C pair in an RNA duplex through major-groove L·G-C base triple formation at physiological pH, with reduced pH dependence as observed for C+·G-C base triple formation. Currently, an unmodified T base is often incorporated into PNAs to recognize a Watson-Crick A-U pair through major-groove T·A-U base triple formation. A substitution of the 5-methyl group in T by hydrogen and halogen atoms (F, Cl, Br, and I) causes a decrease of the pKa of N3 nitrogen atom, which may result in improved hydrogen bonding in addition to enhanced base stacking interactions. Here, we synthesized a series of PNAs incorporating uracil and halouracils, followed by binding studies by non-denaturing polyacrylamide gel electrophoresis, circular dichroism, and thermal melting. Our results suggest that replacing T with uracil and halouracils may enhance the recognition of an A-U pair by PNA·RNA2 triplex formation in a sequence-dependent manner, underscoring the importance of local stacking interactions. Incorporating bromouracils and chlorouracils into a PNA results in a significantly reduced pH dependence of triplex formation even for PNAs containing C bases, likely due to an upshift of the apparent pKa of N3 atoms of C bases. Thus, halogenation and other chemical modifications may be utilized to enhance hydrogen bonding of the adjacent base triples and thus triplex formation. Furthermore, our experimental and computational modelling data suggest that PNA·RNA2 triplexes may be stabilized by incorporating a BrUL step but not an LBrU step, in dsRNA-binding PNAs.

Cellular Small Molecules Contribute to Twister Ribozyme Catalysis

The number of self-cleaving small ribozymes has increased sharply in recent years. Advances have been made in describing these ribozymes in terms of four catalytic strategies: α describes in-line attack, β describes neutralization of the nonbridging oxygens, γ describes activation of the nucleophile, and δ describes stabilization of the leaving group. Current literature presents the rapid self-cleavage of the twister ribozyme in terms of all four of these classic catalytic strategies. Herein, we describe the nonspecific contribution of small molecules to ribozyme catalysis. At biological pH, the rate of the wild-type twister ribozyme is enhanced up to 5-fold in the presence of moderate buffer concentrations, similar to the 3-5-fold effects reported previously for buffer catalysis for protein enzymes. We observe this catalytic enhancement not only with standard laboratory buffers, but also with diverse biological small molecules, including imidazole, amino acids, and amino sugars. Brønsted plots suggest that small molecules assist in proton transfer, most likely with δ catalysis. Cellular small molecules provide a simple way to overcome the limited functional diversity of RNA and have the potential to participate in the catalytic mechanisms of many ribozymes in vivo.

Cellular conditions of weakly chelated magnesium ions strongly promote RNA stability and catalysis

Most RNA folding studies have been performed under non-physiological conditions of high concentrations (≥10 mM) of Mg²⁺_free, while actual cellular concentrations of Mg²⁺_free are only ~1 mM in a background of greater than 50 mM Mg²⁺_total. To uncover cellular behavior of RNA, we devised cytoplasm mimic systems that include biological concentrations of amino acids, which weakly chelate Mg²⁺. Amino acid-chelated Mg²⁺ (aaCM) of ~15 mM dramatically increases RNA folding and prevents RNA degradation. Furthermore, aaCM enhance self-cleavage of several different ribozymes, up to 100,000-fold at Mg²⁺_free of just 0.5 mM, indirectly through RNA compaction. Other metabolites that weakly chelate magnesium offer similar beneficial effects, which implies chelated magnesium may enhance RNA function in the cell in the same way. Overall, these results indicate that the states of Mg²⁺ should not be limited to free and bound only, as weakly bound Mg²⁺ strongly promotes RNA function under cellular conditions.

Impact of macromolecular crowding on RNA/spermine complex coacervation and oligonucleotide compartmentalization

We report the effect of neutral macromolecular crowders poly(ethylene glycol) (PEG) (8 kDa) and Ficoll (70 kDa) on liquid-liquid phase separation in a polyuridylic acid (polyU)/spermine complex coacervate system. The addition of PEG decreased both the amount of spermine required for phase separation and the coacervation temperature (T_C). We interpret these effects on phase behavior as arising due to excluded volume and preferential interactions on both the secondary structure/condensation of spermine-associated polyU molecules and on the association of soluble polyU/spermine polyelectrolyte complexes to form coacervate droplets. Examination of coacervates formed in the presence of fluorescently-labeled PEG or Ficoll crowders indicated that Ficoll is accumulated while PEG is excluded from the coacervate phase, which provides further insight into the differences in phase behavior. Crowding agents impact distribution of a biomolecular solute: partitioning of a fluorescently-labeled U15 RNA oligomer into the polyU/spermine coacervates was increased approximately two-fold by 20 wt% Ficoll 70 kDa and by more than two orders of magnitude by 20 wt% PEG 8 kDa. The volume of the coacervate phase decreased in the presence of crowder relative to a dilute buffer solution. These findings indicate that potential impacts of macromolecular crowding on phase behavior and solute partitioning should be considered in model systems for intracellular membraneless organelles.

Elucidation of Catalytic Strategies of Small Nucleolytic Ribozymes From Comparative Analysis of Active Sites

A number of small, self-cleaving ribozyme classes have been identified including the hammerhead, hairpin, hepatitis delta virus (HDV), Varkud satellite (VS), glmS, twister, hatchet, pistol, and twister sister ribozymes. Within the active sites of these ribozymes, myriad functional groups contribute to catalysis. There has been extensive structure-function analysis of individual ribozymes, but the extent to which catalytic devices are shared across different ribozyme classes is unclear. As such, emergent catalytic principles for ribozymes may await discovery. Identification of conserved catalytic devices can deepen our understanding of RNA catalysis specifically and of enzymic catalysis generally. To probe similarities and differences amongst ribozyme classes, active sites from more than 80 high-resolution crystal structures of self-cleaving ribozymes were compared computationally. We identify commonalities amongst ribozyme classes pertaining to four classic catalytic devices: deprotonation of the 2'OH nucleophile (γ), neutralization of the non-bridging oxygens of the scissile phosphate (β), neutralization of the O5' leaving group (δ), and in-line nucleophilic attack (α). In addition, we uncover conservation of two catalytic devices, each of which centers on the activation of the 2'OH nucleophile by a guanine: one to acidify the 2'OH by hydrogen bond donation to it (γ') and one to acidify the 2'OH by releasing it from non-productive interactions by competitive hydrogen bonding (γ''). Our findings reveal that the amidine functionalities of G, A, and C are especially important for these strategies, and help explain absence of U at ribozyme active sites. The identified γ' and γ'' catalytic strategies help unify the catalytic strategies shared amongst catalytic RNAs and may be important for large ribozymes, as well as protein enzymes that act on nucleic acids.

Corroboration of Zn(ii)–Mg(ii)-tertiary structure interplays essential for the optimal catalysis of a phosphorothiolate thiolesterase ribozyme

The TW17 ribozyme, a catalytic RNA selected from a pool of artificial RNA, is specific for the Zn²⁺-dependent hydrolysis of a phosphorothiolate thiolester bond. Here, we describe the organic synthesis of both guanosine α-thio-monophosphate and the substrates required for selecting and characterizing the TW17 ribozyme, and for deciphering the catalytic mechanism of the ribozyme. By successively substituting the substrate originally conjugated to the RNA pool with structurally modified substrates, we demonstrated that the TW17 ribozyme specifically catalyzes phosphorothiolate thiolester hydrolysis. Metal titration studies of TW17 ribozyme catalysis in the presence of Zn²⁺ alone, Zn²⁺ and Mg²⁺, and Zn²⁺ and [Co(NH₃)₆]³⁺ supported our findings that Zn²⁺ is absolutely required for ribozyme catalysis, and indicated that optimal ribozyme catalysis involves the presence of outer-sphere and one inner-sphere Mg²⁺. A survey of the TW17 ribozyme activity at various pHs revealed that the activity of the ribozyme critically depends on the alkaline conditions. Moreover, a GNRA tetraloop-containing ribozyme constructed with active catalysis in trans provided catalysis and multiple substrate turnover efficiencies significantly higher than ribozymes lacking a GNRA tetraloop. This research supports the essential roles of Zn²⁺, Mg²⁺, and a GNRA tetraloop in modulating the TW17 ribozyme structure for optimal ribozyme catalysis, leading also to the formulation of a proposed reaction mechanism for TW17 ribozyme catalysis.

Zn(ii) and Mg(ii) and GAGA tetraloop in the ion atmosphere of the TW17 ribozyme is critical to optimal ribozyme catalysis at alkaline pH.

Complexity in pH-Dependent Ribozyme Kinetics: Dark pKa Shifts and Wavy Rate-pH Profiles

Charged bases occur in RNA enzymes, or ribozymes, where they play key roles in catalysis. Cationic bases donate protons and perform electrostatic catalysis, while anionic bases accept protons. We previously published simulations of rate-pH profiles for ribozymes in terms of species plots for the general acid and general base that have been useful for understanding how ribozymes respond to pH. In that study, we did not consider interaction between the general acid and general base or interaction with other species on the RNA. Since that report, diverse small ribozyme classes have been discovered, many of which have charged nucleobases or metal ions in the active site that can either directly interact and participate in catalysis or indirectly interact as "influencers". Herein, we simulate experimental rate-pH profiles in terms of species plots in which reverse protonated charged nucleobases interact. These analyses uncover two surprising features of pH-dependent enzyme kinetics. (1) Cooperativity between the general acid and general base enhances population of the functional forms of a ribozyme and manifests itself as hidden or "dark" pK_a shifts, real pK_a shifts that accelerate the reaction but are not readily observed by standard experimental approaches, and (2) influencers favorably shift the pK_as of proton-transferring nucleobases and manifest themselves as "wavy" rate-pH profiles. We identify parallels with the protein enzyme literature, including reverse protonation and wavelike behavior, while pointing out that RNA is more prone to reverse protonation. The complexities uncovered, which arise from simple pairwise interactions, should aid deconvolution of complex rate-pH profiles for RNA and protein enzymes and suggest veiled catalytic devices for promoting catalysis that can be tested by experiment and calculation.

Structure-based insights into self-cleavage by a four-way junctional twister-sister ribozyme

Here we report on the crystal structure and cleavage assays of a four-way junctional twister-sister self-cleaving ribozyme. Notably, 11 conserved spatially separated loop nucleotides are brought into close proximity at the ribozyme core through long-range interactions mediated by hydrated Mg²⁺ cations. The C62–A63 step at the cleavage site adopts a splayed-apart orientation, with flexible C62 directed outwards, whereas A63 is directed inwards and anchored by stacking and hydrogen-bonding interactions. Structure-guided studies of key base, sugar, and phosphate mutations in the twister-sister ribozyme, suggest contributions to the cleavage chemistry from interactions between a guanine at the active site and the non-bridging oxygen of the scissile phosphate, a feature found previously also for the related twister ribozyme. Our four-way junctional pre-catalytic structure differs significantly in the alignment at the cleavage step (splayed-apart vs. base-stacked) and surrounding residues and hydrated Mg²⁺ ions relative to a reported three-way junctional pre-catalytic structure of the twister-sister ribozyme.

Twister-sister is a self-cleaving ribozyme. Here, the authors report the 2.0 Å crystal structure of the four-way junctional twister-sister ribozyme in the pre-catalytic state and discuss mechanistic implications based on their mutagenesis experiments and comparisons with other ribozyme structures.

Divalent Metal Ion Activation of a Guanine General Base in the Hammerhead Ribozyme: Insights from Molecular Simulations

The hammerhead ribozyme is a well-studied nucleolytic ribozyme that catalyzes the self-cleavage of the RNA phosphodiester backbone. Despite experimental and theoretical efforts, key questions remain about details of the mechanism with regard to the activation of the nucleophile by the putative general base guanine (G12). Straightforward interpretation of the measured activity-pH data implies the pK_a value of the N1 position in the G12 nucleobase is significantly shifted by the ribozyme environment. Recent crystallographic and biochemical work has identified pH-dependent divalent metal ion binding at the N7/O6 position of G12, leading to the hypothesis that this binding mode could induce a pK_a shift of G12 toward neutrality. We present computational results that support this hypothesis and provide a model that unifies the interpretation of available structural and biochemical data and paints a detailed mechanistic picture of the general base step of the reaction. Experimentally testable predictions are made for mutational and rescue effects on G12, which will give further insights into the catalytic mechanism. These results contribute to our growing knowledge of the potential roles of divalent metal ions in RNA catalysis.