Nucleic acid catalysis: metals, nucleobases, and other cofactors

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

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

Surface-Accelerated String Method for Locating Minimum Free Energy Paths

We present a surface-accelerated string method (SASM) to efficiently optimize low-dimensional reaction pathways from the sampling performed with expensive quantum mechanical/molecular mechanical (QM/MM) Hamiltonians. The SASM accelerates the convergence of the path using the aggregate sampling obtained from the current and previous string iterations, whereas approaches like the string method in collective variables (SMCV) or the modified string method in collective variables (MSMCV) update the path only from the sampling obtained from the current iteration. Furthermore, the SASM decouples the number of images used to perform sampling from the number of synthetic images used to represent the path. The path is optimized on the current best estimate of the free energy surface obtained from all available sampling, and the proposed set of new simulations is not restricted to being located along the optimized path. Instead, the umbrella potential placement is chosen to extend the range of the free energy surface and improve the quality of the free energy estimates near the path. In this manner, the SASM is shown to improve the exploration for a minimum free energy pathway in regions where the free energy surface is relatively flat. Furthermore, it improves the quality of the free energy profile when the string is discretized with too few images. We compare the SASM, SMCV, and MSMCV using 3 QM/MM applications: a ribozyme methyltransferase reaction using 2 reaction coordinates, the 2'-O-transphosphorylation reaction of Hammerhead ribozyme using 3 reaction coordinates, and a tautomeric reaction in B-DNA using 5 reaction coordinates. We show that SASM converges the paths using roughly 3 times less sampling than the SMCV and MSMCV methods. All three algorithms have been implemented in the FE-ToolKit package made freely available.

An RNA-Cleaving DNAzyme That Requires an Organic Solvent to Function

Engineering functional nucleic acids that are active under unusual conditions will not only reveal their hidden abilities but also lay the groundwork for pursuing them for unique applications. Although many DNAzymes have been derived to catalyze diverse chemical reactions in aqueous solutions, no prior study has been set up to purposely derive DNAzymes that require an organic solvent to function. Herein, we utilized in vitro selection to isolate RNA-cleaving DNAzymes from a random-sequence DNA pool that were "compelled" to accept 35 % dimethyl sulfoxide (DMSO) as a cosolvent, via counter selection in a purely aqueous solution followed by positive selection in the same solution containing 35 % DMSO. This experiment led to the discovery of a new DNAzyme that requires 35 % DMSO for its catalytic activity and exhibits drastically reduced activity without DMSO. This DNAzyme also requires divalent metal ions for catalysis, and its activity is enhanced by monovalent ions. A minimized, more efficient DNAzyme was also derived. This work demonstrates that highly functional, organic solvent-dependent DNAzymes can be isolated from random-sequence DNA libraries via forced in vitro selection, thus expanding the capability and potential utility of catalytic DNA.

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

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

Assessment of Density Functional Theory Methods for the Structural Prediction of Transition and Post-Transition Metal-Nucleic Acid Complexes

Understanding the structure of metal-nucleic acid systems is important for many applications such as the design of new pharmaceuticals, metal detection platforms, and nanomaterials. Herein, we explore the ability of 20 density functional theory (DFT) functionals to reproduce the crystal structure geometry of transition and post-transition metal-nucleic acid complexes identified in the Protein Data Bank and Cambridge Structural Database. The environmental extremes of the gas phase and implicit water were considered, and analysis focused on the global and inner coordination geometry, including the coordination distances. Although gas-phase calculations were unable to describe the structure of 12 out of the 53 complexes in our test set regardless of the DFT functional considered, accounting for the broader environment through implicit solvation or constraining the model truncation points to crystallographic coordinates generally afforded agreement with the experimental structure, suggesting that functional performance for these systems is likely due to the models rather than the methods. For the remaining 41 complexes, our results show that the reliability of functionals depends on the metal identity, with the magnitude of error varying across the periodic table. Furthermore, minimal changes in the geometries of these metal-nucleic acid complexes occur upon use of the Stuttgart-Dresden effective core potential and/or inclusion of an implicit water environment. The overall top three performing functionals are ωB97X-V, ωB97X-D3(BJ), and MN15, which reliably describe the structure of a broad range of metal-nucleic acid systems. Other suitable functionals include MN15-L, which is a cheaper alternative to MN15, and PBEh-3c, which is commonly used in QM/MM calculations of biomolecules. In fact, these five methods were the only functionals tested to reproduce the coordination sphere of Cu²⁺-containing complexes. For metal-nucleic acid systems that do not contain Cu²⁺, ωB97X and ωB97X-D are also suitable choices. These top-performing methods can be utilized in future investigations of diverse metal-nucleic acid complexes of relevance to biology and material science.

Rapid Kinetics of Pistol Ribozyme: Insights into Limits to RNA Catalysis

