1
|
An insight into interaction of the uracil, thymine and cytosine biomolecules with methimazole anti-thyroid drug: DFT and GD3‑DFT approaches. Struct Chem 2022. [DOI: 10.1007/s11224-022-02059-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
2
|
Vecerova A, Triskova I, Trnkova L. The oxidation double peak of reduced guanine residues in short oligodeoxynucleotides: A study of its origin. J Electroanal Chem (Lausanne) 2018. [DOI: 10.1016/j.jelechem.2018.03.013] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
3
|
Nakarada Đ, Etinski M, Petković M. Using Density Functional Theory To Study Neutral and Ionized Stacked Thymine Dimers. J Phys Chem A 2016; 120:7704-7713. [PMID: 27626138 DOI: 10.1021/acs.jpca.6b06493] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Stacking interactions in thymine dimers are studied with density functional theory. According to our calculations, six dimers of comparable stability can be prepared at low temperature, but dimerization is impossible at room temperature due to the large entropy contribution that accompanies it. Analysis of vibrational anharmonic coupling terms shows that each of the dimers exhibits distinct vibrational dynamics. Properties of electron density in the intermolecular region are used to analyze neutral stacked species and their ionized forms. Bond paths and critical points in the intermolecular region are identified, but a simple relationship between binding energy and total electron density in the intermolecular critical points could not be found due to an uneven electron distribution in the binding region. The reduced density gradient was confirmed to be a useful tool for analysis of weak stacking interactions. Those interactions also affect vertical and adiabatic ionization energies, which are computed to be slightly lower for the dimers compared to the monomer.
Collapse
Affiliation(s)
- Đura Nakarada
- Faculty of Physical Chemistry, University of Belgrade , Studentski trg 12-16, 11 158 Belgrade, Serbia
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade , Studentski trg 12-16, 11 158 Belgrade, Serbia
| | - Milena Petković
- Faculty of Physical Chemistry, University of Belgrade , Studentski trg 12-16, 11 158 Belgrade, Serbia
| |
Collapse
|
4
|
Interplay between π···π stacking and cation···π interaction: a theoretical NMR study. JOURNAL OF THE IRANIAN CHEMICAL SOCIETY 2015. [DOI: 10.1007/s13738-015-0663-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
5
|
Cysewski P, Oliński R. Structural, electronic and energetic consequences of epigenetic cytosine modifications. Phys Chem Chem Phys 2015; 17:19616-24. [PMID: 26151626 DOI: 10.1039/c5cp02188j] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The hydrogen bonding patterns of cytosine and its seven C5-modifed analogues paired with canonical guanine were studied using the first principle approach. Both global minima and biologically relevant conformations were studied. The former resulted from full gradient geometry optimizations of hydrogen bonded pairs, while the latter were obtained based on 125 d(GpC) dinucleotides found in the PDB database. The obtained energetic, electronic and structural data lead to the conclusion that the epigenetically relevant modification of cytosine may have serious consequences on hydrogen bonding with guanine. First of all, the significant substituent effects were observed for such trends as charges on sites involved in hydrogen bonding, the total intermolecular interaction energy or electron densities at bond critical points. Moreover, the molecular orbital polarization contribution resulting from energy decomposition expressed in terms of absolutely localized molecular orbitals exhibited an inverse linear correlation with frozen density contributions. A substituent effect on the amount of charge transfer from pyrimidine toward guanine was also observed. The increase of intermolecular interactions of guanine with modified cytosine is associated with the increase of the electro-donating character of the C5-substituent. However, only pairs involving 5-methylcytosine are more stable than those formed by canonical cytosine. Furthermore, the energy differences observed for global minima also remain important for a broad range of displacement and angular parameters defining pair conformations in model d(GpC) dinucleotides. Due to the sensitivities of intermolecular interactions to mutual arrangements of monomers the modification of cytosine at the C5 site can significantly alter the actual energy profiles. Consequently, it may be anticipated that the modified dinucleotides will adopt different conformations than a standard G-C pair in a B-DNA double helix.
Collapse
Affiliation(s)
- P Cysewski
- Collegium Medicum of Bydgoszcz, Nicolaus Copernicus University in Toruń, Physical Chemistry Department, Kurpińskiego 5, Bydgoszcz, Poland.
| | | |
Collapse
|
6
|
Mládek A, Banáš P, Jurečka P, Otyepka M, Zgarbová M, Šponer J. Energies and 2'-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations. J Chem Theory Comput 2013; 10:463-80. [PMID: 26579924 DOI: 10.1021/ct400837p] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Sugar-phosphate backbone is an electronically complex molecular segment imparting RNA molecules high flexibility and architectonic heterogeneity necessary for their biological functions. The structural variability of RNA molecules is amplified by the presence of the 2'-hydroxyl group, capable of forming multitude of intra- and intermolecular interactions. Bioinformatics studies based on X-ray structure database revealed that RNA backbone samples at least 46 substates known as rotameric families. The present study provides a comprehensive analysis of RNA backbone conformational preferences and 2'-hydroxyl group orientations. First, we create a benchmark database of estimated CCSD(T)/CBS relative energies of all rotameric families and test performance of dispersion-corrected DFT-D3 methods and molecular mechanics in vacuum and in continuum solvent. The performance of the DFT-D3 methods is in general quite satisfactory. The B-LYP-D3 method provides the best trade-off between accuracy and computational demands. B3-LYP-D3 slightly outperforms the new PW6B95-D3 and MPW1B95-D3 and is the second most accurate density functional of the study. The best agreement with CCSD(T)/CBS is provided by DSD-B-LYP-D3 double-hybrid functional, although its large-scale applications may be limited by high computational costs. Molecular mechanics does not reproduce the fine energy differences between the RNA backbone substates. We also demonstrate that the differences in the magnitude of the hyperconjugation effect do not correlate with the energy ranking of the backbone conformations. Further, we investigated the 2'-hydroxyl group orientation preferences. For all families, we conducted a QM and MM hydroxyl group rigid scan in gas phase and solvent. We then carried out set of explicit solvent MD simulations of folded RNAs and analyze 2'-hydroxyl group orientations of different backbone families in MD. The solvent energy profiles determined primarily by the sugar pucker match well with the distribution data derived from the simulations. The QM and MM energy profiles predict the same 2'-hydroxyl group orientation preferences. Finally, we demonstrate that the high energy of unfavorable and rarely sampled 2'-hydroxyl group orientations can be attributed to clashes between occupied orbitals.
Collapse
Affiliation(s)
- Arnošt Mládek
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC, Central European Institute of Technology , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
7
|
Copeland KL, Pellock SJ, Cox JR, Cafiero ML, Tschumper GS. Examination of tyrosine/adenine stacking interactions in protein complexes. J Phys Chem B 2013; 117:14001-8. [PMID: 24171662 DOI: 10.1021/jp408027j] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
The π-stacking interactions between tyrosine amino acid side chains and adenine-bearing ligands are examined. Crystalline protein structures from the protein data bank (PDB) exhibiting face-to-face tyrosine/adenine arrangements were used to construct 20 unique 4-methylphenol/N9-methyladenine (p-cresol/9MeA) model systems. Full geometry optimization of the 20 crystal structures with the M06-2X density functional theory method identified 11 unique low-energy conformations. CCSD(T) complete basis set (CBS) limit interaction energies were estimated for all of the structures to determine the magnitude of the interaction between the two ring systems. CCSD(T) computations with double-ζ basis sets (e.g., 6-31G*(0.25) and aug-cc-pVDZ) indicate that the MP2 method overbinds by as much as 3.07 kcal mol(-1) for the crystal structures and 3.90 kcal mol(-1) for the optimized structures. In the 20 crystal structures, the estimated CCSD(T) CBS limit interaction energy ranges from -4.00 to -6.83 kcal mol(-1), with an average interaction energy of -5.47 kcal mol(-1), values remarkably similar to the corresponding data for phenylalanine/adenine stacking interactions. Geometry optimization significantly increases the interaction energies of the p-cresol/9MeA model systems. The average estimated CCSD(T) CBS limit interaction energy of the 11 optimized structures is 3.23 kcal mol(-1) larger than that for the 20 crystal structures.
Collapse
Affiliation(s)
- Kari L Copeland
- Department of Chemistry and Biochemistry, University of Mississippi , University, Mississippi 38677, United States
| | | | | | | | | |
Collapse
|
8
|
Mukherjee S, Kailasam S, Bansal M, Bhattacharyya D. Energy hyperspace for stacking interaction inAU/AUdinucleotide step: Dispersion-corrected density functional theory study. Biopolymers 2013; 101:107-20. [DOI: 10.1002/bip.22289] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2013] [Revised: 05/13/2013] [Accepted: 05/13/2013] [Indexed: 12/20/2022]
Affiliation(s)
- Sanchita Mukherjee
- Biophysics Division; Saha Institute of Nuclear Physics; Kolkata 700064 India
| | - Senthilkumar Kailasam
- Molecular Biophysics Unit; Indian Institute of Science; Bangalore 560012 Karnataka India
| | - Manju Bansal
- Molecular Biophysics Unit; Indian Institute of Science; Bangalore 560012 Karnataka India
| | | |
Collapse
|
9
|
Capobianco A, Caruso T, Celentano M, D'Ursi AM, Scrima M, Peluso A. Stacking interactions between adenines in oxidized oligonucleotides. J Phys Chem B 2013; 117:8947-53. [PMID: 23837863 DOI: 10.1021/jp404133a] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The effects of stacking interactions on the oxidation potentials of single strand oligonucleotides containing up to four consecutive adenines, alternated with thymines and cytosines in different sequences and ratios, have been determined by means of differential pulse voltammetry. Voltammetric measurements point toward the establishment in solution of structured oligonucleotide conformations, in which the nucleobases are well stacked altogether. Molecular dynamics simulations confirm that finding, indicating that single strands assume geometrical parameters characteristic of the B-DNA form. The analysis of the voltammetric signals in terms of a simple effective tight binding quantum model leads one to infer a robust set of parameters for treating hole transfer in one-electron-oxidized DNA containing adenines and thymines.
Collapse
Affiliation(s)
- Amedeo Capobianco
- Dipartimento di Chimica e Biologia, Università di Salerno, I-84084 Fisciano, Salerno, Italy
| | | | | | | | | | | |
Collapse
|
10
|
Villani G. Theoretical investigation of the coupling between hydrogen-atom transfer and stacking interaction in adenine-thymine dimers. Chemphyschem 2013; 14:1256-63. [PMID: 23494877 DOI: 10.1002/cphc.201200971] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2012] [Revised: 02/07/2013] [Indexed: 12/13/2022]
Abstract
Three different dimers of the adenine-thymine (A-T) base pair are studied to point out the changes of important properties (structure, atomic charge, energy and so on) induced by coupling between the movement of the atoms in the hydrogen bonds and the stacking interaction. The comparison of these results with those for the A-T monomer system explains the role of the stacking interaction in the hydrogen-atom transfer in this biologically important base pair. The results support the idea that this coupling depends on the exact dimer considered and is different for the N-N and N-O hydrogen bonds. In particular, the correlation between the hydrogen transfer and the stacking interaction is more relevant for the N-N bridge than for the N-O one. Also, the two different mechanisms of two-hydrogen transfer (step by step and concerted) can be modified by the stacking interaction between the base pairs.
Collapse
Affiliation(s)
- Giovanni Villani
- Istituto di Chimica dei Composti OrganoMetallici, ICCOM-UOS Pisa, Area della Ricerca del CNR, Via G. Moruzzi 1, 56124 Pisa, Italy.
| |
Collapse
|
11
|
Villani G. Theoretical investigation of the coupling between hydrogen atoms transfer and stacking interaction in guanine–cytosine dimers. Phys Chem Chem Phys 2013; 15:19242-52. [DOI: 10.1039/c3cp52855c] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
12
|
Johnson CA, Bloomingdale RJ, Ponnusamy VE, Tillinghast CA, Znosko BM, Lewis M. Computational model for predicting experimental RNA and DNA nearest-neighbor free energy rankings. J Phys Chem B 2011; 115:9244-51. [PMID: 21619071 DOI: 10.1021/jp2012733] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Hydrogen-bonding, intrastrand base-stacking, and interstrand base-stacking energies were calculated for RNA and DNA dimers at the MP2(full)/6-311G** level of theory. Standard A-form RNA and B-form DNA geometries from average fiber diffraction data were employed for all base monomer and dimer geometries, and all dimer binding energies were obtained via single-point calculations. The effects of water solvation were considered using the PCM model. The resulting dimer binding energies were used to calculate the 10 unique RNA and 10 unique DNA computational nearest-neighbor energies, and the ranking of these computational nearest neighbor energies are in excellent agreement with the ranking of the experimental nearest-neighbor free energies. These results dispel the notion that average fiber diffraction geometries are insufficient for calculating RNA and DNA stacking energies.
Collapse
Affiliation(s)
- Charles A Johnson
- Department of Chemistry, Saint Louis University, 3501 Laclede Avenue, Saint Louis, Missouri 63103, USA
| | | | | | | | | | | |
Collapse
|
13
|
Czyżnikowska Ż, Góra RW, Zaleśny R, Lipkowski P, Jarzembska KN, Dominiak PM, Leszczynski J. Structural Variability and the Nature of Intermolecular Interactions in Watson−Crick B-DNA Base Pairs. J Phys Chem B 2010; 114:9629-44. [DOI: 10.1021/jp101258q] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ż. Czyżnikowska
- Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland; Interdisciplinary Center for Nanotoxicity, Department of Chemistry, Jackson State University, 1400 J. R. Lynch St., Jackson, Mississippi 39217; and Department of Chemistry, Warsaw University, Pasteura 1, 02-093 Warsaw, Poland
| | - R. W. Góra
- Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland; Interdisciplinary Center for Nanotoxicity, Department of Chemistry, Jackson State University, 1400 J. R. Lynch St., Jackson, Mississippi 39217; and Department of Chemistry, Warsaw University, Pasteura 1, 02-093 Warsaw, Poland
| | - R. Zaleśny
- Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland; Interdisciplinary Center for Nanotoxicity, Department of Chemistry, Jackson State University, 1400 J. R. Lynch St., Jackson, Mississippi 39217; and Department of Chemistry, Warsaw University, Pasteura 1, 02-093 Warsaw, Poland
| | - P. Lipkowski
- Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland; Interdisciplinary Center for Nanotoxicity, Department of Chemistry, Jackson State University, 1400 J. R. Lynch St., Jackson, Mississippi 39217; and Department of Chemistry, Warsaw University, Pasteura 1, 02-093 Warsaw, Poland
| | - K. N. Jarzembska
- Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland; Interdisciplinary Center for Nanotoxicity, Department of Chemistry, Jackson State University, 1400 J. R. Lynch St., Jackson, Mississippi 39217; and Department of Chemistry, Warsaw University, Pasteura 1, 02-093 Warsaw, Poland
| | - P. M. Dominiak
- Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland; Interdisciplinary Center for Nanotoxicity, Department of Chemistry, Jackson State University, 1400 J. R. Lynch St., Jackson, Mississippi 39217; and Department of Chemistry, Warsaw University, Pasteura 1, 02-093 Warsaw, Poland
| | - J. Leszczynski
- Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland; Interdisciplinary Center for Nanotoxicity, Department of Chemistry, Jackson State University, 1400 J. R. Lynch St., Jackson, Mississippi 39217; and Department of Chemistry, Warsaw University, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
14
|
Svozil D, Hobza P, Sponer J. Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations? J Phys Chem B 2010; 114:1191-203. [PMID: 20000584 DOI: 10.1021/jp910788e] [Citation(s) in RCA: 92] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
High level ab initio methods have been used to study stacking interactions in ten unique base pair steps both in A-RNA and in B-DNA duplexes. The protocol for selection of geometries based on molecular dynamics (MD) simulations is proposed, and its suitability is demonstrated by comparison with stacking in steps at fiber diffraction geometries. It is shown that fiber diffraction geometries are not sufficiently accurate for interaction energy calculations. In addition, the protocol for selection of geometries based on MD simulations allows for the evaluation of the variability of the intrinsic stacking energies along the MD trajectories. The uncertainty in stacking energies (difference between the most and least stable geometry) due to the dynamical nature of systems can be, in some cases, as large as 3.0 kcal x mol(-1), which is almost 50% of the actual sequence dependence of base stacking energies (the energy difference between the most and least stable sequences). Thus, assessing the relative magnitude of the gas phase stacking energy using a single geometry for each sequence is insufficient to obtain an unambiguous order of gas phase stacking energies in canonical double helices. Though the ordering of ten unique dinucleotide steps cannot be definitive, some general conclusions were drawn. The stacking energies of base pair steps in A-RNA are more evenly separated compared to B-DNA, and their ordering is less sensitive to the dynamics of the system compared to be B-DNA. The most stable step both in B-DNA and A-RNA is the GC/GC [corrected] step that is well separated from the second most stable step CG/CG. [corrected] Also the least stable step (the CC/GG step) is well separated from the rest of the structures. The calculations further show that B-DNA stacking is favorable only marginally (on average by 1.14 kcal x mol(-1) per base pair step) over A-RNA stacking, and this difference vanishes after subtracting the stabilizing van der Waals effect of the thymine 5-methyl group that is absent in RNA. Basically, no correlation between the sequence dependence of gas phase stacking energies and the sequence dependence of DeltaG degrees(37) free energies used in nearest-neighbor models was found either for B-DNA or for A-RNA. This reflects the complexity of the balance of forces that are responsible for the sequence dependence of thermodynamics stability of nucleic acids, which masks the effect of the intrinsic interactions between the stacked base pairs.
Collapse
Affiliation(s)
- Daniel Svozil
- Faculty of Chemical Technology, Laboratory of Informatics and Chemistry, Institute of Chemical Technology, Technická 3, 166 28, Prague 6, Czech Republic.
| | | | | |
Collapse
|
15
|
Morgado CA, Jurečka P, Svozil D, Hobza P, Šponer J. Reference MP2/CBS and CCSD(T) quantum-chemical calculations on stacked adenine dimers. Comparison with DFT-D, MP2.5, SCS(MI)-MP2, M06-2X, CBS(SCS-D) and force field descriptions. Phys Chem Chem Phys 2010; 12:3522-34. [DOI: 10.1039/b924461a] [Citation(s) in RCA: 77] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
16
|
Ebrahimi A, Habibi-Khorassani M, Gholipour AR, Masoodi HR. Interaction between uracil nucleobase and phenylalanine amino acid: the role of sodium cation in stacking. Theor Chem Acc 2009. [DOI: 10.1007/s00214-009-0588-y] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
17
|
Czyznikowska Z. How does modification of adenine by hydroxyl radical influence the stability and the nature of stacking interactions in adenine-cytosine complex? J Mol Model 2009; 15:615-22. [PMID: 19198902 DOI: 10.1007/s00894-008-0447-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2008] [Accepted: 12/16/2008] [Indexed: 11/28/2022]
Abstract
This study reports on ab initio calculations of adenine-cytosine complexes in two different context alignments appearing in B-DNA. The influence of adenine modification by hydroxyl radical on the stability of the complexes is also discussed. The analysis was performed on over 40 crystallographic structures for each of the sequence contexts. In most cases, modification of adenine by hydroxyl radical leads to less negative intermolecular interaction energies. The issue of the influence of alteration of structural base step parameters on the stability of modified and unmodified adenine-cytosine complexes is also addressed. Analysis of the dependence of intermolecular interaction energy on base step parameters reveals that for twist and shift modification of adenine by hydroxyl radical leads to quite different interaction energy profiles in comparison with unmodified complexes. In order to elucidate the physical origins of this phenomenon, i.e. to analyze how the modification of adenine by hydroxyl radical is reflected in the change of intermolecular interaction energy components, a variational-perturbational decomposition scheme was applied at the MP2/aug-cc-pVDZ level of theory.
Collapse
Affiliation(s)
- Zaneta Czyznikowska
- Faculty of Pharmacy, Collegium Medicum, Nicolaus Copernicus University, Jagiellonska 13-15, 85-067 Bydgoszcz, Poland.
| |
Collapse
|
18
|
On the importance of electrostatics in stabilization of stacked guanine–adenine complexes appearing in B-DNA crystals. ACTA ACUST UNITED AC 2009. [DOI: 10.1016/j.theochem.2008.10.040] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
19
|
Structural and energetic heterogeneities of canonical and oxidized central guanine triad of B-DNA telomeric fragments. J Mol Model 2009; 15:607-13. [PMID: 19132417 DOI: 10.1007/s00894-008-0438-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2008] [Accepted: 12/05/2008] [Indexed: 10/21/2022]
Abstract
The intermolecular interaction energies in central guanine triad of telomeric B-DNA were estimated based on ab initio quantum chemistry calculations on the MP2/aDZ level of theory. The source of structural information was molecular dynamics simulation of both canonical (AGGGTT) and oxidized (AG8oxoGGTT) telomere units. Our calculations demonstrate that significant stiffness of central triad occurs if 8oxoG is present. The origin of such feature is mainly due to the increase of stacking interactions of 8oxoG with neighbouring guanine molecules and stronger hydrogen bonding formation of 8oxoG with cytosine if compared with canonical guanine. Another interesting observation is the context independence of stacking interactions of 8oxoG. Unlike to 5'-G2/G3-3' and 5'-G3/G4-3' sequences which are energetically different, 5'-G2/8oxoG3-3' and 5'-8oxoG3/G4-3' sequences are almost iso-energetic.
Collapse
|
20
|
Cysewski P. Intra-strand stacking interactions in B-DNA crystals characterized by post-SCF quantum chemistry computations. NEW J CHEM 2009. [DOI: 10.1039/b909327c] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
21
|
The post-SCF quantum chemistry characteristics of inter- and intra-strand stacking interactions in d(CpG) and d(GpC) steps found in B-DNA, A-DNA and Z-DNA crystals. J Mol Model 2008; 15:597-606. [PMID: 19039609 DOI: 10.1007/s00894-008-0378-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2008] [Accepted: 10/31/2008] [Indexed: 10/21/2022]
Abstract
The energies of intra- and inter-strand stacking interactions in model d(GpC) and d(CpG) two-base-pair steps were estimated by MP2/aug-cc-pVDZ single point calculations corrected for basis superposition errors. The stacked two-nucleobase pairs were constructed using experimental values of base pair and base step parameters taken from Nucleic Acid Database (http://ndbserver.rutgers.edu/). Three distinct polymorphic forms were analysed, namely A-, B- and Z-DNA. The applied methodology enables statistical analysis of structural and energetic diversities. The structural relationships between polymorphic forms are quite complex and depend on the sequence of pairs. The variability of parameters such as shift and tilt is almost the same irrespective of the polymorphic form and sequence of steps analysed. In contrast, shift and twist distributions easily discriminate all three polymorphic forms of DNA. Interestingly, despite significant structural diversities, the energies of the most frequent energy ranges are comparable irrespective of the polymorphic form and base sequence. There was observed compensation of inter- and intra-strand interactions, especially for d(GpC) and d(CpG) steps found in A- and B-DNA. Thus, among many other roles, these pairs act as a kind of energetic buffer, balancing the double helix.
Collapse
|
22
|
Czyznikowska Z, Zaleśny R. Theoretical insights into the nature of intermolecular interactions in cytosine dimer. Biophys Chem 2008; 139:137-43. [PMID: 19038487 DOI: 10.1016/j.bpc.2008.11.001] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2008] [Revised: 11/03/2008] [Accepted: 11/04/2008] [Indexed: 10/21/2022]
Abstract
In this study we discuss stacking interactions in cytosine dimer in conformations appearing in B-DNA crystals. The variational-perturbational scheme was applied for decomposition of the intermolecular interaction energy at the MP2 level of theory. The significant influence of the mutual orientation of cytosine monomers was observed not only on the total intermolecular interaction energy but also on its components: Different components of intermolecular interaction energy depend in different manner on parameters describing mutual orientation of cytosine monomers.
Collapse
Affiliation(s)
- Zaneta Czyznikowska
- Institute of Organic and Pharmaceutical Chemistry, The National Hellenic Research Foundation, 48 Vas. Constantinou Avenue, 11635 Athens, Greece
| | | |
Collapse
|
23
|
Řezáč J, Jurečka P, Riley KE, Černý J, Valdes H, Pluháčková K, Berka K, Řezáč T, Pitoňák M, Vondrášek J, Hobza P. Quantum Chemical Benchmark Energy and Geometry Database for Molecular Clusters and Complex Molecular Systems (www.begdb.com): A Users Manual and Examples. ACTA ACUST UNITED AC 2008. [DOI: 10.1135/cccc20081261] [Citation(s) in RCA: 129] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Our previous benchmark CCSD(T)/ complete basis set limit calculations were collected into a database named begdb - Benchmark Energy and Geometry DataBase. Web-based interface to this database was prepared and is available at www.begdb.com. Users can browse, search and plot the data online or download structures and energy tables.
Collapse
|