Catalysis by ribonuclease A is limited by the rate of substrate association

Controlled Level of Contamination Coupled to Deep Sequencing (CoLoC-seq) Probes the Global Localisation Topology of Organelle Transcriptomes

Information on RNA localisation is essential for understanding physiological and pathological processes, such as gene expression, cell reprogramming, host-pathogen interactions, and signalling pathways involving RNA transactions at the level of membrane-less or membrane-bounded organelles and extracellular vesicles. In many cases, it is important to assess the topology of RNA localisation, i.e., to distinguish the transcripts encapsulated within an organelle of interest from those merely attached to its surface. This allows establishing which RNAs can, in principle, engage in local molecular interactions and which are prevented from interacting by membranes or other physical barriers. The most widely used techniques interrogating RNA localisation topology are based on the treatment of isolated organelles with RNases with subsequent identification of the surviving transcripts by northern blotting, qRT-PCR, or RNA-seq. However, this approach produces incoherent results and many false positives. Here, we describe Controlled Level of Contamination coupled to deep sequencing (CoLoC-seq), a more refined subcellular transcriptomics approach that overcomes these pitfalls. CoLoC-seq starts by the purification of organelles of interest. They are then either left intact or lysed and subjected to a gradient of RNase concentrations to produce unique RNA degradation dynamics profiles, which can be monitored by northern blotting or RNA-seq. Through straightforward mathematical modelling, CoLoC-seq distinguishes true membrane-enveloped transcripts from degradable and non-degradable contaminants of any abundance. The method has been implemented in the mitochondria of HEK293 cells, where it outperformed alternative subcellular transcriptomics approaches. It is applicable to other membrane-bounded organelles, e.g., plastids, single-membrane organelles of the vesicular system, extracellular vesicles, or viral particles. Key features • Tested on human mitochondria; potentially applicable to cell cultures, non-model organisms, extracellular vesicles, enveloped viruses, tissues; does not require genetic manipulations or highly pure organelles. • In the case of human cells, the required amount of starting material is ~2,500 cm2 of 80% confluent cells (or ~3 × 108 HEK293 cells). • CoLoC-seq implements a special RNA-seq strategy to selectively capture intact transcripts, which requires RNases generating 5'-hydroxyl and 2'/3'-phosphate termini (e.g., RNase A, RNase I). • Relies on nonlinear regression software with customisable exponential functions.

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 SARS-CoV-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 have not been fully defined. Previous studies report SARS-CoV-2 Nsp15 requires Mn²⁺ 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 single-stranded RNA substrates. Our data confirm that divalent ions are dispensable for catalysis and show that Mn²⁺ activates Nsp15 cleavage of two different ssRNA oligonucleotide substrates, but not a dinucleotide. Furthermore, biphasic kinetics of ssRNA substrates demonstrates that Mn²⁺ stabilizes alternative enzyme states that have faster substrate cleavage on the enzyme. However, we did not detect Mn²⁺-induced conformational changes using CD and fluorescence spectroscopy. The pH-rate profiles in the presence and absence of Mn²⁺ are consistent with active site ionizable groups with similar pK_as of ca. 4.8-5.2. We found the Rp stereoisomer phosphorothioate modification at the scissile phosphate had minimal effect on catalysis, which supports a mechanism involving an anionic transition state. In contrast, the Sp stereoisomer is inactive due to weak binding, consistent with models that position the non-bridging phosphoryl oxygen deep in the active site. Together, these kinetic 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.

CoLoC-seq probes the global topology of organelle transcriptomes

Proper RNA localisation is essential for physiological gene expression. Various kinds of genome-wide approaches permit to comprehensively profile subcellular transcriptomes. Among them, cell fractionation methods, that couple RNase treatment of isolated organelles to the sequencing of protected transcripts, remain most widely used, mainly because they do not require genetic modification of the studied system and can be easily implemented in any cells or tissues, including in non-model species. However, they suffer from numerous false-positives since incompletely digested contaminant RNAs can still be captured and erroneously identified as resident transcripts. Here we introduce Controlled Level of Contamination coupled to deep sequencing (CoLoC-seq) as a new subcellular transcriptomics approach that efficiently bypasses this caveat. CoLoC-seq leverages classical enzymatic kinetics and tracks the depletion dynamics of transcripts in a gradient of an exogenously added RNase, with or without organellar membranes. By means of straightforward mathematical modelling, CoLoC-seq infers the localisation topology of RNAs and robustly distinguishes between genuinely resident, luminal transcripts and merely abundant surface-attached contaminants. Our generic approach performed well on human mitochondria and is in principle applicable to other membrane-bounded organelles, including plastids, compartments of the vacuolar system, extracellular vesicles, and viral particles.

