1
|
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
|
2
|
Rozza AM, Bakó I, Oláh J. Theoretical insights into water network of B-DNA duplex with Watson-Crick and Hoogsteen base pairing geometries. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119728] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
3
|
Yousefi R, Lynch GC, Galbraith M, Pettitt BM. Contributions of higher-order proximal distribution functions to solvent structure around proteins. J Chem Phys 2021; 155:104110. [PMID: 34525817 PMCID: PMC8439718 DOI: 10.1063/5.0062580] [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: 07/06/2021] [Accepted: 08/19/2021] [Indexed: 11/14/2022] Open
Abstract
The proximal distribution function (pDF) quantifies the probability of finding a solvent molecule in the vicinity of solutes. The approach constitutes a hierarchically organized theory for constructing approximate solvation structures around solutes. Given the assumption of universality of atom cluster-specific solvation, reconstruction of the solvent distribution around arbitrary molecules provides a computationally convenient route to solvation thermodynamics. Previously, such solvent reconstructions usually considered the contribution of the nearest-neighbor distribution only. We extend the pDF reconstruction algorithm to terms including next-nearest-neighbor contribution. As a test, small molecules (alanine and butane) are examined. The analysis is then extended to include the protein myoglobin in the P6 crystal unit cell. Molecular dynamics simulations are performed, and solvent density distributions around the solute molecules are compared with the results from different pDF reconstruction models. It is shown that the next-nearest-neighbor modification significantly improves the reconstruction of the solvent number density distribution in concave regions and between solute molecules. The probability densities are then used to calculate the solute-solvent non-bonded interaction energies including van der Waals and electrostatic, which are found to be in good agreement with the simulated values.
Collapse
Affiliation(s)
- Razie Yousefi
- Department of Biochemistry and Molecular Biology, University of Texas Medical Branch, 301 University Boulevard, Galveston, Texas 77555-0304, USA
| | - Gillian C. Lynch
- Department of Biochemistry and Molecular Biology, University of Texas Medical Branch, 301 University Boulevard, Galveston, Texas 77555-0304, USA
| | - Madeline Galbraith
- Department of Biochemistry and Molecular Biology, University of Texas Medical Branch, 301 University Boulevard, Galveston, Texas 77555-0304, USA
| | - B. Montgomery Pettitt
- Department of Biochemistry and Molecular Biology, University of Texas Medical Branch, 301 University Boulevard, Galveston, Texas 77555-0304, USA
| |
Collapse
|
4
|
Fick RJ, Liu AY, Nussbaumer F, Kreutz C, Rangadurai A, Xu Y, Sommer RD, Shi H, Scheiner S, Stelling AL. Probing the Hydrogen-Bonding Environment of Individual Bases in DNA Duplexes with Isotope-Edited Infrared Spectroscopy. J Phys Chem B 2021; 125:7613-7627. [PMID: 34236202 PMCID: PMC8311644 DOI: 10.1021/acs.jpcb.1c01351] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
![]()
Measuring the strength
of the hydrogen bonds between DNA base pairs
is of vital importance for understanding how our genetic code is physically
accessed and recognized in cells, particularly during replication
and transcription. Therefore, it is important to develop probes for
these key hydrogen bonds (H-bonds) that dictate events critical to
cellular function, such as the localized melting of DNA. The vibrations
of carbonyl bonds are well-known probes of their H-bonding environment,
and their signals can be observed with infrared (IR) spectroscopy.
Yet, pinpointing a single bond of interest in the complex IR spectrum
of DNA is challenging due to the large number of carbonyl signals
that overlap with each other. Here, we develop a method using isotope
editing and infrared (IR) spectroscopy to isolate IR signals from
the thymine (T) C2=O carbonyl. We use solvatochromatic studies
to show that the TC2=O signal’s position in the IR spectrum
is sensitive to the H-bonding capacity of the solvent. Our results
indicate that C2=O of a single T base within DNA duplexes experiences
weak H-bonding interactions. This finding is consistent with the existence
of a third, noncanonical CH···O H-bond between adenine
and thymine in both Watson–Crick and Hoogsteen base pairs in
DNA.
Collapse
Affiliation(s)
- Robert J Fick
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Amy Y Liu
- Department of Biochemistry, Duke University Medical Center, Durham, North Carolina 27710, United States
| | - Felix Nussbaumer
- Institute of Organic Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck 6020, Austria
| | - Christoph Kreutz
- Institute of Organic Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck 6020, Austria
| | - Atul Rangadurai
- Department of Biochemistry, Duke University Medical Center, Durham, North Carolina 27710, United States
| | - Yu Xu
- Department of Chemistry, Duke University, Durham, North Carolina 27710, United States
| | - Roger D Sommer
- Molecular Education, Technology, and Research Innovation Center, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Honglue Shi
- Department of Chemistry, Duke University, Durham, North Carolina 27710, United States
| | - Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322, United States
| | - Allison L Stelling
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, Richardson, Texas 75080, United States.,Department of Biochemistry, Duke University Medical Center, Durham, North Carolina 27710, United States
| |
Collapse
|
5
|
Yonetani Y, Nakagawa H. Understanding water-mediated DNA damage production by molecular dynamics calculation of solvent accessibility. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137441] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
6
|
Yu B, Montgomery Pettitt B, Iwahara J. Experimental Evidence of Solvent-Separated Ion Pairs as Metastable States in Electrostatic Interactions of Biological Macromolecules. J Phys Chem Lett 2019; 10:7937-7941. [PMID: 31809050 PMCID: PMC6936746 DOI: 10.1021/acs.jpclett.9b03084] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Electrostatic interactions via ion pairs are vital for biological macromolecules. Regarding the free energy of each ion pair as a function of the interionic distance, continuum electrostatic models predict a single energy minimum corresponding to the contact ion-pair (CIP) state, whereas atomically detailed theoretical hydration studies predict multiple energy minima corresponding to the CIP and solvent-separated ion-pair (SIP) states. Through a statistical analysis of high-resolution crystal structures, we present experimental evidence of the SIP as a metastable state. The histogram of interionic distances between protein side-chain NH3+ and DNA phosphate groups clearly shows two major peaks corresponding to the CIP and SIP states. The statistical data are consistent with the probability distribution of the CIP-SIP equilibria previously obtained with molecular dynamics simulations. Spatial distributions of NH3+ ions and water molecules around phosphates reveal preferential sites for CIP and SIP formations and show how the ions compete with water molecules.
Collapse
|
7
|
Zhang H, Fu H, Shao X, Dehez F, Chipot C, Cai W. Changes in Microenvironment Modulate the B- to A-DNA Transition. J Chem Inf Model 2019; 59:2324-2330. [PMID: 30767527 DOI: 10.1021/acs.jcim.8b00885] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
B- to A-DNA transition is known to be sensitive to the macroscopic properties of the solution, such as salt and ethanol concentrations. Microenvironmental effects on DNA conformational transition have been broadly studied. Providing an intuitive picture of how DNA responds to environmental changes is, however, still needed. Analyzing the chemical equilibrium of B-to-A DNA transition at critical concentrations, employing explicit-solvent simulations, is envisioned to help understand such microenvironmental effects. In the present study, free-energy calculations characterizing the B- to A-DNA transition and the distribution of cations were carried out in solvents with different ethanol concentrations. With the addition of ethanol, the most stable structure of DNA changes from the B- to A-form, in agreement with previous experimental observation. In 60% ethanol, a chemical equilibrium is found, showing reversible transition between B- and A-DNA. Analysis of the microenvironment around DNA suggests that with the increase of ethanol concentration, the cations exhibit a significant tendency to move toward the backbone, and mobility of water molecules around the major groove and backbone decreases gradually, leading eventually to a B-to-A transition. The present results provide a free-energy view of DNA microenvironment and of the role of cation motion in the conformational transition.
Collapse
Affiliation(s)
- Hong Zhang
- Research Center for Analytical Sciences, College of Chemistry , Nankai University, Tianjin Key Laboratory of Biosensing and Molecular Recognition , Tianjin 300071 , China
| | - Haohao Fu
- Research Center for Analytical Sciences, College of Chemistry , Nankai University, Tianjin Key Laboratory of Biosensing and Molecular Recognition , Tianjin 300071 , China
| | - Xueguang Shao
- Research Center for Analytical Sciences, College of Chemistry , Nankai University, Tianjin Key Laboratory of Biosensing and Molecular Recognition , Tianjin 300071 , China.,State Key Laboratory of Medicinal Chemical Biology , Tianjin 300071 , China.,Collaborative Innovation Center of Chemical Science and Engineering , Tianjin 300071 , China
| | - François Dehez
- Laboratoire International Associé CNRS and University of Illinois at Urbana-Champaign , Vandoeuvre-lès-Nancy F-54506 , France.,LPCT, UMR 7019 Université de Lorraine CNRS , Vandoeuvre-lès-Nancy F-54500 , France
| | - Christophe Chipot
- Laboratoire International Associé CNRS and University of Illinois at Urbana-Champaign , Vandoeuvre-lès-Nancy F-54506 , France.,LPCT, UMR 7019 Université de Lorraine CNRS , Vandoeuvre-lès-Nancy F-54500 , France.,Department of Physics , University of Illinois at Urbana-Champaign , Urbana , Illinois 61801 , United States
| | - Wensheng Cai
- Research Center for Analytical Sciences, College of Chemistry , Nankai University, Tianjin Key Laboratory of Biosensing and Molecular Recognition , Tianjin 300071 , China.,Collaborative Innovation Center of Chemical Science and Engineering , Tianjin 300071 , China
| |
Collapse
|
8
|
Abstract
The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.
Collapse
Affiliation(s)
- Damien Laage
- École
Normale Supérieure, PSL Research University, UPMC Univ Paris
06, CNRS, Département de Chimie,
PASTEUR, 24 rue Lhomond, 75005 Paris, France
- Sorbonne
Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR, 75005 Paris, France
| | - Thomas Elsaesser
- Max-Born-Institut
für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| | - James T. Hynes
- École
Normale Supérieure, PSL Research University, UPMC Univ Paris
06, CNRS, Département de Chimie,
PASTEUR, 24 rue Lhomond, 75005 Paris, France
- Sorbonne
Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR, 75005 Paris, France
- Department
of Chemistry and Biochemistry, University
of Colorado, Boulder, Colorado 80309, United
States
| |
Collapse
|
9
|
Liu Y, Guchhait B, Siebert T, Fingerhut BP, Elsaesser T. Molecular couplings and energy exchange between DNA and water mapped by femtosecond infrared spectroscopy of backbone vibrations. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2017; 4:044015. [PMID: 28405593 PMCID: PMC5384856 DOI: 10.1063/1.4980075] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/20/2017] [Accepted: 03/31/2017] [Indexed: 05/15/2023]
Abstract
Molecular couplings between DNA and water together with the accompanying processes of energy exchange are mapped via the ultrafast response of DNA backbone vibrations after OH stretch excitation of the water shell. Native salmon testes DNA is studied in femtosecond pump-probe experiments under conditions of full hydration and at a reduced hydration level with two water layers around the double helix. Independent of their local hydration patterns, all backbone vibrations in the frequency range from 940 to 1120 cm-1 display a quasi-instantaneous reshaping of the spectral envelopes of their fundamental absorption bands upon excitation of the water shell. The subsequent reshaping kinetics encompass a one-picosecond component, reflecting the formation of a hot ground state of the water shell, and a slower contribution on a time scale of tens of picoseconds. Such results are benchmarked by measurements with resonant excitation of the backbone modes, resulting in distinctly different absorption changes. We assign the fast changes of DNA absorption after OH stretch excitation to structural changes in the water shell which couple to DNA through the local electric fields. The second slower process is attributed to a flow of excess energy from the water shell into DNA, establishing a common heated ground state in the molecular ensemble. This interpretation is supported by theoretical calculations of the electric fields exerted by the water shell at different temperatures.
Collapse
Affiliation(s)
- Yingliang Liu
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , 12489 Berlin, Germany
| | - Biswajit Guchhait
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , 12489 Berlin, Germany
| | - Torsten Siebert
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , 12489 Berlin, Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , 12489 Berlin, Germany
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , 12489 Berlin, Germany
| |
Collapse
|
10
|
Nakano M, Tateishi-Karimata H, Tanaka S, Tama F, Miyashita O, Nakano SI, Sugimoto N. Local thermodynamics of the water molecules around single- and double-stranded DNA studied by grid inhomogeneous solvation theory. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.08.032] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
11
|
Siebert T, Guchhait B, Liu Y, Fingerhut BP, Elsaesser T. Range, Magnitude, and Ultrafast Dynamics of Electric Fields at the Hydrated DNA Surface. J Phys Chem Lett 2016; 7:3131-6. [PMID: 27468144 DOI: 10.1021/acs.jpclett.6b01369] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Range and magnitude of electric fields at biomolecular interfaces and their fluctuations in a time window down to the subpicosecond regime have remained controversial, calling for electric-field mapping in space and time. Here, we trace fluctuating electric fields at the surface of native salmon DNA via their interactions with backbone vibrations in a wide range of hydration levels by building the water shell layer by layer. Femtosecond two-dimensional infrared spectroscopy and ab initio based theory establish water molecules in the first two layers as the predominant source of interfacial electric fields, which fluctuate on a 300 fs time scale with an amplitude of 25 MV/cm due to thermally excited water motions. The observed subnanometer range of these electric interactions is decisive for biochemical structure and function.
Collapse
Affiliation(s)
- Torsten Siebert
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , D-12489 Berlin, Germany
| | - Biswajit Guchhait
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , D-12489 Berlin, Germany
| | - Yingliang Liu
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , D-12489 Berlin, Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , D-12489 Berlin, Germany
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , D-12489 Berlin, Germany
| |
Collapse
|
12
|
Guchhait B, Liu Y, Siebert T, Elsaesser T. Ultrafast vibrational dynamics of the DNA backbone at different hydration levels mapped by two-dimensional infrared spectroscopy. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2016; 3:043202. [PMID: 26798841 PMCID: PMC4720115 DOI: 10.1063/1.4936567] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/01/2015] [Accepted: 11/12/2015] [Indexed: 05/28/2023]
Abstract
DNA oligomers are studied at 0% and 92% relative humidity, corresponding to N < 2 and N > 20 water molecules per base pair. Two-dimensional (2D) infrared spectroscopy of DNA backbone modes between 920 and 1120 cm(-1) maps fluctuating interactions at the DNA surface. At both hydration levels, a frequency fluctuation correlation function with a 300 fs decay and a slow decay beyond 10 ps is derived from the 2D lineshapes. The fast component reflects motions of DNA helix, counterions, and water shell. Its higher amplitude at high hydration level reveals a significant contribution of water to the fluctuating forces. The slow component reflects disorder-induced inhomogeneous broadening.
Collapse
Affiliation(s)
- Biswajit Guchhait
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , 12489 Berlin, Germany
| | - Yingliang Liu
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , 12489 Berlin, Germany
| | - Torsten Siebert
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , 12489 Berlin, Germany
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , 12489 Berlin, Germany
| |
Collapse
|
13
|
Yamaguchi N, Zouzumi YK, Shimada N, Nakano SI, Sugimoto N, Maruyama A, Miyoshi D. A reversible B–A transition of DNA duplexes induced by synthetic cationic copolymers. Chem Commun (Camb) 2016; 52:7446-9. [DOI: 10.1039/c6cc02237e] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Reversible B–A transitions of DNA duplexes were induced by synthetic cationic and anionic polymers.
Collapse
Affiliation(s)
- Nonoka Yamaguchi
- Faculty of Frontiers of Innovative Research in Science and Technology (FIRST)
- Konan University
- Kobe 650-0047
- Japan
| | - Yu-ki Zouzumi
- Faculty of Frontiers of Innovative Research in Science and Technology (FIRST)
- Konan University
- Kobe 650-0047
- Japan
| | - Naohiko Shimada
- Department of Biomolecular Engineering
- Graduate School of Bioscience and Biotechnology
- Tokyo Institute of Technology
- Yokohama 226-8501
- Japan
| | - Shu-ichi Nakano
- Faculty of Frontiers of Innovative Research in Science and Technology (FIRST)
- Konan University
- Kobe 650-0047
- Japan
| | - Naoki Sugimoto
- Faculty of Frontiers of Innovative Research in Science and Technology (FIRST)
- Konan University
- Kobe 650-0047
- Japan
- Frontier Institute for Biomolecular Engineering Research (FIBER)
| | - Atsushi Maruyama
- Department of Biomolecular Engineering
- Graduate School of Bioscience and Biotechnology
- Tokyo Institute of Technology
- Yokohama 226-8501
- Japan
| | - Daisuke Miyoshi
- Faculty of Frontiers of Innovative Research in Science and Technology (FIRST)
- Konan University
- Kobe 650-0047
- Japan
| |
Collapse
|
14
|
Siebert T, Guchhait B, Liu Y, Costard R, Elsaesser T. Anharmonic Backbone Vibrations in Ultrafast Processes at the DNA–Water Interface. J Phys Chem B 2015; 119:9670-7. [DOI: 10.1021/acs.jpcb.5b04499] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Torsten Siebert
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Max-Born-Strasse 2a, D-12489 Berlin, Germany
| | - Biswajit Guchhait
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Max-Born-Strasse 2a, D-12489 Berlin, Germany
| | - Yingliang Liu
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Max-Born-Strasse 2a, D-12489 Berlin, Germany
| | - Rene Costard
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Max-Born-Strasse 2a, D-12489 Berlin, Germany
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Max-Born-Strasse 2a, D-12489 Berlin, Germany
| |
Collapse
|
15
|
Son I, Shek YL, Dubins DN, Chalikian TV. Hydration Changes Accompanying Helix-to-Coil DNA Transitions. J Am Chem Soc 2014; 136:4040-7. [DOI: 10.1021/ja5004137] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Ikbae Son
- Department of Pharmaceutical
Sciences, Leslie Dan Faculty of Pharmacy, University of Toronto, 144 College Street, Toronto, Ontario M5S 3M2, Canada
| | - Yuen Lai Shek
- Department of Pharmaceutical
Sciences, Leslie Dan Faculty of Pharmacy, University of Toronto, 144 College Street, Toronto, Ontario M5S 3M2, Canada
| | - David N. Dubins
- Department of Pharmaceutical
Sciences, Leslie Dan Faculty of Pharmacy, University of Toronto, 144 College Street, Toronto, Ontario M5S 3M2, Canada
| | - Tigran V. Chalikian
- Department of Pharmaceutical
Sciences, Leslie Dan Faculty of Pharmacy, University of Toronto, 144 College Street, Toronto, Ontario M5S 3M2, Canada
| |
Collapse
|
16
|
Carbon-14 decay as a source of non-canonical bases in DNA. Biochim Biophys Acta Gen Subj 2014; 1840:526-34. [DOI: 10.1016/j.bbagen.2013.10.003] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2013] [Revised: 08/27/2013] [Accepted: 10/01/2013] [Indexed: 11/22/2022]
|
17
|
Nolin F, Michel J, Wortham L, Tchelidze P, Balossier G, Banchet V, Bobichon H, Lalun N, Terryn C, Ploton D. Changes to cellular water and element content induced by nucleolar stress: investigation by a cryo-correlative nano-imaging approach. Cell Mol Life Sci 2013; 70:2383-94. [PMID: 23385351 PMCID: PMC11113571 DOI: 10.1007/s00018-013-1267-7] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2012] [Revised: 01/10/2013] [Accepted: 01/14/2013] [Indexed: 01/01/2023]
Abstract
The cell is a crowded volume, with estimated mean mass percentage of macromolecules and of water ranging from 7.5 to 45 and 55 to 92.5 %, respectively. However, the concentrations of macromolecules and water at the nanoscale within the various cell compartments are unknown. We recently developed a new approach, correlative cryo-analytical scanning transmission electron microscopy, for mapping the quantity of water within compartments previously shown to display GFP-tagged protein fluorescence on the same ultrathin cryosection. Using energy-dispersive X-ray spectrometry (EDXS), we then identified various elements (C, N, O, P, S, K, Cl, Mg) in these compartments and quantified them in mmol/l. Here, we used this new approach to quantify water and elements in the cytosol, mitochondria, condensed chromatin, nucleoplasm, and nucleolar components of control and stressed cancerous cells. The water content of the control cells was between 60 and 83 % (in the mitochondria and nucleolar fibrillar centers, respectively). Potassium was present at concentrations of 128-462 mmol/l in nucleolar fibrillar centers and condensed chromatin, respectively. The induction of nucleolar stress by treatment with a low dose of actinomycin-D to inhibit rRNA synthesis resulted in both an increase in water content and a decrease in the elements content in all cell compartments. We generated a nanoscale map of water and elements within the cell compartments, providing insight into their changes induced by nucleolar stress.
Collapse
Affiliation(s)
- Frédérique Nolin
- Laboratoire de Recherche en Nanosciences, Université de Reims Champagne Ardenne, Reims, France
| | - Jean Michel
- Laboratoire de Recherche en Nanosciences, Université de Reims Champagne Ardenne, Reims, France
| | - Laurence Wortham
- Laboratoire de Recherche en Nanosciences, Université de Reims Champagne Ardenne, Reims, France
| | - Pavel Tchelidze
- CNRS FRE 3481, Université de Reims Champagne Ardenne, Reims, France
| | - Gérard Balossier
- Laboratoire de Recherche en Nanosciences, Université de Reims Champagne Ardenne, Reims, France
| | - Vincent Banchet
- Laboratoire de Recherche en Nanosciences, Université de Reims Champagne Ardenne, Reims, France
| | - Hélène Bobichon
- CNRS FRE 3481, Université de Reims Champagne Ardenne, Reims, France
| | - Nathalie Lalun
- CNRS FRE 3481, Université de Reims Champagne Ardenne, Reims, France
| | - Christine Terryn
- Plate-forme IBISA, SFR CAP-SANTE, Université de Reims Champagne Ardenne, Reims, France
| | - Dominique Ploton
- CNRS FRE 3481, Université de Reims Champagne Ardenne, Reims, France
| |
Collapse
|
18
|
Mills M, Orr BG, Banaszak Holl MM, Andricioaei I. Attractive hydration forces in DNA-dendrimer interactions on the nanometer scale. J Phys Chem B 2013; 117:973-81. [PMID: 23234339 PMCID: PMC3633417 DOI: 10.1021/jp309616t] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
The energetic contribution of attractive hydration forces arising from water ordering is an interesting but often neglected aspect of macromolecular interactions. Ordering effects of water can bring about cooperativity in many intermolecular transactions, in both the short and long range. Given its high charge density, this is of particular importance for DNA. For instance, in nanotechnology, highly charged dendrimers are used for DNA compaction and transfection. Hypothesizing that water ordering and hydration forces should be maximal for DNA complexes that show charge complementarity (positive-negative), we present here an analysis of water ordering from molecular dynamics simulations and free energy calculations of the interaction between DNA and a nanoparticle with a high positive charge density. Our results indicate not only that complexation of the dendrimer with DNA affects the local water structure but also that ordered water molecules facilitate long-range interactions between the molecules. This contributes significantly to the free energy of binding of dendrimers to DNA and extends the interaction well beyond the electrostatic range of the DNA. Such water effects are of potentially substantial importance in cases when molecules appear to recognize each other across sizable distances, or for which kinetic rates are too fast to be due to pure diffusion. Our results are in good agreement with experiments on the role of solvent in DNA condensation by multivalent cations and exemplify a microscopic realization of mean-field phenomenological theories for hydration forces between mesoscopic surfaces.
Collapse
|
19
|
Fujimoto T, Nakano SI, Sugimoto N, Miyoshi D. Thermodynamics-hydration relationships within loops that affect G-quadruplexes under molecular crowding conditions. J Phys Chem B 2012; 117:963-72. [PMID: 23153339 DOI: 10.1021/jp308402v] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We systematically investigated the effects of loop length on the conformation, thermodynamic stability, and hydration of DNA G-quadruplexes under dilute and molecular crowding conditions in the presence of Na(+). Structural analysis showed that molecular crowding induced conformational switches of oligonucleotides with the longer guanine stretch and the shorter thymine loop. Thermodynamic parameters further demonstrated that the thermodynamic stability of G-quadruplexes increased by increasing the loop length from two to four, whereas it decreased by increasing the loop length from four to six. Interestingly, we found by osmotic pressure analysis that the number of water molecules released from the G-quadruplex decreased with increasing thermodynamic stability. We assumed that base-stacking interactions within the loops not only stabilized the whole G-quadruplex structure but also created hydration sites by accumulating nucleotide functional groups. The molecular crowding effects on the stability of G-quadruplexes composed of abasic sites, which reduce the stacking interactions at the loops, further demonstrated that G-quadruplexes with fewer stacking interactions within the loops released a larger number of water molecules upon folding. These results showed that the stacking interactions within the loops determined the thermodynamic stability and hydration of the whole G-quadruplex.
Collapse
Affiliation(s)
- Takeshi Fujimoto
- Faculty of Frontiers of Innovative Research in Science and Technology, Konan University, 7-1-20, Minatojima-minamimachi, Chuo-ku, Kobe 650-0047, Japan
| | | | | | | |
Collapse
|
20
|
Dixit SB, Mezei M, Beveridge DL. Studies of base pair sequence effects on DNA solvation based on all-atom molecular dynamics simulations. J Biosci 2012; 37:399-421. [PMID: 22750979 DOI: 10.1007/s12038-012-9223-5] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Detailed analyses of the sequence-dependent solvation and ion atmosphere of DNA are presented based on molecular dynamics (MD) simulations on all the 136 unique tetranucleotide steps obtained by the ABC consortium using the AMBER suite of programs. Significant sequence effects on solvation and ion localization were observed in these simulations. The results were compared to essentially all known experimental data on the subject. Proximity analysis was employed to highlight the sequence dependent differences in solvation and ion localization properties in the grooves of DNA. Comparison of the MD-calculated DNA structure with canonical A- and B-forms supports the idea that the G/C-rich sequences are closer to canonical A- than B-form structures, while the reverse is true for the poly A sequences, with the exception of the alternating ATAT sequence. Analysis of hydration density maps reveals that the flexibility of solute molecule has a significant effect on the nature of observed hydration. Energetic analysis of solute-solvent interactions based on proximity analysis of solvent reveals that the GC or CG base pairs interact more strongly with water molecules in the minor groove of DNA that the AT or TA base pairs, while the interactions of the AT or TA pairs in the major groove are stronger than those of the GC or CG pairs. Computation of solvent-accessible surface area of the nucleotide units in the simulated trajectories reveals that the similarity with results derived from analysis of a database of crystallographic structures is excellent. The MD trajectories tend to follow Manning's counterion condensation theory, presenting a region of condensed counterions within a radius of about 17 A from the DNA surface independent of sequence. The GC and CG pairs tend to associate with cations in the major groove of the DNA structure to a greater extent than the AT and TA pairs. Cation association is more frequent in the minor groove of AT than the GC pairs. In general, the observed water and ion atmosphere around the DNA sequences is the MD simulation is in good agreement with experimental observations.
Collapse
Affiliation(s)
- Surjit B Dixit
- Chemistry Department and Molecular Biophysics Program, Wesleyan University, Middletown, CT 06457, USA
| | | | | |
Collapse
|
21
|
Yang M, Szyc Ł, Elsaesser T. Vibrational dynamics of the water shell of DNA studied by femtosecond two-dimensional infrared spectroscopy. J Photochem Photobiol A Chem 2012. [DOI: 10.1016/j.jphotochem.2011.11.014] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
22
|
Yang M, Szyc Ł, Elsaesser T. Decelerated water dynamics and vibrational couplings of hydrated DNA mapped by two-dimensional infrared spectroscopy. J Phys Chem B 2011; 115:13093-100. [PMID: 21972952 DOI: 10.1021/jp208166w] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Double-stranded DNA oligomers containing 23 alternating adenine-thymine base pairs are studied at different hydration levels by femtosecond two-dimensional (2D) infrared spectrosopy. Coupled NH stretching modes of the A-T pairs and OH stretching excitations of the water shell are discerned in the 2D spectra. Limited changes of NH stretching frequencies and line shapes with increasing hydration suggest spectral dynamics governed by DNA rather than water fluctuations. In contrast, OH stretching excitations of the water shell around fully hydrated DNA undergo spectral diffusion on a ~500 fs time scale. The center line slopes of the 2D spectra of hydrated DNA demonstrate a slower decay of the frequency-time correlation function (TCF) than that in neat water, as is evident from a comparison with 2D spectra of neat H(2)O and theoretical TCFs. We attribute this behavior to reduced structural fluctuations of the water shell and a reduced rate of resonant OH stretching energy transfer.
Collapse
Affiliation(s)
- Ming Yang
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin, Germany
| | | | | |
Collapse
|
23
|
Lin B, Wong KY, Hu C, Kokubo H, Pettitt BM. Fast Calculations of Electrostatic Solvation Free Energy from Reconstructed Solvent Density using proximal Radial Distribution Functions. J Phys Chem Lett 2011; 2:1626-1632. [PMID: 21765968 PMCID: PMC3134239 DOI: 10.1021/jz200609v] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Although detailed atomic models may be applied for a full description of solvation, simpler phenomenological models are particularly useful to interpret the results for scanning many, large, complex systems where a full atomic model is too computationally expensive to use. Among the most costly are solvation free energy evaluations by simulation. Here we develop a fast way to calculate electrostatic solvation free energy while retaining much of the accuracy of explicit solvent free energy simulation. The basis of our method is to treat the solvent not as a structureless dielectric continuum, but as a structured medium by making use of universal proximal radial distribution functions. Using a deca-alanine peptide as a test case, we compare the use of our theory with free energy simulations and traditional continuum estimates of the electrostatic solvation free energy.
Collapse
Affiliation(s)
- Bin Lin
- Program in Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, TX 77030 USA
| | | | | | | | | |
Collapse
|
24
|
Pino S, Biasiucci M, Scardamaglia M, Gigli G, Betti MG, Mariani C, Di Mauro E. Nonenzymatic ligation of an RNA oligonucleotide analyzed by atomic force microscopy. J Phys Chem B 2011; 115:6296-303. [PMID: 21500793 DOI: 10.1021/jp200832w] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The products of ligation reaction of a 24 nucleotides long PolyA RNA adsorbed on mica were observed by atomic force microscopy. The occurrence of oligonucleotides at different degrees of polymerization has been quantitatively studied before and after ligation reaction. The microscopy images at the nanoscale show that nonenzymatic ligation of pristine RNA monomers results in the formation of supramolecular aggregates, with prevalence of dimers and tetramers. Analytical conditions were defined allowing the identification, the quantitative evaluation, and their distribution after ligation reaction, also providing an estimate of the degree of hydration of the objects. Such investigation is of particular biological relevance and provides the simplest yet model system for direct investigation of RNA reactions by advanced microscopy.
Collapse
Affiliation(s)
- Samanta Pino
- Department of Physics, Sapienza University of Rome, P.le A. Moro 2, 00185 Rome, Italy
| | | | | | | | | | | | | |
Collapse
|
25
|
Howard JJ, Lynch GC, Pettitt BM. Ion and solvent density distributions around canonical B-DNA from integral equations. J Phys Chem B 2010; 115:547-56. [PMID: 21190358 DOI: 10.1021/jp107383s] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
We calculate the water and ion spatial distributions around charged oligonucleotides using a renormalized three-dimensional reference interaction site theory coupled with the HNC closure. Our goal is to understand the balance between inter-DNA strand forces and solvation forces as a function of oligonucleotide length in the short strand limit. The DNA is considered in aqueous electrolyte solutions of 1 M KCl, 0.1 M KCl, or 0.1 M NaCl. The current theoretical results are compared to molecular dynamics (MD) simulations and experiments. It is found that the integral equation (IE) theory replicates the MD and the experimental results for the base-specific hydration patterns in both the major and the minor grooves. We are also able to discern characteristic structural pattern differences between Na(+) and K(+) ions. When compared to Poisson-Boltzmann methods, the IE theory, like simulation, predicts a richly structured ion environment, which is better described as multilayer rather than double layer.
Collapse
Affiliation(s)
- Jesse J Howard
- Department of Chemistry, University of Houston, Houston, Texas 77204-5003, USA
| | | | | |
Collapse
|
26
|
Virtanen JJ, Makowski L, Sosnick TR, Freed KF. Modeling the hydration layer around proteins: HyPred. Biophys J 2010; 99:1611-9. [PMID: 20816074 DOI: 10.1016/j.bpj.2010.06.027] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2010] [Revised: 06/14/2010] [Accepted: 06/15/2010] [Indexed: 11/25/2022] Open
Abstract
Protein hydration plays an integral role in determining protein function and stability. We develop a simple method with atomic level precision for predicting the solvent density near the surface of a protein. A set of proximal radial distribution functions are defined and calculated for a series of different atom types in proteins using all-atom, explicit solvent molecular dynamic simulations for three globular proteins. A major improvement in predicting the hydration layer is found when the protein is held immobile during the simulations. The distribution functions are used to develop a model for predicting the hydration layer with sub-1-Angstrom resolution without the need for additional simulations. The model and the distribution functions for a given protein are tested in their ability to reproduce the hydration layer from the simulations for that protein, as well as those for other proteins and for simulations in which the protein atoms are mobile. Predictions for the density of water in the hydration shells are then compared with high occupancy sites observed in crystal structures. The accuracy of both tests demonstrates that the solvation model provides a basis for quantitatively understanding protein solvation and thereby predicting the hydration layer without additional simulations.
Collapse
Affiliation(s)
- Jouko J Virtanen
- Department of Chemistry, The University of Chicago, Chicago, Illinois, USA
| | | | | | | |
Collapse
|
27
|
Darré L, Machado MR, Dans PD, Herrera FE, Pantano S. Another Coarse Grain Model for Aqueous Solvation: WAT FOUR? J Chem Theory Comput 2010. [DOI: 10.1021/ct100379f] [Citation(s) in RCA: 94] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Leonardo Darré
- Institut Pasteur de Montevideo, Calle Mataojo 2020, CP 11400, Montevideo, Uruguay, and Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Avda. Rivadavia 1917 - CP C1033AAJ - Cdad. de Buenos Aires, Argentina
| | - Matías R. Machado
- Institut Pasteur de Montevideo, Calle Mataojo 2020, CP 11400, Montevideo, Uruguay, and Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Avda. Rivadavia 1917 - CP C1033AAJ - Cdad. de Buenos Aires, Argentina
| | - Pablo D. Dans
- Institut Pasteur de Montevideo, Calle Mataojo 2020, CP 11400, Montevideo, Uruguay, and Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Avda. Rivadavia 1917 - CP C1033AAJ - Cdad. de Buenos Aires, Argentina
| | - Fernando E. Herrera
- Institut Pasteur de Montevideo, Calle Mataojo 2020, CP 11400, Montevideo, Uruguay, and Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Avda. Rivadavia 1917 - CP C1033AAJ - Cdad. de Buenos Aires, Argentina
| | - Sergio Pantano
- Institut Pasteur de Montevideo, Calle Mataojo 2020, CP 11400, Montevideo, Uruguay, and Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Avda. Rivadavia 1917 - CP C1033AAJ - Cdad. de Buenos Aires, Argentina
| |
Collapse
|
28
|
Kang NS, No KT, Jhon MS. A Study of the Hydrogen-bonded Network Around the Left-handed and Right-handed DNA in an Ionic Solution. MOLECULAR SIMULATION 2010. [DOI: 10.1080/0892702031000065764] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Nam Sook Kang
- a Department of Chemistry and CAMD Research Center , Soong Sil University , Seoul , 156-743 , South Korea
| | - Kyung Tai No
- a Department of Chemistry and CAMD Research Center , Soong Sil University , Seoul , 156-743 , South Korea
| | - Mu Shik Jhon
- b Department of Chemistry , Korea Advanced Institute of Science and Technology , 373-1 Kusung-dong Yusung-gu, Taejon , 305-701 , South Korea
| |
Collapse
|
29
|
Naômé A, Schyman P, Laaksonen A, Vercauteren DP. Molecular Dynamics Simulation of 8-Oxoguanine Containing DNA Fragments Reveals Altered Hydration and Ion Binding Patterns. J Phys Chem B 2010; 114:4789-801. [DOI: 10.1021/jp1000539] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Aymeric Naômé
- Laboratoire de Physico-Chimie Informatique, Unité de Chimie Physique Théorique et Structurale, University of Namur (FUNDP), Namur, Belgium, and Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, Sweden
| | - Patric Schyman
- Laboratoire de Physico-Chimie Informatique, Unité de Chimie Physique Théorique et Structurale, University of Namur (FUNDP), Namur, Belgium, and Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, Sweden
| | - Aatto Laaksonen
- Laboratoire de Physico-Chimie Informatique, Unité de Chimie Physique Théorique et Structurale, University of Namur (FUNDP), Namur, Belgium, and Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, Sweden
| | - Daniel P. Vercauteren
- Laboratoire de Physico-Chimie Informatique, Unité de Chimie Physique Théorique et Structurale, University of Namur (FUNDP), Namur, Belgium, and Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, Sweden
| |
Collapse
|
30
|
Yonetani Y, Kono H. Sequence dependencies of DNA deformability and hydration in the minor groove. Biophys J 2009; 97:1138-47. [PMID: 19686662 PMCID: PMC2726331 DOI: 10.1016/j.bpj.2009.05.049] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2009] [Revised: 05/13/2009] [Accepted: 05/29/2009] [Indexed: 10/20/2022] Open
Abstract
DNA deformability and hydration are both sequence-dependent and are essential in specific DNA sequence recognition by proteins. However, the relationship between the two is not well understood. Here, systematic molecular dynamics simulations of 136 DNA sequences that differ from each other in their central tetramer revealed that sequence dependence of hydration is clearly correlated with that of deformability. We show that this correlation can be illustrated by four typical cases. Most rigid basepair steps are highly likely to form an ordered hydration pattern composed of one water molecule forming a bridge between the bases of distinct strands, but a few exceptions favor another ordered hydration composed of two water molecules forming such a bridge. Steps with medium deformability can display both of these hydration patterns with frequent transition. Highly flexible steps do not have any stable hydration pattern. A detailed picture of this correlation demonstrates that motions of hydration water molecules and DNA bases are tightly coupled with each other at the atomic level. These results contribute to our understanding of the entropic contribution from water molecules in protein or drug binding and could be applied for the purpose of predicting binding sites.
Collapse
Affiliation(s)
- Yoshiteru Yonetani
- Computational Biology Group, Quantum Beam Science Directorate, Japan Atomic Energy Agency, Kyoto, Japan
| | - Hidetoshi Kono
- Computational Biology Group, Quantum Beam Science Directorate, Japan Atomic Energy Agency, Kyoto, Japan
- Quantum Bioinformatics Team, Center for Computational Science and e-Systems, Japan Atomic Energy Agency, Kyoto, Japan
- PRESTO, Japan Science and Technology Agency, Saitama, Japan
| |
Collapse
|
31
|
Dong F, Wagoner JA, Baker NA. Assessing the performance of implicit solvation models at a nucleic acid surface. Phys Chem Chem Phys 2008; 10:4889-902. [PMID: 18688533 PMCID: PMC2538626 DOI: 10.1039/b807384h] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Implicit solvation models are popular alternatives to explicit solvent methods due to their ability to "pre-average" solvent behavior and thus reduce the need for computationally-expensive sampling. Previously, we have demonstrated that Poisson-Boltzmann models for polar solvation and integral-based models for nonpolar solvation can reproduce explicit solvation forces in a low-charge density protein system. In the present work, we examine the ability of these continuum models to describe solvation forces at the surface of a RNA hairpin. While these models do not completely describe all of the details of solvent behavior at this highly-charged biomolecular interface, they do provide a reasonable description of average solvation forces and therefore show significant promise for developing more robust implicit descriptions of solvent around nucleic acid systems for use in biomolecular simulation and modeling. Additionally, we observe fairly good transferability in the nonpolar model parameters optimized for protein systems, suggesting its robustness for modeling general nonpolar solvation phenomena in biomolecular systems.
Collapse
Affiliation(s)
- Feng Dong
- Merck & Co., Inc., 770 Sumneytown Pike, P.O. Box 4, WP42-330, West Point, PA 19486, USA. E-mail:
| | - Jason A. Wagoner
- Department of Chemistry, Stanford University, 333 Campus Drive #121, Mailbox 13, Stanford, CA 94305-5080, USA. E-mail:
| | - Nathan A. Baker
- To whom correspondence should be addressed. Department of Biochemistry and Molecular Biophysics, Center for Computational Biology, Washington University in St. Louis, 700. S. Euclid Ave., St. Louis, MO 63110, USA. E-mail:
| |
Collapse
|
32
|
Yonetani Y, Maruyama Y, Hirata F, Kono H. Comparison of DNA hydration patterns obtained using two distinct computational methods, molecular dynamics simulation and three-dimensional reference interaction site model theory. J Chem Phys 2008; 128:185102. [PMID: 18532849 DOI: 10.1063/1.2904865] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Because proteins and DNA interact with each other and with various small molecules in the presence of water molecules, we cannot ignore their hydration when discussing their structural and energetic properties. Although high-resolution crystal structure analyses have given us a view of tightly bound water molecules on their surface, the structural data are still insufficient to capture the detailed configurations of water molecules around the surface of these biomolecules. Thanks to the invention of various computational algorithms, computer simulations can now provide an atomic view of hydration. Here, we describe the apparent patterns of DNA hydration calculated by using two different computational methods: Molecular dynamics (MD) simulation and three-dimensional reference interaction site model (3D-RISM) theory. Both methods are promising for obtaining hydration properties, but until now there have been no thorough comparisons of the calculated three-dimensional distributions of hydrating water. This rigorous comparison showed that MD and 3D-RISM provide essentially similar hydration patterns when there is sufficient sampling time for MD and a sufficient number of conformations to describe molecular flexibility for 3D-RISM. This suggests that these two computational methods can be used to complement one another when evaluating the reliability of the calculated hydration patterns.
Collapse
Affiliation(s)
- Yoshiteru Yonetani
- Computational Biology Group, Quantum Beam Science Directorate, Japan Atomic Energy Agency, 8-1 Umemidai, Kizugawa, Kyoto, Japan
| | | | | | | |
Collapse
|
33
|
Gantchev TG, Hunting DJ. Probing the interactions of the solvated electron with DNA by molecular dynamics simulations: bromodeoxyuridine substituted DNA. J Mol Model 2008; 14:451-64. [PMID: 18414908 DOI: 10.1007/s00894-008-0296-x] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2007] [Accepted: 02/28/2008] [Indexed: 11/26/2022]
Abstract
Solvated electrons ((e-)(aq)) are produced during water radiolysis and can interact with biological substrates, including DNA. To augment DNA damage, radiosensitizers such as bromo-deoxyuridine (BUdR), often referred to as an "electron affinic radiosensitizer", are incorporated in place of isosteric thymidine. However, little is known about the primary interactions of (e-)(aq) with DNA. In the present study we addressed this problem by applying molecular modeling and molecular dynamics (MD) simulations to a system of normal (BUdR.A)-DNA and a hydrated electron, where the excess electron was modeled as a localized (e-)(H2O)6 anionic cluster. Our goals were to evaluate the suitability of the MD simulations for this application; to characterize the motion of (e-)(aq) around DNA (e.g., diffusion coefficients); to identify and describe configurational states of close (e-)(aq) localization to DNA; and to evaluate the structural dynamics of DNA in the presence of (e-)(aq). The results indicate that (e-)(aq) has distinct space-preferences for forming close contacts with DNA and is more likely to interact directly with nucleotides other than BUdR. Several classes of DNA - (e-)(aq) contact sites, all within the major groove, were distinguished depending on the structure of the intermediate water layer H-bonding pattern (or its absence, i.e., a direct H-bonding of (e-)(aq) with DNA bases). Large-scale structural perturbations were identified during and after the (e-)(aq) approached the DNA from the major groove side, coupled with deeper penetration of sodium counterions in the minor groove.
Collapse
Affiliation(s)
- Tsvetan G Gantchev
- Department of Nuclear Medicine & Radiobiology, Faculty of Medicine, Université de Sherbrooke, Sherbrooke, QC J1H 5N4, Canada.
| | | |
Collapse
|
34
|
Kumar A, Sevilla MD, Suhai S. Microhydration of the guanine-cytosine (GC) base pair in the neutral and anionic radical states: a density functional study. J Phys Chem B 2008; 112:5189-98. [PMID: 18380501 DOI: 10.1021/jp710957p] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A density functional study of the effects of microhydration on the guanine-cytosine (GC) base pair and its anion radical is presented. Geometries of the GC base pair in the presence of 6 and 11 water molecules were fully optimized in the neutral (GC-nH2O) and anion radical [(GC-nH2O)*-] (n = 6 and 11) states using the B3LYP method and the 6-31+G** basis set. Further, vibrational frequency analysis at the same level of theory (B3LYP/6-31+G**) was also performed to ensure the existence of local minima in these hydrated structures. It was found that water molecules surrounding the GC base pair have significant effects on the geometry of the GC base pair and promote nonplanarity in the GC base pair. The calculated structures were found to be in good agreement with those observed experimentally and obtained in molecular dynamics (MD) simulation studies. The water molecules in neutral GC-nH2O complexes lie near the ring plane of the GC base pair where they undergo hydrogen bonding with both GC and each other. However, in the GC anion radical complexes (GC-nH2O, n = 6, 11), the water molecules are displaced substantially from the GC ring plane. For GC-11H2O*-, a water molecule is hydrogen-bonded with the C6 atom of the cytosine base. We found that the hydration shell initially destabilizes the GC base pair toward electron capture as a transient anion. Energetically unstable diffuse states in the hydration shell are suggested to provide an intermediate state for the excess electron before molecular reorganization of the water molecules and the base pair results in a stable anion formation. The singly occupied molecular orbital (SOMO) in the anion radical complexes clearly shows that an excess electron localizes into a pi orbital of cytosine. The zero-point-energy (ZPE-) corrected adiabatic electron affinities (AEAs) of the GC-6H2O and GC-11H2O complexes, at the B3LYP/6-31+G** level of theory, were found to be 0.74 and 0.95 eV, respectively. However, the incorporation of bulk water as a solvent using the polarized continuum model (PCM) increases the EAs of these complexes to 1.77 eV.
Collapse
Affiliation(s)
- Anil Kumar
- Department of Chemistry, Oakland University, Rochester, Michigan 48309, USA
| | | | | |
Collapse
|
35
|
Tanizaki S, Clifford J, Connelly BD, Feig M. Conformational sampling of peptides in cellular environments. Biophys J 2007; 94:747-59. [PMID: 17905846 PMCID: PMC2186233 DOI: 10.1529/biophysj.107.116236] [Citation(s) in RCA: 48] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Biological systems provide a complex environment that can be understood in terms of its dielectric properties. High concentrations of macromolecules and cosolvents effectively reduce the dielectric constant of cellular environments, thereby affecting the conformational sampling of biomolecules. To examine this effect in more detail, the conformational preference of alanine dipeptide, poly-alanine, and melittin in different dielectric environments is studied with computer simulations based on recently developed generalized Born methodology. Results from these simulations suggest that extended conformations are favored over alpha-helical conformations at the dipeptide level at and below dielectric constants of 5-10. Furthermore, lower-dielectric environments begin to significantly stabilize helical structures in poly-alanine at epsilon = 20. In the more complex peptide melittin, different dielectric environments shift the equilibrium between two main conformations: a nearly fully extended helix that is most stable in low dielectrics and a compact, V-shaped conformation consisting of two helices that is preferred in higher dielectric environments. An additional conformation is only found to be significantly populated at intermediate dielectric constants. Good agreement with previous studies of different peptides in specific, less-polar solvent environments, suggest that helix stabilization and shifts in conformational preferences in such environments are primarily due to a reduced dielectric environment rather than specific molecular details. The findings presented here make predictions of how peptide sampling may be altered in dense cellular environments with reduced dielectric response.
Collapse
Affiliation(s)
- Seiichiro Tanizaki
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan
| | - Jacob Clifford
- Department of Physics, Michigan State University, East Lansing, Michigan
| | - Brian D. Connelly
- Department of Computer Science and Engineering, Michigan State University, East Lansing, Michigan
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan
- Department of Computer Science and Engineering, Michigan State University, East Lansing, Michigan
- Department of Chemistry, Michigan State University, East Lansing, Michigan
- Address reprint requests to Michael Feig, Tel.: 517-432-7439; Fax: 517-353-9334.
| |
Collapse
|
36
|
Chocholousová J, Feig M. Implicit solvent simulations of DNA and DNA-protein complexes: agreement with explicit solvent vs experiment. J Phys Chem B 2007; 110:17240-51. [PMID: 16928023 DOI: 10.1021/jp0627675] [Citation(s) in RCA: 60] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Molecular dynamics simulations of biomolecules with implicit solvent reduce the computational cost and complexity of such simulations so that longer time scales and larger system sizes can be reached. While implicit solvent simulations of proteins have become well established, the success of implicit solvent in the simulation of nucleic acids has not been fully established to date. Results obtained in this study demonstrate that stable and efficient simulations of DNA and a protein-DNA complex can be achieved with an implicit solvent model based on continuum dielectric electrostatics. Differences in conformational sampling of DNA with two sets of atomic radii that are used to define the dielectric interface between the solute and the continuum dielectric model of the solvent are investigated. Results suggest that depending on the choice of atomic radii agreement is either closer to experimental data or to explicit solvent simulations. Furthermore, partial conformational transitions toward A-DNA conformations when salt is added within the implicit solvent framework are observed.
Collapse
Affiliation(s)
- Jana Chocholousová
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824-1319, USA
| | | |
Collapse
|
37
|
Auffinger P, Hashem Y. Nucleic acid solvation: from outside to insight. Curr Opin Struct Biol 2007; 17:325-33. [PMID: 17574833 DOI: 10.1016/j.sbi.2007.05.008] [Citation(s) in RCA: 112] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2007] [Revised: 03/28/2007] [Accepted: 05/31/2007] [Indexed: 11/18/2022]
Abstract
Nucleic acids are polyanionic molecules that were historically considered to be solely surrounded by a shell of water molecules and a neutralizing cloud of monovalent and divalent cations. In this respect, recent experimental and theoretical reports demonstrate that water molecules within complex nucleic acid structures can display very long residency times, and assist drug binding and catalytic reactions. Finally, anions can also bind to these polyanionic systems. Many of these recent insights are provided by state-of-the-art molecular dynamics simulations of nucleic acid systems, which will be described together with relevant methodological issues.
Collapse
Affiliation(s)
- Pascal Auffinger
- Architecture et réactivité de l'ARN, Université Louis Pasteur de Strasbourg, CNRS, IBMC, 15 rue René Descartes, 67084 Strasbourg, France.
| | | |
Collapse
|
38
|
|
39
|
Manalo MN, Kong X, LiWang A. Sensitivity of hydrogen bonds of DNA and RNA to hydration, as gauged by 1JNH measurements in ethanol-water mixtures. JOURNAL OF BIOMOLECULAR NMR 2007; 37:257-63. [PMID: 17310327 DOI: 10.1007/s10858-006-9132-8] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/22/2006] [Revised: 12/08/2006] [Accepted: 12/08/2006] [Indexed: 05/14/2023]
Abstract
Hydrogen-bond lengths of nucleic acids are (1) longer in DNA than in RNA, and (2) sequence dependent. The physicochemical basis for these variations in hydrogen-bond lengths is unknown, however. Here, the notion that hydration plays a significant role in nucleic acid hydrogen-bond lengths is tested. Watson-Crick N1...N3 hydrogen-bond lengths of several DNA and RNA duplexes are gauged using imino 1J(NH) measurements, and ethanol is used as a cosolvent to lower water activity. We find that 1J(NH) values of DNA and RNA become less negative with added ethanol, which suggests that mild dehydration reduces hydrogen-bond lengths even as the overall thermal stabilities of these duplexes decrease. The 1J(NH) of DNA are increased in 8 mol% ethanol to those of RNA in water, which suggests that the greater hydration of DNA plays a significant role in its longer hydrogen bonds. The data also suggest that ethanol-induced dehydration is greater for the more hydrated G:C base pairs and thereby results in greater hydrogen-bond shortening than for the less hydrated A:T/U base pairs of DNA and RNA.
Collapse
Affiliation(s)
- Marlon N Manalo
- Department of Biochemistry & Biophysics, Texas A&M University, 2128 TAMU, College Station, TX 77843-2128, USA
| | | | | |
Collapse
|
40
|
René B, Masliah G, Antri SE, Fermandjian S, Mauffret O. Conformations and Dynamics of the Phosphodiester Backbone of a DNA Fragment That Bears a Strong Topoisomerase II Cleavage Site. J Phys Chem B 2007; 111:4235-43. [PMID: 17391020 DOI: 10.1021/jp0683115] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The dynamics of the DNA phosphodiester backbone conformations have been studied for a strong topoisomerase II cleavage site (site 22) using molecular dynamics simulations in explicit water and in the presence of sodium ions. We investigated the backbone motions and more particularly the BI/BII transitions involving the epsilon and zeta angles. The consensus cleavage site is adjacent to the phosphate which shows the most important phosphodiester backbone flexibility in the sequence. We infer that these latter properties could be responsible for the preferential cleavage at this site possibly through the perturbation of the cleavage/ligation activities of the topoisomerase II. More generally, the steps pur-pur and pyr-pur are those presenting the highest BII contents. Relations are observed between the backbone phosphodiester BI/BII transitions and the flexibility of the deoxyribose sugar and the helical parameters such as roll. The roll is sequence dependent when the related phosphate is in the BI form, whereas this appears not to be true when it is in the BII form. The BI/BII transitions are associated with water migration, and new relations are observed with counterions. Indeed, it is observed that a strong coupling exists between the BII form and the presence of sodium ions near the adjacent sugar deoxyribose. The presence of sodium ions in the O4' surroundings or their binding could assist the BI to BII transition by furnishing energy. The implications of these new findings and, namely, their importance in the context of the sequence-dependent behavior of BI/BII transitions will be investigated in future studies.
Collapse
Affiliation(s)
- Brigitte René
- UMR 8113 CNRS, Laboratoire de Biotechnologies et Pharmacologie Génétique Appliquée (Ecole Normale Supérieure de Cachan), Institut Gustave Roussy, 39 rue Camille Desmoulins, 94805 Villejuif Cedex, France
| | | | | | | | | |
Collapse
|
41
|
Pal S, Maiti PK, Bagchi B, Hynes JT. Multiple time scales in solvation dynamics of DNA in aqueous solution: the role of water, counterions, and cross-correlations. J Phys Chem B 2007; 110:26396-402. [PMID: 17181299 DOI: 10.1021/jp065690t] [Citation(s) in RCA: 79] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Recent time domain experiments have explored solvation dynamics of a probe located inside a DNA duplex, in an effort to gain information, e.g., on the dynamics of water molecules in the DNA major and minor grooves and their environment. Multiple time constants in the range of a few picoseconds to several nanoseconds were obtained. We have carried out 15 ns long atomistic molecular dynamics simulations to study the solvation dynamics of bases of a 38 base-pair long DNA duplex in an aqueous solution containing counterions. We have computed the energy-energy time correlation function (TCF) of the four individual bases (A, T, G, and C) to characterize the solvation dynamics. All the TCFs display highly nonexponential decay with time. When the trajectories are analyzed with 100 fs time resolution, the TCF of each base shows initial ultrafast decay (with tau1 approximately equal 60-80 fs) followed by two intermediate components (tau2 approximately equal 1 ps, tau3 approximately equal 20-30 ps), in near complete agreement with a recent time domain experiment on DNA solvation. Interestingly, the solvation dynamics of each of the four different nucleotide bases exhibit rather similar time scales. To explore the existence of slow relaxation at longer times reported recently in a series of experiments, we also analyzed the solvation TCFs calculated with longer time trajectories and with a larger time resolution of 1 ps. In this case, an additional slow component with a time constant of the order of 250 ps is observed. Through an analysis of partial solvation TCFs, we find that the slow decay originates mainly from the interaction of the nucleotides with the dipolar water molecules and the counterions. An interesting negative cross-correlation between water and counterions is observed, which makes an important contribution to relaxation at intermediate to longer times.
Collapse
Affiliation(s)
- Subrata Pal
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore-560012, India
| | | | | | | |
Collapse
|
42
|
Wong KY, Vainrub A, Powdrill T, Hogan M, Pettitt BM. A Non-Watson–Crick Motif of Base-pairing on Surfaces for Untethered Oligonucleotides. MOLECULAR SIMULATION 2006. [DOI: 10.1080/0892702031000152181] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
43
|
Feig M, Im W, Brooks CL. Implicit solvation based on generalized Born theory in different dielectric environments. J Chem Phys 2006; 120:903-11. [PMID: 15267926 DOI: 10.1063/1.1631258] [Citation(s) in RCA: 121] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this paper we are investigating the effect of the dielectric environment on atomic Born radii used in generalized Born (GB) methods. Motivated by the Kirkwood expression for the reaction field of a single off-center charge in a spherical cavity, we are proposing extended formalisms for the calculation of Born radii as a function of external and internal dielectric constants. We demonstrate that reaction field energies calculated from environmentally dependent Born radii lead to much improved agreement with Poisson-Boltzmann solutions for low dielectric external environments, such as biological membranes or organic solvent, compared to previous methods where the calculation of Born radii does not depend on the environment. We also examine how this new approach can be applied for the calculation of transfer free energies from vacuum to a given external dielectric for a system with an internal dielectric larger than one. This has not been possible with standard GB theory but is relevant when scoring minimized or average structures with implicit solvent.
Collapse
Affiliation(s)
- Michael Feig
- Department of Molecular Biology, TPC6, The Scripps Research Institute, La Jolla, California 92037, USA
| | | | | |
Collapse
|
44
|
Guzman MR, Liquier J, Taillandier E. Hydration and conformational transitions in DNA, RNA, and mixed DNA-RNA triplexes studied by gravimetry and FTIR spectroscopy. J Biomol Struct Dyn 2005; 23:331-9. [PMID: 16218757 DOI: 10.1080/07391102.2005.10507068] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
We have studied by gravimetric measurements and FTIR spectroscopy the hydration of duplexes and triplexes formed by combinations of dA(n), dT(n), rA(n), and rU(n) strands. Results obtained on hydrated films show important differences in their hydration and in the structural transitions which can be induced by varying the water content of the samples. The number of water molecules per nucleotide (w/n) measured at high relative humidity (98% R.H.) is found to be 21 for dA(n).dT(n) and 15 for rA(n).rU(n). Addition of a third rU(n) strand does not change the number of water molecules per nucleotide: w/n=21 for rU(n)*dA(n).dT(n) and w/n=15 for rU(n)*rA(n).rU(n). On the contrary, the addition of a third dT(n) strand changes the water content but in a different way, depending whether the duplex is DNA or RNA. Thus, a loss of four water molecules per nucleotide is measured for dT(n)*dA(n).dT(n) while an increase of two water molecules per nucleotide is observed for dT(n)*rA(n).rU(n). The final hydration is the same for both triplexes (w/n=17). The desorption profiles obtained by gravimetry and FTIR spectroscopy are similar for the rA(n).rU(n) duplex and the rU(n)*rA(n).rU(n) triplex. On the contrary, the desorption profiles of the dA(n).dT(n) duplex and the triplexes formed with it (rU(n)*dA(n).dT(n) and dT(n)*dA(n).dT(n)) are different from each other. This is correlated with conformational transitions induced by varying the hydration content of the different structures, as shown by FTIR spectroscopy. Modifications of the phosphate group hydration and of the sugar conformation (S to N type repuckering) induced by decrease of the water content are observed in the case of triplexes formed on the dA(n).dT(n) duplex.
Collapse
Affiliation(s)
- M R Guzman
- Laboratoire BioMoCeTi, UMR CNRS 7033, Universite Paris 13, 74 rue Marcel Cachin, F93017 Bobigny Cedex, France
| | | | | |
Collapse
|
45
|
Feig M, Chocholoušová J, Tanizaki S. Extending the horizon: towards the efficient modeling of large biomolecular complexes in atomic detail. Theor Chem Acc 2005. [DOI: 10.1007/s00214-005-0062-4] [Citation(s) in RCA: 48] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|
46
|
Wibowo FR, Rauch C, Trieb M, Liedl KR. M.TaqI facilitates the base flipping via an unusual DNA backbone conformation. Biopolymers 2005; 79:128-38. [PMID: 16047360 DOI: 10.1002/bip.20341] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
MD simulations have been carried out to understand the dynamical behavior of the DNA substrate of the Thermus aquaticus DNA methyltransferase (M.TaqI) in the methylation process at N6 of adenine. As starting structures, an x-ray structure of M.TaqI in complex with DNA and cofactor analogue (PDB code: 1G 38) and free decamer d(GTTCGATGTC)(2) were taken. The x-ray structure shows two consecutive BII substates that are not observed in the free decamer. These consecutive BII substates are also observed during our simulation. Additionally, their facing backbones adopt the same conformations. These double facing BII substates are stable during the last 9 ns of the trajectories and result in a stretched DNA structure. On the other hand, protein-DNA contacts on 5' and 3' phosphodiester groups of the partner thymine of flipped adenine have changed. The sugar and phosphate parts of thymine have moved further into the empty space left by the flipping base without the influence of protein. Furthermore, readily high populated BII substates at the GpA step of palindromic tetrad TCGA rather than CpG step are observed in the free decamer. On the contrary, the BI substate at the GpA step is observed on the flipped adenine strand. A restrained MD simulation, reproducing the BI/BII pattern in the complex, demonstrated the influence of the unusual backbone conformation on the dynamical behavior of the target base. This finding along with the increased nearby interstrand phosphate distance is supportive to the N6-methylation mechanism.
Collapse
Affiliation(s)
- Fajar R Wibowo
- Institute of General, Inorganic, and Theoretical Chemistry, University of Innsbruck, Innrain 52a, A6020 Innsbruck, Austria
| | | | | | | |
Collapse
|
47
|
Fuxreiter M, Mezei M, Simon I, Osman R. Interfacial water as a "hydration fingerprint" in the noncognate complex of BamHI. Biophys J 2005; 89:903-11. [PMID: 15894630 PMCID: PMC1366640 DOI: 10.1529/biophysj.105.063263] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The molecular code of specific DNA recognition by proteins as a paradigm in molecular biology remains an unsolved puzzle primarily because of the subtle interplay between direct protein-DNA interaction and the indirect contribution from water and ions. Transformation of the nonspecific, low affinity complex to a specific, high affinity complex is accompanied by the release of interfacial water molecules. To provide insight into the conversion from the loose to the tight form, we characterized the structure and energetics of water at the protein-DNA interface of the BamHI complex with a noncognate sequence and in the specific complex. The fully hydrated models were produced with Grand Canonical Monte Carlo simulations. Proximity analysis shows that water distributions exhibit sequence dependent variations in both complexes and, in particular, in the noncognate complex they discriminate between the correct and the star site. Variations in water distributions control the number of water molecules released from a given sequence upon transformation from the loose to the tight complex as well as the local entropy contribution to the binding free energy. We propose that interfacial waters can serve as a "hydration fingerprint" of a given DNA sequence.
Collapse
Affiliation(s)
- Monika Fuxreiter
- Institute of Enzymology, Hungarian Academy of Sciences, Budapest.
| | | | | | | |
Collapse
|
48
|
Kumar A, Mishra PC, Suhai S. Adiabatic Electron Affinities of the Polyhydrated Adenine−Thymine Base Pair: A Density Functional Study. J Phys Chem A 2005; 109:3971-9. [PMID: 16833718 DOI: 10.1021/jp0456178] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Adiabatic electron affinities (AEAs) of the adenine-thymine (AT) base pair surrounded by 5 and 13 water molecules have been studied by density functional theory (DFT). Geometries of neutral AT x nH2O and anionic (AT x nH2O)- complexes (n = 5 and 13) were fully optimized, and vibrational frequency analysis was performed at the B3LYP/6-31+G** level of theory. The optimized structures of the neutral (AT x nH2O) and (AT x nH2O)- complexes were found to be somewhat nonplanar. Some of the water molecules are displaced away from the AT ring plane and linked with one another by hydrogen bonds. The optimized structures of the complexes are found to be in a satisfactory agreement with the observed experimental and molecular dynamics simulation results. In the optimized anionic complexes, the thymine (T) moiety was found to be puckered, whereas the adenine (A) moiety remained almost planar. Natural population analysis (NPA) performed using the B3LYP/6-31+G** method shows that the thymine moiety in the anionic (AT x nH2O)- complexes (n = 5 and 13) has most of the excess electronic charge, i.e., approximately -0.87 and approximately -0.81 (in the unit of magnitude of the electronic charge), respectively. The zero-point energy corrected adiabatic electron affinities of the hydrated AT base pair were found to be positive both for n = 5 and 13 and have the values of 0.97 and 0.92 eV, respectively, which are almost three times the AEA of the AT base pair. The results show that the presence of water molecules appreciably enhances the EA of the base pair.
Collapse
Affiliation(s)
- Anil Kumar
- Department of Physics, Banaras Hindu University, Varanasi-221 005, India
| | | | | |
Collapse
|
49
|
Tikhomirova A, Taulier N, Chalikian TV. Energetics of nucleic acid stability: the effect of DeltaCP. J Am Chem Soc 2005; 126:16387-94. [PMID: 15600340 DOI: 10.1021/ja046387d] [Citation(s) in RCA: 49] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report high-resolution differential scanning calorimetric data on the poly(dAdT)poly(dAdT), poly(dA)poly(dT), poly(dIdC)poly(dIdC), poly(dGdC)poly(dGdC), poly(rA)poly(rU), and poly(rI)poly(rC) nucleic acid duplexes. We use these data to evaluate the melting temperatures, TM, enthalpy changes, DeltaHM, and heat capacity changes, DeltaCP, accompanying helix-to-coil transitions of each polymeric duplex studied in this work at different NaCl concentrations. In agreement with previous reports, we have found that DeltaCP exhibits a positive, nonzero value, which, on average, equals 268 +/- 33 J mol(-1) K(-1). With DeltaCP, we have calculated the transition free energies, DeltaG, enthalpies, DeltaH, and entropies, DeltaS, for the duplexes as a function of temperature. Since, DeltaG, DeltaH, and DeltaS all strongly depend on temperature, the thermodynamic comparison between DNA and/or RNA duplexes (that may differ from one another with respect to sequence, composition, conformation, etc.) is physically meaningful only if extrapolated to a common temperature. We have performed such comparative analyses to derive differential thermodynamic parameters of formation of GC versus AT, AU, and IC base pairs as well as B' versus A and B helix conformations. We have proposed some general microscopic interpretations for the observed sequence-specific and conformation-specific thermodynamic differences between the duplexes.
Collapse
Affiliation(s)
- Anna Tikhomirova
- Department of Pharmaceutical Sciences, Leslie Dan Faculty of Pharmacy, University of Toronto, 19 Russell Street, Toronto, Ontario M5S 2S2, Canada
| | | | | |
Collapse
|
50
|
Wibowo FR, Trieb M, Rauch C, Wellenzohn B, Liedl KR. The N6-Methyl Group of Adenine Further Increases the BI Stability of DNA Compared to C5-Methyl Groups. J Phys Chem B 2004; 109:557-64. [PMID: 16851047 DOI: 10.1021/jp048519v] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Methylated DNA bases are natural modifications which play an important role in protein-DNA interactions. Recent experimental and theoretical results have shown an influence of the base modification on the conformational behavior of the DNA backbone. MD simulations of four different B-DNA dodecamers (d(GC)(6), d(AT)(6), d(G(5mCG)(5)C), and d(A(T6mA)(5)T)) have been performed with the aim to examine the influence of methyl groups on the B-DNA backbone behavior. An additional control simulation of d(AU)(6) has also been performed to examine the further influence of the C5-methyl group in thymine. Methyl groups in the major groove (as in C5-methylcytosine, thymine, or N6-methyladenine) decrease the BII substate population of RpY steps. Due to methylation a clearer distinction of the BI substate stability between YpR and RpY (CpG/GpC or TpA/ApT) steps arises. A positive correlation between the BII substate population and base stacking distances is seen only for poly(GC). A methyl group added into the major groove increases mean water residence times around the purine N7 atom, which may stabilize the BI substate by improving the hydration network between the DNA backbone and the major groove. The N6-methyl group also forms a water molecule bridge between the N6 and O4 atoms, and thus further stabilizes the BI substate.
Collapse
Affiliation(s)
- Fajar R Wibowo
- Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 52a, A-6020 Innsbruck, Austria
| | | | | | | | | |
Collapse
|