A minimal RNA substrate with dual fluorescent probes enables rapid kinetics and provides insight into bacterial RNase P active site interactions.

Bacterial ribonuclease P (RNase P) is a tRNA processing endonuclease that occurs primarily as a ribonucleoprotein with a catalytic RNA subunit (P RNA). As one of the first ribozymes discovered, P RNA is a well-studied model system for understanding RNA catalysis and substrate recognition. Extensive structural and biochemical studies have revealed the structure of RNase P bound to precursor tRNA (ptRNA) and product tRNA. These studies also helped to define active site residues and propose the molecular interactions that are involved in substrate binding and catalysis. However, a detailed quantitative model of the reaction cycle that includes the structures of intermediates and the process of positioning active site metal ions for catalysis is lacking. To further this goal, we used a chemically modified minimal RNA duplex substrate (MD1) to establish a kinetic framework for measuring the functional effects of P RNA active site mutations. Substitution of U69, a critical nucleotide involved in active site Mg2+ binding, was found to reduce catalysis >500-fold as expected, but had no measurable effect on ptRNA binding kinetics. In contrast, the same U69 mutations had little effect on catalysis in Ca2+ compared to reactions containing native Mg2+ ions. CryoEM structures and SHAPE mapping suggested increased flexibility of U69 and adjacent nucleotides in Ca2+ compared to Mg2+. These results support a model in which slow catalysis in Ca2+ is due to inability to engage U69. These studies establish a set of experimental tools to analyze RNase P kinetics and mechanism and can be expanded to gain new insights into the assembly of the active RNase P-ptRNA complex.

Rapid Ribonuclease P Kinetics Measured by Stopped-Flow Fluorescence and Fluorescence Anisotropy.

Stopped-flow fluorescence spectroscopy is a highly sensitive method for measuring rapid enzyme kinetics. A wide range of fluorophores can be employed, and fluorescence and fluorescence polarization can be measured. Thus, binding, conformational changes, and catalysis can, in principle, be measured, making it helpful in probing the entire kinetic landscape of a reaction. In this chapter, we use the bacterial RNA processing enzyme ribonuclease P (RNase P) as a model system to illustrate the determination of the kinetic constants for substrate binding and cleavage, thus allowing mechanistic questions regarding the effects of reaction conditions, mutations, or drug binding to be answered.

The specificity landscape of bacterial ribonuclease P.

Developing quantitative models of substrate specificity for RNA processing enzymes is a key step toward understanding their biology and guiding applications in biotechnology and biomedicine. Optimally, models to predict relative rate constants for alternative substrates should integrate an understanding of structures of the enzyme bound to "fast" and "slow" substrates, large datasets of rate constants for alternative substrates, and transcriptomic data identifying in vivo processing sites. Such data are either available or emerging for bacterial ribonucleoprotein RNase P a widespread and essential tRNA 5' processing endonuclease, thus making it a valuable model system for investigating principles of biological specificity. Indeed, the well-established structure and kinetics of bacterial RNase P enabled the development of high throughput measurements of rate constants for tRNA variants and provided the necessary framework for quantitative specificity modeling. Several studies document the importance of conformational changes in the precursor tRNA substrate as well as the RNA and protein subunits of bacterial RNase P during binding, although the functional roles and dynamics are still being resolved. Recently, results from cryo-EM studies of E. coli RNase P with alternative precursor tRNAs are revealing prospective mechanistic relationships between conformational changes and substrate specificity. Yet, extensive uncharted territory remains, including leveraging these advances for drug discovery, achieving a complete accounting of RNase P substrates, and understanding how the cellular context contributes to RNA processing specificity in vivo.

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 pKa 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.

Kinetic analysis of RNA cleavage by coronavirus Nsp15 endonuclease: Evidence for acid-base catalysis and substrate-dependent metal ion activation.

