Stereoselection in the diels–alderase ribozyme: A molecular dynamics study

The Diels‐Alderase ribozyme is an in vitro‐evolved ribonucleic acid enzyme that catalyzes a [4 + 2] cycloaddition reaction between an anthracene diene and a maleimide dienophile. The ribozyme can in principle be used to selectively synthesize only one product enantiomer, depending on which of the two entrances to the catalytic pocket, “front” or “back”, the substrate is permitted to use. Here, we investigate stereoselection and substrate recognition in the ribozyme by means of multiple molecular dynamics simulations, performed on each of the two substrates individually in the pocket, on the reactant state, and on the product state. The results are consistent with a binding mechanism in which the maleimide likely binds first followed by the anthracene, which enters preferentially through the front door. The free energy profiles for anthracene binding indicate that the pre‐(R,R)‐enantiomer conformation is slightly preferred, in agreement with the experimentally observed small enantiomeric excess of the (R,R)‐enantiomer of the product. The reactant state is stabilized by the simultaneous presence of both substrates bound to their binding sites in the hydrophobic pocket as well as by stacking interactions between them. © 2012 Wiley Periodicals, Inc.

[1]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[2]  Nils G Walter,et al.  Molecular dynamics simulations of RNA: an in silico single molecule approach. , 2007, Biopolymers.

[3]  Hans W. Horn,et al.  ELECTRONIC STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS: THE PROGRAM SYSTEM TURBOMOLE , 1989 .

[4]  Dario A. Estrin,et al.  Comparing and combining implicit ligand sampling with multiple steered molecular dynamics to study ligand migration processes in heme proteins , 2011, J. Comput. Chem..

[5]  C. Jarzynski.  Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[6]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[7]  Aaron T. Frank,et al.  Constructing RNA dynamical ensembles by combining MD and motionally decoupled NMR RDCs: new insights into RNA dynamics and adaptive ligand recognition , 2009, Nucleic acids research.

[8]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[9]  T. Cech,et al.  Self-splicing RNA: Autoexcision and autocyclization of the ribosomal RNA intervening sequence of tetrahymena , 1982, Cell.

[10]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[11]  Hongwei Huang,et al.  Molecular Dynamics Simulation and Binding Energy Calculation for Estimation of Oligonucleotide Duplex Thermostability in RNA-Based Therapeutics , 2011, J. Chem. Inf. Model..

[12]  R. Rosenfeld.  Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[13]  Modesto Orozco,et al.  Recent advances in the study of nucleic acid flexibility by molecular dynamics. , 2008, Current opinion in structural biology.

[14]  D. York,et al.  Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation. , 2008, Journal of the American Chemical Society.

[15]  A. Cavalli,et al.  Single-molecule pulling simulations can discern active from inactive enzyme inhibitors. , 2010, Journal of the American Chemical Society.

[16]  J. Doudna,et al.  Creative catalysis: pieces of the RNA world jigsaw. , 2001, Trends in biochemical sciences.

[17]  A. Jäschke,et al.  Characterization of an RNA active site: interactions between a Diels-Alderase ribozyme and its substrates and products. , 2002, Journal of the American Chemical Society.

[18]  Tamsyn Montagnon,et al.  The Diels--Alder reaction in total synthesis. , 2002, Angewandte Chemie.

[19]  Gerhard Hummer,et al.  Free energy surfaces from single-molecule force spectroscopy. , 2005, Accounts of chemical research.

[20]  R. Nussinov,et al.  Folding funnels and binding mechanisms. , 1999, Protein engineering.

[21]  Otto Diels,et al.  Synthesen in der hydroaromatischen Reihe , 1928 .

[22]  W. Gilbert.  Origin of life: The RNA world , 1986, Nature.

[23]  Alexander D. MacKerell,et al.  CHARMM force field parameters for simulation of reactive intermediates in native and thio‐substituted ribozymes , 2007, J. Comput. Chem..

[24]  G. F. Joyce,et al.  Directed evolution of nucleic acid enzymes. , 2003, Annual review of biochemistry.

[25]  J. Elenewski,et al.  Free energy landscape of the retinol/serum retinol binding protein complex: a biological host-guest system. , 2010, The journal of physical chemistry. B.

[26]  A. Serganov,et al.  Control of stereoselectivity in an enzymatic reaction by backdoor access. , 2006, Angewandte Chemie.

[27]  K. Schulten,et al.  Calculating potentials of mean force from steered molecular dynamics simulations. , 2004, The Journal of chemical physics.

