Structural insights into the effects of 2′-5′ linkages on the RNA duplex

Significance The nonenzymatic replication of RNA is thought to have been a critical step in the emergence of simple cellular life from prebiotic chemistry. However, the chemical copying of RNA templates generates product strands that contain 2′-5′ backbone linkages and normal 3′-5′ linkages. Our recent finding that RNAs with such mixed backbones can still fold into functional structures raised the question of how RNA accommodates the presence of 2′-5′ linkages. Here we use X-ray crystallography and molecular dynamics simulations to reveal how 3′-5′–linked RNA duplexes accommodate interspersed 2′-5′ linkages. The diminished thermal and chemical stability of such RNA duplexes reflects local structural changes, but compensatory changes result in a global RNA duplex structure with relatively minor alterations. The mixture of 2′-5′ and 3′-5′ linkages generated during the nonenzymatic replication of RNA has long been regarded as a central problem for the origin of the RNA world. However, we recently observed that both a ribozyme and an RNA aptamer retain considerable functionality in the presence of prebiotically plausible levels of linkage heterogeneity. To better understand the RNA structure and function in the presence of backbone linkage heterogeneity, we obtained high-resolution X-ray crystal structures of a native 10-mer RNA duplex (1.32 Å) and two variants: one containing one 2′-5′ linkage per strand (1.55 Å) and one containing three such linkages per strand (1.20 Å). We found that RNA duplexes adjust their local structures to accommodate the perturbation caused by 2′-5′ linkages, with the flanking nucleotides buffering the disruptive effects of the isomeric linkage and resulting in a minimally altered global structure. Although most 2′-linked sugars were in the expected 2′-endo conformation, some were partially or fully in the 3′-endo conformation, suggesting that the energy difference between these conformations was relatively small. Our structural and molecular dynamic studies also provide insight into the diminished thermal and chemical stability of the duplex state associated with the presence of 2′-5′ linkages. Our results contribute to the view that a low level of 2′-5′ substitution would not have been fatal in an early RNA world and may in contrast have been helpful for both the emergence of nonenzymatic RNA replication and the early evolution of functional RNAs.

[1]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[2]  R. Sarma,et al.  Why do nucleic acids have 3′5′ phosphodiester bonds? , 1978, Nature.

[3]  J. Szostak,et al.  Evolution of functional nucleic acids in the presence of nonheritable backbone heterogeneity , 2011, Proceedings of the National Academy of Sciences.

[4]  M. Viswamitra,et al.  Visualisation of a 2'-5' parallel stranded double helix at atomic resolution: crystal structure of cytidylyl-2',5'-adenosine. , 1991, Nucleic Acids Research.

[5]  D. A. Usher,et al.  Hydrolytic stability of helical RNA: a selective advantage for the natural 3',5'-bond. , 1976, Proceedings of the National Academy of Sciences of the United States of America.

[6]  G. F. Joyce RNA evolution and the origins of life , 1989, Nature.

[7]  H. M. Sobell,et al.  Visualization of drug--nucleic acid interactions at atomic resolution. VI. Structure of two drug--dinucleoside monophosphate crystalline complexes, ellipticine--5-iodocytidylyy (3'-5') guanosine and 3,5,6,8-tetramethyl-N-methyl phenanthrolinium--5-iodocytidylyl (3'-5') guanosine. , 1979, Journal of molecular biology.

[8]  Z. Otwinowski,et al.  Processing of X-ray diffraction data collected in oscillation mode. , 1997, Methods in enzymology.

[9]  J W Szostak,et al.  Nonenzymatic, template-directed ligation of oligoribonucleotides is highly regioselective for the formation of 3'-5' phosphodiester bonds. , 1996, Journal of the American Chemical Society.

[10]  M. Woolfson,et al.  A flexible and efficient procedure for the solution and phase refinement of protein structures. , 2000, Acta crystallographica. Section D, Biological crystallography.

[11]  D. Bartel,et al.  RNA-Catalyzed RNA Polymerization: Accurate and General RNA-Templated Primer Extension , 2001, Science.

[12]  G. Murshudov,et al.  Refinement of macromolecular structures by the maximum-likelihood method. , 1997, Acta crystallographica. Section D, Biological crystallography.

[13]  Garib N. Murshudov,et al.  JLigand: a graphical tool for the CCP4 template-restraint library , 2012, Acta crystallographica. Section D, Biological crystallography.

[14]  A. Noronha,et al.  Physicochemical and biochemical properties of 2',5'-linked RNA and 2',5'-RNA:3',5'-RNA "hybrid" duplexes. , 1998, Biochemistry.

[15]  L. Orgel Evolution of the genetic apparatus. , 1968, Journal of molecular biology.

[16]  Bert L. de Groot,et al.  g_wham—A Free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates , 2010 .

[17]  B. Roux The calculation of the potential of mean force using computer simulations , 1995 .

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

[19]  R. K. Stodola,et al.  Drug-nucleic acid interactions: conformational flexibility at the intercalation site. , 1978, Proceedings of the National Academy of Sciences of the United States of America.

[20]  R Parthasarathy,et al.  X-ray structure of a dinucleoside monophosphate A2'p5'C that contains a 2'-5' link found in (2'-5')oligo(A)s induced by interferons: single-stranded helical conformation of 2'-5'-linked oligonucleotides. , 1982, Proceedings of the National Academy of Sciences of the United States of America.

[21]  J. H. Maddocks,et al.  Conformational analysis of nucleic acids revisited: Curves+ , 2009, Nucleic acids research.

[22]  J. Ferris,et al.  Oligomerization of ribonucleotides on montmorillonite: reaction of the 5'-phosphorimidazolide of adenosine. , 1992, Science.

