1
|
Grille L, Gallego D, Darré L, da Rosa G, Battistini F, Orozco M, Dans PD. The pseudotorsional space of RNA. RNA (NEW YORK, N.Y.) 2023; 29:1896-1909. [PMID: 37793790 PMCID: PMC10653382 DOI: 10.1261/rna.079821.123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 09/01/2023] [Indexed: 10/06/2023]
Abstract
The characterization of the conformational landscape of the RNA backbone is rather complex due to the ability of RNA to assume a large variety of conformations. These backbone conformations can be depicted by pseudotorsional angles linking RNA backbone atoms, from which Ramachandran-like plots can be built. We explore here different definitions of these pseudotorsional angles, finding that the most accurate ones are the traditional η (eta) and θ (theta) angles, which represent the relative position of RNA backbone atoms P and C4'. We explore the distribution of η - θ in known experimental structures, comparing the pseudotorsional space generated with structures determined exclusively by one experimental technique. We found that the complete picture only appears when combining data from different sources. The maps provide a quite comprehensive representation of the RNA accessible space, which can be used in RNA-structural predictions. Finally, our results highlight that protein interactions lead to significant changes in the population of the η - θ space, pointing toward the role of induced-fit mechanisms in protein-RNA recognition.
Collapse
Affiliation(s)
- Leandro Grille
- Computational Biophysics Group, Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República, 50000 Salto, Uruguay
- Bioinformatics Unit, Institute Pasteur of Montevideo, 11400 Montevideo, Uruguay
| | - Diego Gallego
- Molecular Modelling and Bioinformatics Group, Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Department of Biochemistry and Molecular Biology, Faculty of Biology, University of Barcelona, 08028 Barcelona, Spain
| | - Leonardo Darré
- Bioinformatics Unit, Institute Pasteur of Montevideo, 11400 Montevideo, Uruguay
- Molecular Modelling and Bioinformatics Group, Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Gabriela da Rosa
- Computational Biophysics Group, Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República, 50000 Salto, Uruguay
- Bioinformatics Unit, Institute Pasteur of Montevideo, 11400 Montevideo, Uruguay
| | - Federica Battistini
- Molecular Modelling and Bioinformatics Group, Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Department of Biochemistry and Molecular Biology, Faculty of Biology, University of Barcelona, 08028 Barcelona, Spain
| | - Modesto Orozco
- Molecular Modelling and Bioinformatics Group, Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Department of Biochemistry and Molecular Biology, Faculty of Biology, University of Barcelona, 08028 Barcelona, Spain
| | - Pablo D Dans
- Computational Biophysics Group, Department of Biological Sciences, CENUR Litoral Norte, Universidad de la República, 50000 Salto, Uruguay
- Bioinformatics Unit, Institute Pasteur of Montevideo, 11400 Montevideo, Uruguay
- Molecular Modelling and Bioinformatics Group, Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| |
Collapse
|
2
|
Kolesnikov ES, Gushchin IY, Zhilyaev PA, Onufriev AV. Why Na+ has higher propensity than K+ to condense DNA in a crowded environment. J Chem Phys 2023; 159:145103. [PMID: 37815107 DOI: 10.1063/5.0159341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 08/22/2023] [Indexed: 10/11/2023] Open
Abstract
Experimentally, in the presence of the crowding agent polyethylene glycol (PEG), sodium ions compact double-stranded DNA more readily than potassium ions. Here, we have used molecular dynamics simulations and the "ion binding shells model" of DNA condensation to provide an explanation for the observed variations in condensation of short DNA duplexes in solutions containing different monovalent cations and PEG; several predictions are made. According to the model we use, externally bound ions contribute the most to the ion-induced aggregation of DNA duplexes. The simulations reveal that for two adjacent DNA duplexes, the number of externally bound Na+ ions is larger than the number of K+ ions over a wide range of chloride concentrations in the presence of PEG, providing a qualitative explanation for the higher propensity of sodium ions to compact DNA under crowded conditions. The qualitative picture is confirmed by an estimate of the corresponding free energy of DNA aggregation that is at least 0.2kBT per base pair more favorable in solution with NaCl than with KCl at the same ion concentration. The estimated attraction free energy of DNA duplexes in the presence of Na+ depends noticeably on the DNA sequence; we predict that AT-rich DNA duplexes are more readily condensed than GC-rich ones in the presence of Na+. Counter-intuitively, the addition of a small amount of a crowding agent with high affinity for the specific condensing ion may lead to the weakening of the ion-mediated DNA-DNA attraction, shifting the equilibrium away from the DNA condensed phase.
Collapse
Affiliation(s)
- Egor S Kolesnikov
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, 141700 Dolgoprudny, Russia
- Department of Physics, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Ivan Yu Gushchin
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, 141700 Dolgoprudny, Russia
| | - Petr A Zhilyaev
- The Center for Materials Technologies, Skolkovo Institute of Science and Technology, Bolshoy Boulevard 30, bld. 1, Moscow 121205, Russia
| | - Alexey V Onufriev
- Department of Physics, Virginia Tech, Blacksburg, Virginia 24061, USA
- Department of Computer Science, Virginia Tech, 2160C Torgersen Hall, Blacksburg, Virginia 24061, USA
- Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
3
|
Mu ZC, Tan YL, Liu J, Zhang BG, Shi YZ. Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding. Molecules 2023; 28:4833. [PMID: 37375388 DOI: 10.3390/molecules28124833] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 06/11/2023] [Accepted: 06/14/2023] [Indexed: 06/29/2023] Open
Abstract
DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.
Collapse
Affiliation(s)
- Zi-Chun Mu
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan 430073, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Jie Liu
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Ben-Gong Zhang
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan 430073, China
| |
Collapse
|
4
|
Muniz MI, Bustos AH, Slott S, Astakhova K, Weber G. Cation valence dependence of hydrogen bond and stacking potentials in DNA mesoscopic models. Biophys Chem 2023; 294:106949. [PMID: 36706510 DOI: 10.1016/j.bpc.2022.106949] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Revised: 12/23/2022] [Accepted: 12/25/2022] [Indexed: 12/30/2022]
Abstract
Monovalent and divalent cations play a crucial role in living cells and for molecular techniques such as PCR. Here we evaluate DNA melting temperatures in magnesium (Mg2+) and magnesium‑potassium (Mg2++ K+) buffers with a mesoscopic model that allows us to estimate hydrogen bonds and stacking interaction potentials. The Mg2+ and Mg2++ K+ results are compared to previous calculations for sodium ions (Na+), in terms of equivalent sodium concentration and ionic strength. Morse potentials, related to hydrogen bonding, were found to be essentially constant and unaffected by cation conditions. However, for stacking interactions we find a clear dependence with ionic strength and cation valence. The highest ionic strength variations, for both hydrogen bonds and stacking interactions, was found at the sequence terminals. This suggests that end-to-end interactions in DNA will be strongly dependent on cation valence and ionic strength.
Collapse
Affiliation(s)
- Maria Izabel Muniz
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil
| | - Adrian H Bustos
- Department of Chemistry, Technical University of Denmark, 206-207 Kemitorvet, 2800 Kongens Lyngby, Denmark
| | - Sofie Slott
- Department of Chemistry, Technical University of Denmark, 206-207 Kemitorvet, 2800 Kongens Lyngby, Denmark
| | - Kira Astakhova
- Department of Chemistry, Technical University of Denmark, 206-207 Kemitorvet, 2800 Kongens Lyngby, Denmark
| | - Gerald Weber
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil.
| |
Collapse
|
5
|
Pant P. Harmonizing Interstrand Electrostatic Repulsion by Conformational Rigidity in Counterion-Deprived Z-DNA: A Molecular Dynamics Study. J Phys Chem B 2022; 126:9956-9963. [PMID: 36412276 DOI: 10.1021/acs.jpcb.2c04527] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Deoxyribonucleic acid (DNA) is a vital biomacromolecule. Although the right-handed B-DNA type helical structure is the most abundant and extensively studied form of DNA, several noncanonical forms, such as triplex, quadruplex, Z-DNA, A-DNA, and ss-DNA, have been probed from time to time to gain insights into the DNA's function. Z-DNA was recently found to be involved in cancer and several autoimmune diseases. In the present Article, we evaluated the conformational stability of locked-sugar-based Z-DNA via all-atom explicit-solvent molecular dynamics simulations and found that the modified DNA maintained the left-handed conformation even in the absence of counterions, wherein the structural rigidity dominates over the electrostatic repulsion between the complementary strands. The control Z-DNA without counterions, as expected, instantaneously resulted in unfolded states. The remarkable stability of the conformationally locked model system was thoroughly investigated via structural and energetic perspectives and was probably the result of the backbone widening in tandem with enhanced electrostatics between complementary strands. We believe that the design of the proposed modified Z-DNA construct could help understand the otherwise delicate Z-DNA conformation even in salt-deprived conditions. The design could also motivate the medicinal use of short segments of such modified nucleotides and could be utilized in more advanced modeling techniques, such as DNA origami which has gained popularity in recent years.
Collapse
Affiliation(s)
- Pradeep Pant
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi 110016, India
| |
Collapse
|
6
|
Patra P, Gao YQ. Sequence-Specific Structural Features and Solvation Properties of Transcription Factor Binding DNA Motifs: Insights from Molecular Dynamics Simulation. J Phys Chem B 2022; 126:9187-9206. [PMID: 36322688 DOI: 10.1021/acs.jpcb.2c05749] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Sequence-specific recognition of transcription factor (TF) binding motifs in the target site of DNA over the vast amount of non-target DNA is of primary importance for the transcriptional regulation of gene expression by the TFs. Binding of TFs to the target site of DNA relies not only on the direct contact formation but also on the structural and conformational features of DNA. Recognition of DNA structural features or shape readout by proteins is an important factor in the context of TF-DNA interaction. Based on the atomistic molecular simulation, here we report the sequence-dependent unique structural features, solvation, and ion-binding properties of biologically relevant AT- and GC-rich human TF binding motifs in DNA. Counterion and water distribution around the motif is found to be sensitive to the motif sequence, which is accompanied with the DNA shape features. The motif sequence affects the electrostatic potential along the grooves, and cytosine methylation alters the DNA shape features. Characteristic solvation properties of TF binding motif DNA fragments infer that an ionic environment and hydration influences are essential to describe TF-DNA interactions.
Collapse
Affiliation(s)
- Piya Patra
- Shenzhen Bay Laboratory, Institute of Systems and Physical Biology, Shenzhen 518107, China
| | - Yi Qin Gao
- Shenzhen Bay Laboratory, Institute of Systems and Physical Biology, Shenzhen 518107, China.,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China.,Biomedical Pioneering Innovation Center, Peking University, Beijing 100871, China
| |
Collapse
|
7
|
Cruz-León S, Schwierz N. RNA Captures More Cations than DNA: Insights from Molecular Dynamics Simulations. J Phys Chem B 2022; 126:8646-8654. [PMID: 36260822 PMCID: PMC9639116 DOI: 10.1021/acs.jpcb.2c04488] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
The distribution of cations around nucleic acids is essential for a broad variety of processes ranging from DNA condensation and RNA folding to the detection of biomolecules in biosensors. Predicting the exact distribution of ions remains challenging since the distribution and, hence, a broad variety of nucleic acid properties depend on the salt concentration, the valency of the ions, and the ion type. Despite the importance, a general theory to quantify ion-specific effects for highly charged biomolecules is still lacking. Moreover, recent experiments reveal that despite their similar building blocks, DNA and RNA duplexes can react differently to the same ionic conditions. The aim of our current work is to provide a comprehensive set of molecular dynamics simulations using more than 180 μs of simulation time. For the mono- and divalent cations Li+, Na+, K+, Cs+, Ca2+, Sr2+, and Ba2+, the simulations allow us to reveal the ion-specific distributions and binding patterns for DNA and RNA duplexes. The microscopic insights from the simulations display the origin of ion-specificity and shed light on the question of why DNA and RNA show opposing behavior in the same ionic conditions. Finally, the detailed binding patterns from the simulations reveal why RNA can capture more cations than DNA.
Collapse
Affiliation(s)
- Sergio Cruz-León
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Str. 3, 60438Frankfurt am Main, Germany
| | - Nadine Schwierz
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Str. 3, 60438Frankfurt am Main, Germany,Institute
of Physics, University of Augsburg, Universitätsstraße 1, 86159Augsburg, Germany,E-mail:
| |
Collapse
|
8
|
Singh A, Maity A, Singh N. Structure and Dynamics of dsDNA in Cell-like Environments. ENTROPY (BASEL, SWITZERLAND) 2022; 24:1587. [PMID: 36359677 PMCID: PMC9689892 DOI: 10.3390/e24111587] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/25/2022] [Revised: 10/28/2022] [Accepted: 10/31/2022] [Indexed: 06/01/2023]
Abstract
Deoxyribonucleic acid (DNA) is a fundamental biomolecule for correct cellular functioning and regulation of biological processes. DNA's structure is dynamic and has the ability to adopt a variety of structural conformations in addition to its most widely known double-stranded DNA (dsDNA) helix structure. Stability and structural dynamics of dsDNA play an important role in molecular biology. In vivo, DNA molecules are folded in a tightly confined space, such as a cell chamber or a channel, and are highly dense in solution; their conformational properties are restricted, which affects their thermodynamics and mechanical properties. There are also many technical medical purposes for which DNA is placed in a confined space, such as gene therapy, DNA encapsulation, DNA mapping, etc. Physiological conditions and the nature of confined spaces have a significant influence on the opening or denaturation of DNA base pairs. In this review, we summarize the progress of research on the stability and dynamics of dsDNA in cell-like environments and discuss current challenges and future directions. We include studies on various thermal and mechanical properties of dsDNA in ionic solutions, molecular crowded environments, and confined spaces. By providing a better understanding of melting and unzipping of dsDNA in different environments, this review provides valuable guidelines for predicting DNA thermodynamic quantities and for designing DNA/RNA nanostructures.
Collapse
|
9
|
Abdelgawwad AMA, Monari A, Tuñón I, Francés-Monerris A. Spatial and Temporal Resolution of the Oxygen-Independent Photoinduced DNA Interstrand Cross-Linking by a Nitroimidazole Derivative. J Chem Inf Model 2022; 62:3239-3252. [PMID: 35771238 PMCID: PMC9277591 DOI: 10.1021/acs.jcim.2c00460] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
DNA damage is ubiquitous in nature and is at the basis of emergent treatments such as photodynamic therapy, which is based on the activation of highly oxidative reactive oxygen species by photosensitizing O2. However, hypoxia observed in solid tumors imposes the necessity to devise oxygen-independent modes of action able to induce DNA damage under a low oxygen concentration. The complexity of these DNA damage mechanisms in realistic environments grows exponentially when taking into account light absorption and subsequent excited-state population, photochemical and (photo)-redox reactions, the multiple species involved in different electronic states, noncovalent interactions, multiple reaction steps, and the large number of DNA reactive sites. This work tackles all the intricate reactivity of a photosensitizer based on a nitroimidazole derivative reacting toward DNA in solution under UV light exposition. This is performed through a combination of ground- and excited-state quantum chemistry, classical molecular dynamics, and hybrid QM/MM simulations to rationalize in detail the formation of DNA interstrand cross-links (ICLs) exerted by the noncanonical noncovalent photosensitizer. Unprecedented spatial and temporal resolution of these phenomena is achieved, revealing that the ICL is sequence-specific and that the fastest reactions take place at AT, GC, and GT steps involving either the opposite nucleobases or adjacent Watson-Crick base pairs. The N7 and O6 positions of guanine, the N7 and N3 sites of adenine, the N4 position of cytosine, and the O2 atom of thymine are deemed as the most nucleophile sites and are positively identified to participate in the ICL productions. This work provides a multiscale computational protocol to study DNA reactivity with noncovalent photosensitizers, and contributes to the understanding of therapies based on photoinduced DNA damage at molecular and electronic levels. In addition, we believe the depth understanding of these processes should assist the design of new photosensitizers considering their molecular size, electronic properties, and the observed regioselectivity toward nucleic acids.
Collapse
Affiliation(s)
| | - Antonio Monari
- Université Paris Cité, CNRS, ITODYS, F-75006 Paris, France.,Université de Lorraine and CNRS, UMR 7019 LPCT, F-5400 Nancy, France
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | | |
Collapse
|
10
|
Cruz-León S, Vanderlinden W, Müller P, Forster T, Staudt G, Lin YY, Lipfert J, Schwierz N. Twisting DNA by salt. Nucleic Acids Res 2022; 50:5726-5738. [PMID: 35640616 PMCID: PMC9177979 DOI: 10.1093/nar/gkac445] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2021] [Revised: 05/06/2022] [Accepted: 05/10/2022] [Indexed: 12/11/2022] Open
Abstract
The structure and properties of DNA depend on the environment, in particular the ion atmosphere. Here, we investigate how DNA twist -one of the central properties of DNA- changes with concentration and identity of the surrounding ions. To resolve how cations influence the twist, we combine single-molecule magnetic tweezer experiments and extensive all-atom molecular dynamics simulations. Two interconnected trends are observed for monovalent alkali and divalent alkaline earth cations. First, DNA twist increases monotonously with increasing concentration for all ions investigated. Second, for a given salt concentration, DNA twist strongly depends on cation identity. At 100 mM concentration, DNA twist increases as Na+ < K+ < Rb+ < Ba2+ < Li+ ≈ Cs+ < Sr2+ < Mg2+ < Ca2+. Our molecular dynamics simulations reveal that preferential binding of the cations to the DNA backbone or the nucleobases has opposing effects on DNA twist and provides the microscopic explanation of the observed ion specificity. However, the simulations also reveal shortcomings of existing force field parameters for Cs+ and Sr2+. The comprehensive view gained from our combined approach provides a foundation for understanding and predicting cation-induced structural changes both in nature and in DNA nanotechnology.
Collapse
Affiliation(s)
- Sergio Cruz-León
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| | - Willem Vanderlinden
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Peter Müller
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Tobias Forster
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Georgina Staudt
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Yi-Yun Lin
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Jan Lipfert
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| |
Collapse
|
11
|
Zhang C, Tian F, Lu Y, Yuan B, Tan ZJ, Zhang XH, Dai L. Twist-diameter coupling drives DNA twist changes with salt and temperature. SCIENCE ADVANCES 2022; 8:eabn1384. [PMID: 35319990 PMCID: PMC8942373 DOI: 10.1126/sciadv.abn1384] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/05/2021] [Accepted: 01/31/2022] [Indexed: 06/14/2023]
Abstract
DNA deformations upon environmental changes, e.g., salt and temperature, play crucial roles in many biological processes and material applications. Here, our magnetic tweezers experiments observed that the increase in NaCl, KCl, or RbCl concentration leads to substantial DNA overwinding. Our simulations and theoretical calculation quantitatively explain the salt-induced twist change through the mechanism: More salt enhances the screening of interstrand electrostatic repulsion and hence reduces DNA diameter, which is transduced to twist increase through twist-diameter coupling. We determined that the coupling constant is 4.5 ± 0.8 kBT/(degrees∙nm) for one base pair. The coupling comes from the restraint of the contour length of DNA backbone. On the basis of this coupling constant and diameter-dependent DNA conformational entropy, we predict the temperature dependence of DNA twist Δωbp/ΔT ≈ -0.01 degree/°C, which agrees with our and previous experimental results. Our analysis suggests that twist-diameter coupling is a common driving force for salt- and temperature-induced DNA twist changes.
Collapse
Affiliation(s)
- Chen Zhang
- College of Life Sciences, The Institute for Advanced Studies, State Key Laboratory of Virology, Hubei Key Laboratory of Cell Homeostasis, Wuhan University, Wuhan 430072, China
| | - Fujia Tian
- Department of Physics, City University of Hong Kong, Hong Kong 999077, China
| | - Ying Lu
- Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Bing Yuan
- Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, China
| | - Zhi-Jie Tan
- School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Xing-Hua Zhang
- College of Life Sciences, The Institute for Advanced Studies, State Key Laboratory of Virology, Hubei Key Laboratory of Cell Homeostasis, Wuhan University, Wuhan 430072, China
| | - Liang Dai
- Department of Physics, City University of Hong Kong, Hong Kong 999077, China
- Shenzhen Research Institute, City University of Hong Kong, Shenzhen 518057, China
| |
Collapse
|
12
|
Pant P, Aggarwal L. Assessing the DNA structural integrity via selective annihilation of Watson-Crick hydrogen bonds: Insights from molecular dynamics simulations. Biophys Chem 2022; 282:106758. [DOI: 10.1016/j.bpc.2021.106758] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Revised: 12/09/2021] [Accepted: 12/30/2021] [Indexed: 01/17/2023]
|
13
|
Vasiliu T, Mocci F, Laaksonen A, Engelbrecht LDV, Perepelytsya S. Caging Polycations: Effect of Increasing Confinement on the Modes of Interaction of Spermidine3+ With DNA Double Helices. Front Chem 2022; 10:836994. [PMID: 35281557 PMCID: PMC8915389 DOI: 10.3389/fchem.2022.836994] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Accepted: 02/04/2022] [Indexed: 11/23/2022] Open
Abstract
Polyamines have important roles in the modulation of the cellular function and are ubiquitous in cells. The polyamines putrescine2+, spermidine3+, and spermine4+ represent the most abundant organic counterions of the negatively charged DNA in the cellular nucleus. These polyamines are known to stabilize the DNA structure and, depending on their concentration and additional salt composition, to induce DNA aggregation, which is often referred to as condensation. However, the modes of interactions of these elongated polycations with DNA and how they promote condensation are still not clear. In the present work, atomistic molecular dynamics (MD) computer simulations of two DNA fragments surrounded by spermidine3+ (Spd3+) cations were performed to study the structuring of Spd3+ “caged” between DNA molecules. Microsecond time scale simulations, in which the parallel DNA fragments were constrained at three different separations, but allowed to rotate axially and move naturally, provided information on the conformations and relative orientations of surrounding Spm3+ cations as a function of DNA-DNA separation. Novel geometric criteria allowed for the classification of DNA-Spd3+ interaction modes, with special attention given to Spd3+ conformational changes in the space between the two DNA molecules (caged Spd3+). This work shows how changes in the accessible space, or confinement, around DNA affect DNA-Spd3+ interactions, information fundamental to understanding the interactions between DNA and its counterions in environments where DNA is compacted, e.g. in the cellular nucleus.
Collapse
Affiliation(s)
- Tudor Vasiliu
- Centre of Advanced Research in Bionanoconjugates and Biopolymers “Petru Poni” Institute of Macromolecular Chemistry, Iasi, Romania
| | - Francesca Mocci
- Dipartimento di Scienze Chimiche e Geologiche, Cagliari University, Cagliari, Italy
- *Correspondence: Francesca Mocci, ; Aatto Laaksonen, ; Sergiy Perepelytsya,
| | - Aatto Laaksonen
- Centre of Advanced Research in Bionanoconjugates and Biopolymers “Petru Poni” Institute of Macromolecular Chemistry, Iasi, Romania
- Dipartimento di Scienze Chimiche e Geologiche, Cagliari University, Cagliari, Italy
- Division of Energy Science, Energy Engineering, Luleå University of Technology, Luleå, Sweden
- Division of Physical Chemistry, Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, Sweden
- State Key Laboratory of Materials-Oriented and Chemical Engineering, Nanjing Tech University, Nanjing, China
- *Correspondence: Francesca Mocci, ; Aatto Laaksonen, ; Sergiy Perepelytsya,
| | | | - Sergiy Perepelytsya
- Bogolyubov Institute for Theoretical Physics of the NAS of Ukraine, Kyiv, Ukraine
- *Correspondence: Francesca Mocci, ; Aatto Laaksonen, ; Sergiy Perepelytsya,
| |
Collapse
|
14
|
Kristoffersen EL, Burman M, Noy A, Holliger P. Rolling circle RNA synthesis catalysed by RNA. eLife 2022; 11:75186. [PMID: 35108196 PMCID: PMC8937235 DOI: 10.7554/elife.75186] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2021] [Accepted: 02/01/2022] [Indexed: 11/13/2022] Open
Abstract
RNA-catalyzed RNA replication is widely considered a key step in the emergence of life’s first genetic system. However, RNA replication can be impeded by the extraordinary stability of duplex RNA products, which must be dissociated for re-initiation of the next replication cycle. Here, we have explored rolling circle synthesis (RCS) as a potential solution to this strand separation problem. We observe sustained RCS by a triplet polymerase ribozyme beyond full-length circle synthesis with strand displacement yielding concatemeric RNA products. Furthermore, we show RCS of a circular Hammerhead ribozyme capable of self-cleavage and re-circularization. Thus, all steps of a viroid-like RNA replication pathway can be catalyzed by RNA alone. Finally, we explore potential RCS mechanisms by molecular dynamics simulations, which indicate a progressive build-up of conformational strain upon RCS with destabilization of nascent strand 5′- and 3′-ends. Our results have implications for the emergence of RNA replication and for understanding the potential of RNA to support complex genetic processes. Many organisms today rely on a trio of molecules for their survival: DNA, to store their genetic information; proteins, to conduct the biological processes required for growth or replication; and RNA, to mainly act as an intermediary between DNA and proteins. Yet, how these inanimate molecules first came together to form a living system remains unclear. Circumstantial evidence suggests that the first lifeforms relied to a much greater exrtent on RNA to conduct all necessary biological processes. There is no trace of this ‘RNA world’ today, but molecular ‘fossils’ may exist in current biology. Viroids, for example, are agents which can infect and replicate inside plant cells. They are formed of nothing but a circular strand of RNA that serves not only as genetic storage but also as ribozymes (RNA-based enzymes). Viroids need proteins from the host plant to replicate, but scientists have been able to engineer ribozymes that can copy complex RNA strands. This suggests that viroid-like replication could be achieved using only RNA. Kristoffersen et al. put this idea to the test and showed that it is possible to use RNA enzymatic activity alone to carry out all the steps of a viroid-like copying mechanism. This process included copying a viroid-like RNA circle with RNA, followed by trimming the copy to the right size and reforming the circle. These two latter steps could be carried out by a ribozyme that could itself be encoded on the RNA circle. A computer simulation indicated that RNA synthesis on the circle caused increasing tension that could ease some of the barriers to replication. These results increase our understanding of how RNA copying by RNA could be possible. This may lead to developing molecular models of a primordial RNA-based replication, which could be used to investigate early genetic systems and may have potential applications in synthetic biology.
Collapse
Affiliation(s)
| | - Matthew Burman
- Department of Physics, University of York, York, United Kingdom
| | - Agnes Noy
- Department of Physics, University of York, York, United Kingdom
| | | |
Collapse
|
15
|
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
|
16
|
Fingerhut BP, Schauss J, Kundu A, Elsaesser T. Contact pairs of RNA with magnesium ions-electrostatics beyond the Poisson-Boltzmann equation. Biophys J 2021; 120:5322-5332. [PMID: 34715079 PMCID: PMC8715182 DOI: 10.1016/j.bpj.2021.10.029] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2020] [Revised: 08/24/2021] [Accepted: 10/22/2021] [Indexed: 11/18/2022] Open
Abstract
The electrostatic interaction of RNA with its aqueous environment is most relevant for defining macromolecular structure and biological function. The attractive interaction of phosphate groups in the RNA backbone with ions in the water environment leads to the accumulation of positively charged ions in the first few hydration layers around RNA. Electrostatics of this ion atmosphere and the resulting ion concentration profiles have been described by solutions of the nonlinear Poisson-Boltzmann equation and atomistic molecular dynamics (MD) simulations. Much less is known on contact pairs of RNA phosphate groups with ions at the RNA surface, regarding their abundance, molecular geometry, and role in defining RNA structure. Here, we present a combined theoretical and experimental study of interactions of a short RNA duplex with magnesium (Mg2+) ions. MD simulations covering a microsecond time range give detailed hydration geometries as well as electrostatics and spatial arrangements of phosphate-Mg2+ pairs, including both pairs in direct contact and separated by a single water layer. The theoretical predictions are benchmarked by linear infrared absorption and nonlinear two-dimensional infrared spectra of the asymmetric phosphate stretch vibration which probes both local interaction geometries and electric fields. Contact pairs of phosphate groups and Mg2+ ions are identified via their impact on the vibrational frequency position and line shape. A quantitative analysis of infrared spectra for a range of Mg2+-excess concentrations and comparison with fluorescence titration measurements shows that on average 20-30% of the Mg2+ ions interacting with the RNA duplex form contact pairs. The experimental and MD results are in good agreement. In contrast, calculations based on the nonlinear Poisson-Boltzmann equation fail in describing the ion arrangement, molecular electrostatic potential, and local electric field strengths correctly. Our results underline the importance of local electric field mapping and molecular-level simulations to correctly account for the electrostatics at the RNA-water interface.
Collapse
|
17
|
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
|
18
|
Kolesnikov ES, Gushchin IY, Zhilyaev PA, Onufriev AV. Similarities and Differences between Na + and K + Distributions around DNA Obtained with Three Popular Water Models. J Chem Theory Comput 2021; 17:7246-7259. [PMID: 34633813 DOI: 10.1021/acs.jctc.1c00332] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We have compared distributions of sodium and potassium ions around double-stranded DNA, simulated using fixed charge SPC/E, TIP3P, and OPC water models and the Joung/Cheatham (J/C) ion parameter set, as well as the Li/Merz HFE 6-12 (L/M HFE) ion parameters for OPC water. In all the simulations, the ion distributions are in qualitative agreement with Manning's condensation theory and the Debye-Hückel theory, where expected. In agreement with experiment, binding affinity of monovalent ions to DNA does not depend on ion type in every solvent model. However, behavior of deeply bound ions, including ions bound to specific sites, depends strongly on the solvent model. In particular, the number of potassium ions in the minor groove of AT-tracts differs at least 3-fold between the solvent models tested. The number of sodium ions associated with the DNA agrees quantitatively with the experiment for the OPC water model, followed closely by TIP3P+J/C; the largest deviation from the experiment, ∼10%, is seen for SPC/E+J/C. On the other hand, SPC/E+J/C model is most consistent (67%) with the experimental potassium binding sites, followed by OPC+J/C (60%), TIP3P+J/C (53%), and OPC+L/M HFE (27%). The use of NBFIX correction with TIP3P+J/C improves its consistency with the experiment. In summary, the choice of the solvent model matters little for simulating the diffuse atmosphere of sodium and potassium ions around DNA, but ion distributions become increasingly sensitive to the solvent model near the helical axis. We offer an explanation for these trends. There is no single gold standard solvent model, although OPC water with J/C ions or TIP3P with J/C + NBFIX may offer an imperfect compromise for practical simulations of ionic atmospheres around DNA.
Collapse
Affiliation(s)
- Egor S Kolesnikov
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, Dolgoprudny 141700, Russia
| | - Ivan Yu Gushchin
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, Dolgoprudny 141700, Russia
| | - Petr A Zhilyaev
- Center for Design, Manufacturing and Materials, Skolkovo Institute of Science and Technology, Bolshoy Boulevard 30, bld. 1, Moscow 121205, Russia
| | - Alexey V Onufriev
- Department of Computer Science, Virginia Tech, Blacksburg 24061-0131, United States.,Department of Physics, Virginia Tech, Blacksburg 24061-0131, United States.,Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg 24061-0131, United States
| |
Collapse
|
19
|
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
|
20
|
Dohnalová H, Lankaš F. Deciphering the mechanical properties of
B‐DNA
duplex. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1575] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Hana Dohnalová
- Department of Informatics and Chemistry University of Chemistry and Technology Prague Praha 6 Czech Republic
| | - Filip Lankaš
- Department of Informatics and Chemistry University of Chemistry and Technology Prague Praha 6 Czech Republic
| |
Collapse
|
21
|
Galindo-Murillo R, Cheatham TE. Ethidium bromide interactions with DNA: an exploration of a classic DNA-ligand complex with unbiased molecular dynamics simulations. Nucleic Acids Res 2021; 49:3735-3747. [PMID: 33764383 PMCID: PMC8053101 DOI: 10.1093/nar/gkab143] [Citation(s) in RCA: 42] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Revised: 02/03/2021] [Accepted: 03/16/2021] [Indexed: 01/27/2023] Open
Abstract
Visualization of double stranded DNA in gels with the binding of the fluorescent dye ethidium bromide has been a basic experimental technique in any molecular biology laboratory for >40 years. The interaction between ethidium and double stranded DNA has been observed to be an intercalation between base pairs with strong experimental evidence. This presents a unique opportunity for computational chemistry and biomolecular simulation techniques to benchmark and assess their models in order to see if the theory can reproduce experiments and ultimately provide new insights. We present molecular dynamics simulations of the interaction of ethidium with two different double stranded DNA models. The first model system is the classic sequence d(CGCGAATTCGCG)2 also known as the Drew–Dickerson dodecamer. We found that the ethidium ligand binds mainly stacked on, or intercalated between, the terminal base pairs of the DNA with little to no interaction with the inner base pairs. As the intercalation at the terminal CpG steps is relatively rapid, the resultant DNA unwinding, rigidification, and increased stability of the internal base pair steps inhibits further intercalation. In order to reduce these interactions and to provide a larger groove space, a second 18-mer DNA duplex system with the sequence d(GCATGAACGAACGAACGC) was tested. We computed molecular dynamics simulations for 20 independent replicas with this sequence, each with ∼27 μs of sampling time. Results show several spontaneous intercalation and base-pair eversion events that are consistent with experimental observations. The present work suggests that extended MD simulations with modern DNA force fields and optimized simulation codes are allowing the ability to reproduce unbiased intercalation events that we were not able to previously reach due to limits in computing power and the lack of extensively tested force fields and analysis tools.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 306, Salt Lake City, UT 84112, USA
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 306, Salt Lake City, UT 84112, USA
| |
Collapse
|
22
|
Elsaesser T, Schauss J, Kundu A, Fingerhut BP. Phosphate Vibrations Probe Electric Fields in Hydrated Biomolecules: Spectroscopy, Dynamics, and Interactions. J Phys Chem B 2021; 125:3899-3908. [PMID: 33834783 PMCID: PMC8154594 DOI: 10.1021/acs.jpcb.1c01502] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Electric interactions
have a strong impact on the structure and
dynamics of biomolecules in their native water environment. Given
the variety of water arrangements in hydration shells and the femto-
to subnanosecond time range of structural fluctuations, there is a
strong quest for sensitive noninvasive probes of local electric fields.
The stretching vibrations of phosphate groups, in particular the asymmetric
(PO2)− stretching vibration νAS(PO2)−, allow for a quantitative
mapping of dynamic electric fields in aqueous environments via a field-induced
redshift of their transition frequencies and concomitant changes of
vibrational line shapes. We present a systematic study of νAS(PO2)− excitations in molecular
systems of increasing complexity, including dimethyl phosphate (DMP),
short DNA and RNA duplex structures, and transfer RNA (tRNA) in water.
A combination of linear infrared absorption, two-dimensional infrared
(2D-IR) spectroscopy, and molecular dynamics (MD) simulations gives
quantitative insight in electric-field tuning rates of vibrational
frequencies, electric field and fluctuation amplitudes, and molecular
interaction geometries. Beyond neat water environments, the formation
of contact ion pairs of phosphate groups with Mg2+ ions
is demonstrated via frequency upshifts of the νAS(PO2)− vibration, resulting in a distinct
vibrational band. The frequency positions of contact geometries are
determined by an interplay of attractive electric and repulsive exchange
interactions.
Collapse
Affiliation(s)
- Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Jakob Schauss
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Achintya Kundu
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| |
Collapse
|
23
|
Schauss J, Kundu A, Fingerhut BP, Elsaesser T. Magnesium Contact Ions Stabilize the Tertiary Structure of Transfer RNA: Electrostatics Mapped by Two-Dimensional Infrared Spectra and Theoretical Simulations. J Phys Chem B 2021; 125:740-747. [PMID: 33284610 PMCID: PMC7848891 DOI: 10.1021/acs.jpcb.0c08966] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
![]()
Ions interacting with hydrated RNA
play a central role in defining
its secondary and tertiary structure. While spatial arrangements of
ions, water molecules, and phosphate groups have been inferred from
X-ray studies, the role of electrostatic and other noncovalent interactions
in stabilizing compact folded RNA structures is not fully understood
at the molecular level. Here, we demonstrate that contact ion pairs
of magnesium (Mg2+) and phosphate groups embedded in local
water shells stabilize the tertiary equilibrium structure of transfer
RNA (tRNA). Employing dialyzed tRNAPhe from yeast and tRNA
from Escherichia coli, we follow the
population of Mg2+ sites close to phosphate groups of the
ribose-phosphodiester backbone step by step, combining linear and
nonlinear infrared spectroscopy of phosphate vibrations with molecular
dynamics simulations and ab initio vibrational frequency calculations.
The formation of up to six Mg2+/phosphate contact pairs
per tRNA and local field-induced reorientations of water molecules
balance the phosphate–phosphate repulsion in nonhelical parts
of tRNA, thus stabilizing the folded structure electrostatically.
Such geometries display limited sub-picosecond fluctuations in the
arrangement of water molecules and ion residence times longer than
1 μs. At higher Mg2+ excess, the number of contact
ion pairs per tRNA saturates around 6 and weakly interacting ions
prevail. Our results suggest a predominance of contact ion pairs over
long-range coupling of the ion atmosphere and the biomolecule in defining
and stabilizing the tertiary structure of tRNA.
Collapse
Affiliation(s)
- Jakob Schauss
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Achintya Kundu
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| |
Collapse
|
24
|
Zdorevskyi OO, Perepelytsya SM. Dynamics of K + counterions around DNA double helix in the external electric field: A molecular dynamics study. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2020; 43:77. [PMID: 33306165 DOI: 10.1140/epje/i2020-12000-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2020] [Accepted: 11/05/2020] [Indexed: 06/12/2023]
Abstract
The structure of DNA double helix is stabilized by metal counterions condensed to a diffuse layer around the macromolecule. The dynamics of counterions in real conditions is governed by the electric fields from DNA and other biological macromolecules. In the present work the molecular dynamics study was performed for the system of DNA double helix with neutralizing K+ counterions and for the system of KCl salt solution in an external electric field of different strength (up to 32mV/Å). The analysis of ionic conductivities of these systems has shown that the counterions around the DNA double helix are slowed down compared with the KCl salt solution. The calculated values of ion mobility are within (0.05-0.4)mS/cm depending on the orientation of the external electric field relatively to the double helix. Under the electric field parallel to the macromolecule K+ counterions move along the grooves of the double helix staying longer in the places with narrower minor groove. Under the electric field perpendicular to the macromolecule the dynamics of counterions is less affected by DNA atoms, and starting with the electric field values about 30mV/Å the double helix undergoes a phase transition from a double-stranded to a single-strand state.
Collapse
Affiliation(s)
- O O Zdorevskyi
- Bogolyubov Institute for Theoretical Physics of the National Academy of Sciences of Ukraine, 14-b, Metrolohichna Str., 03143, Kiev, Ukraine.
| | - S M Perepelytsya
- Bogolyubov Institute for Theoretical Physics of the National Academy of Sciences of Ukraine, 14-b, Metrolohichna Str., 03143, Kiev, Ukraine
| |
Collapse
|
25
|
Soni A, Augsburg M, Buchholz F, Pisabarro MT. Nearest-neighbor amino acids of specificity-determining residues influence the activity of engineered Cre-type recombinases. Sci Rep 2020; 10:13985. [PMID: 32814809 PMCID: PMC7438526 DOI: 10.1038/s41598-020-70867-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2020] [Accepted: 08/03/2020] [Indexed: 11/29/2022] Open
Abstract
The tyrosine-type site-specific DNA recombinase Cre recombines its target site, loxP, with high activity and specificity without cross-recombining the target sites of highly related recombinases. Understanding how Cre achieves this precision is key to be able to rationally engineer site-specific recombinases (SSRs) for genome editing applications. Previous work has revealed key residues for target site selectivity in the Cre/loxP and the related Dre/rox recombinase systems. However, enzymes in which these residues were changed to the respective counterpart only showed weak activity on the foreign target site. Here, we use molecular modeling and dynamics simulation techniques to comprehensively explore the mechanisms by which these residues determine target recognition in the context of their flanking regions in the protein-DNA interface, and we establish a structure-based rationale for the design of improved recombination activities. Our theoretical models reveal that nearest-neighbors to the specificity-determining residues are important players for enhancing SSR activity on the foreign target site. Based on the established rationale, we design new Cre variants with improved rox recombination activities, which we validate experimentally. Our work provides new insights into the target recognition mechanisms of Cre-like recombinases and represents an important step towards the rational design of SSRs for applied genome engineering.
Collapse
Affiliation(s)
- Anjali Soni
- Structural Bioinformatics, BIOTEC, TU Dresden, Tatzberg 47-51, 01307, Dresden, Germany
| | - Martina Augsburg
- University Carl Gustav Carus and Medical Faculty, UCC, Medical Systems Biology, TU Dresden, Fetscherstrasse 74, Dresden, Germany
| | - Frank Buchholz
- University Carl Gustav Carus and Medical Faculty, UCC, Medical Systems Biology, TU Dresden, Fetscherstrasse 74, Dresden, Germany
| | - M Teresa Pisabarro
- Structural Bioinformatics, BIOTEC, TU Dresden, Tatzberg 47-51, 01307, Dresden, Germany.
| |
Collapse
|
26
|
|
27
|
Salsbury AM, Dean TJ, Lemkul JA. Polarizable Molecular Dynamics Simulations of Two c-kit Oncogene Promoter G-Quadruplexes: Effect of Primary and Secondary Structure on Loop and Ion Sampling. J Chem Theory Comput 2020; 16:3430-3444. [PMID: 32307997 DOI: 10.1021/acs.jctc.0c00191] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
G-quadruplexes (GQs) are highly ordered nucleic acid structures that play fundamental roles in regulating gene expression and maintaining genomic stability. GQs are topologically diverse and enriched in promoter sequences of growth regulatory genes and proto-oncogenes, suggesting that they may serve as attractive targets for drug design at the level of transcription rather than inhibiting the activity of the protein products of these genes. The c-kit promoter contains three adjacent GQ-forming sequences that have proposed antagonistic effects on gene expression and thus are promising drug targets for diseases such as gastrointestinal stromal tumors, mast cell disease, and leukemia. Because GQ stability is influenced by primary structure, secondary structure, and ion interactions, a greater understanding of GQ structure, dynamics, and ion binding properties is needed to develop novel, GQ-targeting therapeutics. Here, we performed molecular dynamics simulations to systematically study the c-kit2 and c-kit* GQs, evaluating nonpolarizable and polarizable force fields (FFs) and examining the effects of base substitutions and cation type (K+, Na+, and Li+) on the dynamics of their isolated and linked structures. We found that the Drude polarizable FF outperformed the additive CHARMM36 FF in two- and three-tetrad GQs and solutions of KCl, NaCl, and LiCl. Drude simulations with different cations agreed with the known GQ stabilization preference (K+ > Na+ > Li+) and illustrated that tetrad core-ion coordination differs as a function of cation type. Finally, we showed that differences in primary and secondary structure influence loop sampling, ion binding, and core-ion energetics of GQs.
Collapse
Affiliation(s)
- Alexa M Salsbury
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Tanner J Dean
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Justin A Lemkul
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
28
|
Aqueous Contact Ion Pairs of Phosphate Groups with Na+, Ca2+ and Mg2+ – Structural Discrimination by Femtosecond Infrared Spectroscopy and Molecular Dynamics Simulations. Z PHYS CHEM 2020. [DOI: 10.1515/zpch-2020-1614] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Abstract
The extent of contact and solvent shared ion pairs of phosphate groups with Na+, Ca2+ and Mg2+ ions in aqueous environment and their relevance for the stability of polyanionic DNA and RNA structures is highly debated. Employing the asymmetric phosphate stretching vibration of dimethyl phosphate (DMP), a model system of the sugar-phosphate backbone of DNA and RNA, we present linear infrared, femtosecond infrared pump-probe and absorptive 2D-IR spectra that report on contact ion pair formation via the presence of blue shifted spectral signatures. Compared to the linear infrared spectra, the nonlinear spectra reveal contact ion pairs with increased sensitivity because the spectra accentuate differences in peak frequency, transition dipole moment strength, and excited state lifetime. The experimental results are corroborated by long time scale MD simulations, benchmarked by density functional simulations on phosphate-ion-water clusters. The microscopic interpretation reveals subtle structural differences of ion pairs formed by the phosphate group and the ions Na+, Ca2+ and Mg2+. Intricate properties of the solvation shell around the phosphate group and the ion are essential to explain the experimental observations. The present work addresses a challenging to probe topic with the help of a model system and establishes new experimental data of contact ion pair formation, thereby underlining the potential of nonlinear 2D-IR spectroscopy as an analytical probe of phosphate-ion interactions in complex biological systems.
Collapse
|
29
|
Perepelytsya S, Uličný J, Laaksonen A, Mocci F. Pattern preferences of DNA nucleotide motifs by polyamines putrescine2+, spermidine3+ and spermine4. Nucleic Acids Res 2020; 47:6084-6097. [PMID: 31114917 PMCID: PMC6614828 DOI: 10.1093/nar/gkz434] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2018] [Revised: 04/10/2019] [Accepted: 05/08/2019] [Indexed: 12/25/2022] Open
Abstract
The interactions of natural polyamines (putrescine2+, spermidine3+ and spermine4+) with DNA double helix are studied to characterize their nucleotide sequence pattern preference. Atomistic Molecular Dynamics simulations have been carried out for three systems consisting of the same DNA fragment d(CGCGAATTCGCGAATTCGCG) with different polyamines. The results show that polyamine molecules are localized with well-recognized patterns along the double helix with different residence times. We observed a clear hierarchy in the residence times of the polyamines, with the longest residence time (ca 100ns) in the minor groove. The analysis of the sequence dependence shows that polyamine molecules prefer the A-tract regions of the minor groove - in its narrowest part. The preferable localization of putrescine2+, spermidine3+ and spermine4+ in the minor groove with A-tract motifs is correlated with modulation of the groove width by a specific nucleotide sequences. We did develop a theoretical model pointing to the electrostatic interactions as the main driving force in this phenomenon, making it even more prominent for polyamines with higher charges. The results of the study explain the specificity of polyamine interactions with A-tract region of the DNA double helix which is also observed in experiments.
Collapse
Affiliation(s)
- Sergiy Perepelytsya
- Bogolyubov Institute for Theoretical Physics of the National Academy of Sciences of Ukraine, 03143 Kyiv, Ukraine.,Department of Theoretical and Mathematical Physics, Kyiv Academic University, 03142 Kyiv, Ukraine
| | - Jozef Uličný
- Department of Biophysics, Institute of Physics, P. J. Šafárik University, 041 54 Košice, Slovakia
| | - Aatto Laaksonen
- State Key Laboratory of Materials-Oriented and Chemical Engineering, Nanjing Tech University, 210009 Nanjing, China.,Division of Physical Chemistry, Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, 10691 Stockholm, Sweden.,Centre of Advanced Research in Bionanoconjugates and Biopolymers, Petru Poni Institute of Macromolecular Chemistry, Iasi, 700487, Romania
| | - Francesca Mocci
- Centre of Advanced Research in Bionanoconjugates and Biopolymers, Petru Poni Institute of Macromolecular Chemistry, Iasi, 700487, Romania.,Department of Chemical and Geological Sciences, University of Cagliari, I-09042 Monserrato, Italy
| |
Collapse
|
30
|
Schauss J, Kundu A, Fingerhut BP, Elsaesser T. Contact Ion Pairs of Phosphate Groups in Water: Two-Dimensional Infrared Spectroscopy of Dimethyl Phosphate and ab Initio Simulations. J Phys Chem Lett 2019; 10:6281-6286. [PMID: 31560211 DOI: 10.1021/acs.jpclett.9b02157] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The interaction of phosphate groups with ions in an aqueous environment has a strong impact on the structure and folding processes of DNA and RNA. The dynamic variety of ionic arrangements, including both contact pairs and water separated ions, and the molecular coupling mechanisms are far from being understood. In a combined experimental and theoretical approach, we address the properties of contact ion pairs of the prototypical system dimethyl phosphate with Na+, Ca2+, and Mg2+ ions in water. Linear and femtosecond two-dimensional infrared (2D-IR) spectroscopy of the asymmetric (PO2)- stretching vibration separates and characterizes the different species via their blue-shifted vibrational signatures and 2D-IR line shapes. Phosphate-magnesium contact pairs stand out as the most compact geometry while the contact pairs with Ca2+ and Na+ display a wider structural variation. Microscopic density functional theory simulations rationalize the observed frequency shifts and reveal distinct differences between the contact geometries.
Collapse
Affiliation(s)
- Jakob Schauss
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| | - Achintya Kundu
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| |
Collapse
|
31
|
Galindo-Murillo R, Cheatham TE. Lessons learned in atomistic simulation of double-stranded DNA: Solvation and salt concerns [Article v1.0]. ACTA ACUST UNITED AC 2019; 1. [PMID: 33073182 DOI: 10.33011/livecoms.1.2.9974] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Nucleic acids are highly charged macromolecules sensitive to their surroundings of water, salt, and other biomolecules. Molecular dynamics simulations with accurate biomolecular force fields provide a detailed atomistic view into DNA and RNA that has been useful to study the structure and dynamics of these molecules and their biological relevance. In this work we study the Drew-Dickerson dodecamer duplex with the sequence d(GCGCAATTGCGC)2 in three different salt concentrations and using different monvalent salt types to detect possible structural influence. Overall, the DNA shows no major structural changes regardless of amount or type of monovalent ions used. Our results show that only at very high salt conditions (5M) is a small structural effect observed in the DNA duplex, which mainly consist of narrowing of the grooves due to increased residence of ions. We also present the importance of sampling time to achieve a converged ensemble, which is of major relevance in any simulation to avoid biased or non-meaningful results.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, L. S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, L. S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112
| |
Collapse
|
32
|
Ghoshdastidar D, Senapati S. Dehydrated DNA in B-form: ionic liquids in rescue. Nucleic Acids Res 2019; 46:4344-4353. [PMID: 29669113 PMCID: PMC5961314 DOI: 10.1093/nar/gky253] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2017] [Accepted: 04/10/2018] [Indexed: 12/19/2022] Open
Abstract
The functional B-conformation of DNA succumbs to the A-form at low water activity. Methods for room temperature DNA storage that rely upon ‘anhydrobiosis’, thus, often encounter the loss of DNA activity due to the B→A-DNA transition. Here, we show that ionic liquids, an emerging class of green solvents, can induce conformational transitions in DNA and even rescue the dehydrated DNA in the functional B-form. CD spectroscopic analyses not only reveal rapid transition of A-DNA in 78% ethanol medium to B-conformation in presence of ILs, but also the high resistance of IL-bound B-form to transit to A-DNA under dehydration. Molecular dynamics simulations show the unique ability of ILs to disrupt Na+ ion condensation and form ‘IL spine’ in DNA minor groove to drive the A→B transition. Implications of these findings range from the plausible use of ILs as novel anhydrobiotic DNA storage medium to a switch for modulating DNA conformational transitions.
Collapse
Affiliation(s)
- Debostuti Ghoshdastidar
- Department of Biotechnology, BJM School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Sanjib Senapati
- Department of Biotechnology, BJM School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
33
|
|
34
|
Carter EK, Laughlin-Toth S, Dodd T, Wilson WD, Ivanov I. Small molecule binders recognize DNA microstructural variations via an induced fit mechanism. Phys Chem Chem Phys 2019; 21:1841-1851. [PMID: 30629058 PMCID: PMC6497476 DOI: 10.1039/c8cp05537h] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Regulation of gene-expression by specific targeting of protein-nucleic acid interactions has been a long-standing goal in medicinal chemistry. Transcription factors are considered "undruggable" because they lack binding sites well suited for binding small-molecules. In order to overcome this obstacle, we are interested in designing small molecules that bind to the corresponding promoter sequences and either prevent or modulate transcription factor association via an allosteric mechanism. To achieve this, we must design small molecules that are both sequence-specific and able to target G/C base pair sites. A thorough understanding of the relationship between binding affinity and the structural aspects of the small molecule-DNA complex would greatly aid in rational design of such compounds. Here we present a comprehensive analysis of sequence-specific DNA association of a synthetic minor groove binder using long timescale molecular dynamics. We show how binding selectivity arises from a combination of structural factors. Our results provide a framework for the rational design and optimization of synthetic small molecules in order to improve site-specific targeting of DNA for therapeutic uses in the design of selective DNA binders targeting transcription regulation.
Collapse
Affiliation(s)
- E. Kathleen Carter
- Department of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA. ;
| | - Sarah Laughlin-Toth
- Department of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA. ;
- Department of Pediatrics, Aflac Cancer and Blood Disorders Center, Emory University School of Medicine, Atlanta, GA, USA
| | - Thomas Dodd
- Department of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA. ;
| | - W. David Wilson
- Department of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA. ;
| | - Ivaylo Ivanov
- Department of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA. ;
| |
Collapse
|
35
|
Schauss J, Dahms F, Fingerhut BP, Elsaesser T. Phosphate-Magnesium Ion Interactions in Water Probed by Ultrafast Two-Dimensional Infrared Spectroscopy. J Phys Chem Lett 2019; 10:238-243. [PMID: 30599134 DOI: 10.1021/acs.jpclett.8b03568] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Electric interactions between ions and ionic molecular groups in aqueous solution play a fundamental role in chemistry and biology. While Mg2+ ions are known to strongly affect the structure and folding dynamics of biomolecules, the relevance of different solvation geometries and the underlying interactions are mainly unresolved. We study dynamics and couplings between the hydrated Mg2+ and the dimethylphosphate anion, an established model system for the DNA and RNA backbone. The asymmetric (PO2-) stretching vibration serves as a sensitive noninvasive probe of phosphate-ion interactions. Femtosecond two-dimensional infrared (2D-IR) spectroscopy directly maps Mg2+ ions in contact with the phosphate groups via a distinct blue-shifted signature in the 2D spectrum. Data for different Mg2+ concentrations are analyzed by microscopic density functional theory modeling of cluster geometries and associated spectroscopic features, providing spatial assignments of the observed 2D-IR signatures. Phosphate-ion interactions arising from electrostatic Coulomb forces and exchange repulsion are the predominant origin of the observed frequency shifts.
Collapse
Affiliation(s)
- Jakob Schauss
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| | - Fabian Dahms
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| |
Collapse
|
36
|
Bochicchio A, Krepl M, Yang F, Varani G, Sponer J, Carloni P. Molecular basis for the increased affinity of an RNA recognition motif with re-engineered specificity: A molecular dynamics and enhanced sampling simulations study. PLoS Comput Biol 2018; 14:e1006642. [PMID: 30521520 PMCID: PMC6307825 DOI: 10.1371/journal.pcbi.1006642] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2018] [Revised: 12/27/2018] [Accepted: 11/13/2018] [Indexed: 12/25/2022] Open
Abstract
The RNA recognition motif (RRM) is the most common RNA binding domain across eukaryotic proteins. It is therefore of great value to engineer its specificity to target RNAs of arbitrary sequence. This was recently achieved for the RRM in Rbfox protein, where four mutations R118D, E147R, N151S, and E152T were designed to target the precursor to the oncogenic miRNA 21. Here, we used a variety of molecular dynamics-based approaches to predict specific interactions at the binding interface. Overall, we have run approximately 50 microseconds of enhanced sampling and plain molecular dynamics simulations on the engineered complex as well as on the wild-type Rbfox·pre-miRNA 20b from which the mutated systems were designed. Comparison with the available NMR data on the wild type molecules (protein, RNA, and their complex) served to establish the accuracy of the calculations. Free energy calculations suggest that further improvements in affinity and selectivity are achieved by the S151T replacement. RNA is an outstanding target for oncological intervention. Engineering the most common RNA binding motif in human proteins (called RRM) so as to bind to a specific RNA has an enormous pharmacological potential. Yet, it is highly non trivial to design RRM-bearing protein variants with RNA selectivity and affinity sufficiently high for clinical applications. Here we present an extensive molecular simulation study which shed light on the exquisite molecular recognition of the empirically-engineered complex between the RRM-bearing protein Rbfox and its RNA target pre-miR21. The simulations allow predicting a variant, the S151T, which may lead to further enhancement of selectivity and affinity for pre-miR21.
Collapse
Affiliation(s)
- Anna Bochicchio
- Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich, Jülich, Germany
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Brno, Czech Republic
- * E-mail: (MK); (PC)
| | - Fan Yang
- Department of Chemistry, University of Washington, Seattle, Washington, United States of America
| | - Gabriele Varani
- Department of Chemistry, University of Washington, Seattle, Washington, United States of America
| | - Jiri Sponer
- Institute of Biophysics of the Czech Academy of Sciences, Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, Olomouc, Czech Republic
| | - Paolo Carloni
- Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich, Jülich, Germany
- JARA-HPC, Jülich Supercomputing Centre, Forschungszentrum Jülich GmbH, Jülich, Germany
- * E-mail: (MK); (PC)
| |
Collapse
|
37
|
Jaiswal AK, Srivastava R, Pandey P, Bandyopadhyay P. Microscopic picture of water-ethylene glycol interaction near a model DNA by computer simulation: Concentration dependence, structure, and localized thermodynamics. PLoS One 2018; 13:e0206359. [PMID: 30427849 PMCID: PMC6235303 DOI: 10.1371/journal.pone.0206359] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2018] [Accepted: 10/11/2018] [Indexed: 01/03/2023] Open
Abstract
It is known that crowded molecular environment affects the structure, thermodynamics, and dynamics of macromolecules. Most of the previous works on molecular crowding have majorly focused on the behavior of the macromolecule with less emphasis on the behavior of the crowder and water molecules. In the current study, we have precisely focused on the behavior of the crowder, (ethylene glycol (EG)), salt ions, and water in the presence of a DNA with the increase of the EG concentration. We have probed the behavior of water and crowder using molecular dynamics (MD) simulation and by calculating localized thermodynamic properties. Our results show an interesting competition between EG and water molecules to make hydrogen bonds (H-bond) with DNA. Although the total number of H-bonds involving DNA with both EG and water remains essentially same irrespective of the increase in EG concentration, there is a proportional change in the H-bonding pattern between water-water, EG-EG, and EG-water near DNA and in bulk. At low concentrations of EG, the displacement of water molecules near DNA is relatively easy. However, the displacement of water becomes more difficult as the concentration of EG increases. The density of Na+ (Cl-) near DNA increases (decreases) as the concentration of EG is increased. The density of Cl- near Na+ increases with the increase in EG concentration. It was also found that the average free energy per water in the first solvation shell increases with the increase in EG concentration. Putting all these together, a microscopic picture of EG, water, salt interaction in the presence of DNA, as a function of EG concentration, has emerged.
Collapse
Affiliation(s)
- Atul Kumar Jaiswal
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Rakesh Srivastava
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Preeti Pandey
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
- * E-mail:
| |
Collapse
|
38
|
Francés-Monerris A, Gattuso H, Roca-Sanjuán D, Tuñón I, Marazzi M, Dumont E, Monari A. Dynamics of the excited-state hydrogen transfer in a (dG)·(dC) homopolymer: intrinsic photostability of DNA. Chem Sci 2018; 9:7902-7911. [PMID: 30450180 PMCID: PMC6202918 DOI: 10.1039/c8sc03252a] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2018] [Accepted: 09/17/2018] [Indexed: 12/12/2022] Open
Abstract
The intrinsic photostability of nucleic acids is intimately related to evolution of life, while its understanding at the molecular and electronic levels remains a challenge for modern science. Among the different decay pathways proposed in the last two decades, the excited-state hydrogen transfer between guanine-cytosine base pairs has been identified as an efficient non-reactive channel to dissipate the excess of energy provided by light absorption. The present work studies the dynamics of such phenomena taking place in a (dG)·(dC) B-DNA homopolymer in water solution using state-of-the-art molecular modelling and simulation methods. A dynamic effect that boosts the photostability of the inter-strand hydrogen atom transfers, inherent to the Watson-Crick base pairing, is unveiled and ascribed to the energy released during the proton transfer step. Our results also reveal a novel mechanism of DNA decay named four proton transfer (FPT), in which two protons of two adjacent G-C base pairs are transferred to form a biradical zwitterionic intermediate. Decay of the latter intermediate to the ground state triggers the transfer of the protons back to the guanine molecules recovering the Watson-Crick structure of the tetramer. This FPT process is activated by the close interaction of a nearby Na+ counterion with the oxygen atoms of the guanine nucleobases and hence represents a photostable channel operative in natural nucleic acids.
Collapse
Affiliation(s)
| | - Hugo Gattuso
- Theoretical Physical Chemistry , Research Unit Molecular Systems (UR MOLSYS) , University of Liège , 4000 Liège , Belgium
| | - Daniel Roca-Sanjuán
- Instituto de Ciencia Molecular , Universitat de València , Apartado 22085 , ES-46071 Paterna , Spain
| | - Iñaki Tuñón
- Departamento de Química Física , Universitat de València , 46100 Burjassot , Spain
| | - Marco Marazzi
- Departamento de Química , Centro de Investigación en Síntesis Química (CISQ) , Universidad de La Rioja , 26006 Logroño , Spain
| | - Elise Dumont
- Univ. Lyon , ENS de Lyon , CNRS UMR 5182 , Université Claude Bernard Lyon 1 , Laboratoire de Chimie , F69342 , Lyon , France
| | - Antonio Monari
- Université de Lorraine , CNRS , LPCT , Nancy F-54000 , France . ;
| |
Collapse
|
39
|
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
|
40
|
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]
|
41
|
Lowe BM, Skylaris CK, Green NG, Shibuta Y, Sakata T. Molecular dynamics simulation of potentiometric sensor response: the effect of biomolecules, surface morphology and surface charge. NANOSCALE 2018; 10:8650-8666. [PMID: 29700545 DOI: 10.1039/c8nr00776d] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The silica-water interface is critical to many modern technologies in chemical engineering and biosensing. One technology used commonly in biosensors, the potentiometric sensor, operates by measuring the changes in electric potential due to changes in the interfacial electric field. Predictive modelling of this response caused by surface binding of biomolecules remains highly challenging. In this work, through the most extensive molecular dynamics simulation of the silica-water interfacial potential and electric field to date, we report a novel prediction and explanation of the effects of nano-morphology on sensor response. Amorphous silica demonstrated a larger potentiometric response than an equivalent crystalline silica model due to increased sodium adsorption, in agreement with experiments showing improved sensor response with nano-texturing. We provide proof-of-concept that molecular dynamics can be used as a complementary tool for potentiometric biosensor response prediction. Effects that are conventionally neglected, such as surface morphology, water polarisation, biomolecule dynamics and finite-size effects, are explicitly modelled.
Collapse
Affiliation(s)
- B M Lowe
- Department of Materials Engineering, School of Engineering, The University of Tokyo 7-3-1 Hongo, Bunkyo, Tokyo 113-8656, Japan.
| | | | | | | | | |
Collapse
|
42
|
O' Lee DJ. Introducing a model of pairing based on base pair specific interactions between identical DNA sequences. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:075102. [PMID: 29219116 DOI: 10.1088/1361-648x/aaa043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
At present, there have been suggested two types of physical mechanism that may facilitate preferential pairing between DNA molecules, with identical or similar base pair texts, without separation of base pairs. One mechanism solely relies on base pair specific patterns of helix distortion being the same on the two molecules, discussed extensively in the past. The other mechanism proposes that there are preferential interactions between base pairs of the same composition. We introduce a model, built on this second mechanism, where both thermal stretching and twisting fluctuations are included, as well as the base pair specific helix distortions. Firstly, we consider an approximation for weak pairing interactions, or short molecules. This yields a dependence of the energy on the square root of the molecular length, which could explain recent experimental data. However, analysis suggests that this approximation is no longer valid at large DNA lengths. In a second approximation, for long molecules, we define two adaptation lengths for twisting and stretching, over which the pairing interaction can limit the accumulation of helix disorder. When the pairing interaction is sufficiently strong, both adaptation lengths are finite; however, as we reduce pairing strength, the stretching adaptation length remains finite but the torsional one becomes infinite. This second state persists to arbitrarily weak values of the pairing strength; suggesting that, if the molecules are long enough, the pairing energy scales as length. To probe differences between the two pairing mechanisms, we also construct a model of similar form. However, now, pairing between identical sequences solely relies on the intrinsic helix distortion patterns. Between the two models, we see interesting qualitative differences. We discuss our findings, and suggest new work to distinguish between the two mechanisms.
Collapse
Affiliation(s)
- Dominic J O' Lee
- Department of Chemistry, Imperial College London, SW7 2AZ, London, United Kingdom
| |
Collapse
|
43
|
Grad JN, Gigante A, Wilms C, Dybowski JN, Ohl L, Ottmann C, Schmuck C, Hoffmann D. Locating Large, Flexible Ligands on Proteins. J Chem Inf Model 2018; 58:315-327. [PMID: 29266929 DOI: 10.1021/acs.jcim.7b00413] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Many biologically important ligands of proteins are large, flexible, and in many cases charged molecules that bind to extended regions on the protein surface. It is infeasible or expensive to locate such ligands on proteins with standard methods such as docking or molecular dynamics (MD) simulation. The alternative approach proposed here is scanning of a spatial and angular grid around the protein with smaller fragments of the large ligand. Energy values for complete grids can be computed efficiently with a well-known fast Fourier transform-accelerated algorithm and a physically meaningful interaction model. We show that the approach can readily incorporate flexibility of the protein and ligand. The energy grids (EGs) resulting from the ligand fragment scans can be transformed into probability distributions and then directly compared to probability distributions estimated from MD simulations and experimental structural data. We test the approach on a diverse set of complexes between proteins and large, flexible ligands, including a complex of sonic hedgehog protein and heparin, three heparin sulfate substrates or nonsubstrates of an epimerase, a multibranched supramolecular ligand that stabilizes a protein-peptide complex, a flexible zwitterionic ligand that binds to a surface basin of a Kringle domain, and binding of ATP to a flexible site of an ion channel. In all cases, the EG approach gives results that are in good agreement with experimental data or MD simulations.
Collapse
Affiliation(s)
- Jean-Noël Grad
- Bioinformatics and Computational Biophysics, Faculty of Biology, University of Duisburg-Essen , Universitätstraße 7, 45141 Essen, Germany
| | - Alba Gigante
- Institute of Organic Chemistry, University of Duisburg-Essen , Universitätstraße 7, 45141 Essen, Germany
| | - Christoph Wilms
- Bioinformatics and Computational Biophysics, Faculty of Biology, University of Duisburg-Essen , Universitätstraße 7, 45141 Essen, Germany
| | - Jan Nikolaj Dybowski
- Bioinformatics and Computational Biophysics, Faculty of Biology, University of Duisburg-Essen , Universitätstraße 7, 45141 Essen, Germany
| | - Ludwig Ohl
- Bioinformatics and Computational Biophysics, Faculty of Biology, University of Duisburg-Essen , Universitätstraße 7, 45141 Essen, Germany
| | - Christian Ottmann
- Laboratory of Chemical Biology, Department of Biomedical Engineering, and Institute for Complex Molecular Systems, Eindhoven University of Technology , Den Dolech 2, 5612 AZ Eindhoven, The Netherlands
| | - Carsten Schmuck
- Institute of Organic Chemistry, University of Duisburg-Essen , Universitätstraße 7, 45141 Essen, Germany
| | - Daniel Hoffmann
- Bioinformatics and Computational Biophysics, Faculty of Biology, University of Duisburg-Essen , Universitätstraße 7, 45141 Essen, Germany
| |
Collapse
|
44
|
Etheve L, Martin J, Lavery R. Decomposing protein-DNA binding and recognition using simplified protein models. Nucleic Acids Res 2017; 45:10270-10283. [PMID: 28973439 PMCID: PMC5622342 DOI: 10.1093/nar/gkx627] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2017] [Accepted: 07/11/2017] [Indexed: 12/30/2022] Open
Abstract
We analyze the role of different physicochemical factors in protein/DNA binding and recognition by comparing the results from all-atom molecular dynamics simulations with simulations using simplified protein models. These models enable us to separate the role of specific amino acid side chains, formal amino acid charges and hydrogen bonding from the effects of the low-dielectric volume occupied by the protein. Comparisons are made on the basis of the conformation of DNA after protein binding, the ionic distribution around the complex and the sequence specificity. The results for four transcription factors, binding in either the minor or major grooves of DNA, show that the protein volume and formal charges, with one exception, play a predominant role in binding. Adding hydrogen bonding and a very small number of key amino acid side chains at the all-atom level yields results in DNA conformations and sequence recognition close to those seen in the reference all-atom simulations.
Collapse
Affiliation(s)
- Loïc Etheve
- MMSB UMR 5086 CNRS / Univ. Lyon I, Institut de Biologie et Chimie des Protéines, 7 passage du Vercors, Lyon 69367, France
| | - Juliette Martin
- MMSB UMR 5086 CNRS / Univ. Lyon I, Institut de Biologie et Chimie des Protéines, 7 passage du Vercors, Lyon 69367, France
| | - Richard Lavery
- MMSB UMR 5086 CNRS / Univ. Lyon I, Institut de Biologie et Chimie des Protéines, 7 passage du Vercors, Lyon 69367, France
| |
Collapse
|
45
|
Pasi M, Zakrzewska K, Maddocks JH, Lavery R. Analyzing DNA curvature and its impact on the ionic environment: application to molecular dynamics simulations of minicircles. Nucleic Acids Res 2017; 45:4269-4277. [PMID: 28180333 PMCID: PMC5397150 DOI: 10.1093/nar/gkx092] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2016] [Accepted: 02/07/2017] [Indexed: 01/16/2023] Open
Abstract
We propose a method for analyzing the magnitude and direction of curvature within nucleic acids, based on the curvilinear helical axis calculated by Curves+. The method is applied to analyzing curvature within minicircles constructed with varying degrees of over- or under-twisting. Using the molecular dynamics trajectories of three different minicircles, we are able to quantify how curvature varies locally both in space and in time. We also analyze how curvature influences the local environment of the minicircles, notably via increased heterogeneity in the ionic distributions surrounding the double helix. The approach we propose has been integrated into Curves+ and the utilities Canal (time trajectory analysis) and Canion (environmental analysis) and can be used to study a wide variety of static and dynamic structural data on nucleic acids.
Collapse
Affiliation(s)
- Marco Pasi
- MMSB, University of Lyon I/CNRS UMR 5086, Institut de Biologie et Chimie des Protéines, 7 Passage du Vercors, 69367 Lyon, France
| | - Krystyna Zakrzewska
- MMSB, University of Lyon I/CNRS UMR 5086, Institut de Biologie et Chimie des Protéines, 7 Passage du Vercors, 69367 Lyon, France
| | - John H Maddocks
- Section de Mathématiques, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Richard Lavery
- MMSB, University of Lyon I/CNRS UMR 5086, Institut de Biologie et Chimie des Protéines, 7 Passage du Vercors, 69367 Lyon, France
| |
Collapse
|
46
|
Abstract
The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.
Collapse
Affiliation(s)
- Damien Laage
- École
Normale Supérieure, PSL Research University, UPMC Univ Paris
06, CNRS, Département de Chimie,
PASTEUR, 24 rue Lhomond, 75005 Paris, France
- Sorbonne
Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR, 75005 Paris, France
| | - Thomas Elsaesser
- Max-Born-Institut
für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| | - James T. Hynes
- École
Normale Supérieure, PSL Research University, UPMC Univ Paris
06, CNRS, Département de Chimie,
PASTEUR, 24 rue Lhomond, 75005 Paris, France
- Sorbonne
Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR, 75005 Paris, France
- Department
of Chemistry and Biochemistry, University
of Colorado, Boulder, Colorado 80309, United
States
| |
Collapse
|
47
|
Lee DJO. Statistical mechanical model for a closed loop plectoneme with weak helix specific forces. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:145101. [PMID: 28251958 DOI: 10.1088/1361-648x/aa521c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We develop a statistical mechanical framework, based on a variational approximation, to describe closed loop plectonemes. This framework incorporates weak helix structure dependent forces into the determination of the free energy and average structure of a plectoneme. Notably, due to their chiral nature, helix structure dependent forces break the symmetry between left and right handed supercoiling. The theoretical approach, presented here, also provides a systematic way of enforcing the topological constraint of closed loop supercoiling in the variational approximation. At large plectoneme lengths, by considering correlation functions in an expansion in terms of the spatial mean twist density about its thermally averaged value, it can be argued that topological constraint may be approximated by replacing twist and writhe by their thermal averages. A Lagrange multiplier, containing the sum of average twist and writhe, can be added to the free energy to conveniently inforce this result. The average writhe can be calculated through the thermal average of the Gauss' integral in the variational approximation. Furthermore, this approach allows for a possible way to calculate finite size corrections due to the topological constraint. Using interaction energy terms from the mean-field Kornyshev-Leikin theory, for parameter values that correspond to weak helix dependent forces, we calculate the free energy, fluctuation magnitudes and mean geometric parameters for the plectoneme. We see a slight asymmetry, where interestingly, left handed supercoils have a looser structure than right handed ones, although with a lower free energy, unlike what the previous ground state calculations would suggest.
Collapse
Affiliation(s)
- Dominic J O' Lee
- Department of Chemistry, Imperial College London, SW7 2AZ, London, United Kingdom
| |
Collapse
|
48
|
Andrews CT, Campbell BA, Elcock AH. Direct Comparison of Amino Acid and Salt Interactions with Double-Stranded and Single-Stranded DNA from Explicit-Solvent Molecular Dynamics Simulations. J Chem Theory Comput 2017; 13:1794-1811. [PMID: 28288277 DOI: 10.1021/acs.jctc.6b00883] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Given the ubiquitous nature of protein-DNA interactions, it is important to understand the interaction thermodynamics of individual amino acid side chains for DNA. One way to assess these preferences is to perform molecular dynamics (MD) simulations. Here we report MD simulations of 20 amino acid side chain analogs interacting simultaneously with both a 70-base-pair double-stranded DNA and with a 70-nucleotide single-stranded DNA. The relative preferences of the amino acid side chains for dsDNA and ssDNA match well with values deduced from crystallographic analyses of protein-DNA complexes. The estimated apparent free energies of interaction for ssDNA, on the other hand, correlate well with previous simulation values reported for interactions with isolated nucleobases, and with experimental values reported for interactions with guanosine. Comparisons of the interactions with dsDNA and ssDNA indicate that, with the exception of the positively charged side chains, all types of amino acid side chain interact more favorably with ssDNA, with intercalation of aromatic and aliphatic side chains being especially notable. Analysis of the data on a base-by-base basis indicates that positively charged side chains, as well as sodium ions, preferentially bind to cytosine in ssDNA, and that negatively charged side chains, and chloride ions, preferentially bind to guanine in ssDNA. These latter observations provide a novel explanation for the lower salt dependence of DNA duplex stability in GC-rich sequences relative to AT-rich sequences.
Collapse
Affiliation(s)
- Casey T Andrews
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Brady A Campbell
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Adrian H Elcock
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| |
Collapse
|
49
|
Pasi M, Angelov D, Bednar J, Dimitrov S, Lavery R. Extra views on structure and dynamics of DNA loops on nucleosomes studied with molecular simulations. Nucleus 2016; 7:554-559. [PMID: 27874316 PMCID: PMC5214536 DOI: 10.1080/19491034.2016.1260800] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2016] [Accepted: 11/10/2016] [Indexed: 02/07/2023] Open
Abstract
It has been shown experimentally that the action of the RSC chromatin remodeler leads to the formation of an irregular, partially remodeled nucleosome, termed a remosome. The remosome contains an extra 30-40 base pairs of DNA compared to a canonical nucleosome. Large-scale molecular simulations have provided information on the probable structure of remosomes and have explained why they remain stable in the absence of RSC. Here we explain how these simulations were carried out and what the resulting remosome models imply in terms of the mechanism of action of RSC. We notably show that local kinks within DNA are key in explaining how extra DNA can be in added to nucleosomes without unduly disturbing DNA-histone binding.
Collapse
Affiliation(s)
- Marco Pasi
- Université de Lyon/CNRS UMR 5086, MMSB, Institut de Biologie et Chimie des Protéines, Lyon, France
| | - Dimitar Angelov
- Université de Lyon/CNRS/INRA UMR 5239, LBMC, Ecole Normale Supérieure de Lyon, Lyon, France
| | - Jan Bednar
- Institut Albert Bonniot, Université de Grenoble Alpes/INSERM U1209/CNRS UMR 5309, Grenoble, France
| | - Stefan Dimitrov
- Institut Albert Bonniot, Université de Grenoble Alpes/INSERM U1209/CNRS UMR 5309, Grenoble, France
| | - Richard Lavery
- Université de Lyon/CNRS UMR 5086, MMSB, Institut de Biologie et Chimie des Protéines, Lyon, France
| |
Collapse
|
50
|
Etheve L, Martin J, Lavery R. Protein-DNA interfaces: a molecular dynamics analysis of time-dependent recognition processes for three transcription factors. Nucleic Acids Res 2016; 44:9990-10002. [PMID: 27658967 PMCID: PMC5175364 DOI: 10.1093/nar/gkw841] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2016] [Revised: 09/12/2016] [Accepted: 09/13/2016] [Indexed: 12/14/2022] Open
Abstract
We have studied the dynamics of three transcription factor-DNA complexes using all-atom, microsecond-scale MD simulations. In each case, the salt bridges and hydrogen bond interactions formed at the protein-DNA interface are found to be dynamic, with lifetimes typically in the range of tens to hundreds of picoseconds, although some interactions, notably those involving specific binding to DNA bases, can be a hundred times longer lived. Depending on the complex studied, this dynamics may or may not lead to the existence of distinct conformational substates. Using a sequence threading technique, it has been possible to determine whether DNA sequence recognition is sensitive or not to such conformational changes, and, in one case, to show that recognition appears to be locally dependent on protein-mediated cation distributions.
Collapse
Affiliation(s)
- Loïc Etheve
- MMSB UMR 5086 CNRS/University of Lyon I, Institut de Biologie et Chimie des Protéines, 7 passage du Vercors, Lyon 69367, France
| | - Juliette Martin
- MMSB UMR 5086 CNRS/University of Lyon I, Institut de Biologie et Chimie des Protéines, 7 passage du Vercors, Lyon 69367, France
| | - Richard Lavery
- MMSB UMR 5086 CNRS/University of Lyon I, Institut de Biologie et Chimie des Protéines, 7 passage du Vercors, Lyon 69367, France
| |
Collapse
|