Pistol ribozyme (Psr) is a distinct class of small endonucleolytic ribozymes, which are important experimental systems for defining fundamental principles of RNA catalysis and designing valuable tools in biotechnology. High-resolution structures of Psr, extensive structure-function studies, and computation support a mechanism involving one or more catalytic guanosine nucleobases acting as a general base and divalent metal ion-bound water acting as an acid to catalyze RNA 2'-O-transphosphorylation. Yet, for a wide range of pH and metal ion concentrations, the rate of Psr catalysis is too fast to measure manually and the reaction steps that limit catalysis are not well understood. Here, we use stopped-flow fluorescence spectroscopy to evaluate Psr temperature dependence, solvent H/D isotope effects, and divalent metal ion affinity and specificity unconstrained by limitations due to fast kinetics. The results show that Psr catalysis is characterized by small apparent activation enthalpy and entropy changes and minimal transition state H/D fractionation, suggesting that one or more pre-equilibrium steps rather than chemistry is rate limiting. Quantitative analyses of divalent ion dependence confirm that metal aquo ion pK_a correlates with higher rates of catalysis independent of differences in ion binding affinity. However, ambiguity regarding the rate-limiting step and similar correlation with related attributes such as ionic radius and hydration free energy complicate a definitive mechanistic interpretation. These new data provide a framework for further interrogation of Psr transition state stabilization and show how thermal instability, metal ion insolubility at optimal pH, and pre-equilibrium steps such as ion binding and folding limit the catalytic power of Psr suggesting potential strategies for further optimization.

A Smart Metal-Polyphenol-DNAzyme nanoplatform for Gene-Chemodynamic Synergistic Tumor therapy

DNAzyme-based gene regulation shows great potential for the therapy of many cancers. However, ineffective delivery and insufficient cofactor supply pose challenges for potent gene therapy. In this study, we constructed a smart metal-polyphenol-DNAzyme nanoplatform (TA-Mn@Dz NPs) with intrinsic stability, effective delivery, and cofactor self-supply ability for gene-chemodynamic synergistic tumor therapy. Tannic acid, a plant-derived polyphenol, acts as an intermediate structural unit to mediate the assembly of Mn²⁺/DNAzyme and tumor acid environment-responsive nanocarriers. Intracellularly, the acidic environment triggers the decomposition of TA-Mn@Dz NPs to release DNAzyme and Mn²⁺. The Mn²⁺ ion not only boosts the catalytic cleavage of surviving mRNA for effective gene therapy but also activates chemodynamic therapy (CDT), generating highly toxic ·OH from endogenous H₂O₂. When tail intravenously injected into MCF-7 tumor-bearing mice, the TA-Mn@Dz NPs display desirable synergistic gene-chemodynamic antitumor effects, paving the way for developing DNAzyme-based multifunctional theranostic platforms for biomedical applications. STATEMENT OF SIGNIFICANCE: 1. A smart metal-polyphenol-DNAzyme nanoplatform was constructed for gene-chemodynamic synergistic tumor therapy. 2. Tannic acid act as intermediate structural units to mediate the assembly of Mn²⁺/DNAzyme and tumor acid environment-responsive nanocarriers. 3. The Mn²⁺-ion could not only boost the catalytic cleavage of surviving mRNA for effective gene therapy, but also catalyze endogenous H₂O₂ to form cytotoxic hydroxyl radicals for chemodynamic therapy. 4. Our work paves an extremely simple way to integrate gene therapy with CDT for the dual-catalytic tumor treatment.

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

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

A Practical Approach to Bicyclic Carbamoyl Pyridones with Application to the Synthesis of HIV-1 Integrase Strand Transfer Inhibitors

An efficient one-pot synthetic method has been developed for the preparation of bicyclic carbamoyl pyridones from the known common intermediate methyl 5-((2,4-difluorobenzyl)carbamoyl)-1-(2,2-dimethoxyethyl)-3-methoxy-4-oxo-1,4-dihydropyridine-2-carboxylate (8). The scalable protocol is facile and employs readily available reagents, needing only a single purification as the final step. The utility of the approach was demonstrated by preparing a library of HIV-1 integrase strand transfer inhibitors (INSTIs) that differ by the presence or absence of a double bond in the B-ring of the bicyclic carbamoyl pyridines 6 and 7. Several of the analogs show good antiviral potencies in single-round HIV-1 replication antiviral assays and show no cytotoxicity in cell culture assays. In general, the compounds with a B-ring double bond have higher antiviral potencies than their saturated congeners. Our methodology should be applicable to the synthesis of a range of new metal-chelating analogs.

Structural and functional characterisation of HepTH1-5 peptide as a potential hepcidin replacement

