51
|
Molecular Dynamics Simulation of Homo-DNA: The Role of Crystal Packing in Duplex Conformation. CRYSTALS 2019. [DOI: 10.3390/cryst9100532] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
The (4′→6′)-linked DNA homolog 2′,3′-dideoxy-β-D-glucopyranosyl nucleic acid (dideoxy-glucose nucleic acid or homo-DNA) exhibits stable self-pairing of the Watson–Crick and reverse-Hoogsteen types, but does not cross-pair with DNA. Molecular modeling and NMR solution studies of homo-DNA duplexes pointed to a conformation that was nearly devoid of a twist and a stacking distance in excess of 4.5 Å. By contrast, the crystal structure of the homo-DNA octamer dd(CGAATTCG) revealed a right-handed duplex with average values for helical twist and rise of ca. 15° and 3.8 Å, respectively. Other key features of the structure were strongly inclined base-pair and backbone axes in the duplex with concomitant base-pair slide and cross-strand stacking, and the formation of a dimer across a crystallographic dyad with inter-duplex base swapping. To investigate the conformational flexibility of the homo-DNA duplex and a potential influence of lattice interactions on its geometry, we used molecular dynamics (MD) simulations of the crystallographically observed dimer of duplexes and an isolated duplex in the solution state. The dimer of duplexes showed limited conformational flexibility, and key parameters such as helical rise, twist, and base-pair slide exhibited only minor fluctuations. The single duplex was clearly more flexible by comparison and underwent partial unwinding, albeit without significant lengthening. Thus, base stacking was preserved in the isolated duplex and two adenosines extruded from the stack in the dimer of duplexes were reinserted into the duplex and pair with Ts in a Hoogsteen mode. Our results confirmed that efficient stacking in homo-DNA seen in the crystal structure of a dimer of duplexes was maintained in the separate duplex. Therefore, lattice interactions did not account for the different geometries of the homo-DNA duplex in the crystal and earlier models that resembled inclined ladders with large base-pair separations that precluded efficient stacking.
Collapse
|
52
|
Galindo-Murillo R, Barroso-Flores J. Hydrophobic unnatural base pairs show a Watson-Crick pairing in micro-second molecular dynamics simulations. J Biomol Struct Dyn 2019; 38:4098-4106. [PMID: 31542995 DOI: 10.1080/07391102.2019.1671898] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Two unnatural hydrophobic nucleotides, d5SICS (2,6-dimethyl-2H-isoquiniline-1-thione) and dNaM (2-methoxy-3-methylnaphthalene), were previously replicated in vivo by a modified E. coli strand, however, a consistent structure for their pairing in terms of specific and selective directional interactions remains elusive, as data from spectroscopy experiments and simulations are inconsistent. The proposed d5SICS-dNaM pairing has been suggested to be a stacked configuration as suggested by NMR data; simulations have failed to reproduce this configuration and a Watson-Crick like pairing is observed. Previously, we focused on reproducing the d5SICS-dNaM Unnatural Base Pair (UBP) paring using an older (bsc0) AMBER force field, which was not able to correctly reproduce the experimental data. We present our efforts to reproduce the experimental pairing using the current version of the AMBER DNA force fields (OL15 and bsc1), two water models (OPC and TIP3P) and external electrostatic stabilization by Mg2+ ions. Opposite to previously reported simulations, a Watson-Crick-like pairing with no hydrogen bonds persists throughout all our results. Despite our efforts to replicate the reported stacked conformation, we cannot confirm its plausibility nor obtain a consistent structure that is independent of the neighboring nucleotides. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT, USA
| | - Joaquín Barroso-Flores
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM, Estado de México, C.P., México.,Instituto de Química, Universidad Nacional Autónoma de México, Ciudad Universitaria, Ciudad de México, México
| |
Collapse
|
53
|
Hognon C, Gebus A, Barone G, Monari A. Human DNA Telomeres in Presence of Oxidative Lesions: The Crucial Role of Electrostatic Interactions on the Stability of Guanine Quadruplexes. Antioxidants (Basel) 2019; 8:antiox8090337. [PMID: 31443537 PMCID: PMC6770428 DOI: 10.3390/antiox8090337] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2019] [Revised: 08/14/2019] [Accepted: 08/20/2019] [Indexed: 01/23/2023] Open
Abstract
By using all atom molecular dynamics simulations, we studied the behavior of human DNA telomere sequences in guanine quadruplex (G4) conformation and in the presence of oxidative lesions, namely abasic sites. In particular, we evidenced that while removing one guanine base induces a significant alteration and destabilization of the involved leaflet, human telomere oligomers tend, in most cases, to maintain at least a partial quadruplex structure, eventually by replacing the empty site with undamaged guanines of different leaflets. This study shows that (i) the disruption of the quadruplex leaflets induces the release of at least one of the potassium cations embedded in the quadruplex channel and that (ii) the electrostatic interactions of the DNA sequence with the aforementioned cations are fundamental to the maintenance of the global quadruplex structure.
Collapse
Affiliation(s)
- Cecilia Hognon
- Université de Lorraine, CNRS, LPCT UMR 7019, F54000 Nancy, France.
| | - Adrien Gebus
- Université de Lorraine, CNRS, LPCT UMR 7019, F54000 Nancy, France
| | - Giampaolo Barone
- Department of Biological, Chenical and Pharmaceutical Sciences and Technologies, Università di Palermo, Viale delle Scienze, 90128 Palermo, Italy
| | - Antonio Monari
- Université de Lorraine, CNRS, LPCT UMR 7019, F54000 Nancy, France.
| |
Collapse
|
54
|
Hognon C, Besancenot V, Gruez A, Grandemange S, Monari A. Cooperative Effects of Cytosine Methylation on DNA Structure and Dynamics. J Phys Chem B 2019; 123:7365-7371. [PMID: 31365827 DOI: 10.1021/acs.jpcb.9b05835] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The behavior of the structural parameters of DNA considering different levels of methylation in CpG islands is studied by means of full-atom molecular dynamics simulations and electronic circular dichroism, both in an artificial model system and in a gene promoter sequence. It is demonstrated that methylation although intrinsically brings quite local perturbations may, if its level is high enough, induce cooperative effects that strongly modify the DNA backbone torsional parameters altering the helicity as compared to the nonmethylated case. Because methylation of the CpG island is correlated with the regulation of gene expression, understanding the structural modifications induced in DNA is crucial to characterize all the fine equilibria into play in epigenetics phenomena.
Collapse
Affiliation(s)
- Cécilia Hognon
- Université de Lorraine and CNRS, UMR 7019 LPCT , F-54000 Nancy , France
| | | | - Arnaud Gruez
- Université de Lorraine and CNRS, UMR 7356 IMOPA , F-54000 Nancy , France
| | | | - Antonio Monari
- Université de Lorraine and CNRS, UMR 7019 LPCT , F-54000 Nancy , France
| |
Collapse
|
55
|
Krupa P, Quoc Huy PD, Li MS. Properties of monomeric Aβ42 probed by different sampling methods and force fields: Role of energy components. J Chem Phys 2019. [DOI: 10.1063/1.5093184] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Pawel Krupa
- Institute of Physics Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Pham Dinh Quoc Huy
- Institute of Physics Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Mai Suan Li
- Institute of Physics Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
- Institute for Computational Science and Technology, SBI Building, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam
| |
Collapse
|
56
|
Satange R, Chang CK, Hou MH. A survey of recent unusual high-resolution DNA structures provoked by mismatches, repeats and ligand binding. Nucleic Acids Res 2019; 46:6416-6434. [PMID: 29945186 PMCID: PMC6061790 DOI: 10.1093/nar/gky561] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2018] [Accepted: 06/08/2018] [Indexed: 12/20/2022] Open
Abstract
The structure of the DNA duplex is arguably one of the most important biological structures elucidated in modern history. DNA duplex structure is closely associated with essential biological functions such as DNA replication and RNA transcription. In addition to the classical A-, B- and Z-DNA conformations, DNA duplexes are capable of assuming a variety of alternative conformations depending on the sequence and environmental context. A considerable number of these unusual DNA duplex structures have been identified in the past decade, and some of them have been found to be closely associated with different biological functions and pathological conditions. In this manuscript, we review a selection of unusual DNA duplex structures, particularly those originating from base pair mismatch, repetitive sequence motifs and ligand-induced structures. Although the biological significance of these novel structures has not yet been established in most cases, the illustrated conformational versatility of DNA could have relevance for pharmaceutical or nanotechnology development. A perspective on the future directions of this field is also presented.
Collapse
Affiliation(s)
- Roshan Satange
- Institute of Genomics and Bioinformatics, National Chung Hsing University, Taichung, Taiwan.,Ph.D. Program in Medical Biotechnology, National Chung Hsing University, Taichung, Taiwan
| | - Chung-Ke Chang
- Institute of Biomedical Sciences, Academia Sinica, Taipei, Taiwan
| | - Ming-Hon Hou
- Institute of Genomics and Bioinformatics, National Chung Hsing University, Taichung, Taiwan.,Ph.D. Program in Medical Biotechnology, National Chung Hsing University, Taichung, Taiwan
| |
Collapse
|
57
|
Ekesan Ş, York DM. Framework for Conducting and Analyzing Crystal Simulations of Nucleic Acids to Aid in Modern Force Field Evaluation. J Phys Chem B 2019; 123:4611-4624. [PMID: 31002511 PMCID: PMC6614744 DOI: 10.1021/acs.jpcb.8b11923] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Crystal simulations provide useful tools, along with solution simulations, to test nucleic acid force fields, but should be interpreted with care owing to the difficulty of establishing the environmental conditions needed to reproduce experimental crystal packing. These challenges underscore the need to construct proper protocols for carrying out crystal simulations and analyzing results to identify the origin of deviations from crystallographic data. Toward this end, we introduce a novel framework for B-factor decomposition into additive intramolecular, rotational, and translational atomic fluctuation components and partitioning of each of these components into individual asymmetric unit and lattice contributions. We apply the framework to a benchmark set of A-DNA, Z-DNA, and B-DNA double helix systems of various chain lengths. Overall, the intramolecular deviations from the crystal were quite small (≤1.0 Å), suggesting high accuracy of the force field, whereas crystal packing was not well reproduced. The present work establishes a framework to conduct and analyze crystal simulations that ultimately take on issues of crystal packing and can provide insight into nucleic acid force fields.
Collapse
Affiliation(s)
- Şölen Ekesan
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854 , United States
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854 , United States
| |
Collapse
|
58
|
Wang X, Sun Z. Determination of Base-Flipping Free-Energy Landscapes from Nonequilibrium Stratification. J Chem Inf Model 2019; 59:2980-2994. [PMID: 31124677 DOI: 10.1021/acs.jcim.9b00263] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Correct calculation of the variation of free energy upon base flipping is crucial in understanding the dynamics of DNA systems. The free-energy landscape along the flipping pathway gives the thermodynamic stability and the flexibility of base-paired states. Although numerous free-energy simulations are performed in the base flipping cases, no theoretically rigorous nonequilibrium techniques are devised and employed to investigate the thermodynamics of base flipping. In the current work, we report a general nonequilibrium stratification scheme for the efficient calculation of the free-energy landscape of base flipping in DNA duplex. We carefully monitor the convergence behavior of the equilibrium sampling based free-energy simulation and the nonequilibrium stratification and determine the empirical length of time blocks required for converged sampling. Comparison between the performances of the equilibrium umbrella sampling and the nonequilibrium stratification is given. The results show that nonequilibrium free-energy simulation achieves similar accuracy and efficiency compared with the equilibrium enhanced sampling technique in the base flipping cases. We further test a convergence criterion we previously proposed and it comes out that the convergence determined by this criterion agrees with those given by the time-invariant behavior of PMF and the nonlinear dependence of standard deviation on the sample size.
Collapse
Affiliation(s)
- Xiaohui Wang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering , East China Normal University , Shanghai 200062 , China.,Institute of Computational Science , Università della Svizzera Italiana (USI) , Via Giuseppe Buffi 13 , CH-6900 , Lugano , Ticino , Switzerland
| | - Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering , East China Normal University , Shanghai 200062 , China.,Computational Biomedicine (IAS-5/INM-9) , Forschungszentrum Jülich , Jülich 52425 , Germany
| |
Collapse
|
59
|
Balasubramanian R, Pal S, Joshi H, Rao A, Naik A, Varma M, Chakraborty B, Maiti PK. DNA Translocation through Hybrid Bilayer Nanopores. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2019; 123:11908-11916. [PMID: 31316684 PMCID: PMC6636640 DOI: 10.1021/acs.jpcc.9b00399] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Pore functionalization has been explored by several groups as a strategy to control DNA translocation through solid-state nanopores. Here we present a hybrid nanopore system consisting of single-layer graphene and a DNA origami layer to achieve base-selective control of DNA translocation rate through aligned nanopores of the two layers. This is achieved by incorporating unpaired dangling bases called overhangs to the origami near the pore region. Molecular dynamics simulations were used to optimize the design of the origami nanopore and the overhangs. Specifically, we considered the influence of the number and spatial distribution of overhangs on translocation times. The simulations revealed that specific interactions between the overhangs and the translocating single-stranded DNA resulted in base-specific residence times.
Collapse
Affiliation(s)
| | - Sohini Pal
- Centre for Nano Science and Engineering, Indian Institute of Science, Bangalore 560012, India
| | - Himanshu Joshi
- Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Anjana Rao
- Division of Signaling and Gene Expression, La Jolla Institute for Allergy and Immunology, La Jolla, California 92037, United States
| | - Akshay Naik
- Centre for Nano Science and Engineering, Indian Institute of Science, Bangalore 560012, India
| | - Manoj Varma
- Centre for Nano Science and Engineering, Indian Institute of Science, Bangalore 560012, India
- Corresponding Authors: Dr. Banani Chakraborty: ., Prof. Manoj Varma: ., Prof. Prabal K. Maiti:
| | - Banani Chakraborty
- Department of Chemical Engineering, Indian Institute of Science, Bangalore 560012, India
- Corresponding Authors: Dr. Banani Chakraborty: ., Prof. Manoj Varma: ., Prof. Prabal K. Maiti:
| | - Prabal K. Maiti
- Department of Physics, Indian Institute of Science, Bangalore 560012, India
- Corresponding Authors: Dr. Banani Chakraborty: ., Prof. Manoj Varma: ., Prof. Prabal K. Maiti:
| |
Collapse
|
60
|
Igaev M, Kutzner C, Bock LV, Vaiana AC, Grubmüller H. Automated cryo-EM structure refinement using correlation-driven molecular dynamics. eLife 2019; 8:e43542. [PMID: 30829573 PMCID: PMC6424565 DOI: 10.7554/elife.43542] [Citation(s) in RCA: 59] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2018] [Accepted: 03/03/2019] [Indexed: 12/22/2022] Open
Abstract
We present a correlation-driven molecular dynamics (CDMD) method for automated refinement of atomistic models into cryo-electron microscopy (cryo-EM) maps at resolutions ranging from near-atomic to subnanometer. It utilizes a chemically accurate force field and thermodynamic sampling to improve the real-space correlation between the modeled structure and the cryo-EM map. Our framework employs a gradual increase in resolution and map-model agreement as well as simulated annealing, and allows fully automated refinement without manual intervention or any additional rotamer- and backbone-specific restraints. Using multiple challenging systems covering a wide range of map resolutions, system sizes, starting model geometries and distances from the target state, we assess the quality of generated models in terms of both model accuracy and potential of overfitting. To provide an objective comparison, we apply several well-established methods across all examples and demonstrate that CDMD performs best in most cases.
Collapse
Affiliation(s)
- Maxim Igaev
- Department of Theoretical and Computational BiophysicsMax Planck Institute for Biophysical ChemistryGöttingenGermany
| | - Carsten Kutzner
- Department of Theoretical and Computational BiophysicsMax Planck Institute for Biophysical ChemistryGöttingenGermany
| | - Lars V Bock
- Department of Theoretical and Computational BiophysicsMax Planck Institute for Biophysical ChemistryGöttingenGermany
| | - Andrea C Vaiana
- Department of Theoretical and Computational BiophysicsMax Planck Institute for Biophysical ChemistryGöttingenGermany
| | - Helmut Grubmüller
- Department of Theoretical and Computational BiophysicsMax Planck Institute for Biophysical ChemistryGöttingenGermany
| |
Collapse
|
61
|
Nilsson L, Villa A. Modeling and Simulation of Oligonucleotide Hybrids: Outlining a Strategy. Methods Mol Biol 2019; 2036:113-126. [PMID: 31410793 DOI: 10.1007/978-1-4939-9670-4_6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Molecular dynamics simulations with a state-of-the-art force field provide an atomistic detailed description of the structural and thermodynamic features of biomolecules. Effects of chemical modifications and of the environment such as sequence, solvent, and ionic strength can explicitly be taken into account. Molecular simulation techniques can also provide insight in change in binding affinity, in protonation (pKa shift) and tautomeric propensity due to changes in the environment or in the molecular system. The quality and reliability of a simulation depend on the quality of the force field and on the reproducibility of the data, and validation depends on the availability of suitable experimental data. Here, we describe the workflow to investigate oligonucleotide hybrids using molecular simulation including hardware and software information.
Collapse
Affiliation(s)
- Lennart Nilsson
- Department of Biosciences and Nutrition, Karolinska Institutet, Huddinge, Sweden
| | - Alessandra Villa
- Department of Biosciences and Nutrition, Karolinska Institutet, Huddinge, Sweden.
| |
Collapse
|
62
|
Lemkul JA, MacKerell AD. Polarizable force field for RNA based on the classical drude oscillator. J Comput Chem 2018; 39:2624-2646. [PMID: 30515902 PMCID: PMC6284239 DOI: 10.1002/jcc.25709] [Citation(s) in RCA: 59] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2018] [Revised: 08/01/2018] [Accepted: 09/23/2018] [Indexed: 12/15/2022]
Abstract
RNA molecules are highly dynamic and capable of adopting a wide range of complex, folded structures. The factors driving the folding and dynamics of these structures are dependent on a balance of base pairing, hydration, base stacking, ion interactions, and the conformational sampling of the 2'-hydroxyl group in the ribose sugar. The representation of these features is a challenge for empirical force fields used in molecular dynamics simulations. Toward meeting this challenge, the inclusion of explicit electronic polarization is important in accurately modeling RNA structure. In this work, we present a polarizable force field for RNA based on the classical Drude oscillator model, which represents electronic degrees of freedom via negatively charged particles attached to their parent atoms by harmonic springs. Beginning with parametrization against quantum mechanical base stacking interaction energy and conformational energy data, we have extended the Drude-2017 nucleic acid force field to include RNA. The conformational sampling of a range of RNA sequences were used to validate the force field, including canonical A-form RNA duplexes, stem-loops, and complex tertiary folds that bind multiple Mg2+ ions. Overall, the Drude-2017 RNA force field reproduces important properties of these structures, including the conformational sampling of the 2'-hydroxyl and key interactions with Mg2+ ions. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| |
Collapse
|
63
|
Galindo-Murillo R, Cheatham TE, Hopkins RC. Exploring potentially alternative non-canonical DNA duplex structures through simulation. J Biomol Struct Dyn 2018; 37:2201-2210. [PMID: 30047316 DOI: 10.1080/07391102.2018.1483839] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Hopkins proposed an alternative and chirally distinct family of double-stranded DNA (dsDNA) models that have antiparallel chains with 5'→3' senses opposite to those of the right-handed Watson-Crick (WC) family. Termed configuration II, this family of dsDNA models contains both right-handed (II-R) and left-handed (II-L) forms, with Z-DNA as an example of the latter. Relative interstrand binding energies for six DNA duplex models, two each of configuration I-R (standard WC canonical B-DNA), II-R, and II-L for the duplex d(CGCGAATTCGCG), have been estimated under identical conditions using MM-PBSA analysis from molecular dynamics trajectories using three different AMBER force fields. These simulations support the stereo chemical soundness of configuration II dsDNA forms. Recent force fields (Barcelona Supercomputing Center [BSC] [bsc1] and Olomouc 2015 [OL15]) successfully render stable II-L structures, whereas the previous force field, bsc0, generated stable II-R structures, although with an energy difference between II-R and II-L of ∼30 kcal/mol. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- a Department of Medicinal Chemistry , L. S. Skaggs Pharmacy Institute, University of Utah , Salt Lake City , UT , USA
| | - Thomas E Cheatham
- a Department of Medicinal Chemistry , L. S. Skaggs Pharmacy Institute, University of Utah , Salt Lake City , UT , USA
| | - Robert C Hopkins
- a Department of Medicinal Chemistry , L. S. Skaggs Pharmacy Institute, University of Utah , Salt Lake City , UT , USA.,b Department of Chemistry, University of Houston-Clear Lake , Houston , TX , USACommunicated by Ramaswamy H. Sarma
| |
Collapse
|
64
|
Karwowski BT. The AT Interstrand Cross-Link: Structure, Electronic Properties, and Influence on Charge Transfer in dsDNA. MOLECULAR THERAPY-NUCLEIC ACIDS 2018; 13:665-685. [PMID: 30500729 PMCID: PMC6258832 DOI: 10.1016/j.omtn.2018.10.014] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/01/2018] [Revised: 10/22/2018] [Accepted: 10/22/2018] [Indexed: 12/01/2022]
Abstract
The interaction of chemical and physical agents with genetic material can lead to almost 80 different DNA damage formations. The targeted intentional DNA damage by radiotherapy or chemotherapy is a front-line anticancer therapy. An interstrand cross-link can result from ionization radiation or specific chemical agents, such as trans-/cisplatin activity. Here, the influence of the adenine and thymidine (AT) interstrand linkage, the covalent bond between the adenine N6 and thymidine C5 methylene group, on the isolated base pair as well as double-stranded DNA (dsDNA) was taken into quantum mechanical/molecular mechanical (QM/MM) consideration at the m062x/6-31+G*:UFF level of theory in the aqueous phase. All the results presented in this article, for the first time, show that an AT-interstrand cross-link (ICL) changes the positive and negative charge migration process due to a higher activation energy forced by the cross-link’s presence. However, the final radical cation destination in cross-linked DNA is left in the same place as in a native double-stranded-deoxyoligonucleotide. Additionally, the direction of the radical anion transfer was found to be opposite to that of native dsDNA. Therefore, it can be postulated that the appearance of the AT-ICL does not disturb the hole migration in the double helix, with subsequent effective changes in the electron migration process.
Collapse
Affiliation(s)
- Boleslaw T Karwowski
- DNA Damage Laboratory of Food Science Department, Faculty of Pharmacy, Medical University of Lodz, ul. Muszynskiego 1, 90-151 Lodz, Poland.
| |
Collapse
|
65
|
Affiliation(s)
- Vito Genna
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| | - Elisa Donati
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| | - Marco De Vivo
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| |
Collapse
|
66
|
Grotz KK, Nueesch MF, Holmstrom ED, Heinz M, Stelzl LS, Schuler B, Hummer G. Dispersion Correction Alleviates Dye Stacking of Single-Stranded DNA and RNA in Simulations of Single-Molecule Fluorescence Experiments. J Phys Chem B 2018; 122:11626-11639. [PMID: 30285443 DOI: 10.1021/acs.jpcb.8b07537] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
We combine single-molecule Förster resonance energy transfer (single-molecule FRET) experiments with extensive all-atom molecular dynamics (MD) simulations (>100 μs) to characterize the conformational ensembles of single-stranded (ss) DNA and RNA in solution. From MD simulations with explicit dyes attached to single-stranded nucleic acids via flexible linkers, we calculate FRET efficiencies and fluorescence anisotropy decays. We find that dispersion-corrected water models alleviate the problem of overly abundant interactions between fluorescent dyes and the aromatic ring systems of nucleobases. To model dye motions in a computationally efficient and conformationally exhaustive manner, we introduce a dye-conformer library, built from simulations of dinucleotides with covalently attached dye molecules. We use this library to calculate FRET efficiencies for dT19, dA19, and rA19 simulated without explicit labels over a wide range of salt concentrations. For end-labeled homopolymeric pyrimidine ssDNA, MD simulations with the parmBSC1 force field capture the overall trend in salt-dependence of single-molecule FRET based distance measurements. For homopolymeric purine ssRNA and ssDNA, the DESRES and parmBSC1 force fields, respectively, provide useful starting points, even though our comparison also identifies clear deviations from experiment.
Collapse
Affiliation(s)
- Kara K Grotz
- Department of Theoretical Biophysics , Max Planck Institute of Biophysics , 60438 Frankfurt am Main , Germany
| | - Mark F Nueesch
- Department of Biochemistry , University of Zurich , 8057 Zurich , Switzerland
| | - Erik D Holmstrom
- Department of Biochemistry , University of Zurich , 8057 Zurich , Switzerland
| | - Marcel Heinz
- Department of Theoretical Biophysics , Max Planck Institute of Biophysics , 60438 Frankfurt am Main , Germany
| | - Lukas S Stelzl
- Department of Theoretical Biophysics , Max Planck Institute of Biophysics , 60438 Frankfurt am Main , Germany
| | - Benjamin Schuler
- Department of Biochemistry , University of Zurich , 8057 Zurich , Switzerland.,Department of Physics , University of Zurich , 8057 Zurich , Switzerland
| | - Gerhard Hummer
- Department of Theoretical Biophysics , Max Planck Institute of Biophysics , 60438 Frankfurt am Main , Germany.,Institute of Biophysics , Goethe University Frankfurt , 60438 Frankfurt am Main , Germany
| |
Collapse
|
67
|
Azuma R, Kishi S, Gutmann G, Konagaya A. <b>All-atom molecular dynamics of film supported flat-shaped DNA origami in water</b>. CHEM-BIO INFORMATICS JOURNAL 2018. [DOI: 10.1273/cbij.18.96] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
|
68
|
Mao P, Brown AJ, Esaki S, Lockwood S, Poon GMK, Smerdon MJ, Roberts SA, Wyrick JJ. ETS transcription factors induce a unique UV damage signature that drives recurrent mutagenesis in melanoma. Nat Commun 2018; 9:2626. [PMID: 29980679 PMCID: PMC6035183 DOI: 10.1038/s41467-018-05064-0] [Citation(s) in RCA: 84] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2018] [Accepted: 06/07/2018] [Indexed: 11/12/2022] Open
Abstract
Recurrent mutations are frequently associated with transcription factor (TF) binding sites (TFBS) in melanoma, but the mechanism driving mutagenesis at TFBS is unclear. Here, we use a method called CPD-seq to map the distribution of UV-induced cyclobutane pyrimidine dimers (CPDs) across the human genome at single nucleotide resolution. Our results indicate that CPD lesions are elevated at active TFBS, an effect that is primarily due to E26 transformation-specific (ETS) TFs. We show that ETS TFs induce a unique signature of CPD hotspots that are highly correlated with recurrent mutations in melanomas, despite high repair activity at these sites. ETS1 protein renders its DNA binding targets extremely susceptible to UV damage in vitro, due to binding-induced perturbations in the DNA structure that favor CPD formation. These findings define a mechanism responsible for recurrent mutations in melanoma and reveal that DNA binding by ETS TFs is inherently mutagenic in UV-exposed cells. Many factors contribute to mutation hotspots in cancer cells. Here the authors map UV damage at single-nucleotide resolution across the human genome and find that binding sites of ETS transcription factors are especially prone to forming UV lesions, leading to mutation hotspots in melanoma.
Collapse
Affiliation(s)
- Peng Mao
- School of Molecular Biosciences, Washington State University, Pullman, WA, 99164, USA
| | - Alexander J Brown
- School of Molecular Biosciences, Washington State University, Pullman, WA, 99164, USA
| | - Shingo Esaki
- Department of Chemistry, Georgia State University, Atlanta, GA, 30303, USA
| | - Svetlana Lockwood
- Paul G. Allen School for Global Animal Health, Washington State University, Pullman, WA, 99164, USA
| | - Gregory M K Poon
- Department of Chemistry, Georgia State University, Atlanta, GA, 30303, USA.,Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, 30303, USA
| | - Michael J Smerdon
- School of Molecular Biosciences, Washington State University, Pullman, WA, 99164, USA
| | - Steven A Roberts
- School of Molecular Biosciences, Washington State University, Pullman, WA, 99164, USA. .,Center for Reproductive Biology, Washington State University, Pullman, WA, 99164, USA.
| | - John J Wyrick
- School of Molecular Biosciences, Washington State University, Pullman, WA, 99164, USA. .,Center for Reproductive Biology, Washington State University, Pullman, WA, 99164, USA.
| |
Collapse
|
69
|
Marazzi M, Gattuso H, Monari A, Assfeld X. Steady-State Linear and Non-linear Optical Spectroscopy of Organic Chromophores and Bio-macromolecules. Front Chem 2018; 6:86. [PMID: 29666792 PMCID: PMC5891624 DOI: 10.3389/fchem.2018.00086] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2018] [Accepted: 03/12/2018] [Indexed: 01/05/2023] Open
Abstract
Bio-macromolecules as DNA, lipid membranes and (poly)peptides are essential compounds at the core of biological systems. The development of techniques and methodologies for their characterization is therefore necessary and of utmost interest, even though difficulties can be experienced due to their intrinsic complex nature. Among these methods, spectroscopies, relying on optical properties are especially important to determine their macromolecular structures and behaviors, as well as the possible interactions and reactivity with external dyes—often drugs or pollutants—that can (photo)sensitize the bio-macromolecule leading to eventual chemical modifications, thus damages. In this review, we will focus on the theoretical simulation of electronic spectroscopies of bio-macromolecules, considering their secondary structure and including their interaction with different kind of (photo)sensitizers. Namely, absorption, emission and electronic circular dichroism (CD) spectra are calculated and compared with the available experimental data. Non-linear properties will be also taken into account by two-photon absorption, a highly promising technique (i) to enhance absorption in the red and infra-red windows and (ii) to enhance spatial resolution. Methodologically, the implications of using implicit and explicit solvent, coupled to quantum and thermal samplings of the phase space, will be addressed. Especially, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are explored for a comparison with solely QM methods, in order to address the necessity to consider an accurate description of environmental effects on spectroscopic properties of biological systems.
Collapse
Affiliation(s)
- Marco Marazzi
- Laboratoire de Physique et Chimie Théoriques, Université de Lorraine-Nancy, UMR 7019, Vandoeuvre-lés-Nancy, France.,Laboratoire de Physique et Chimie Théoriques, Centre National de la Recherche Scientifique, UMR 7019, Vandoeuvre-lès-Nancy, France.,Departamento de Química, Centro de Investigacíon en Síntesis Química (CISQ), Universidad de La Rioja, Logroño, Spain
| | - Hugo Gattuso
- Laboratoire de Physique et Chimie Théoriques, Université de Lorraine-Nancy, UMR 7019, Vandoeuvre-lés-Nancy, France.,Laboratoire de Physique et Chimie Théoriques, Centre National de la Recherche Scientifique, UMR 7019, Vandoeuvre-lès-Nancy, France
| | - Antonio Monari
- Laboratoire de Physique et Chimie Théoriques, Université de Lorraine-Nancy, UMR 7019, Vandoeuvre-lés-Nancy, France.,Laboratoire de Physique et Chimie Théoriques, Centre National de la Recherche Scientifique, UMR 7019, Vandoeuvre-lès-Nancy, France
| | - Xavier Assfeld
- Laboratoire de Physique et Chimie Théoriques, Université de Lorraine-Nancy, UMR 7019, Vandoeuvre-lés-Nancy, France.,Laboratoire de Physique et Chimie Théoriques, Centre National de la Recherche Scientifique, UMR 7019, Vandoeuvre-lès-Nancy, France
| |
Collapse
|
70
|
Reshetnikov RV, Stolyarova AV, Zalevsky AO, Panteleev DY, Pavlova GV, Klinov DV, Golovin AV, Protopopova AD. A coarse-grained model for DNA origami. Nucleic Acids Res 2018; 46:1102-1112. [PMID: 29267876 PMCID: PMC5814798 DOI: 10.1093/nar/gkx1262] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2017] [Revised: 11/14/2017] [Accepted: 12/07/2017] [Indexed: 01/20/2023] Open
Abstract
Modeling tools provide a valuable support for DNA origami design. However, current solutions have limited application for conformational analysis of the designs. In this work we present a tool for a thorough study of DNA origami structure and dynamics. The tool is based on a novel coarse-grained model dedicated to geometry optimization and conformational analysis of DNA origami. We explored the ability of the model to predict dynamic behavior, global shapes, and fine details of two single-layer systems designed in hexagonal and square lattices using atomic force microscopy, Förster resonance energy transfer spectroscopy, and all-atom molecular dynamic simulations for validation of the results. We also examined the performance of the model for multilayer systems by simulation of DNA origami with published cryo-electron microscopy and atomic force microscopy structures. A good agreement between the simulated and experimental data makes the model suitable for conformational analysis of DNA origami objects. The tool is available at http://vsb.fbb.msu.ru/cosm as a web-service and as a standalone version.
Collapse
Affiliation(s)
- Roman V Reshetnikov
- Institute of Gene Biology, Russian Academy of Sciences, Vavilova str., 34/5, 119334 Moscow, Russia
- A.N.Belozersky Institute of Physical and Chemical Biology, Lomonosov Moscow State University, Leninskye gory, 1-40, 119992 Moscow, Russia
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, GSP-1, Leninskiye Gory, 1-73, 119234 Moscow, Russia
| | - Anastasia V Stolyarova
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, GSP-1, Leninskiye Gory, 1-73, 119234 Moscow, Russia
- Skolkovo Institute of Science and Technology, Nobel Street 3, 143026 Moscow, Russia
| | - Arthur O Zalevsky
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, GSP-1, Leninskiye Gory, 1-73, 119234 Moscow, Russia
| | - Dmitry Y Panteleev
- Institute of Gene Biology, Russian Academy of Sciences, Vavilova str., 34/5, 119334 Moscow, Russia
| | - Galina V Pavlova
- Institute of Gene Biology, Russian Academy of Sciences, Vavilova str., 34/5, 119334 Moscow, Russia
| | - Dmitry V Klinov
- Federal Research and Clinical Center of Physical-Chemical Medicine, Malaya Pirogovskaya str. 1a, 119435 Moscow, Russia
- Moscow Institute of Physics and Technology (State University), 9 Institutskiy per. Dolgoprudny, 141700 Moscow Region, Russia
| | - Andrey V Golovin
- A.N.Belozersky Institute of Physical and Chemical Biology, Lomonosov Moscow State University, Leninskye gory, 1-40, 119992 Moscow, Russia
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, GSP-1, Leninskiye Gory, 1-73, 119234 Moscow, Russia
- Sechenov First Moscow State Medical University, Institute of Molecular Medicine, Trubetskaya str. 8-2, 119991 Moscow, Russia
| | - Anna D Protopopova
- Federal Research and Clinical Center of Physical-Chemical Medicine, Malaya Pirogovskaya str. 1a, 119435 Moscow, Russia
- Department of Cell & Developmental Biology, Perelman School of Medicine, University of Pennsylvania, BRB II/III 421 Curie Boulevard, Philadelphia, PA 19104, USA
| |
Collapse
|
71
|
Chakraborty D, Wales DJ. Energy Landscape and Pathways for Transitions between Watson-Crick and Hoogsteen Base Pairing in DNA. J Phys Chem Lett 2018; 9:229-241. [PMID: 29240425 DOI: 10.1021/acs.jpclett.7b01933] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The recent discovery that Hoogsteen (HG) base pairs are widespread in DNA across diverse sequences and positional contexts could have important implications for understanding DNA replication and DNA-protein recognition. While evidence is emerging that the Hoogsteen conformation could be a thermodynamically accessible conformation of the DNA duplex and provide a means to expand its functionality, relatively little is known about the molecular mechanism underlying the Watson-Crick (WC) to HG transition. In this Perspective, we describe pathways and kinetics for this transition at an atomic level of detail, using the energy landscape perspective. We show that competition between the duplex conformations results in a double funnel landscape, which explains some recent experimental observations. The interconversion pathways feature a number of intermediates, with a variable number of WC and HG base pairs. The relatively slow kinetics, with possible deviations from two-state behavior, suggest that this conformational switch is likely to be a challenging target for both simulation and experiment.
Collapse
Affiliation(s)
- Debayan Chakraborty
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Department of Chemistry, The University of Texas at Austin , 24th Street Stop A5300, Austin, Texas 78712, United States
| | - David J Wales
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
72
|
Yang C, Kulkarni M, Lim M, Pak Y. Insilico direct folding of thrombin-binding aptamer G-quadruplex at all-atom level. Nucleic Acids Res 2017; 45:12648-12656. [PMID: 29112755 PMCID: PMC5728390 DOI: 10.1093/nar/gkx1079] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2017] [Revised: 10/13/2017] [Accepted: 10/20/2017] [Indexed: 12/25/2022] Open
Abstract
The reversible folding of the thrombin-binding DNA aptamer G-quadruplexes (GQs) (TBA-15) starting from fully unfolded states was demonstrated using a prolonged time scale (10-12 μs) parallel tempering metadynamics (PTMetaD) simulation method in conjunction with a modified version of the AMBER bsc1 force field. For unbiased descriptions of the folding free energy landscape of TBA-15, this force field was minimally modified. From this direct folding simulation using the modified bsc1 force field, reasonably converged free energy landscapes were obtained in K+-rich aqueous solution (150 mM), providing detailed atomistic pictures of GQ folding mechanisms for TBA-15. This study found that the TBA folding occurred via multiple folding pathways with two major free energy barriers of 13 and 15 kcal/mol in the presence of several intermediate states of G-triplex variants. The early formation of these intermediates was associated with a single K+ ion capturing. Interestingly, these intermediate states appear to undergo facile transitions among themselves through relatively small energy barriers.
Collapse
Affiliation(s)
- Changwon Yang
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 609-735, South Korea
| | - Mandar Kulkarni
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 609-735, South Korea
| | - Manho Lim
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 609-735, South Korea
| | - Youngshang Pak
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 609-735, South Korea
| |
Collapse
|
73
|
Effect of single-residue bulges on RNA double-helical structures: crystallographic database analysis and molecular dynamics simulation studies. J Mol Model 2017; 23:311. [DOI: 10.1007/s00894-017-3480-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2017] [Accepted: 09/19/2017] [Indexed: 11/26/2022]
|
74
|
Stelzl LS, Erlenbach N, Heinz M, Prisner TF, Hummer G. Resolving the Conformational Dynamics of DNA with Ångstrom Resolution by Pulsed Electron-Electron Double Resonance and Molecular Dynamics. J Am Chem Soc 2017; 139:11674-11677. [PMID: 28777549 DOI: 10.1021/jacs.7b05363] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Pulsed electron-electron double resonance (PELDOR/DEER) experiments of nucleic acids with rigid spin labels provide highly accurate distance and orientation information. Here we combine PELDOR experiments with molecular dynamics (MD) simulations to arrive at an atomistic view of the conformational dynamics of DNA. The MD simulations closely reproduce the PELDOR time traces, and demonstrate that bending, in addition to twist-stretch motions, underpin the sub-μs dynamics of DNA. PELDOR experiments correctly rank DNA force fields and resolve subtle differences in the conformational ensembles of nucleic acids, on the order of 1-2 Å. Long-range distance and angle measurements with rigid spin labels provide critical input for the refinement of computer models and the elucidation of the structure and dynamics of complex biomolecules.
Collapse
Affiliation(s)
- Lukas S Stelzl
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics , 60438 Frankfurt am Main, Germany
| | - Nicole Erlenbach
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt , 60438 Frankfurt am Main, Germany
| | - Marcel Heinz
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics , 60438 Frankfurt am Main, Germany
| | - Thomas F Prisner
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt , 60438 Frankfurt am Main, Germany
| | - Gerhard Hummer
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics , 60438 Frankfurt am Main, Germany.,Institute for Biophysics, Goethe University Frankfurt , 60438 Frankfurt am Main, Germany
| |
Collapse
|
75
|
Andrabi M, Hutchins AP, Miranda-Saavedra D, Kono H, Nussinov R, Mizuguchi K, Ahmad S. Predicting conformational ensembles and genome-wide transcription factor binding sites from DNA sequences. Sci Rep 2017; 7:4071. [PMID: 28642456 PMCID: PMC5481346 DOI: 10.1038/s41598-017-03199-6] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2016] [Accepted: 04/26/2017] [Indexed: 12/24/2022] Open
Abstract
DNA shape is emerging as an important determinant of transcription factor binding beyond just the DNA sequence. The only tool for large scale DNA shape estimates, DNAshape was derived from Monte-Carlo simulations and predicts four broad and static DNA shape features, Propeller twist, Helical twist, Minor groove width and Roll. The contributions of other shape features e.g. Shift, Slide and Opening cannot be evaluated using DNAshape. Here, we report a novel method DynaSeq, which predicts molecular dynamics-derived ensembles of a more exhaustive set of DNA shape features. We compared the DNAshape and DynaSeq predictions for the common features and applied both to predict the genome-wide binding sites of 1312 TFs available from protein interaction quantification (PIQ) data. The results indicate a good agreement between the two methods for the common shape features and point to advantages in using DynaSeq. Predictive models employing ensembles from individual conformational parameters revealed that base-pair opening - known to be important in strand separation - was the best predictor of transcription factor-binding sites (TFBS) followed by features employed by DNAshape. Of note, TFBS could be predicted not only from the features at the target motif sites, but also from those as far as 200 nucleotides away from the motif.
Collapse
Affiliation(s)
- Munazah Andrabi
- National Institutes of Biomedical Innovation Health and Nutrition, 7-6-8, Saito-Asagi, Ibaraki, Osaka, 5670085, Japan
- Faculty of Biology,Medicine and Health, Michael Smith Building, The University of Manchester, Dover Street, Manchester, M13 9PT, UK
| | - Andrew Paul Hutchins
- Department of Biology, Southern University of Science and Technology of China, Shenzhen, 518055, China
| | - Diego Miranda-Saavedra
- World Premier International (WPI) Immunology Frontier Research Center (IFReC), Osaka University, 3-1 Yamadaoka, Suita, 565-0871, Osaka, Japan
- Centro de Biología Molecular Severo Ochoa, CSIC/Universidad Autónoma de Madrid, 28049, Madrid, Spain
- Department of Computer Science, University of Oxford Wolfson Building, Parks Road, OXFORD, OX1 3QD, United Kingdom
| | - Hidetoshi Kono
- Molecular Modeling and Simulation (MMS) Group, National Institutes for Quantum and Radiological Science and Technology, 8-1-7, Umemidai, Kizugawa, Kyoto, 619-0215, Japan
| | - Ruth Nussinov
- National Cancer Institute, Cancer and Inflammation Program, Leidos Biomedical Research, Inc. Frederick, Maryland, USA
- Department of Biochemistry and Human Genetics, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Kenji Mizuguchi
- National Institutes of Biomedical Innovation Health and Nutrition, 7-6-8, Saito-Asagi, Ibaraki, Osaka, 5670085, Japan
| | - Shandar Ahmad
- National Institutes of Biomedical Innovation Health and Nutrition, 7-6-8, Saito-Asagi, Ibaraki, Osaka, 5670085, Japan.
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Mehrauli Road, New Delhi, 110067, India.
| |
Collapse
|
76
|
Lemkul JA, MacKerell AD. Polarizable Force Field for DNA Based on the Classical Drude Oscillator: II. Microsecond Molecular Dynamics Simulations of Duplex DNA. J Chem Theory Comput 2017; 13:2072-2085. [PMID: 28398748 PMCID: PMC5485260 DOI: 10.1021/acs.jctc.7b00068] [Citation(s) in RCA: 61] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The structure and dynamics of DNA are governed by a sensitive balance between base stacking and pairing, hydration, and interactions with ions. Force-field models that include explicit representations of electronic polarization are capable of more accurately modeling the subtle details of these interactions versus commonly used additive force fields. In this work, we validate our recently refined polarizable force field for DNA based on the classical Drude oscillator model, in which electronic degrees of freedom are represented as negatively charged particles attached to their parent atoms via harmonic springs. The previous version of the force field, called Drude-2013, produced stable A- and B-DNA trajectories on the order of hundreds of nanoseconds, but deficiencies were identified that included weak base stacking ultimately leading to distortion of B-DNA duplexes and unstable Z-DNA. As a result of extensive refinement of base nonbonded terms and bonded parameters in the deoxyribofuranose sugar and phosphodiester backbone, we demonstrate that the new version of the Drude DNA force field is capable of simulating A- and B-forms of DNA on the microsecond time scale and the resulting conformational ensembles agree well with a broad set of experimental properties, including solution X-ray scattering profiles. In addition, simulations of Z-form duplex DNA in its crystal environment are stable on the order of 100 ns. The revised force field is to be called Drude-2017.
Collapse
Affiliation(s)
- Justin A. Lemkul
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| |
Collapse
|
77
|
Lemkul JA, MacKerell AD. Polarizable Force Field for DNA Based on the Classical Drude Oscillator: I. Refinement Using Quantum Mechanical Base Stacking and Conformational Energetics. J Chem Theory Comput 2017; 13:2053-2071. [PMID: 28399366 DOI: 10.1021/acs.jctc.7b00067] [Citation(s) in RCA: 59] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Empirical force fields seek to relate the configuration of a set of atoms to its energy, thus yielding the forces governing its dynamics, using classical physics rather than more expensive quantum mechanical calculations that are computationally intractable for large systems. Most force fields used to simulate biomolecular systems use fixed atomic partial charges, neglecting the influence of electronic polarization, instead making use of a mean-field approximation that may not be transferable across environments. Recent hardware and software developments make polarizable simulations feasible, and to this end, polarizable force fields represent the next generation of molecular dynamics simulation technology. In this work, we describe the refinement of a polarizable force field for DNA based on the classical Drude oscillator model by targeting quantum mechanical interaction energies and conformational energy profiles of model compounds necessary to build a complete DNA force field. The parametrization strategy employed in the present work seeks to correct weak base stacking in A- and B-DNA and the unwinding of Z-DNA observed in the previous version of the force field, called Drude-2013. Refinement of base nonbonded terms and reparametrization of dihedral terms in the glycosidic linkage, deoxyribofuranose rings, and important backbone torsions resulted in improved agreement with quantum mechanical potential energy surfaces. Notably, we expand on previous efforts by explicitly including Z-DNA conformational energetics in the refinement.
Collapse
Affiliation(s)
- Justin A Lemkul
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| |
Collapse
|
78
|
Balaceanu A, Pasi M, Dans PD, Hospital A, Lavery R, Orozco M. The Role of Unconventional Hydrogen Bonds in Determining BII Propensities in B-DNA. J Phys Chem Lett 2017; 8:21-28. [PMID: 27935717 DOI: 10.1021/acs.jpclett.6b02451] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
An accurate understanding of DNA backbone transitions is likely to be the key for elucidating the puzzle of the intricate sequence-dependent mechanical properties that govern most of the biologically relevant functions of the double helix. One factor believed to be important in indirect recognition within protein-DNA complexes is the combined effect of two DNA backbone torsions (ε and ζ) which give rise to the well-known BI/BII conformational equilibrium. In this work we explain the sequence-dependent BII propensity observed in RpY steps (R = purine; Y = pyrimidine) at the tetranucleotide level with the help of a previously undetected C-H···O contact between atoms belonging to adjacent bases. Our results are supported by extensive multimicrosecond molecular dynamics simulations from the Ascona B-DNA Consortium, high-level quantum mechanical calculations, and data mining of the experimental structures deposited in the Protein Data Bank.
Collapse
Affiliation(s)
- Alexandra Balaceanu
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology , Baldiri Reixac 10-12, Barcelona 08028, Spain
- Joint BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine , Baldiri Reixac 10-12, Barcelona 08028, Spain
| | - Marco Pasi
- MMSB, Univ. Lyon I/CNRS UMR 5086, Institut de Biologie et Chimie des Protéines, 7 Passage du Vercors, Lyon 69367, France
- School of Pharmacy and Centre for Biomolecular Sciences, University of Nottingham , University Park NG7 2RD, U.K
| | - Pablo D Dans
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology , Baldiri Reixac 10-12, Barcelona 08028, Spain
- Joint BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine , Baldiri Reixac 10-12, Barcelona 08028, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology , Baldiri Reixac 10-12, Barcelona 08028, Spain
- Joint BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine , Baldiri Reixac 10-12, Barcelona 08028, Spain
| | - Richard Lavery
- MMSB, Univ. Lyon I/CNRS UMR 5086, Institut de Biologie et Chimie des Protéines, 7 Passage du Vercors, Lyon 69367, France
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology , Baldiri Reixac 10-12, Barcelona 08028, Spain
- Joint BSC-IRB Program in Computational Biology, Institute for Research in Biomedicine , Baldiri Reixac 10-12, Barcelona 08028, Spain
- Department of Biochemistry and Biomedicine, Faculty of Biology, University of Barcelona . Diagonal 643, Barcelona 08028, Spain
| |
Collapse
|