Understanding the functional properties of severe acute respiratory syndrome coronavirus 2 nonstructural proteins is essential for defining their roles in the viral life cycle, developing improved therapeutics and diagnostics, and countering future variants. Coronavirus nonstructural protein Nsp15 is a hexameric U-specific endonuclease whose functions, substrate specificity, mechanism, and dynamics are not fully defined. Previous studies report that Nsp15 requires Mn2+ ions for optimal activity; however, the effects of divalent ions on Nsp15 reaction kinetics have not been investigated in detail. Here, we analyzed the single- and multiple-turnover kinetics for model ssRNA substrates. Our data confirm that divalent ions are dispensable for catalysis and show that Mn2+ activates Nsp15 cleavage of two different ssRNA oligonucleotide substrates but not a dinucleotide. Biphasic kinetics of ssRNA substrates demonstrates that Mn2+ stabilizes alternative enzyme states that have faster substrate cleavage on the enzyme. However, we did not detect Mn2+-induced conformational changes using CD and fluorescence spectroscopy. The pH-rate profiles in the presence and absence of Mn2+ reveal active-site ionizable groups with similar pKas of ca. 4.8 to 5.2. An Rp stereoisomer phosphorothioate modification at the scissile phosphate had minimal effect on catalysis supporting a mechanism involving an anionic transition state. However, the Sp stereoisomer is inactive because of weak binding, consistent with models that position the nonbridging phosphoryl oxygen deep in the active site. Together, these data demonstrate that Nsp15 employs a conventional acid-base catalytic mechanism passing through an anionic transition state, and that divalent ion activation is substrate dependent.

A general and efficient approach to synthesize the phosphoramidites of 5'-18O labeled purine nucleosides.

5'-18O labeled RNA oligos are important probes to investigate the mechanism of 2'-O-transphosphorylation reactions. Here we describe a general and efficient synthetic approach to the phosphoramidite derivatives of 5'-18O labeled nucleosides starting from the corresponding commercially available 5'-O-DMT protected nucleosides. Using this method, we prepared 5'-18O-guanosine phosphoramidite in 8 steps (13.2% overall yield), 5'-18O-adenosine phosphoramidite in 9 steps (10.1% overall yield) and 5'-18O-2'-deoxyguanosine phosphoramidite in 6 steps (12.8% overall yield). These 5'-18O labeled phosphoramidites can be incorporated into RNA oligos by solid phase synthesis for determination of heavy atom isotope effects in RNA 2'-O-transphosphorylation reactions.

Altered Mechanisms for Acid-Catalyzed RNA Cleavage and Isomerization Reactions Models.

RNA strand cleavage by 2'-O-transphosphorylation is catalyzed not only by numerous nucleolytic RNA enzymes (ribozymes) but also by hydroxide or hydronium ions. In experiments, both cleavage of the 5'-linked nucleoside and isomerization between 3',5'- and 2',5'-phosphodiesters occur under acidic conditions, while only the cleavage reaction is observed under basic conditions. An ab initio path-integral approach for simulating kinetic isotope effects is used to reveal the reaction mechanisms for RNA cleavage and isomerization reactions under acidic conditions. Moreover, the proposed mechanisms can also be combined through the experimental pH-rate profiles.

Dissociative Transition State in Hepatitis Delta Virus Ribozyme Catalysis.

Ribonucleases and small nucleolytic ribozymes are both able to catalyze RNA strand cleavage through 2'-O-transphosphorylation, provoking the question of whether protein and RNA enzymes facilitate mechanisms that pass through the same or distinct transition states. Here, we report the primary and secondary 18O kinetic isotope effects for hepatitis delta virus ribozyme catalysis that reveal a dissociative, metaphosphate-like transition state in stark contrast to the late, associative transition states observed for reactions catalyzed by specific base, Zn2+ ions, or ribonuclease A. This new information provides evidence for a discrete ribozyme active site design that modulates the RNA cleavage pathway to pass through an altered transition state.

Structural and mechanistic basis for recognition of alternative tRNA precursor substrates by bacterial ribonuclease P.