Hepcidin is a principal regulator of iron homeostasis and its dysregulation has been recognised as a causative factor in cancers and iron disorders. The strategy of manipulating the presence of hepcidin peptide has been used for cancer treatment. However, this has demonstrated poor efficiency and has been short-lived in patients. Many studies reported using minihepcidin therapy as an alternative way to treat hepcidin dysregulation, but this was only applied to non-cancer patients. Highly conserved fish hepcidin protein, HepTH1-5, was investigated to determine its potential use in developing a hepcidin replacement for human hepcidin (Hepc25) and as a therapeutic agent by targeting the tumour suppressor protein, p53, through structure-function analysis. The authors found that HepTH1-5 is stably bound to ferroportin, compared to Hepc25, by triggering the ferroportin internalisation via Lys42 and Lys270 ubiquitination, in a similar manner to the Hepc25 activity. Moreover, the residues Ile24 and Gly24, along with copper and zinc ligands, interacted with similar residues, Lys24 and Asp1 of Hepc25, respectively, showing that those molecules are crucial to the hepcidin replacement strategy. HepTH1-5 interacts with p53 and activates its function through phosphorylation. This finding shows that HepTH1-5 might be involved in the apoptosis signalling pathway upon a DNA damage response. This study will be very helpful for understanding the mechanism of the hepcidin replacement and providing insights into the HepTH1-5 peptide as a new target for hepcidin and cancer therapeutics.Communicated by Ramaswamy H. Sarma.

Tellurium-Modified Nucleosides, Nucleotides, and Nucleic Acids with Potential Applications

Tellurium was successfully incorporated into proteins and applied to protein structure determination through X-ray crystallography. However, studies on tellurium modification of DNA and RNA are limited. This review highlights the recent development of Te-modified nucleosides, nucleotides, and nucleic acids, and summarizes the main synthetic approaches for the preparation of 5-PhTe, 2'-MeTe, and 2'-PhTe modifications. Those modifications are compatible with solid-phase synthesis and stable during Te-oligonucleotide purification. Moreover, the ideal electronic and atomic properties of tellurium for generating clear isomorphous signals give Te-modified DNA and RNA great potential applications in 3D crystal structure determination through X-ray diffraction. STM study also shows that Te-modified DNA has strong topographic and current peaks, which immediately suggests potential applications in nucleic acid direct imaging, nanomaterials, molecular electronics, and diagnostics. Theoretical studies indicate the potential application of Te-modified nucleosides in cancer therapy.

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.

RNA Electrostatics: How Ribozymes Engineer Active Sites to Enable Catalysis

Electrostatic interactions are fundamental to RNA structure and function, and intimately influenced by solvation and the ion atmosphere. RNA enzymes, or ribozymes, are catalytic RNAs that are able to enhance reaction rates over a million-fold, despite having only a limited repertoire of building blocks and available set of chemical functional groups. Ribozyme active sites usually occur at junctions where negatively charged helices come together, and in many cases leverage this strained electrostatic environment to recruit metal ions in solution that can assist in catalysis. Similar strategies have been implicated in related artificially engineered DNA enzymes. Herein, we apply Poisson-Boltzmann, 3D-RISM, and molecular simulations to study a set of metal-dependent small self-cleaving ribozymes (hammerhead, pistol, and Varkud satellite) as well as an artificially engineered DNAzyme (8-17) to examine electrostatic features and their relation to the recruitment of monovalent and divalent metal ions important for activity. We examine several fundamental roles for these ions that include: (1) structural integrity of the catalytically active state, (2) pKa tuning of residues involved in acid-base catalysis, and (3) direct electrostatic stabilization of the transition state via Lewis acid catalysis. Taken together, these examples demonstrate how RNA electrostatics orchestrates the site-specific and territorial binding of metal ions to play important roles in catalysis.

Who stole the proton? Suspect general base guanine found with a smoking gun in the pistol ribozyme

The pistol ribozyme (Psr) is one among the most recently discovered classes of small nucleolytic ribozymes that catalyze site-specific RNA self-cleavage through 2'-O-transphosphorylation. The Psr contains a conserved guanine (G40) that in crystal structures is in a position suggesting it plays the role of the general base to abstract a proton from the nucleophile to activate the reaction. Although some functional data is consistent with this mechanistic role, a notable exception is 2-aminopurine (2AP) substitution which has no effect on the rate, unlike similar substitutions across other so-called "G + M" and "G + A" ribozyme classes. Herein we postulate that an alternate conserved guanine, G42, is the primary general base, and provide evidence from molecular simulations that the active site of Psr can undergo local refolding into a structure that is consistent with the common "L-platform/L-scaffold" architecture identified in G + M and G + A ribozyme classes with Psr currently the notable exception. We summarize the key currently available experimental data and present new classical and combined quantum mechanical/molecular mechanical simulation results that collectively suggest a new hypothesis. We hypothesize that there are two available catalytic pathways supported by different conformational states connected by a local refolding of the active site: (1) a primary pathway with an active site architecture aligned with the L-platform/L-scaffold framework where G42 acts as a general base, and (2) a secondary pathway with the crystallographic active site architecture where G40 acts as a general base. We go on to make several experimentally testable predictions, and suggest specific experiments that would ultimately bring closure to the mystery as to "who stole the proton in the pistol ribozyme?".

Combined QM/MM, Machine Learning Path Integral Approach to Compute Free Energy Profiles and Kinetic Isotope Effects in RNA Cleavage Reactions