[28]  Alexander D. MacKerell,et al.  Atomic detail investigation of the structure and dynamics of DNA.RNA hybrids: a molecular dynamics study. , 2008, The journal of physical chemistry. B.

[29]  D. York,et al.  Insight into the role of Mg in hammerhead ribozyme catalysis from X-ray crystallography and molecular dynamics simulation. , 2007, Journal of chemical theory and computation.

[30]  A. Serganov,et al.  Structural basis for Diels-Alder ribozyme-catalyzed carbon-carbon bond formation , 2005, Nature Structural &Molecular Biology.

[31]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[32]  Alexander D. MacKerell,et al.  Development and current status of the CHARMM force field for nucleic acids , 2000, Biopolymers.

[33]  Jeremy C. Smith,et al.  Molecular mechanics force field parameterization of the fluorescent probe rhodamine 6G using automated frequency matching , 2003, J. Comput. Chem..

[34]  E. Westhof,et al.  Molecular dynamics investigations of hammerhead ribozyme RNA , 1998, European Biophysics Journal.

[35]  Klaus Schulten,et al.  Sugar transport across lactose permease probed by steered molecular dynamics. , 2007, Biophysical journal.

[36]  P. Tavan,et al.  Ligand Binding: Molecular Mechanics Calculation of the Streptavidin-Biotin Rupture Force , 1996, Science.

[37]  P. Auffinger,et al.  A short guide for molecular dynamics simulations of RNA systems. , 2009, Methods.

[38]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[39]  A. Jäschke,et al.  Nucleic acid enzymes. , 2005, Current opinion in biotechnology.

[40]  A. Jäschke,et al.  Evolution of DNA and RNA as catalysts for chemical reactions. , 2000, Current opinion in chemical biology.

[41]  Alexander D. MacKerell,et al.  Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening. , 2003, Nucleic acids research.

[42]  Alexander D. MacKerell,et al.  Impact of 2′‐hydroxyl sampling on the conformational properties of RNA: Update of the CHARMM all‐atom additive force field for RNA , 2011, J. Comput. Chem..

[43]  D. Bartel,et al.  RNA-catalysed nucleotide synthesis , 1998, Nature.

[44]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[45]  M. Vogel,et al.  A structural basis for the activity of retro-Diels–Alder catalytic antibodies: Evidence for a catalytic aromatic residue , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[46]  Alexander D. MacKerell,et al.  All‐atom empirical force field for nucleic acids: II. Application to molecular dynamics simulations of DNA and RNA in solution , 2000 .

[47]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[48]  K. Schulten,et al.  Free energy calculation from steered molecular dynamics simulations using Jarzynski's equality , 2003 .

[49]  G. Hummer,et al.  Free energy reconstruction from nonequilibrium single-molecule pulling experiments , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[50]  A. Jäschke,et al.  A small catalytic RNA motif with Diels-Alderase activity. , 1999, Chemistry & biology.

[51]  Jeremy C. Smith,et al.  AFMM: A molecular mechanics force field vibrational parametrization program , 2005, Comput. Phys. Commun..

[52]  Gerhard Hummer,et al.  Free energy profiles from single-molecule pulling experiments , 2010, Proceedings of the National Academy of Sciences.

[53]  Andres Jäschke,et al.  Enantioselective Ribozyme Catalysis of a Bimolecular Cycloaddition Reaction , 2000 .

[54]  K. Réblová,et al.  Conformational transitions of flanking purines in HIV‐1 RNA dimerization initiation site kissing complexes studied by CHARMM explicit solvent molecular dynamics , 2008, Biopolymers.

[55]  T. Tarasow,et al.  RNA-catalysed carbon–carbon bond formation , 1997, Nature.

[56]  S A Benner,et al.  Modern metabolism as a palimpsest of the RNA world. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[57]  Jeremy C. Smith,et al.  Magnesium-dependent active-site conformational selection in the Diels-Alderase ribozyme. , 2010, Journal of the American Chemical Society.

[58]  P. Kollman,et al.  Encyclopedia of computational chemistry , 1998 .

[59]  J Andrew McCammon,et al.  Potentials of mean force for acetylcholine unbinding from the alpha7 nicotinic acetylcholine receptor ligand-binding domain. , 2006, Journal of the American Chemical Society.

[60]  Jack W. Szostak,et al.  In vitro evolution of a self-alkylatlng ribozyme , 1995, Nature.

[61]  Weihua Li,et al.  Exploring coumarin egress channels in human cytochrome p450 2a6 by random acceleration and steered molecular dynamics simulations , 2011, Proteins.