Binding of precursor tRNAs (ptRNAs) by bacterial ribonuclease P (RNase P) involves an encounter complex (ES) that isomerizes to a catalytic conformation (ES*). However, the structures of intermediates and the conformational changes that occur during binding are poorly understood. Here, we show that pairing between the 5' leader and 3'RCCA extending the acceptor stem of ptRNA inhibits ES* formation. Cryo-electron microscopy single particle analysis reveals a dynamic enzyme that becomes ordered upon formation of ES* in which extended acceptor stem pairing is unwound. Comparisons of structures with alternative ptRNAs reveals that once unwinding is completed RNase P primarily uses stacking interactions and shape complementarity to accommodate alternative sequences at its cleavage site. Our study reveals active site interactions and conformational changes that drive molecular recognition by RNase P and lays the foundation for understanding how binding interactions are linked to helix unwinding and catalysis.

Atomic-Scale View of Protein-PEG Interactions that Redirect the Thermal Unfolding Pathway of PEGylated Human Galectin-3.

PEGylation is a promising approach to address the central challenge of applying biologics, i.e., lack of protein stability in the demanding environment of the human body. Wider application is hindered by lack of atomic level understanding of protein-PEG interactions, preventing design of conjugates with predicted properties. We deployed an integrative structural and biophysical approach to address this critical challenge with the PEGylated carbohydrate recognition domain of human galectin-3 (Gal3C), a lectin essential for cell adhesion and potential biologic. PEGylation dramatically increased Gal3C thermal stability, forming a stable intermediate and redirecting its unfolding pathway. Structural details revealed by NMR pointed to a potential role of PEG localization facilitated by charged residues. Replacing these residues subtly altered the protein-PEG interface and thermal unfolding behavior, providing insight into rationally designing conjugates while preserving PEGylation benefits.

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 D2O to H2O 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 D2O ratio (n) versus pL versus kn/k0. 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.

An Ontology for Facilitating Discussion of Catalytic Strategies of RNA-Cleaving Enzymes.

A predictive understanding of the mechanisms of RNA cleavage is important for the design of emerging technology built from biological and synthetic molecules that have promise for new biochemical and medicinal applications. Over the past 15 years, RNA cleavage reactions involving 2'-O-transphosphorylation have been discussed using a simplified framework introduced by Breaker that consists of four fundamental catalytic strategies (designated α, ß, γ, and δ) that contribute to rate enhancement. As more detailed mechanistic data emerge, there is need for the framework to evolve and keep pace. We develop an ontology for discussion of strategies of enzymes that catalyze RNA cleavage via 2'-O-transphosphorylation that stratifies Breaker's framework into primary (1°), secondary (2°), and tertiary (3°) contributions to enable more precise interpretation of mechanism in the context of structure and bonding. Further, we point out instances where atomic-level changes give rise to changes in more than one catalytic contribution, a phenomenon we refer to as "functional blurring". We hope that this ontology will help clarify our conversations and pave the path forward toward a consensus view of these fundamental and fascinating mechanisms. The insight gained will deepen our understanding of RNA cleavage reactions catalyzed by natural protein and RNA enzymes, as well as aid in the design of new engineered DNA and synthetic enzymes.

Structure-guided design of anti-cancer ribonucleotide reductase inhibitors.

Ribonucleotide reductase (RR) catalyses the rate-limiting step of dNTP synthesis, establishing it as an important cancer target. While RR is traditionally inhibited by nucleoside-based antimetabolites, we recently discovered a naphthyl salicyl acyl hydrazone-based inhibitor (NSAH) that binds reversibly to the catalytic site (C-site). Here we report the synthesis and in vitro evaluation of 13 distinct compounds (TP1-13) with improved binding to hRR over NSAH (TP8), with lower KD's and more predicted residue interactions. Moreover, TP6 displayed the greatest growth inhibiting effect in the Panc1 pancreatic cancer cell line with an IC50 of 0.393 µM. This represents more than a 2-fold improvement over NSAH, making TP6 the most potent compound against pancreatic cancer emerging from the hydrazone inhibitors. NSAH was optimised by the addition of cyclic and polar groups replacing the naphthyl moiety, which occupies the phosphate-binding pocket in the C-site, establishing a new direction in inhibitor design.

Distributive enzyme binding controlled by local RNA context results in 3' to 5' directional processing of dicistronic tRNA precursors by Escherichia coli ribonuclease P.