A method for quantifying how the activity of an enzyme is affected by the net charge of its nearest crowded neighbor

The electrostatic effects of protein crowding have not been systematically explored. Rather, protein crowding is generally studied with co‐solvents or crowders that are electrostatically neutral, with no methods to measure how the net charge (Z) of a crowder affects protein function. For example, can the activity of an enzyme be affected electrostatically by the net charge of its neighbor in crowded milieu? This paper reports a method for crowding proteins of different net charge to an enzyme via semi‐random chemical crosslinking. As a proof of concept, RNase A was crowded (at distances ≤ the Debye length) via crosslinking to different heme proteins with Z = +8.50 ± 0.04, Z = +6.39 ± 0.12, or Z = −10.30 ± 1.32. Crosslinking did not disrupt the structure of proteins, according to amide H/D exchange, and did not inhibit RNase A activity. For RNase A, we found that the electrostatic environment of each crowded neighbor had significant effects on rates of RNA hydrolysis. Crowding with cationic cytochrome c led to increases in activity, while crowding with anionic “supercharged” cytochrome c or myoglobin diminished activity. Surprisingly, electrostatic crowding effects were amplified at high ionic strength (I = 0.201 M) and attenuated at low ionic strength (I = 0.011 M). This salt dependence might be caused by a unique set of electric double layers at the dimer interspace (maximum distance of 8 Å, which cannot accommodate four layers). This new method of crowding via crosslinking can be used to search for electrostatic effects in protein crowding.

Hofmeister Effects of Group II Cations as Seen in the Unfolding of Ribonuclease A

This work studies the effects of alkaline-earth cation addition upon the unfolding free energy of a model protein, pancreatic Ribonuclease A (RNase A) by DSC analysis.  RNase A was chosen because it: a) does not specifically bind Mg 2+ , Ca 2+ and Sr 2+ cations and b) maintains its structural integrity throughout a large pH range. We have measured and compared the effects of NaCl, MgCl 2 , CaCl 2 and SrCl 2 addition on the melting point of RNase A.  Our results show that even though the addition of group II cations to aqueous solvent reduces the solubility of nonpolar residues (and enhances the hydrophobic effect), their interactions with the amide moieties are strong enough to "salt-them-in" the solvent, thereby causing an overall reduction in protein stability. We demonstrate that amide-cation interactions are a major contributor to the observed "Hofmeister Effects" of group II cations in protein folding. Our analysis suggests that protein folding "Hofmeister Effects" of group II cations, are mostly the aggregate sum of how cation addition simultaneously salts-out hydrophobic moieties through increasing the cavitation free energy, while promoting the salting-in of amide moieties through contact pair formation.

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.

Cleaning Up Mechanistic Debris Generated by Twister Ribozymes Using Computational RNA Enzymology

The catalytic properties of RNA have been a subject of fascination and intense research since their discovery over 30 years ago. Very recently, several classes of nucleolytic ribozymes have emerged and been characterized structurally. Among these, the twister ribozyme has been center-stage, and a topic of debate about its architecture and mechanism owing to conflicting interpretations of different crystal structures, and in some cases conflicting interpretations of the same functional data. In the present work, we attempt to clean up the mechanistic "debris" generated by twister ribozymes using a comprehensive computational RNA enzymology approach aimed to provide a unified interpretation of existing structural and functional data. Simulations in the crystalline environment and in solution provide insight into the origins of observed differences in crystal structures, and coalesce on a common active site architecture, and dynamical ensemble in solution. We use GPU-accelerated free energy methods with enhanced sampling to ascertain microscopic nucleobase pK a values of the implicated general acid and base, from which predicted activity-pH profiles can be compared directly with experiments. Next, ab initio quantum mechanical/molecular mechanical (QM/MM) simulations with full dynamic solvation under periodic boundary conditions are used to determine mechanistic pathways through multi-dimensional free energy landscapes for the reaction. We then characterize the rate-controlling transition state, and make predictions about kinetic isotope effects and linear free energy relations. Computational mutagenesis is performed to explain the origin of rate effects caused by chemical modifications and make experimentally testable predictions. Finally, we provide evidence that helps to resolve conflicting issues related to the role of metal ions in catalysis. Throughout each stage, we highlight how a conserved L-platform structural motif, to- gether with a key L-anchor residue, forms the characteristic active site scaffold enabling each of the catalytic strategies to come together not only for the twister ribozyme, but the majority of the known small nucleolytic ribozyme classes.

Do soft anions promote protein denaturation through binding interactions? A case study using ribonuclease A

