51
|
Šponer J, Bussi G, Stadlbauer P, Kührová P, Banáš P, Islam B, Haider S, Neidle S, Otyepka M. Folding of guanine quadruplex molecules-funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies. Biochim Biophys Acta Gen Subj 2016; 1861:1246-1263. [PMID: 27979677 DOI: 10.1016/j.bbagen.2016.12.008] [Citation(s) in RCA: 74] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2016] [Revised: 12/04/2016] [Accepted: 12/11/2016] [Indexed: 01/18/2023]
Abstract
BACKGROUND Guanine quadruplexes (GQs) play vital roles in many cellular processes and are of much interest as drug targets. In contrast to the availability of many structural studies, there is still limited knowledge on GQ folding. SCOPE OF REVIEW We review recent molecular dynamics (MD) simulation studies of the folding of GQs, with an emphasis paid to the human telomeric DNA GQ. We explain the basic principles and limitations of all types of MD methods used to study unfolding and folding in a way accessible to non-specialists. We discuss the potential role of G-hairpin, G-triplex and alternative GQ intermediates in the folding process. We argue that, in general, folding of GQs is fundamentally different from funneled folding of small fast-folding proteins, and can be best described by a kinetic partitioning (KP) mechanism. KP is a competition between at least two (but often many) well-separated and structurally different conformational ensembles. MAJOR CONCLUSIONS The KP mechanism is the only plausible way to explain experiments reporting long time-scales of GQ folding and the existence of long-lived sub-states. A significant part of the natural partitioning of the free energy landscape of GQs comes from the ability of the GQ-forming sequences to populate a large number of syn-anti patterns in their G-tracts. The extreme complexity of the KP of GQs typically prevents an appropriate description of the folding landscape using just a few order parameters or collective variables. GENERAL SIGNIFICANCE We reconcile available computational and experimental studies of GQ folding and formulate basic principles characterizing GQ folding landscapes. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic; Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic.
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy
| | - Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic; Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Barira Islam
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Shozeb Haider
- UCL School of Pharmacy, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Stephen Neidle
- UCL School of Pharmacy, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
52
|
Karg B, Haase L, Funke A, Dickerhoff J, Weisz K. Observation of a Dynamic G-Tetrad Flip in Intramolecular G-Quadruplexes. Biochemistry 2016; 55:6949-6955. [PMID: 27951645 DOI: 10.1021/acs.biochem.6b00925] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A MYC sequence forming an intramolecular G-quadruplex with a parallel topology was modified by the incorporation of 8-bromoguanosine (BrG) analogues in one of its outer G-tetrads. The propensity of the BrG analogues to adopt a syn glycosidic torsion angle results in an exceptional monomolecular quadruplex conformation featuring a complete flip of one tetrad while keeping a parallel orientation of all G-tracts as shown by circular dichroism and nuclear magnetic resonance spectroscopic studies. When substituting three of the four G-tetrad residues with BrG analogues, two coexisting quadruplex conformational isomers with an all-syn and all-anti outer G-quartet are approximately equally populated in solution. A dynamic interconversion of the two quadruplexes with an exchange rate (kex) of 0.2 s-1 is demonstrated through the observation of exchange crosspeaks in rotating frame Overhauser effect spectroscopy and nuclear Overhauser effect spectroscopy experiments at 50 °C. The kinetic properties suggest disruption of the corresponding outer G-tetrad but not of the whole quadruplex core during the tetrad flip. Conformational syn-anti isomers with homopolar and heteropolar stacking interactions are nearly isoenergetic with a transition enthalpy of 18.2 kJ/mol in favor of the all-syn isomer.
Collapse
Affiliation(s)
- Beatrice Karg
- Institute of Biochemistry, Ernst-Moritz-Arndt University Greifswald , Felix-Hausdorff-Strasse 4, D-17487 Greifswald, Germany
| | - Linn Haase
- Institute of Biochemistry, Ernst-Moritz-Arndt University Greifswald , Felix-Hausdorff-Strasse 4, D-17487 Greifswald, Germany
| | - Andrea Funke
- Institute of Biochemistry, Ernst-Moritz-Arndt University Greifswald , Felix-Hausdorff-Strasse 4, D-17487 Greifswald, Germany
| | - Jonathan Dickerhoff
- Institute of Biochemistry, Ernst-Moritz-Arndt University Greifswald , Felix-Hausdorff-Strasse 4, D-17487 Greifswald, Germany
| | - Klaus Weisz
- Institute of Biochemistry, Ernst-Moritz-Arndt University Greifswald , Felix-Hausdorff-Strasse 4, D-17487 Greifswald, Germany
| |
Collapse
|
53
|
ALS and FTD linked GGGGCC-repeat containing DNA oligonucleotide folds into two distinct G-quadruplexes. Biochim Biophys Acta Gen Subj 2016; 1861:1237-1245. [PMID: 27856299 DOI: 10.1016/j.bbagen.2016.11.018] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2016] [Revised: 11/11/2016] [Accepted: 11/13/2016] [Indexed: 01/24/2023]
Abstract
BACKGROUND The most common genetic cause of neurological disorders ALS and FTD is a largely increased number of GGGGCC repeats in C9orf72 gene. Non-canonical structures including G-quadruplexes adopted by expanded repeats are hypothesized to be crucial in pathogenesis. Recently, we have shown that structural polymorphism of oligonucleotide d(G4C2)3G4 is reduced by dG to 8Br-dG substitution. High-resolution structure of one of the two major G-quadruplexes adopts antiparallel topology comprising of four G-quartets. Herein, we describe the topology of the second major G-quadruplex structure and influence of folding conditions on relative populations of the two folds. METHODS Influence of folding conditions was explored by 1H 1D NMR. Determination of topology was achieved by 2D NMR complemented with PAGE and CD. UV melting experiment was used to explore thermal stability of structures. RESULTS Two structures adopted by oligonucleotide d(G4C2)3GGBrGG denoted AQU and NAN coexist in solution and ratio of their populations is determined by pH and rate of cooling when folding from thermally denatured state in the presence of K+ ions. CONCLUSIONS AQU is kinetically favored and forms by folding at low pH, while NAN is favored thermodynamically and at neutral pH. AQU and NAN share similar antiparallel topology with four G-quartets and three edgewise loops, however they exhibit distinct structural and dynamic properties. GENERAL SIGNIFICANCE Novel G-quadruplex topology adds insight into diverse polymorphism of DNA sequences comprising potentially pathological GGGGCC repeat. Relative populations of the two structures and their dependence on folding conditions contribute to understanding of factors that govern G-quadruplex folding. This article is part of a Special Issue entitled "Gquadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.
Collapse
|
54
|
Stadlbauer P, Mazzanti L, Cragnolini T, Wales DJ, Derreumaux P, Pasquali S, Šponer J. Coarse-Grained Simulations Complemented by Atomistic Molecular Dynamics Provide New Insights into Folding and Unfolding of Human Telomeric G-Quadruplexes. J Chem Theory Comput 2016; 12:6077-6097. [DOI: 10.1021/acs.jctc.6b00667] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Petr Stadlbauer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
- Regional
Centre of Advanced Technologies and Materials, Departments of Physical
Chemistry, Faculty of Science, Palacký University, 17. listopadu
1192/12, 771 46 Olomouc, Czech Republic
| | - Liuba Mazzanti
- Laboratoire
de Biochimie Théorique, IBPC, CNRS UPR9080, Université Sorbonne Paris Cite, Paris Diderot, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Tristan Cragnolini
- Department
of Chemistry, Cambridge University, Lensfield Road, Cambridge CB2 1EW, U.K
| | - David J. Wales
- Department
of Chemistry, Cambridge University, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Philippe Derreumaux
- Laboratoire
de Biochimie Théorique, IBPC, CNRS UPR9080, Université Sorbonne Paris Cite, Paris Diderot, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Samuela Pasquali
- Laboratoire
de Biochimie Théorique, IBPC, CNRS UPR9080, Université Sorbonne Paris Cite, Paris Diderot, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Jiří Šponer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| |
Collapse
|
55
|
Kogut M, Kleist C, Czub J. Molecular dynamics simulations reveal the balance of forces governing the formation of a guanine tetrad-a common structural unit of G-quadruplex DNA. Nucleic Acids Res 2016; 44:3020-30. [PMID: 26980278 PMCID: PMC4838382 DOI: 10.1093/nar/gkw160] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2016] [Accepted: 03/01/2016] [Indexed: 12/13/2022] Open
Abstract
G-quadruplexes (G4) are nucleic acid conformations of guanine-rich sequences, in which guanines are arranged in the square-planar G-tetrads, stacked on one another. G4 motifs form in vivo and are implicated in regulation of such processes as gene expression and chromosome maintenance. The structure and stability of various G4 topologies were determined experimentally; however, the driving forces for their formation are not fully understood at the molecular level. Here, we used all-atom molecular dynamics to probe the microscopic origin of the G4 motif stability. By computing the free energy profiles governing the dissociation of the 3′-terminal G-tetrad in the telomeric parallel-stranded G4, we examined the thermodynamic and kinetic stability of a single G-tetrad, as a common structural unit of G4 DNA. Our results indicate that the energetics of guanine association alone does not explain the overall stability of the G-tetrad and that interactions involving sugar–phosphate backbone, in particular, the constrained minimization of the phosphate–phosphate repulsion energy, are crucial in providing the observed enthalpic stabilization. This enthalpic gain is largely compensated by the unfavorable entropy change due to guanine association and optimization of the backbone topology.
Collapse
Affiliation(s)
- Mateusz Kogut
- Department of Physical Chemistry, Gdansk University of Technology, ul. Narutowicza 11/12, 80-233 Gdansk, Poland
| | - Cyprian Kleist
- Department of Physical Chemistry, Gdansk University of Technology, ul. Narutowicza 11/12, 80-233 Gdansk, Poland
| | - Jacek Czub
- Department of Physical Chemistry, Gdansk University of Technology, ul. Narutowicza 11/12, 80-233 Gdansk, Poland
| |
Collapse
|
56
|
Islam B, Stadlbauer P, Neidle S, Haider S, Sponer J. Can We Execute Reliable MM-PBSA Free Energy Computations of Relative Stabilities of Different Guanine Quadruplex Folds? J Phys Chem B 2016; 120:2899-912. [DOI: 10.1021/acs.jpcb.6b01059] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Barira Islam
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Stephen Neidle
- UCL School of Pharmacy, 29-39
Brunswick Square, London WC1N 1AX, U.K
| | - Shozeb Haider
- UCL School of Pharmacy, 29-39
Brunswick Square, London WC1N 1AX, U.K
| | - Jiri Sponer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
- CEITEC
- Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
57
|
Gkionis K, Kruse H, Šponer J. Derivation of Reliable Geometries in QM Calculations of DNA Structures: Explicit Solvent QM/MM and Restrained Implicit Solvent QM Optimizations of G-Quadruplexes. J Chem Theory Comput 2016; 12:2000-16. [DOI: 10.1021/acs.jctc.5b01025] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Konstantinos Gkionis
- Institute of Biophysics,
Academy of Sciences of the Czech Republic, Královopolská 135, 612
65 Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics,
Academy of Sciences of the Czech Republic, Královopolská 135, 612
65 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics,
Academy of Sciences of the Czech Republic, Královopolská 135, 612
65 Brno, Czech Republic
- CEITEC
− Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
58
|
Largy E, Marchand A, Amrane S, Gabelica V, Mergny JL. Quadruplex Turncoats: Cation-Dependent Folding and Stability of Quadruplex-DNA Double Switches. J Am Chem Soc 2016; 138:2780-92. [PMID: 26837276 DOI: 10.1021/jacs.5b13130] [Citation(s) in RCA: 103] [Impact Index Per Article: 12.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Quadruplex (G4) nucleic acids, a family of secondary structures formed by guanine-rich sequences, exhibit an important structural polymorphism. We demonstrate here that G-rich DNA sequences may function as a double switch based on different triggers, provided that their quadruplex structures and stability display a high dependence on cation nature and concentration. A first switch is based on a remarkable antiparallel-to-parallel conversion, taking place in a few seconds at room temperature by addition of low KCl amounts to a sodium-rich sample. The second switch involves the conversion of alternative antiparallel quadruplex structures binding only one cation, formed in the presence of sub-millimolar potassium or strontium concentrations, to parallel structures by increasing the cation concentration. Incidentally, extremely low K(+) or Sr(2+) concentrations (≤5 equiv) are sufficient to induce G4 formation in a buffer devoid of other G4-promoting cations, and we suggest that the alternative structures observed contain only two tetrads. Such DNA systems are biological relevant targets, can be used in nanotechnology applications, and are valuable methodological tools for understanding DNA quadruplex folding, notably at low cation concentrations. We demonstrate that this behavior is not restricted to a narrow set of sequences but can also be found for other G-quadruplex-forming motifs, arguing for widespread applications.
Collapse
Affiliation(s)
- Eric Largy
- U1212, ARNA Laboratory, Inserm , F-33000 Bordeaux, France.,IECB, ARNA Laboratory, Université de Bordeaux , F-33600 Pessac, France.,UMR 5320, ARNA Laboratory, CNRS , F-33600 Pessac, France
| | - Adrien Marchand
- U1212, ARNA Laboratory, Inserm , F-33000 Bordeaux, France.,IECB, ARNA Laboratory, Université de Bordeaux , F-33600 Pessac, France.,UMR 5320, ARNA Laboratory, CNRS , F-33600 Pessac, France
| | - Samir Amrane
- U1212, ARNA Laboratory, Inserm , F-33000 Bordeaux, France.,IECB, ARNA Laboratory, Université de Bordeaux , F-33600 Pessac, France.,UMR 5320, ARNA Laboratory, CNRS , F-33600 Pessac, France
| | - Valérie Gabelica
- U1212, ARNA Laboratory, Inserm , F-33000 Bordeaux, France.,IECB, ARNA Laboratory, Université de Bordeaux , F-33600 Pessac, France.,UMR 5320, ARNA Laboratory, CNRS , F-33600 Pessac, France
| | - Jean-Louis Mergny
- U1212, ARNA Laboratory, Inserm , F-33000 Bordeaux, France.,IECB, ARNA Laboratory, Université de Bordeaux , F-33600 Pessac, France.,UMR 5320, ARNA Laboratory, CNRS , F-33600 Pessac, France
| |
Collapse
|
59
|
Largy E, Mergny JL, Gabelica V. Role of Alkali Metal Ions in G-Quadruplex Nucleic Acid Structure and Stability. Met Ions Life Sci 2016; 16:203-58. [PMID: 26860303 DOI: 10.1007/978-3-319-21756-7_7] [Citation(s) in RCA: 105] [Impact Index Per Article: 13.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
G-quadruplexes are guanine-rich nucleic acids that fold by forming successive quartets of guanines (the G-tetrads), stabilized by intra-quartet hydrogen bonds, inter-quartet stacking, and cation coordination. This specific although highly polymorphic type of secondary structure deviates significantly from the classical B-DNA duplex. G-quadruplexes are detectable in human cells and are strongly suspected to be involved in a number of biological processes at the DNA and RNA levels. The vast structural polymorphism exhibited by G-quadruplexes, together with their putative biological relevance, makes them attractive therapeutic targets compared to canonical duplex DNA. This chapter focuses on the essential and specific coordination of alkali metal cations by G-quadruplex nucleic acids, and most notably on studies highlighting cation-dependent dissimilarities in their stability, structure, formation, and interconversion. Section 1 surveys G-quadruplex structures and their interactions with alkali metal ions while Section 2 presents analytical methods used to study G-quadruplexes. The influence of alkali cations on the stability, structure, and kinetics of formation of G-quadruplex structures of quadruplexes will be discussed in Sections 3 and 4. Section 5 focuses on the cation-induced interconversion of G-quadruplex structures. In Sections 3 to 5, we will particularly emphasize the comparisons between cations, most often K(+) and Na(+) because of their prevalence in the literature and in cells.
Collapse
Affiliation(s)
- Eric Largy
- ARNA Laboratory, Université Bordeaux, IECB, 2, rue Robert Escarpit, F-33600, Pessac, France.,ARNA Laboratory, INSERM, U869, F-33000, Bordeaux, France
| | - Jean-Louis Mergny
- ARNA Laboratory, Université Bordeaux, IECB, 2, rue Robert Escarpit, F-33600, Pessac, France. .,ARNA Laboratory, INSERM, U869, F-33000, Bordeaux, France.
| | - Valérie Gabelica
- ARNA Laboratory, Université Bordeaux, IECB, 2, rue Robert Escarpit, F-33600, Pessac, France. .,ARNA Laboratory, INSERM, U869, F-33000, Bordeaux, France.
| |
Collapse
|
60
|
Zgarbová M, Šponer J, Otyepka M, Cheatham TE, Galindo-Murillo R, Jurečka P. Refinement of the Sugar-Phosphate Backbone Torsion Beta for AMBER Force Fields Improves the Description of Z- and B-DNA. J Chem Theory Comput 2015; 11:5723-36. [PMID: 26588601 DOI: 10.1021/acs.jctc.5b00716] [Citation(s) in RCA: 336] [Impact Index Per Article: 37.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Z-DNA duplexes are a particularly complicated test case for current force fields. We performed a set of explicit solvent molecular dynamics (MD) simulations with various AMBER force field parametrizations including our recent refinements of the ε/ζ and glycosidic torsions. None of these force fields described the ZI/ZII and other backbone substates correctly, and all of them underpredicted the population of the important ZI substate. We show that this underprediction can be attributed to an inaccurate potential for the sugar-phosphate backbone torsion angle β. We suggest a refinement of this potential, β(OL1), which was derived using our recently introduced methodology that includes conformation-dependent solvation effects. The new potential significantly increases the stability of the dominant ZI backbone substate and improves the overall description of the Z-DNA backbone. It also has a positive (albeit small) impact on another important DNA form, the antiparallel guanine quadruplex (G-DNA), and improves the description of the canonical B-DNA backbone by increasing the population of BII backbone substates, providing a better agreement with experiment. We recommend using β(OL1) in combination with our previously introduced corrections, εζ(OL1) and χ(OL4), (the combination being named OL15) as a possible alternative to the current β torsion potential for more accurate modeling of nucleic acids.
Collapse
Affiliation(s)
- Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic.,Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, University of Utah , 30 South 2000 East, Skaggs 105, Salt Lake City, Utah 84112, United States
| | - Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, University of Utah , 30 South 2000 East, Skaggs 105, Salt Lake City, Utah 84112, United States
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic
| |
Collapse
|
61
|
Stadlbauer P, Kührová P, Banáš P, Koča J, Bussi G, Trantírek L, Otyepka M, Šponer J. Hairpins participating in folding of human telomeric sequence quadruplexes studied by standard and T-REMD simulations. Nucleic Acids Res 2015; 43:9626-44. [PMID: 26433223 PMCID: PMC4787745 DOI: 10.1093/nar/gkv994] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2015] [Accepted: 09/22/2015] [Indexed: 11/13/2022] Open
Abstract
DNA G-hairpins are potential key structures participating in folding of human telomeric guanine quadruplexes (GQ). We examined their properties by standard MD simulations starting from the folded state and long T-REMD starting from the unfolded state, accumulating ∼130 μs of atomistic simulations. Antiparallel G-hairpins should spontaneously form in all stages of the folding to support lateral and diagonal loops, with sub-μs scale rearrangements between them. We found no clear predisposition for direct folding into specific GQ topologies with specific syn/anti patterns. Our key prediction stemming from the T-REMD is that an ideal unfolded ensemble of the full GQ sequence populates all 4096 syn/anti combinations of its four G-stretches. The simulations can propose idealized folding pathways but we explain that such few-state pathways may be misleading. In the context of the available experimental data, the simulations strongly suggest that the GQ folding could be best understood by the kinetic partitioning mechanism with a set of deep competing minima on the folding landscape, with only a small fraction of molecules directly folding to the native fold. The landscape should further include non-specific collapse processes where the molecules move via diffusion and consecutive random rare transitions, which could, e.g. structure the propeller loops.
Collapse
Affiliation(s)
- Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jaroslav Koča
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic National Center for Biomolecular Research, Faculty of Science, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy
| | - Lukáš Trantírek
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
62
|
Marušič M, Plavec J. The Effect of DNA Sequence Directionality on G-Quadruplex Folding. Angew Chem Int Ed Engl 2015. [DOI: 10.1002/ange.201505348] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
|
63
|
Marušič M, Plavec J. The Effect of DNA Sequence Directionality on G‐Quadruplex Folding. Angew Chem Int Ed Engl 2015; 54:11716-9. [DOI: 10.1002/anie.201505348] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2015] [Revised: 07/20/2015] [Indexed: 11/07/2022]
Affiliation(s)
- Maja Marušič
- Slovenian NMR Center, National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana (Slovenia)
| | - Janez Plavec
- Slovenian NMR Center, National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana (Slovenia)
- EN‐FIST Center of Excellence, 1000 Ljubljana (Slovenia)
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, 1000 Ljubljana (Slovenia)
| |
Collapse
|
64
|
Islam B, Stadlbauer P, Krepl M, Koca J, Neidle S, Haider S, Sponer J. Extended molecular dynamics of a c-kit promoter quadruplex. Nucleic Acids Res 2015; 43:8673-93. [PMID: 26245347 PMCID: PMC4605300 DOI: 10.1093/nar/gkv785] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2015] [Accepted: 07/21/2015] [Indexed: 01/29/2023] Open
Abstract
The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.
Collapse
Affiliation(s)
- Barira Islam
- Central European Institute of Technology (CEITEC), Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Jaroslav Koca
- Central European Institute of Technology (CEITEC), Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic National Center for Biomolecular Research, Faculty of Science, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Stephen Neidle
- UCL School of Pharmacy, University College London, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Shozeb Haider
- UCL School of Pharmacy, University College London, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Jiri Sponer
- Central European Institute of Technology (CEITEC), Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| |
Collapse
|
65
|
Tsvetkov V, Pozmogova G, Varizhuk A. The systematic approach to describing conformational rearrangements in G-quadruplexes. J Biomol Struct Dyn 2015; 34:705-15. [PMID: 26017012 PMCID: PMC4867883 DOI: 10.1080/07391102.2015.1055303] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Conformational changes in DNA G-quadruplex (GQ)-forming regions affect genome function and, thus, compose an interesting research topic. Computer modelling may yield insight into quadruplex folding and rearrangement, particularly molecular dynamics simulations. Here, we show that specific parameters, which are distinct from those commonly used in DNA conformational analyses, must be introduced for adequate interpretation and, most importantly, convenient visual representation of the quadruplex modelling results. We report a set of parameters that comprehensively and systematically describe GQ geometry in dynamics. The parameters include those related to quartet planarity, quadruplex twist, and quartet stacking; they are used to quantitatively characterise various types of quadruplexes and rearrangements, such as quartet distortion/disruption or deviation/bulging of a single nucleotide from the quartet plane. Our approach to describing conformational changes in quadruplexes using the new parameters is exemplified by telomeric quadruplex rearrangement, and the benefits of applying this approach to analyse other structures are discussed.
Collapse
Affiliation(s)
- Vladimir Tsvetkov
- a Department of Molecular Biology and Genetics , SRI of Physical-Chemical Medicine , Moscow , 119435 , Russia.,b Department of Polyelectrolytes and Surface-active Polymers , Topchiev Institute of Petrochemical Synthesis , Moscow , 119991 , Russia
| | - Galina Pozmogova
- a Department of Molecular Biology and Genetics , SRI of Physical-Chemical Medicine , Moscow , 119435 , Russia
| | - Anna Varizhuk
- a Department of Molecular Biology and Genetics , SRI of Physical-Chemical Medicine , Moscow , 119435 , Russia.,c Department of Structure-Functional Analysis of Biopolymers , Engelhardt Institute of Molecular Biology , Vavilov str. 32, Moscow , 119991 , Russia
| |
Collapse
|
66
|
Collie GW, Campbell NH, Neidle S. Loop flexibility in human telomeric quadruplex small-molecule complexes. Nucleic Acids Res 2015; 43:4785-99. [PMID: 25940631 PMCID: PMC4446451 DOI: 10.1093/nar/gkv427] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2015] [Accepted: 04/20/2015] [Indexed: 01/25/2023] Open
Abstract
Quadruplex nucleic acids can be formed at the ends of eukaryotic chromosomes. Their formation and stabilisation by appropriate small molecules can be used as a means of inhibiting the telomere maintenance functions of telomerase in human cancer cells. The crystal structures have been determined for a number of complexes between these small molecules and human telomeric DNA and RNA quadruplexes. The detailed structural characteristics of these complexes have been surveyed here and the variations in conformation for the TTA and UUA loops have been explored. Loop conformations have been classified in terms of a number of discrete types and their distribution among the crystal structures. Sugar conformation and backbone angles have also been examined and trends highlighted. One particular loop class has been found to be most prevalent. Implications for in particular, rational drug design, are discussed.
Collapse
Affiliation(s)
- Gavin W Collie
- UCL School of Pharmacy, University College London, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Nancy H Campbell
- UCL School of Pharmacy, University College London, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Stephen Neidle
- UCL School of Pharmacy, University College London, 29-39 Brunswick Square, London WC1N 1AX, UK
| |
Collapse
|
67
|
Dickerhoff J, Weisz K. Flipping a G-tetrad in a unimolecular quadruplex without affecting its global fold. Angew Chem Int Ed Engl 2015; 54:5588-91. [PMID: 25775974 DOI: 10.1002/anie.201411887] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2014] [Revised: 02/22/2015] [Indexed: 11/07/2022]
Abstract
A unimolecular G-quadruplex with a hybrid-type topology and propeller, diagonal, and lateral loops was examined for its ability to undergo structural changes upon specific modifications. Substituting 2'-deoxy-2'-fluoro analogues with a propensity to adopt an anti glycosidic conformation for two or three guanine deoxyribonucleosides in syn positions of the 5'-terminal G-tetrad significantly alters the CD spectral signature of the quadruplex. An NMR analysis reveals a polarity switch of the whole tetrad with glycosidic conformational changes detected for all four guanine nucleosides in the modified sequence. As no additional rearrangement of the overall fold occurs, a novel type of G-quadruplex is formed with guanosines in the four columnar G-tracts lined up in either an all-syn or an all-anti glycosidic conformation.
Collapse
Affiliation(s)
- Jonathan Dickerhoff
- Institut für Biochemie, Ernst-Moritz-Arndt-Universität Greifswald, Felix-Hausdorff-Strasse 4, 17487 Greifswald (Germany)
| | | |
Collapse
|
68
|
Dickerhoff J, Weisz K. Flipping a G-Tetrad in a Unimolecular Quadruplex Without Affecting Its Global Fold. Angew Chem Int Ed Engl 2015. [DOI: 10.1002/ange.201411887] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
69
|
Ghosh S, Jana J, Kar RK, Chatterjee S, Dasgupta D. Plant alkaloid chelerythrine induced aggregation of human telomere sequence--a unique mode of association between a small molecule and a quadruplex. Biochemistry 2015; 54:974-86. [PMID: 25566806 DOI: 10.1021/bi501117x] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Small molecules that interact with G-quadruplex structures formed by the human telomeric region and stabilize them have the potential to evolve as anticancer therapeutic agents. Herein we report the interaction of a putative anticancer agent from a plant source, chelerythrine, with the human telomeric DNA sequence. It has telomerase inhibitory potential as demonstrated from telomerase repeat amplification assay in cancer cell line extract. We have attributed this to the quadruplex binding potential of the molecule and characterized the molecular details of the interaction by means of optical spectroscopy such as absorbance and circular dichroism and calorimetric techniques such as isothermal titration calorimetry and differential scanning calorimetry. The results show that chelerythrine binds with micromolar dissociation constant and 2:1 binding stoichiometry to the human telomeric DNA sequence. Chelerythrine association stabilizes the G-quadruplex. Nuclear magnetic resonance spectroscopy ((1)H and (31)P) shows that chelerythrine binds to both G-quartet and phosphate backbone of the quadruplex leading to quadruplex aggregation. Molecular dynamics simulation studies support the above inferences and provide further insight into the mechanism of ligand binding. The specificity toward quartet binding for chelerythrine is higher compared to that of groove binding. MM-PBSA calculation mines out the energy penalty for quartet binding to be -4.7 kcal/mol, whereas that of the groove binding is -1.7 kcal/mol. We propose that the first chelerythrine molecule binds to the quartet followed by a second molecule which binds to the groove. This second molecule might bring about aggregation of the quadruplex structure which is evident from the results of nuclear magnetic resonance.
Collapse
Affiliation(s)
- Saptaparni Ghosh
- Biophysics & Structural Genomics Division, Saha Institute of Nuclear Physics , Block-AF, Sector-I, Bidhannagar, Kolkata-700064, India
| | | | | | | | | |
Collapse
|
70
|
Wei D, Husby J, Neidle S. Flexibility and structural conservation in a c-KIT G-quadruplex. Nucleic Acids Res 2014; 43:629-44. [PMID: 25452341 PMCID: PMC4288176 DOI: 10.1093/nar/gku1282] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
A quadruplex sequence from the promoter region of the c-KIT gene forms a stable quadruplex, as characterized by crystallographic and NMR methods. Two new crystal structures are reported here, together with molecular dynamics simulation studies on these quadruplex crystal structures and an NMR structure. The new crystal structures, each in a distinct space group and lattice packing arrangement, together with the existing structures, demonstrate that the c-KIT quadruplex fold does not change with differing environments, suggesting that quadruplex topological dynamism is not a general phenomenon. The single and dinucleotide loops in these structures show a high degree of conformational flexibility within the three crystal forms and the NMR ensemble, with no evidence of clustering to particular conformers. This is in accord with the findings of high loop flexibility from the molecular dynamics studies. It is suggested that intramolecular quadruplexes can be grouped into two broad classes (i) those with at least one single-nucleotide loop, often showing singular topologies even though loops are highly flexible, and (ii) with all loops comprising at least two nucleotides, leading to topological dynamism. The loops can have more stable and less dynamic base-stacked secondary structures.
Collapse
Affiliation(s)
- Dengguo Wei
- UCL School of Pharmacy, University College London, London WC1N 1AX, UK
| | - Jarmila Husby
- UCL School of Pharmacy, University College London, London WC1N 1AX, UK
| | - Stephen Neidle
- UCL School of Pharmacy, University College London, London WC1N 1AX, UK
| |
Collapse
|
71
|
Triplex intermediates in folding of human telomeric quadruplexes probed by microsecond-scale molecular dynamics simulations. Biochimie 2014; 105:22-35. [DOI: 10.1016/j.biochi.2014.07.009] [Citation(s) in RCA: 62] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2014] [Accepted: 07/01/2014] [Indexed: 11/23/2022]
|
72
|
Parrotta L, Ortuso F, Moraca F, Rocca R, Costa G, Alcaro S, Artese A. Targeting unimolecular G-quadruplex nucleic acids: a new paradigm for the drug discovery? Expert Opin Drug Discov 2014; 9:1167-87. [PMID: 25109710 DOI: 10.1517/17460441.2014.941353] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
INTRODUCTION G-quadruplexes (G4s) are targets of great interest because of their roles in crucial biological processes, such as aging and cancer. G4s are based on the formation of G-quartets, stabilised by Hoogsteen-type hydrogen bonds and by interaction with cations between the tetrads. These biologically relevant conformations were first discovered in eukaryotic chromosomal telomeric DNA, but have also been found in the proximal location of promoters in a number of human genes. Therefore, the extensive analysis of an intriguing target could move towards the rational drug design of new selective anticancer agents. AREAS COVERED The authors review G4 structural characterisation, with detailed insight related to the polymorphism issue. The authors describe the topologically distinct G4 structural forms and the factors involved in their interconversion mechanisms, such as the sequence of the oligonucleotides, the strand stoichiometry and orientation, the syn-anti conformation of the guanine glycosidic bonds and the G4 loop types and the environmental factors. Furthermore, the authors report several studies related to folding and unfolding kinetic profiles in order to understand the conformational view of monomolecular G4 formations. EXPERT OPINION G4 unimolecular nucleic acids can be considered as valid targets for the rational drug development of novel anticancer agents. Structural biology represents an essential link between the biology and medicinal chemistry knowledge in this field. In silico methods have already been demonstrated to be useful, especially if well integrated with biophysical tests. If this proves successful, the G4-targeting paradigm could also be extended to drug discovery beyond neoplastic pathologies.
Collapse
Affiliation(s)
- Lucia Parrotta
- Università degli Studi "Magna Græcia", Dipartimento di Scienze della Salute , Campus "S. Venuta", Viale Europa, Germaneto, 88100, Catanzaro , Italy
| | | | | | | | | | | | | |
Collapse
|
73
|
Zhu H, Xiao S, Wang L, Liang H. Communication: Asymmetrical cation movements through G-quadruplex DNA. J Chem Phys 2014; 141:041103. [DOI: 10.1063/1.4891218] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hong Zhu
- CAS Key Laboratory of Soft Matter Chemistry, Department of Polymer Science and Engineering, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| | - Shiyan Xiao
- CAS Key Laboratory of Soft Matter Chemistry, Department of Polymer Science and Engineering, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| | - Lei Wang
- CAS Key Laboratory of Soft Matter Chemistry, Department of Polymer Science and Engineering, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| | - Haojun Liang
- CAS Key Laboratory of Soft Matter Chemistry, Department of Polymer Science and Engineering, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
- Hefei National Laboratory for Physical Sciences at Microscale, Collaborative Innovation Center of Chemistry for Energy Materials, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| |
Collapse
|
74
|
Zgarbová M, Otyepka M, Šponer J, Lankaš F, Jurečka P. Base Pair Fraying in Molecular Dynamics Simulations of DNA and RNA. J Chem Theory Comput 2014; 10:3177-89. [DOI: 10.1021/ct500120v] [Citation(s) in RCA: 114] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Marie Zgarbová
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacky University, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Michal Otyepka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacky University, 17. listopadu 12, 77146 Olomouc, Czech Republic
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacky University, 17. listopadu 12, 77146 Olomouc, Czech Republic
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Filip Lankaš
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nam. 2, 166 10 Prague, Czech Republic
| | - Petr Jurečka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacky University, 17. listopadu 12, 77146 Olomouc, Czech Republic
| |
Collapse
|
75
|
Šponer J, Banáš P, Jurečka P, Zgarbová M, Kührová P, Havrila M, Krepl M, Stadlbauer P, Otyepka M. Molecular Dynamics Simulations of Nucleic Acids. From Tetranucleotides to the Ribosome. J Phys Chem Lett 2014; 5:1771-82. [PMID: 26270382 DOI: 10.1021/jz500557y] [Citation(s) in RCA: 113] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
We present a brief overview of explicit solvent molecular dynamics (MD) simulations of nucleic acids. We explain physical chemistry limitations of the simulations, namely, the molecular mechanics (MM) force field (FF) approximation and limited time scale. Further, we discuss relations and differences between simulations and experiments, compare standard and enhanced sampling simulations, discuss the role of starting structures, comment on different versions of nucleic acid FFs, and relate MM computations with contemporary quantum chemistry. Despite its limitations, we show that MD is a powerful technique for studying the structural dynamics of nucleic acids with a fast growing potential that substantially complements experimental results and aids their interpretation.
Collapse
Affiliation(s)
- Jiří Šponer
- †Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- ‡CEITEC - Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Pavel Banáš
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petr Jurečka
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Marie Zgarbová
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petra Kührová
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Marek Havrila
- †Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- ‡CEITEC - Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Miroslav Krepl
- †Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petr Stadlbauer
- †Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- §Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
76
|
Improta R. Quantum mechanical calculations unveil the structure and properties of the absorbing and emitting excited electronic states of guanine quadruplex. Chemistry 2014; 20:8106-15. [PMID: 24828154 DOI: 10.1002/chem.201400065] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2014] [Indexed: 12/20/2022]
Abstract
Herein, a full quantum mechanical study, in solution, of several models of guanine-quadruplex helices, both parallel and antiparallel, containing up to eight guanine residues, in their electronic excited state is reported. By exploiting TD-DFT calculations and including solvent effects by the polarizable continuum model, we provide the first atomistic description of the processes triggered by the absorption of UV light, reproducing and assigning the experimental optical and electronic circular dichroism spectra. The absorbing excited states are delocalized over multiple bases, whereas emission involves a stacked guanine dimer or a monomer. Several states, with a varying degree of localization and charge-transfer character, rule the photoexcited dynamics, which are deeply affected by the quadruplex topology. The lowest excited-state minimum for parallel quadruplex is an asymmetric excimer involving two stacked guanines, with a small charge transfer character, whereas for the anti-parallel structure, with the same topology of the thrombin binding aptamer, it is a fully symmetric excimer, characterized by a strong decrease of the stacking distance. A monomer-like decay path is the most relevant nonradiative decay pathway. Insights on the effect of the ions (K(+) or Na(+)) on the excited state decay are also provided.
Collapse
Affiliation(s)
- Roberto Improta
- Istituto Biostrutture e Bioimmagini, Consiglio Nazionale delle Ricerche, Via Mezzocannone 16, I-80134, Napoli (Italy).
| |
Collapse
|
77
|
Kruse H, Havrila M, Šponer J. QM Computations on Complete Nucleic Acids Building Blocks: Analysis of the Sarcin–Ricin RNA Motif Using DFT-D3, HF-3c, PM6-D3H, and MM Approaches. J Chem Theory Comput 2014; 10:2615-29. [DOI: 10.1021/ct500183w] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Holger Kruse
- CEITEC
− Central European Institute of Technology, Campus Bohunice, Kamenice
5, 625 00 Brno, Czech Republic
| | - Marek Havrila
- CEITEC
− Central European Institute of Technology, Campus Bohunice, Kamenice
5, 625 00 Brno, Czech Republic
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Jiřı́ Šponer
- CEITEC
− Central European Institute of Technology, Campus Bohunice, Kamenice
5, 625 00 Brno, Czech Republic
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| |
Collapse
|
78
|
Sun L, Jin H, Zhao X, Liu Z, Guan Y, Yang Z, Zhang L, Zhang L. Unfolding and conformational variations of thrombin-binding DNA aptamers: synthesis, circular dichroism and molecular dynamics simulations. ChemMedChem 2014; 9:993-1001. [PMID: 24715713 DOI: 10.1002/cmdc.201300564] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2014] [Revised: 03/14/2014] [Indexed: 01/08/2023]
Abstract
Thrombin-binding DNA aptamer (TBA), with a consensus 15-base sequence: d(GGTTGGTGTGGTTGG), can fold into an antiparallel unimolecular G-quadruplex structure that is necessary for its interaction with thrombin. For the first time, using steered molecular dynamics (SMD) simulations, we have successfully simulated the unfolding process of native TBA G-quadruplex. The unfolding pathway proposed is in agreement with previously reported experimental NMR data. Moreover, the critical intermediate structure in the unfolding pathway, predicted by the NMR results, was identified. The structural characteristics of several TBA oligonucleotides modified with locked nucleoside (LNA) or 2'-O-methyl-nucleoside (MNA) at different positions and number were also investigated by CD spectroscopy. An oligonucleotide substituted with either LNA or MNA at position 2 folds into a native-like G-quadruplex, while doubly substituted derivatives of TBA where LNA or MNA is incorporated at positions 11 and 14 are no longer able to form a G-quadruplex. Starting from the same initial intermediate structure, we successfully overcame sampling limitations, and simulated the large conformational variations in structures of native TBA and modified TBAs by classic MD. Analysis of the models showed that inversion of the glycosyl orientation at position 14 contributes significantly to the disruption of G-quadruplex formation in both of the di-substituted modified TBA systems. Our calculations provide a simple and reliable theoretical model that can be used to investigate and predict the effects of the modifications of an oligonucleotide on the resultant G-quadruplex structure. In addition, the computational protocol described can also be applied for other G-quadruplex systems.
Collapse
Affiliation(s)
- Lidan Sun
- State Key Laboratory of Natural and Biomimetic Drugs, School of Pharmaceutical Sciences, Peking University, Beijing 100191 (P. R. China)
| | | | | | | | | | | | | | | |
Collapse
|
79
|
Gkionis K, Kruse H, Platts JA, Mládek A, Koča J, Šponer J. Ion Binding to Quadruplex DNA Stems. Comparison of MM and QM Descriptions Reveals Sizable Polarization Effects Not Included in Contemporary Simulations. J Chem Theory Comput 2014; 10:1326-40. [DOI: 10.1021/ct4009969] [Citation(s) in RCA: 67] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Konstantinos Gkionis
- CEITEC
- Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Holger Kruse
- CEITEC
- Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - James A. Platts
- School
of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Arnošt Mládek
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Jaroslav Koča
- CEITEC
- Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jiří Šponer
- CEITEC
- Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| |
Collapse
|
80
|
Lai H, Xiao Y, Yan S, Tian F, Zhong C, Liu Y, Weng X, Zhou X. Symmetric cyanovinyl-pyridinium triphenylamine: a novel fluorescent switch-on probe for an antiparallel G-quadruplex. Analyst 2014; 139:1834-8. [DOI: 10.1039/c3an02269b] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
We describe a probe based on a cyanovinyl pyridinium triphenylamine (CPT) derivative, which showed fluorescent switch-on properties toward an antiparallel G-quadruplex.
Collapse
Affiliation(s)
- Hai Lai
- College of Chemistry and Molecular Sciences
- Wuhan University
- Wuhan 430072, P. R. China
| | - Yijie Xiao
- College of Chemistry and Molecular Sciences
- Wuhan University
- Wuhan 430072, P. R. China
| | - Shengyong Yan
- College of Chemistry and Molecular Sciences
- Wuhan University
- Wuhan 430072, P. R. China
| | - Fangfang Tian
- College of Chemistry and Molecular Sciences
- Wuhan University
- Wuhan 430072, P. R. China
| | - Cheng Zhong
- College of Chemistry and Molecular Sciences
- Wuhan University
- Wuhan 430072, P. R. China
| | - Yi Liu
- College of Chemistry and Molecular Sciences
- Wuhan University
- Wuhan 430072, P. R. China
| | - Xiaocheng Weng
- College of Chemistry and Molecular Sciences
- Wuhan University
- Wuhan 430072, P. R. China
- The State Key Laboratory of Natural and Biomimetic Drugs
- State
| | - Xiang Zhou
- College of Chemistry and Molecular Sciences
- Wuhan University
- Wuhan 430072, P. R. China
- The State Key Laboratory of Natural and Biomimetic Drugs
- State
| |
Collapse
|
81
|
Lech CJ, Phan AT, Michel-Beyerle ME, Voityuk AA. Electron-Hole Transfer in G-Quadruplexes with Different Tetrad Stacking Geometries: A Combined QM and MD Study. J Phys Chem B 2013; 117:9851-6. [DOI: 10.1021/jp404788t] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Christopher J. Lech
- Division of Physics and Applied Physics,
School of
Physical and Mathematical Sciences, Nanyang Technological University, Singapore
| | - Anh Tuân Phan
- Division of Physics and Applied Physics,
School of
Physical and Mathematical Sciences, Nanyang Technological University, Singapore
| | - Maria-Elisabeth Michel-Beyerle
- Division of Physics and Applied Physics,
School of
Physical and Mathematical Sciences, Nanyang Technological University, Singapore
| | - Alexander A. Voityuk
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Barcelona, Spain
and Institute of Computational Chemistry, Universitat de Girona, 17071 Girona, Spain
| |
Collapse
|
82
|
Zhu H, Xiao S, Liang H. Structural dynamics of human telomeric G-quadruplex loops studied by molecular dynamics simulations. PLoS One 2013; 8:e71380. [PMID: 23951152 PMCID: PMC3738534 DOI: 10.1371/journal.pone.0071380] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2012] [Accepted: 06/30/2013] [Indexed: 11/18/2022] Open
Abstract
Loops which are linkers connecting G-strands and supporting the G-tetrad core in G-quadruplex are important for biological roles of G-quadruplexes. TTA loop is a common sequence which mainly resides in human telomeric DNA (hTel) G-quadruplex. A series of molecular dynamics (MD) simulations were carried out to investigate the structural dynamics of TTA loops. We found that (1) the TA base pair formed in TTA loops are very stable, the occupied of all hydrogen bonds are more than 0.95. (2) The TA base pair makes the adjacent G-quartet more stable than others. (3) For the edgewise loop and the diagonal loop, most loop bases are stacking with others, only few bases have considerable freedom. (4) The stabilities of these stacking structures are distinct. Part of the loops, especially TA base pairs, and bases stacking with the G-quartet, maintain certain stable conformations in the simulation, but other parts, like TT and TA stacking structures, are not stable enough. For the first time, spontaneous conformational switches of TTA edgewise loops were observed in our long time MD simulations. (5) For double chain reversal loop, it is really hard to maintain a stable conformation in the long time simulation under present force fields (parm99 and parmbsc0), as it has multiple conformations with similar free energies.
Collapse
Affiliation(s)
- Hong Zhu
- CAS Key Laboratory of Soft Matter Chemistry, Department of Polymer Science and Engineering, University of Science and Technology of China, Hefei, Anhui, People's Republic of China
| | - Shiyan Xiao
- CAS Key Laboratory of Soft Matter Chemistry, Department of Polymer Science and Engineering, University of Science and Technology of China, Hefei, Anhui, People's Republic of China
| | - Haojun Liang
- CAS Key Laboratory of Soft Matter Chemistry, Department of Polymer Science and Engineering, University of Science and Technology of China, Hefei, Anhui, People's Republic of China
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui, People's Republic of China
| |
Collapse
|
83
|
Stadlbauer P, Krepl M, Cheatham TE, Koča J, Šponer J. Structural dynamics of possible late-stage intermediates in folding of quadruplex DNA studied by molecular simulations. Nucleic Acids Res 2013; 41:7128-43. [PMID: 23700306 PMCID: PMC3737530 DOI: 10.1093/nar/gkt412] [Citation(s) in RCA: 95] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2013] [Revised: 04/18/2013] [Accepted: 04/24/2013] [Indexed: 12/31/2022] Open
Abstract
Explicit solvent molecular dynamics simulations have been used to complement preceding experimental and computational studies of folding of guanine quadruplexes (G-DNA). We initiate early stages of unfolding of several G-DNAs by simulating them under no-salt conditions and then try to fold them back using standard excess salt simulations. There is a significant difference between G-DNAs with all-anti parallel stranded stems and those with stems containing mixtures of syn and anti guanosines. The most natural rearrangement for all-anti stems is a vertical mutual slippage of the strands. This leads to stems with reduced numbers of tetrads during unfolding and a reduction of strand slippage during refolding. The presence of syn nucleotides prevents mutual strand slippage; therefore, the antiparallel and hybrid quadruplexes initiate unfolding via separation of the individual strands. The simulations confirm the capability of G-DNA molecules to adopt numerous stable locally and globally misfolded structures. The key point for a proper individual folding attempt appears to be correct prior distribution of syn and anti nucleotides in all four G-strands. The results suggest that at the level of individual molecules, G-DNA folding is an extremely multi-pathway process that is slowed by numerous misfolding arrangements stabilized on highly variable timescales.
Collapse
Affiliation(s)
- Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT 84124, USA and CEITEC – Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT 84124, USA and CEITEC – Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Thomas E. Cheatham
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT 84124, USA and CEITEC – Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jaroslav Koča
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT 84124, USA and CEITEC – Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT 84124, USA and CEITEC – Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
84
|
Šponer J, Mládek A, Špačková N, Cang X, Cheatham TE, Grimme S. Relative stability of different DNA guanine quadruplex stem topologies derived using large-scale quantum-chemical computations. J Am Chem Soc 2013; 135:9785-96. [PMID: 23742743 PMCID: PMC3775466 DOI: 10.1021/ja402525c] [Citation(s) in RCA: 99] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
We provide theoretical predictions of the intrinsic stability of different arrangements of guanine quadruplex (G-DNA) stems. Most computational studies of nucleic acids have applied Molecular Mechanics (MM) approaches using simple pairwise-additive force fields. The principle limitation of such calculations is the highly approximate nature of the force fields. In this study, we for the first time apply accurate QM computations (DFT-D3 with large atomic orbital basis sets) to essentially complete DNA building blocks, seven different folds of the cation-stabilized two-quartet G-DNA stem, each having more than 250 atoms. The solvent effects are approximated by COSMO continuum solvent. We reveal sizable differences between MM and QM descriptions of relative energies of different G-DNA stems, which apparently reflect approximations of the DNA force field. Using the QM energy data, we propose correction to earlier free energy estimates of relative stabilities of different parallel, hybrid, and antiparallel G-stem folds based on classical simulations. The new energy ranking visibly improves the agreement between theory and experiment. We predict the 5'-anti-anti-3' GpG dinucleotide step to be the most stable one, closely followed by the 5'-syn-anti-3' step. The results are in good agreement with known experimental structures of 2-, 3-, and 4-quartet G-DNA stems. Besides providing specific results for G-DNA, our study highlights basic limitations of force field modeling of nucleic acids. Although QM computations have their own limitations, mainly the lack of conformational sampling and the approximate description of the solvent, they can substantially improve the quality of calculations currently relying exclusively on force fields.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Arnošt Mládek
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Nad’a Špačková
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- Department of Condensed Matter Physics, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic
| | - Xiaohui Cang
- Institute of Genetics, School of Life Science, Zhejiang University, Hangzhou, China 310058
| | - Thomas E. Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT 84124, USA
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute of Physical & Theoretical Chemistry, University of Bonn, Beringstrasse. 4, D-53115 Bonn, Germany
| |
Collapse
|
85
|
Karsisiotis AI, O'Kane C, Webba da Silva M. DNA quadruplex folding formalism--a tutorial on quadruplex topologies. Methods 2013; 64:28-35. [PMID: 23791747 DOI: 10.1016/j.ymeth.2013.06.004] [Citation(s) in RCA: 80] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2013] [Revised: 06/04/2013] [Accepted: 06/10/2013] [Indexed: 10/26/2022] Open
Abstract
Quadruplexes of DNA adopt a large variety of topologies that are dependent on their environment. We have been developing a formalism for quadruplex folding based on the relationship between base and its sugar--as defined by the glycosidic bond angle. By reducing the quadruplex stem to a description based on two finite states of the range of angles the glycosidic bond angle may adopt, the description of the relationships of type of loop and groove widths of a quadruplex stem are possible. In its current form this formalism has allowed for the prediction of some unimolecular quadruplex topologies. Its rules, whilst developed for unimolecular quadruplexes of three loops, are of general utility in understanding the interdependency of structural characteristics of multimolecular folds, as well as unimolecular quadruplexes of more than three loops. Here we describe current understanding of the interdependent structural features that define the quadruplex fold, and provide a tutorial for the use and application of this formalism.
Collapse
Affiliation(s)
- Andreas Ioannis Karsisiotis
- School of Pharmacy & Pharmaceutical Sciences, Biomedical Sciences Research Institute, University of Ulster, Cromore Road, BT52 1SA, Coleraine, UK
| | | | | |
Collapse
|
86
|
Šponer J, Šponer JE, Mládek A, Banáš P, Jurečka P, Otyepka M. How to understand quantum chemical computations on DNA and RNA systems? A practical guide for non-specialists. Methods 2013; 64:3-11. [PMID: 23747334 DOI: 10.1016/j.ymeth.2013.05.025] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2013] [Revised: 05/29/2013] [Accepted: 05/30/2013] [Indexed: 12/30/2022] Open
Abstract
In this review primarily written for non-experts we explain basic methodological aspects and interpretation of modern quantum chemical (QM) computations applied to nucleic acids. We introduce current reference QM computations on small model systems consisting of dozens of atoms. Then we comment on recent advance of fast and accurate dispersion-corrected density functional theory methods, which will allow computations of small but complete nucleic acids building blocks in the near future. The qualitative difference between QM and molecular mechanics (MM, force field) computations is discussed. We also explain relation of QM and molecular simulation computations to experiments.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolská 135, 612 65 Brno, Czech Republic; CEITEC - Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic.
| | | | | | | | | | | |
Collapse
|
87
|
Islam B, Sgobba M, Laughton C, Orozco M, Sponer J, Neidle S, Haider S. Conformational dynamics of the human propeller telomeric DNA quadruplex on a microsecond time scale. Nucleic Acids Res 2013; 41:2723-35. [PMID: 23293000 PMCID: PMC3575793 DOI: 10.1093/nar/gks1331] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2012] [Revised: 11/26/2012] [Accepted: 11/27/2012] [Indexed: 01/12/2023] Open
Abstract
The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. Although G-quartet stems have been well characterized, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 μs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.
Collapse
Affiliation(s)
- Barira Islam
- Centre for Cancer Research and Cell Biology, Queen’s University of Belfast, Belfast BT9 7BL, UK, School of Pharmacy, Nottingham University, University Park, Nottingham NG7 2RD, UK, Institute of Research in Biomedicine, Barcelona 08028, Spain, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovoplka 135, Brno 612 65, Czech Republic, Central European Institute of Technology, Campus Bohunice, Kamenice 5, Brno 625 00, Czech Republic and University College London, School of Pharmacy, Brunswick Square, London WC1N 1AX, UK
| | - Miriam Sgobba
- Centre for Cancer Research and Cell Biology, Queen’s University of Belfast, Belfast BT9 7BL, UK, School of Pharmacy, Nottingham University, University Park, Nottingham NG7 2RD, UK, Institute of Research in Biomedicine, Barcelona 08028, Spain, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovoplka 135, Brno 612 65, Czech Republic, Central European Institute of Technology, Campus Bohunice, Kamenice 5, Brno 625 00, Czech Republic and University College London, School of Pharmacy, Brunswick Square, London WC1N 1AX, UK
| | - Charlie Laughton
- Centre for Cancer Research and Cell Biology, Queen’s University of Belfast, Belfast BT9 7BL, UK, School of Pharmacy, Nottingham University, University Park, Nottingham NG7 2RD, UK, Institute of Research in Biomedicine, Barcelona 08028, Spain, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovoplka 135, Brno 612 65, Czech Republic, Central European Institute of Technology, Campus Bohunice, Kamenice 5, Brno 625 00, Czech Republic and University College London, School of Pharmacy, Brunswick Square, London WC1N 1AX, UK
| | - Modesto Orozco
- Centre for Cancer Research and Cell Biology, Queen’s University of Belfast, Belfast BT9 7BL, UK, School of Pharmacy, Nottingham University, University Park, Nottingham NG7 2RD, UK, Institute of Research in Biomedicine, Barcelona 08028, Spain, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovoplka 135, Brno 612 65, Czech Republic, Central European Institute of Technology, Campus Bohunice, Kamenice 5, Brno 625 00, Czech Republic and University College London, School of Pharmacy, Brunswick Square, London WC1N 1AX, UK
| | - Jiri Sponer
- Centre for Cancer Research and Cell Biology, Queen’s University of Belfast, Belfast BT9 7BL, UK, School of Pharmacy, Nottingham University, University Park, Nottingham NG7 2RD, UK, Institute of Research in Biomedicine, Barcelona 08028, Spain, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovoplka 135, Brno 612 65, Czech Republic, Central European Institute of Technology, Campus Bohunice, Kamenice 5, Brno 625 00, Czech Republic and University College London, School of Pharmacy, Brunswick Square, London WC1N 1AX, UK
| | - Stephen Neidle
- Centre for Cancer Research and Cell Biology, Queen’s University of Belfast, Belfast BT9 7BL, UK, School of Pharmacy, Nottingham University, University Park, Nottingham NG7 2RD, UK, Institute of Research in Biomedicine, Barcelona 08028, Spain, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovoplka 135, Brno 612 65, Czech Republic, Central European Institute of Technology, Campus Bohunice, Kamenice 5, Brno 625 00, Czech Republic and University College London, School of Pharmacy, Brunswick Square, London WC1N 1AX, UK
| | - Shozeb Haider
- Centre for Cancer Research and Cell Biology, Queen’s University of Belfast, Belfast BT9 7BL, UK, School of Pharmacy, Nottingham University, University Park, Nottingham NG7 2RD, UK, Institute of Research in Biomedicine, Barcelona 08028, Spain, Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovoplka 135, Brno 612 65, Czech Republic, Central European Institute of Technology, Campus Bohunice, Kamenice 5, Brno 625 00, Czech Republic and University College London, School of Pharmacy, Brunswick Square, London WC1N 1AX, UK
| |
Collapse
|
88
|
Smestad J, Maher LJ. Ion-dependent conformational switching by a DNA aptamer that induces remyelination in a mouse model of multiple sclerosis. Nucleic Acids Res 2013; 41:1329-42. [PMID: 23175609 PMCID: PMC3553947 DOI: 10.1093/nar/gks1093] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
We recently reported that a guanosine-rich 40-mer DNA aptamer (LJM-3064) mediates remyelination in the Theiler's murine encephalomyelitis virus mouse model of multiple sclerosis. Here, we characterize the G-quadruplex forms of this aptamer in vitro, and demonstrate using circular dichroism spectroscopy that LJM-3064 undergoes a monovalent ion-dependent conformational switch. In the presence of sodium ions and no potassium ions, LJM-3064 adopts an antiparallel-stranded G-quadruplex structure. When presented with low concentrations of potassium ions in a buffer that mimics the composition of interstitial fluid and blood plasma, LJM-3064 rapidly switches to a parallel-stranded G-quadruplex conformation, which is presumably the physiologically active folded form. We characterize these conformational states using dimethyl sulfate reactivity studies and Bal 31 nuclease probing. Our analysis indicates that only the 5'-terminal 26 nucleotides are involved in G-quadruplex formation. Thermodynamic characterization of LJM-3064 at physiologically relevant ion concentrations reveals the G-quadruplex to be metastable at human body temperature. These data provide important structural and thermodynamic insights that may be valuable in optimizing LJM-3064 as a therapeutic remyelinating agent.
Collapse
Affiliation(s)
- John Smestad
- Department of Biochemistry and Molecular Biology, Mayo Clinic College of Medicine, Rochester, MN 55905, USA
| | | |
Collapse
|
89
|
Abstract
G-quadruplexes constitute a class of nucleic acid structures defined by stacked guanine tetrads (or G-tetrads) with guanine bases from neighboring tetrads stacking with one another within the G-tetrad core. Individual G-quadruplexes can also stack with one another at their G-tetrad interface leading to higher-order structures as observed in telomeric repeat-containing DNA and RNA. In this study, we investigate how guanine base stacking influences the stability of G-quadruplexes and their stacked higher-order structures. A structural survey of the Protein Data Bank is conducted to characterize experimentally observed guanine base stacking geometries within the core of G-quadruplexes and at the interface between stacked G-quadruplex structures. We couple this survey with a systematic computational examination of stacked G-tetrad energy landscapes using quantum mechanical computations. Energy calculations of stacked G-tetrads reveal large energy differences of up to 12 kcal/mol between experimentally observed geometries at the interface of stacked G-quadruplexes. Energy landscapes are also computed using an AMBER molecular mechanics description of stacking energy and are shown to agree quite well with quantum mechanical calculated landscapes. Molecular dynamics simulations provide a structural explanation for the experimentally observed preference of parallel G-quadruplexes to stack in a 5′–5′ manner based on different accessible tetrad stacking modes at the stacking interfaces of 5′–5′ and 3′–3′ stacked G-quadruplexes.
Collapse
Affiliation(s)
- Christopher Jacques Lech
- Division of Physics and Applied Physics, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371, Singapore
| | | | | |
Collapse
|
90
|
New insights from molecular dynamic simulation studies of the multiple binding modes of a ligand with G-quadruplex DNA. J Comput Aided Mol Des 2012; 26:1355-68. [PMID: 23239169 DOI: 10.1007/s10822-012-9619-1] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2012] [Accepted: 11/27/2012] [Indexed: 01/06/2023]
Abstract
G-quadruplexes are higher-order DNA and RNA structures formed from guanine-rich sequences. These structures have recently emerged as a new class of potential molecular targets for anticancer drugs. An understanding of the three-dimensional interactions between small molecular ligands and their G-quadruplex targets in solution is crucial for rational drug design and the effective optimization of G-quadruplex ligands. Thus far, rational ligand design has been focused mainly on the G-quartet platform. It should be noted that small molecules can also bind to loop nucleotides, as observed in crystallography studies. Hence, it would be interesting to elucidate the mechanism underlying how ligands in distinct binding modes influence the flexibility of G-quadruplex. In the present study, based on a crystal structure analysis, the models of a tetra-substituted naphthalene diimide ligand bound to a telomeric G-quadruplex with different modes were built and simulated with a molecular dynamics simulation method. Based on a series of computational analyses, the structures, dynamics, and interactions of ligand-quadruplex complexes were studied. Our results suggest that the binding of the ligand to the loop is viable in aqueous solutions but dependent on the particular arrangement of the loop. The binding of the ligand to the loop enhances the flexibility of the G-quadruplex, while the binding of the ligand simultaneously to both the quartet and the loop diminishes its flexibility. These results add to our understanding of the effect of a ligand with different binding modes on G-quadruplex flexibility. Such an understanding will aid in the rational design of more selective and effective G-quadruplex binding ligands.
Collapse
|
91
|
Dršata T, Pérez A, Orozco M, Morozov AV, Sponer J, Lankaš F. Structure, Stiffness and Substates of the Dickerson-Drew Dodecamer. J Chem Theory Comput 2012; 9:707-721. [PMID: 23976886 DOI: 10.1021/ct300671y] [Citation(s) in RCA: 68] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The Dickerson-Drew dodecamer (DD) d-[CGCGAATTCGCG]2 is a prototypic B-DNA molecule whose sequence-specific structure and dynamics have been investigated by many experimental and computational studies. Here, we present an analysis of DD properties based on extensive atomistic molecular dynamics (MD) simulations using different ionic conditions and water models. The 0.6-2.4-µs-long MD trajectories are compared to modern crystallographic and NMR data. In the simulations, the duplex ends can adopt an alternative base-pairing, which influences the oligomer structure. A clear relationship between the BI/BII backbone substates and the basepair step conformation has been identified, extending previous findings and exposing an interesting structural polymorphism in the helix. For a given end pairing, distributions of the basepair step coordinates can be decomposed into Gaussian-like components associated with the BI/BII backbone states. The nonlocal stiffness matrices for a rigid-base mechanical model of DD are reported for the first time, suggesting salient stiffness features of the central A-tract. The Riemann distance and Kullback-Leibler divergence are used for stiffness matrix comparison. The basic structural parameters converge very well within 300 ns, convergence of the BI/BII populations and stiffness matrices is less sharp. Our work presents new findings about the DD structural dynamics, mechanical properties, and the coupling between basepair and backbone configurations, including their statistical reliability. The results may also be useful for optimizing future force fields for DNA.
Collapse
Affiliation(s)
- Tomáš Dršata
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 166 10 Prague, Czech Republic ; Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, Brno, Czech Republic
| | | | | | | | | | | |
Collapse
|
92
|
Kim E, Yang C, Pak Y. Free-Energy Landscape of a Thrombin-Binding DNA Aptamer in Aqueous Environment. J Chem Theory Comput 2012; 8:4845-51. [DOI: 10.1021/ct300714u] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Eunae Kim
- BK-21 Project
Team, College of Pharmacy, Chosun University, Gwangju 501-759, Republic of Korea
- Supercomputing Center, Korea Institute of Science and Technology Information, Daejeon 305-806, Republic of
Korea
| | - Changwon Yang
- Department of Chemistry
and Institute
of Functional Materials, Pusan National University, Busan 609-735, Republic of Korea
| | - Youngshang Pak
- Department of Chemistry
and Institute
of Functional Materials, Pusan National University, Busan 609-735, Republic of Korea
| |
Collapse
|
93
|
Šket P, Virgilio A, Esposito V, Galeone A, Plavec J. Strand directionality affects cation binding and movement within tetramolecular G-quadruplexes. Nucleic Acids Res 2012; 40:11047-57. [PMID: 22977177 PMCID: PMC3510487 DOI: 10.1093/nar/gks851] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Nuclear magnetic resonance study of G-quadruplex structures formed by d(TG3T) and its modified analogs containing a 5′-5′ or 3′-3′ inversion of polarity sites, namely d(3′TG5′-5′G2T3′), d(3′T5′-5′G3T3′) and d(5′TG3′-3′G2T5’) demonstrates formation of G-quadruplex structures with tetrameric topology and distinct cation-binding preferences. All oligonucleotides are able to form quadruplex structures with two binding sites, although the modified oligonucleotides also form, in variable amounts, quadruplex structures with only one bound cation. The inter-quartet cavities at the inversion of polarity sites bind ammonium ions less tightly than a naturally occurring 5′-3′ backbone. Exchange of 15 ions between G-quadruplex and bulk solution is faster at the 3′-end in comparison to the 5′-end. In addition to strand directionality, cation movement is influenced by formation of an all-syn G-quartet. Formation of such quartet has been observed also for the parent d(TG3T) that besides the canonical quadruplex with only all-anti G-quartets, forms a tetramolecular parallel quadruplex containing one all-syn G-quartet, never observed before in unmodified quadruplex structures.
Collapse
Affiliation(s)
- Primoz Šket
- Slovenian NMR Center, National Institute of Chemistry, Hajdrihova 19, Ljubljana SI-1000, Slovenia
| | | | | | | | | |
Collapse
|
94
|
Akhshi P, Acton G, Wu G. Molecular dynamics simulations to provide new insights into the asymmetrical ammonium ion movement inside of the [d(G3T4G4)]2 G-quadruplex DNA structure. J Phys Chem B 2012; 116:9363-70. [PMID: 22780684 DOI: 10.1021/jp304760k] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
We have used both adaptive biasing force (ABF) and regular molecular dynamics (MD) simulations to investigate the asymmetrical NH(4)(+) ion movement inside of a bimolecular G-quadruplex DNA structure [d(G(3)T(4)G(4))](2). The free-energy landscapes obtained from ABF MD simulations suggest that the NH(4)(+) ion exiting the [d(G(3)T(4)G(4))](2) G-quadruplex stem in the direction toward the edge-type loop (denoted as the upper direction) experiences a lower free-energy barrier than that toward the diagonal loop (denoted as the lower direction) by approximately 3-4 kcal mol(-1). This result is in qualitative agreement with the previous discovery made by Šket and Plavec on the same G-quadruplex structure from (15)N NMR experiments (J. Am. Chem. Soc. 2007, 129, 8794). In the Na(+) form of the same G-quadruplex, Na(+) ion movement was found to be symmetrical, with a free-energy barrier of only 5-7 kcal mol(-1) to cross all three G-quartets, that is, [d(G(3)T(4)G(4))](2) still exhibits ion-channel-like behaviors for Na(+) ions. On the basis of the new computational results, we hypothesize that the stiffness of a G-quartet is primarily determined by the base stacking interactions within the G-quadruplex stem. Therefore, the structural origin for the asymmetrical NH(4)(+) ion movement in [d(G(3)T(4)G(4))](2) is the presence of two different modes of base stacking around the NH(4)(+) binding sites, a more stable 5'-syn-anti mode between lower and central G-quartets and a less stable 5'-anti-anti mode between upper and central G-quartets. Simulations also suggest that loop topology at the end of a G-quadruplex stem only controls the direction at which an exiting NH(4)(+) ion reaches bulk solution but does not impose significant free-energy barriers.
Collapse
Affiliation(s)
- Parisa Akhshi
- Department of Chemistry, Queen's University, 90 Bader Lane, Kingston, Ontario, Canada K7L 3N6
| | | | | |
Collapse
|
95
|
šponer J, Cang X, Cheatham TE. Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures. Methods 2012; 57:25-39. [PMID: 22525788 PMCID: PMC3775459 DOI: 10.1016/j.ymeth.2012.04.005] [Citation(s) in RCA: 105] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2012] [Revised: 04/04/2012] [Accepted: 04/06/2012] [Indexed: 11/29/2022] Open
Abstract
The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.
Collapse
Affiliation(s)
- Jiří šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- CEITEC – Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Xiaohui Cang
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Thomas E. Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, Skaggs Hall 201, 2000 East 30 South, University of Utah, Salt Lake City, UT 84112, United States
| |
Collapse
|
96
|
Clark GR, Pytel PD, Squire CJ. The high-resolution crystal structure of a parallel intermolecular DNA G-4 quadruplex/drug complex employing syn glycosyl linkages. Nucleic Acids Res 2012; 40:5731-8. [PMID: 22373921 PMCID: PMC3384316 DOI: 10.1093/nar/gks193] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
We have determined the X-ray structure of the complex between the DNA quadruplex d(5′-GGGG-3′)4 and daunomycin, as a potential model for studying drug–telomere interactions. The structure was solved at 1.08 Å by direct methods in space group I4. The asymmetric unit comprises a linear arrangement of one d(GGGG) strand, four daunomycin molecules, a second d(GGGG) strand facing in the opposite direction to the first, and Na and Mg cations. The crystallographic 4-fold axis generates the biological unit, which is a 12-layered structure comprising two sets of four guanine layers, with four layers each of four daunomycins stacked between the 5′ faces of the two quadruplexes. The daunomycin layers fall into two groups which are novel in their mode of self assembly. The only contacts between daunomycin molecules within any one of these layers are van der Waals interactions, however there is substantial π–π stacking between successive daunomycin layers and also with adjacent guanine layers. The structure differs significantly from all other parallel d(TGGGGT)4 quadruplexes in that the 5′ guanine adopts the unusual syn glycosyl linkage, refuting the widespread belief that such conformations should all be anti. In contrast to the related d(TGGGGT)/daunomycin complex, there are no ligand–quadruplex groove insertion interactions.
Collapse
Affiliation(s)
- George R Clark
- Chemistry Department and School of Biological Sciences, The University of Auckland, Private Bag 92019, Auckland 1044, New Zealand.
| | | | | |
Collapse
|
97
|
Li R, Xiong C, Xiao Z, Ling L. Colorimetric detection of cholesterol with G-quadruplex-based DNAzymes and ABTS2-. Anal Chim Acta 2012; 724:80-5. [PMID: 22483213 DOI: 10.1016/j.aca.2012.02.015] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2011] [Revised: 02/02/2012] [Accepted: 02/09/2012] [Indexed: 12/17/2022]
Abstract
A novel colorimetric method for detection of cholesterol was developed with hemin-G-quadruplex DNAzyme by transducing oxidation of cholesterol into the color change of 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt (ABTS(2-)). Oligonucleotide 5'-GTGGGTAGGGCGGGTTGG-3' (Oligo-1) formed G-quadruplex structure in the presence of K(+), it acted as a horseradish peroxidase (HRP) mimicking DNAzyme when binding hemin and catalyzed the oxidation of colorless ABTS(2-) to green ABTS(·-) by H(2)O(2), which was produced by the reaction of cholesterol and oxygen that catalyzed by cholesterol oxidase. Therefore, the oxidation of cholesterol could be transduced into the color change of ABTS(2-) by combining these two reactions. Under the optimum conditions, the absorbance was proportional to the concentration of cholesterol over the range of 1.0-30 μM, with a linear regression equation of A=0.362+0.0256C (C: μM, R=0.998) and a detection limit of 0.10 μM (3σ/slope). Moreover, the practicability of the assay in the detection of cholesterol in human serum was studied as well.
Collapse
Affiliation(s)
- Ruimin Li
- School of Chemistry and Chemical Engineering, Sun Yat-Sen University, Guangzhou, PR China
| | | | | | | |
Collapse
|
98
|
Wei D, Parkinson GN, Reszka AP, Neidle S. Crystal structure of a c-kit promoter quadruplex reveals the structural role of metal ions and water molecules in maintaining loop conformation. Nucleic Acids Res 2012; 40:4691-700. [PMID: 22287624 PMCID: PMC3378867 DOI: 10.1093/nar/gks023] [Citation(s) in RCA: 105] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
We report here the 1.62 Å crystal structure of an intramolecular quadruplex DNA formed from a sequence in the promoter region of the c-kit gene. This is the first reported crystal structure of a promoter quadruplex and the first observation of localized magnesium ions in a quadruplex structure. The structure reveals that potassium and magnesium ions have an unexpected yet significant structural role in stabilizing particular quadruplex loops and grooves that is distinct from but in addition to the role of potassium ions in the ion channel at the centre of all quadruplex structures. The analysis also shows how ions cluster together with structured water molecules to stabilize the quadruplex arrangement. This particular quadruplex has been previously studied by NMR methods, and the present X-ray structure is in accord with the earlier topology assignment. However, as well as the observations of potassium and magnesium ions, the crystal structure has revealed a highly significant difference in the dimensions of the large cleft in the structure, which is a plausible target for small molecules. This difference can be understood by the stabilizing role of structured water networks.
Collapse
Affiliation(s)
- Dengguo Wei
- CRUK Biomolecular Structure Group, UCL School of Pharmacy, University College London, 29-39 Brunswick Square, WC1N 1AX, London, UK
| | | | | | | |
Collapse
|
99
|
Virgilio A, Esposito V, Citarella G, Pepe A, Mayol L, Galeone A. The insertion of two 8-methyl-2'-deoxyguanosine residues in tetramolecular quadruplex structures: trying to orientate the strands. Nucleic Acids Res 2011; 40:461-75. [PMID: 21908403 PMCID: PMC3245916 DOI: 10.1093/nar/gkr670] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this article, we report a structural study, based on NMR and CD spectroscopies, and molecular modelling of all possible d(TG3T) and d(TG4T) analogues containing two 8-methyl-2′-deoxyguanosine residues (M). Particularly, the potential ability of these modified residues to orientate the strands and then to affect the folding topology of tetramolecular quadruplex structures has been investigated. Oligodeoxynucleotides (ODNs) TMMGT (T12) and TMMGGT (F12) form parallel tetramolecular quadruplexes, characterized by an all-syn M-tetrad at the 5′-side stacked to all-anti M- and G-tetrads. ODNs TMGMT (T13) and TMGGMT (F14) form parallel tetramolecular quadruplexes, in which an all-anti G core is sandwiched between two all-syn M-tetrads at the 5′- and the 3′-side. Notably, the quadruplex formed by T13 corresponds to an unprecedented structure in which the syn residues exceed in number the anti ones. Conversely, ODN TGMGMT (F24) adopts a parallel arrangement in which all-anti G-tetrads alternate with all-syn M-tetrads. Most importantly, all data strongly suggest that ODN TMGMGT (F13) forms an unprecedented anti-parallel tetramolecular quadruplex in which G and M residues adopt anti and syn glycosidic conformations, respectively. This article opens up new understandings and perspectives about the intricate relationship between the quadruplex strands orientation and the glycosidic conformation of the residues.
Collapse
Affiliation(s)
- Antonella Virgilio
- Dipartimento di Chimica delle Sostanze Naturali, Università degli Studi di Napoli Federico II, Via D. Montesano 49, I-80131 Napoli, Italy
| | | | | | | | | | | |
Collapse
|
100
|
Reshetnikov RV, Sponer J, Rassokhina OI, Kopylov AM, Tsvetkov PO, Makarov AA, Golovin AV. Cation binding to 15-TBA quadruplex DNA is a multiple-pathway cation-dependent process. Nucleic Acids Res 2011; 39:9789-802. [PMID: 21893589 PMCID: PMC3239185 DOI: 10.1093/nar/gkr639] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
A combination of explicit solvent molecular dynamics simulation (30 simulations reaching 4 µs in total), hybrid quantum mechanics/molecular mechanics approach and isothermal titration calorimetry was used to investigate the atomistic picture of ion binding to 15-mer thrombin-binding quadruplex DNA (G-DNA) aptamer. Binding of ions to G-DNA is complex multiple pathway process, which is strongly affected by the type of the cation. The individual ion-binding events are substantially modulated by the connecting loops of the aptamer, which play several roles. They stabilize the molecule during time periods when the bound ions are not present, they modulate the route of the ion into the stem and they also stabilize the internal ions by closing the gates through which the ions enter the quadruplex. Using our extensive simulations, we for the first time observed full spontaneous exchange of internal cation between quadruplex molecule and bulk solvent at atomistic resolution. The simulation suggests that expulsion of the internally bound ion is correlated with initial binding of the incoming ion. The incoming ion then readily replaces the bound ion while minimizing any destabilization of the solute molecule during the exchange.
Collapse
Affiliation(s)
- Roman V Reshetnikov
- Department of Boiengineering and Bioinformatics, Lomonosov Moscow State University, GSP-1, Leninskie Gory, Moscow, 119991, Russian Federation
| | | | | | | | | | | | | |
Collapse
|