We present a fast, accurate, and robust approach for determination of free energy profiles and kinetic isotope effects for RNA 2'-O-transphosphorylation reactions with inclusion of nuclear quantum effects. We apply a deep potential range correction (DPRc) for combined quantum mechanical/molecular mechanical (QM/MM) simulations of reactions in the condensed phase. The method uses the second-order density-functional tight-binding method (DFTB2) as a fast, approximate base QM model. The DPRc model modifies the DFTB2 QM interactions and applies short-range corrections to the QM/MM interactions to reproduce ab initio DFT (PBE0/6-31G*) QM/MM energies and forces. The DPRc thus enables both QM and QM/MM interactions to be tuned to high accuracy, and the QM/MM corrections are designed to smoothly vanish at a specified cutoff boundary (6 Å in the present work). The computational speed-up afforded by the QM/MM+DPRc model enables free energy profiles to be calculated that include rigorous long-range QM/MM interactions under periodic boundary conditions and nuclear quantum effects through a path integral approach using a new interface between the AMBER and i-PI software. The approach is demonstrated through the calculation of free energy profiles of a native RNA cleavage model reaction and reactions involving thio-substitutions, which are important experimental probes of the mechanism. The DFTB2+DPRc QM/MM free energy surfaces agree very closely with the PBE0/6-31G* QM/MM results, and it is vastly superior to the DFTB2 QM/MM surfaces with and without weighted thermodynamic perturbation corrections. ¹⁸O and ³⁴S primary kinetic isotope effects are compared, and the influence of nuclear quantum effects on the free energy profiles is examined.

Phosphorescent Ir(III) complexes derived from purine nucleobases

We report the preparation and the study of new types of neutral and cationic phosphorescent heteroleptic Ir(III) complexes derived from 6-phenylpurine nucleosides and nucleotides. Neutral complexes of general formula Ir(C^N)₂(acac) 7, and 8a-c (HC^N = 9-substituted-6-phenyl purine) are orange-red emissive upon photoexcitation, with short lifetimes and good quantum yields (0.42-0.65) in both PMMA films and 2-MeTHF at room temperature. In turn, cationic complexes [Ir(C^N)₂(dtb-bpy)][PF₆] 9, 12a and 12c (dtb-bpy = 4,4'-di-tert-butyl-2,2'-dipyridine) are yellow-green emitters with moderate quantum yields (0.24-0.32).

Not making the cut: Techniques to prevent RNA cleavage in structural studies of RNase-RNA complexes

RNases are varied in the RNA structures and sequences they target for cleavage and are an important type of enzyme in cells. Despite the numerous examples of RNases known, and of those with determined three-dimensional structures, relatively few examples exist with the RNase bound to intact cognate RNA substrate prior to cleavage. To better understand RNase structure and sequence specificity for RNA targets, in vitro methods used to assemble these enzyme complexes trapped in a pre-cleaved state have been developed for a number of different RNases. We have surveyed the Protein Data Bank for such structures and in this review detail methodologies that have successfully been used and relate them to the corresponding structures. We also offer ideas and suggestions for future method development. Many strategies within this review can be used in combination with X-ray crystallography, as well as cryo-EM, and other structure-solving techniques. Our hope is that this review will be used as a guide to resolve future yet-to-be-determined RNase-substrate complex structures.

A DNA Nanocomplex Containing Cascade DNAzymes and Promoter-Like Zn-Mn-Ferrite for Combined Gene/Chemo-dynamic Therapy

Sequential control of exogenous chemical events inside cells is a promising way to regulate cell functions and fate. Herein we report a DNA nanocomplex containing cascade DNAzymes and promoter-like Zn-Mn-Ferrite (ZMF), achieving combined gene/chemo-dynamic therapy. The promoter-like ZMF decomposed in response to intratumoral glutathione to release a sufficient quantity of metal ions, thus promoting cascade DNA/RNA cleavage and free radical generation. Two kinds of DNAzymes were designed for sequential cascade enzymatic reaction, in which metal ions functioned as cofactors. The primary DNAzyme self-cleaved the DNA chain with Zn²⁺ as cofactor, and produced the secondary DNAzyme; the secondary DNAzyme afterwards cleaved the EGR-1 mRNA, and thus downregulated the expression of target EGR-1 protein, achieving DNAzyme-based gene therapy. Meanwhile, the released Zn²⁺ , Mn²⁺ and Fe²⁺ induced Fenton/Fenton-like reactions, during which free radicals were catalytically generated and efficient chemo-dynamic therapy was achieved. In a breast cancer mouse model, the administration of DNA nanocomplex led to a significant therapeutic efficacy of tumor growth suppression.

Bioinspired Artificial Photosynthetic Systems