It has long been known that large soft anions like bromide, iodide and thiocyanate are protein denaturing agents, but their mechanism of action is still unclear. In this work we have investigated the protein denaturing properties of these anions using Ribonuclease A (RNase A) as a model protein system. Salt-induced perturbations to the protein folding free energy were determined using differential scanning calorimetry and the results demonstrate that the addition of sodium iodide and sodium thiocyanate significantly decreases the melting temperature of the protein. In order to account for this reduction in protein stability, we show that the introduction of salts that contain soft anions to the aqueous solvent perturbs the protein unfolding free energy through three mechanisms: (a) screening Coulomb interactions that exist between charged protein residues, (b) Hofmeister effects, and (c) specific anion binding to CH and CH2 moieties in the protein polypeptide backbone. Using the micellization of 1,2-hexanediol as a ruler for hydrophobicity, we have devised a practical methodology that separates the Coulomb and Hofmeister contributions of salts to the protein unfolding free energy. This allowing us to isolate the contribution of soft anion binding interactions to the unfolding process. The analysis shows that binding contributions have the largest magnitude, confirming that it is the binding of soft anions to the polypeptide backbone that is the main promoter of protein unfolding.

Structural and mechanistic basis for preferential deadenylation of U6 snRNA by Usb1

Post-transcriptional modification of snRNA is central to spliceosome function. Usb1 is an exoribonuclease that shortens the oligo-uridine tail of U6 snRNA, resulting in a terminal 2',3' cyclic phosphate group in most eukaryotes, including humans. Loss of function mutations in human Usb1 cause the rare disorder poikiloderma with neutropenia (PN), and result in U6 snRNAs with elongated 3' ends that are aberrantly adenylated. Here, we show that human Usb1 removes 3' adenosines with 20-fold greater efficiency than uridines, which explains the presence of adenylated U6 snRNAs in cells lacking Usb1. We determined three high-resolution co-crystal structures of Usb1: wild-type Usb1 bound to the substrate analog adenosine 5'-monophosphate, and an inactive mutant bound to RNAs with a 3' terminal adenosine and uridine. These structures, along with QM/MM MD simulations of the catalytic mechanism, illuminate the molecular basis for preferential deadenylation of U6 snRNA. The extent of Usb1 processing is influenced by the secondary structure of U6 snRNA.

Self-Transfecting Micellar RNA: Modulating Nanoparticle Cell Interactions via High Density Display of Small Molecule Ligands on Micelle Coronas

The intracellular delivery of synthetic nucleic acids represents a major challenge in biotechnology and in biomedicine. Methods to deliver short, double-stranded RNA to living cells are of particular interest because of the potential to engage the RNA interference machinery and to regulate mRNA expression. In this work, we describe novel RNA-polymer amphiphiles that assemble into spherical micellar nanoparticles with diameters of ca. 15-30 nm and efficiently enter live cells without transfection reagents. Each micelle consists of approximately 100 RNA strands forming a densely packed corona around a polymeric core. Importantly, the surface-displayed RNA remains accessible for hybridization with complementary RNA. Chemical modification of the termini of hybridized RNA strands enabled the display of small organic moieties on the outer surface of the micelle corona. We found that some of these modifications can have a tremendous impact on cellular internalization efficiencies. The display of hydrophobic dabcyl or stilbene units dramatically increased cell uptake, whereas hydrophilic neutral hydroxy or anionic phosphate residues were ineffective. Interestingly, neither of these modifications mediated noticeable uptake of free RNA oligonucleotides. We infer that their high density display on micellar nanoparticle surfaces is key for the observed effect; achieved with local effective surface concentrations in the millimolar range. We speculate that weak interactions with cell surface receptors that are amplified by the multivalent presentation of such modifications may be responsible. The installation of small molecule ligands on nanomaterial surfaces via hybridization of chemically modified oligonucleotides offers a simple and straightforward way to modulate cellular uptake of nanoparticles. Biological functionality of micellar RNA was demonstrated through the sequence-specific regulation of mRNA expression in HeLa cells.

Integration of kinetic isotope effect analyses to elucidate ribonuclease mechanism

