1
|
Affiliation(s)
- A. Subha Mahadevi
- Centre for Molecular Modelling, CSIR-Indian Institute of Chemical Technology, Tarnaka, Hyderabad, India 500607
| | - G. Narahari Sastry
- Centre for Molecular Modelling, CSIR-Indian Institute of Chemical Technology, Tarnaka, Hyderabad, India 500607
| |
Collapse
|
2
|
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
|
3
|
Banáš P, Mládek A, Otyepka M, Zgarbová M, Jurečka P, Svozil D, Lankaš F, Šponer J. Can We Accurately Describe the Structure of Adenine Tracts in B-DNA? Reference Quantum-Chemical Computations Reveal Overstabilization of Stacking by Molecular Mechanics. J Chem Theory Comput 2012; 8:2448-60. [PMID: 26588974 DOI: 10.1021/ct3001238] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Sequence-dependent local variations of helical parameters, structure, and flexibility are crucial for molecular recognition processes involving B-DNA. A-tracts, i.e., stretches of several consecutive adenines in one strand that are in phase with the DNA helical repeat, mediate significant DNA bending. During the past few decades, there have been intense efforts to understand the sequence dependence of helical parameters in DNA. Molecular dynamics (MD) simulations can provide valuable insights into the molecular mechanism behind the relationship between sequence and structure. However, although recent improvements in empirical force fields have helped to capture many sequence-dependent B-DNA properties, several problems remain, such as underestimation of the helical twist and suspected underestimation of the propeller twist in A-tracts. Here, we employ reference quantum mechanical (QM) calculations, explicit solvent MD, and bioinformatics to analyze the underestimation of propeller twisting of A-tracts in simulations. Although we did not identify a straightforward explanation, we discovered two imbalances in the empirical force fields. The first was overestimation of stacking interactions accompanied by underestimation of base-pairing energy, which we attribute to anisotropic polarizabilities that are not reflected by the isotropic force fields. This may lead to overstacking with potentially important consequences for MD simulations of nucleic acids. The second observed imbalance was steric clash between A(N1) and T(N3) nitrogens of AT base pairs in force-field descriptions, resulting in overestimation of the AT pair stretch in MD simulations. We also substantially extend the available set of benchmark estimated CCSD(T)/CBS data for B-DNA base stacking and provide a code that allows the generation of diverse base-stacking geometries suitable for QM computations with predefined intra- and interbase pair parameters.
Collapse
Affiliation(s)
- 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
| | - Arnošt Mládek
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, 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
| | - 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
| | - Daniel Svozil
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic.,Laboratory of Informatics and Chemistry, Faculty of Chemical Technology, Institute of Chemical Technology, Technicka 5, 166 28 Prague, Czech Republic
| | - Filip Lankaš
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nam. 6, 166 10 Prague, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, 625 00 Brno, Czech Republic
| |
Collapse
|
4
|
Rutledge LR, Wetmore SD. A computational proposal for the experimentally observed discriminatory behavior of hypoxanthine, a weak universal nucleobase. Phys Chem Chem Phys 2012; 14:2743-53. [PMID: 22270716 DOI: 10.1039/c2cp23600a] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
A computational model composed of six nucleobases was used to investigate why hypoxanthine does not yield duplexes of equal stability when paired opposite each of the natural DNA nucleobases. The magnitudes of all nearest-neighbor interactions in a DNA helix were calculated, including hydrogen-bonding, intra- and interstrand stacking interactions, as well as 1-3 intrastrand stacking interactions. Although the stacking interactions in DNA relevant arrangements are significant and account for at least one third of the total stabilization energy in our nucleobase complexes, the trends in the magnitude of the stacking interactions cannot explain the relative experimental melting temperatures previously reported in the literature. Furthermore, although the total hydrogen-bonding interactions explain why hypoxanthine preferentially pairs with cytosine, the experimental trend for the remaining nucleobases (A, T, G) is not explained. In fact, the calculated pairing preference of hypoxanthine matches that determined experimentally only when the sum of all types of nearest-neighbor interactions is considered. This finding highlights a strong correlation between the relative magnitude of the total nucleobase-nucleobase interactions and measured melting temperatures for DNA strands containing hypoxanthine despite the potential role of other factors (including hydration, temperature, sugar-phosphate backbone). By considering a large range of sequence combinations, we reveal that the binding preference of hypoxanthine is strongly dependent on the nucleobase sequence, which may explain the varied ability of hypoxanthine to universally bind to the natural bases. As a result, we propose that future work should closely examine the interplay between the dominant nucleobase-nucleobase interactions and the overall strand stability to fully understand how sequence context affects the universal binding properties of modified bases and to aid the design of new molecules with ambiguous pairing properties.
Collapse
Affiliation(s)
- Lesley R Rutledge
- Department of Chemistry & Biochemistry, University of Lethbridge, 4401 University Drive, Lethbridge, Alberta, Canada T1K 3M4
| | | |
Collapse
|