Mimicking photosynthesis using artificial systems, as a means for solar energy conversion and green fuel generation, is one of the holy grails of modern science. This perspective presents recent advances towards developing artificial photosynthetic systems. In one approach, native photosystems are interfaced with electrodes to yield photobioelectrochemical cells that transform light energy into electrical power. This is exemplified by interfacing photosystem I (PSI) and photosystem II (PSII) as an electrically contacted assembly mimicking the native Z-scheme, and by the assembly of an electrically wired PSI/glucose oxidase biocatalytic conjugate on an electrode support. Illumination of the functionalized electrodes led to light-induced generation of electrical power, or to the generation of photocurrents using glucose as the fuel. The second approach introduces supramolecular photosensitizer nucleic acid/electron acceptor complexes as functional modules for effective photoinduced electron transfer stimulating the subsequent biocatalyzed generation of NADPH or the Pt-nanoparticle-catalyzed evolution of molecular hydrogen. Application of the DNA machineries for scaling-up the photosystems is demonstrated. A third approach presents the integration of artificial photosynthetic modules into dynamic nucleic acid networks undergoing reversible reconfiguration or dissipative transient operation in the presence of auxiliary triggers. Control over photoinduced electron transfer reactions and photosynthetic transformations by means of the dynamic networks is demonstrated.

Zn2+ -Coordination-Driven RNA Assembly with Retained Integrity and Biological Functions

Metal-coordination-directed biomolecule crosslinking in nature has been used for synthesizing various biopolymers, including DNA, peptides, proteins, and polysaccharides. However, the RNA biopolymer has been avoided so far, as due to the poor stability of the RNA molecules, the formation of a biopolymer may alter the biological function of the molecules. Herein, for the first time, we report Zn²⁺ -driven RNA self-assembly forming spherical nanoparticles while retaining the integrity and biological function of RNA. Various functional RNAs of different compositions, shapes, and lengths from 20 to nearly 1000 nucleotides were used, highlighting the versatility of this approach. The assembled nanospheres possess a superior RNA-loading efficiency, pharmacokinetics, and bioavailability. In-vitro and in-vivo evaluation demonstrated mRNA delivery for expressing GFP proteins, and microRNA delivery to triple-negative breast cancer. This coordination-directed self-assembly behavior amplifies the horizons of RNA coordination chemistry and the application scope of RNA-based therapeutics.

Biosensing with DNAzymes

This article provides a comprehensive review of biosensing with DNAzymes, providing an overview of different sensing applications while highlighting major progress and seminal contributions to the field of portable biosensor devices and point-of-care diagnostics. Specifically, the field of functional nucleic acids is introduced, with a specific focus on DNAzymes. The incorporation of DNAzymes into bioassays is then described, followed by a detailed overview of recent advances in the development of in vivo sensing platforms and portable sensors incorporating DNAzymes for molecular recognition. Finally, a critical perspective on the field, and a summary of where DNAzyme-based devices may make the biggest impact are provided.

Protein cofactors and substrate influence Mg2+-dependent structural changes in the catalytic RNA of archaeal RNase P

The ribonucleoprotein (RNP) form of archaeal RNase P comprises one catalytic RNA and five protein cofactors. To catalyze Mg²⁺-dependent cleavage of the 5′ leader from pre-tRNAs, the catalytic (C) and specificity (S) domains of the RNase P RNA (RPR) cooperate to recognize different parts of the pre-tRNA. While ∼250–500 mM Mg²⁺ renders the archaeal RPR active without RNase P proteins (RPPs), addition of all RPPs lowers the Mg²⁺ requirement to ∼10–20 mM and improves the rate and fidelity of cleavage. To understand the Mg²⁺- and RPP-dependent structural changes that increase activity, we used pre-tRNA cleavage and ensemble FRET assays to characterize inter-domain interactions in Pyrococcus furiosus (Pfu) RPR, either alone or with RPPs ± pre-tRNA. Following splint ligation to doubly label the RPR (Cy3-RPR_{C domain} and Cy5-RPR_{S domain}), we used native mass spectrometry to verify the final product. We found that FRET correlates closely with activity, the Pfu RPR and RNase P holoenzyme (RPR + 5 RPPs) traverse different Mg²⁺-dependent paths to converge on similar functional states, and binding of the pre-tRNA by the holoenzyme influences Mg²⁺ cooperativity. Our findings highlight how Mg²⁺ and proteins in multi-subunit RNPs together favor RNA conformations in a dynamic ensemble for functional gains.

Enhancement of the Catalytic Activity of Hammerhead Ribozymes by Organic Cations

Catalytic turnover is important for the application of ribozymes to biotechnology. However, the turnover is often impaired because of the intrinsic high stability of base pairs with cleaved RNA products. Here, organic cations were used as additives to improve the catalytic performance of hammerhead ribozyme constructs that exhibit different kinetic behaviors. Kinetic analysis of substrate cleavage demonstrated that bulky cations, specifically tetra-substituted ammonium ions containing pentyl groups or a benzyl group, have the ability to greatly increase the turnover rate of the ribozymes. Thermal stability analysis of RNA structures revealed that the bulky cations promote the dissociation of cleaved products and refolding of incorrectly folded structures with small disruption of the catalytic structure. The use of bulky cations is a convenient method for enhancing the catalytic activity of hammerhead ribozymes, and the approach may be useful for advancing ribozyme technologies.

A Synergistic DNA-polydopamine-MnO2 Nanocomplex for Near-Infrared-Light-Powered DNAzyme-Mediated Gene Therapy