The well-studied mechanism of ribonuclease A is believed to involve concerted general acid-base catalysis by two histidine residues, His12 and His119. The basic features of this mechanism are often cited to explain rate enhancement by both protein and RNA enzymes that catalyze RNA 2'-O-transphosphorylation. Recent kinetic isotope effect analyses and computational studies are providing a more chemically detailed description of the mechanism of RNase A and the rate limiting transition state. Overall, the results support an asynchronous mechanism for both solution and ribonuclease catalyzed reactions in which breakdown of a transient dianoinic phosphorane intermediate by 5'OP bond cleavage is rate limiting. Relative to non-enzymatic reactions catalyzed by specific base, a smaller KIE on the 5'O leaving group and a less negative βLG are observed for RNase A catalysis. Quantum mechanical calculations consistent with these data support a model in which electrostatic and H-bonding interactions with the non-bridging oxygens and proton transfer from His119 render departure of the 5'O less advanced and stabilize charge buildup in the transition state. Both experiment and computation indicate advanced 2'OP bond formation in the rate limiting transition state. However, this feature makes it difficult to resolve the chemical steps involved in 2'O activation. Thus, modeling the transition state for RNase A catalysis underscores those elements of its chemical mechanism that are well resolved, as well as highlighting those where ambiguity remains. This article is part of a Special Issue entitled: Enzyme Transition States from Theory and Experiment.

Interpretation of pH-activity profiles for acid-base catalysis from molecular simulations

The measurement of reaction rate as a function of pH provides essential information about mechanism. These rates are sensitive to the pK(a) values of amino acids directly involved in catalysis that are often shifted by the enzyme active site environment. Experimentally observed pH-rate profiles are usually interpreted using simple kinetic models that allow estimation of "apparent pK(a)" values of presumed general acid and base catalysts. One of the underlying assumptions in these models is that the protonation states are uncorrelated. In this work, we introduce the use of constant pH molecular dynamics simulations in explicit solvent (CpHMD) with replica exchange in the pH-dimension (pH-REMD) as a tool to aid in the interpretation of pH-activity data of enzymes and to test the validity of different kinetic models. We apply the methods to RNase A, a prototype acid-base catalyst, to predict the macroscopic and microscopic pK(a) values, as well as the shape of the pH-rate profile. Results for apo and cCMP-bound RNase A agree well with available experimental data and suggest that deprotonation of the general acid and protonation of the general base are not strongly coupled in transphosphorylation and hydrolysis steps. Stronger coupling, however, is predicted for the Lys41 and His119 protonation states in apo RNase A, leading to the requirement for a microscopic kinetic model. This type of analysis may be important for other catalytic systems where the active forms of the implicated general acid and base are oppositely charged and more highly correlated. These results suggest a new way for CpHMD/pH-REMD simulations to bridge the gap with experiments to provide a molecular-level interpretation of pH-activity data in studies of enzyme mechanisms.

DNA crystals as vehicles for biocatalysis

Here we demonstrate that protein enzymes captured in the solvent channels of three-dimensional DNA crystals are catalytically active. Using RNase A as a model enzyme system, we show that crystals infused with enzyme can cleave a dinucleotide substrate with similar kinetic restrictions as other immobilized enzyme systems. This new vehicle for immobilized enzymes, created entirely from biomolecules, opens possibilities for developing modular solid-state catalysts that could be both biocompatible and biodegradable.

A simple assay for the ribonuclease activity of ribonucleases in the presence of ethidium bromide

The ribonuclease (RNase) activity of ribonucleases has been assayed by observing the change in fluorescence intensity of ethidium bromide on binding with yeast RNA. The binding of EtBr with RNA was monitored via UV-vis and fluorimetric methods. The degradation of RNA by RNase A was monitored by the change in fluorescence emission intensity of ethidium bromide at 600nm on excitation at 510nm. The ribonucleolytic activity of RNase A and angiogenin at various pH values was determined by this method. From this technique we have also determined the macroscopic pKa values of active site residues of these enzymes. This assay permits the evaluation of the catalytic efficiency of enzymatic proteins ranging from high ribonucleolytic activity to low ribonucleolytic activity toward the natural substrate RNA.

A measure of the broad substrate specificity of enzymes based on 'duplicate' catalytic residues

The ability of an enzyme to select and act upon a specific class of compounds with unerring precision and efficiency is an essential feature of life. Simultaneously, these enzymes often catalyze the reaction of a range of similar substrates of the same class, and also have promiscuous activities on unrelated substrates. Previously, we have established a methodology to quantify promiscuous activities in a wide range of proteins. In the current work, we quantitatively characterize the active site for the ability to catalyze distinct, yet related, substrates (BRASS). A protein with known structure and active site residues provides the framework for computing ‘duplicate’ residues, each of which results in slightly modified replicas of the active site scaffold. Such spatial congruence is supplemented by Finite difference Poisson Boltzmann analysis which filters out electrostatically unfavorable configurations. The congruent configurations are used to compute an index (BrassIndex), which reflects the broad substrate profile of the active site. We identify an acetylhydrolase and a methyltransferase as having the lowest and highest BrassIndex, respectively, from a set of non-homologous proteins extracted from the Catalytic Site Atlas. The acetylhydrolase, a regulatory enzyme, is known to be highly specific for platelet-activating factor. In the methyltransferase (PDB: 1QAM), various combinations of glycine (Gly38/40/42), asparagine (Asn101/11) and glutamic acid (Glu59/36) residues having similar spatial and electrostatic profiles with the specified scaffold (Gly38, Asn101 and Glu59) exemplifies the broad substrate profile such an active site may provide. ‘Duplicate’ residues identified by relaxing the spatial and/or electrostatic constraints can be the target of directed evolution methodologies, like saturation mutagenesis, for modulating the substrate specificity of proteins.

