1
|
Franco LR, Toledo KCF, Matias TA, Araujo CM, Araki K, Coutinho K. Theoretical investigation of solvent and oxidation/deprotonation effects on the electronic structure of a mononuclear Ru-aqua-polypyridine complex in aqueous solution. Phys Chem Chem Phys 2023; 25:24475-24494. [PMID: 37655780 DOI: 10.1039/d3cp02154h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Mononuclear polypyridine ruthenium (Ru) complexes can catalyze various reactions, including water splitting, and can also serve as photosensitizers in solar cells. Despite recent progress in their synthesis, accurately modeling their physicochemical properties, particularly in solution, remains challenging. Herein, we conduct a theoretical investigation of the structural and electronic properties of a mononuclear Ru-aqua polypyridine complex in aqueous solution, considering five of its possible oxidation/protonation states species: [RuII(H2O)(py)(bpy)2]2+, [RuII(OH)(py)(bpy)2]+, [RuIII(H2O)(py)(bpy)2]3+, [RuIII(OH)(py)(bpy)2]2+ and [RuIV(O)(py)(bpy)2]2+, where py = pyridine and bpy = 2,2'-bipyridine. At first, we investigate the impact of proton-coupled and non-coupled electron transfer reactions on the geometry and electronic structure of the complexes in vacuum and in solution, using an implicit solvent model. Then, using a sequential multiscale approach that combines quantum mechanics and molecular mechanics (S-QM/MM), we examine the explicit solvent effects on the electronic excitations of the complexes, and compare them with the experimental results. The complexes were synthesized, and their absorption spectra measured in aqueous solution. To accurately describe the QM interactions between the metal center and the aqueous ligand in the MM simulations, we developed new force field parameters for the Ru atom. We analyze the solvent structure around the complexes and account for its explicit influence on the polarization and electronic excitations of the complexes. Notably, accounting for the explicit solvent polarization effects of the first solvation shells is essential to correctly describe the energy of the electronic transitions, and the explicit treatment of the hydrogen bonds at the QM level in the excitation calculations improves the accuracy of the description of the metal-to-ligand charge-transfer bands. Transition density matrix analysis is used to characterize all electronic transitions in the visible and ultraviolet ranges according to their charge-transfer (CT) character. This study elucidates the electronic structure of those ruthenium polypyridyl complexes in aqueous solution and underscores the importance of precisely describing solvent effects, which can be achieved employing the S-QM/MM method.
Collapse
Affiliation(s)
- Leandro Rezende Franco
- Instituto de Fisica, Universidade de Sao Paulo, Cidade Universitaria, 05508-090 Sao Paulo, SP, Brazil
- Department of Engineering and Physics, Karlstad University, 65188 Karlstad, Sweden.
| | | | - Tiago Araujo Matias
- Instituto de Quimica, Universidade de Sao Paulo, Av. Lineu Prestes 748, Butanta, 05508-000 Sao Paulo, SP, Brazil
| | - C Moyses Araujo
- Department of Engineering and Physics, Karlstad University, 65188 Karlstad, Sweden.
- Materials Theory Division, Department of Physics and Astronomy, Ångström Laboratory, Uppsala University, 75120 Uppsala, Sweden
| | - Koiti Araki
- Instituto de Quimica, Universidade de Sao Paulo, Av. Lineu Prestes 748, Butanta, 05508-000 Sao Paulo, SP, Brazil
| | - Kaline Coutinho
- Instituto de Fisica, Universidade de Sao Paulo, Cidade Universitaria, 05508-090 Sao Paulo, SP, Brazil
| |
Collapse
|
2
|
Pederson K, Meints GA, Drobny GP. Base Dynamics in the HhaI Protein Binding Site. J Phys Chem B 2023; 127:7266-7275. [PMID: 37561575 PMCID: PMC10461302 DOI: 10.1021/acs.jpcb.3c03687] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Revised: 07/18/2023] [Indexed: 08/12/2023]
Abstract
Protein-DNA interactions play an important role in numerous biological functions within the living cell. In many of these interactions, the DNA helix is significantly distorted upon protein-DNA complex formation. The HhaI restriction-modification system is one such system, where the methylation target is flipped out of the helix when bound to the methyltransferase. However, the base flipping mechanism is not well understood. The dynamics of the binding site of the HhaI methyltransferase and endonuclease (underlined) within the DNA oligomer [d(G1A2T3A4G5C6G7C8T9A10T11C12)]2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs deuterated on the base of nucleotides within and flanking the [5'-GCGC-3']2 sequence indicate that all of these positions are structurally flexible. Previously, conformational flexibility within the phosphodiester backbone and furanose ring within the target sequence has been observed and hypothesized to play a role in the distortion mechanism. However, whether that distortion was occurring through an active or passive mechanism remained unclear. These NMR data demonstrate that although the [5'-GCGC-3']2 sequence is dynamic, the target cytosine is not passively flipping out of the double-helix on the millisecond-picosecond time scale. Additionally, although previous studies have shown that both the furanose ring and phosphodiester backbone experience a change in dynamics upon methylation, which may play a role in recognition and cleavage by the endonuclease, our observations here indicate that methylation has no effect on the dynamics of the base itself.
Collapse
Affiliation(s)
- Kari Pederson
- Department
of Chemistry & Biochemistry, California
State University at Dominguez Hills, Carson, California 90747, United States
| | - Gary A. Meints
- Department
of Chemistry, Missouri State University, Springfield, Missouri 65897, United States
| | - Gary P. Drobny
- Department
of Chemistry, University of Washington, Seattle, Washington 98195-1700, United
States
| |
Collapse
|
3
|
Parisis N, Dans PD, Jbara M, Singh B, Schausi-Tiffoche D, Molina-Serrano D, Brun-Heath I, Hendrychová D, Maity SK, Buitrago D, Lema R, Nait Achour T, Giunta S, Girardot M, Talarek N, Rofidal V, Danezi K, Coudreuse D, Prioleau MN, Feil R, Orozco M, Brik A, Wu PYJ, Krasinska L, Fisher D. Histone H3 serine-57 is a CHK1 substrate whose phosphorylation affects DNA repair. Nat Commun 2023; 14:5104. [PMID: 37607906 PMCID: PMC10444856 DOI: 10.1038/s41467-023-40843-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2018] [Accepted: 08/12/2023] [Indexed: 08/24/2023] Open
Abstract
Histone post-translational modifications promote a chromatin environment that controls transcription, DNA replication and repair, but surprisingly few phosphorylations have been documented. We report the discovery of histone H3 serine-57 phosphorylation (H3S57ph) and show that it is implicated in different DNA repair pathways from fungi to vertebrates. We identified CHK1 as a major human H3S57 kinase, and disrupting or constitutively mimicking H3S57ph had opposing effects on rate of recovery from replication stress, 53BP1 chromatin binding, and dependency on RAD52. In fission yeast, mutation of all H3 alleles to S57A abrogated DNA repair by both non-homologous end-joining and homologous recombination, while cells with phospho-mimicking S57D alleles were partly compromised for both repair pathways, presented aberrant Rad52 foci and were strongly sensitised to replication stress. Mechanistically, H3S57ph loosens DNA-histone contacts, increasing nucleosome mobility, and interacts with H3K56. Our results suggest that dynamic phosphorylation of H3S57 is required for DNA repair and recovery from replication stress, opening avenues for investigating the role of this modification in other DNA-related processes.
Collapse
Affiliation(s)
- Nikolaos Parisis
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
- Equipe labellisée Ligue contre le Cancer, Paris, France
- BPMP, CNRS, INRA, Montpellier SupAgro, University of Montpellier, Montpellier, France
- Institut Jacques Monod, CNRS, University Paris Diderot, Paris, France
| | - Pablo D Dans
- IRB Barcelona, BIST, Barcelona, Spain
- Bioinformatics Unit, Institute Pasteur of Montevideo, Montevideo, Uruguay
- Department of Biological Sciences, CENUR North Riverside, University of the Republic (UdelaR), Salto, Uruguay
| | - Muhammad Jbara
- Schulich Faculty of Chemistry, Technion Israel Institute of Technology, Haifa, Israel
- School of Chemistry, Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv, Israel
| | | | | | | | | | - Denisa Hendrychová
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
- Equipe labellisée Ligue contre le Cancer, Paris, France
- Department of Experimental Biology, Faculty of Science, Palacký University Olomouc, Olomouc, Czech Republic
| | - Suman Kumar Maity
- Schulich Faculty of Chemistry, Technion Israel Institute of Technology, Haifa, Israel
| | | | | | - Thiziri Nait Achour
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
- Equipe labellisée Ligue contre le Cancer, Paris, France
| | - Simona Giunta
- The Rockefeller University, New York, NY, USA
- Laboratory of Genome Evolution, Department of Biology and Biotechnology "Charles Darwin", University of Rome Sapienza, Rome, Italy
| | - Michael Girardot
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
| | - Nicolas Talarek
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
| | - Valérie Rofidal
- BPMP, CNRS, INRA, Montpellier SupAgro, University of Montpellier, Montpellier, France
| | - Katerina Danezi
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
- Equipe labellisée Ligue contre le Cancer, Paris, France
| | - Damien Coudreuse
- IGDR, CNRS, University of Rennes, Rennes, France
- IBGC, CNRS, University of Bordeaux, Bordeaux, France
| | | | - Robert Feil
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France
| | | | - Ashraf Brik
- Schulich Faculty of Chemistry, Technion Israel Institute of Technology, Haifa, Israel
| | - Pei-Yun Jenny Wu
- IGDR, CNRS, University of Rennes, Rennes, France
- IBGC, CNRS, University of Bordeaux, Bordeaux, France
| | - Liliana Krasinska
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France.
- Equipe labellisée Ligue contre le Cancer, Paris, France.
| | - Daniel Fisher
- IGMM, CNRS, INSERM, University of Montpellier, Montpellier, France.
- Equipe labellisée Ligue contre le Cancer, Paris, France.
| |
Collapse
|
4
|
Rodriguez J, Battistini F, Learte-Aymamí S, Orozco M, Mascareñas JL. Molecular dynamics modelling of the interaction of a synthetic zinc-finger miniprotein with DNA. RSC Chem Biol 2023; 4:486-493. [PMID: 37415868 PMCID: PMC10320839 DOI: 10.1039/d3cb00053b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2023] [Accepted: 05/09/2023] [Indexed: 07/08/2023] Open
Abstract
We report the modelling of the DNA complex of an artificial miniprotein composed of two zinc finger modules and an AT-hook linking peptide. The computational study provides for the first time a structural view of these types of complexes, dissecting interactions that are key to modulate their stability. The relevance of these interactions was validated experimentally. These results confirm the potential of this type of computational approach for studying peptide-DNA complexes and suggest that they could be very useful for the rational design of non-natural, DNA binding miniproteins.
Collapse
Affiliation(s)
- Jessica Rodriguez
- Centro Singular de Investigación en Química Biolóxica e Materiais Moleculares (CIQUS), and Departamento de Química Orgánica, Universidade de Santiago de Compostela Rúa Jenaro de la Fuente s/n 15782 Santiago de Compostela Spain
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST) Baldiri Reixac 10-12 08028 Barcelona Spain
| | - Federica Battistini
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST) Baldiri Reixac 10-12 08028 Barcelona Spain
- Department of Biochemistry and Molecular Biology, University of Barcelona 08028 Barcelona Spain
| | - Soraya Learte-Aymamí
- Centro Singular de Investigación en Química Biolóxica e Materiais Moleculares (CIQUS), and Departamento de Química Orgánica, Universidade de Santiago de Compostela Rúa Jenaro de la Fuente s/n 15782 Santiago de Compostela Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST) Baldiri Reixac 10-12 08028 Barcelona Spain
- Department of Biochemistry and Molecular Biology, University of Barcelona 08028 Barcelona Spain
| | - José L Mascareñas
- Centro Singular de Investigación en Química Biolóxica e Materiais Moleculares (CIQUS), and Departamento de Química Orgánica, Universidade de Santiago de Compostela Rúa Jenaro de la Fuente s/n 15782 Santiago de Compostela Spain
| |
Collapse
|
5
|
Bubble Relaxation Dynamics in Homopolymer DNA Sequences. Molecules 2023; 28:molecules28031041. [PMID: 36770707 PMCID: PMC9920605 DOI: 10.3390/molecules28031041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2022] [Revised: 01/11/2023] [Accepted: 01/15/2023] [Indexed: 01/22/2023] Open
Abstract
Understanding the inherent timescales of large bubbles in DNA is critical to a thorough comprehension of its physicochemical characteristics, as well as their potential role on helix opening and biological function. In this work, we employ the coarse-grained Peyrard-Bishop-Dauxois model of DNA to study relaxation dynamics of large bubbles in homopolymer DNA, using simulations up to the microsecond time scale. By studying energy autocorrelation functions of relatively large bubbles inserted into thermalised DNA molecules, we extract characteristic relaxation times from the equilibration process for both adenine-thymine (AT) and guanine-cytosine (GC) homopolymers. Bubbles of different amplitudes and widths are investigated through extensive statistics and appropriate fittings of their relaxation. Characteristic relaxation times increase with bubble amplitude and width. We show that, within the model, relaxation times are two orders of magnitude longer in GC sequences than in AT sequences. Overall, our results confirm that large bubbles leave a lasting impact on the molecule's dynamics, for times between 0.5-500 ns depending on the homopolymer type and bubble shape, thus clearly affecting long-time evolutions of the molecule.
Collapse
|
6
|
Rozza AM, Bakó I, Oláh J. Theoretical insights into water network of B-DNA duplex with Watson-Crick and Hoogsteen base pairing geometries. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119728] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
7
|
Berryman JT, Taghavi A, Mazur F, Tkatchenko A. Quantum machine learning corrects classical forcefields: Stretching DNA base pairs in explicit solvent. J Chem Phys 2022; 157:064107. [PMID: 35963717 DOI: 10.1063/5.0094727] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In order to improve the accuracy of molecular dynamics simulations, classical forcefields are supplemented with a kernel-based machine learning method trained on quantum-mechanical fragment energies. As an example application, a potential-energy surface is generalized for a small DNA duplex, taking into account explicit solvation and long-range electron exchange-correlation effects. A long-standing problem in molecular science is that experimental studies of the structural and thermodynamic behavior of DNA under tension are not well confirmed by simulation; study of the potential energy vs extension taking into account a novel correction shows that leading classical DNA models have excessive stiffness with respect to stretching. This discrepancy is found to be common across multiple forcefields. The quantum correction is in qualitative agreement with the experimental thermodynamics for larger DNA double helices, providing a candidate explanation for the general and long-standing discrepancy between single molecule stretching experiments and classical calculations of DNA stretching. The new dataset of quantum calculations should facilitate multiple types of nucleic acid simulation, and the associated Kernel Modified Molecular Dynamics method (KMMD) is applicable to biomolecular simulations in general. KMMD is made available as part of the AMBER22 simulation software.
Collapse
Affiliation(s)
- Joshua T Berryman
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Amirhossein Taghavi
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Florian Mazur
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
8
|
Diaz A, Jothiraman HB, Ramakrishnan V. Effect of glycerol on free DNA: A molecular dynamics simulation study. J Mol Graph Model 2022; 114:108169. [DOI: 10.1016/j.jmgm.2022.108169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2021] [Revised: 03/16/2022] [Accepted: 03/16/2022] [Indexed: 11/29/2022]
|
9
|
Schneider AL, Albrecht AV, Huang K, Germann MW, Poon GMK. Self-Consistent Parameterization of DNA Residues for the Non-Polarizable AMBER Force Fields. Life (Basel) 2022; 12:life12050666. [PMID: 35629334 PMCID: PMC9143812 DOI: 10.3390/life12050666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Revised: 04/27/2022] [Accepted: 04/27/2022] [Indexed: 11/22/2022] Open
Abstract
Fixed-charge (non-polarizable) forcefields are accurate and computationally efficient tools for modeling the molecular dynamics of nucleic acid polymers, particularly DNA, well into the µs timescale. The continued utility of these forcefields depends in part on expanding the residue set in step with advancing nucleic acid chemistry and biology. A key step in parameterizing new residues is charge derivation which is self-consistent with the existing residues. As atomic charges are derived by fitting against molecular electrostatic potentials, appropriate structural models are critical. Benchmarking against the existing charge set used in current AMBER nucleic acid forcefields, we report that quantum mechanical models of deoxynucleosides, even at a high level of theory, are not optimal structures for charge derivation. Instead, structures from molecular mechanics minimization yield charges with up to 6-fold lower RMS deviation from the published values, due to the choice of such an approach in the derivation of the original charge set. We present a contemporary protocol for rendering self-consistent charges as well as optimized charges for a panel of nine non-canonical residues that will permit comparison with literature as well as studying the dynamics of novel DNA polymers.
Collapse
Affiliation(s)
- Amelia L. Schneider
- Department of Chemistry, Georgia State University, Atlanta, GA 30303, USA; (A.L.S.); (A.V.A.); (K.H.)
| | - Amanda V. Albrecht
- Department of Chemistry, Georgia State University, Atlanta, GA 30303, USA; (A.L.S.); (A.V.A.); (K.H.)
| | - Kenneth Huang
- Department of Chemistry, Georgia State University, Atlanta, GA 30303, USA; (A.L.S.); (A.V.A.); (K.H.)
| | - Markus W. Germann
- Department of Chemistry, Georgia State University, Atlanta, GA 30303, USA; (A.L.S.); (A.V.A.); (K.H.)
- Department of Biology, Georgia State University, Atlanta, GA 30303, USA
- Correspondence: (M.W.G.); (G.M.K.P.)
| | - Gregory M. K. Poon
- Department of Chemistry, Georgia State University, Atlanta, GA 30303, USA; (A.L.S.); (A.V.A.); (K.H.)
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA 30303, USA
- Correspondence: (M.W.G.); (G.M.K.P.)
| |
Collapse
|
10
|
Martin B, Dans PD, Wieczór M, Villegas N, Brun-Heath I, Battistini F, Terrazas M, Orozco M. Molecular basis of Arginine and Lysine DNA sequence-dependent thermo-stability modulation. PLoS Comput Biol 2022; 18:e1009749. [PMID: 35007284 PMCID: PMC8782489 DOI: 10.1371/journal.pcbi.1009749] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Revised: 01/21/2022] [Accepted: 12/13/2021] [Indexed: 11/19/2022] Open
Abstract
We have used a variety of theoretical and experimental techniques to study the role of four basic amino acids-Arginine, Lysine, Ornithine and L-2,4-Diaminobutyric acid-on the structure, flexibility and sequence-dependent stability of DNA. We found that the presence of organic ions stabilizes the duplexes and significantly reduces the difference in stability between AT- and GC-rich duplexes with respect to the control conditions. This suggests that these amino acids, ingredients of the primordial soup during abiogenesis, could have helped to equalize the stability of AT- and GC-rich DNA oligomers, facilitating a general non-catalysed self-replication of DNA. Experiments and simulations demonstrate that organic ions have an effect that goes beyond the general electrostatic screening, involving specific interactions along the grooves of the double helix. We conclude that organic ions, largely ignored in the DNA world, should be reconsidered as crucial structural elements far from mimics of small inorganic cations.
Collapse
Affiliation(s)
- Benjamin Martin
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
- Institute of Bioengineering, School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Pablo D. Dans
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
- Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República (UdelaR), Salto, Uruguay
- Functional Genomics Unit, Institut Pasteur de Montevideo, Montevideo, Uruguay
| | - Milosz Wieczór
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Nuria Villegas
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Isabelle Brun-Heath
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Federica Battistini
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Montserrat Terrazas
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), the Barcelona Institute of Science and Technology, Barcelona, Spain
- Department of Biochemistry and Biomedicine, Faculty of Biology, University of Barcelona, Barcelona, Spain
| |
Collapse
|
11
|
Abstract
:
Clustered regularly interspaced short palindromic repeats along with CRISPR-associated protein
mechanisms preserve the memory of previous experiences with DNA invaders, in particular spacers
that are embedded in CRISPR arrays between coordinate repeats. There has been a fast progression in
the comprehension of this immune system and its implementations; however, there are numerous points
of view that anticipate explanations to make the field an energetic research zone. The efficiency of
CRISPR-Cas depends upon well-considered single guide RNA; for this purpose, many bioinformatics
methods and tools are created to support the design of greatly active and precise single guide RNA. Insilico
single guide RNA architecture is a crucial point for effective gene editing by means of the
CRISPR technique. Persistent attempts have been made to improve in-silico single guide RNA formulation
having great on-target effectiveness and decreased off-target effects. This review offers a summary
of the CRISPR computational tools to help different researchers pick a specific tool for their work according
to pros and cons, along with new thoughts to make new computational tools to overcome all existing
limitations.
Collapse
Affiliation(s)
- Mohsin Ali Nasir
- Center for Informational Biology, University of Electronic Science and Technology of China, No. 2006, Xiyuan Ave,
West Hi-Tech Zone, Chengdu 611731, China
| | - Samia Nawaz
- Center for Informational Biology, University of Electronic Science and Technology of China, No. 2006, Xiyuan Ave,
West Hi-Tech Zone, Chengdu 611731, China
| | - Jian Huang
- Center for Informational Biology, University of Electronic Science and Technology of China, No. 2006, Xiyuan Ave,
West Hi-Tech Zone, Chengdu 611731, China
| |
Collapse
|
12
|
Sequence-dependent structural properties of B-DNA: what have we learned in 40 years? Biophys Rev 2021; 13:995-1005. [DOI: 10.1007/s12551-021-00893-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2021] [Accepted: 11/01/2021] [Indexed: 11/27/2022] Open
|
13
|
Battistini F, Dans PD, Terrazas M, Castellazzi CL, Portella G, Labrador M, Villegas N, Brun-Heath I, González C, Orozco M. The Impact of the HydroxyMethylCytosine epigenetic signature on DNA structure and function. PLoS Comput Biol 2021; 17:e1009547. [PMID: 34748533 PMCID: PMC8601608 DOI: 10.1371/journal.pcbi.1009547] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2021] [Revised: 11/18/2021] [Accepted: 10/10/2021] [Indexed: 12/30/2022] Open
Abstract
We present a comprehensive, experimental and theoretical study of the impact of 5-hydroxymethylation of DNA cytosine. Using molecular dynamics, biophysical experiments and NMR spectroscopy, we found that Ten-Eleven translocation (TET) dioxygenases generate an epigenetic variant with structural and physical properties similar to those of 5-methylcytosine. Experiments and simulations demonstrate that 5-methylcytosine (mC) and 5-hydroxymethylcytosine (hmC) generally lead to stiffer DNA than normal cytosine, with poorer circularization efficiencies and lower ability to form nucleosomes. In particular, we can rule out the hypothesis that hydroxymethylation reverts to unmodified cytosine physical properties, as hmC is even more rigid than mC. Thus, we do not expect dramatic changes in the chromatin structure induced by differences in physical properties between d(mCpG) and d(hmCpG). Conversely, our simulations suggest that methylated-DNA binding domains (MBDs), associated with repression activities, are sensitive to the substitution d(mCpG) ➔ d(hmCpG), while MBD3 which has a dual activation/repression activity is not sensitive to the d(mCpG) d(hmCpG) change. Overall, while gene activity changes due to cytosine methylation are the result of the combination of stiffness-related chromatin reorganization and MBD binding, those associated to 5-hydroxylation of methylcytosine could be explained by a change in the balance of repression/activation pathways related to differential MBD binding. In Eukaryotic cells, DNA epigenetic modifications play an important role in gene expression and regulation, and protein recognition. In this work we investigate the physical implications of cytosine 5-hydroxymethylation on DNA, its structural and flexibility differences with methylated and unmodified cytosine using molecular dynamics, biophysical experiments and NMR spectroscopy. In particular the effect of hydroxyl group on free energy of nucleosome and Methyl binding Protein (MBD) binding, comparing in silico and experimental data to shed light on the effect of the reduced flexibility and the direct protein-DNA recognition.
Collapse
Affiliation(s)
- Federica Battistini
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Pablo D. Dans
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
- Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República (UdelaR), Salto, Uruguay
- Functional Genomics Lab., Institut Pasteur of Montevideo, Montevideo, Uruguay
| | - Montserrat Terrazas
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Chiara L. Castellazzi
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Guillem Portella
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
- Chemistry Department, University of Cambridge, Cambridge, United Kingdom
| | - Mireia Labrador
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Núria Villegas
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Isabelle Brun-Heath
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Carlos González
- Instituto Química Física Rocasolano, Consejo Superior de Investigaciones Científicas (CSIC), Madrid, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology, Barcelona, Spain
- Department of Biochemistry and Molecular Biology, University of Barcelona, Barcelona, Spain
- * E-mail:
| |
Collapse
|
14
|
Schlick T, Portillo-Ledesma S, Myers CG, Beljak L, Chen J, Dakhel S, Darling D, Ghosh S, Hall J, Jan M, Liang E, Saju S, Vohr M, Wu C, Xu Y, Xue E. Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field. Annu Rev Biophys 2021; 50:267-301. [PMID: 33606945 PMCID: PMC8105287 DOI: 10.1146/annurev-biophys-091720-102019] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, New York, New York 10003, USA;
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200122, China
| | | | - Christopher G Myers
- Department of Chemistry, New York University, New York, New York 10003, USA;
| | - Lauren Beljak
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Justin Chen
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sami Dakhel
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Daniel Darling
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sayak Ghosh
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Joseph Hall
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Mikaeel Jan
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Emily Liang
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sera Saju
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Mackenzie Vohr
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Chris Wu
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Yifan Xu
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Eva Xue
- College of Arts and Science, New York University, New York, New York 10003, USA
| |
Collapse
|
15
|
Nierzwicki Ł, Arantes PR, Saha A, Palermo G. Establishing the allosteric mechanism in CRISPR-Cas9. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2021; 11:e1503. [PMID: 34322166 PMCID: PMC8315640 DOI: 10.1002/wcms.1503] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Accepted: 08/26/2020] [Indexed: 12/18/2022]
Abstract
Allostery is a fundamental property of proteins, which regulates biochemical information transfer between spatially distant sites. Here, we report on the critical role of molecular dynamics (MD) simulations in discovering the mechanism of allosteric communication within CRISPR-Cas9, a leading genome editing machinery with enormous promises for medicine and biotechnology. MD revealed how allostery intervenes during at least three steps of the CRISPR-Cas9 function: affecting DNA recognition, mediating the cleavage and interfering with the off-target activity. An allosteric communication that activates concerted DNA cleavages was found to led through the L1/L2 loops, which connect the HNH and RuvC catalytic domains. The identification of these "allosteric transducers" inspired the development of novel variants of the Cas9 protein with improved specificity, opening a new avenue for controlling the CRISPR-Cas9 activity. Discussed studies also highlight the critical role of the recognition lobe in the conformational activation of the catalytic HNH domain. Specifically, the REC3 region was found to modulate the dynamics of HNH by sensing the formation of the RNA:DNA hybrid. The role of REC3 was revealed to be particularly relevant in the presence of DNA mismatches. Indeed, interference of REC3 with the RNA:DNA hybrid containing mismatched pairs at specific positions resulted in locking HNH in an inactive "conformational checkpoint" conformation, thereby hampering off-target cleavages. Overall, MD simulations established the fundamental mechanisms underlying the allosterism of CRISPR-Cas9, aiding engineering strategies to develop new CRISPR-Cas9 variants for improved genome editing.
Collapse
Affiliation(s)
- Łukasz Nierzwicki
- Department of Bioengineering, University of California Riverside, Riverside, California
| | - Pablo Ricardo Arantes
- Department of Bioengineering, University of California Riverside, Riverside, California
| | - Aakash Saha
- Department of Bioengineering, University of California Riverside, Riverside, California
| | - Giulia Palermo
- Department of Bioengineering and Department of Chemistry, University of California Riverside, Riverside, California
| |
Collapse
|
16
|
Schlick T, Portillo-Ledesma S. Biomolecular modeling thrives in the age of technology. NATURE COMPUTATIONAL SCIENCE 2021; 1:321-331. [PMID: 34423314 PMCID: PMC8378674 DOI: 10.1038/s43588-021-00060-9] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 03/22/2021] [Indexed: 12/12/2022]
Abstract
The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, New York, NY, USA
- Courant Institute of Mathematical Sciences, New York University, New York, NY, USA
- New York University–East China Normal University Center for Computational Chemistry at New York University Shanghai, Shanghai, China
| | | |
Collapse
|
17
|
Ibrahim M, Saeed T, Hekmatifar M, Sabetvand R, Chu YM, Toghraie D, Iran TG. The atomic interactions between Histone and 3LPT protein using an equilibrium molecular dynamics simulation. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.115397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
18
|
Mondal M, Yang L, Cai Z, Patra P, Gao YQ. A perspective on the molecular simulation of DNA from structural and functional aspects. Chem Sci 2021; 12:5390-5409. [PMID: 34168783 PMCID: PMC8179617 DOI: 10.1039/d0sc05329e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.
Collapse
Affiliation(s)
- Manas Mondal
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory 518055 Shenzhen China
| | - Lijiang Yang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China
| | - Zhicheng Cai
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China.,Biomedical Pioneering Innovation Center, Peking University 100871 Beijing China
| | - Piya Patra
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory 518055 Shenzhen China .,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China
| | - Yi Qin Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory 518055 Shenzhen China .,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University 100871 Beijing China.,Biomedical Pioneering Innovation Center, Peking University 100871 Beijing China.,Beijing Advanced Innovation Center for Genomics, Peking University 100871 Beijing China
| |
Collapse
|
19
|
Menssen RJ, Kimmel GJ, Tokmakoff A. Investigation into the mechanism and dynamics of DNA association and dissociation utilizing kinetic Monte Carlo simulations. J Chem Phys 2021; 154:045101. [PMID: 33514113 DOI: 10.1063/5.0035187] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
In this work, we present a kinetic Markov state Monte Carlo model designed to complement temperature-jump (T-jump) infrared spectroscopy experiments probing the kinetics and dynamics of short DNA oligonucleotides. The model is designed to be accessible to experimental researchers in terms of both computational simplicity and expense while providing detailed insights beyond those provided by experimental methods. The model is an extension of a thermodynamic lattice model for DNA hybridization utilizing the formalism of the nucleation-zipper mechanism. Association and dissociation trajectories were generated utilizing the Gillespie algorithm and parameters determined via fitting the association and dissociation timescales to previously published experimental data. Terminal end fraying, experimentally observed following a rapid T-jump, in the sequence 5'-ATATGCATAT-3' was replicated by the model that also demonstrated that experimentally observed fast dynamics in the sequences 5'-C(AT)nG-3', where n = 2-6, were also due to terminal end fraying. The dominant association pathways, isolated by transition pathway theory, showed two primary motifs: initiating at or next to a G:C base pair, which is enthalpically favorable and related to the increased strength of G:C base pairs, and initiating in the center of the sequence, which is entropically favorable and related to minimizing the penalty associated with the decrease in configurational entropy due to hybridization.
Collapse
Affiliation(s)
- Ryan J Menssen
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 929 East 57th Street, Chicago, Illinois 60637, USA
| | - Gregory J Kimmel
- Moffitt Cancer Center, 12902 USF Magnolia Drive, Tampa, Florida 33612, USA
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 929 East 57th Street, Chicago, Illinois 60637, USA
| |
Collapse
|
20
|
Lai CL, Chen C, Ou SC, Prentiss M, Pettitt BM. Interactions between identical DNA double helices. Phys Rev E 2021; 101:032414. [PMID: 32289903 DOI: 10.1103/physreve.101.032414] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2019] [Accepted: 01/02/2020] [Indexed: 01/09/2023]
Abstract
The molecular mechanism of specific interactions between double stranded DNA molecules has been investigated for many years. Problems remain in how confinement, ions, and condensing agents change the interactions. We consider how the orientational alignment of DNAs contributes to the interactions via free energy simulations. Here we report on the effective interactions between two parallel DNA double helices in 150-mM NaCl solution using all atom models. We calculate the potential of mean force (PMF) of DNA-DNA interactions as a function of two coordinates, interhelical separation of parallel double helices and relative rotation of a DNA molecule with respect to the other about the helical axis. We generate the two-dimensional PMF to better understand the effective interactions when a DNA molecule is in juxtaposition with another. The analysis of the ion and solvent distributions around the DNA and particularly in the interface region shows that certain alignments of the DNA pair enhance the interactions. At local free energy minima in distance and alignment, water molecules and Na^{+} ions form a hydrogen bonded network with the phosphates from each DNA. This network contributes an attractive energy component to the DNA-DNA interactions. Our results provide a molecular mechanism whereby local DNA-DNA interactions, depending on the helical orientation, give a potential mechanism for stabilizing pairing of much larger lengths of homologous DNA that have been seen experimentally. The study suggests an atomically detailed local picture of relevance to certain aspects of DNA condensation or aggregation.
Collapse
Affiliation(s)
- Chun-Liang Lai
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555, USA
| | - Chuanying Chen
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555, USA
| | - Shu-Ching Ou
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555, USA
| | - Mara Prentiss
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA
| | - B Montgomery Pettitt
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555, USA
| |
Collapse
|
21
|
Afek A, Shi H, Rangadurai A, Sahay H, Senitzki A, Xhani S, Fang M, Salinas R, Mielko Z, Pufall MA, Poon GMK, Haran TE, Schumacher MA, Al-Hashimi HM, Gordân R. DNA mismatches reveal conformational penalties in protein-DNA recognition. Nature 2020; 587:291-296. [PMID: 33087930 PMCID: PMC7666076 DOI: 10.1038/s41586-020-2843-2] [Citation(s) in RCA: 59] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2019] [Accepted: 09/17/2020] [Indexed: 12/17/2022]
Abstract
Transcription factors recognize specific genomic sequences to regulate complex gene-expression programs. Although it is well-established that transcription factors bind to specific DNA sequences using a combination of base readout and shape recognition, some fundamental aspects of protein-DNA binding remain poorly understood1,2. Many DNA-binding proteins induce changes in the structure of the DNA outside the intrinsic B-DNA envelope. However, how the energetic cost that is associated with distorting the DNA contributes to recognition has proven difficult to study, because the distorted DNA exists in low abundance in the unbound ensemble3-9. Here we use a high-throughput assay that we term SaMBA (saturation mismatch-binding assay) to investigate the role of DNA conformational penalties in transcription factor-DNA recognition. In SaMBA, mismatched base pairs are introduced to pre-induce structural distortions in the DNA that are much larger than those induced by changes in the Watson-Crick sequence. Notably, approximately 10% of mismatches increased transcription factor binding, and for each of the 22 transcription factors that were examined, at least one mismatch was found that increased the binding affinity. Mismatches also converted non-specific sites into high-affinity sites, and high-affinity sites into 'super sites' that exhibit stronger affinity than any known canonical binding site. Determination of high-resolution X-ray structures, combined with nuclear magnetic resonance measurements and structural analyses, showed that many of the DNA mismatches that increase binding induce distortions that are similar to those induced by protein binding-thus prepaying some of the energetic cost incurred from deforming the DNA. Our work indicates that conformational penalties are a major determinant of protein-DNA recognition, and reveals mechanisms by which mismatches can recruit transcription factors and thus modulate replication and repair activities in the cell10,11.
Collapse
Affiliation(s)
- Ariel Afek
- Center for Genomic and Computational Biology, Duke University School of Medicine, Durham, NC, USA
- Department of Biostatistics and Bioinformatics, Duke University School of Medicine, Durham, NC, USA
| | - Honglue Shi
- Department of Chemistry, Duke University, Durham, NC, USA
| | - Atul Rangadurai
- Department of Biochemistry, Duke University School of Medicine, Durham, NC, USA
| | - Harshit Sahay
- Center for Genomic and Computational Biology, Duke University School of Medicine, Durham, NC, USA
- Program in Computational Biology and Bioinformatics, Duke University School of Medicine, Durham, NC, USA
| | - Alon Senitzki
- Department of Biology, Technion-Israel Institute of Technology, Haifa, Israel
| | - Suela Xhani
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
| | - Mimi Fang
- Department of Biochemistry, Carver College of Medicine, University of Iowa, Iowa City, IA, USA
- Holden Comprehensive Cancer Center, University of Iowa, Iowa City, IA, USA
| | - Raul Salinas
- Department of Biochemistry, Duke University School of Medicine, Durham, NC, USA
| | - Zachery Mielko
- Center for Genomic and Computational Biology, Duke University School of Medicine, Durham, NC, USA
- Program in Genetics and Genomics, Duke University School of Medicine, Durham, NC, USA
| | - Miles A Pufall
- Department of Biochemistry, Carver College of Medicine, University of Iowa, Iowa City, IA, USA
- Holden Comprehensive Cancer Center, University of Iowa, Iowa City, IA, USA
| | - Gregory M K Poon
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| | - Tali E Haran
- Department of Biology, Technion-Israel Institute of Technology, Haifa, Israel
| | - Maria A Schumacher
- Department of Biochemistry, Duke University School of Medicine, Durham, NC, USA
| | - Hashim M Al-Hashimi
- Department of Chemistry, Duke University, Durham, NC, USA.
- Department of Biochemistry, Duke University School of Medicine, Durham, NC, USA.
| | - Raluca Gordân
- Center for Genomic and Computational Biology, Duke University School of Medicine, Durham, NC, USA.
- Department of Biostatistics and Bioinformatics, Duke University School of Medicine, Durham, NC, USA.
- Department of Computer Science, Duke University, Durham, NC, USA.
- Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, NC, USA.
| |
Collapse
|
22
|
Ray D, Andricioaei I. Free Energy Landscape and Conformational Kinetics of Hoogsteen Base Pairing in DNA vs. RNA. Biophys J 2020; 119:1568-1579. [PMID: 32946766 DOI: 10.1016/j.bpj.2020.08.031] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2019] [Revised: 05/10/2020] [Accepted: 08/25/2020] [Indexed: 10/23/2022] Open
Abstract
Genetic information is encoded in the DNA double helix, which, in its physiological milieu, is characterized by the iconical Watson-Crick nucleo-base pairing. Recent NMR relaxation experiments revealed the transient presence of an alternative, Hoogsteen (HG) base pairing pattern in naked DNA duplexes, and estimated its relative stability and lifetime. In contrast with DNA, such structures were not observed in RNA duplexes. Understanding HG base pairing is important because the underlying "breathing" motion between the two conformations can significantly modulate protein binding. However, a detailed mechanistic insight into the transition pathways and kinetics is still missing. We performed enhanced sampling simulation (with combined metadynamics and adaptive force-bias method) and Markov state modeling to obtain accurate free energy, kinetics, and the intermediates in the transition pathway between Watson-Crick and HG base pairs for both naked B-DNA and A-RNA duplexes. The Markov state model constructed from our unbiased MD simulation data revealed previously unknown complex extrahelical intermediates in the seemingly simple process of base flipping in B-DNA. Extending our calculation to A-RNA, for which HG base pairing is not observed experimentally, resulted in relatively unstable, single-hydrogen-bonded, distorted Hoogsteen-like bases. Unlike B-DNA, the transition pathway primarily involved base paired and intrahelical intermediates with transition timescales much longer than that of B-DNA. The seemingly obvious flip-over reaction coordinate (i.e., the glycosidic torsion angle) is unable to resolve the intermediates. Instead, a multidimensional picture involving backbone dihedral angles and distance between hydrogen bond donor and acceptor atoms is required to gain insight into the molecular mechanism.
Collapse
Affiliation(s)
| | - Ioan Andricioaei
- Department of Chemistry; Department of Physics and Astronomy, University of California Irvine, Irvine, California.
| |
Collapse
|
23
|
The assessment of three dimensional modelling design for single strand DNA aptamers for computational chemistry application. Biophys Chem 2020; 267:106492. [PMID: 33035750 DOI: 10.1016/j.bpc.2020.106492] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Revised: 09/17/2020] [Accepted: 09/29/2020] [Indexed: 11/23/2022]
Abstract
Aptamers are oligonucleotides and peptides around 15-100 bases in length and are suitable as detection probes or as therapeutics molecules. There are growing interests in the aptamer screening approach through computational simulation methods. DNA and RNA modelling lacks of validation on their predicted 3D structures due to less number of validation tools, unlike protein structures. We suggest an approach to design the stem-loop/hairpin for the three dimensional structure of DNA aptamers through serial applications of computational prediction methods by comparing the simulated structures with the experimental data deposited in PDB Data bank, followed by MD simulations. The result shows minimal structural differences were observed between the designed and the original NMR aptamers, and the stem-loop conformational structures were also retained during the MD thus suggesting the proposed aptamers designing methods are able to synthesize a high quality molecular structure of hairpin aptamers, comparable to the NMR structures.
Collapse
|
24
|
Alexiou TS, Mintis DG, Mavrantzas VG. Molecular Dynamics Simulation of the Diffusion Dynamics of Linear DNA Fragments in Dilute Solution with the Parmbsc1 Force Field and Comparison with Experimental Data and Theoretical Models. Macromolecules 2020. [DOI: 10.1021/acs.macromol.0c00222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Terpsichori S. Alexiou
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
| | - Dimitris G. Mintis
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
| | - Vlasis G. Mavrantzas
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
- Department of Mechanical and Process Engineering, Particle Technology Laboratory, ETH Zürich, CH-8092 Zürich, Switzerland
| |
Collapse
|
25
|
Colibactin DNA-damage signature indicates mutational impact in colorectal cancer. Nat Med 2020; 26:1063-1069. [PMID: 32483361 DOI: 10.1038/s41591-020-0908-2] [Citation(s) in RCA: 140] [Impact Index Per Article: 35.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2020] [Accepted: 04/23/2020] [Indexed: 02/06/2023]
Abstract
The mucosal epithelium is a common target of damage by chronic bacterial infections and the accompanying toxins, and most cancers originate from this tissue. We investigated whether colibactin, a potent genotoxin1 associated with certain strains of Escherichia coli2, creates a specific DNA-damage signature in infected human colorectal cells. Notably, the genomic contexts of colibactin-induced DNA double-strand breaks were enriched for an AT-rich hexameric sequence motif, associated with distinct DNA-shape characteristics. A survey of somatic mutations at colibactin target sites of several thousand cancer genomes revealed notable enrichment of this motif in colorectal cancers. Moreover, the exact double-strand-break loci corresponded with mutational hot spots in cancer genomes, reminiscent of a trinucleotide signature previously identified in healthy colorectal epithelial cells3. The present study provides evidence for the etiological role of colibactin in human cancer.
Collapse
|
26
|
Mitchell BP, Hsu RV, Medrano MA, Zewde NT, Narkhede YB, Palermo G. Spontaneous Embedding of DNA Mismatches Within the RNA:DNA Hybrid of CRISPR-Cas9. Front Mol Biosci 2020; 7:39. [PMID: 32258048 PMCID: PMC7093078 DOI: 10.3389/fmolb.2020.00039] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2019] [Accepted: 02/19/2020] [Indexed: 12/15/2022] Open
Abstract
CRISPR-Cas9 is the forefront technology for editing the genome. In this system, the Cas9 protein is programmed with guide RNAs to process DNA sequences that match the guide RNA forming an RNA:DNA hybrid structure. However, the binding of DNA sequences that do not fully match the guide RNA can limit the applicability of CRISPR-Cas9 for genome editing, resulting in the so-called off-target effects. Here, molecular dynamics is used to probe the effect of DNA base pair mismatches within the RNA:DNA hybrid in CRISPR-Cas9. Molecular simulations revealed that the presence of mismatched pairs in the DNA at distal sites with respect to the Protospacer Adjacent Motif (PAM) recognition sequence induces an extended opening of the RNA:DNA hybrid, leading to novel interactions established by the unwound nucleic acids and the protein counterpart. On the contrary, mismatched pairs upstream of the RNA:DNA hybrid are rapidly incorporated within the heteroduplex, with minor effect on the protein-nucleic acid interactions. As a result, mismatched pairs at PAM distal ends interfere with the activation of the catalytic HNH domain, while mismatches fully embedded in the RNA:DNA do not affect the HNH dynamics and enable its activation to cleave the DNA. These findings provide a mechanistic understanding to the intriguing experimental evidence that PAM distal mismatches hamper a proper function of HNH, explaining also why mismatches within the heteroduplex are much more tolerated. This constitutes a step forward in understanding off-target effects in CRISPR-Cas9, which encourages novel structure-based engineering efforts aimed at preventing the onset of off-target effects.
Collapse
Affiliation(s)
- Brandon P. Mitchell
- Department of Bioengineering, University of California, Riverside, Riverside, CA, United States
| | - Rohaine V. Hsu
- Department of Bioengineering, University of California, Riverside, Riverside, CA, United States
| | - Marco A. Medrano
- Department of Bioengineering, University of California, Riverside, Riverside, CA, United States
| | - Nehemiah T. Zewde
- Department of Bioengineering, University of California, Riverside, Riverside, CA, United States
| | - Yogesh B. Narkhede
- Department of Bioengineering, University of California, Riverside, Riverside, CA, United States
| | - Giulia Palermo
- Department of Bioengineering, University of California, Riverside, Riverside, CA, United States
- Department of Chemistry, University of California, Riverside, Riverside, CA, United States
| |
Collapse
|
27
|
Jolfaei NA, Jolfaei NA, Hekmatifar M, Piranfar A, Toghraie D, Sabetvand R, Rostami S. Investigation of thermal properties of DNA structure with precise atomic arrangement via equilibrium and non-equilibrium molecular dynamics approaches. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2020; 185:105169. [PMID: 31715331 DOI: 10.1016/j.cmpb.2019.105169] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/26/2019] [Revised: 10/26/2019] [Accepted: 10/28/2019] [Indexed: 06/10/2023]
Abstract
BACKGROUND AND OBJECTIVE Thermal conductivity of Deoxyribonucleic acid molecules is important for nanotechnology applications. Theoretical simulations based on simple models predict thermal conductivity for these molecular structures. METHODS In this work, we calculate the thermal properties of Deoxyribonucleic acid with precise atomic arrangement via equilibrium and non-equilibrium molecular dynamics approaches. In these methods, each Deoxyribonucleic acid molecule is represented by C, N, O, and P atoms and implemented dreidng potential to describe their atomic interactions. RESULTS Our calculated rate for thermal conductivity via equilibrium and non-equilibrium molecular dynamics methods is 0.381 W/m K and 0.373 W/m K, respectively. By comparing results from these two methods, it was found that the results from equilibrium and non-equilibrium molecular dynamics methods are identical, approximately. On the other hand, the number of DNA molecules and the equilibrium temperature of the simulated structures were important factors in their thermal conductivity rates, and their thermal conductivity was calculated at 0.323 W/m K-0.381 W/m K intervals for equilibrium and 0.303 W/m K-0.373 W/m K interval for non-equilibrium calculations. CONCLUSIONS These results are in good agreement with thermal conductivity calculation with other research groups.
Collapse
Affiliation(s)
| | - Niyusha Adavoodi Jolfaei
- Department of Pharmaceutical Sciences, KLE College of Pharmacy, 2nd block, Rajajinagar, Bengaluru, Karnatka, India
| | - Maboud Hekmatifar
- Department of Mechanical Engineering, Khomeinishahr Branch, Islamic Azad University, Khomeinishahr, Iran
| | - Anahita Piranfar
- Biomechanic Department, Biomedical Engineering Faculty, Mashhad Branch, Islamic Azad University, Mashhad, Iran
| | - Davood Toghraie
- Department of Mechanical Engineering, Khomeinishahr Branch, Islamic Azad University, Khomeinishahr, Iran
| | - Roozbeh Sabetvand
- Department of Energy Engineering and Physics, Faculty of Condensed Matter Physics, Amirkabir University of Technology, Tehran, Iran
| | - Sara Rostami
- Laboratory of Magnetism and Magnetic Materials, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Vietnam; Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam.
| |
Collapse
|
28
|
Sarkar S, Singh PC. Alteration of the groove width of DNA induced by the multimodal hydrogen bonding of denaturants with DNA bases in its grooves affects their stability. Biochim Biophys Acta Gen Subj 2020; 1864:129498. [DOI: 10.1016/j.bbagen.2019.129498] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2019] [Revised: 11/20/2019] [Accepted: 11/25/2019] [Indexed: 02/08/2023]
|
29
|
Cuppari A, Fernández-Millán P, Battistini F, Tarrés-Solé A, Lyonnais S, Iruela G, Ruiz-López E, Enciso Y, Rubio-Cosials A, Prohens R, Pons M, Alfonso C, Tóth K, Rivas G, Orozco M, Solà M. DNA specificities modulate the binding of human transcription factor A to mitochondrial DNA control region. Nucleic Acids Res 2020; 47:6519-6537. [PMID: 31114891 PMCID: PMC6614842 DOI: 10.1093/nar/gkz406] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2018] [Revised: 04/30/2019] [Accepted: 05/15/2019] [Indexed: 12/18/2022] Open
Abstract
Human mitochondrial DNA (h-mtDNA) codes for 13 subunits of the oxidative phosphorylation pathway, the essential route that produces ATP. H-mtDNA transcription and replication depends on the transcription factor TFAM, which also maintains and compacts this genome. It is well-established that TFAM activates the mtDNA promoters LSP and HSP1 at the mtDNA control region where DNA regulatory elements cluster. Previous studies identified still uncharacterized, additional binding sites at the control region downstream from and slightly similar to LSP, namely sequences X and Y (Site-X and Site-Y) (Fisher et al., Cell 50, pp 247-258, 1987). Here, we explore TFAM binding at these two sites and compare them to LSP by multiple experimental and in silico methods. Our results show that TFAM binding is strongly modulated by the sequence-dependent properties of Site-X, Site-Y and LSP. The high binding versatility of Site-Y or the considerable stiffness of Site-X tune TFAM interactions. In addition, we show that increase in TFAM/DNA complex concentration induces multimerization, which at a very high concentration triggers disruption of preformed complexes. Therefore, our results suggest that mtDNA sequences induce non-uniform TFAM binding and, consequently, direct an uneven distribution of TFAM aggregation sites during the essential process of mtDNA compaction.
Collapse
Affiliation(s)
- Anna Cuppari
- Structural MitoLab, Structural Biology Department, Maria de Maeztu Unit of Excellence, Molecular Biology Institute Barcelona (IBMB-CSIC), 08028 Barcelona, Spain
| | - Pablo Fernández-Millán
- Structural MitoLab, Structural Biology Department, Maria de Maeztu Unit of Excellence, Molecular Biology Institute Barcelona (IBMB-CSIC), 08028 Barcelona, Spain
| | - Federica Battistini
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Aleix Tarrés-Solé
- Structural MitoLab, Structural Biology Department, Maria de Maeztu Unit of Excellence, Molecular Biology Institute Barcelona (IBMB-CSIC), 08028 Barcelona, Spain
| | - Sébastien Lyonnais
- Structural MitoLab, Structural Biology Department, Maria de Maeztu Unit of Excellence, Molecular Biology Institute Barcelona (IBMB-CSIC), 08028 Barcelona, Spain
| | - Guillermo Iruela
- BioNMR Laboratory, Inorganic and Organic Chemistry Department, Universitat de Barcelona, 08028 Barcelona, Spain
| | - Elena Ruiz-López
- Structural MitoLab, Structural Biology Department, Maria de Maeztu Unit of Excellence, Molecular Biology Institute Barcelona (IBMB-CSIC), 08028 Barcelona, Spain
| | - Yuliana Enciso
- Structural MitoLab, Structural Biology Department, Maria de Maeztu Unit of Excellence, Molecular Biology Institute Barcelona (IBMB-CSIC), 08028 Barcelona, Spain
| | - Anna Rubio-Cosials
- Structural MitoLab, Structural Biology Department, Maria de Maeztu Unit of Excellence, Molecular Biology Institute Barcelona (IBMB-CSIC), 08028 Barcelona, Spain
| | - Rafel Prohens
- Unitat de Polimorfisme i Calorimetria, Centres Científics i Tecnològics, University of Barcelona, 08028 Barcelona, Spain
| | - Miquel Pons
- BioNMR Laboratory, Inorganic and Organic Chemistry Department, Universitat de Barcelona, 08028 Barcelona, Spain
| | - Carlos Alfonso
- Centro de Investigaciones Biológicas, Consejo Superior de Investigaciones Científicas (CSIC), 28040 Madrid, Spain
| | - Katalin Tóth
- Deutsches Krebsforschungszentrum, Division Biophysics of Macromolecules, Heidelberg, Germany
| | - Germán Rivas
- Centro de Investigaciones Biológicas, Consejo Superior de Investigaciones Científicas (CSIC), 28040 Madrid, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain.,Department of Biochemistry and Biomedicine, University of Barcelona, Barcelona 08028, Spain
| | - Maria Solà
- Structural MitoLab, Structural Biology Department, Maria de Maeztu Unit of Excellence, Molecular Biology Institute Barcelona (IBMB-CSIC), 08028 Barcelona, Spain
| |
Collapse
|
30
|
Paredes-Rosan CA, Valencia DE, Barazorda-Ccahuana HL, Aguilar-Pineda JA, Gómez B. Amyloid beta oligomers: how pH influences over trimer and pentamer structures? J Mol Model 2019; 26:1. [PMID: 31834477 DOI: 10.1007/s00894-019-4247-5] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2019] [Accepted: 10/30/2019] [Indexed: 02/07/2023]
Abstract
The aggregation of proteins in the brain is one of the main features of neurodegenerative diseases. In Alzheimer's disease, the abnormal aggregation of Aβ-42 is due to intrinsic and extrinsic factors. The latter is due to variations in the environment, such as temperature, salt concentration, and pH. We evaluated the effect of protonation/deprotonation of residues that are part of trimeric and pentameric oligomers at pH 5, pH 6, and pH 7. Molecular dynamics simulation at 200 ns in the canonical ensemble was implemented. The results have revealed that histidine, glutamic acid, and aspartic acid residues showed a protonation/deprotonation effect in oligomers. The root mean square deviation analysis was used to analyze the structural stability at different pHs. We found an increase in hydrophobicity in the side chains of the trimer, while in the pentamer, the structural instability of a compact structure at pH 5 caused the hydrophobic core to open, revealing the hydrophobic region to the environment. At this point, we believe that conformational changes mediated by pH are essential in the aggregation of Aβ-42 oligomers.
Collapse
Affiliation(s)
- Carla A Paredes-Rosan
- Centro de Investigación en Ingeniería Molecular - CIIM, Vicerrectorado de Investigación, Universidad Católica de Santa María, Samuel Velarde 320, Arequipa, Peru
| | - Diego E Valencia
- Centro de Investigación en Ingeniería Molecular - CIIM, Vicerrectorado de Investigación, Universidad Católica de Santa María, Samuel Velarde 320, Arequipa, Peru
| | - Haruna L Barazorda-Ccahuana
- Centro de Investigación en Ingeniería Molecular - CIIM, Vicerrectorado de Investigación, Universidad Católica de Santa María, Samuel Velarde 320, Arequipa, Peru
| | - Jorge A Aguilar-Pineda
- Centro de Investigación en Ingeniería Molecular - CIIM, Vicerrectorado de Investigación, Universidad Católica de Santa María, Samuel Velarde 320, Arequipa, Peru
| | - Badhin Gómez
- Centro de Investigación en Ingeniería Molecular - CIIM, Vicerrectorado de Investigación, Universidad Católica de Santa María, Samuel Velarde 320, Arequipa, Peru. .,Facultad de Ciencias Farmacéuticas, Bioquímicas y Biotecnológicas, Departamento de Farmacia, Bioquímica y Biotecnología, Universidad Católica de Santa María, Samuel Velarde 320, Arequipa, Peru.
| |
Collapse
|
31
|
Dans PD, Balaceanu A, Pasi M, Patelli AS, Petkevičiūtė D, Walther J, Hospital A, Bayarri G, Lavery R, Maddocks JH, Orozco M. The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules. Nucleic Acids Res 2019; 47:11090-11102. [PMID: 31624840 PMCID: PMC6868377 DOI: 10.1093/nar/gkz905] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2019] [Revised: 09/25/2019] [Accepted: 10/06/2019] [Indexed: 12/12/2022] Open
Abstract
We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence. We characterize the sequence-dependent choreography of backbone and base movements modulating the non-Gaussian or anharmonic effects manifested in the higher moments of the dynamics of the duplex when sampling the equilibrium distribution. Contrary to prior assumptions, such anharmonic deformations are not rare in DNA and can play a significant role in determining DNA conformation within complexes. Polymorphisms in helical geometries are particularly prevalent for certain tetranucleotide sequence contexts and are always coupled to a complex network of coordinated changes in the backbone. The analysis of our simulations, which contain instances of all tetranucleotide sequences, allow us to extend Calladine-Dickerson rules used for decades to interpret the average geometry of DNA, leading to a set of rules with quantitative predictive power that encompass nonlocal sequence-dependence and anharmonic fluctuations.
Collapse
Affiliation(s)
- Pablo D Dans
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
- Department of Biological Sciences, University of the Republic (UdelaR), CENUR Gral. Rivera 1350, 50000 Salto, Uruguay
| | - Alexandra Balaceanu
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
| | - Marco Pasi
- LBPA, École normale supérieure Paris-Saclay, 61 Av. du Pdt Wilson, Cachan 94235, France
- Bases Moléculaires et Structurales des Systèmes Infectieux, Univ. Lyon I/CNRS UMR 5086, IBCP, 7 Passage du Vercors, Lyon 69367, France
| | - Alessandro S Patelli
- Institute of Mathematics, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Daiva Petkevičiūtė
- Institute of Mathematics, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
- Faculty of Mathematics and Natural Sciences, Kaunas University of Technology, Studentų g. 50, 51368 Kaunas, Lithuania
| | - Jürgen Walther
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
| | - Genís Bayarri
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
| | - Richard Lavery
- Bases Moléculaires et Structurales des Systèmes Infectieux, Univ. Lyon I/CNRS UMR 5086, IBCP, 7 Passage du Vercors, Lyon 69367, France
| | - John H Maddocks
- Institute of Mathematics, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
- Department of Biochemistry and Molecular Biology. University of Barcelona, 08028 Barcelona, Spain
| |
Collapse
|
32
|
Owens EA, Huynh HT, Stroeva EM, Barman A, Ziabrev K, Paul A, Nguyen SV, Laramie M, Hamelberg D, Germann MW, Wilson WD, Henary M. Second Generation G-Quadruplex Stabilizing Trimethine Cyanines. Bioconjug Chem 2019; 30:2647-2663. [PMID: 31518105 DOI: 10.1021/acs.bioconjchem.9b00571] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
G-Quadruplex DNA has been recognized as a highly appealing target for the development of new selective chemotherapeutics, which could result in markedly reduced toxicity toward normal cells. In particular, the cyanine dyes that bind selectively to G-quadruplex structures without targeting duplex DNA have attracted attention due to their high amenability to structural modifications that allows fine-tuning of their biomolecular interactions. We have previously reported pentamethine and symmetric trimethine cyanines designed to effectively bind G-quadruplexes through end stacking interactions. Herein, we are reporting a second generation of drug candidates, the asymmetric trimethine cyanines. These have been synthesized and evaluated for their quadruplex binding properties. Incorporating a benz[c,d]indolenine heterocyclic unit increased overall quadruplex binding, and elongating the alkyl length increases the quadruplex-to-duplex binding specificity.
Collapse
Affiliation(s)
- Eric A Owens
- Department of Chemistry , Georgia State University , Petit Science Center, 100 Piedmont Ave SE. Atlanta Georgia 30303 , United States.,Center for Diagnostics and Therapeutics , Georgia State University , Petit Science Center, 100 Piedmont Ave SE , Atlanta , Georgia 30303 , United States
| | - Hang T Huynh
- Department of Chemistry , Georgia State University , Petit Science Center, 100 Piedmont Ave SE. Atlanta Georgia 30303 , United States
| | - Ekaterina M Stroeva
- Department of Chemistry , Georgia State University , Petit Science Center, 100 Piedmont Ave SE. Atlanta Georgia 30303 , United States
| | | | - Kostiantyn Ziabrev
- Department of Chemistry , Georgia State University , Petit Science Center, 100 Piedmont Ave SE. Atlanta Georgia 30303 , United States
| | | | | | | | - Donald Hamelberg
- Department of Chemistry , Georgia State University , Petit Science Center, 100 Piedmont Ave SE. Atlanta Georgia 30303 , United States.,Center for Diagnostics and Therapeutics , Georgia State University , Petit Science Center, 100 Piedmont Ave SE , Atlanta , Georgia 30303 , United States
| | - Markus W Germann
- Department of Chemistry , Georgia State University , Petit Science Center, 100 Piedmont Ave SE. Atlanta Georgia 30303 , United States.,Department of Biology , Georgia State University , Petit Science Center, 100 Piedmont Ave. , Atlanta , Georgia 30303 , United States
| | - W David Wilson
- Department of Chemistry , Georgia State University , Petit Science Center, 100 Piedmont Ave SE. Atlanta Georgia 30303 , United States.,Center for Diagnostics and Therapeutics , Georgia State University , Petit Science Center, 100 Piedmont Ave SE , Atlanta , Georgia 30303 , United States
| | - Maged Henary
- Department of Chemistry , Georgia State University , Petit Science Center, 100 Piedmont Ave SE. Atlanta Georgia 30303 , United States.,Center for Diagnostics and Therapeutics , Georgia State University , Petit Science Center, 100 Piedmont Ave SE , Atlanta , Georgia 30303 , United States
| |
Collapse
|
33
|
Sarkar S, Singh PC. Mechanistic Aspects of Fungicide-Induced DNA Damage: Spectroscopic and Molecular Dynamics Simulation Studies. J Phys Chem B 2019; 123:8653-8661. [DOI: 10.1021/acs.jpcb.9b06009] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Sunipa Sarkar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Kolkata 700032, India
| | - Prashant Chandra Singh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Kolkata 700032, India
| |
Collapse
|
34
|
How B-DNA Dynamics Decipher Sequence-Selective Protein Recognition. J Mol Biol 2019; 431:3845-3859. [DOI: 10.1016/j.jmb.2019.07.021] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2019] [Revised: 07/09/2019] [Accepted: 07/10/2019] [Indexed: 11/23/2022]
|
35
|
Sun Z, Wang X, Zhang JZH, He Q. Sulfur-substitution-induced base flipping in the DNA duplex. Phys Chem Chem Phys 2019; 21:14923-14940. [PMID: 31233058 DOI: 10.1039/c9cp01989h] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Base flipping is widely observed in a number of important biological processes. The genetic codes deposited inside the DNA duplex become accessible to external agents upon base flipping. The sulfur substitution of guanine leads to thioguanine, which alters the thermodynamic stability of the GC base pairs and the GT mismatches. Experimental studies conclude that the sulfur substitution decreases the lifetime of the GC base pair. In this work, under three AMBER force fields for nucleotide systems, we firstly performed equilibrium and nonequilibrium free energy simulations to investigate the variation of the thermodynamic profiles in base flipping upon sulfur substitution. It is found that the bsc0 modification, the bsc1 modification and the OL15 modification of AMBER force fields are able to qualitatively describe the sulfur-substitution dependent behavior of the thermodynamics. However, only the two last-generation AMBER force fields are able to provide quantitatively correct predictions. The second computational study on the sulfur substitutions focused on the relative stability of the S6G-C base pair and the S6G-T mismatch. Two conflicting experimental observations were reported by the same authors. One suggested that the S6G-C base pair was more stable, while the other concludes that the S6G-T mismatch was more stable. We answered this question by constructing the free energy profiles along the base flipping pathway computationally.
Collapse
Affiliation(s)
- Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China and Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich 52425, Germany.
| | - Xiaohui Wang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China and Institute of Computational Science, Università della Svizzera italiana (USI), Via Giuseppe Buffi 13, CH-6900, Lugano, Ticino, Switzerland
| | - John Z H Zhang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China and NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China and Department of Chemistry, New York University, NY, NY 10003, USA
| | - Qiaole He
- Forschungszentrum Jülich GmbH, IBG-1: Biotechnology, Wilhelm-Johnen-Str. 1, 52425 Jülich, Germany. and State Key Laboratory of Bioreactor Engineering, R&D Center of Separation and Extraction Technology in Fermentation Industry, East China University of Science and Technology, Shanghai 200237, China
| |
Collapse
|
36
|
Wagner JR, Demir Ö, Carpenter MA, Aihara H, Harki DA, Harris RS, Amaro RE. Determinants of Oligonucleotide Selectivity of APOBEC3B. J Chem Inf Model 2019; 59:2264-2273. [PMID: 30130104 PMCID: PMC6644697 DOI: 10.1021/acs.jcim.8b00427] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
APOBEC3B (A3B) is a prominent source of mutation in many cancers. To date, it has been difficult to capture the native protein-DNA interactions that confer A3B's substrate specificity by crystallography due to the highly dynamic nature of wild-type A3B active site. We use computational tools to restore a recent crystal structure of a DNA-bound A3B C-terminal domain mutant construct to its wild type sequence, and run molecular dynamics simulations to study its substrate recognition mechanisms. Analysis of these simulations reveal dynamics of the native A3Bctd-oligonucleotide interactions, including the experimentally inaccessible loop 1-oligonucleotide interactions. A second series of simulations in which the target cytosine nucleotide was computationally mutated from a deoxyribose to a ribose show a change in sugar ring pucker, leading to a rearrangement of the binding site and revealing a potential intermediate in the binding pathway. Finally, apo simulations of A3B, starting from the DNA-bound open state, experience a rapid and consistent closure of the binding site, reaching conformations incompatible with substrate binding. This study reveals a more realistic and dynamic view of the wild type A3B binding site and provides novel insights for structure-guided design efforts for A3B.
Collapse
Affiliation(s)
- Jeffrey R Wagner
- Department of Chemistry and Biochemistry , University of California, San Diego , La Jolla , California 92093-0340 , United States
| | - Özlem Demir
- Department of Chemistry and Biochemistry , University of California, San Diego , La Jolla , California 92093-0340 , United States
| | - Michael A Carpenter
- Department of Biochemistry, Molecular Biology and Biophysics , University of Minnesota , Minneapolis , Minnesota 55455 , United States
- Masonic Cancer Center , University of Minnesota , Minneapolis , Minnesota 55455 , United States
- Institute for Molecular Virology , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| | - Hideki Aihara
- Department of Biochemistry, Molecular Biology and Biophysics , University of Minnesota , Minneapolis , Minnesota 55455 , United States
- Masonic Cancer Center , University of Minnesota , Minneapolis , Minnesota 55455 , United States
- Institute for Molecular Virology , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| | - Daniel A Harki
- Department of Medicinal Chemistry , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| | - Reuben S Harris
- Department of Biochemistry, Molecular Biology and Biophysics , University of Minnesota , Minneapolis , Minnesota 55455 , United States
- Masonic Cancer Center , University of Minnesota , Minneapolis , Minnesota 55455 , United States
- Institute for Molecular Virology , University of Minnesota , Minneapolis , Minnesota 55455 , United States
- Howard Hughes Medical Institute , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry , University of California, San Diego , La Jolla , California 92093-0340 , United States
| |
Collapse
|
37
|
Ricci C, Chen JS, Miao Y, Jinek M, Doudna JA, McCammon JA, Palermo G. Deciphering Off-Target Effects in CRISPR-Cas9 through Accelerated Molecular Dynamics. ACS CENTRAL SCIENCE 2019; 5:651-662. [PMID: 31041385 PMCID: PMC6487449 DOI: 10.1021/acscentsci.9b00020] [Citation(s) in RCA: 81] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/09/2019] [Indexed: 05/02/2023]
Abstract
CRISPR-Cas9 is the state-of-the-art technology for editing and manipulating nucleic acids. However, the occurrence of off-target mutations can limit its applicability. Here, all-atom enhanced molecular dynamics (MD) simulations-using Gaussian accelerated MD (GaMD)-are used to decipher the mechanism of off-target binding at the molecular level. GaMD reveals that base pair mismatches in the target DNA at distal sites with respect to the protospacer adjacent motif (PAM) can induce an extended opening of the RNA:DNA heteroduplex, which leads to newly formed interactions between the unwound DNA and the L2 loop of the catalytic HNH domain. These conserved interactions constitute a "lock" effectively decreasing the conformational freedom of the HNH domain and hampering its activation for cleavage. Remarkably, depending on their positions at PAM distal sites, DNA mismatches responsible for off-target cleavages are unable to "lock" the HNH domain, thereby leading to the unselective cleavage of DNA sequences. In consistency with the available experimental data, the ability to "lock" the catalytic HNH domain in an inactive "conformational checkpoint" is shown to be a key determinant in the onset of off-target effects. This mechanistic rationale contributes in clarifying a long lasting open issue in the CRISPR-Cas9 function and poses the foundation for designing novel and more specific Cas9 variants, which could be obtained by magnifying the "locking" interactions between HNH and the target DNA in the presence of any incorrect off-target sequence, thus preventing undesired cleavages.
Collapse
Affiliation(s)
- Clarisse
G. Ricci
- Department
of Pharmacology, Department of Chemistry and Biochemistry, and National Biomedical
Computation Resource, University of California
San Diego, La Jolla, California 92093, United States
| | - Janice S. Chen
- Department of Molecular
and Cell Biology, Department of Chemistry, Howard Hughes Medical Institute, Innovative Genomics
Institute, and Molecular Biophysics and Integrated Bioimaging Division, Lawrence
Berkeley National Laboratory, University
of California Berkeley, Berkeley, California 94720, United States
| | - Yinglong Miao
- Center
for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Martin Jinek
- Department
of Biochemistry, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Jennifer A. Doudna
- Department of Molecular
and Cell Biology, Department of Chemistry, Howard Hughes Medical Institute, Innovative Genomics
Institute, and Molecular Biophysics and Integrated Bioimaging Division, Lawrence
Berkeley National Laboratory, University
of California Berkeley, Berkeley, California 94720, United States
| | - J. Andrew McCammon
- Department
of Pharmacology, Department of Chemistry and Biochemistry, and National Biomedical
Computation Resource, University of California
San Diego, La Jolla, California 92093, United States
| | - Giulia Palermo
- Department
of Bioengineering, Bourns College of Engineering, University of California Riverside, 900 University Avenue, Riverside, California 92521, United States
| |
Collapse
|
38
|
Pant P, Pathak A, Jayaram B. Symmetrization of the backbone of nucleic acids: a molecular dynamics study. J Biomol Struct Dyn 2019; 38:673-681. [DOI: 10.1080/07391102.2019.1585292] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Pradeep Pant
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi, India
- Supercomputing Facility for Bioinformatics & Computational Biology, New Delhi, India
| | - Amita Pathak
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi, India
- Supercomputing Facility for Bioinformatics & Computational Biology, New Delhi, India
| | - B. Jayaram
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi, India
- Supercomputing Facility for Bioinformatics & Computational Biology, New Delhi, India
- Kusuma School of Biological Sciences, Indian Institute of Technology Delhi, New Delhi, India
| |
Collapse
|
39
|
Jakubec D, Vondrášek J. Can All-Atom Molecular Dynamics Simulations Quantitatively Describe Homeodomain-DNA Binding Equilibria? J Chem Theory Comput 2019; 15:2635-2648. [PMID: 30807142 DOI: 10.1021/acs.jctc.8b01144] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
We systematically investigate the applicability of a molecular dynamics-based setup for the calculations of standard binding free energies of biologically relevant protein-DNA complexes. The free energies are extracted from a potential of mean force calculated using umbrella sampling simulations. Two protein-DNA systems derived from a homeodomain transcription factor complex are studied in order to investigate the binding of both disordered and globular proteins. Free energies and trajectories obtained using two modern molecular mechanical force fields are compared to each other and to experimental data. The temperature dependence of the calculated standard binding free energies is investigated by performing all simulations over a range of temperatures. We show that the values of standard binding free energies obtained from these simulations are overestimated compared to experimental results. Significant differences are observed between the two protein-DNA systems and between the two force fields, which are explained by different propensities to form inter- and intramolecular contacts. The number of protein-DNA contacts increases with increasing temperature, in agreement with the experimentally known temperature dependence of enthalpies of binding. However, conclusions about the temperature dependence of the standard binding free energies cannot be made with confidence, as the differences among the values are on the order of statistical uncertainty.
Collapse
Affiliation(s)
- David Jakubec
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry , Czech Academy of Sciences , 166 10 Praha 6, Czech Republic.,Department of Physical and Macromolecular Chemistry, Faculty of Science , Charles University , 128 43 Praha 2, Czech Republic
| | - Jiří Vondrášek
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry , Czech Academy of Sciences , 166 10 Praha 6, Czech Republic
| |
Collapse
|
40
|
Affiliation(s)
- Saumyak Mukherjee
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| | - Sayantan Mondal
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| | - Subhajit Acharya
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| | - Biman Bagchi
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
41
|
Lippens JL, Ranganathan SV, D'Esposito RJ, Fabris D. Modular calibrant sets for the structural analysis of nucleic acids by ion mobility spectrometry mass spectrometry. Analyst 2018; 141:4084-99. [PMID: 27152369 DOI: 10.1039/c6an00453a] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
This study explored the use of modular nucleic acid (NA) standards to generate calibration curves capable of translating primary ion mobility readouts into corresponding collision cross section (CCS) data. Putative calibrants consisted of single- (ss) and double-stranded (ds) oligo-deoxynucleotides reaching up to ∼40 kDa in size (i.e., 64 bp) and ∼5700 Å(2) in CCS. To ensure self-consistency among reference CCS values, computational data obtained in house were preferred to any experimental or computational data from disparate sources. Such values were obtained by molecular dynamics (MD) simulations and either the exact hard sphere scattering (EHSS) or the projection superposition approximation (PSA) methods, and then plotted against the corresponding experimental values to generate separate calibration curves. Their performance was evaluated on the basis of their correlation coefficients and ability to provide values that matched the CCS of selected test samples mimicking typical unknowns. The results indicated that the predictive power benefited from the exclusion of higher charged species that were more susceptible to the destabilizing effects of Coulombic repulsion. The results revealed discrepancies between EHSS and PSA data that were ascribable to the different approximations used to describe the ion mobility process. Within the boundaries defined by these approximations and the challenges of modeling NA structure in a solvent-free environment, the calibrant sets enabled the experimental determination of CCS with excellent reproducibility (precision) and error (accuracy), which will support the analysis of progressively larger NA samples of biological significance.
Collapse
Affiliation(s)
| | | | | | - Daniele Fabris
- University at Albany, Albany, New York, USA. and SUNY, Albany, The RNA Institute, 1400 Washington Avenue, Albany, New York, USA
| |
Collapse
|
42
|
Rubio‐Magnieto J, Kajouj S, Di Meo F, Fossépré M, Trouillas P, Norman P, Linares M, Moucheron C, Surin M. Binding Modes and Selectivity of Ruthenium Complexes to Human Telomeric DNA G‐Quadruplexes. Chemistry 2018; 24:15577-15588. [DOI: 10.1002/chem.201802147] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2018] [Revised: 07/12/2018] [Indexed: 02/05/2023]
Affiliation(s)
- Jenifer Rubio‐Magnieto
- Laboratory for Chemistry of Novel Materials Center for Innovation in Materials and Polymers University of Mons-UMONS 20 Place du Parc 7000 Mons Belgium
- Current address: Bioinspired Supramolecular Chemistry and Materials group Departament de Química Inorgànica i Orgànica Universitat Jaume I Avda Sos Baynat s/n E-12071 Castelló Spain
| | - Sofia Kajouj
- Chimie Organique et Photochimie CP160/08 Université libre de Bruxelles 50 avenue F. D. Roosevelt 1050 Bruxelles Belgium
| | - Florent Di Meo
- INSERM U1248 IPPRITT University of Limoges School of Pharmacy 2 rue du Dr. Marcland 87025 Limoges France
| | - Mathieu Fossépré
- Laboratory for Chemistry of Novel Materials Center for Innovation in Materials and Polymers University of Mons-UMONS 20 Place du Parc 7000 Mons Belgium
| | - Patrick Trouillas
- INSERM U1248 IPPRITT University of Limoges School of Pharmacy 2 rue du Dr. Marcland 87025 Limoges France
- RCPTM Palacký University Faculty of Sciences Šlechtitelů 27 78371 Olomouc Czech Republic
| | - Patrick Norman
- Department of Theoretical Chemistry and Biology School of Engineering Sciences in Chemistry Biotechnology and Health KTH Royal Institute of Technology SE-106 91 Stockholm Sweden
| | - Mathieu Linares
- Department of Theoretical Chemistry and Biology School of Engineering Sciences in Chemistry Biotechnology and Health KTH Royal Institute of Technology SE-106 91 Stockholm Sweden
- Swedish e-Science Research Centre (SeRC) KTH Royal Institute of Technology 104 50 Stockholm Sweden
| | - Cécile Moucheron
- Chimie Organique et Photochimie CP160/08 Université libre de Bruxelles 50 avenue F. D. Roosevelt 1050 Bruxelles Belgium
| | - Mathieu Surin
- Laboratory for Chemistry of Novel Materials Center for Innovation in Materials and Polymers University of Mons-UMONS 20 Place du Parc 7000 Mons Belgium
| |
Collapse
|
43
|
Balaceanu A, Pérez A, Dans PD, Orozco M. Allosterism and signal transfer in DNA. Nucleic Acids Res 2018; 46:7554-7565. [PMID: 29905860 PMCID: PMC6125689 DOI: 10.1093/nar/gky549] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2018] [Revised: 05/11/2018] [Accepted: 06/06/2018] [Indexed: 12/14/2022] Open
Abstract
We analysed the basic mechanisms of signal transmission in DNA and the origins of the allostery exhibited by systems such as the ternary complex BAMHI-DNA-GRDBD. We found that perturbation information generated by a primary protein binding event travels as a wave to distant regions of DNA following a hopping mechanism. However, such a structural perturbation is transient and does not lead to permanent changes in the DNA geometry and interaction properties at the secondary binding site. The BAMHI-DNA-GRDBD allosteric mechanism does not occur through any traditional models: direct (protein-protein), indirect (reorganization of the secondary site) readout or solvent-release. On the contrary, it is generated by a subtle and less common entropy-mediated mechanism, which might have an important role to explain other DNA-mediated cooperative effects.
Collapse
Affiliation(s)
- Alexandra Balaceanu
- Joint IRB-BSC Program on Computational Biology, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028 Barcelona, Spain
| | - Alberto Pérez
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA
| | - Pablo D Dans
- Joint IRB-BSC Program on Computational Biology, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028 Barcelona, Spain
| | - Modesto Orozco
- Joint IRB-BSC Program on Computational Biology, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028 Barcelona, Spain
- Department of Biochemistry and Biomedicine, University of Barcelona, 08028 Barcelona, Spain
| |
Collapse
|
44
|
Garai A, Ghoshdastidar D, Senapati S, Maiti PK. Ionic liquids make DNA rigid. J Chem Phys 2018; 149:045104. [PMID: 30068211 DOI: 10.1063/1.5026640] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
Persistence length of double-stranded DNA (dsDNA) is known to decrease with an increase in ionic concentration of the solution. In contrast to this, here we show that the persistence length of dsDNA increases dramatically as a function of ionic liquid (IL) concentration. Using all atom explicit solvent molecular dynamics simulations and theoretical models, we present, for the first time, a systematic study to determine the mechanical properties of dsDNA in various hydrated ILs at different concentrations. We find that dsDNA in 50 wt % ILs have lower persistence length and stretch modulus in comparison to 80 wt % ILs. We further observe that both the persistence length and stretch modulus of dsDNA increase as we increase the concentration of ILs. The present trend of the stretch modulus and persistence length of dsDNA with IL concentration supports the predictions of the macroscopic elastic theory, in contrast to the behavior exhibited by dsDNA in monovalent salt. Our study further suggests the preferable ILs that can be used for maintaining DNA stability during long-term storage.
Collapse
Affiliation(s)
- Ashok Garai
- Department of Physics, Centre for Condensed Matter Theory, Indian Institute of Science, Bangalore 560012, India
| | - Debostuti Ghoshdastidar
- Bhupat and Jyoti Mehta School of Biosciences, Department of Biotechnology, Indian Institute of Technology Madras, Chennai 600036, India
| | - Sanjib Senapati
- Bhupat and Jyoti Mehta School of Biosciences, Department of Biotechnology, Indian Institute of Technology Madras, Chennai 600036, India
| | - Prabal K Maiti
- Department of Physics, Centre for Condensed Matter Theory, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
45
|
Ghoshdastidar D, Bansal M. Dynamics of physiologically relevant noncanonical DNA structures: an overview from experimental and theoretical studies. Brief Funct Genomics 2018; 18:192-204. [DOI: 10.1093/bfgp/ely026] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2018] [Revised: 06/23/2018] [Accepted: 07/09/2018] [Indexed: 12/23/2022] Open
Abstract
Abstract
DNA is a complex molecule with phenomenal inherent plasticity and the ability to form different hydrogen bonding patterns of varying stabilities. These properties enable DNA to attain a variety of structural and conformational polymorphic forms. Structurally, DNA can exist in single-stranded form or as higher-order structures, which include the canonical double helix as well as the noncanonical duplex, triplex and quadruplex species. Each of these structural forms in turn encompasses an ensemble of dynamically heterogeneous conformers depending on the sequence composition and environmental context. In vivo, the widely populated canonical B-DNA attains these noncanonical polymorphs during important cellular processes. While several investigations have focused on the structure of these noncanonical DNA, studying their dynamics has remained nontrivial. Here, we outline findings from some recent advanced experimental and molecular simulation techniques that have significantly contributed toward understanding the complex dynamics of physiologically relevant noncanonical forms of DNA.
Collapse
Affiliation(s)
| | - Manju Bansal
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore 560 012, India
| |
Collapse
|
46
|
|
47
|
Sweeny EA, Singh AB, Chakravarti R, Martinez-Guzman O, Saini A, Haque MM, Garee G, Dans PD, Hannibal L, Reddi AR, Stuehr DJ. Glyceraldehyde-3-phosphate dehydrogenase is a chaperone that allocates labile heme in cells. J Biol Chem 2018; 293:14557-14568. [PMID: 30012884 DOI: 10.1074/jbc.ra118.004169] [Citation(s) in RCA: 84] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2018] [Revised: 07/05/2018] [Indexed: 11/06/2022] Open
Abstract
Cellular heme is thought to be distributed between a pool of sequestered heme that is tightly bound within hemeproteins and a labile heme pool required for signaling and transfer into proteins. A heme chaperone that can hold and allocate labile heme within cells has long been proposed but never been identified. Here, we show that the glycolytic protein glyceraldehyde-3-phosphate dehydrogenase (GAPDH) fulfills this role by acting as an essential repository and allocator of bioavailable heme to downstream protein targets. We identified a conserved histidine in GAPDH that is needed for its robust heme binding both in vitro and in mammalian cells. Substitution of this histidine, and the consequent decreases in GAPDH heme binding, antagonized heme delivery to both cytosolic and nuclear hemeprotein targets, including inducible nitric-oxide synthase (iNOS) in murine macrophages and the nuclear transcription factor Hap1 in yeast, even though this GAPDH variant caused cellular levels of labile heme to rise dramatically. We conclude that by virtue of its heme-binding property, GAPDH binds and chaperones labile heme to create a heme pool that is bioavailable to downstream proteins. Our finding solves a fundamental question in cell biology and provides a new foundation for exploring heme homeostasis in health and disease.
Collapse
Affiliation(s)
- Elizabeth A Sweeny
- From the Department of Inflammation and Immunity, Lerner Research Institute, The Cleveland Clinic, Cleveland, Ohio 44195
| | - Anuradha Bharara Singh
- From the Department of Inflammation and Immunity, Lerner Research Institute, The Cleveland Clinic, Cleveland, Ohio 44195
| | - Ritu Chakravarti
- From the Department of Inflammation and Immunity, Lerner Research Institute, The Cleveland Clinic, Cleveland, Ohio 44195
| | - Osiris Martinez-Guzman
- the School of Chemistry and Biochemistry and Parker Petit Institute for Bioengineering and Biosciences, Georgia Institute of Technology, Atlanta, Georgia 30332
| | - Arushi Saini
- the School of Chemistry and Biochemistry and Parker Petit Institute for Bioengineering and Biosciences, Georgia Institute of Technology, Atlanta, Georgia 30332
| | - Mohammad Mahfuzul Haque
- From the Department of Inflammation and Immunity, Lerner Research Institute, The Cleveland Clinic, Cleveland, Ohio 44195
| | - Greer Garee
- From the Department of Inflammation and Immunity, Lerner Research Institute, The Cleveland Clinic, Cleveland, Ohio 44195
| | - Pablo D Dans
- the Institute for Research in Biomedicine (IRB Barcelona), Barcelona Institute of Science and Technology, Balidiri Reixac 10-12, Barcelona 08028, Spain, and
| | - Luciana Hannibal
- the Laboratory of Clinical Biochemistry and Metabolism, Center of Pediatrics, Medical Center, University of Freiburg, D-79106 Freiburg, Germany
| | - Amit R Reddi
- the School of Chemistry and Biochemistry and Parker Petit Institute for Bioengineering and Biosciences, Georgia Institute of Technology, Atlanta, Georgia 30332
| | - Dennis J Stuehr
- From the Department of Inflammation and Immunity, Lerner Research Institute, The Cleveland Clinic, Cleveland, Ohio 44195,
| |
Collapse
|
48
|
Hydration of counterions interacting with DNA double helix: a molecular dynamics study. J Mol Model 2018; 24:171. [DOI: 10.1007/s00894-018-3704-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2017] [Accepted: 06/06/2018] [Indexed: 12/12/2022]
|
49
|
Rodríguez-Decuadro S, Barraco-Vega M, Dans PD, Pandolfi V, Benko-Iseppon AM, Cecchetto G. Antimicrobial and structural insights of a new snakin-like peptide isolated from Peltophorum dubium (Fabaceae). Amino Acids 2018; 50:1245-1259. [PMID: 29948342 DOI: 10.1007/s00726-018-2598-3] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2018] [Accepted: 05/31/2018] [Indexed: 02/02/2023]
Abstract
Snakins are antimicrobial peptides (AMPs) found, so far, exclusively in plants, and known to be important in the defense against a wide range of pathogens. Like other plant AMPs, they contain several positively charged amino acids, and an even number of cysteine residues forming disulfide bridges which are considered important for their usual function. Despite its importance, studies on snakin tertiary structure and mode of action are still scarce. In this study, a new snakin-like gene was isolated from the native plant Peltophorum dubium, and its expression was verified in seedlings and adult leaves. The deduced peptide (PdSN1) shows 84% sequence identity with potato snakin-1 mature peptide, with the 12 cysteines characteristic from this peptide family at the GASA domain. The mature PdSN1 coding sequence was successfully expressed in Escherichia coli. The purified recombinant peptide inhibits the growth of important plant and human pathogens, like the economically relevant potato pathogen Streptomyces scabies and the opportunistic fungi Candida albicans and Aspergillus niger. Finally, homology and ab initio modeling techniques coupled to extensive molecular dynamics simulations were used to gain insight on the 3D structure of PdSN1, which exhibited a helix-turn-helix motif conserved in both native and recombinant peptides. We found this motif to be strongly coded in the sequence of PdSN1, as it is stable under different patterns of disulfide bonds connectivity, and even when the 12 cysteines are considered in their reduced form, explaining the previous experimental evidences.
Collapse
Affiliation(s)
- Susana Rodríguez-Decuadro
- Departamento de Biología Vegetal, Facultad de Agronomía, Universidad de la República, Garzón 780, 12900, Montevideo, Uruguay
| | - Mariana Barraco-Vega
- Departamento de Biociencias, Facultad de Química, Universidad de la República, General Flores 2124, 11800, Montevideo, Uruguay
| | - Pablo D Dans
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, 08028, Barcelona, Spain.,Joint BSC-IRB Research Program in Computational Biology, Baldiri Reixac 10-12, 08028, Barcelona, Spain
| | - Valesca Pandolfi
- Universidade Federal de Pernambuco, Centro de Biociências, Av. Prof. Moraes Rego, 1235, Recife, PE, CEP 50.670-420, Brazil
| | - Ana Maria Benko-Iseppon
- Universidade Federal de Pernambuco, Centro de Biociências, Av. Prof. Moraes Rego, 1235, Recife, PE, CEP 50.670-420, Brazil
| | - Gianna Cecchetto
- Departamento de Biociencias, Facultad de Química, Universidad de la República, General Flores 2124, 11800, Montevideo, Uruguay. .,Instituto de Química Biológica, Facultad de Ciencias, Facultad de Química, Universidad de la República, General Flores 2124, 11800, Montevideo, Uruguay.
| |
Collapse
|
50
|
Brovarets' OO, Tsiupa KS, Hovorun DM. Surprising Conformers of the Biologically Important A·T DNA Base Pairs: QM/QTAIM Proofs. Front Chem 2018; 6:8. [PMID: 29536003 PMCID: PMC5835050 DOI: 10.3389/fchem.2018.00008] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2017] [Accepted: 01/11/2018] [Indexed: 11/22/2022] Open
Abstract
For the first time novel high-energy conformers-A·T(wWC) (5.36), A·T(wrWC) (5.97), A·T(wH) (5.78), and A·T(wrH) (ΔG = 5.82 kcal·mol-1) (See Graphical Abstract) were revealed for each of the four biologically important A·T DNA base pairs - Watson-Crick A·T(WC), reverse Watson-Crick A·T(rWC), Hoogsteen A·T(H) and reverse Hoogsteen A·T(rH) at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of quantum-mechanical theory in the continuum with ε = 4 under normal conditions. Each of these conformers possesses substantially non-planar wobble (w) structure and is stabilized by the participation of the two anti-parallel N6H/N6H'…O4/O2 and N3H…N6 H-bonds, involving the pyramidalized amino group of the A DNA base as an acceptor and a donor of the H-bonding. The transition states - TSA·T(WC)↔A·T(wWC), TSA·T(rWC)↔A·T(wrWC), TSA·T(H)↔A·T(wH), and TSA·T(rH)↔A·T(wrH), controlling the dipole-active transformations of the conformers from the main plane-symmetric state into the high-energy, significantly non-planar state and vice versa, were localized. They also possess wobble structures similarly to the high-energy conformers and are stabilized by the participation of the N6H/N6H'…O4/O2 and N3H…N6 H-bonds. Discovered conformers of the A·T DNA base pairs are dynamically stable short-lived structures [lifetime τ = (1.4-3.9) ps]. Their possible biological significance and future perspectives have been briefly discussed.
Collapse
Affiliation(s)
- Ol'ha O. Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| | - Kostiantyn S. Tsiupa
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| |
Collapse
|