RNA processing by ribonucleases and RNA modifying enzymes often involves sequential reactions of the same enzyme on a single precursor transcript. In Escherichia coli, processing of polycistronic tRNA precursors involves separation into individual pre-tRNAs by one of several ribonucleases followed by 5' end maturation by ribonuclease P. A notable exception are valine and lysine tRNAs encoded by three polycistronic precursors that follow a recently discovered pathway involving initial 3' to 5' directional processing by RNase P. Here, we show that the dicistronic precursor containing tRNAvalV and tRNAvalW undergoes accurate and efficient 3' to 5' directional processing by RNase P in vitro. Kinetic analyses reveal a distributive mechanism involving dissociation of the enzyme between the two cleavage steps. Directional processing is maintained despite swapping or duplicating the two tRNAs consistent with inhibition of processing by 3' trailer sequences. Structure-function studies identify a stem-loop in 5' leader of tRNAvalV that inhibits RNase P cleavage and further enforces directional processing. The results demonstrate that directional processing is an intrinsic property of RNase P and show how RNA sequence and structure context can modulate reaction rates in order to direct precursors along specific pathways.

Evidence That Nucleophile Deprotonation Exceeds Bond Formation in the HDV Ribozyme Transition State.

Steric constraints imposed by the active sites of protein and RNA enzymes pose major challenges to the investigation of structure-function relationships within these systems. As a strategy to circumvent such constraints in the HDV ribozyme, we have synthesized phosphoramidites from propanediol derivatives and incorporated them at the 5'-termini of RNA and DNA oligonucleotides to generate a series of novel substrates with nucleophiles perturbed electronically through geminal fluorination. In nonenzymatic, hydroxide-catalyzed intramolecular transphosphorylation of the DNA substrates, pH-rate profiles revealed that fluorine substitution reduces the maximal rate and the kinetic p Ka, consistent with the expected electron-withdrawing effect. In HDV ribozyme reactions, we observed that the RNA substrates undergo transphosphorylation relatively efficiently, suggesting that the conformational constraints imposed by a ribofuranose ring are not strictly required for ribozyme catalysis. In contrast to the nonenzymatic reactions, however, substrate fluorination modestly increases the ribozyme reaction rate, consistent with a mechanism in which (1) the 2'-hydroxyl nucleophile exists predominantly in its neutral, protonated form in the ground state and (2) the 2'-hydroxyl bears some negative charge in the rate-determining step, consistent with a transition state in which the extent of 2'-OH deprotonation exceeds the extent of P-O bond formation.

Structure-Guided Synthesis and Mechanistic Studies Reveal Sweetspots on Naphthyl Salicyl Hydrazone Scaffold as Non-Nucleosidic Competitive, Reversible Inhibitors of Human Ribonucleotide Reductase.

Ribonucleotide reductase (RR), an established cancer target, is usually inhibited by antimetabolites, which display multiple cross-reactive effects. Recently, we discovered a naphthyl salicyl acyl hydrazone-based inhibitor (NSAH or E-3a) of human RR (hRR) binding at the catalytic site (C-site) and inhibiting hRR reversibly. We herein report the synthesis and biochemical characterization of 25 distinct analogs. We designed each analog through docking to the C-site of hRR based on our 2.7 Å X-ray crystal structure (PDB ID: 5TUS). Broad tolerance to minor structural variations preserving inhibitory potency is observed. E-3f (82% yield) displayed an in vitro IC50 of 5.3 ± 1.8 µM against hRR, making it the most potent in this series. Kinetic assays reveal that E-3a, E-3c, E-3t, and E-3w bind and inhibit hRR through a reversible and competitive mode. Target selectivity toward the R1 subunit of hRR is established, providing a novel way of inhibition of this crucial enzyme.

Kinetic Isotope Effect Analysis of RNA 2'-O-Transphosphorylation.