DNAzyme is emerging for gene therapy. The administration of the in vivo catalytic activity of DNAzyme has proven important but challenging for clinical applications. Herein, we report a synergistic DNA-polydopamine-MnO₂ nanocomplex, which enables near-infrared (NIR)-light-powered catalytic activity of DNAzyme in vivo. The nanocomplex has a hierarchical structure: a DNA nanoframework as the scaffold and polydopamine-MnO₂ (PM) as the coating layer. The DNA nanoframework contains repeated DNAzyme sequences. PM assembles on the surface of the DNA nanoframework. When the nanocomplex accumulates at tumor sites, upon NIR-light radiation, polydopamine induces a temperature elevation at tumor sites via photothermal conversion; meanwhile, glutathione triggers decomposition of PM to release Mn²⁺ to activate DNAzyme in the cytoplasm for gene regulation. In vitro and in vivo experiments show that the PM-induced temperature elevation enhances the Egr-1 mRNA cleavage activity of DNAzyme, promoting downregulation of the Egr-1 protein in tumor cells. In addition, the temperature elevation induces heat stress, achieving a synergistic tumor ablation effect.

Probing Metal-Dependent Phosphate Binding for the Catalysis of the 17E DNAzyme

The RNA-cleaving 17E DNAzyme exhibits different levels of cleavage activity in the presence of various divalent metal ions, with Pb²⁺ giving the fastest cleavage. In this study, the metal-phosphate interaction is probed to understand the trend of activity with different metal ions. For the first-row transition metals, the lowest activity shown by Ni²⁺ correlates with the inhibition by the inorganic phosphate and its water ligand exchange rate, suggesting inner-sphere metal coordination. Cleavage activity with the two stereoisomers of the phosphorothioate-modified substrates, R_p and S_p, indicated that Mg²⁺, Mn²⁺, Fe²⁺, and Co²⁺ had the highest S_p:R_p activity ratio of >900. Comparatively, the activity was much less affected using the thiophilic metals, including Pb²⁺, suggesting inner-sphere coordination. The pH-rate profiles showed that Pb²⁺ was different than the rest of the metal ions in having a smaller slope and a similar fitted apparent pK_a and the pK_a of metal-bound water. Combining previous reports and our current results, we propose that Pb²⁺ most likely plays the role of a general acid while the other metal ions are Lewis acid catalysts interacting with the scissile phosphate.

Developing and Assessing Nonbonded Dummy Models of Magnesium Ion with Different Hydration Free Energy References

A large diversity in the targeted hydration free energies (HFEs) during model parameterization of metal ions was reported in the literature with a difference by dozens of kcal/mol. Here, we developed a series of nonbonded dummy models of the Mg²⁺ ion targeting different HFE references in TIP3P water, followed by assessments of the designed models in the simulations of MgCl₂ solution and biological systems. Together with the comparison of existing models, we conclude that the difference in the targeted HFEs has a limited influence on the model performance, while the usability of these models differs from case to case. The feasibility of reproducing more properties of Mg²⁺ such as diffusion constants and water exchange rates using a nonbonded dummy model is demonstrated. Underestimated activity derivative and osmotic coefficient of MgCl₂ solutions in high concentration reveal a necessity for further optimization of ion-pair interactions. The developed dummy models are applicable to metal coordination with Asp, Glu, and His residues in metalloenzymes, and the performance in predicting monodentate or bidentate binding modes of Asp/Glu residues depends on the complexity of metal centers and the choice of protein force fields. When both the binding modes coexist, the nonbonded dummy models outperform point charge models, probably in need of considering polarization of metal-binding residues by, for instance, charge calibration in classical force fields. This work is valuable for the use and further development of magnesium ion models for simulations of metal-containing systems with good accuracy.

Biochemical analysis of cleavage and ligation activities of the pistol ribozyme from Paenibacillus polymyxa

Nine distinct classes of self-cleaving ribozymes are known to date, of which the pistol ribozyme class was discovered only 5 years ago. Self-cleaving ribozymes are able to cleave their own phosphodiester backbone at a specific site with rates much higher than those of spontaneous RNA degradation. Our study focuses on a bioinformatically predicted pistol ribozyme from the bacterium Paenibacillus polymyxa. We provide a biochemical characterization of this ribozyme, which includes an investigation of the effect of various metal ions on ribozyme cleavage and a kinetic analysis of ribozyme activity under increasing Mg²⁺ concentrations and pH. Based on the obtained results, we discuss a possible catalytic role of divalent metal ions. Moreover, we investigated the ligation activity of the P. polymyxa pistol ribozyme - an aspect that has not been previously analysed for this ribozyme class. We determined that the P. polymyxa pistol ribozyme is almost fully cleaved at equilibrium with the ligation rate constant being nearly 30-fold lower than the cleavage rate constant. In summary, we have characterized an additional representative of this recently discovered ribozyme class isolated from P. polymyxa. We expect that our biochemical characterization of a pistol representative in a cultivatable, genetically tractable organism will support our future investigation of the biological roles of this ribozyme class in bacteria.