A measure of the promiscuity of proteins and characteristics of residues in the vicinity of the catalytic site that regulate promiscuity

Promiscuity, the basis for the evolution of new functions through 'tinkering' of residues in the vicinity of the catalytic site, is yet to be quantitatively defined. We present a computational method Promiscuity Indices Estimator (PROMISE)--based on signatures derived from the spatial and electrostatic properties of the catalytic residues, to estimate the promiscuity (PromIndex) of proteins with known active site residues and 3D structure. PromIndex reflects the number of different active site signatures that have congruent matches in close proximity of its native catalytic site, the quality of the matches and difference in the enzymatic activity. Promiscuity in proteins is observed to follow a lognormal distribution (μ = 0.28, σ = 1.1 reduced chi-square = 3.0E-5). The PROMISE predicted promiscuous functions in any protein can serve as the starting point for directed evolution experiments. PROMISE ranks carboxypeptidase A and ribonuclease A amongst the more promiscuous proteins. We have also investigated the properties of the residues in the vicinity of the catalytic site that regulates its promiscuity. Linear regression establishes a weak correlation (R(2)∼0.1) between certain properties of the residues (charge, polar, etc) in the neighborhood of the catalytic residues and PromIndex. A stronger relationship states that most proteins with high promiscuity have high percentages of charged and polar residues within a radius of 3 Å of the catalytic site, which is validated using one-tailed hypothesis tests (P-values∼0.05). Since it is known that these characteristics are key factors in catalysis, their relationship with the promiscuity index cross validates the methodology of PROMISE.

Reengineering rate-limiting, millisecond enzyme motions by introduction of an unnatural amino acid

Rate-limiting millisecond motions in wild-type (WT) Ribonuclease A (RNase A) are modulated by histidine 48. Here, we incorporate an unnatural amino acid, thia-methylimidazole, at this site (H48C-4MI) to investigate the effects of a single residue on protein motions over multiple timescales and on enzyme catalytic turnover. Molecular dynamics simulations reveal that H48C-4MI retains some crucial WT-like hydrogen bonding interactions but the extent of protein-wide correlated motions in the nanosecond regime is decreased relative to WT. NMR Carr-Purcell-Meiboom-Gill relaxation dispersion experiments demonstrate that millisecond conformational motions in H48C-4MI are present over a similar pH range compared to WT. Furthermore, incorporation of this nonnatural amino acid allows retention of WT-like catalytic activity over the full pH range. These studies demonstrate that the complexity of the protein energy landscape during the catalytic cycle can be maintained using unnatural amino acids, which may prove useful in enzyme design efforts.

Alteration of hydrogen bonding in the vicinity of histidine 48 disrupts millisecond motions in RNase A

The motion of amino acid residues on the millisecond (ms) time scale is involved in the tight regulation of catalytic function in numerous enzyme systems. Using a combination of mutational, enzymological, and relaxation-compensated (15)N Carr-Purcell-Meiboom-Gill (CPMG) methods, we have previously established the conformational significance of the distant His48 residue and the neighboring loop 1 in RNase A function. These studies suggested that RNase A relies on an intricate network of hydrogen bonding interactions involved in propagating functionally relevant, long-range ms motions to the catalytic site of the enzyme. To further investigate the dynamic importance of this H-bonding network, this study focuses on the individual replacement of Thr17 and Thr82 with alanine, effectively altering the key H-bonding interactions that connect loop 1 and His48 to the rest of the protein. (15)N CPMG dispersion studies, nuclear magnetic resonance (NMR) chemical shift analysis, and NMR line shape analysis of point mutants T17A and T82A demonstrate that the evolutionarily conserved single H-bond linking His48 to Thr82 is essential for propagating ms motions from His48 to the active site of RNase A on the time scale of catalytic turnover, whereas the T17A mutation increases the off rate and conformational exchange motions in loop 1. Accumulating evidence from our mutational studies indicates that residues experiencing conformational exchange in RNase A can be grouped into two separate clusters displaying distinct dynamical features, which appear to be independently affected by mutation. Overall, this study illuminates how tightly controlled and finely tuned ms motions are in RNase A, suggesting that designed modulation of protein motions may be possible.