The breaking of RNA strands by 2'-O-transphosphorylation is a ubiquitous reaction in biology, and enzymes that catalyze this reaction play key roles in RNA metabolism. The mechanisms of 2'-O-transphosphorylation in solution are relatively well studied, but complex and can involve different transition states depending on how the reaction is catalyzed. Because of this complexity and the lack of experimental information on transition-state structure, pinning down the chemical details of enzyme-catalyzed RNA strand cleavage has been difficult. Kinetic isotope effects (KIEs) provide information about changes in bonding as a reaction proceeds from ground state to transition state, and therefore they provide a powerful tool for revealing mechanistic detail. Application of kinetic isotope analyses to RNA 2'-O-transphosphorylation faces three fundamental challenges: synthesis of RNA substrate isotopomers with 18O substitutions at the 2'-O, 5'-O and nonbridging phosphoryl oxygens; determination of the 18O/16O ratios in the residual unreacted substrate or product RNAs; and analyzing these data to allow calculation of the KIEs for use in evaluating different mechanistic scenarios. In this chapter, we outline methods for surmounting these challenges for solution RNA 2'-O-transphosphorylation reactions, and we describe their initial application to understand nonenzymatic solution reactions and reactions catalyzed by the enzyme ribonuclease A.

Phylogenetic sequence analysis and functional studies reveal compensatory amino acid substitutions in loop 2 of human ribonucleotide reductase.

Eukaryotic class I ribonucleotide reductases (RRs) generate deoxyribonucleotides for DNA synthesis. Binding of dNTP effectors is coupled to the formation of active dimers and induces conformational changes in a short loop (loop 2) to regulate RR specificity among its nucleoside diphosphate substrates. Moreover, ATP and dATP bind at an additional allosteric site 40 Å away from loop 2 and thereby drive formation of activated or inactive hexamers, respectively. To better understand how dNTP binding influences specificity, activity, and oligomerization of human RR, we aligned >300 eukaryotic RR sequences to examine natural sequence variation in loop 2. We found that most amino acids in eukaryotic loop 2 were nearly invariant in this sample; however, two positions co-varied as nonconservative substitutions (N291G and P294K; human numbering). We also found that the individual N291G and P294K substitutions in human RR additively affect substrate specificity. The P294K substitution significantly impaired effector-induced oligomerization required for enzyme activity, and oligomerization was rescued in the N291G/P294K enzyme. None of the other mutants exhibited altered ATP-mediated hexamerization; however, certain combinations of loop 2 mutations and dNTP effectors perturbed ATP's role as an allosteric activator. Our results demonstrate that the observed compensatory covariation of amino acids in eukaryotic loop 2 is essential for its role in dNTP-induced dimerization. In contrast, defects in substrate specificity are not rescued in the double mutant, implying that functional sequence variation elsewhere in the protein is necessary. These findings yield insight into loop 2's roles in regulating RR specificity, allostery, and oligomerization.

Potent competitive inhibition of human ribonucleotide reductase by a nonnucleoside small molecule.

Human ribonucleotide reductase (hRR) is crucial for DNA replication and maintenance of a balanced dNTP pool, and is an established cancer target. Nucleoside analogs such as gemcitabine diphosphate and clofarabine nucleotides target the large subunit (hRRM1) of hRR. These drugs have a poor therapeutic index due to toxicity caused by additional effects, including DNA chain termination. The discovery of nonnucleoside, reversible, small-molecule inhibitors with greater specificity against hRRM1 is a key step in the development of more effective treatments for cancer. Here, we report the identification and characterization of a unique nonnucleoside small-molecule hRR inhibitor, naphthyl salicylic acyl hydrazone (NSAH), using virtual screening, binding affinity, inhibition, and cell toxicity assays. NSAH binds to hRRM1 with an apparent dissociation constant of 37 µM, and steady-state kinetics reveal a competitive mode of inhibition. A 2.66-Å resolution crystal structure of NSAH in complex with hRRM1 demonstrates that NSAH functions by binding at the catalytic site (C-site) where it makes both common and unique contacts with the enzyme compared with NDP substrates. Importantly, the IC50 for NSAH is within twofold of gemcitabine for growth inhibition of multiple cancer cell lines, while demonstrating little cytotoxicity against normal mobilized peripheral blood progenitor cells. NSAH depresses dGTP and dATP levels in the dNTP pool causing S-phase arrest, providing evidence for RR inhibition in cells. This report of a nonnucleoside reversible inhibitor binding at the catalytic site of hRRM1 provides a starting point for the design of a unique class of hRR inhibitors.