Effects of Ionic Liquids on Metalloproteins

In the past decade, innovative protein therapies and bio-similar industries have grown rapidly. Additionally, ionic liquids (ILs) have been an area of great interest and rapid development in industrial processes over a similar timeline. Therefore, there is a pressing need to understand the structure and function of proteins in novel environments with ILs. Understanding the short-term and long-term stability of protein molecules in IL formulations will be key to using ILs for protein technologies. Similarly, ILs have been investigated as part of therapeutic delivery systems and implicated in numerous studies in which ILs impact the activity and/or stability of protein molecules. Notably, many of the proteins used in industrial applications are involved in redox chemistry, and thus often contain metal ions or metal-associated cofactors. In this review article, we focus on the current understanding of protein structure-function relationship in the presence of ILs, specifically focusing on the effect of ILs on metal containing proteins.

A new Pb2+-specific DNAzyme by revisiting the catalytic core of 10-23 DNAzyme

10-23 DNAzyme is a catalytic DNA molecule from in vitro selection, the 15-mer catalytic core was investigated for more DNAzyme variants by block deletions. DNAzyme DZM01 was selected with metal ion dependence of Pb²⁺ ≫ Mn²⁺, with no activity in the presence of Mg²⁺ (20 mM), Ca²⁺ (20 mM), Zn²⁺ (20 mM, pH 6). The unique binding properties of Pb²⁺ with nucleic acids might be responsible for the formation of the catalytic core, which is different from that of other divalent metal ions. More DNAzyme variants are expected to be derived for specific metal ion dependence by various nucleobase sequences and modifications.

Zn2+ -Dependent DNAzymes: From Solution Chemistry to Analytical, Materials and Therapeutic Applications

Since 1994, deoxyribozymes or DNAzymes have been in vitro selected to catalyze various types of reactions. Metal ions play a critical role in DNAzyme catalysis, and Zn²⁺ is a very important one among them. Zn²⁺ has good biocompatibility and can be used for intracellular applications. Chemically, Zn²⁺ is a Lewis acid and it can bind to both the phosphate backbone and the nucleobases of DNA. Zn²⁺ undergoes hydrolysis even at neutral pH, and the partially hydrolyzed polynuclear complexes can affect the interactions with DNA. These features have made Zn²⁺ a unique cofactor for DNAzyme reactions. This review summarizes Zn²⁺ -dependent DNAzymes with an emphasis on RNA-/DNA-cleaving reactions. A key feature is the sharp Zn²⁺ concentration and pH-dependent activity for many of the DNAzymes. The applications of these DNAzymes as biosensors for Zn²⁺ , as therapeutic agents to cleave intracellular RNA, and as chemical biology tools to manipulate DNA are discussed. Future studies can focus on the selection of new DNAzymes with improved performance and detailed biochemical characterizations to understand the role of Zn²⁺ , which can facilitate practical applications of Zn²⁺ -dependent DNAzymes.

Cutting in-line with iron: ribosomal function and non-oxidative RNA cleavage

Divalent metal cations are essential to the structure and function of the ribosome. Previous characterizations of the ribosome performed under standard laboratory conditions have implicated Mg²⁺ as a primary mediator of ribosomal structure and function. Possible contributions of Fe²⁺ as a ribosomal cofactor have been largely overlooked, despite the ribosome's early evolution in a high Fe²⁺ environment, and the continued use of Fe²⁺ by obligate anaerobes inhabiting high Fe²⁺ niches. Here, we show that (i) Fe²⁺ cleaves RNA by in-line cleavage, a non-oxidative mechanism that has not previously been shown experimentally for this metal, (ii) the first-order in-line rate constant with respect to divalent cations is >200 times greater with Fe²⁺ than with Mg²⁺, (iii) functional ribosomes are associated with Fe²⁺ after purification from cells grown under low O₂ and high Fe²⁺ and (iv) a small fraction of Fe²⁺ that is associated with the ribosome is not exchangeable with surrounding divalent cations, presumably because those ions are tightly coordinated by rRNA and deeply buried in the ribosome. In total, these results expand the ancient role of iron in biochemistry and highlight a possible new mechanism of iron toxicity.

In vitro Selection of Chemically Modified DNAzymes

DNAzymes are in vitro selected DNA oligonucleotides with catalytic activities. RNA cleavage is one of the most extensively studied DNAzyme reactions. To expand the chemical functionality of DNA, various chemical modifications have been made during and after selection. In this review, we summarize examples of RNA-cleaving DNAzymes and focus on those modifications introduced during in vitro selection. By incorporating various modified nucleotides via polymerase chain reaction (PCR) or primer extension, a few DNAzymes were obtained that can be specifically activated by metal ions such as Zn2+ and Hg2+. In addition, some modifications were introduced to mimic RNase A that can cleave RNA substrates in the absence of divalent metal ions. In addition, single modifications at the fixed regions of DNA libraries, especially at the cleavage junctions, have been tested, and examples of DNAzymes with phosphorothioate and histidine-glycine modified tertiary amine were successfully obtained specific for Cu2+, Cd2+, Zn2+, and Ni2+. Labeling fluorophore/quencher pair right next to the cleavage junction was also used to obtain signaling DNAzymes for detecting various metal ions and cells. Furthermore, we reviewed work on the cleavage of 2'-5' linked RNA and L-RNA substrates. Finally, applications of these modified DNAzymes as biosensors, RNases, and biochemical probes are briefly described with a few future research opportunities outlined at the end.