The crystal structure of ribonuclease A in complex with thymidine-3'-monophosphate provides further insight into ligand binding

Thymidine-3'-monophosphate (3'-TMP) is a competitive inhibitor analogue of the 3'-CMP and 3'-UMP natural product inhibitors of bovine pancreatic ribonuclease A (RNase A). Isothermal titration calorimetry experiments show that 3'-TMP binds the enzyme with a dissociation constant (K(d)) of 15 microM making it one of the strongest binding members of the five natural bases found in nucleic acids (A, C, G, T, and U). To further investigate the molecular properties of this potent natural affinity, we have determined the crystal structure of bovine pancreatic RNase A in complex with 3'-TMP at 1.55 A resolution and we have performed NMR binding experiments with 3'-CMP and 3'-TMP. Our results show that binding of 3'-TMP is very similar to other natural and non-natural pyrimidine ligands, demonstrating that single nucleotide affinity is independent of the presence or absence of a 2'-hydroxyl on the ribose moiety of pyrimidines and suggesting that the pyrimidine binding subsite of RNase A is not a significant contributor of inhibitor discrimination. Accumulating evidence suggests that very subtle structural, chemical, and potentially motional variations contribute to ligand discrimination in this enzyme.

Investigation of the catalytic mechanism of a synthetic DNAzyme with protein-like functionality: an RNaseA mimic?

The protein enzyme ribonuclease A (RNaseA) cleaves RNA with catalytic perfection, although with little sequence specificity, by a divalent metal ion (M(2+))-independent mechanism in which a pair of imidazoles provides general acid and base catalysis, while a cationic amine provides electrostatic stabilization of the transition state. Synthetic imitation of this remarkable organo-catalyst ("RNaseA mimicry") has been a longstanding goal in biomimetic chemistry. The 9(25)-11 DNAzyme contains synthetically modified nucleotides presenting both imidazole and cationic amine side chains, and catalyzes RNA cleavage with turnover in the absence of M(2+) similarly to RNaseA. Nevertheless, the catalytic roles, if any, of the "protein-like" functional groups have not been defined, and hence the question remains whether 9(25)-11 engages any of these functionalities to mimic aspects of the mechanism of RNaseA. To address this question, we report a mechanistic investigation of 9(25)-11 catalysis wherein we have employed a variety of experiments, such as DNAzyme functional group deletion, mechanism-based affinity labeling, and bridging and nonbridging phosphorothioate substitution of the scissile phosphate. Several striking parallels exist between the results presented here for 9(25)-11 and the results of analogous experiments applied previously to RNaseA. Specifically, our results implicate two particular imidazoles in general acid and base catalysis and suggest that a specific cationic amine stabilizes the transition state via diastereoselective interaction with the scissile phosphate. Overall, 9(25)-11 appears to meet the minimal criteria of an RNaseA mimic; this demonstrates how added synthetic functionality can expand the mechanistic repertoire available to a synthetic DNA-based catalyst.

Molecular dynamics simulation of bovine pancreatic ribonuclease A-CpA and transition state-like complexes

The mechanisms of enzymes are intimately connected with their overall structure and dynamics in solution. Experimentally, it is considerably challenging to provide detailed atomic level information about the conformational events that occur at different stages along the chemical reaction path. Here, theoretical tools may offer new potential insights that complement those obtained from experiments that may not yield an unambiguous mechanistic interpretation. In this study, we apply molecular dynamics simulations of bovine pancreatic ribonuclease A, an archetype ribonuclease, to study the conformational dynamics, structural relaxation, and differential solvation that occur at discrete stages of the transesterification and cleavage reaction. Simulations were performed with explicit solvation with rigorous electrostatics and utilize recently developed molecular mechanical force field parameters for transphosphorylation and hydrolysis transition state analogues. Herein, we present results for the enzyme complexed with the dinucleotide substrate cytidilyl-3',5'-adenosine (CpA) in the reactant, and transphosphorylation and hydrolysis transition states. A detailed analysis of active site structures and hydrogen-bond patterns is presented and compared. The integrity of the overall backbone structure is preserved in the simulations and supports a mechanism whereby His12 stabilizes accumulating negative charge at the transition states through hydrogen-bond donation to the nonbridge oxygens. Lys41 is shown to be highly versatile along the reaction coordinate and can aid in the stabilization of the dianionic transition state, while being poised to act as a general acid catalyst in the hydrolysis step.

The flexibility of a distant loop modulates active site motion and product release in ribonuclease A