[23]  H M Berman,et al.  A standard reference frame for the description of nucleic acid base-pair geometry. , 2001, Journal of molecular biology.

[24]  N. Hud,et al.  Assembly of an antiparallel homo-adenine DNA duplex by small-molecule binding. , 2004, Journal of the American Chemical Society.

[25]  F. H. C. CRICK,et al.  Origin of the Genetic Code , 1967, Nature.

[26]  P. A. Giannaris,et al.  Oligoribonucleotides containing 2',5'-phosphodiester linkages exhibit binding selectivity for 3',5'-RNA over 3',5'-ssDNA. , 1993, Nucleic acids research.

[27]  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..

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

[29]  R R Breaker,et al.  Relationship between internucleotide linkage geometry and the stability of RNA. , 1999, RNA.

[30]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[31]  E. Kool,et al.  Encoding phenotype in bacteria with an alternative genetic set. , 2011, Journal of the American Chemical Society.

[32]  B. Ganem RNA world , 1987, Nature.

[33]  Z. Otwinowski,et al.  [20] Processing of X-ray diffraction data collected in oscillation mode. , 1997, Methods in enzymology.

[34]  J. Szostak,et al.  Functional RNAs exhibit tolerance for non-heritable 2'-5' versus 3'-5' backbone heterogeneity. , 2013, Nature chemistry.

[35]  J. Plavec,et al.  Solution structure of a modified 2′,5′-linked RNA hairpin involved in an equilibrium with duplex , 2005, Nucleic acids research.

[36]  M. Sundaralingam,et al.  Solution structure of 2',5' d(G4C4). Relevance to topological restrictions and nature's choice of phosphodiester links. , 2004, European journal of biochemistry.

[37]  R J Read,et al.  Crystallography & NMR system: A new software suite for macromolecular structure determination. , 1998, Acta crystallographica. Section D, Biological crystallography.

[38]  R. Dickerson,et al.  Definitions and nomenclature of nucleic acid structure components. , 1989, Nucleic acids research.

[39]  E. Ohtsuka,et al.  Identification of I:A mismatch base-pairing structure in DNA. , 1987, The Journal of biological chemistry.

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

[41]  L. Orgel,et al.  Catalysis of accurate poly(C)-directed synthesis of 3'-5'-linked oligoguanylates by Zn2+. , 1980, Journal of molecular biology.

[42]  David P. Bartel,et al.  RNA-catalysed RNA polymerization using nucleoside triphosphates , 1996, Nature.

[43]  L. Orgel,et al.  Oligomerization of (guanosine 5'-phosphor)-2-methylimidazolide on poly(C). An RNA polymerase model. , 1982, Journal of molecular biology.

[44]  S. Grzesiek,et al.  NMRPipe: A multidimensional spectral processing system based on UNIX pipes , 1995, Journal of biomolecular NMR.

[45]  Zaida Luthey-Schulten,et al.  Dynamics of Recognition between tRNA and elongation factor Tu. , 2008, Journal of molecular biology.

[46]  C. Sherrill,et al.  Quantum-mechanical analysis of the energetic contributions to π stacking in nucleic acids versus rise, twist, and slide. , 2013, Journal of the American Chemical Society.

[47]  D. Langreth,et al.  An application of the van der Waals density functional: Hydrogen bonding and stacking interactions between nucleobases. , 2008, The Journal of chemical physics.

[48]  H. Berman,et al.  New parameters for the refinement of nucleic acid-containing structures. , 1996, Acta crystallographica. Section D, Biological crystallography.

[49]  A. H. Wang,et al.  DNA with 2'-5' Phosphodiester Bonds Forms a Duplex Structure in the A-type Conformation , 1995 .

[50]  N. Yathindra,et al.  Structural basis for the unusual properties of 2',5' nucleic acids and their complexes with RNA and DNA. , 2002, Biophysical chemistry.

[51]  C. Switzer,et al.  Activity of 2′,5′‐Linked RNA in the Template‐Directed Oligomerization of Mononucleotides , 1997 .

[52]  Christopher K. W. Chan,et al.  Prebiotically plausible oligoribonucleotide ligation facilitated by chemoselective acetylation. , 2013, Nature chemistry.

[53]  E J Dodson,et al.  ACORN2: new developments of the ACORN concept. , 2009, Acta crystallographica. Section D, Biological crystallography.

[54]  Kevin Cowtan,et al.  research papers Acta Crystallographica Section D Biological , 2005 .

[55]  S. Agrawal,et al.  NMR structure of a 2',5' RNA favors A type duplex with compact C2'endo nucleotide repeat. , 2001, Biochemical and biophysical research communications.

[56]  D. Turner,et al.  Association of 2'-5' oligoribonucleotides. , 1992, Nucleic acids research.

[57]  Xiang-Jun Lu,et al.  3DNA: a versatile, integrated software system for the analysis, rebuilding and visualization of three-dimensional nucleic-acid structures , 2008, Nature Protocols.

[58]  N. Hud,et al.  Solution structure and thermodynamics of 2',5' RNA intercalation. , 2009, Journal of the American Chemical Society.

[59]  R. Read Improved Fourier Coefficients for Maps Using Phases from Partial Structures with Errors , 1986 .

[60]  Randy J. Read,et al.  Phaser crystallographic software , 2007, Journal of applied crystallography.

[61]  Shri Jain,et al.  Visualization of drug-nucleic acid interactions at atomic resolution. III. Unifying structural concepts in understanding drug-DNA interactions and their broader implications in understanding protein-DNA interactions. , 1977, Journal of molecular biology.