1
|
Ekesan Ş, York DM. Dynamical ensemble of the active state and transition state mimic for the RNA-cleaving 8-17 DNAzyme in solution. Nucleic Acids Res 2019; 47:10282-10295. [PMID: 31511899 PMCID: PMC6821293 DOI: 10.1093/nar/gkz773] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2019] [Revised: 08/20/2019] [Accepted: 09/03/2019] [Indexed: 02/01/2023] Open
Abstract
We perform molecular dynamics simulations, based on recent crystallographic data, on the 8-17 DNAzyme at four states along the reaction pathway to determine the dynamical ensemble for the active state and transition state mimic in solution. A striking finding is the diverse roles played by Na+ and Pb2+ ions in the electrostatically strained active site that impact all four fundamental catalytic strategies, and share commonality with some features recently inferred for naturally occurring hammerhead and pistol ribozymes. The active site Pb2+ ion helps to stabilize in-line nucleophilic attack, provides direct electrostatic transition state stabilization, and facilitates leaving group departure. A conserved guanine residue is positioned to act as the general base, and is assisted by a bridging Na+ ion that tunes the pKa and facilitates in-line fitness. The present work provides insight into how DNA molecules are able to solve the RNA-cleavage problem, and establishes functional relationships between the mechanism of these engineered DNA enzymes with their naturally evolved RNA counterparts. This adds valuable information to our growing body of knowledge on general mechanisms of phosphoryl transfer reactions catalyzed by RNA, proteins and DNA.
Collapse
Affiliation(s)
- Şölen Ekesan
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| |
Collapse
|
2
|
Wierzchowski J, Antosiewicz JM, Shugar D. 8-Azapurines as isosteric purine fluorescent probes for nucleic acid and enzymatic research. MOLECULAR BIOSYSTEMS 2015; 10:2756-74. [PMID: 25124808 DOI: 10.1039/c4mb00233d] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The 8-azapurines, and their 7-deaza and 9-deaza congeners, represent a unique class of isosteric (isomorphic) analogues of the natural purines, frequently capable of substituting for the latter in many biochemical processes. Particularly interesting is their propensity to exhibit pH-dependent room-temperature fluorescence in aqueous medium, and in non-polar media. We herein review the physico-chemical properties of this class of compounds, with particular emphasis on the fluorescence emission properties of their neutral and/or ionic species, which has led to their widespread use as fluorescent probes in enzymology, including enzymes involved in purine metabolism, agonists/antagonists of adenosine receptors, mechanisms of catalytic RNAs, RNA editing, etc. They are also exceptionally useful fluorescent probes for analytical and clinical applications in crude cell homogenates.
Collapse
Affiliation(s)
- Jacek Wierzchowski
- Department of Biophysics, University of Varmia & Masuria, Oczapowskiego 4, 10-719 Olsztyn, Poland.
| | | | | |
Collapse
|
3
|
Panteva MT, Dissanayake T, Chen H, Radak BK, Kuechler ER, Giambaşu GM, Lee TS, York DM. Multiscale methods for computational RNA enzymology. Methods Enzymol 2015; 553:335-74. [PMID: 25726472 PMCID: PMC4739856 DOI: 10.1016/bs.mie.2014.10.064] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
RNA catalysis is of fundamental importance to biology and yet remains ill-understood due to its complex nature. The multidimensional "problem space" of RNA catalysis includes both local and global conformational rearrangements, changes in the ion atmosphere around nucleic acids and metal ion binding, dependence on potentially correlated protonation states of key residues, and bond breaking/forming in the chemical steps of the reaction. The goal of this chapter is to summarize and apply multiscale modeling methods in an effort to target the different parts of the RNA catalysis problem space while also addressing the limitations and pitfalls of these methods. Classical molecular dynamics simulations, reference interaction site model calculations, constant pH molecular dynamics (CpHMD) simulations, Hamiltonian replica exchange molecular dynamics, and quantum mechanical/molecular mechanical simulations will be discussed in the context of the study of RNA backbone cleavage transesterification. This reaction is catalyzed by both RNA and protein enzymes, and here we examine the different mechanistic strategies taken by the hepatitis delta virus ribozyme and RNase A.
Collapse
Affiliation(s)
- Maria T Panteva
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Thakshila Dissanayake
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Haoyuan Chen
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Brian K Radak
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Erich R Kuechler
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - George M Giambaşu
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Tai-Sung Lee
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA
| | - Darrin M York
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, USA.
| |
Collapse
|
4
|
Radak BK, Harris ME, York DM. Molecular simulations of RNA 2'-O-transesterification reaction models in solution. J Phys Chem B 2013; 117:94-103. [PMID: 23214417 PMCID: PMC3574632 DOI: 10.1021/jp3084277] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
We employ quantum mechanical/molecular mechanical umbrella sampling simulations to probe the free energy surfaces of a series of increasingly complex reaction models of RNA 2'-O-transesterification in aqueous solution under alkaline conditions. Such models are valuable for understanding the uncatalyzed processes underlying catalytic cleavage of the phosphodiester backbone of RNA, a reaction of fundamental importance in biology. The chemically reactive atoms are modeled by the AM1/d-PhoT quantum model for phosphoryl transfer, whereas the aqueous solvation environment is modeled with a molecular mechanics force field. Several simulation protocols were compared that used different ionic conditions and force field models. The results provide insight into how variation of the structural environment of the nucleophile and leaving group affects the free energy profile for the transesterification reaction. Results for a simple RNA backbone model are compared with recent experiments by Harris et al. on the specific base-catalyzed cleavage of a UpG dinucleotide. The calculated and measured free energies of activation match extremely well (ΔF(‡) = 19.9-20.8 vs 19.9 kcal/mol). Solvation is seen to play a crucial role and is characterized by a network of hydrogen bonds that envelopes the pentacoordinate dianionic phosphorane transition state and provides preferential stabilization relative to the reactant state.
Collapse
Affiliation(s)
- Brian K. Radak
- BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087 USA
- Department of Chemistry, University of Minnesota, Minneapolis, MN 55455-0431 USA
| | - Michael E. Harris
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, OH 44106 USA
| | - Darrin M. York
- BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087 USA
| |
Collapse
|
5
|
Abstract
Metal ions are inextricably involved with nucleic acids due to their polyanionic nature. In order to understand the structure and function of RNAs and DNAs, one needs to have detailed pictures on the structural, thermodynamic, and kinetic properties of metal ion interactions with these biomacromolecules. In this review we first compile the physicochemical properties of metal ions found and used in combination with nucleic acids in solution. The main part then describes the various methods developed over the past decades to investigate metal ion binding by nucleic acids in solution. This includes for example hydrolytic and radical cleavage experiments, mutational approaches, as well as kinetic isotope effects. In addition, spectroscopic techniques like EPR, lanthanide(III) luminescence, IR and Raman as well as various NMR methods are summarized. Aside from gaining knowledge about the thermodynamic properties on the metal ion-nucleic acid interactions, especially NMR can be used to extract information on the kinetics of ligand exchange rates of the metal ions applied. The final section deals with the influence of anions, buffers, and the solvent permittivity on the binding equilibria between metal ions and nucleic acids. Little is known on some of these aspects, but it is clear that these three factors have a large influence on the interaction between metal ions and nucleic acids.
Collapse
Affiliation(s)
- Maria Pechlaner
- Institute of Inorganic Chemistry, University of Zürich, Zürich, Switzerland
| | | |
Collapse
|
6
|
Lassila JK, Zalatan JG, Herschlag D. Biological phosphoryl-transfer reactions: understanding mechanism and catalysis. Annu Rev Biochem 2011; 80:669-702. [PMID: 21513457 DOI: 10.1146/annurev-biochem-060409-092741] [Citation(s) in RCA: 286] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Phosphoryl-transfer reactions are central to biology. These reactions also have some of the slowest nonenzymatic rates and thus require enormous rate accelerations from biological catalysts. Despite the central importance of phosphoryl transfer and the fascinating catalytic challenges it presents, substantial confusion persists about the properties of these reactions. This confusion exists despite decades of research on the chemical mechanisms underlying these reactions. Here we review phosphoryl-transfer reactions with the goal of providing the reader with the conceptual and experimental background to understand this body of work, to evaluate new results and proposals, and to apply this understanding to enzymes. We describe likely resolutions to some controversies, while emphasizing the limits of our current approaches and understanding. We apply this understanding to enzyme-catalyzed phosphoryl transfer and provide illustrative examples of how this mechanistic background can guide and deepen our understanding of enzymes and their mechanisms of action. Finally, we present important future challenges for this field.
Collapse
Affiliation(s)
- Jonathan K Lassila
- Department of Biochemistry, Stanford University, Stanford, California 94305, USA.
| | | | | |
Collapse
|
7
|
Georgianna WE, Young DD. Development and utilization of non-coding RNA–small molecule interactions. Org Biomol Chem 2011; 9:7969-78. [DOI: 10.1039/c1ob06324c] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
8
|
Harris ME, Dai Q, Gu H, Kellerman DL, Piccirilli JA, Anderson VE. Kinetic isotope effects for RNA cleavage by 2'-O- transphosphorylation: nucleophilic activation by specific base. J Am Chem Soc 2010; 132:11613-21. [PMID: 20669950 DOI: 10.1021/ja103550e] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To better understand the interactions between catalysts and transition states during RNA strand cleavage, primary (18)O kinetic isotope effects (KIEs) and solvent D(2)O isotope effects were measured to probe the mechanism of base-catalyzed 2'-O-transphosphorylation of the RNA dinucleotide 5'-UpG-3'. The observed (18)O KIEs for the nucleophilic 2'-O and in the 5'-O leaving group at pH 14 are both large relative to reactions of phosphodiesters with good leaving groups, indicating that the reaction catalyzed by hydroxide has a transition state (TS) with advanced phosphorus-oxygen bond fission to the leaving group ((18)k(LG) = 1.034 +/- 0.004) and phosphorus-nucleophile bond formation ((18)k(NUC) = 0.984 +/- 0.004). A breakpoint in the pH dependence of the 2'-O-transphosphorylation rate to a pH independent phase above pH 13 has been attributed to the pK(a) of the 2'-OH nucleophile. A smaller nucleophile KIE is observed at pH 12 ((18)k(NUC) = 0.995 +/- 0.004) that is interpreted as the combined effect of the equilibrium isotope effect (ca. 1.02) on deprotonation of the 2'-hydroxyl nucleophile and the intrinsic KIE on the nucleophilic addition step (ca. 0.981). An alternative mechanism in which the hydroxide ion acts as a general base is considered unlikely given the lack of a solvent deuterium isotope effect above the breakpoint in the pH versus rate profile. These results represent the first direct analysis of the transition state for RNA strand cleavage. The primary (18)O KIE results and the lack of a kinetic solvent deuterium isotope effect together provide strong evidence for a late transition state and 2'-O nucleophile activation by specific base catalysis.
Collapse
Affiliation(s)
- Michael E Harris
- RNA Center and Departmet of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, Ohio 44118, USA.
| | | | | | | | | | | |
Collapse
|
9
|
Rios AC, Tor Y. Model systems: how chemical biologists study RNA. Curr Opin Chem Biol 2009; 13:660-8. [PMID: 19879179 DOI: 10.1016/j.cbpa.2009.09.028] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2009] [Revised: 09/14/2009] [Accepted: 09/25/2009] [Indexed: 11/30/2022]
Abstract
Ribonucleic acids are structurally and functionally sophisticated biomolecules and the use of models, frequently truncated or modified sequences representing functional domains of the natural systems, is essential to their exploration. Functional noncoding RNAs such as miRNAs, riboswitches, and, in particular, ribozymes, have changed the view of RNA's role in biology and its catalytic potential. The well-known truncated hammerhead model has recently been refined and new data provide a clearer molecular picture of the elements responsible for its catalytic power. A model for the spliceosome, a massive and highly intricate ribonucleoprotein, is also emerging, although its true utility is yet to be cemented. Such catalytic model systems could also serve as 'chemo-paleontological' tools, further refining the RNA world hypothesis and its relevance to the origin and evolution of life.
Collapse
Affiliation(s)
- Andro C Rios
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, CA 92093-0358, USA
| | | |
Collapse
|
10
|
Banáš P, Jurečka P, Walter NG, Šponer J, Otyepka M. Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM. Methods 2009; 49:202-16. [PMID: 19398008 PMCID: PMC2753711 DOI: 10.1016/j.ymeth.2009.04.007] [Citation(s) in RCA: 75] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2009] [Revised: 04/07/2009] [Accepted: 04/07/2009] [Indexed: 11/28/2022] Open
Abstract
Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.
Collapse
Affiliation(s)
- Pavel Banáš
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. Svobody 26, 771 46 Olomouc, Czech Republic
| | - Petr Jurečka
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. Svobody 26, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Nils G. Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 N. University Ave., Ann Arbor, MI 48109-1055, USA
| | - Jiří Šponer
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. Svobody 26, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. Svobody 26, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| |
Collapse
|
11
|
Frederiksen JK, Piccirilli JA. Identification of catalytic metal ion ligands in ribozymes. Methods 2009; 49:148-66. [PMID: 19651216 DOI: 10.1016/j.ymeth.2009.07.005] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2009] [Revised: 07/20/2009] [Accepted: 07/29/2009] [Indexed: 01/05/2023] Open
Abstract
Site-bound metal ions participate in the catalytic mechanisms of many ribozymes. Understanding these mechanisms therefore requires knowledge of the specific ligands on both substrate and ribozyme that coordinate these catalytic metal ions. A number of different structural and biochemical strategies have been developed and refined for identifying metal ion binding sites within ribozymes, and for assessing the catalytic contributions of the metal ions bound at those sites. We review these approaches and provide examples of their application, focusing in particular on metal ion rescue experiments and their roles in the construction of the transition state models for the Tetrahymena group I and RNase P ribozymes.
Collapse
Affiliation(s)
- John K Frederiksen
- The Pritzker School of Medicine, The University of Chicago, Chicago, IL 60637, USA
| | | |
Collapse
|