The role of the flexible loop 1 in protein conformational motion and in the dissociation of enzymatic product from ribonuclease A (RNase A) was investigated by creation of a chimeric enzyme in which a 6-residue loop 1 from the RNase A homologue, eosinophil cationic protein (ECP), replaced the 12-residue loop 1 in RNase A. The chimera (RNase A(ECP)) experiences only local perturbations in NMR backbone chemical shifts compared to WT RNase A. Many of the flexible residues that were previously identified in WT as involved in an important conformational change now experience no NMR-detected millisecond motions in the chimera. Likewise, binding of the product analogue, 3'-CMP, to RNase A(ECP) results in only minor chemical shift changes in the enzyme similar to what is observed for the H48A mutant of RNase A and in contrast to WT enzyme. For both RNase A(ECP) and H48A there is a 10-fold decrease in the product release rate constant, k(off), compared to WT, in agreement with previous studies indicating the importance of flexibility in RNase A in the overall rate-limiting product release step. Together, these NMR and biochemical experiments provide additional insight into the mechanism of millisecond motions in the RNase A catalytic cycle.

The 3'-flap pocket of human flap endonuclease 1 is critical for substrate binding and catalysis

Flap endonuclease 1 (FEN1) proteins, which are present in all kingdoms of life, catalyze the sequence-independent hydrolysis of the bifurcated nucleic acid intermediates formed during DNA replication and repair. How FEN1s have evolved to preferentially cleave flap structures is of great interest especially in light of studies wherein mice carrying a catalytically deficient FEN1 were predisposed to cancer. Structural studies of FEN1s from phage to human have shown that, although they share similar folds, the FEN1s of higher organisms contain a 3'-extrahelical nucleotide (3'-flap) binding pocket. When presented with 5'-flap substrates having a 3'-flap, archaeal and eukaryotic FEN1s display enhanced reaction rates and cleavage site specificity. To investigate the role of this interaction, a kinetic study of human FEN1 (hFEN1) employing well defined DNA substrates was conducted. The presence of a 3'-flap on substrates reduced Km and increased multiple- and single turnover rates of endonucleolytic hydrolysis at near physiological salt concentrations. Exonucleolytic and fork-gap-endonucleolytic reactions were also stimulated by the presence of a 3'-flap, and the absence of a 3'-flap from a 5'-flap substrate was more detrimental to hFEN1 activity than removal of the 5'-flap or introduction of a hairpin into the 5'-flap structure. hFEN1 reactions were predominantly rate-limited by product release regardless of the presence or absence of a 3'-flap. Furthermore, the identity of the stable enzyme product species was deduced from inhibition studies to be the 5'-phosphorylated product. Together the results indicate that the presence of a 3'-flap is the critical feature for efficient hFEN1 substrate recognition and catalysis.

An in silico model for HIF-alpha regulation and hypoxia response in tumor cells

