1
|
Zheng CC, Chen YL, Dong HL, Zhang XH, Tan ZJ. Effect of ethanol on the elasticities of double-stranded RNA and DNA revealed by magnetic tweezers and simulations. J Chem Phys 2024; 161:075101. [PMID: 39145565 DOI: 10.1063/5.0211869] [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: 04/02/2024] [Accepted: 07/29/2024] [Indexed: 08/16/2024] Open
Abstract
The elasticities of double-stranded (ds) DNA and RNA, which are critical to their biological functions and applications in materials science, can be significantly modulated by solution conditions such as ions and temperature. However, there is still a lack of a comprehensive understanding of the role of solvents in the elasticities of dsRNA and dsDNA in a comparative way. In this work, we explored the effect of ethanol solvent on the elasticities of dsRNA and dsDNA by magnetic tweezers and all-atom molecular dynamics simulations. We found that the bending persistence lengths and contour lengths of dsRNA and dsDNA decrease monotonically with the increase in ethanol concentration. Furthermore, the addition of ethanol weakens the positive twist-stretch coupling of dsRNA, while promotes the negative twist-stretch coupling of dsDNA. Counter-intuitively, the lower dielectric environment of ethanol causes a significant re-distribution of counterions and enhanced ion neutralization, which overwhelms the enhanced repulsion along dsRNA/dsDNA, ultimately leading to the softening in bending for dsRNA and dsDNA. Moreover, for dsRNA, ethanol causes slight ion-clamping across the major groove, which weakens the major groove-mediated twist-stretch coupling, while for dsDNA, ethanol promotes the stretch-radius correlation due to enhanced ion binding and consequently enhances the helical radius-mediated twist-stretch coupling.
Collapse
Affiliation(s)
- Chen-Chen Zheng
- School of Physics and Technology, College of Life Sciences, Wuhan University, Wuhan 430072, China
| | - Yun-Long Chen
- School of Physics and Technology, College of Life Sciences, Wuhan University, Wuhan 430072, China
| | - Hai-Long Dong
- School of Physics and Technology, College of Life Sciences, Wuhan University, Wuhan 430072, China
| | - Xing-Hua Zhang
- School of Physics and Technology, College of Life Sciences, Wuhan University, Wuhan 430072, China
| | - Zhi-Jie Tan
- School of Physics and Technology, College of Life Sciences, Wuhan University, Wuhan 430072, China
| |
Collapse
|
2
|
Mamun Y, Tse-Dinh YC, Chapagain P. Insights into the DNA and RNA Interactions of Human Topoisomerase III Beta Using Molecular Dynamics Simulations. J Chem Inf Model 2024; 64:6062-6071. [PMID: 39024468 PMCID: PMC11323020 DOI: 10.1021/acs.jcim.4c00472] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2024] [Revised: 06/27/2024] [Accepted: 06/28/2024] [Indexed: 07/20/2024]
Abstract
Human topoisomerase III beta (hTOP3B) is the only topoisomerase in the human cell that can act on both DNA and RNA substrates. Recent findings have emphasized the physiological importance of hTOP3B and consolidated it as a valuable drug target for antiviral and anticancer therapeutics. Although type IA topoisomerases of different organisms have been studied over the years, the step-by-step interaction of hTOP3B and nucleic acid substrates is still not well understood. Due to the lack of hTOP3B-RNA structures as well as DNA/RNA covalent complexes, computational investigations have been limited. In our study, we utilized molecular dynamics (MD) simulations to study the interactions between hTOP3B and nucleic acids to get a closer look into the residues that play a role in binding DNA or RNA and facilitate catalysis, along with the differences and similarities when hTOP3B interacts with DNA compared to RNA. For this, we generated multiple models of hTOP3B complexed with DNA and RNA sequences using the hTOP3B crystal structure and 8-mer single-stranded DNA and RNA sequences. These models include both covalent and noncovalent complexes, which are then subjected to MD simulations and analyzed. Our findings highlight the complexes' stability, sequence preference, and interactions of the binding pocket residues with different nucleotides. Our work demonstrates that hTOP3B forms stable complexes with both DNA and RNA and provides a better understanding of the enzyme's interaction with different nucleic acid substrate sequences.
Collapse
Affiliation(s)
- Yasir Mamun
- Department
of Chemistry and Biochemistry, Florida International
University, Miami, Florida 33199, United States
| | - Yuk-Ching Tse-Dinh
- Department
of Chemistry and Biochemistry, Florida International
University, Miami, Florida 33199, United States
- Biomolecular
Sciences Institute, Florida International
University, Miami, Florida 33199, United States
| | - Prem Chapagain
- Department
of Physics, Florida International University, Miami, Florida 33199, United States
- Biomolecular
Sciences Institute, Florida International
University, Miami, Florida 33199, United States
| |
Collapse
|
3
|
Mak CH. Hydration Waters Make Up for the Missing Third Hydrogen Bond in the A·T Base Pair. ACS PHYSICAL CHEMISTRY AU 2024; 4:180-190. [PMID: 38560756 PMCID: PMC10979491 DOI: 10.1021/acsphyschemau.3c00058] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Revised: 01/05/2024] [Accepted: 01/08/2024] [Indexed: 04/04/2024]
Abstract
Base pairing complementarity is central to DNA function. G·C and A·T pair specificity is thought to originate from the different number of hydrogen bonds the pairs make. Quantifying how many hydrogen bonds exist can be difficult because water molecules in the surrounding can make up for or disrupt direct hydrogen bonds, and the hydration structures around A·T and G·C pairs on duplex DNA are distinct. Large-scale computer simulations have been used here to create a detailed map for the hydration structure on A·T and G·C base pairs in water. The contributions of specific hydration waters to the free energy of each of the hydrogen bonds in the A·T and G·C pairs were computed. Using the equilibrium fractions of hydrated versus unhydrated states from the hydration profiles, the impact of specific bound waters on each hydrogen bond can be uniquely quantified using a thermodynamic construction. The findings suggest that hydration water in the minor groove of an A·T pair can provide up to about 2 kcal/mol of free energy advantage, effectively making up for the missing third hydrogen bond in the A·T pair compared to G·C, rendering the intrinsic thermodynamic stability of the A·T pair almost synonymous with G·C.
Collapse
Affiliation(s)
- Chi H. Mak
- Departments of Chemistry
and Quantitative and Computational Biology, and Center of Applied
Mathematical Sciences, University of Southern
California, Los Angeles, California 90089, United States
| |
Collapse
|
4
|
Kılıç M, Diamantis P, Johnson SK, Toth O, Rothlisberger U. Redox-Based Defect Detection in Packed DNA: Insights from Hybrid Quantum Mechanical/Molecular Mechanics Molecular Dynamics Simulations. J Chem Theory Comput 2023; 19:8434-8445. [PMID: 37963372 PMCID: PMC10687876 DOI: 10.1021/acs.jctc.3c01013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Revised: 10/30/2023] [Accepted: 10/31/2023] [Indexed: 11/16/2023]
Abstract
The impact of an 8-oxoguanine (8oxoG) defect on the redox properties of DNA within the nucleosome core particle (NCP) was investigated employing hybrid quantum mechanical/molecular mechanics (QM/MM) molecular dynamics simulations of native and 8oxoG-containing NCP systems with an explicit representation of a biologically relevant environment. Two distinct NCP positions with varying solvent accessibility were considered for 8oxoG insertion. In both cases, it is found that the presence of 8oxoG drastically decreases the redox free energy of oxidation by roughly 1 eV, which is very similar to what was recently reported for free native and 8oxoG-containing DNA. In contrast, the effect of 8oxoG on the reorganization free energy is even smaller for packed DNA (decrease of 0.13 and 0.01 eV for defect-free and defect-containing systems, respectively) compared to the one for free DNA (0.25 eV), consistent with the increased rigidity of the NCP as compared to free DNA. Furthermore, the presence of an 8oxoG defect does not yield any significant changes in the packed DNA structure. Such a conclusion favors the idea that in the case of chromatin, defect-induced changes in DNA redox chemistry can also be exploited to detect damaged bases via DNA-mediated hole transfer.
Collapse
Affiliation(s)
| | | | - Sophia K. Johnson
- Laboratory of Computational Chemistry
and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale
de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Oliver Toth
- Laboratory of Computational Chemistry
and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale
de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry
and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale
de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
5
|
Templeton C, Hamilton I, Russell R, Elber R. Impact of Ion-Mixing Entropy on Orientational Preferences of DNA Helices: FRET Measurements and Computer Simulations. J Phys Chem B 2023; 127:8796-8808. [PMID: 37815452 PMCID: PMC11341850 DOI: 10.1021/acs.jpcb.3c04354] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/11/2023]
Abstract
Biological processes require DNA and RNA helices to pack together in specific interhelical orientations. While electrostatic repulsion between backbone charges is expected to be maximized when helices are in parallel alignment, such orientations are commonplace in nature. To better understand how the repulsion is overcome, we used experimental and computational approaches to investigate how the orientational preferences of DNA helices depend on the concentration and valence of mobile cations. We used Förster resonance energy transfer (FRET) to probe the relative orientations of two 24-bp helices held together via a freely rotating PEG linker. At low cation concentrations, the helices preferred more "cross"-like orientations over those closer to parallel, and this preference was reduced with increasing salt concentrations. The results were in good quantitative agreement with Poisson-Boltzmann (PB) calculations for monovalent salt (Na+). However, PB underestimated the ability of mixtures of monovalent and divalent ions (Mg2+) to reduce the conformational preference. As a complementary approach, we performed all-atom molecular dynamics (MD) simulations and found better agreement with the experimental results. While MD and PB predict similar electrostatic forces, MD predicts a greater accumulation of Mg2+ in the ion atmosphere surrounding the DNA. Mg2+ occupancy is predicted to be greater in conformations close to the parallel orientation than in conformations close to the crossed orientation, enabling a greater release of Na+ ions and providing an entropic gain (one bound ion for two released). MD predicts an entropy gain larger than that of PB because of the increased Mg2+ occupancy. The entropy changes have a negligible effect at low Mg2+ concentrations because the free energies are dominated by electrostatic repulsion. However, as the Mg2+ concentration increases, charge screening is more effective and the mixing entropy produces readily detectable changes in packing preferences. Our results underline the importance of mixing entropy of counterions in nucleic acid interactions and provide a new understanding on the impact of a mixed ion atmosphere on the packing of DNA helices.
Collapse
Affiliation(s)
- Clark Templeton
- Department of Chemical Engineering, University of Texas at Austin, Austin, TX 78712, USA
- Department of Physics, FU Berlin, 14195 Berlin, Germany
| | - Ian Hamilton
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712, USA
| | - Rick Russell
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712, USA
| | - Ron Elber
- Institute for Computational Engineering and Science, Department of Chemistry, University of Texas at Austin, Austin, TX 78712, USA
| |
Collapse
|
6
|
Zubova EA, Strelnikov IA. Experimental detection of conformational transitions between forms of DNA: problems and prospects. Biophys Rev 2023; 15:1053-1078. [PMID: 37974981 PMCID: PMC10643659 DOI: 10.1007/s12551-023-01143-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2023] [Accepted: 09/06/2023] [Indexed: 11/19/2023] Open
Abstract
Under different conditions, the DNA double helix can take different geometric forms. Of the large number of its conformations, in addition to the "canonical" B form, the A, C, and Z forms are widely known, and the D, Hoogsteen, and X forms are less known. DNA locally takes the A, C, and Z forms in the cell, in complexes with proteins. We compare different methods for detecting non-canonical DNA conformations: X-ray, IR, and Raman spectroscopy, linear and circular dichroism in both the infrared and ultraviolet regions, as well as NMR (measurement of chemical shifts and their anisotropy, scalar and residual dipolar couplings and inter-proton distances from NOESY (nuclear Overhauser effect spectroscopy) data). We discuss the difficulties in applying these methods, the problems of theoretical interpretation of the experimental results, and the prospects for reliable identification of non-canonical DNA conformations.
Collapse
Affiliation(s)
- Elena A. Zubova
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 4 Kosygin St., Moscow, 119991 Russia
| | - Ivan A. Strelnikov
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 4 Kosygin St., Moscow, 119991 Russia
| |
Collapse
|
7
|
Battistini F, Sala A, Hospital A, Orozco M. Sequence-Dependent Properties of the RNA Duplex. J Chem Inf Model 2023; 63:5259-5271. [PMID: 37577978 DOI: 10.1021/acs.jcim.3c00741] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/15/2023]
Abstract
Sequence-dependent properties of the DNA duplex have been accurately described using extensive molecular dynamics simulations. The RNA duplex meanwhile─which is typically represented as a sequence-averaged rigid rod─does not benefit from having equivalent molecular dynamics simulations. In this paper, we present a massive simulation effort using a set of ABC-optimized duplexes from which we derived tetramer-resolution properties of the RNA duplex and a simple mesoscopic model that can represent elastic properties of long RNA duplexes. Despite the extreme chemical similarity between DNA and RNA, the local and global elastic properties of the duplexes are very different. DNA duplexes show a complex and nonelastic pattern of flexibility, for instance, while RNA duplexes behave as an elastic system whose deformations can be represented by simple harmonic potentials. In RNA duplexes (RNA2), not only are intra- and interbase pair parameters (equilibrium and mechanical) different from those in the equivalent DNA duplex sequences (DNA2) but the correlations between movements also differ. Simple statements on the relative flexibility or stability of both polymers are meaningless and should be substituted by a more detailed description depending on the sequence and the type of deformation considered.
Collapse
Affiliation(s)
- Federica Battistini
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac 10, Barcelona 08028, Spain
- Departament de Bioquímica i Biomedicina. Facultat de Biologia, Universitat de Barcelona, Avgda Diagonal 647, Barcelona 08028, Spain
| | - Alba Sala
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac 10, Barcelona 08028, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac 10, Barcelona 08028, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac 10, Barcelona 08028, Spain
- Departament de Bioquímica i Biomedicina. Facultat de Biologia, Universitat de Barcelona, Avgda Diagonal 647, Barcelona 08028, Spain
| |
Collapse
|
8
|
Alexiou TS, Likos CN. Effective Interactions between Double-Stranded DNA Molecules in Aqueous Electrolyte Solutions: Effects of Molecular Architecture and Counterion Valency. J Phys Chem B 2023; 127:6969-6981. [PMID: 37493448 PMCID: PMC10424236 DOI: 10.1021/acs.jpcb.3c02216] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2023] [Revised: 06/12/2023] [Indexed: 07/27/2023]
Abstract
A computational investigation of the effects of molecular topology, namely, linear and circular, as well as counterion valency, on the ensuing pairwise effective interactions between DNA molecules in an unlinked state is presented. Umbrella sampling simulations have been performed through the introduction of bias potential along a reaction coordinate defined as the distance between the centers-of-mass of pairs of DNA molecules, and effective pair interaction potentials have been computed by employing the weighted histogram analysis method. An interesting comparison can be drawn between the different DNA topologies studied here, especially with regard to the contrasting effects of divalent counterions on the effective pair potentials: while DNA-DNA repulsion in short center-of-mass distances decreases significantly in the presence of divalent counterion-ions (as compared to monovalent ions) for linear DNA, the opposite effect occurs for the DNA minicircles. This can be attributed to the fact that linear DNA fragments can easily adopt relative orientations that minimize electrostatic and steric repulsions by rotating relative to one another and by exhibiting more pronounced bending due to the presence of free ends.
Collapse
Affiliation(s)
| | - Christos N Likos
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| |
Collapse
|
9
|
Genna V, Iglesias-Fernández J, Reyes-Fraile L, Villegas N, Guckian K, Seth P, Wan B, Cabrero C, Terrazas M, Brun-Heath I, González C, Sciabola S, Villalobos A, Orozco M. Controlled sulfur-based engineering confers mouldability to phosphorothioate antisense oligonucleotides. Nucleic Acids Res 2023; 51:4713-4725. [PMID: 37099382 PMCID: PMC10250214 DOI: 10.1093/nar/gkad309] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Revised: 04/04/2023] [Accepted: 04/24/2023] [Indexed: 04/27/2023] Open
Abstract
Phosphorothioates (PS) have proven their effectiveness in the area of therapeutic oligonucleotides with applications spanning from cancer treatment to neurodegenerative disorders. Initially, PS substitution was introduced for the antisense oligonucleotides (PS ASOs) because it confers an increased nuclease resistance meanwhile ameliorates cellular uptake and in-vivo bioavailability. Thus, PS oligonucleotides have been elevated to a fundamental asset in the realm of gene silencing therapeutic methodologies. But, despite their wide use, little is known on the possibly different structural changes PS-substitutions may provoke in DNA·RNA hybrids. Additionally, scarce information and significant controversy exists on the role of phosphorothioate chirality in modulating PS properties. Here, through comprehensive computational investigations and experimental measurements, we shed light on the impact of PS chirality in DNA-based antisense oligonucleotides; how the different phosphorothioate diastereomers impact DNA topology, stability and flexibility to ultimately disclose pro-Sp S and pro-Rp S roles at the catalytic core of DNA Exonuclease and Human Ribonuclease H; two major obstacles in ASOs-based therapies. Altogether, our results provide full-atom and mechanistic insights on the structural aberrations PS-substitutions provoke and explain the origin of nuclease resistance PS-linkages confer to DNA·RNA hybrids; crucial information to improve current ASOs-based therapies.
Collapse
Affiliation(s)
- Vito Genna
- Mechanisms of Diseases, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, Barcelona 08028, Spain
- NBD | Nostrum Biodiscovery, Baldiri Reixac 10, Barcelona 08028, Spain
| | | | - Laura Reyes-Fraile
- Mechanisms of Diseases, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, Barcelona 08028, Spain
| | - Nuria Villegas
- Mechanisms of Diseases, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, Barcelona 08028, Spain
| | | | - Punit Seth
- Ionis Pharmaceuticals Inc., 2855 Gazelle Court, Carlsbad, CA 92010, USA
| | - Brad Wan
- Ionis Pharmaceuticals Inc., 2855 Gazelle Court, Carlsbad, CA 92010, USA
| | - Cristina Cabrero
- Instituto de Química Física Rocasolano, C/ Serrano 119, Madrid 28006, Spain
| | - Montserrat Terrazas
- Mechanisms of Diseases, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, Barcelona 08028, Spain
- Department of Inorganic and Organic Chemistry, Section of Organic Chemistry, IBUB, University of Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain
| | - Isabelle Brun-Heath
- Mechanisms of Diseases, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, Barcelona 08028, Spain
| | - Carlos González
- Instituto de Química Física Rocasolano, C/ Serrano 119, Madrid 28006, Spain
| | | | | | - Modesto Orozco
- Mechanisms of Diseases, Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, Barcelona 08028, Spain
- Department of Biochemistry and Biomedicine, University of Barcelona, Barcelona 08028, Spain
| |
Collapse
|
10
|
Bignon E, Grandemange S, Dumont E, Monari A. How SARS-CoV-2 Alters the Regulation of Gene Expression in Infected Cells. J Phys Chem Lett 2023; 14:3199-3207. [PMID: 36971439 PMCID: PMC10068877 DOI: 10.1021/acs.jpclett.3c00582] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2023] [Accepted: 03/23/2023] [Indexed: 06/18/2023]
Abstract
Nonstructural accessory proteins in viruses play a key role in hijacking the basic cellular mechanisms, which is essential to promote the virus survival and evasion of the immune system. The immonuglobulin-like open reading frame 8 (ORF8) protein expressed by SARS-CoV-2 accumulates in the nucleus and may influence the regulation of the gene expression in infected cells. In this contribution, by using microsecond time-scale all-atom molecular dynamics simulations, we unravel the structural bases behind the epigenetic action of ORF8. In particular, we highlight how the protein is able to form stable aggregates with DNA through a histone tail-like motif, and how this interaction is influenced by post-translational modifications, such as acetylation and methylation, which are known epigenetic markers in histones. Our work not only clarifies the molecular mechanisms behind the perturbation of the epigenetic regulation caused by the viral infection but also offers an unusual perspective which may foster the development of original antivirals.
Collapse
Affiliation(s)
- Emmanuelle Bignon
- Université
de Lorraine and CNRS, UMR 7019 LPCT, F-54000 Nancy, France
| | | | - Elise Dumont
- Université
Côte d’Azur, Institut de Chimie de Nice, UMR 7272, Parc Valrose, 28 avenue Valrose, F-06108 Nice, France
- Institut
Universitaire de France, 5 rue Descartes, F-75005 Paris, France
| | - Antonio Monari
- Université
Paris Cité and CNRS, ITODYS, F-75006 Paris, France
| |
Collapse
|
11
|
Liebl K, Zacharias M. Toward Force Fields with Improved Base Stacking Descriptions. J Chem Theory Comput 2023; 19:1529-1536. [PMID: 36795949 DOI: 10.1021/acs.jctc.2c01121] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
Recent DNA force fields indicate good performance in describing flexibility and structural stability of double-stranded B-DNA. However, it is not clear how accurately base stacking interactions are represented that are critical for simulating structure formation processes and conformational changes. Based on the equilibrium nucleoside association and base pair nicking, we find that the recent Tumuc1 force field improves the description of base stacking compared to previous state-of-the-art force fields. Nevertheless, base pair stacking is still overstabilized compared to experiment. We propose a rapid method to reweight calculated free energies of stacking upon force field modifications in order to generate improved parameters. A decrease of the Lennard-Jones attraction between nucleo-bases alone appears insufficient; however, adjustments in the partial charge distribution on base atoms could help to further improve the force field description of base stacking.
Collapse
Affiliation(s)
- Korbinian Liebl
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - Martin Zacharias
- Physics Department and Center of Protein Assemblies, Technical University of Munich, Garching 85748, Germany
| |
Collapse
|
12
|
Nakagawa S, Kimura A, Okamoto Y. Polarizable Molecular Block Model: Toward the Development of an Induced Dipole Force Field for DNA. J Phys Chem B 2022; 126:10646-10661. [PMID: 36512703 DOI: 10.1021/acs.jpcb.2c06227] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
For flexible and highly ionized macromolecules such as DNA, it is important to correctly evaluate the intramolecular polarization in an induced dipole force field. In a proposed polarizable molecular block (PMB) model, a large molecule is divided into several molecular blocks. The atomic charges of the blocks are optimized by using the respective electrostatic potentials (ESPs) on the molecular surface. By using the capped hydrogen removal operation, the total charge of the blocks is controlled exactly to have an integer charge. The atomic polarizabilities of the blocks are optimized by using the respective polarized one-electron potentials that are the differences between ESPs with and without an external test charge. Induced dipole-charge interactions between the blocks are all included, but those interactions within the blocks are strictly excluded. All dipole-dipole interactions are included, but the damping functions are applied to the close dipole-dipole pairs. Several types of damping (simple scaling, exponential, linear, and Gaussian) are evaluated. The validity of the PMB model was verified by using trinucleotide duplexes which have A-, B-, and Z-DNA forms. The reference energies of trinucleotide duplexes including counterions (GGT3Na-ACC3Na, GAC3Na-GTC3Na, and GCG3Na-CGC3Na) are calculated using ωB97XD/aug-cc-pVDZ. All damping types reproduced well the reference interaction energies, dipole moments, and ESPs. Among them, the simple scaling with strong attenuation to 1-2 atomic pairs showed the highest stability against the polarization catastrophe. This study shows that it is possible to develop a high-quality polarizable force field by treating the intramolecular polarization on a block-by-block basis.
Collapse
Affiliation(s)
- Setsuko Nakagawa
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi464-8602, Japan.,Kinjo Gakuin University, Nagoya, Aichi463-8521, Japan
| | - Akihiro Kimura
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi464-8602, Japan
| | - Yuko Okamoto
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi464-8602, Japan.,High Performance Computing Division, Information Technology Center, Nagoya University, Nagoya, Aichi464-8601, Japan.,Global Engagement Center, International Affairs, Nagoya University, Nagoya, Aichi464-8601, Japan
| |
Collapse
|
13
|
Berryman JT, Taghavi A, Mazur F, Tkatchenko A. Quantum machine learning corrects classical forcefields: Stretching DNA base pairs in explicit solvent. J Chem Phys 2022; 157:064107. [PMID: 35963717 DOI: 10.1063/5.0094727] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In order to improve the accuracy of molecular dynamics simulations, classical forcefields are supplemented with a kernel-based machine learning method trained on quantum-mechanical fragment energies. As an example application, a potential-energy surface is generalized for a small DNA duplex, taking into account explicit solvation and long-range electron exchange-correlation effects. A long-standing problem in molecular science is that experimental studies of the structural and thermodynamic behavior of DNA under tension are not well confirmed by simulation; study of the potential energy vs extension taking into account a novel correction shows that leading classical DNA models have excessive stiffness with respect to stretching. This discrepancy is found to be common across multiple forcefields. The quantum correction is in qualitative agreement with the experimental thermodynamics for larger DNA double helices, providing a candidate explanation for the general and long-standing discrepancy between single molecule stretching experiments and classical calculations of DNA stretching. The new dataset of quantum calculations should facilitate multiple types of nucleic acid simulation, and the associated Kernel Modified Molecular Dynamics method (KMMD) is applicable to biomolecular simulations in general. KMMD is made available as part of the AMBER22 simulation software.
Collapse
Affiliation(s)
- Joshua T Berryman
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Amirhossein Taghavi
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Florian Mazur
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
14
|
Galindo-Murillo R, Cheatham TE. Transient Hoogsteen Base Pairs Observed in Unbiased Molecular Dynamics Simulations of DNA. J Phys Chem Lett 2022; 13:6283-6287. [PMID: 35775742 PMCID: PMC9466961 DOI: 10.1021/acs.jpclett.2c01348] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Duplex DNA is modeled as canonical B-DNA displaying the characteristic Watson-Crick base pairs. A less common and short-lived pairing of the nucleobases is the Hoogsteen (HG) conformation. The low population of the HG base pairs (<1%) necessitates extended sampling times in order to analyze through unbiased molecular dynamics (MD) simulations. We have discovered that with extended sampling times using multiple independent copies of an 18-mer sequence, the MD trajectories reproduce the expected and transient HG base pairing. Consistent with experiment, the percentage of the HG events are within the range of ∼0.1-1.0% over the combined aggregate sampling time of more than 3.6 ms. We present the reliability of the current AMBER set of nucleic acid force fields and tools to accurately simulate naturally occurring base-pairing and opening events without any bias or restraints. The mechanism consists of base pair fraying, flipping of the purine, and reformation with HG base pairs.
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, Utah 84112, United States
| | - Thomas E. Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 306, Salt Lake City, Utah 84112, United States
| |
Collapse
|
15
|
Tucker MR, Piana S, Tan D, LeVine MV, Shaw DE. Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA-Protein Complexes. J Phys Chem B 2022; 126:4442-4457. [PMID: 35694853 PMCID: PMC9234960 DOI: 10.1021/acs.jpcb.1c10971] [Citation(s) in RCA: 25] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
![]()
Although molecular
dynamics (MD) simulations have been used extensively
to study the structural dynamics of proteins, the role of MD simulation
in studies of nucleic acid based systems has been more limited. One
contributing factor to this disparity is the historically lower level
of accuracy of the physical models used in such simulations to describe
interactions involving nucleic acids. By modifying nonbonded and torsion
parameters of a force field from the Amber family of models, we recently
developed force field parameters for RNA that achieve a level of accuracy
comparable to that of state-of-the-art protein force fields. Here
we report force field parameters for DNA, which we developed by transferring
nonbonded parameters from our recently reported RNA force field and
making subsequent adjustments to torsion parameters. We have also
modified the backbone charges in both the RNA and DNA parameter sets
to make the treatment of electrostatics compatible with our recently
developed variant of the Amber protein and ion force field. We name
the force field resulting from the union of these three parameter
sets (the new DNA parameters, the revised RNA parameters, and the
existing protein and ion parameters) DES-Amber. Extensive
testing of DES-Amber indicates that it can describe the thermal stability
and conformational flexibility of single- and double-stranded DNA
systems with a level of accuracy comparable to or, especially for
disordered systems, exceeding that of state-of-the-art nucleic acid
force fields. Finally, we show that, in certain favorable cases, DES-Amber
can be used for long-timescale simulations of protein–nucleic
acid complexes.
Collapse
Affiliation(s)
| | - Stefano Piana
- D. E. Shaw Research, New York, New York 10036, United States
| | - Dazhi Tan
- D. E. Shaw Research, New York, New York 10036, United States
| | | | - David E Shaw
- D. E. Shaw Research, New York, New York 10036, United States.,Department of Biochemistry and Molecular Biophysics, Columbia University, New York, New York 10032, United States
| |
Collapse
|
16
|
Segers M, Voorspoels A, Sakaue T, Carlon E. Mechanical properties of nucleic acids and the non-local twistable wormlike chain model. J Chem Phys 2022; 156:234105. [PMID: 35732531 DOI: 10.1063/5.0089166] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Mechanical properties of nucleic acids play an important role in many biological processes that often involve physical deformations of these molecules. At sufficiently long length scales (say, above ∼20-30 base pairs), the mechanics of DNA and RNA double helices is described by a homogeneous Twistable Wormlike Chain (TWLC), a semiflexible polymer model characterized by twist and bending stiffnesses. At shorter scales, this model breaks down for two reasons: the elastic properties become sequence-dependent and the mechanical deformations at distal sites get coupled. We discuss in this paper the origin of the latter effect using the framework of a non-local Twistable Wormlike Chain (nlTWLC). We show, by comparing all-atom simulations data for DNA and RNA double helices, that the non-local couplings are of very similar nature in these two molecules: couplings between distal sites are strong for tilt and twist degrees of freedom and weak for roll. We introduce and analyze a simple double-stranded polymer model that clarifies the origin of this universal distal couplings behavior. In this model, referred to as the ladder model, a nlTWLC description emerges from the coarsening of local (atomic) degrees of freedom into angular variables that describe the twist and bending of the molecule. Different from its local counterpart, the nlTWLC is characterized by a length-scale-dependent elasticity. Our analysis predicts that nucleic acids are mechanically softer at the scale of a few base pairs and are asymptotically stiffer at longer length scales, a behavior that matches experimental data.
Collapse
Affiliation(s)
- Midas Segers
- Soft Matter and Biophysics Unit, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
| | - Aderik Voorspoels
- Soft Matter and Biophysics Unit, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
| | - Takahiro Sakaue
- Department of Physical Sciences, Aoyama Gakuin University, 5-10-1 Fuchinobe, Chuo-ku, Sagamihara 252-5258, Kanagawa, Japan
| | - Enrico Carlon
- Soft Matter and Biophysics Unit, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
| |
Collapse
|
17
|
Wieczór M, Genna V, Aranda J, Badia RM, Gelpí JL, Gapsys V, de Groot BL, Lindahl E, Municoy M, Hospital A, Orozco M. Pre-exascale HPC approaches for molecular dynamics simulations. Covid-19 research: A use case. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2022; 13:e1622. [PMID: 35935573 PMCID: PMC9347456 DOI: 10.1002/wcms.1622] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/02/2022] [Revised: 04/25/2022] [Accepted: 04/28/2022] [Indexed: 06/15/2023]
Abstract
Exascale computing has been a dream for ages and is close to becoming a reality that will impact how molecular simulations are being performed, as well as the quantity and quality of the information derived for them. We review how the biomolecular simulations field is anticipating these new architectures, making emphasis on recent work from groups in the BioExcel Center of Excellence for High Performance Computing. We exemplified the power of these simulation strategies with the work done by the HPC simulation community to fight Covid-19 pandemics. This article is categorized under:Data Science > Computer Algorithms and ProgrammingData Science > Databases and Expert SystemsMolecular and Statistical Mechanics > Molecular Dynamics and Monte-Carlo Methods.
Collapse
Affiliation(s)
- Miłosz Wieczór
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and TechnologyBarcelonaSpain
- Department of Physical ChemistryGdansk University of TechnologyGdańskPoland
| | - Vito Genna
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and TechnologyBarcelonaSpain
| | - Juan Aranda
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and TechnologyBarcelonaSpain
| | | | - Josep Lluís Gelpí
- Barcelona Supercomputing CenterBarcelonaSpain
- Department of Biochemistry and BiomedicineUniversity of BarcelonaBarcelonaSpain
| | - Vytautas Gapsys
- Max Planck Institute for Multidisciplinary SciencesComputational Biomolecular Dynamics GroupGoettingenGermany
| | - Bert L. de Groot
- Max Planck Institute for Multidisciplinary SciencesComputational Biomolecular Dynamics GroupGoettingenGermany
| | - Erik Lindahl
- Department of Applied PhysicsSwedish e‐Science Research Center, KTH Royal Institute of TechnologyStockholmSweden
- Department of Biochemistry and Biophysics, Science for Life LaboratoryStockholm UniversityStockholmSweden
| | | | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and TechnologyBarcelonaSpain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and TechnologyBarcelonaSpain
- Department of Biochemistry and BiomedicineUniversity of BarcelonaBarcelonaSpain
| |
Collapse
|
18
|
Miclot T, Hognon C, Bignon E, Terenzi A, Grandemange S, Barone G, Monari A. Never Cared for What They Do: High Structural Stability of Guanine-Quadruplexes in the Presence of Strand-Break Damage. Molecules 2022; 27:3256. [PMID: 35630732 PMCID: PMC9146567 DOI: 10.3390/molecules27103256] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2022] [Revised: 05/17/2022] [Accepted: 05/18/2022] [Indexed: 11/16/2022] Open
Abstract
DNA integrity is an important factor that assures genome stability and, more generally, the viability of cells and organisms. In the presence of DNA damage, the normal cell cycle is perturbed when cells activate their repair processes. Although efficient, the repair system is not always able to ensure complete restoration of gene integrity. In these cases, mutations not only may occur, but the accumulation of lesions can either lead to carcinogenesis or reach a threshold that induces apoptosis and programmed cell death. Among the different types of DNA lesions, strand breaks produced by ionizing radiation are the most toxic due to the inherent difficultly of repair, which may lead to genomic instability. In this article we show, by using classical molecular simulation techniques, that compared to canonical double-helical B-DNA, guanine-quadruplex (G4) arrangements show remarkable structural stability, even in the presence of two strand breaks. Since G4-DNA is recognized for its regulatory roles in cell senescence and gene expression, including oncogenes, this stability may be related to an evolutionary cellular response aimed at minimizing the effects of ionizing radiation.
Collapse
Affiliation(s)
- Tom Miclot
- Department of Biological, Chemical and Pharmaceutical Sciences, University of Palermo, viale delle Scienze, Ed. 17, 90128 Palermo, Italy; (T.M.); (A.T.)
- Université de Lorraine and CNRS, LPCT UMR 7019, F-54000 Nancy, France; (C.H.); (E.B.)
| | - Cécilia Hognon
- Université de Lorraine and CNRS, LPCT UMR 7019, F-54000 Nancy, France; (C.H.); (E.B.)
| | - Emmanuelle Bignon
- Université de Lorraine and CNRS, LPCT UMR 7019, F-54000 Nancy, France; (C.H.); (E.B.)
| | - Alessio Terenzi
- Department of Biological, Chemical and Pharmaceutical Sciences, University of Palermo, viale delle Scienze, Ed. 17, 90128 Palermo, Italy; (T.M.); (A.T.)
| | | | - Giampaolo Barone
- Department of Biological, Chemical and Pharmaceutical Sciences, University of Palermo, viale delle Scienze, Ed. 17, 90128 Palermo, Italy; (T.M.); (A.T.)
| | - Antonio Monari
- Université Paris Cité and CNRS, ITODYS, F-75006 Paris, France
| |
Collapse
|
19
|
Schneider AL, Albrecht AV, Huang K, Germann MW, Poon GMK. Self-Consistent Parameterization of DNA Residues for the Non-Polarizable AMBER Force Fields. Life (Basel) 2022; 12:life12050666. [PMID: 35629334 PMCID: PMC9143812 DOI: 10.3390/life12050666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Revised: 04/27/2022] [Accepted: 04/27/2022] [Indexed: 11/22/2022] Open
Abstract
Fixed-charge (non-polarizable) forcefields are accurate and computationally efficient tools for modeling the molecular dynamics of nucleic acid polymers, particularly DNA, well into the µs timescale. The continued utility of these forcefields depends in part on expanding the residue set in step with advancing nucleic acid chemistry and biology. A key step in parameterizing new residues is charge derivation which is self-consistent with the existing residues. As atomic charges are derived by fitting against molecular electrostatic potentials, appropriate structural models are critical. Benchmarking against the existing charge set used in current AMBER nucleic acid forcefields, we report that quantum mechanical models of deoxynucleosides, even at a high level of theory, are not optimal structures for charge derivation. Instead, structures from molecular mechanics minimization yield charges with up to 6-fold lower RMS deviation from the published values, due to the choice of such an approach in the derivation of the original charge set. We present a contemporary protocol for rendering self-consistent charges as well as optimized charges for a panel of nine non-canonical residues that will permit comparison with literature as well as studying the dynamics of novel DNA polymers.
Collapse
Affiliation(s)
- Amelia L. Schneider
- Department of Chemistry, Georgia State University, Atlanta, GA 30303, USA; (A.L.S.); (A.V.A.); (K.H.)
| | - Amanda V. Albrecht
- Department of Chemistry, Georgia State University, Atlanta, GA 30303, USA; (A.L.S.); (A.V.A.); (K.H.)
| | - Kenneth Huang
- Department of Chemistry, Georgia State University, Atlanta, GA 30303, USA; (A.L.S.); (A.V.A.); (K.H.)
| | - Markus W. Germann
- Department of Chemistry, Georgia State University, Atlanta, GA 30303, USA; (A.L.S.); (A.V.A.); (K.H.)
- Department of Biology, Georgia State University, Atlanta, GA 30303, USA
- Correspondence: (M.W.G.); (G.M.K.P.)
| | - Gregory M. K. Poon
- Department of Chemistry, Georgia State University, Atlanta, GA 30303, USA; (A.L.S.); (A.V.A.); (K.H.)
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA 30303, USA
- Correspondence: (M.W.G.); (G.M.K.P.)
| |
Collapse
|
20
|
Whitford PC. Overview of the Biomolecular Association and Dynamics session at the 20th IUPAB congress, 45th Brazilian congress of SBBF, and the 50th annual meeting of SBBq. Biophys Rev 2022; 13:863-865. [PMID: 35059009 DOI: 10.1007/s12551-021-00884-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Accepted: 11/01/2021] [Indexed: 05/31/2023] Open
Abstract
In this session, experts in molecular biophysics described the dynamics of biopolymers across a wide range of length and time scales. This discussion highlighted numerous techniques that span from highly detailed simulations, to coarse-grained theoretical models, as well as high-resolution structural analysis. The topics were equally diverse, where there was discussion of biological processes at small (individual atoms), intermediate (assemblies) and very large scales (phase separation).
Collapse
Affiliation(s)
- Paul Charles Whitford
- Department of Physics and Center for Theoretical Biological Physics, Northeastern University, Boston, MA USA
| |
Collapse
|
21
|
Rajasekaran T, Freestone GC, Galindo-Murillo R, Lugato B, Rico L, Salinas JC, Gaus H, Migawa MT, Swayze EE, Cheatham TE, Hanessian S, Seth PP. Backbone Hydrocarbon-Constrained Nucleic Acids Modulate Hybridization Kinetics for RNA. J Am Chem Soc 2022; 144:1941-1950. [PMID: 35041415 DOI: 10.1021/jacs.1c12323] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
The binding affinity of therapeutic oligonucleotides (ONs) for their cognate RNA is determined by the rates of association (ka) and dissociation (kd). Single-stranded ONs are highly flexible and can adopt multiple conformations in solution, some of which may not be conducive for hybridization. We investigated if restricting rotation around the sugar-phosphate backbone, by tethering two adjacent backbone phosphonate esters using hydrocarbon bridges, can modulate hybridization kinetics of the modified ONs for complementary RNA. Given the large number of possible analogues with different tether lengths and configurations at the phosphorus atoms, we employed molecular dynamic simulations to optimize the size of the hydrocarbon bridge to guide the synthetic efforts. The backbone-constrained nucleotide trimers with stereodefined configurations at the contiguous backbone phosphorus atoms were assembled using a ring-closing metathesis reaction, then incorporated into oligonucleotides by an in situ synthesis of the phosphoramidites followed by coupling to solid supports. Evaluation of the modified oligonucleotides revealed that 15-membered macrocyclic-constrained analogues displayed similar or slightly improved on-rates but significantly increased off-rates compared to unmodified DNA ONs, resulting in reduced duplex stability. In contrast, LNA ONs with conformationally preorganized furanose rings showed similar on-rates to DNA ONs but very slow off-rates, resulting in net improvement in duplex stability. Furthermore, the experimental data generally supported the molecular dynamics simulation results, suggesting that this strategy can be used as a predictive tool for designing the next generation of constrained backbone ON analogues with improved hybridization properties.
Collapse
Affiliation(s)
| | - Graeme C Freestone
- Department of Medicinal Chemistry, Ionis Pharmaceuticals, 2855 Gazelle Court, Carlsbad, California 92010, United States
| | - Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 201, Salt Lake City, Utah 84112, United States
| | - Barbara Lugato
- Department of Chemistry, Université de Montréal, Quebec H3C 3J7, Canada
| | - Lorena Rico
- Department of Chemistry, Université de Montréal, Quebec H3C 3J7, Canada
| | - Juan C Salinas
- Department of Chemistry, Université de Montréal, Quebec H3C 3J7, Canada
| | - Hans Gaus
- Department of Medicinal Chemistry, Ionis Pharmaceuticals, 2855 Gazelle Court, Carlsbad, California 92010, United States
| | - Michael T Migawa
- Department of Medicinal Chemistry, Ionis Pharmaceuticals, 2855 Gazelle Court, Carlsbad, California 92010, United States
| | - Eric E Swayze
- Department of Medicinal Chemistry, Ionis Pharmaceuticals, 2855 Gazelle Court, Carlsbad, California 92010, United States
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 201, Salt Lake City, Utah 84112, United States
| | - Stephen Hanessian
- Department of Chemistry, Université de Montréal, Quebec H3C 3J7, Canada
| | - Punit P Seth
- Department of Medicinal Chemistry, Ionis Pharmaceuticals, 2855 Gazelle Court, Carlsbad, California 92010, United States
| |
Collapse
|
22
|
Liebl K, Zacharias M. Tumuc1: A New Accurate DNA Force Field Consistent with High-Level Quantum Chemistry. J Chem Theory Comput 2021; 17:7096-7105. [PMID: 34662102 DOI: 10.1021/acs.jctc.1c00682] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An accurate molecular mechanics force field forms the basis of Molecular Dynamics simulations to obtain a realistic view of the structure and dynamics of biomolecules such as DNA. Although frequently updated to improve agreement with available experimental data, DNA force fields still rely in part on parameters introduced more than 20 years ago. We have developed an entirely new DNA force field, Tumuc1, derived from quantum mechanical calculations to obtain a consistent set of bonded parameters and partial atomic charges. The performance of the force field was extensively tested on a variety of DNA molecules. It excels in accuracy of B-DNA simulations but also performs very well on other types of DNA structures and structure formation processes such as hairpin folding, duplex formation, and dynamics of DNA-protein complexes. It can complement existing force fields in order to provide an increasingly accurate description of the structure and dynamics of DNA during simulation studies.
Collapse
Affiliation(s)
- Korbinian Liebl
- Physics Department and Center of Protein Assemblies, Technical University of Munich, James-Franck-Str. 1, Garching 85748, Germany
| | - Martin Zacharias
- Physics Department and Center of Protein Assemblies, Technical University of Munich, James-Franck-Str. 1, Garching 85748, Germany
| |
Collapse
|
23
|
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
|
24
|
Donati E, Vidossich P, De Vivo M. Molecular Mechanism of Phosphate Steering for DNA Binding, Cleavage Localization, and Substrate Release in Nucleases. ACS Catal 2021. [DOI: 10.1021/acscatal.1c03346] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Elisa Donati
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, via Morego 30, 16163 Genoa, Italy
| | - Pietro Vidossich
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, via Morego 30, 16163 Genoa, Italy
| | - Marco De Vivo
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, via Morego 30, 16163 Genoa, Italy
| |
Collapse
|
25
|
Strelnikov IA, Kovaleva NA, Zubova EA. Variability of the DNA Backbone Geometry in DNA-Protein Complexes: Experimental Data Analysis. J Chem Inf Model 2021; 61:4783-4794. [PMID: 34529915 DOI: 10.1021/acs.jcim.1c00506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We have analyzed and compared the available experimental data (PDB) on the backbone geometry of the DNA in solution (NMR), in crystals (X-rays), and in complexes with proteins (X-rays and cryo-electron microscopy). The deoxyribose (pseudorotational angle τ0) and ε/ζ (BI-BII transition in phosphates) flexibilities are practically the same in the four samples. The α/γ mobility is minimal in crystalline DNA: on the histograms, there is one canonical and one noncanonical t/t peak. The α/γ mobility increases in DNA solutions (three more noncanonical peaks) and is maximal in DNA-protein complexes (another additional peak). On a large amount of data, we have confirmed that the three main degrees of freedom of the sugar-phosphate backbone are "orthogonal": changes in any of the angles τ0, (ζ-ε), and (γ-α) occur, as a rule, at a constant (usually canonical) value of any other. In the DNA-protein complexes, none of the geometrical parameters commonly used to distinguish the A and B forms of DNA, except for Zp and its simpler analog Zp', show an unambiguous correlation with τ0. Proteins, binding to DNA, in 59% of cases change the local shape of the helix up to the characteristic of the A-form without switching the deoxyribose conformation from south to north. However, we have found simple local characteristics of one nucleotide that correlate with the angles τ0 and (ζ-ε). These are the angles C3'C1'N* and C4'C3'P(2), respectively. They are orthogonal in DNA-protein complexes exactly as the pair τ0 and (ζ-ε). Most characteristics of DNA in complexes with proteins are the same in X-ray and in cryo-EM data, except for the histogram for the angle τ0. We offer a possible explanation for this difference. We also discuss the artifacts on the ε/ζ histogram for DNA in solutions caused by the currently used NMR refinement protocols.
Collapse
Affiliation(s)
- Ivan A Strelnikov
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 4 Kosygin Street, Moscow 119991, Russia
| | - Natalya A Kovaleva
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 4 Kosygin Street, Moscow 119991, Russia
| | - Elena A Zubova
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 4 Kosygin Street, Moscow 119991, Russia
| |
Collapse
|
26
|
Hognon C, Bignon E, Harle G, Touche N, Grandemange S, Monari A. The Iron Maiden. Cytosolic Aconitase/IRP1 Conformational Transition in the Regulation of Ferritin Translation and Iron Hemostasis. Biomolecules 2021; 11:biom11091329. [PMID: 34572542 PMCID: PMC8469783 DOI: 10.3390/biom11091329] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2021] [Revised: 08/29/2021] [Accepted: 09/07/2021] [Indexed: 01/16/2023] Open
Abstract
Maintaining iron homeostasis is fundamental for almost all living beings, and its deregulation correlates with severe and debilitating pathologies. The process is made more complicated by the omnipresence of iron and by its role as a fundamental component of a number of crucial metallo proteins. The response to modifications in the amount of the free-iron pool is performed via the inhibition of ferritin translation by sequestering consensus messenger RNA (mRNA) sequences. In turn, this is regulated by the iron-sensitive conformational equilibrium between cytosolic aconitase and IRP1, mediated by the presence of an iron-sulfur cluster. In this contribution, we analyze by full-atom molecular dynamics simulation, the factors leading to both the interaction with mRNA and the conformational transition. Furthermore, the role of the iron-sulfur cluster in driving the conformational transition is assessed by obtaining the related free energy profile via enhanced sampling molecular dynamics simulations.
Collapse
Affiliation(s)
- Cécilia Hognon
- Université de Lorraine and CNRS, UMR 7019 LPCT, F-54000 Nancy, France; (C.H.); (E.B.)
| | - Emmanuelle Bignon
- Université de Lorraine and CNRS, UMR 7019 LPCT, F-54000 Nancy, France; (C.H.); (E.B.)
| | - Guillaume Harle
- Université de Lorraine and CNRS, UMR 7039 CRAN, F-54000 Nancy, France; (G.H.); (N.T.)
| | - Nadège Touche
- Université de Lorraine and CNRS, UMR 7039 CRAN, F-54000 Nancy, France; (G.H.); (N.T.)
| | - Stéphanie Grandemange
- Université de Lorraine and CNRS, UMR 7039 CRAN, F-54000 Nancy, France; (G.H.); (N.T.)
- Correspondence: (S.G.); (A.M.)
| | - Antonio Monari
- Université de Lorraine and CNRS, UMR 7019 LPCT, F-54000 Nancy, France; (C.H.); (E.B.)
- Université de Paris and CNRS, ITODYS, F-75006 Paris, France
- Correspondence: (S.G.); (A.M.)
| |
Collapse
|
27
|
Bie LH, Fei JW, Gao J. Molecular mechanism of methyl-dependent and spatial-specific DNA recognition of c-Jun homodimer. J Mol Model 2021; 27:227. [PMID: 34264385 DOI: 10.1007/s00894-021-04840-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Accepted: 07/01/2021] [Indexed: 10/20/2022]
Abstract
DNA methylation is important in regulation of gene expression and normal development because it alters the interplay between protein and DNA. Experiments have shown that a single 5-methylcytosine at different CpG sites (mCpG) might have different effects on specific recognition, but the atomistic origin and dynamic details are largely unclear. In this work, we investigated the mechanism of monomethylation at different CpG sites in the cognate motif and the cooperativity of full methylation. By constructing four models of c-Jun/Jun protein binding to the 5[Formula: see text]-XGAGTCA-3[Formula: see text] (X represents C or methylated C) motif, we characterized the dynamics of the contact interface using the all-atom molecular dynamics method. Free energy analysis of MM/GBSA suggests that regardless of whether the C12pG13 site of the bottom strand is methylated, the effects from mC25 of the top strand are dominant and can moderately enhance the binding by [Formula: see text] 31 kcal/mol, whereas mC12 showed a relatively small contribution, in agreement with the experimental data. Remarkably, we found that this spatial-specific influence was induced by different regulatory rules. The influence of the mC25 site is mainly mediated by steric hindrance. The additional methyl group leads to the conformational changes in nearby residues and triggers an obvious structural bending in the protein, which results in the formation of a new T-Asn-C triad that enhances the specific recognition of TCA half-sites. The substitution of the methyl group at the mC12 site of the bottom strand breaks the original H-bonds directly. Such changes in electrostatic interactions also lead to the remote allosteric effects of protein by multifaceted interactions but have negligible contributions to binding. Although these two influence modes are different, they can both fine-tune the local environment, which might produce remote allosteric effects through protein-protein interactions. Further analysis reveals that the discrepancies in these two modes are primarily due to their location. Moreover, when both sites are methylated, the major determinant of binding specificity depends on the context and the location of the methylation site, which is the result of crosstalk and cooperativity.
Collapse
Affiliation(s)
- Li-Hua Bie
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, 430070, People's Republic of China
| | - Jun-Wen Fei
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, 430070, People's Republic of China
| | - Jun Gao
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, 430070, People's Republic of China.
| |
Collapse
|
28
|
Hognon C, Monari A. Staring at the Naked Goddess: Unraveling the Structure and Reactivity of Artemis Endonuclease Interacting with a DNA Double Strand. Molecules 2021; 26:molecules26133986. [PMID: 34210101 PMCID: PMC8271620 DOI: 10.3390/molecules26133986] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Revised: 06/23/2021] [Accepted: 06/26/2021] [Indexed: 11/16/2022] Open
Abstract
Artemis is an endonuclease responsible for breaking hairpin DNA strands during immune system adaptation and maturation as well as the processing of potentially toxic DNA lesions. Thus, Artemis may be an important target in the development of anticancer therapy, both for the sensitization of radiotherapy and for immunotherapy. Despite its importance, its structure has been resolved only recently, and important questions concerning the arrangement of its active center, the interaction with the DNA substrate, and the catalytic mechanism remain unanswered. In this contribution, by performing extensive molecular dynamic simulations, both classically and at the hybrid quantum mechanics/molecular mechanics level, we evidenced the stable interaction modes of Artemis with a model DNA strand. We also analyzed the catalytic cycle providing the free energy profile and key transition states for the DNA cleavage reaction.
Collapse
|
29
|
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: 51] [Impact Index Per Article: 17.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
|
30
|
Beyond the double helix: DNA structural diversity and the PDB. J Biol Chem 2021; 296:100553. [PMID: 33744292 PMCID: PMC8063756 DOI: 10.1016/j.jbc.2021.100553] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Revised: 01/15/2021] [Accepted: 03/16/2021] [Indexed: 12/11/2022] Open
Abstract
The determination of the double helical structure of DNA in 1953 remains the landmark event in the development of modern biological and biomedical science. This structure has also been the starting point for the determination of some 2000 DNA crystal structures in the subsequent 68 years. Their structural diversity has extended to the demonstration of sequence-dependent local structure in duplex DNA, to DNA bending in short and long sequences and in the DNA wound round the nucleosome, and to left-handed duplex DNAs. Beyond the double helix itself, in circumstances where DNA sequences are or can be induced to unwind from being duplex, a wide variety of topologies and forms can exist. Quadruplex structures, based on four-stranded cores of stacked G-quartets, are prevalent though not randomly distributed in the human and other genomes and can play roles in transcription, translation, and replication. Yet more complex folds can result in DNAs with extended tertiary structures and enzymatic/catalytic activity. The Protein Data Bank is the depository of all these structures, and the resource where structures can be critically examined and validated, as well as compared one with another to facilitate analysis of conformational and base morphology features. This review will briefly survey the major structural classes of DNAs and illustrate their significance, together with some examples of how the use of the Protein Data Bank by for example, data mining, has illuminated DNA structural concepts.
Collapse
|
31
|
Golyshev VM, Pyshnyi DV, Lomzov AA. Effects of Phosphoryl Guanidine Modification of Phosphate Residues on the Structure and Hybridization of Oligodeoxyribonucleotides. J Phys Chem B 2021; 125:2841-2855. [PMID: 33724825 DOI: 10.1021/acs.jpcb.0c10214] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Phosphoryl guanidine oligonucleotides (PGOs) are promising tools for biological research and development of biosensors and therapeutics. We performed structural and hybridization analyses of octa-, deca-, and dodecamers with all phosphate residues modified by 1,3-dimethylimidazolidine-2-imine moieties. Similarity of the B-form double helix between native and modified duplexes was noted. In PGO duplexes, we detected a decrease in the proportion of C2'-endo and an increased proportion of C1'-exo sugar conformations of the modified chain. Applicability of the two-state model to denaturation transition of all studied duplexes was proved for the first time. Sequence-dependent effects of this modification on hybridization properties were observed. The thermal stability of PGO complexes is almost native at 100 mM NaCl and slightly increases with decreasing ionic strength. An increase in water activity and dramatic changes in interaction with cations and in solvation of PGOs and their duplexes were noted, resulting in slight elevation of the melting temperature after an ionic-strength decrease from 1 M NaCl down to deionized water. Decreased binding of sodium ions and decreased water solvation were documented for PGOs and their duplexes. In contrast to DNA, the PGO duplex formation leads to a release of several cations. The water shell is significantly more disordered near PGOs and their complexes. Nevertheless, changes in solvation during the formation of native and PGO complexes are similar and indicate that it is possible to develop models for predictive calculations of the thermodynamic properties of phosphoryl guanidine oligomers. Our results may help devise an approach for the rational design of PGOs as novel improved molecular probes and tools for many modern methods involving oligonucleotides.
Collapse
Affiliation(s)
- Victor M Golyshev
- Institute of Chemical Biology and Fundamental Medicine, Novosibirsk 630090, Russia.,Novosibirsk State University, Novosibirsk 630090, Russia
| | - Dmitrii V Pyshnyi
- Institute of Chemical Biology and Fundamental Medicine, Novosibirsk 630090, Russia.,Novosibirsk State University, Novosibirsk 630090, Russia
| | - Alexander A Lomzov
- Institute of Chemical Biology and Fundamental Medicine, Novosibirsk 630090, Russia.,Novosibirsk State University, Novosibirsk 630090, Russia
| |
Collapse
|
32
|
Bouchal T, Durník I, Illík V, Réblová K, Kulhánek P. Importance of base-pair opening for mismatch recognition. Nucleic Acids Res 2020; 48:11322-11334. [PMID: 33080020 PMCID: PMC7672436 DOI: 10.1093/nar/gkaa896] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2020] [Revised: 09/09/2020] [Accepted: 09/30/2020] [Indexed: 01/04/2023] Open
Abstract
Mismatch repair is a highly conserved cellular pathway responsible for repairing mismatched dsDNA. Errors are detected by the MutS enzyme, which most likely senses altered mechanical property of damaged dsDNA rather than a specific molecular pattern. While the curved shape of dsDNA in crystallographic MutS/DNA structures suggests the role of DNA bending, the theoretical support is not fully convincing. Here, we present a computational study focused on a base-pair opening into the minor groove, a specific base-pair motion observed upon interaction with MutS. Propensities for the opening were evaluated in terms of two base-pair parameters: Opening and Shear. We tested all possible base pairs in anti/anti, anti/syn and syn/anti orientations and found clear discrimination between mismatches and canonical base-pairs only for the opening into the minor groove. Besides, the discrimination gap was also confirmed in hotspot and coldspot sequences, indicating that the opening could play a more significant role in the mismatch recognition than previously recognized. Our findings can be helpful for a better understanding of sequence-dependent mutability. Further, detailed structural characterization of mismatches can serve for designing anti-cancer drugs targeting mismatched base pairs.
Collapse
Affiliation(s)
- Tomáš Bouchal
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Ivo Durník
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Viktor Illík
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Kamila Réblová
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Petr Kulhánek
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
33
|
Hörberg J, Reymer A. Specifically bound BZIP transcription factors modulate DNA supercoiling transitions. Sci Rep 2020; 10:18795. [PMID: 33139763 PMCID: PMC7606469 DOI: 10.1038/s41598-020-75711-4] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2020] [Accepted: 10/18/2020] [Indexed: 01/01/2023] Open
Abstract
Torsional stress on DNA, introduced by molecular motors, constitutes an important regulatory mechanism of transcriptional control. Torsional stress can modulate specific binding of transcription factors to DNA and introduce local conformational changes that facilitate the opening of promoters and nucleosome remodelling. Using all-atom microsecond scale molecular dynamics simulations together with a torsional restraint that controls the total twist of a DNA fragment, we address the impact of torsional stress on DNA complexation with a human BZIP transcription factor, MafB. We gradually over- and underwind DNA alone and in complex with MafB by 0.5° per dinucleotide step, starting from the relaxed state to a maximum of 5° per dinucleotide step, monitoring the evolution of the protein-DNA contacts at different degrees of torsional strain. Our computations show that MafB changes the DNA sequence-specific response to torsional stress. The dinucleotide steps that are susceptible to absorbing most of the torsional stress become more torsionally rigid, as they are involved in protein-DNA contacts. Also, the protein undergoes substantial conformational changes to follow the stress-induced DNA deformation, but mostly maintains the specific contacts with DNA. This results in a significant asymmetric increase of free energy of DNA twisting transitions, relative to free DNA, where overtwisting is more energetically unfavourable. Our data suggest that specifically bound BZIP factors could act as torsional stress insulators, modulating the propagation of torsional stress along the chromatin fibre, which might promote cooperative binding of collaborative DNA-binding factors.
Collapse
Affiliation(s)
- Johanna Hörberg
- Department of Chemistry and Molecular Biology, University of Gothenburg, 40530, Gothenburg, Sweden
| | - Anna Reymer
- Department of Chemistry and Molecular Biology, University of Gothenburg, 40530, Gothenburg, Sweden.
| |
Collapse
|
34
|
Peterhans L, Nicolaidou E, Diamantis P, Alloa E, Leclerc M, Surin M, Clément S, Rothlisberger U, Banerji N, Hayes SC. Structural and Photophysical Templating of Conjugated Polyelectrolytes with Single-Stranded DNA. CHEMISTRY OF MATERIALS : A PUBLICATION OF THE AMERICAN CHEMICAL SOCIETY 2020; 32:7347-7362. [PMID: 33122875 PMCID: PMC7587141 DOI: 10.1021/acs.chemmater.0c02251] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/29/2020] [Revised: 08/06/2020] [Indexed: 06/11/2023]
Abstract
A promising approach to influence and control the photophysical properties of conjugated polymers is directing their molecular conformation by templating. We explore here the templating effect of single-stranded DNA oligomers (ssDNAs) on cationic polythiophenes with the goal to uncover the intermolecular interactions that direct the polymer backbone conformation. We have comprehensively characterized the optical behavior and structure of the polythiophenes in conformationally distinct complexes depending on the sequence of nucleic bases and addressed the effect on the ultrafast excited-state relaxation. This, in combination with molecular dynamics simulations, allowed us a detailed atomistic-level understanding of the structure-property correlations. We find that electrostatic and other noncovalent interactions direct the assembly with the polymer, and we identify that optimal templating is achieved with (ideally 10-20) consecutive cytosine bases through numerous π-stacking interactions with the thiophene rings and side groups of the polymer, leading to a rigid assembly with ssDNA, with highly ordered chains and unique optical signatures. Our insights are an important step forward in an effective approach to structural templating and optoelectronic control of conjugated polymers and organic materials in general.
Collapse
Affiliation(s)
- Lisa Peterhans
- Department
of Chemistry and Biochemistry, University
of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Eliana Nicolaidou
- Department
of Chemistry, University of Cyprus, P.O. Box 20537, 1678, Nicosia, Cyprus
| | - Polydefkis Diamantis
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Elisa Alloa
- Department
of Chemistry, University of Cyprus, P.O. Box 20537, 1678, Nicosia, Cyprus
| | - Mario Leclerc
- Department
of Chemistry, Université Laval, G1K 7P4 Quebec
City, Quebec, Canada
| | - Mathieu Surin
- Laboratory
for Chemistry of Novel Materials, Center for Innovation in Materials
and Polymers, University of Mons −
UMONS, 20 Place du Parc, B-7000 Mons, Belgium
| | - Sébastien Clément
- Institut
Charles Gerhardt Montpellier, ICGM, UMR 5253 CNRS, Université de Montpellier, Place Eugène Bataillon, F-34095 Montpellier, Cedex
05, France
| | - Ursula Rothlisberger
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Natalie Banerji
- Department
of Chemistry and Biochemistry, University
of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Sophia C. Hayes
- Department
of Chemistry, University of Cyprus, P.O. Box 20537, 1678, Nicosia, Cyprus
| |
Collapse
|
35
|
Matoušková E, Bignon E, Claerbout VEP, Dršata T, Gillet N, Monari A, Dumont E, Lankaš F. Impact of the Nucleosome Histone Core on the Structure and Dynamics of DNA-Containing Pyrimidine-Pyrimidone (6-4) Photoproduct. J Chem Theory Comput 2020; 16:5972-5981. [PMID: 32810397 DOI: 10.1021/acs.jctc.0c00593] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The pyrimidine-pyrimidone (6-4) photoproduct (64-PP) is an important photoinduced DNA lesion constituting a mutational signature for melanoma. The structural impact of 64-PP on DNA complexed with histones affects the lesion mutagenicity and repair but remains poorly understood. Here we investigate the conformational dynamics of DNA-containing 64-PP within the nucleosome core particle by atomic-resolution molecular dynamics simulations and multiscale data analysis. We demonstrate that the histone core exerts important mechanical restraints that largely decrease global DNA structural fluctuations. However, the local DNA flexibility at the damaged site is enhanced due to imperfect structural adaptation to restraints imposed by the histone core. If 64-PP faces the histone core and is therefore not directly accessible by the repair protein, the complementary strand facing the solvent is deformed and exhibits higher flexibility than the corresponding strand in a naked, undamaged DNA. This may serve as an initial recognition signal for repair. Our simulations also pinpoint the structural role of proximal residues from the truncated histone tails.
Collapse
Affiliation(s)
- Eva Matoušková
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| | - Emmanuelle Bignon
- Université de Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratoire de Chimie, F69342, Lyon, France
| | - Victor E P Claerbout
- Université de Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratoire de Chimie, F69342, Lyon, France
| | - Tomáš Dršata
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| | - Natacha Gillet
- Université de Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratoire de Chimie, F69342, Lyon, France
| | - Antonio Monari
- Université de Lorraine and CNRS, LPCT UMR 7019, F-54000 Nancy, France
| | - Elise Dumont
- Université de Lyon, ENS de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratoire de Chimie, F69342, Lyon, France.,Institut Universitaire de France, 5 rue Descartes, 75005 Paris, France
| | - Filip Lankaš
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| |
Collapse
|
36
|
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
|
37
|
Sanchez JC, Zhang L, Evoli S, Schnicker NJ, Nunez-Hernandez M, Yu L, Wereszczynski J, Pufall MA, Musselman CA. The molecular basis of selective DNA binding by the BRG1 AT-hook and bromodomain. BIOCHIMICA ET BIOPHYSICA ACTA. GENE REGULATORY MECHANISMS 2020; 1863:194566. [PMID: 32376391 PMCID: PMC7350285 DOI: 10.1016/j.bbagrm.2020.194566] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/03/2020] [Revised: 04/21/2020] [Accepted: 04/21/2020] [Indexed: 12/18/2022]
Abstract
The ATP-dependent BAF chromatin remodeling complex plays a critical role in gene regulation by modulating chromatin architecture, and is frequently mutated in cancer. Indeed, subunits of the BAF complex are found to be mutated in >20% of human tumors. The mechanism by which BAF properly navigates chromatin is not fully understood, but is thought to involve a multivalent network of histone and DNA contacts. We previously identified a composite domain in the BRG1 ATPase subunit that is capable of associating with both histones and DNA in a multivalent manner. Mapping the DNA binding pocket revealed that it contains several cancer mutations. Here, we utilize SELEX-seq to investigate the DNA specificity of this composite domain and NMR spectroscopy and molecular modelling to determine the structural basis of DNA binding. Finally, we demonstrate that cancer mutations in this domain alter the mode of DNA association.
Collapse
Affiliation(s)
- Julio C Sanchez
- Department of Biochemistry, Carver College of Medicine, University of Iowa, Iowa City, IA 52242, United States
| | - Liyang Zhang
- Department of Biochemistry, Carver College of Medicine, University of Iowa, Iowa City, IA 52242, United States; Integrated DNA Technologies IDT, Coralville, IA 52241, United States
| | - Stefania Evoli
- Department of Physics and The Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, IL, United States
| | - Nicholas J Schnicker
- Protein & Crystallography Facility, Carver College of Medicine, University of Iowa, Iowa City, IA 52242, United States
| | - Maria Nunez-Hernandez
- Department of Biochemistry, Carver College of Medicine, University of Iowa, Iowa City, IA 52242, United States
| | - Liping Yu
- Department of Biochemistry, Carver College of Medicine NMR Core Facility, University of Iowa, Iowa City, IA 52242, United States; The Iowa City Veterans Affairs Medical Center, Iowa City, IA 52242, United States
| | - Jeff Wereszczynski
- Department of Physics and The Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, IL, United States.
| | - Miles A Pufall
- Department of Biochemistry, Carver College of Medicine, University of Iowa, Iowa City, IA 52242, United States.
| | - Catherine A Musselman
- Department of Biochemistry, Carver College of Medicine, University of Iowa, Iowa City, IA 52242, United States; Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, CO 80045, United States.
| |
Collapse
|
38
|
Alexiou TS, Alatas PV, Tsalikis DG, Mavrantzas VG. Conformational and Dynamic Properties of Short DNA Minicircles in Aqueous Solution from Atomistic Molecular Dynamics Simulations. Macromolecules 2020. [DOI: 10.1021/acs.macromol.0c00821] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [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
| | - Panagiotis V. Alatas
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
| | - Dimitrios G. Tsalikis
- 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
|
39
|
Craggs TD, Sustarsic M, Plochowietz A, Mosayebi M, Kaju H, Cuthbert A, Hohlbein J, Domicevica L, Biggin PC, Doye JPK, Kapanidis AN. Substrate conformational dynamics facilitate structure-specific recognition of gapped DNA by DNA polymerase. Nucleic Acids Res 2020; 47:10788-10800. [PMID: 31544938 PMCID: PMC6846080 DOI: 10.1093/nar/gkz797] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2019] [Revised: 09/02/2019] [Accepted: 09/18/2019] [Indexed: 01/23/2023] Open
Abstract
DNA-binding proteins utilise different recognition mechanisms to locate their DNA targets; some proteins recognise specific DNA sequences, while others interact with specific DNA structures. While sequence-specific DNA binding has been studied extensively, structure-specific recognition mechanisms remain unclear. Here, we study structure-specific DNA recognition by examining the structure and dynamics of DNA polymerase I Klenow Fragment (Pol) substrates both alone and in DNA–Pol complexes. Using a docking approach based on a network of 73 distances collected using single-molecule FRET, we determined a novel solution structure of the single-nucleotide-gapped DNA–Pol binary complex. The structure resembled existing crystal structures with regards to the downstream primer-template DNA substrate, and revealed a previously unobserved sharp bend (∼120°) in the DNA substrate; this pronounced bend was present in living cells. MD simulations and single-molecule assays also revealed that 4–5 nt of downstream gap-proximal DNA are unwound in the binary complex. Further, experiments and coarse-grained modelling showed the substrate alone frequently adopts bent conformations with 1–2 nt fraying around the gap, suggesting a mechanism wherein Pol recognises a pre-bent, partially-melted conformation of gapped DNA. We propose a general mechanism for substrate recognition by structure-specific enzymes driven by protein sensing of the conformational dynamics of their DNA substrates.
Collapse
Affiliation(s)
- Timothy D Craggs
- Biological Physics Research Group, Clarendon Laboratory, Department of Physics, University of Oxford, Oxford OX1 3PU, UK
| | - Marko Sustarsic
- Biological Physics Research Group, Clarendon Laboratory, Department of Physics, University of Oxford, Oxford OX1 3PU, UK
| | - Anne Plochowietz
- Biological Physics Research Group, Clarendon Laboratory, Department of Physics, University of Oxford, Oxford OX1 3PU, UK
| | - Majid Mosayebi
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, UK.,School of Physics, Institute for Research in Fundamental Sciences (IPM), Tehran 19538-33511, Iran
| | - Hendrik Kaju
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, UK
| | - Andrew Cuthbert
- Biological Physics Research Group, Clarendon Laboratory, Department of Physics, University of Oxford, Oxford OX1 3PU, UK
| | - Johannes Hohlbein
- Laboratory of Biophysics, Wageningen University & Research, Wageningen 6708 WE, The Netherlands.,Microspectroscopy Research Facility Wageningen, Wageningen University & Research, Wageningen 6708 WE, The Netherlands
| | - Laura Domicevica
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, UK
| | - Philip C Biggin
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, UK
| | - Jonathan P K Doye
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, UK
| | - Achillefs N Kapanidis
- Biological Physics Research Group, Clarendon Laboratory, Department of Physics, University of Oxford, Oxford OX1 3PU, UK
| |
Collapse
|
40
|
Dohnalová H, Dršata T, Šponer J, Zacharias M, Lipfert J, Lankaš F. Compensatory Mechanisms in Temperature Dependence of DNA Double Helical Structure: Bending and Elongation. J Chem Theory Comput 2020; 16:2857-2863. [DOI: 10.1021/acs.jctc.0c00037] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hana Dohnalová
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| | - Tomáš Dršata
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Martin Zacharias
- Physics-Department T38, Technical University of Munich, James-Franck-Strasse 1, 85748 Garching, Germany
| | - Jan Lipfert
- Department of Physics and Center for Nanoscience, LMU Munich, Amalienstrasse 54, 80799 Munich, Germany
| | - Filip Lankaš
- Department of Informatics and Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28 Prague, Czech Republic
| |
Collapse
|
41
|
Walther J, Dans PD, Balaceanu A, Hospital A, Bayarri G, Orozco M. A multi-modal coarse grained model of DNA flexibility mappable to the atomistic level. Nucleic Acids Res 2020; 48:e29. [PMID: 31956910 PMCID: PMC7049737 DOI: 10.1093/nar/gkaa015] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2019] [Revised: 12/10/2019] [Accepted: 01/04/2020] [Indexed: 12/13/2022] Open
Abstract
We present a new coarse grained method for the simulation of duplex DNA. The algorithm uses a generalized multi-harmonic model that can represent any multi-normal distribution of helical parameters, thus avoiding caveats of current mesoscopic models for DNA simulation and representing a breakthrough in the field. The method has been parameterized from accurate parmbsc1 atomistic molecular dynamics simulations of all unique tetranucleotide sequences of DNA embedded in long duplexes and takes advantage of the correlation between helical states and backbone configurations to derive atomistic representations of DNA. The algorithm, which is implemented in a simple web interface and in a standalone package reproduces with high computational efficiency the structural landscape of long segments of DNA untreatable by atomistic molecular dynamics simulations.
Collapse
Affiliation(s)
- Jürgen Walther
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Pablo D Dans
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Department of Biological Sciences, CENUR North Coast, University of the Republic, 50000 Salto, Uruguay
| | - Alexandra Balaceanu
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Genís Bayarri
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Department of Biochemistry and Biomedicine, The University of Barcelona, 08028 Barcelona, Spain
| |
Collapse
|
42
|
Bendandi A, Dante S, Zia SR, Diaspro A, Rocchia W. Chromatin Compaction Multiscale Modeling: A Complex Synergy Between Theory, Simulation, and Experiment. Front Mol Biosci 2020; 7:15. [PMID: 32158765 PMCID: PMC7051991 DOI: 10.3389/fmolb.2020.00015] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2019] [Accepted: 01/27/2020] [Indexed: 12/24/2022] Open
Abstract
Understanding the mechanisms that trigger chromatin compaction, its patterns, and the factors they depend on, is a fundamental and still open question in Biology. Chromatin compacts and reinforces DNA and is a stable but dynamic structure, to make DNA accessible to proteins. In recent years, computational advances have provided larger amounts of data and have made large-scale simulations more viable. Experimental techniques for the extraction and reconstitution of chromatin fibers have improved, reinvigorating theoretical and experimental interest in the topic and stimulating debate on points previously considered as certainties regarding chromatin. A great assortment of approaches has emerged, from all-atom single-nucleosome or oligonucleosome simulations to various degrees of coarse graining, to polymer models, to fractal-like structures and purely topological models. Different fiber-start patterns have been studied in theory and experiment, as well as different linker DNA lengths. DNA is a highly charged macromolecule, making ionic and electrostatic interactions extremely important for chromatin topology and dynamics. Indeed, the repercussions of varying ionic concentration have been extensively examined at the computational level, using all-atom, coarse-grained, and continuum techniques. The presence of high-curvature AT-rich segments in DNA can cause conformational variations, attesting to the fact that the role of DNA is both structural and electrostatic. There have been some tentative attempts to describe the force fields governing chromatin conformational changes and the energy landscapes of these transitions, but the intricacy of the system has hampered reaching a consensus. The study of chromatin conformations is an intrinsically multiscale topic, influenced by a wide range of biological and physical interactions, spanning from the atomic to the chromosome level. Therefore, powerful modeling techniques and carefully planned experiments are required for an overview of the most relevant phenomena and interactions. The topic provides fertile ground for interdisciplinary studies featuring a synergy between theoretical and experimental scientists from different fields and the cross-validation of respective results, with a multi-scale perspective. Here, we summarize some of the most representative approaches, and focus on the importance of electrostatics and solvation, often overlooked aspects of chromatin modeling.
Collapse
Affiliation(s)
- Artemi Bendandi
- Physics Department, University of Genoa, Genoa, Italy
- Nanophysics & NIC@IIT, Istituto Italiano di Tecnologia, Genoa, Italy
| | - Silvia Dante
- Nanophysics & NIC@IIT, Istituto Italiano di Tecnologia, Genoa, Italy
| | - Syeda Rehana Zia
- Dr. Panjwani Center for Molecular Medicine and Drug Research, International Center for Chemical and Biological Sciences, University of Karachi, Karachi, Pakistan
| | - Alberto Diaspro
- Physics Department, University of Genoa, Genoa, Italy
- Nanophysics & NIC@IIT, Istituto Italiano di Tecnologia, Genoa, Italy
| | - Walter Rocchia
- Concept Lab, Istituto Italiano di Tecnologia, Genoa, Italy
| |
Collapse
|
43
|
Donati E, Genna V, De Vivo M. Recruiting Mechanism and Functional Role of a Third Metal Ion in the Enzymatic Activity of 5' Structure-Specific Nucleases. J Am Chem Soc 2020; 142:2823-2834. [PMID: 31939291 PMCID: PMC7993637 DOI: 10.1021/jacs.9b10656] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
![]()
Enzymes of the 5′ structure-specific
nuclease family are crucial for DNA repair, replication, and recombination.
One such enzyme is the human exonuclease 1 (hExo1) metalloenzyme,
which cleaves DNA strands, acting primarily as a processive 5′-3′
exonuclease and secondarily as a 5′-flap endonuclease. Recently,
in crystallo reaction intermediates have elucidated how hExo1 exerts
hydrolysis of DNA phosphodiester bonds. These hExo1 structures show
a third metal ion intermittently bound close to the two-metal-ion
active site, to which recessed ends or 5′-flap substrates bind.
Evidence of this third ion has been observed in several nucleic-acid-processing
metalloenzymes. However, there is still debate over what triggers
the (un)binding of this transient third ion during catalysis and whether
this ion has a catalytic function. Using extended molecular dynamics
and enhanced sampling free-energy simulations, we observed that the
carboxyl side chain of Glu89 (located along the arch motif in hExo1)
flips frequently from the reactant state to the product state. The
conformational flipping of Glu89 allows one metal ion to be recruited
from the bulk and promptly positioned near the catalytic center. This
is in line with the structural evidence. Additionally, our simulations
show that the third metal ion assists the departure, through the mobile
arch, of the nucleotide monophosphate product from the catalytic site.
Structural comparisons of nuclease enzymes suggest that this Glu(Asp)-mediated
mechanism for third ion recruitment and nucleic acid hydrolysis may
be shared by other 5′ structure-specific nucleases.
Collapse
Affiliation(s)
- Elisa Donati
- Laboratory of Molecular Modelling & Drug Discovery , Istituto Italiano di Tecnologia , Via Morego 30 , 16163 Genoa , Italy
| | - Vito Genna
- Laboratory of Molecular Modelling & Drug Discovery , Istituto Italiano di Tecnologia , Via Morego 30 , 16163 Genoa , Italy
| | - Marco De Vivo
- Laboratory of Molecular Modelling & Drug Discovery , Istituto Italiano di Tecnologia , Via Morego 30 , 16163 Genoa , Italy
| |
Collapse
|
44
|
Heinz M, Erlenbach N, Stelzl LS, Thierolf G, Kamble NR, Sigurdsson ST, Prisner TF, Hummer G. High-resolution EPR distance measurements on RNA and DNA with the non-covalent Ǵ spin label. Nucleic Acids Res 2020; 48:924-933. [PMID: 31777925 PMCID: PMC6954412 DOI: 10.1093/nar/gkz1096] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2019] [Revised: 11/01/2019] [Accepted: 11/20/2019] [Indexed: 12/25/2022] Open
Abstract
Pulsed electron paramagnetic resonance (EPR) experiments, among them most prominently pulsed electron-electron double resonance experiments (PELDOR/DEER), resolve the conformational dynamics of nucleic acids with high resolution. The wide application of these powerful experiments is limited by the synthetic complexity of some of the best-performing spin labels. The recently developed $\bf\acute{G}$ (G-spin) label, an isoindoline-nitroxide derivative of guanine, can be incorporated non-covalently into DNA and RNA duplexes via Watson-Crick base pairing in an abasic site. We used PELDOR and molecular dynamics (MD) simulations to characterize $\bf\acute{G}$, obtaining excellent agreement between experiments and time traces calculated from MD simulations of RNA and DNA double helices with explicitly modeled $\bf\acute{G}$ bound in two abasic sites. The MD simulations reveal stable hydrogen bonds between the spin labels and the paired cytosines. The abasic sites do not significantly perturb the helical structure. $\bf\acute{G}$ remains rigidly bound to helical RNA and DNA. The distance distributions between the two bound $\bf\acute{G}$ labels are not substantially broadened by spin-label motions in the abasic site and agree well between experiment and MD. $\bf\acute{G}$ and similar non-covalently attached spin labels promise high-quality distance and orientation information, also of complexes of nucleic acids and proteins.
Collapse
Affiliation(s)
- Marcel Heinz
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
| | - Nicole Erlenbach
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt, Max-von-Laue-Straße 7, 60438 Frankfurt am Main, Germany
| | - Lukas S Stelzl
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
| | - Grace Thierolf
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt, Max-von-Laue-Straße 7, 60438 Frankfurt am Main, Germany
| | - Nilesh R Kamble
- Department of Chemistry, Science Institute, University of Iceland, Dunhaga 3, 107 Reykjavk, Iceland
| | - Snorri Th Sigurdsson
- Department of Chemistry, Science Institute, University of Iceland, Dunhaga 3, 107 Reykjavk, Iceland
| | - Thomas F Prisner
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University Frankfurt, Max-von-Laue-Straße 7, 60438 Frankfurt am Main, Germany
| | - Gerhard Hummer
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
- Institute for Biophysics, Goethe University Frankfurt, 60438 Frankfurt am Main, Germany
| |
Collapse
|
45
|
Babbitt GA, Fokoue EP, Evans JR, Diller KI, Adams LE. DROIDS 3.0-Detecting Genetic and Drug Class Variant Impact on Conserved Protein Binding Dynamics. Biophys J 2019; 118:541-551. [PMID: 31928763 PMCID: PMC7002913 DOI: 10.1016/j.bpj.2019.12.008] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2019] [Revised: 12/09/2019] [Accepted: 12/10/2019] [Indexed: 01/07/2023] Open
Abstract
The application of statistical methods to comparatively framed questions about the molecular dynamics (MD) of proteins can potentially enable investigations of biomolecular function beyond the current sequence and structural methods in bioinformatics. However, the chaotic behavior in single MD trajectories requires statistical inference that is derived from large ensembles of simulations representing the comparative functional states of a protein under investigation. Meaningful interpretation of such complex forms of big data poses serious challenges to users of MD. Here, we announce Detecting Relative Outlier Impacts from Molecular Dynamic Simulation (DROIDS) 3.0, a method and software package for comparative protein dynamics that includes maxDemon 1.0, a multimethod machine learning application that trains on large ensemble comparisons of concerted protein motions in opposing functional states generated by DROIDS and deploys learned classifications of these states onto newly generated MD simulations. Local canonical correlations in learning patterns generated from independent, yet identically prepared, MD validation runs are used to identify regions of functionally conserved protein dynamics. The subsequent impacts of genetic and/or drug class variants on conserved dynamics can also be analyzed by deploying the classifiers on variant MD simulations and quantifying how often these altered protein systems display opposing functional states. Here, we present several case studies of complex changes in functional protein dynamics caused by temperature, genetic mutation, and binding interactions with nucleic acids and small molecules. We demonstrate that our machine learning algorithm can properly identify regions of functionally conserved dynamics in ubiquitin and TATA-binding protein (TBP). We quantify the impact of genetic variation in TBP and drug class variation targeting the ATP-binding region of Hsp90 on conserved dynamics. We identify regions of conserved dynamics in Hsp90 that connect the ATP binding pocket to other functional regions. We also demonstrate that dynamic impacts of various Hsp90 inhibitors rank accordingly with how closely they mimic natural ATP binding.
Collapse
Affiliation(s)
- Gregory A Babbitt
- Thomas H. Gosnell School of Life Sciences, Rochester Institute of Technology, Rochester, New York.
| | - Ernest P Fokoue
- School of Mathematical Sciences, Rochester Institute of Technology, Rochester, New York
| | - Joshua R Evans
- Thomas H. Gosnell School of Life Sciences, Rochester Institute of Technology, Rochester, New York
| | - Kyle I Diller
- Thomas H. Gosnell School of Life Sciences, Rochester Institute of Technology, Rochester, New York; Golisano College for Computing and Information Science, Rochester, New York
| | - Lily E Adams
- Thomas H. Gosnell School of Life Sciences, Rochester Institute of Technology, Rochester, New York
| |
Collapse
|
46
|
Minhas V, Sun T, Mirzoev A, Korolev N, Lyubartsev AP, Nordenskiöld L. Modeling DNA Flexibility: Comparison of Force Fields from Atomistic to Multiscale Levels. J Phys Chem B 2019; 124:38-49. [DOI: 10.1021/acs.jpcb.9b09106] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Vishal Minhas
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Tiedong Sun
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Alexander Mirzoev
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Nikolay Korolev
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Alexander P. Lyubartsev
- Department of Materials and Environmental Chemistry, Stockholm University, 10691 Stockholm, Sweden
| | - Lars Nordenskiöld
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| |
Collapse
|
47
|
Dans PD, Balaceanu A, Pasi M, Patelli AS, Petkevičiūtė D, Walther J, Hospital A, Bayarri G, Lavery R, Maddocks JH, Orozco M. The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules. Nucleic Acids Res 2019; 47:11090-11102. [PMID: 31624840 PMCID: PMC6868377 DOI: 10.1093/nar/gkz905] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2019] [Revised: 09/25/2019] [Accepted: 10/06/2019] [Indexed: 12/12/2022] Open
Abstract
We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence. We characterize the sequence-dependent choreography of backbone and base movements modulating the non-Gaussian or anharmonic effects manifested in the higher moments of the dynamics of the duplex when sampling the equilibrium distribution. Contrary to prior assumptions, such anharmonic deformations are not rare in DNA and can play a significant role in determining DNA conformation within complexes. Polymorphisms in helical geometries are particularly prevalent for certain tetranucleotide sequence contexts and are always coupled to a complex network of coordinated changes in the backbone. The analysis of our simulations, which contain instances of all tetranucleotide sequences, allow us to extend Calladine-Dickerson rules used for decades to interpret the average geometry of DNA, leading to a set of rules with quantitative predictive power that encompass nonlocal sequence-dependence and anharmonic fluctuations.
Collapse
Affiliation(s)
- Pablo D Dans
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
- Department of Biological Sciences, University of the Republic (UdelaR), CENUR Gral. Rivera 1350, 50000 Salto, Uruguay
| | - Alexandra Balaceanu
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
| | - Marco Pasi
- LBPA, École normale supérieure Paris-Saclay, 61 Av. du Pdt Wilson, Cachan 94235, France
- Bases Moléculaires et Structurales des Systèmes Infectieux, Univ. Lyon I/CNRS UMR 5086, IBCP, 7 Passage du Vercors, Lyon 69367, France
| | - Alessandro S Patelli
- Institute of Mathematics, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Daiva Petkevičiūtė
- Institute of Mathematics, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
- Faculty of Mathematics and Natural Sciences, Kaunas University of Technology, Studentų g. 50, 51368 Kaunas, Lithuania
| | - Jürgen Walther
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
| | - Genís Bayarri
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
| | - Richard Lavery
- Bases Moléculaires et Structurales des Systèmes Infectieux, Univ. Lyon I/CNRS UMR 5086, IBCP, 7 Passage du Vercors, Lyon 69367, France
| | - John H Maddocks
- Institute of Mathematics, Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10–12, 08028 Barcelona, Spain
- Department of Biochemistry and Molecular Biology. University of Barcelona, 08028 Barcelona, Spain
| |
Collapse
|
48
|
Hospital A, Battistini F, Soliva R, Gelpí JL, Orozco M. Surviving the deluge of biosimulation data. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1449] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Adam Hospital
- Institut de Recerca Biomèdica, IRB Barcelona, The Barcelona Institute of Science and Technology Joint IRB‐BSC Program in Computational Biology Barcelona Spain
| | - Federica Battistini
- Institut de Recerca Biomèdica, IRB Barcelona, The Barcelona Institute of Science and Technology Joint IRB‐BSC Program in Computational Biology Barcelona Spain
| | | | - Josep Lluis Gelpí
- Barcelona Supercomputing Center Join IRB‐BSC Program in Computational Biology Barcelona Spain
- Departament de Bioquímica i Biomedicina Universitat de Barcelona Barcelona Spain
| | - Modesto Orozco
- Institut de Recerca Biomèdica, IRB Barcelona, The Barcelona Institute of Science and Technology Joint IRB‐BSC Program in Computational Biology Barcelona Spain
- Departament de Bioquímica i Biomedicina Universitat de Barcelona Barcelona Spain
| |
Collapse
|
49
|
Balaceanu A, Buitrago D, Walther J, Hospital A, Dans PD, Orozco M. Modulation of the helical properties of DNA: next-to-nearest neighbour effects and beyond. Nucleic Acids Res 2019; 47:4418-4430. [PMID: 30957854 PMCID: PMC6511876 DOI: 10.1093/nar/gkz255] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2019] [Revised: 03/23/2019] [Accepted: 03/30/2019] [Indexed: 12/21/2022] Open
Abstract
We used extensive molecular dynamics simulations to study the structural and dynamic properties of the central d(TpA) step in the highly polymorphic d(CpTpApG) tetranucleotide. Contrary to the assumption of the dinucleotide-model and its nearest neighbours (tetranucleotide-model), the properties of the central d(TpA) step change quite significantly dependent on the next-to-nearest (hexanucleotide) sequence context and in a few cases are modulated by even remote neighbours (beyond next-to-nearest from the central TpA). Our results highlight the existence of previously undescribed dynamical mechanisms for the transmission of structural information into the DNA and demonstrate the existence of certain sequences with special physical properties that can impact on the global DNA structure and dynamics.
Collapse
Affiliation(s)
- Alexandra Balaceanu
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028 Barcelona, Spain
| | - Diana Buitrago
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028 Barcelona, Spain
| | - Jürgen Walther
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028 Barcelona, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028 Barcelona, Spain
| | - Pablo D Dans
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), 08028 Barcelona, Spain
| | - Modesto Orozco
- 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
|
50
|
Francés-Monerris A, Tuñón I, Monari A. Hypoxia-Selective Dissociation Mechanism of a Nitroimidazole Nucleoside in a DNA Environment. J Phys Chem Lett 2019; 10:6750-6754. [PMID: 31609626 DOI: 10.1021/acs.jpclett.9b02760] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Photodynamic therapy is a promising approach to treat a variety of superficial tumors and other diseases. One of its major limitations arises from its dependence on molecular oxygen, which decreases the efficiency of the therapy in hypoxia conditions commonly developed by solid tumors. The present contribution reveals the molecular mechanism of a modified thymine bearing a nitroimidazole substituent, a photosensitizer able to produce highly harmful interstrand cross-links in the DNA double strand after irradiation selectively in absence of oxygen. The mechanism is resolved at a fully atomistic and electronic level relying on quantum mechanics (CASPT2, coupled-cluster, DFT, and TD-DFT methods), classical molecular dynamics, and advanced biased QM/MM simulations, revealing an energy penalty of ∼8 kcal/mol for the anionic nitromidazole release. Our findings indicate that the global interstrand cross-link production is driven by a combination of multiple factors, namely, the reverse energy penalty, the diffusion of the nitroimidazole anion, and the further reactivity of the formed thymine radical. On the basis of these results, we also suggest some possible strategies to improve the efficiency of interstrand cross-link production.
Collapse
Affiliation(s)
- Antonio Francés-Monerris
- Université de Lorraine, CNRS, LPCT UMR 7019 , F-54000 Nancy , France
- Departamento de Química Física , Universitat de València , 46100 Burjassot , Spain
| | - Iñaki Tuñón
- Departamento de Química Física , Universitat de València , 46100 Burjassot , Spain
| | - Antonio Monari
- Université de Lorraine, CNRS, LPCT UMR 7019 , F-54000 Nancy , France
| |
Collapse
|