1
|
Tripathi M, Thakur Y, Syed R, Asatkar AK, Alqahtani MS, Das D, Agrawal R, Verma B, Pande R. In-vitro and in-silico analysis and antitumor studies of novel Cu(II) and V(V) complexes of N-p-Tolylbenzohydroxamic acid. Int J Biol Macromol 2024; 268:131768. [PMID: 38663706 DOI: 10.1016/j.ijbiomac.2024.131768] [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: 09/21/2023] [Revised: 03/07/2024] [Accepted: 04/20/2024] [Indexed: 05/04/2024]
Abstract
Copper(L2Cu) and vanadium(L2VOCl) complexes of N-p-tolylbenzohydroxamic acid (LH) ligand have been investigated for DNA binding efficacy by multiple analytical, spectral, and computational techniques. The results revealed that complexes as groove binders as evidenced by UV absorption. Fluorescence studies including displacement assay using classical intercalator ethidium bromide as fluorescent probe also confirmed as groove binders. The viscometric analysis too supports the inferences as strong groove binders for both the complexes. Molecular docking too exposed DNA as a target to the complexes which precisely binds L2Cu, in the minor groove region while L2VOCl in major groove region. Molecular dynamic simulation performed on L2Cu complex revealing the interaction of complex with DNA within 20 ns time. The complex stacked into the nitrogen bases of oligonucleotides and the bonding features were intrinsically preserved for longer simulation times. In-vitro cytotoxicity study was undertaken employing MTT assay against the breast cancer cell line (MCF-7). Potential cytotoxic activities were observed for L2Cu and L2VOCl complexes with IC50 values of showing 71 % and 74 % of inhibition respectively.
Collapse
Affiliation(s)
- Mamta Tripathi
- School of Studies in Chemistry, Pt. Ravishankar Shukla University, Raipur 492010, Chhattisgarh, India.
| | - Yamini Thakur
- School of Studies in Chemistry, Pt. Ravishankar Shukla University, Raipur 492010, Chhattisgarh, India; Department of Chemistry, Govt. J. Yoganandam Chhattisgarh College, Raipur C.G. - 492001
| | - Rabbani Syed
- Department of Pharmaceutics, College of Pharmacy, King Saud University, P.O. Box 2457, Riyadh 11451, Saudi Arabia
| | - Ashish Kumar Asatkar
- Department of Chemistry, Satya Narayan Agrawal Govt. Arts and Commerce College, Kohka-Neora, Dist. Raipur, CG 493114, India
| | - Mohammad S Alqahtani
- Department of Pharmaceutics, College of Pharmacy, King Saud University, P.O. Box 2457, Riyadh 11451, Saudi Arabia
| | - Devashish Das
- Department of Chemical Engineering, Konkuk University, Seoul, South Korea
| | - Rainy Agrawal
- School of Studies in Chemistry, Pt. Ravishankar Shukla University, Raipur 492010, Chhattisgarh, India
| | - Bharati Verma
- School of Studies in Chemistry, Pt. Ravishankar Shukla University, Raipur 492010, Chhattisgarh, India
| | - Rama Pande
- School of Studies in Chemistry, Pt. Ravishankar Shukla University, Raipur 492010, Chhattisgarh, India
| |
Collapse
|
2
|
Afanasyev AY, Onufriev AV. Stretching of Long Double-Stranded DNA and RNA Described by the Same Approach. J Chem Theory Comput 2022; 18:3911-3920. [PMID: 35544776 DOI: 10.1021/acs.jctc.1c01221] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We propose an approach to help interpret polymer force-extension curves that exhibit plateau regimes. When coupled to a bead-spring dynamic model, the approach accurately reproduces a variety of experimental force-extension curves of long double-stranded DNA and RNA, including torsionally constrained and unconstrained DNA and negatively supercoiled DNA. A key feature of the model is a specific nonconvex energy function of the spring. We provide an algorithm to obtain the five required parameters of the model from experimental force-extension curves. The applicability of the approach to the force-extension curves of double-stranded (ds) DNA of variable GC content as well as to a DNA/RNA hybrid structure is explored and confirmed. We use the approach to explain counterintuitive sequence-dependent trends and make predictions. In the plateau region of the force-extension curves, our molecular dynamics simulations show that the polymer separates into a mix of weakly and strongly stretched states without forming macroscopically distinct phases. The distribution of these states is predicted to depend on the sequence.
Collapse
Affiliation(s)
- Alexander Y Afanasyev
- Department of Biomedical Engineering and Mechanics, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
| | - Alexey V Onufriev
- Departments of Computer Science and Physics, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
| |
Collapse
|
3
|
Sun LZ, Qian JL, Cai P, Hu HX, Xu X, Luo MB. Mg2+ effects on the single-stranded DNA conformations and nanopore translocation dynamics. POLYMER 2022. [DOI: 10.1016/j.polymer.2022.124895] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
|
4
|
Lim W, Randisi F, Doye JPK, Louis AA. The interplay of supercoiling and thymine dimers in DNA. Nucleic Acids Res 2022; 50:2480-2492. [PMID: 35188542 PMCID: PMC8934635 DOI: 10.1093/nar/gkac082] [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] [Received: 09/27/2021] [Revised: 01/21/2022] [Accepted: 02/04/2022] [Indexed: 11/16/2022] Open
Abstract
Thymine dimers are a major mutagenic photoproduct induced by UV radiation. While they have been the subject of extensive theoretical and experimental investigations, questions of how DNA supercoiling affects local defect properties, or, conversely, how the presence of such defects changes global supercoiled structure, are largely unexplored. Here, we introduce a model of thymine dimers in the oxDNA forcefield, parametrized by comparison to melting experiments and structural measurements of the thymine dimer induced bend angle. We performed extensive molecular dynamics simulations of double-stranded DNA as a function of external twist and force. Compared to undamaged DNA, the presence of a thymine dimer lowers the supercoiling densities at which plectonemes and bubbles occur. For biologically relevant supercoiling densities and forces, thymine dimers can preferentially segregate to the tips of the plectonemes, where they enhance the probability of a localized tip-bubble. This mechanism increases the probability of highly bent and denatured states at the thymine dimer site, which may facilitate repair enzyme binding. Thymine dimer-induced tip-bubbles also pin plectonemes, which may help repair enzymes to locate damage. We hypothesize that the interplay of supercoiling and local defects plays an important role for a wider set of DNA damage repair systems.
Collapse
Affiliation(s)
- Wilber Lim
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Parks Road, Oxford OX1 3PU, UK
| | - Ferdinando Randisi
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Parks Road, Oxford OX1 3PU, UK
- FabricNano Limited, 192 Drummond St, London NW1 3HP, UK
| | - Jonathan P K Doye
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford, OX1 3QZ, UK
| | - Ard A Louis
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Parks Road, Oxford OX1 3PU, UK
| |
Collapse
|
5
|
Gravina NM, Gumbart JC, Kim HD. Coarse-Grained Simulations of DNA Reveal Angular Dependence of Sticky-End Binding. J Phys Chem B 2021; 125:4016-4024. [PMID: 33870695 DOI: 10.1021/acs.jpcb.1c00432] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Annealing between sticky ends of DNA is an intermediate step in ligation. It can also be utilized to program specific binding sites for DNA tile and origami assembly. This reaction is generally understood as a bimolecular reaction dictated by the local concentration of the sticky ends. Its dependence on the relative orientation between the sticky ends, however, is less understood. Here we report on the interactions between DNA sticky ends using the coarse-grained oxDNA model; specifically, we consider how the orientational alignment of the double-stranded DNA (dsDNA) segments affects the time required for the sticky ends to bind, τb. We specify the orientation of the dsDNA segments with three parameters: θ, which measures the angle between the helical axes, and ϕ1 and ϕ2, which measure rotations of each strand around the helical axis. We find that the binding time depends strongly on both θ and ϕ2: ∼20-fold change with θ and 10-fold change with ϕ2. The binding time is the fastest when the helical axes of duplexes are pointing toward each other and the sticky ends protrude from the farthest two points. Our result is relevant for predicting hybridization efficiency of sticky ends that are rotationally restricted.
Collapse
Affiliation(s)
- Nicholas M Gravina
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - James C Gumbart
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Harold D Kim
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| |
Collapse
|
6
|
Alexiou TS, Mintis DG, Mavrantzas VG. Molecular Dynamics Simulation of the Diffusion Dynamics of Linear DNA Fragments in Dilute Solution with the Parmbsc1 Force Field and Comparison with Experimental Data and Theoretical Models. Macromolecules 2020. [DOI: 10.1021/acs.macromol.0c00222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Terpsichori S. Alexiou
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
| | - Dimitris G. Mintis
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
| | - Vlasis G. Mavrantzas
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
- Department of Mechanical and Process Engineering, Particle Technology Laboratory, ETH Zürich, CH-8092 Zürich, Switzerland
| |
Collapse
|
7
|
Das AK, Demerdash ON, Head-Gordon T. Improvements to the AMOEBA Force Field by Introducing Anisotropic Atomic Polarizability of the Water Molecule. J Chem Theory Comput 2018; 14:6722-6733. [DOI: 10.1021/acs.jctc.8b00978] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
|
8
|
Henrich O, Gutiérrez Fosado YA, Curk T, Ouldridge TE. Coarse-grained simulation of DNA using LAMMPS : An implementation of the oxDNA model and its applications. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2018; 41:57. [PMID: 29748779 DOI: 10.1140/epje/i2018-11669-8] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/20/2018] [Accepted: 04/19/2018] [Indexed: 05/25/2023]
Abstract
During the last decade coarse-grained nucleotide models have emerged that allow us to study DNA and RNA on unprecedented time and length scales. Among them is oxDNA, a coarse-grained, sequence-specific model that captures the hybridisation transition of DNA and many structural properties of single- and double-stranded DNA. oxDNA was previously only available as standalone software, but has now been implemented into the popular LAMMPS molecular dynamics code. This article describes the new implementation and analyses its parallel performance. Practical applications are presented that focus on single-stranded DNA, an area of research which has been so far under-investigated. The LAMMPS implementation of oxDNA lowers the entry barrier for using the oxDNA model significantly, facilitates future code development and interfacing with existing LAMMPS functionality as well as other coarse-grained and atomistic DNA models.
Collapse
Affiliation(s)
- Oliver Henrich
- Department of Physics, SUPA, University of Strathclyde, G4 0NG, Glasgow, Scotland, UK.
| | | | - Tine Curk
- CAS Key Laboratory of Soft Matter Physics, Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, 100190, Beijing, China
- Department of Chemistry, University of Cambridge, CB2 1EW, Cambridge, UK
| | - Thomas E Ouldridge
- Department of Bioengineering, Imperial College London, SW7 2AZ, London, UK
- Centre of Synthetic Biology, Imperial College London, SW7 2AZ, London, UK
| |
Collapse
|
9
|
|
10
|
Zavadlav J, Podgornik R, Praprotnik M. Order and interactions in DNA arrays: Multiscale molecular dynamics simulation. Sci Rep 2017; 7:4775. [PMID: 28684875 PMCID: PMC5500594 DOI: 10.1038/s41598-017-05109-2] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2017] [Accepted: 06/06/2017] [Indexed: 11/21/2022] Open
Abstract
While densely packed DNA arrays are known to exhibit hexagonal and orthorhombic local packings, the detailed mechanism governing the associated phase transition remains rather elusive. Furthermore, at high densities the atomistic resolution is paramount to properly account for fine details, encompassing the DNA molecular order, the contingent ordering of counterions and the induced molecular ordering of the bathing solvent, bringing together electrostatic, steric, thermal and direct hydrogen-bonding interactions, resulting in the observed osmotic equation of state. We perform a multiscale simulation of dense DNA arrays by enclosing a set of 16 atomistically resolved DNA molecules within a semi-permeable membrane, allowing the passage of water and salt ions, and thus mimicking the behavior of DNA arrays subjected to external osmotic stress in a bathing solution of monovalent salt and multivalent counterions. By varying the DNA density, local packing symmetry, and counterion type, we obtain osmotic equation of state together with the hexagonal-orthorhombic phase transition, and full structural characterization of the DNA subphase in terms of its positional and angular orientational fluctuations, counterion distributions, and the solvent local dielectric response profile with its order parameters that allow us to identify the hydration force as the primary interaction mechanism at high DNA densities.
Collapse
Affiliation(s)
- Julija Zavadlav
- Department of Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, SI-1001, Ljubljana, Slovenia.,Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, SI-1000, Ljubljana, Slovenia.,Chair of Computational Science, ETH Zurich, Clausiusstrasse 33, CH-8092, Zurich, Switzerland
| | - Rudolf Podgornik
- Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, SI-1000, Ljubljana, Slovenia. .,Theoretical Physics Department, J. Stefan Institute, Jamova c. 39, SI-1000, Ljubljana, Slovenia.
| | - Matej Praprotnik
- Department of Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, SI-1001, Ljubljana, Slovenia. .,Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, SI-1000, Ljubljana, Slovenia.
| |
Collapse
|
11
|
Hilder TA, Hodgkiss JM. Molecular Mechanism of Binding between 17β-Estradiol and DNA. Comput Struct Biotechnol J 2016; 15:91-97. [PMID: 28066533 PMCID: PMC5196241 DOI: 10.1016/j.csbj.2016.12.001] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2016] [Revised: 11/30/2016] [Accepted: 12/05/2016] [Indexed: 02/05/2023] Open
Abstract
Although 17β-estradiol (E2) is a natural molecule involved in the endocrine system, its widespread use in various applications has resulted in its accumulation in the environment and its classification as an endocrine-disrupting molecule. These molecules can interfere with the hormonal system, and have been linked to various adverse effects such as the proliferation of breast cancer. It has been proposed that E2 could contribute to breast cancer by the induction of DNA damage. Mass spectrometry has demonstrated that E2 can bind to DNA but the mechanism by which E2 interacts with DNA has yet to be elucidated. Using all-atom molecular dynamics simulations, we demonstrate that E2 intercalates (inserts between two successive DNA base pairs) in DNA at the location specific to estrogen receptor binding, known as the estrogen response element (ERE), and to other random sequences of DNA. Our results suggest that excess E2 has the potential to disrupt processes in the body which rely on binding to DNA, such as the binding of the estrogen receptor to the ERE and the activity of enzymes that bind DNA, and could lead to DNA damage.
Collapse
Affiliation(s)
- Tamsyn A. Hilder
- School of Chemical and Physical Sciences, Victoria University of Wellington, Wellington 6040, New Zealand
- Computational Biophysics Group, Research School of Biology, Canberra, ACT 0200, Australia
- Correspondence to: T.A. Hilder, School of Chemical and Physical Sciences, Victoria University of Wellington, Wellington 6040, New Zealand.Computational Biophysics GroupResearch School of BiologyCanberraACT0200Australia
| | - Justin M. Hodgkiss
- School of Chemical and Physical Sciences, Victoria University of Wellington, Wellington 6040, New Zealand
- The MacDiarmid Institute of Advanced Materials and Nanotechnology, New Zealand
- Correspondence to: J.M. Hodgkiss, The MacDiarmid Institute of Advanced Materials and Nanotechnology, New Zealand.The MacDiarmid Institute of Advanced Materials and NanotechnologyNew Zealand
| |
Collapse
|
12
|
Lyubartsev AP, Naômé A, Vercauteren DP, Laaksonen A. Systematic hierarchical coarse-graining with the inverse Monte Carlo method. J Chem Phys 2016; 143:243120. [PMID: 26723605 DOI: 10.1063/1.4934095] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
We outline our coarse-graining strategy for linking micro- and mesoscales of soft matter and biological systems. The method is based on effective pairwise interaction potentials obtained in detailed ab initio or classical atomistic Molecular Dynamics (MD) simulations, which can be used in simulations at less accurate level after scaling up the size. The effective potentials are obtained by applying the inverse Monte Carlo (IMC) method [A. P. Lyubartsev and A. Laaksonen, Phys. Rev. E 52(4), 3730-3737 (1995)] on a chosen subset of degrees of freedom described in terms of radial distribution functions. An in-house software package MagiC is developed to obtain the effective potentials for arbitrary molecular systems. In this work we compute effective potentials to model DNA-protein interactions (bacterial LiaR regulator bound to a 26 base pairs DNA fragment) at physiological salt concentration at a coarse-grained (CG) level. Normally the IMC CG pair-potentials are used directly as look-up tables but here we have fitted them to five Gaussians and a repulsive wall. Results show stable association between DNA and the model protein as well as similar position fluctuation profile.
Collapse
Affiliation(s)
- Alexander P Lyubartsev
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, S 106 91 Stockholm, Sweden
| | - Aymeric Naômé
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, S 106 91 Stockholm, Sweden
| | | | - Aatto Laaksonen
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, S 106 91 Stockholm, Sweden
| |
Collapse
|
13
|
Nick TA, de Oliveira TE, Pilat DW, Spenkuch F, Butt HJ, Helm M, Netz PA, Berger R. Stability of a Split Streptomycin Binding Aptamer. J Phys Chem B 2016; 120:6479-89. [PMID: 27281393 DOI: 10.1021/acs.jpcb.6b02440] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Thomas A Nick
- Max Planck Institute for Polymer Research , 55128 Mainz, Germany
| | - Tiago E de Oliveira
- Instituto de Química, Universidade Federal do Rio Grande do Sul , Avenida Bento Gonçalves, 9500, 91501-970 Porto Alegre-RS, Brazil
| | - Dominik W Pilat
- Max Planck Institute for Polymer Research , 55128 Mainz, Germany
| | - Felix Spenkuch
- Johannes Gutenberg Universität Mainz , Institute of Pharmacy and Biochemistry, 55128 Mainz, Germany
| | - Hans-Jürgen Butt
- Max Planck Institute for Polymer Research , 55128 Mainz, Germany
| | - Mark Helm
- Johannes Gutenberg Universität Mainz , Institute of Pharmacy and Biochemistry, 55128 Mainz, Germany
| | - Paulo A Netz
- Instituto de Química, Universidade Federal do Rio Grande do Sul , Avenida Bento Gonçalves, 9500, 91501-970 Porto Alegre-RS, Brazil
| | - Rüdiger Berger
- Max Planck Institute for Polymer Research , 55128 Mainz, Germany
| |
Collapse
|
14
|
Snodin BEK, Randisi F, Mosayebi M, Šulc P, Schreck JS, Romano F, Ouldridge TE, Tsukanov R, Nir E, Louis AA, Doye JPK. Introducing improved structural properties and salt dependence into a coarse-grained model of DNA. J Chem Phys 2016; 142:234901. [PMID: 26093573 DOI: 10.1063/1.4921957] [Citation(s) in RCA: 217] [Impact Index Per Article: 27.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
We introduce an extended version of oxDNA, a coarse-grained model of deoxyribonucleic acid (DNA) designed to capture the thermodynamic, structural, and mechanical properties of single- and double-stranded DNA. By including explicit major and minor grooves and by slightly modifying the coaxial stacking and backbone-backbone interactions, we improve the ability of the model to treat large (kilobase-pair) structures, such as DNA origami, which are sensitive to these geometric features. Further, we extend the model, which was previously parameterised to just one salt concentration ([Na(+)] = 0.5M), so that it can be used for a range of salt concentrations including those corresponding to physiological conditions. Finally, we use new experimental data to parameterise the oxDNA potential so that consecutive adenine bases stack with a different strength to consecutive thymine bases, a feature which allows a more accurate treatment of systems where the flexibility of single-stranded regions is important. We illustrate the new possibilities opened up by the updated model, oxDNA2, by presenting results from simulations of the structure of large DNA objects and by using the model to investigate some salt-dependent properties of DNA.
Collapse
Affiliation(s)
- Benedict E K Snodin
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Ferdinando Randisi
- Life Sciences Interface Doctoral Training Center, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Majid Mosayebi
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Petr Šulc
- Center for Studies in Physics and Biology, The Rockefeller University, 1230 York Avenue, New York, New York 10065, USA
| | - John S Schreck
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Flavio Romano
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Thomas E Ouldridge
- Department of Mathematics, Imperial College, 180 Queen's Gate, London SW7 2AZ, United Kingdom
| | - Roman Tsukanov
- Department of Chemistry and the Ilse Katz Institute for Nanoscale Science and Technology, Ben-Gurion University of the Negev, Beer Sheva, Israel
| | - Eyal Nir
- Department of Chemistry and the Ilse Katz Institute for Nanoscale Science and Technology, Ben-Gurion University of the Negev, Beer Sheva, Israel
| | - Ard A Louis
- Rudolf Peierls Centre for Theoretical Physics, 1 Keble Road, Oxford OX1 3NP, United Kingdom
| | - Jonathan P K Doye
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
15
|
Kocak A, Yilmaz H, Faiz O, Andac O. Experimental and theoretical studies on Cu(II) complex of N,N′-disalicylidene-2,3-diaminopyridine ligand reveal indirect evidence for DNA intercalation. Polyhedron 2016. [DOI: 10.1016/j.poly.2015.11.037] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
16
|
Sedeh RS, Pan K, Adendorff MR, Hallatschek O, Bathe KJ, Bathe M. Computing Nonequilibrium Conformational Dynamics of Structured Nucleic Acid Assemblies. J Chem Theory Comput 2015; 12:261-73. [PMID: 26636351 DOI: 10.1021/acs.jctc.5b00965] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Synthetic nucleic acids can be programmed to form precise three-dimensional structures on the nanometer-scale. These thermodynamically stable complexes can serve as structural scaffolds to spatially organize functional molecules including multiple enzymes, chromophores, and force-sensing elements with internal dynamics that include substrate reaction-diffusion, excitonic energy transfer, and force-displacement response that often depend critically on both the local and global conformational dynamics of the nucleic acid assembly. However, high molecular weight assemblies exhibit long time-scale and large length-scale motions that cannot easily be sampled using all-atom computational procedures such as molecular dynamics. As an alternative, here we present a computational framework to compute the overdamped conformational dynamics of structured nucleic acid assemblies and apply it to a DNA-based tweezer, a nine-layer DNA origami ring, and a pointer-shaped DNA origami object, which consist of 204, 3,600, and over 7,000 basepairs, respectively. The framework employs a mechanical finite element model for the DNA nanostructure combined with an implicit solvent model to either simulate the Brownian dynamics of the assembly or alternatively compute its Brownian modes. Computational results are compared with an all-atom molecular dynamics simulation of the DNA-based tweezer. Several hundred microseconds of Brownian dynamics are simulated for the nine-layer ring origami object to reveal its long time-scale conformational dynamics, and the first ten Brownian modes of the pointer-shaped structure are predicted.
Collapse
Affiliation(s)
| | | | | | - Oskar Hallatschek
- Department of Physics, University of California, Berkeley , Berkeley, California 94720, United States
| | | | | |
Collapse
|
17
|
Zgarbová M, Šponer J, Otyepka M, Cheatham TE, Galindo-Murillo R, Jurečka P. Refinement of the Sugar-Phosphate Backbone Torsion Beta for AMBER Force Fields Improves the Description of Z- and B-DNA. J Chem Theory Comput 2015; 11:5723-36. [PMID: 26588601 DOI: 10.1021/acs.jctc.5b00716] [Citation(s) in RCA: 347] [Impact Index Per Article: 38.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Z-DNA duplexes are a particularly complicated test case for current force fields. We performed a set of explicit solvent molecular dynamics (MD) simulations with various AMBER force field parametrizations including our recent refinements of the ε/ζ and glycosidic torsions. None of these force fields described the ZI/ZII and other backbone substates correctly, and all of them underpredicted the population of the important ZI substate. We show that this underprediction can be attributed to an inaccurate potential for the sugar-phosphate backbone torsion angle β. We suggest a refinement of this potential, β(OL1), which was derived using our recently introduced methodology that includes conformation-dependent solvation effects. The new potential significantly increases the stability of the dominant ZI backbone substate and improves the overall description of the Z-DNA backbone. It also has a positive (albeit small) impact on another important DNA form, the antiparallel guanine quadruplex (G-DNA), and improves the description of the canonical B-DNA backbone by increasing the population of BII backbone substates, providing a better agreement with experiment. We recommend using β(OL1) in combination with our previously introduced corrections, εζ(OL1) and χ(OL4), (the combination being named OL15) as a possible alternative to the current β torsion potential for more accurate modeling of nucleic acids.
Collapse
Affiliation(s)
- Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic.,Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, University of Utah , 30 South 2000 East, Skaggs 105, Salt Lake City, Utah 84112, United States
| | - Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, University of Utah , 30 South 2000 East, Skaggs 105, Salt Lake City, Utah 84112, United States
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic
| |
Collapse
|
18
|
Iacovelli F, Falconi M. Decoding the conformation-linked functional properties of nucleic acids by the use of computational tools. FEBS J 2015; 282:3298-310. [DOI: 10.1111/febs.13315] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2015] [Revised: 04/16/2015] [Accepted: 04/29/2015] [Indexed: 12/25/2022]
Affiliation(s)
| | - Mattia Falconi
- Department of Biology; University of Rome “Tor Vergata”; Italy
| |
Collapse
|
19
|
Baker CM. Polarizable force fields for molecular dynamics simulations of biomolecules. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2015. [DOI: 10.1002/wcms.1215] [Citation(s) in RCA: 100] [Impact Index Per Article: 11.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
20
|
Argudo D, Purohit PK. Equilibrium and kinetics of DNA overstretching modeled with a quartic energy landscape. Biophys J 2014; 107:2151-63. [PMID: 25418100 DOI: 10.1016/j.bpj.2014.09.014] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2014] [Revised: 08/22/2014] [Accepted: 09/16/2014] [Indexed: 11/18/2022] Open
Abstract
It is well known that the dsDNA molecule undergoes a phase transition from B-DNA into an overstretched state at high forces. For some time, the structure of the overstretched state remained unknown and highly debated, but recent advances in experimental techniques have presented evidence of more than one possible phase (or even a mixed phase) depending on ionic conditions, temperature, and basepair sequence. Here, we present a theoretical model to study the overstretching transition with the possibility that the overstretched state is a mixture of two phases: a structure with portions of inner strand separation (melted or M-DNA), and an extended phase that retains the basepair structure (S-DNA). We model the double-stranded DNA as a chain composed of n segments of length l, where the transition is studied by means of a Landau quartic potential with statistical fluctuations. The length l is a measure of cooperativity of the transition and is key to characterizing the overstretched phase. By analyzing the different values of l corresponding to a wide spectrum of experiments, we find that for a range of temperatures and ionic conditions, the overstretched form is likely to be a mix of M-DNA and S-DNA. For a transition close to a pure S-DNA state, where the change in extension is close to 1.7 times the original B-DNA length, we find l ? 25 basepairs regardless of temperature and ionic concentration. Our model is fully analytical, yet it accurately reproduces the force-extension curves, as well as the transient kinetic behavior, seen in DNA overstretching experiments.
Collapse
Affiliation(s)
- David Argudo
- Department of Mechanical Engineering and Applied Mechanics, University of Pennsylvania, Philadelphia, Pennsylvania
| | - Prashant K Purohit
- Department of Mechanical Engineering and Applied Mechanics, University of Pennsylvania, Philadelphia, Pennsylvania.
| |
Collapse
|
21
|
Li NK, Kim HS, Nash JA, Lim M, Yingling YG. Progress in molecular modelling of DNA materials. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.913792] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
22
|
Guardiani C, Cencini M, Cecconi F. Coarse-grained modeling of protein unspecifically bound to DNA. Phys Biol 2014; 11:026003. [PMID: 24685517 DOI: 10.1088/1478-3975/11/2/026003] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
There is now a certain consensus that transcription factors (TFs) reach their target sites, where they regulate gene transcription, via a mechanism dubbed facilitated diffusion (FD). In FD, the TF cycles between events of 3D diffusion in solution (jumps), 1D diffusion along DNA (sliding), and small jumps (hopping), achieving association rates higher than for 3D diffusion alone. We investigate the FD phenomenology through molecular dynamics simulations in the framework of coarse-grained modeling. We show that, despite the crude approximations, the model generates, upon varying the equilibrium distance of the DNA-TF interaction, a phenomenology matching a number of experimental and numerical results obtained with more refined models. In particular, focusing on the kinematics of the process, we characterize the geometrical properties of TF trajectories during sliding. We find that sliding occurs via helical paths around the DNA helix, leading to a coupling of translation along the DNA axis with rotation around it. The 1D diffusion constant measured in simulations is found to be interwoven with the geometrical properties of sliding and we develop a simple argument that can be used to quantitatively reproduce the measured values.
Collapse
Affiliation(s)
- Carlo Guardiani
- Dipartimento di Fisica, University 'Sapienza', Piazzale A. Moro 2, 00185 Rome, Italy
| | | | | |
Collapse
|
23
|
DNA-binding of nickel(II), copper(II) and zinc(II) complexes: Structure–affinity relationships. Coord Chem Rev 2013. [DOI: 10.1016/j.ccr.2013.02.023] [Citation(s) in RCA: 203] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
24
|
Abstract
The physiological properties of biological soft matter are the product of collective interactions, which span many time and length scales. Recent computational modeling efforts have helped illuminate experiments that characterize the ways in which proteins modulate membrane physics. Linking these models across time and length scales in a multiscale model explains how atomistic information propagates to larger scales. This paper reviews continuum modeling and coarse-grained molecular dynamics methods, which connect atomistic simulations and single-molecule experiments with the observed microscopic or mesoscale properties of soft-matter systems essential to our understanding of cells, particularly those involved in sculpting and remodeling cell membranes.
Collapse
Affiliation(s)
- Ryan Bradley
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Ravi Radhakrishnan
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA 19104, USA
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA 19104, USA
- Author to whom correspondence should be addressed; ; Tel.: +1-215-898-0487; Fax: +1-215-573-2071
| |
Collapse
|
25
|
Metal complex–DNA binding: Insights from molecular dynamics and DFT/MM calculations. J Inorg Biochem 2013; 124:63-9. [DOI: 10.1016/j.jinorgbio.2013.03.010] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2012] [Revised: 03/19/2013] [Accepted: 03/20/2013] [Indexed: 11/19/2022]
|
26
|
Hospital A, Faustino I, Collepardo-Guevara R, González C, Gelpí JL, Orozco M. NAFlex: a web server for the study of nucleic acid flexibility. Nucleic Acids Res 2013; 41:W47-55. [PMID: 23685436 PMCID: PMC3692121 DOI: 10.1093/nar/gkt378] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
We present NAFlex, a new web tool to study the flexibility of nucleic acids, either isolated or bound to other molecules. The server allows the user to incorporate structures from protein data banks, completing gaps and removing structural inconsistencies. It is also possible to define canonical (average or sequence-adapted) nucleic acid structures using a variety of predefined internal libraries, as well to create specific nucleic acid conformations from the sequence. The server offers a variety of methods to explore nucleic acid flexibility, such as a colorless wormlike-chain model, a base-pair resolution mesoscopic model and atomistic molecular dynamics simulations with a wide variety of protocols and force fields. The trajectories obtained by simulations, or imported externally, can be visualized and analyzed using a large number of tools, including standard Cartesian analysis, essential dynamics, helical analysis, local and global stiffness, energy decomposition, principal components and in silico NMR spectra. The server is accessible free of charge from the mmb.irbbarcelona.org/NAFlex webpage.
Collapse
Affiliation(s)
- Adam Hospital
- Institute for Research in Biomedicine, Molecular Modelling and Bioinformatics Department, 08028, Barcelona, Spain
| | | | | | | | | | | |
Collapse
|
27
|
Abstract
Small DNA circles can occur in Nature, for example as protein-constrained loops, and can be synthesized by a number of methods. Such small circles provide tractable systems for the study of the structure, thermodynamics and molecular dynamics of closed-circular DNA. In the present article, we review the occurrence and synthesis of small DNA circles, and examine their utility in studying the properties of DNA and DNA–protein interactions. In particular, we highlight the analysis of small circles using atomistic simulations.
Collapse
|
28
|
Potoyan DA, Savelyev A, Papoian GA. Recent successes in coarse-grained modeling of DNA. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2012. [DOI: 10.1002/wcms.1114] [Citation(s) in RCA: 66] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
|
29
|
Abstract
This article provides a retrospective on the ABC initiative in the area of all-atom molecular dynamics (MD) simulations including explicit solvent on all tetranucleotide steps of duplex B-form DNA duplex, ca. 2012. The ABC consortium has completed two phases of simulations, the most current being a set of 50-100 trajectories based on the AMBER ff99 force field together with the parmbsc0 modification. Some general perspectives on the field of MD on DNA and sequence effects on DNA structure are provided, followed by an overview our MD results, including a detailed comparison of the ff99/parmbsc0 results with crystal and NMR structures available for d(CGCGAATTCGCG). Some projects inspired by or related to the ABC initiative and database are also reviewed, including methods for the trajectory analyses, informatics of dealing with the large database of results, compressions of trajectories for efficacy of distribution, DNA solvation by water and ions, parameterization of coarse-grained models with applications and gene finding and genome annotation.
Collapse
Affiliation(s)
- David L Beveridge
- Department of Chemistry and Molecular Biophysics Program, Wesleyan University Middletown, CT 06459, USA.
| | | | | |
Collapse
|
30
|
Abstract
It has been known for decades that DNA is extremely flexible and polymorphic, but our knowledge of its accessible conformational space remains limited. Structural data, primarily from X-ray diffraction studies, is sparse in comparison to the manifold configurations possible, and direct experimental examinations of DNA's flexibility still suffer from many limitations. In the face of these shortcomings, molecular dynamics (MD) is now an essential tool in the study of DNA. It affords detailed structural and dynamical insights, which explains its recent transition from a small number of highly specialized laboratories to a large variety of groups dealing with challenging biological problems. MD is now making an irreversible journey to the mainstream of research in biology, with the attendant opportunities and challenges. But given the speed with which MD studies of DNA have spread, the roots remain somewhat shallow: in many cases, there is a lack of deep knowledge about the foundations, strengths, and limits of the technique. In this Account, we discuss how MD has become the most important source of structural and flexibility data on DNA, focusing on advances since 2007 of atomistic MD in the description of DNA under near-physiological conditions and highlighting the possibilities and shortcomings of the technique. The evolution in the field over the past four years is a prelude to the ongoing revolution. The technique has gained in robustness and predictive power, which when coupled with the spectacular improvements in software and hardware has enabled the tackling of systems of increasing complexity. Simulation times of microseconds have now been achieved, with even longer times when specialized hardware is used. As a result, we have seen the first real-time simulation of large conformational transitions, including folding and unfolding of short DNA duplexes. Noteworthy advances have also been made in the study of DNA-ligand interactions, and we predict that a global thermodynamic and kinetic picture of the binding landscape of DNA will become available in a few years. MD will become a crucial tool in areas such as biomolecular engineering and synthetic biology. MD has also been shown to be an excellent source of parameters for mesoscopic models of DNA flexibility. Such models can be refined through atomistic MD simulations on small duplexes and then applied to the study of entire chromosomes. Recent evidence suggests that MD-derived elastic models can successfully predict the position of regulatory regions in DNA and can help advance our understanding of nucleosome positioning and chromatin plasticity. If these results are confirmed, MD simulations can become the ultimate tool to decipher a physical code that can contribute to gene regulation. We are entering the golden age of MD simulations of DNA. Undoubtedly, the expectations are high, but the challenges are also enormous. These include the need for more accurate potential energy functionals and for longer and more complex simulations in more realistic systems. The joint research effort of several groups will be crucial for adapting the technique to the requirements of the coming decade.
Collapse
Affiliation(s)
- Alberto Pérez
- Joint IRB-BSC Program in Computational Biology, Institute of Research in Biomedicine Barcelona, Baldiri i Reixac 10, Barcelona 08028, Spain
| | - F. Javier Luque
- Department de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona 08028, Spain
| | - Modesto Orozco
- Joint IRB-BSC Program in Computational Biology, Institute of Research in Biomedicine Barcelona, Baldiri i Reixac 10, Barcelona 08028, Spain
- Departament de Bioquímica, Universitat de Barcelona, Avgda Diagonal 647, Barcelona 08028, Spain, and Instituto Nacional de Bioinformàtica, Parc Científic de Barcelona, Baldiri i Reixac 10, Barcelona 08028, Spain
| |
Collapse
|