The dependency of the growth and metastasis of tumors on the new blood vessel formation, or angiogenesis, has opened up new potentials to tumor therapy, nevertheless understanding the molecular mechanisms involved in angiogenesis is crucial in the bioengineering of novel anti-angiogenic drugs. The key component in hypoxia sensing in tumor cells is the hypoxia-inducible factor, HIF-1alpha, which is inactivated through proteosome-mediated degradation under normoxic conditions. Two enzymes have been reported to hydroxylate HIF-1alpha, namely prolyl hydroxylase (PH), recruiting the proetasome complex and degrading cytoplasmic HIF-1alpha, and asparaginyl hydroxylase/factor inhibiting HIF-1alpha (FIH-1), downregulating the recruitment of p300 to the promoter, thereby reducing the transcriptional activity of HIF-1alpha. In this study, we have constructed an in silico model of a tumor cell using the GEPASI 3.30 biochemical simulation software (http://www.gepasi.org) and studied the performances of PH and FIH-1 on HIF-1alpha degradation and inactivation, respectively, as monitored by expression of the vascular endothelial growth factor, VEGF, during hypoxia. In our biochemical models, FIH-1 can successfully increase hypoxic transcription of VEGF, however FIH-1 on its own is not sufficient to inactivate HIF-1 completely, leading to background VEGF transcription under normoxic conditions. On the other hand, PH is necessary to increase the hypoxic transcriptional response, and can effectively shut off normoxic transcription. We therefore propose that regulating PH activity can be a primary target for anti-angiogenic bioengineering research.

Back to the future: Ribonuclease A

Pancreatic ribonuclease A (EC 3.1.27.5, RNase) is, perhaps, the best-studied enzyme of the 20th century. It was isolated by René Dubos, crystallized by Moses Kunitz, sequenced by Stanford Moore and William Stein, and synthesized in the laboratory of Bruce Merrifield, all at the Rockefeller Institute/University. It has proven to be an excellent model system for many different types of experiments, both as an enzyme and as a well-characterized protein for biophysical studies. Of major significance was the demonstration by Chris Anfinsen at NIH that the primary sequence of RNase encoded the three-dimensional structure of the enzyme. Many other prominent protein chemists/enzymologists have utilized RNase as a dominant theme in their research. In this review, the history of RNase and its offspring, RNase S (S-protein/S-peptide), will be considered, especially the work in the Merrifield group, as a preface to preliminary data and proposed experiments addressing topics of current interest. These include entropy-enthalpy compensation, entropy of ligand binding, the impact of protein modification on thermal stability, and the role of protein dynamics in enzyme action. In continuing to use RNase as a prototypical enzyme, we stand on the shoulders of the giants of protein chemistry to survey the future.

The mechanism of rate-limiting motions in enzyme function

The ability to use conformational flexibility is a hallmark of enzyme function. Here we show that protein motions and catalytic activity in a RNase are coupled and display identical solvent isotope effects. Solution NMR relaxation experiments identify a cluster of residues, some distant from the active site, that are integral to this motion. These studies implicate a single residue, histidine-48, as the key modulator in coupling protein motion with enzyme function. Mutation of H48 to alanine results in loss of protein motion in the isotope-sensitive region of the enzyme. In addition, k(cat) decreases for this mutant and the kinetic solvent isotope effect on k(cat), which was 2.0 in WT, is near unity in H48A. Despite being located 18 A from the enzyme active site, H48 is essential in coordinating the motions involved in the rate-limiting enzymatic step. These studies have identified, of approximately 160 potential exchangeable protons, a single site that is integral in the rate-limiting step in RNase A enzyme function.

Inhibition of human pancreatic ribonuclease by the human ribonuclease inhibitor protein

The ribonuclease inhibitor protein (RI) binds to members of the bovine pancreatic ribonuclease (RNase A) superfamily with an affinity in the femtomolar range. Here, we report on structural and energetic aspects of the interaction between human RI (hRI) and human pancreatic ribonuclease (RNase 1). The structure of the crystalline hRI x RNase 1 complex was determined at a resolution of 1.95 A, revealing the formation of 19 intermolecular hydrogen bonds involving 13 residues of RNase 1. In contrast, only nine such hydrogen bonds are apparent in the structure of the complex between porcine RI and RNase A. hRI, which is anionic, also appears to use its horseshoe-shaped structure to engender long-range Coulombic interactions with RNase 1, which is cationic. In accordance with the structural data, the hRI.RNase 1 complex was found to be extremely stable (t(1/2)=81 days; K(d)=2.9 x 10(-16) M). Site-directed mutagenesis experiments enabled the identification of two cationic residues in RNase 1, Arg39 and Arg91, that are especially important for both the formation and stability of the complex, and are thus termed "electrostatic targeting residues". Disturbing the electrostatic attraction between hRI and RNase 1 yielded a variant of RNase 1 that maintained ribonucleolytic activity and conformational stability but had a 2.8 x 10(3)-fold lower association rate for complex formation and 5.9 x 10(9)-fold lower affinity for hRI. This variant of RNase 1, which exhibits the largest decrease in RI affinity of any engineered ribonuclease, is also toxic to human erythroleukemia cells. Together, these results provide new insight into an unusual and important protein-protein interaction, and could expedite the development of human ribonucleases as chemotherapeutic agents.

Potent inhibition of ribonuclease A by oligo(vinylsulfonic acid)

Ribonuclease A (RNase A) can make multiple contacts with an RNA substrate. In particular, the enzymatic active site and adjacent subsites bind sequential phosphoryl groups in the RNA backbone through Coulombic interactions. Here, oligomers of vinylsulfonic acid (OVS) are shown to be potent inhibitors of RNase A that exploit these interactions. Inhibition is competitive with substrate and has Ki = 11 pm in assays at low salt concentration. The effect of salt concentration on inhibition indicates that nearly eight favorable Coulombic interactions occur in the RNase A.OVS complex. The phosphonic acid and sulfuric acid analogs of OVS are also potent inhibitors although slightly less effective. OVS is also shown to be a contaminant of MES and other buffers that contain sulfonylethyl groups. Oligomers greater than nine units in length can be isolated from commercial MES buffer. Inhibition by contaminating OVS is responsible for the apparent decrease in catalytic activity that has been observed in assays of RNase A at low salt concentration. Thus, OVS is both a useful inhibitor of RNase A and a potential bane to chemists and biochemists who use ethanesulfonic acid buffers.