Occurrence and stability of lone pair-π and OH-π interactions between water and nucleobases in functional RNAs

We identified over 1000 instances of water-nucleobase stacking contacts in a variety of RNA molecules from a non-redundant set of crystal structures with resolution ≤3.0 Å. Such contacts may be of either the lone pair-π (lp-π) or the OH-π type, in nature. The distribution of the distances of the water oxygen from the nucleobase plane peaks at 3.5 Å for A, G and C, and approximately at 3.1-3.2 Å for U. Quantum mechanics (QM) calculations confirm, as expected, that the optimal energy is reached at a shorter distance for the lp-π interaction as compared to the OH-π one (3.0 versus 3.5 Å). The preference of each nucleobase for either type of interaction closely correlates with its electrostatic potential map. Furthermore, QM calculations show that for all the nucleobases a favorable interaction, of either the lp-π or the OH-π type, can be established at virtually any position of the water molecule above the nucleobase skeleton, which is consistent with the uniform projection of the OW atoms over the nucleobases ring we observed in the experimental occurrences. Finally, molecular dynamics simulations of a model system for the characterization of water-nucleobase stacking contacts confirm the stability of these interactions also under dynamic conditions.

Frequency and hydrogen bonding of nucleobase homopairs in small molecule crystals

We used the high resolution and accuracy of the Cambridge Structural Database (CSD) to provide detailed information regarding base pairing interactions of selected nucleobases. We searched for base pairs in which nucleobases interact with each other through two or more hydrogen bonds and form more or less planar structures. The investigated compounds were either free forms or derivatives of adenine, guanine, hypoxanthine, thymine, uracil and cytosine. We divided our findings into categories including types of pairs, protonation patterns and whether they are formed by free bases or substituted ones. We found base pair types that are exclusive to small molecule crystal structures, some that can be found only in RNA containing crystal structures and many that are native to both environments. With a few exceptions, nucleobase protonation generally followed a standard pattern governed by pKa values. The lengths of hydrogen bonds did not depend on whether the nucleobases forming a base pair were charged or not. The reasons why particular nucleobases formed base pairs in a certain way varied significantly.

Insights into Catalytic Gas-Phase Hydrolysis of Organophosphate Chemical Warfare Agents by MOF-Supported Bimetallic Metal-Oxo Clusters

Zirconium-based metal-organic frameworks (Zr-MOFs) have been reported to be efficient catalysts for the hydrolysis of organophosphate chemical warfare agents (CWAs) in buffered solutions. However, for the gas-phase reaction, which is more relevant to the situation in a battlefield gas mask application, the kinetics of Zr-MOF catalysts may be severely hindered by strong product inhibition. To improve the catalytic performance, we computationally screened a series of synthetically accessible Zr-MOF-supported bimetallic metal-oxo clusters in which the metal-oxygen-metal active motif is preserved, aiming to find catalysts that have lower binding affinities to the hydrolysis product. For the promising catalyst Al₂O₂(OH)₂@NU-1000 identified from the screening using density functional theory, we mapped out the full reaction pathway of gas-phase dimethyl p-nitrophenolphosphate (DMNP) hydrolysis and analyzed the free energy profile as well as the turnover frequency (TOF). We found that the catalytic mechanism on the new catalyst is slightly different from the one on NU-1000, which also led to a different TOF-limiting step. Additional factors that can affect the overall catalytic performance in practical application, such as the amount of ambient moisture and the existence of acid gases that may poison the catalyst, have also been evaluated.

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.

Catalytic Nucleic Acids: Biochemistry, Chemical Biology, Biosensors, and Nanotechnology

Since the initial discovery of ribozymes in the early 1980s, catalytic nucleic acids have been used in different areas. Compared with protein enzymes, catalytic nucleic acids are programmable in structure, easy to modify, and more stable especially for DNA. We take a historic view to summarize a few main interdisciplinary areas of research on nucleic acid enzymes that may have broader impacts. Early efforts on ribozymes in the 1980s have broken the notion that all enzymes are proteins, supplying new evidence for the RNA world hypothesis. In 1994, the first catalytic DNA (DNAzyme) was reported. Since 2000, the biosensor applications of DNAzymes have emerged and DNAzymes are particularly useful for detecting metal ions, a challenging task for enzymes and antibodies. Combined with nanotechnology, DNAzymes are key building elements for switches allowing dynamic control of materials assembly. The search for new DNAzymes and ribozymes is facilitated by developments in DNA sequencing and computational algorithms, further broadening our fundamental understanding of their biochemistry.