Free Energy Landscape and Conformational Kinetics of Hoogsteen Base-Pairing in DNA Vs RNA

Total Page:16

File Type:pdf, Size:1020Kb

Free Energy Landscape and Conformational Kinetics of Hoogsteen Base-Pairing in DNA Vs RNA bioRxiv preprint doi: https://doi.org/10.1101/2020.01.09.868018; this version posted January 9, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Manuscript submitted to BiophysicalJournal Article Free Energy Landscape and Conformational Kinetics of Hoogsteen Base-Pairing in DNA vs RNA D. Ray and I. Andricioaei ABSTRACT Genetic information is encoded in the DNA double helix which, in its physiological milieu, is characterized by the iconical Watson- Crick nucleobase pairing. Recent NMR relaxation experiments revealed the transient presence of an alternative, Hoogsteen base pairing pattern in naked DNA duplexes and estimated its relative stability and lifetime. In contrast, HG transitions in RNA were not observed. Understanding Hoogsteen (HG) base pairing is important because the underlying "breathing" can modulate significantly DNA/RNA recognition by proteins. However, a detailed mechanistic insight into the transition pathways and kinetics is still missing. We performed enhanced sampling simulation (with combined metadynamics and adaptive force bias method) and Markov State modeling to obtain accurate free energy, kinetics and the intermediates in the transition pathway between WC and HG base pair for both naked B-DNA and A-RNA duplexes. The Markov state model constructed from our unbiased MD simulation data revealed previously unknown complex extra-helical intermediates in this seemingly simple process of base pair conformation switching in B-DNA. Extending our calculation to A-RNA, for which HG base pair is not observed experimentally, resulted in relatively unstable single hydrogen bonded distorted Hoogsteen like base pair. Unlike B-DNA the transition pathway primarily involved base paired and intra-helical intermediates with transition timescales much higher than that of B-DNA. The seemingly obvious flip-over reaction coordinate, i.e., the glycosidic torsion angle is unable to resolve the intermediates; so a multidimensional picture, involving backbone dihedral angles and distance between atoms participating in hydrogen bonds, is required to gain insight into the molecular mechanism. SIGNIFICANCE Formation of unconventional Hoogsteen (HG) base pairing is an important problem in DNA biophysics owing to its key role in facilitating the binding of DNA repairing enzymes, proteins and drugs to damaged DNA. X-ray crystallography and NMR relaxation experiments revealed the presence of HG base pair in naked DNA duplex and protein- DNA complex but no HG base pair was observed in RNA. Molecular dynamics simulations could reproduce the experimental free energy cost of HG base pairing in DNA although a detailed mechanistic insight is still missing. We performed enhanced sampling simulation and Markov state modeling to obtain accurate free energy, kinetics and the intermediates in the transition pathway between WC and HG base pair for both B-DNA and A-RNA. INTRODUCTION The one-dimensional genetic information encoded in the sequence of DNA base pairs is intrinsically related to its three- dimensional structure described by the iconical Watson-Crick (WC) helix (1), with its specific hydrogen-bonding pattern between purine and pyrimidine (A-T, C-G) complementary nucleobases. Recently, a series of exciting NMR relaxation studies (2) on free DNA in solution have shown that the predominant WC base-pairing pattern between A-T and G-C nucleotides is in dynamic equilibrium with a less common and transient (µs-ms lifetime) Hoogsteen (HG) form (3), where the purine nucleotide in the base pair flips over. The WC to HG base pair transition, a form of “DNA breathing” (4), hints to the existence of a “secondary genetic code” signaled by base pair flipping, apart from the primary encoding in the nucleotide sequence. Crystallographic experiments have reported the presence of Hoogsteen base pairs only in specific contexts, such as when DNA binds to proteins like p53 tumor suppressor protein (5) or integration host factor (IHF) in E. coli (6). Detailed biochemical studies showed that HG base pairing plays a key role in circumventing interference in DNA replication mechanism by DNA polymerase by avoiding appearance of lesions (7, 8). In p53 binding sequence HG conformation exposes negative charged regions of nucleic acid to the positive arginine residue of the protein (5) leading to stability of the complex and playing crucial role in DNA recognition (9, 10). Examples of specific binding to DNA guided by HG base pair extends beyond proteins and even includes small molecules and drugs, which in turn facilitates its recognition by proteins and enzymes (11, 12). In a nutshell, HG base pairing therefore expands the structural and functional repertoire of Manuscript submitted to Biophysical Journal 1 bioRxiv preprint doi: https://doi.org/10.1101/2020.01.09.868018; this version posted January 9, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. D. Ray and I. Andricioaei duplex DNA beyond what can be achieved by Watson-Crick base-pairing alone. For example, HG base pairs are recognized by DNA repairing enzymes resulting in selective binding in the damaged regions of DNA, rich in syn-anti configuration over the anti-anti WC base pair rich normal DNA. Thus, apart from providing structural integrity to the damaged and deformed DNA, HG base pairs assist in the DNA repair mechanism (13). Apart from multiple appearance of HG base pairing in damaged or distorted DNA (14–16) several studies have confirmed presence of Hoogsteen conformation in undistorted DNA duplex (17, 18); the first X-ray crystallographic evidence of Hoogsteen base pairing in a simple stretch of DNA sequence composed exclusively of A-T base pairs was provided by Abrescia and coworkers (19). Early computational work by Gould and Kollman had showed that HG base pairing between adenine and thymine nucleobases (i.e., without sugar-phosphate backbone) in vacuum was energetically more favourable (∼1 kcal/mol) than its WC counterpart (20) in agreement with the structure in crystalline environment (3). Subsequent ab-initio quantum chemistry calculations for gas phase A-T base pairs have painted a more complex picture, showing that neither HG or WC forms are the global minimum (21). Furthermore, NMR studies have found that N-methylation of adenine favors syn HG base pairing over the anti WC form by preventing formation of a WC specific hydrogen bond by inducing steric clash (16). NMR chemical shift perturbation, NOESY experiments and constant pH MD simulation show that the HG base pairing among G-C base pairs is ∼20 times less abundant than that in A-T base pairs arguably because a protonation of cytosine is necessary for the transition to HG form (11, 22). Experiments of bare DNA duplex like A-T rich A6-DNA segment have shown presence of Hoogsteen base pair in the A16-T9 position from carbon R1ρ relaxation dispersion NMR signal (2). Such experiments on DNA bound to base intercalating drugs like echinomycin and actinomycin showed formation of Hoogsteen base pair (23). Since RNA duplexes also involve Watson-Crick like base pairs, it was natural to look for the existence of Hoogsteen transitions there, as well. However, HG base pairing was not observed in NMR relaxation experiment in a common A-form RNA hairpin with same sequence as A6-DNA. In contrast to DNA N-methylation of RNA adenine base to induce HG base pair formation resulted in base pair melting and consequent dissociation of the RNA helical region (24). In the present paper, we perform a computational comparison of the free energy psurface and transition paths of DNA vs RNA WC to HG transitions, seeking to understand the relative stability of the HG state and the factors that modulate the dynamical equilibrium in DNA vs RNA. The first report (2) of the NMR measured WC-HG dynamic equilibrium also included our preliminary work on searching pathways of conformational transition using a biased molecular dynamics simulation. Further computational work from our group and others (12, 23–26) have since uncovered several other pieces of information, including energetic details, transition pathways and kinetics of the transition between WC and HG base pairs in naked A-T rich DNA helical duplexes. The paths predicted in our initial study were also recently sampled by Vreede et al.(27), who have calculated the rate constants of back and forth transitions between WC and HG base pairs using Transition Path Sampling (TPS) (28). Experimental estimates put the WC-HG transition in DNA on a timescale of 50-250 ms (2). Therefore, direct simulation of Hoogsteen base pair formation in explicit water for a biologically relevant system is still beyond the reach of present day computational power. In previous work, we showed that formation of Hoogsteen base pair can be achieved in molecular dynamics simulation by applying an artificial biasing force on the the glycosidic torsion angle for free DNA (2, 24) and echinomycin-bound DNA (23). Enhanced sampling methods like umbrella sampling (29) and metadynamics (30, 31) can be utilized to obtain the free energy surface along slow coordinates. For example, Yang et. al. have calculated the free energy landscape of DNA breathing motion for a A-T rich duplex DNA segment from umbrella sampling (25) and multiple walker well tempered metadynamics simulation(26). Even for a short DNA segment, it required 6-40 µs long enhanced sampling simulations to obtain a converged physically meaningful free energy surface. Adaptive biasing force (ABF) (32) has recently become a good alternative to the traditional free energy sampling methods (33) for a wide range of biophysical problems including membrane permeation (34), ion transport (35), ligand diffusion in protein (36) and binding free energy calculation (37).
Recommended publications
  • 1 History for Canonical and Non-Canonical Structures of Nucleic Acids
    1 1 History for Canonical and Non-canonical Structures of Nucleic Acids The main points of the learning: Understand canonical and non-canonical structures of nucleic acids and think of historical scientists in the research field of nucleic acids. 1.1 Introduction This book is to interpret the non-canonical structures and their stabilities of nucleic acids from the viewpoint of the chemistry and study their biological significances. There is more than 60 years’ history after the discovery of the double helix DNA structure by James Dewey Watson and Francis Harry Compton Crick in 1953, and chemical biology of nucleic acids is facing a new aspect today. Through this book, I expect that readers understand how the uncommon structure of nucleic acids became one of the common structures that fascinate us now. In this chapter, I introduce the history of nucleic acid structures and the perspective of research for non-canonical nucleic acid structures (see also Chapter 15). 1.2 History of Duplex The opening of the history of genetics was mainly done by three researchers. Charles Robert Darwin, who was a scientist of natural science, pioneered genetics. The proposition of genetic concept is indicated in his book On the Origin of Species published in 1859. He indicated the theory of biological evolution, which is the basic scientific hypothesis of natural diversity. In other words, he proposed biological evolution, which changed among individuals by adapting to the environment and be passed on to the next generation. However, that was still a primitive idea for the genetic concept. After that, Gregor Johann Mendel, who was a priest in Brno, Czech Republic, confirmed the mechanism of gene evolution by using “factor” inherited from parent to children using pea plant in 1865.
    [Show full text]
  • Accurate Energies of Hydrogen Bonded Nucleic Acid Base Pairs and Triplets in Trna Tertiary Interactions Romina Oliva1,2,*, Luigi Cavallo3 and Anna Tramontano2
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by PubMed Central Nucleic Acids Research, 2006, Vol. 34, No. 3 865–879 doi:10.1093/nar/gkj491 Accurate energies of hydrogen bonded nucleic acid base pairs and triplets in tRNA tertiary interactions Romina Oliva1,2,*, Luigi Cavallo3 and Anna Tramontano2 1Centro Linceo Interdisciplinare ‘Beniamino Segre’, Accademia dei Lincei, I-00165 Rome, Italy, 2Dipartimento di Scienze Biochimiche ‘A. Rossi Fanelli’, Universita` di Roma ‘La Sapienza’, I-00185 Rome, Italy and 3Dipartimento di Chimica, Universita` di Salerno, Via Salvador Allende, Baronissi (SA), I-84081, Italy Received December 4, 2005; Revised and Accepted January 19, 2006 ABSTRACT intercalated A58–G18–A/G57–G19 purine bases. All these interactions are highly conserved in cytosolic tRNAs, indic- Tertiary interactions are crucial in maintaining the ating their importance for the tRNA function. Available tRNA structure and functionality. We used a com- experimental data also support the functional importance of bined sequence analysis and quantum mechanics such a region. Mutants of the Escherichia coli tRNAAla(CUA) approach to calculate accurate energies of the most obtained by random mutagenesis and having at least one of frequent tRNA tertiary base pairing interactions. Our the above-mentioned interactions disrupted were shown to be analysis indicates that six out of the nine classical non-functional (6), and insertion of additional nucleotides into tertiary interactions are held in place mainly by the T-loop or deletion of G19 from the D-loop affected the H-bonds between the bases. In the remaining three accuracy of the codon–anticodon recognition, resulting in a cases other effects have to be considered.
    [Show full text]
  • Plant G-Quadruplex (G4) Motifs in DNA and RNA; Abundant, Intriguing Sequences of Unknown Function T
    Plant Science 269 (2018) 143–147 Contents lists available at ScienceDirect Plant Science journal homepage: www.elsevier.com/locate/plantsci Review article Review: Plant G-quadruplex (G4) motifs in DNA and RNA; abundant, intriguing sequences of unknown function T ⁎ Brianna D. Griffin, Hank W. Bass Department of Biological Science, 319 Stadium Drive, Florida State University, Tallahassee, FL, 32306-4295, USA ARTICLE INFO ABSTRACT Keywords: DNA sequences capable of forming G-quadruplex (G4) structures can be predicted and mapped in plant genomes G-quadruplex using computerized pattern search programs. Non-telomeric G4 motifs have recently been found to number in G4 the thousands across many plant species and enriched around gene promoters, prompting speculation that they Gene regulation may represent a newly uncovered and ubiquitous family of cis-acting elements. Comparative analysis shows that cis-regulatory element monocots exhibit five to ten times higher G4 motif density than eudicots, but the significance of this difference has not been determined. The vast scale and complexity of G4 functions, actual or theoretical, are reviewed in relation to the multiple modes of action and myriad genetic functions for which G4s have been implicated in DNA and RNA. Future experimental strategies and opportunities include identifying plant G4-interactomes, resolving the structures of G4s with and without their binding partners, and defining molecular mechanisms through reporter gene, genetic, or genome editing approaches. Given the global importance of plants for food, clothing, medicine, and energy, together with the potential role of G4 motifs as a widely conserved set of DNA sequences that could coordinate gene regulation, future plant G4 research holds great potential for use in plant improvement strategies.
    [Show full text]
  • De Novo Nucleic Acids: a Review of Synthetic Alternatives to DNA and RNA That Could Act As † Bio-Information Storage Molecules
    life Review De Novo Nucleic Acids: A Review of Synthetic Alternatives to DNA and RNA That Could Act as y Bio-Information Storage Molecules Kevin G Devine 1 and Sohan Jheeta 2,* 1 School of Human Sciences, London Metropolitan University, 166-220 Holloway Rd, London N7 8BD, UK; [email protected] 2 Network of Researchers on the Chemical Evolution of Life (NoR CEL), Leeds LS7 3RB, UK * Correspondence: [email protected] This paper is dedicated to Professor Colin B Reese, Daniell Professor of Chemistry, Kings College London, y on the occasion of his 90th Birthday. Received: 17 November 2020; Accepted: 9 December 2020; Published: 11 December 2020 Abstract: Modern terran life uses several essential biopolymers like nucleic acids, proteins and polysaccharides. The nucleic acids, DNA and RNA are arguably life’s most important, acting as the stores and translators of genetic information contained in their base sequences, which ultimately manifest themselves in the amino acid sequences of proteins. But just what is it about their structures; an aromatic heterocyclic base appended to a (five-atom ring) sugar-phosphate backbone that enables them to carry out these functions with such high fidelity? In the past three decades, leading chemists have created in their laboratories synthetic analogues of nucleic acids which differ from their natural counterparts in three key areas as follows: (a) replacement of the phosphate moiety with an uncharged analogue, (b) replacement of the pentose sugars ribose and deoxyribose with alternative acyclic, pentose and hexose derivatives and, finally, (c) replacement of the two heterocyclic base pairs adenine/thymine and guanine/cytosine with non-standard analogues that obey the Watson–Crick pairing rules.
    [Show full text]
  • Conformation and Polymorphism in PNA-DNA and PNA-RNA Hybrids (Chimera/Duplex/Triple Helix/Molecular Mechanics)
    Proc. Natl. Acad. Sci. USA Vol. 90, pp. 9542-9546, October 1993 Biochemistry Peptide nucleic acid (PNA) conformation and polymorphism in PNA-DNA and PNA-RNA hybrids (chimera/duplex/triple helix/molecular mechanics) ORN ALMARSSON AND THOMAS C. BRUICE Department of Chemistry, University of California at Santa Barbara, Santa Barbara, CA 93106 Contributed by Thomas C. Bruice, July 19, 1993 ABSTRACT Two hydrogen-bonding motifs have been pro- nucleic acids, is their stability toward nucleases. Although the posed to account for the extraordinary stability of polyamide designing of the PNA analogues of DNA has been described "peptide" nucleic acid (PNA) hybrids with nucleic acids. These by the inventors (4), the design rationale does not explain the interresidue- and intraresidue-hydrogen-bond motifs were in- great stability of the PNA hybrids with nucleic acids, as vestigated by molecular mechanics calculations. Energy- exemplified by PNA-strand invasion of double-helical DNA. iinimized structures of Watson-Crick base-paired decameric Our recent molecular mechanics study (5) with the B-helix duplexes of PNA with A-, B-, and Z-DNA and A-RNA poly- structure of the DNA*PNA duplex (dAp)lo'(pnaT)lo [where morphs indicate that the inherent stability of the complemen- (dAp)10 is the decanucleotide of 2'-deoxyadenosine 5'- tary PNA helical structures is derived from interresidue, rather phosphate and (pnaT)lo is the decamer of thyminyl PNA with than from intraresidue, hydrogen bonds in all hybrids studied. lysine amide attached to the C terminus] strongly suggests that Intraresidue-hydrogen-bond lengths are consistently longer the stability of DNA*PNA duplexes is due to interresidue than interresidue hydrogen bonds.
    [Show full text]
  • Binary Complex Crystal Structure of DNA Polymerase Β Reveals Multiple Conformations of the Templating 8-Oxoguanine Lesion
    Binary complex crystal structure of DNA polymerase β reveals multiple conformations of the templating 8-oxoguanine lesion Vinod K. Batraa, David D. Shocka, William A. Bearda, Charles E. McKennab, and Samuel H. Wilsona,1 aLaboratory of Structural Biology, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709-2233; and bDepartment of Chemistry, University of Southern California, Los Angeles, CA 90089-0744 Edited by Philip C. Hanawalt, Stanford University, Stanford, CA, and approved October 31, 2011 (received for review July 27, 2011) Oxidation of genomic DNA forms the guanine lesion 7,8-dihydro-8- well-extended by DNA polymerases facilitating G to T mutagen- oxoguanine (8-oxoG). When in the template base position during esis. The structural basis for DNA polymerase preferences in DNA synthesis the 8-oxoG lesion has dual coding potential by vir- mutagenic bypass of the 8-oxoG lesion is emerging across the tue of its anti- and syn-conformations, base pairing with cytosine DNA polymerase families (8–12) and is the subject of the present and adenine, respectively. This impacts mutagenesis, because inser- report. tion of adenine opposite template 8-oxoG can result in a G to T Repair of the 8-oxoG-DNA lesion in mammalian cells pro- transversion. DNA polymerases vary by orders of magnitude in ceeds by several DNA repair pathways (13–15), but the major their preferences for mutagenic vs. error-free 8-oxoG lesion bypass. pathway is considered to be a form of the base excision repair Yet, the structural basis for lesion bypass specificity is not well un- (BER) pathway referred to as the “GO system” (16).
    [Show full text]
  • Structural Basis of Error-Prone Replication and Stalling at a Thymine Base by Human DNA Polymerase Iota Kevin N
    Western University Scholarship@Western Biochemistry Publications Biochemistry Department 6-3-2009 Structural Basis of Error-prone Replication and Stalling at a Thymine Base by Human DNA Polymerase iota Kevin N. Kirouac University of Western Ontario, [email protected] Hong Ling University of Western Ontario, [email protected] Follow this and additional works at: https://ir.lib.uwo.ca/biochempub Part of the Biochemistry Commons Citation of this paper: Kirouac, Kevin N. and Ling, Hong, "Structural Basis of Error-prone Replication and Stalling at a Thymine Base by Human DNA Polymerase iota" (2009). Biochemistry Publications. 66. https://ir.lib.uwo.ca/biochempub/66 Structural basis of error-prone replication and stalling at a thymine base by human DNA polymerase ι Kevin N Kirouac & Hong Ling* Department of Biochemistry, University of Western Ontario, London, Ontario, Canada, N6A 5C1 *Correspondence should be addressed to H. L . ( [email protected] ) Subject Category: Structural Biology; Genome Stability and Dynamics Characters (with spaces) 48 073 Figures 5 Table 1 Running Title: Polymerase ι misinserts and stalls at T bases Competing interest statement: The authors declare no competing financial interests 1 ABSTRACT Human DNA polymerase ι (pol ιιι) is a unique member of Y-family polymerases, which preferentially misincorporates nucleotides opposite thymines (T), and halts replication at T bases. The structural basis of the high error rates remains elusive. We present three crystal structures of pol ι complexed with DNA containing a thymine base, paired with correct or incorrect incoming nucleotides. A narrowed active site supports a pyrimidine:pyrimidine mismatch and excludes Watson-Crick base pairing in pol ιιι.
    [Show full text]
  • How Does Human 8-Oxoguanine DNA Glycosylase 1 (OGG1) Find 8-Oxoguanine in the Genome? Ostiane D’Augustin, Sébastien Huet, Anna Campalans, Juan Pablo Radicella
    Lost in the Crowd: How Does Human 8-Oxoguanine DNA Glycosylase 1 (OGG1) Find 8-Oxoguanine in the Genome? Ostiane D’augustin, Sébastien Huet, Anna Campalans, Juan Pablo Radicella To cite this version: Ostiane D’augustin, Sébastien Huet, Anna Campalans, Juan Pablo Radicella. Lost in the Crowd: How Does Human 8-Oxoguanine DNA Glycosylase 1 (OGG1) Find 8-Oxoguanine in the Genome?. International Journal of Molecular Sciences, MDPI, 2020, 21 (21), 10.3390/ijms21218360. hal- 03007117 HAL Id: hal-03007117 https://hal.archives-ouvertes.fr/hal-03007117 Submitted on 16 Nov 2020 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Distributed under a Creative Commons Attribution| 4.0 International License International Journal of Molecular Sciences Review Lost in the Crowd: How Does Human 8-Oxoguanine DNA Glycosylase 1 (OGG1) Find 8-Oxoguanine in the Genome? Ostiane D’Augustin 1,2 ,Sébastien Huet 2,3,* , Anna Campalans 1,* and Juan Pablo Radicella 1,* 1 Institute of Cellular and Molecular Radiobiology, Institut de Biologie François Jacob, CEA, Université Paris-Saclay,
    [Show full text]
  • Reverse Watson-Crick Base-Pairs in Theá D(GCGCGCG)/D(TCGCGCG)Á and D(GCGCGCG)/ D(CCGCGCG) Heteroduplexes, Respectively
    J. Mol. Biol. (1997) 269, 796±810 The Structures and Relative Stabilities of d(G G) Reverse Hoogsteen, d(G T) Reverse Wobble,Á and d(G C) Reverse Watson-CrickÁ Base-pairs in DNAÁ Crystals BlaineH.M.Mooers,BrandtF.EichmanandP.ShingHo* 5 Department of Biochemistry We have solved the structures of the homoduplex d(Gm CGCGCG)2, and and Biophysics, ALS 2011 the heteroduplexes d(GCGCGCG)/d(TCGCGCG) and d(GCGCGCG)/ Oregon State University d(CCGCGCG). The structures form six base-pairs of identical Z-DNA Corvallis, OR 97331, USA duplexes with single nucleotides overhanging at the 50-ends. The over- hanging nucleotide from one strand remains stacked and sandwiched between the blunt-ends of two adjacent Z-DNA duplexes, while the over- hanging base of the opposing strand is extra-helical. The stacked and the extra-helical bases from adjacent duplexes pair to form a distorted 5 d(G G) reverse Hoogsteen base-pair in the d(Gm CGCGCG)2 homo- duplex,Á and d(G T) reverse wobble and d(G C) reverse Watson-Crick base-pairs in theÁ d(GCGCGCG)/d(TCGCGCG)Á and d(GCGCGCG)/ d(CCGCGCG) heteroduplexes, respectively. Interestingly, only the d(G,T) and d(G C) base-pairs were observed in the heteroduplexes, suggesting that bothÁ the d(G T) reverse wobble and d(G C) reverse Watson-Crick base-pairs are moreÁ stable in this crystal environmentÁ than the d(G G) reverse Hoogsteen base-pair. To estimate the relative stability of the threeÁ types of reverse base-pairs, crystals were grown using various mixtures of sequences and their strand compositions analyzed by mass spec- trometry.
    [Show full text]
  • A Unified Computational View of DNA Duplex, Triplex, Quadruplex and Their
    Published online 24 April 2021 Nucleic Acids Research, 2021, Vol. 49, No. 9 4919–4933 doi: 10.1093/nar/gkab285 A unified computational view of DNA duplex, triplex, quadruplex and their donor–acceptor interactions Gyuri Park1, Byunghwa Kang1, Soyeon V. Park1, Donghwa Lee1,2,3,* and Seung Soo Oh 1,3,4,* 1Department of Materials Science and Engineering, Pohang University of Science and Technology (POSTECH), Pohang 37673, South Korea, 2Division of Advanced Materials Science, Pohang University of Science and Technology (POSTECH), Pohang 37673, South Korea, 3Institute for Convergence Research and Education in Advanced Technology (I-CREATE), Yonsei University, Incheon 21983, South Korea and 4School of Interdisciplinary Bioscience and Bioengineering, Pohang University of Science and Technology (POSTECH), Pohang 37673, South Korea Received October 15, 2020; Revised April 07, 2021; Editorial Decision April 07, 2021; Accepted April 14, 2021 ABSTRACT INTRODUCTION DNA can assume various structures as a result Deoxyribonucleic acid (DNA) is a unique material as com- of interactions at atomic and molecular levels posed of nitrogenous bases (adenine (A), thymine (T), (e.g., hydrogen bonds, ␲–␲ stacking interactions, guanine (G) or cytosine (C)), sugar rings, and phosphate and electrostatic potentials), so understanding of groups. DNA is programmable, so it can be rationally the consequences of these interactions could guide designed into molecular structures ranging from simple Watson–Crick base-pairing primers to DNA origami-based development of ways to produce elaborate pro- complex 3D constructs (1,2). Even spatial and temporal grammable DNA for applications in bio- and nan- control of DNA nanostructures is achievable; sophisticated otechnology. We conducted advanced ab initio cal- DNA molecular machines can perform a series of nanome- culations to investigate nucleobase model struc- chanical motions in a controllable manner, so this abil- tures by componentizing their donor-acceptor inter- ity provides unprecedented applications in bio- and nan- actions.
    [Show full text]
  • A DNA G-Quadruplex/I-Motif Hybrid
    bioRxiv preprint doi: https://doi.org/10.1101/737296; this version posted August 15, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. A DNA G-quadruplex/i-motif hybrid Betty Chu, Daoning Zhang, and Paul J. Paukstelis* Department of Chemistry and Biochemistry, Center for Biomolecular Structure and Organization, University of Maryland, College Park, Maryland, 20742, United States Correspondence: Paul J. Paukstelis [email protected] 1 bioRxiv preprint doi: https://doi.org/10.1101/737296; this version posted August 15, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Abstract DNA can form many structures beyond the canonical Watson-Crick double helix. It is now clear that noncanonical structures are present in genomic DNA and have biological functions. G-rich G-quadruplexes and C-rich i-motifs are the most well-characterized noncanonical DNA motifs that have been detected in vivo with either proscribed or postulated biological roles. Because of their independent sequence requirements, these structures have largely been considered distinct types of quadruplexes. Here, we describe the crystal structure of the DNA oligonucleotide, d(CCAGGCTGCAA), that self-associates to form a quadruplex structure containing two central antiparallel G-tetrads and six i-motif C-C+ base pairs. Solution studies suggest a robust structural motif capable of assembling as a tetramer of individual strands or as a dimer when composed of tandem repeats. This hybrid structure highlights the growing structural diversity of DNA and suggests that biological systems may harbor many functionally important non-duplex structures.
    [Show full text]
  • Canonical and Non-Canonical Base Pairs in DNA Or RNA: Structure, Function and Dynamics
    Canonical and Non-canonical Base Pairs in DNA or RNA: Structure, Function and Dynamics Dhananjay Bhattacharyya1 and Abhijit Mitra2 1 Computational Science Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata 700064, INDIA 2 Center for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology (IIIT-H), Gachibowli, Hyderabad 500032, INDIA ABSTRACT: Non-canonical base pairs are planar hydrogen bonded pairs of nucleobases, having hydrogen bonding patterns which differ from the patterns observed in Watson-Crick base pairs, as in the classic double helical DNA. The structures of polynucleotide strands of both DNA and RNA molecules can be understood in terms of sugar-phosphate backbones consisting of phosphodiester-linked D-2‘-deoxyribofuranose (D-ribofuranose in RNA) sugar moieties, with purine or pyrimidine nucleobases covalently linked to them. Here, the N9 atoms of the purines, guanine and adenine, and the N1 atoms of the pyrimidines, cytosine and thymine (uracil in RNA), respectively, form glycosidic linkages with the C1‘ atom of the sugars. These nucleobases can be schematically represented as triangles with one of their vertices linked to the sugar, and the three sides accounting for three edges through which they can form hydrogen bonds with other moieties, including with other nucleobases. As also explained in greater details later in this article, the side opposite to the sugar linked vertex is traditionally called the Watson-Crick edge, since they are involved in forming the Watson- Crick base pairs which constitute building blocks of double helical DNA. The two sides adjacent to the sugar-linked vertex are referred to, respectively, as the sugar and Hoogsteen (C-H for pyrimidines) edges.
    [Show full text]