1
|
Jayachandran P, Angamuthu A, Gopalan P. Quantum Chemical Study on the Structure and Energetics of Binary Ionic Porphyrin Complexes. J CHIN CHEM SOC-TAIP 2018. [DOI: 10.1002/jccs.201700383] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
| | - Abiram Angamuthu
- Department of Physics; Karunya Institute of Technology and Sciences; Coimbatore 641114 India
| | - Praveena Gopalan
- Department of Physics; PSGR Krishnammal College for Women; Coimbatore 641004 India
| |
Collapse
|
2
|
Chawla M, Credendino R, Chermak E, Oliva R, Cavallo L. Theoretical Characterization of the H-Bonding and Stacking Potential of Two Nonstandard Nucleobases Expanding the Genetic Alphabet. J Phys Chem B 2016; 120:2216-24. [DOI: 10.1021/acs.jpcb.6b00125] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Mohit Chawla
- Physical
Sciences and Engineering Division (PSE), Kaust Catalysis Center (KCC), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Raffaele Credendino
- Physical
Sciences and Engineering Division (PSE), Kaust Catalysis Center (KCC), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Edrisse Chermak
- Physical
Sciences and Engineering Division (PSE), Kaust Catalysis Center (KCC), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Romina Oliva
- Department
of Sciences and Technologies, University Parthenope of Naples, Centro Direzionale Isola C4, I-80143 Naples, Italy
| | - Luigi Cavallo
- Physical
Sciences and Engineering Division (PSE), Kaust Catalysis Center (KCC), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
3
|
Cybulski H, Sadlej J. Symmetry-Adapted Perturbation-Theory Interaction-Energy Decomposition for Hydrogen-Bonded and Stacking Structures. J Chem Theory Comput 2015; 4:892-7. [PMID: 26621230 DOI: 10.1021/ct800067m] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
This letter reports the computational ab initio studies on the stacked and hydrogen-bonded geometries of the uracil dimer and pyrimidine···p-benzoquinone complex with a special regard to the ratios of different interaction-energy terms calculated by means of the symmetry-adapted perturbation theory (SAPT). In the hydrogen-bonded systems the absolute value of the dispersion term constitutes approximately half of the absolute value of the total SAPT0 interaction energy, while in the stacking complexes the ratio of the dispersion to the total interaction energy is much larger, ca. 1.2-2.0. Our SAPT results are compared with the DFT-SAPT results published recently by the Hobza group (J. Chem. Phys. 2007, 127, 075104), and the role of the dispersion contribution in stacking and hydrogen-bonded arrangements is discussed. The methodological part of this letter presents the influence of counterpoise corrections in the optimization procedure on the geometries of the systems and the calculated SAPT contributions.
Collapse
Affiliation(s)
- Hubert Cybulski
- Department of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Joanna Sadlej
- Department of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
4
|
Kruse H, Mladek A, Gkionis K, Hansen A, Grimme S, Sponer J. Quantum chemical benchmark study on 46 RNA backbone families using a dinucleotide unit. J Chem Theory Comput 2015; 11:4972-91. [PMID: 26574283 DOI: 10.1021/acs.jctc.5b00515] [Citation(s) in RCA: 80] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
We have created a benchmark set of quantum chemical structure-energy data denoted as UpU46, which consists of 46 uracil dinucleotides (UpU), representing all known 46 RNA backbone conformational families. Penalty-function-based restrained optimizations with COSMO TPSS-D3/def2-TZVP ensure a balance between keeping the target conformation and geometry relaxation. The backbone geometries are close to the clustering-means of their respective RNA bioinformatics family classification. High-level wave function methods (DLPNO-CCSD(T) as reference) and a wide-range of dispersion-corrected or inclusive DFT methods (DFT-D3, VV10, LC-BOP-LRD, M06-2X, M11, and more) are used to evaluate the conformational energies. The results are compared to the Amber RNA bsc0χOL3 force field. Most dispersion-corrected DFT methods surpass the Amber force field significantly in accuracy and yield mean absolute deviations (MADs) for relative conformational energies of ∼0.4-0.6 kcal/mol. Double-hybrid density functionals represent the most accurate class of density functionals. Low-cost quantum chemical methods such as PM6-D3H+, HF-3c, DFTB3-D3, as well as small basis set calculations corrected for basis set superposition errors (BSSEs) by the gCP procedure are also tested. Unfortunately, the presently available low-cost methods are struggling to describe the UpU conformational energies with satisfactory accuracy. The UpU46 benchmark is an ideal test for benchmarking and development of fast methods to describe nucleic acids, including force fields.
Collapse
Affiliation(s)
- Holger Kruse
- 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
| | - Arnost Mladek
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Konstantinos Gkionis
- 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
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn , Beringstr. 4, D-53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn , Beringstr. 4, D-53115 Bonn, Germany
| | - Jiri Sponer
- 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
|
5
|
Anticancer drug IUdR and other 5-halogen derivatives of 2′-deoxyuridine: conformers, hydrates, and structure–activity relationships. Struct Chem 2013. [DOI: 10.1007/s11224-013-0225-1] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
6
|
Sharma P, Lait LA, Wetmore SD. yDNA versus yyDNA pyrimidines: computational analysis of the effects of unidirectional ring expansion on the preferred sugar-base orientation, hydrogen-bonding interactions and stacking abilities. Phys Chem Chem Phys 2013; 15:2435-48. [PMID: 23303174 DOI: 10.1039/c2cp43910g] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The properties of natural, y- and yy-pyrimidines are compared using computational (B3LYP, MP2) methods. Ring expansion upon incorporation of benzene or naphthalene into the natural pyrimidines affects the preferred orientation of the base about the glycosidic bond in the corresponding nucleoside to a similar extent. Specifically, although the natural pyrimidines preferentially adopt the anti orientation with respect to the 2'-deoxyribose moiety, the expanded analogues will likely display (anti/syn) conformational heterogeneity, which may lead to alternate hydrogen-bonding modes in double-stranded duplexes. Nevertheless, the A:T Watson-Crick hydrogen-bond strengths do not significantly change upon base expansion, while the G:C interaction energy is slightly strengthened upon incorporation of either expanded pyrimidine. The largest effect of base expansion occurs in the stacking energies. Specifically, the maximum (most negative) stacking energies in isolated dimers formed by aligning the nucleobase centers of mass can be increased up to 45% by inclusion of a single y-pyrimidine and up to 55% by consideration of a yy-pyrimidine. Similar increases in the stacking interactions are found when a simplified duplex model composed of two stacked (hydrogen-bonded) base pairs is considered, where both the intrastrand and interstrand stacking interactions can be increased and the effects are more pronounced for the yy-pyrimidines. Moreover, the total stability (sum of all hydrogen-bonding and stacking interactions) is greater for duplexes containing expanded yy-pyrimidines compared to y-pyrimidines, which is mainly due to enhanced stacking interactions. Thus, our calculations suggest that multiple unidirectional increases in the size of the nucleobase spacer can continuously enhance the stability of expanded duplexes.
Collapse
Affiliation(s)
- Purshotam Sharma
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta, Canada T1K 3M4
| | | | | |
Collapse
|
7
|
Deepa P, Kolandaivel P, Senthilkumar K. Structural properties and the effect of platinum drugs with DNA base pairs. Struct Chem 2012. [DOI: 10.1007/s11224-012-0087-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
8
|
Hunter RS, van Mourik T. DNA base stacking: The stacked uracil/uracil and thymine/thymine minima. J Comput Chem 2012; 33:2161-72. [DOI: 10.1002/jcc.23052] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2012] [Revised: 06/01/2012] [Accepted: 06/04/2012] [Indexed: 01/24/2023]
|
9
|
Morgado CA, Svozil D, Turner DH, Šponer J. Understanding the role of base stacking in nucleic acids. MD and QM analysis of tandem GA base pairs in RNA duplexes. Phys Chem Chem Phys 2012; 14:12580-91. [PMID: 22722325 DOI: 10.1039/c2cp40556c] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Preceding NMR experiments show that the conformation of tandem GA base pairs, an important recurrent non-canonical building block in RNA duplexes, is context dependent. The GA base pairs adopt "sheared" N3(G)-N6(A), N2(G)-N7(A) geometry in the r(CGAG)(2) and r(iGGAiC)(2) contexts while switching to "imino" N1(G)-N1(A), O6(G)-N6(A) geometry in the r(GGAC)(2) and r(iCGAiG)(2) contexts (iC and iG stand for isocytosine and isoguanine, respectively). As base stacking is likely to be one of the key sources of the context dependence of the conformation of GA base pairs, we calculated base stacking energies in duplexes containing such base pairs, to see if this dependence can be predicted by stacking energy calculations. When investigating the context dependence of the GA geometry two different conformations of the same duplex were compared (imino vs. sheared). The geometries were generated via explicit solvent MD simulations of the respective RNA duplexes, while the subsequent QM energy calculations focused on base stacking interactions of the four internal base pairs. Geometrical relaxation of nucleobase atoms prior to the stacking energy computations has a non-negligible effect on the results. The stacking energies were derived at the DFT-D/6-311++G(3df,3pd) level. We show a rather good correspondence between the intrinsic gas-phase stacking energies and the NMR-determined GA geometries. The conformation with more favorable gas-phase stacking is in most cases the one observed in experiments. This correlation is not improved when including solvent effects via the COSMO method. On the other side, the stacking calculations do not predict the relative thermodynamic stability of duplex formation for different sequences.
Collapse
Affiliation(s)
- Claudio A Morgado
- Universidad Técnica Federico Santa María, Departamento de Química, Casilla 110-V, Valparaíso, Chile.
| | | | | | | |
Collapse
|
10
|
Sychrovský V, Sochorová Vokáčová Z, Trantírek L. Guanine bases in DNA G-quadruplex adopt nonplanar geometries owing to solvation and base pairing. J Phys Chem A 2012; 116:4144-51. [PMID: 22471881 DOI: 10.1021/jp2110049] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The effect of base pairing and solvation on pyramidalization of the glycosidic nitrogen found in the residues of parallel G-quadruplex with NDB ID UDF062 is analyzed and explained with theoretical calculations. The extent of the pyramidalization depends on the local geometry of the 2'-deoxyguanosine residues, namely on their glycosidic torsion and sugar pucker, which are predetermined by the 3D-architecture of G-quadruplex. Pyramidal inversion of the glycosidic nitrogen found in 2'-deoxyguanosines of G-quadruplex is induced owing to site-specifically coordinated solvent. Different adiabatic structural constraints used for fixing the base-to-sugar orientation of 2'-deoxyguanosine in geometry optimizations result in different extents of pyramidalization and induce pyramidal inversion of the glycosidic nitrogen. These model geometry constraints helped us analyze the effect of real constraints represented by explicit molecular environment of selected residues of the G-quadruplex. The maximal extent of the glycosidic nitrogen pyramidalization found in the high-resolution crystal structure corresponds to the calculation to deformation energy of only 1 kcal mol(-1). The out-of-plane deformations of nucleobases thus provide a way for compensating the site-specific external environmental stress on the G-quadruplex.
Collapse
Affiliation(s)
- Vladimír Sychrovský
- Institute of Organic Chemistry and Biochemistry, v.v.i., Academy of Sciences of the Czech Republic, Flemingovo square 2, 166 10 Prague 6, Czech Republic.
| | | | | |
Collapse
|
11
|
Nikolaienko TY, Bulavin LA, Hovorun DM. How Flexible are DNA Constituents? The Quantum-Mechanical Study. J Biomol Struct Dyn 2011; 29:563-75. [DOI: 10.1080/07391102.2011.10507406] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
12
|
|
13
|
Shishkin OV, Gorb L, Leszczynski J. Conformational flexibility of pyrimidine rings of nucleic acid bases in polar environment: PCM study. Struct Chem 2009. [DOI: 10.1007/s11224-009-9477-1] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
14
|
Morgado CA, Jurečka P, Svozil D, Hobza P, Šponer J. Balance of Attraction and Repulsion in Nucleic-Acid Base Stacking: CCSD(T)/Complete-Basis-Set-Limit Calculations on Uracil Dimer and a Comparison with the Force-Field Description. J Chem Theory Comput 2009; 5:1524-44. [DOI: 10.1021/ct9000125] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Claudio A. Morgado
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Prague 6, Czech Republic, and Department of Physical Chemistry, Palacky University, tr. Svobody 26, 771 46, Olomouc, Czech Republic
| | - Petr Jurečka
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Prague 6, Czech Republic, and Department of Physical Chemistry, Palacky University, tr. Svobody 26, 771 46, Olomouc, Czech Republic
| | - Daniel Svozil
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Prague 6, Czech Republic, and Department of Physical Chemistry, Palacky University, tr. Svobody 26, 771 46, Olomouc, Czech Republic
| | - Pavel Hobza
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Prague 6, Czech Republic, and Department of Physical Chemistry, Palacky University, tr. Svobody 26, 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, Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Prague 6, Czech Republic, and Department of Physical Chemistry, Palacky University, tr. Svobody 26, 771 46, Olomouc, Czech Republic
| |
Collapse
|
15
|
He X, Fusti-Molnar L, Cui G, Merz KM. Importance of dispersion and electron correlation in ab initio protein folding. J Phys Chem B 2009; 113:5290-300. [PMID: 19320454 PMCID: PMC2737261 DOI: 10.1021/jp8106952] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Dispersion is well-known to be important in biological systems, but the effect of electron correlation in such systems remains unclear. In order to assess the relationship between the structure of a protein and its electron correlation energy, we employed both full system Hartree-Fock (HF) and second-order Møller-Plesset perturbation (MP2) calculations in conjunction with the Polarizable Continuum Model (PCM) on the native structures of two proteins and their corresponding computer-generated decoy sets. Because of the expense of the MP2 calculation, we have utilized the fragment molecular orbital method (FMO) in this study. We show that the sum of the Hartree-Fock (HF) energy and force field (LJ6)-derived dispersion energy (HF + LJ6) is well correlated with the energies obtained using second-order Møller-Plesset perturbation (MP2) theory. In one of the two examples studied, the correlation energy as well as the empirical dispersive energy term was able to discriminate between native and decoy structures. On the other hand, for the second protein we studied, neither the correlation energy nor dispersion energy showed discrimination capabilities; however, the ab initio MP2 energy and the HF+LJ6 both ranked the native structure correctly. Furthermore, when we randomly scrambled the Lennard-Jones parameters, the correlation between the MP2 energy and the sum of the HF energy and dispersive energy (HF+LJ6) significantly drops, which indicates that the choice of Lennard-Jones parameters is important.
Collapse
Affiliation(s)
- Xiao He
- Department of Chemistry and the Quantum Theory Project, 2328 New Physics Building, P.O. Box 118435, University of Florida, Gainesville, Florida 32611-8435
| | - Laszlo Fusti-Molnar
- Department of Chemistry and the Quantum Theory Project, 2328 New Physics Building, P.O. Box 118435, University of Florida, Gainesville, Florida 32611-8435
| | - Guanglei Cui
- Department of Chemistry and the Quantum Theory Project, 2328 New Physics Building, P.O. Box 118435, University of Florida, Gainesville, Florida 32611-8435
| | - Kenneth M. Merz
- Department of Chemistry and the Quantum Theory Project, 2328 New Physics Building, P.O. Box 118435, University of Florida, Gainesville, Florida 32611-8435
| |
Collapse
|
16
|
Interaction of psoralens with DNA-bases (II): An ab initio quantum chemical, density functional theory and second-order MØller-Plesset perturbational study. ACTA ACUST UNITED AC 2009. [DOI: 10.1016/j.theochem.2008.10.012] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
|
17
|
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
|
18
|
The post-SCF quantum chemistry characteristics of the energetic heterogeneity of stacked guanine–guanine pairs found in B-DNA and A-DNA crystals. ACTA ACUST UNITED AC 2008. [DOI: 10.1016/j.theochem.2008.06.019] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
19
|
Wenska G, Filipiak P, Asmus KD, Bobrowski K, Koput J, Marciniak B. Formation of a Sandwich-Structure Assisted, Relatively Long-Lived Sulfur-Centered Three-Electron Bonded Radical Anion in the Reduction of a Bis(1-substituted-uracilyl) Disulfide in Aqueous Solution. J Phys Chem B 2008; 112:10045-53. [DOI: 10.1021/jp8041928] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Grazyna Wenska
- Faculty of Chemistry, Adam Mickiewicz University, 60-780 Poznan, Poland, and Institute of Nuclear Chemistry and Technology, 03-195 Warszawa, Poland
| | - Piotr Filipiak
- Faculty of Chemistry, Adam Mickiewicz University, 60-780 Poznan, Poland, and Institute of Nuclear Chemistry and Technology, 03-195 Warszawa, Poland
| | - Klaus-Dieter Asmus
- Faculty of Chemistry, Adam Mickiewicz University, 60-780 Poznan, Poland, and Institute of Nuclear Chemistry and Technology, 03-195 Warszawa, Poland
| | - Krzysztof Bobrowski
- Faculty of Chemistry, Adam Mickiewicz University, 60-780 Poznan, Poland, and Institute of Nuclear Chemistry and Technology, 03-195 Warszawa, Poland
| | - Jacek Koput
- Faculty of Chemistry, Adam Mickiewicz University, 60-780 Poznan, Poland, and Institute of Nuclear Chemistry and Technology, 03-195 Warszawa, Poland
| | - Bronislaw Marciniak
- Faculty of Chemistry, Adam Mickiewicz University, 60-780 Poznan, Poland, and Institute of Nuclear Chemistry and Technology, 03-195 Warszawa, Poland
| |
Collapse
|
20
|
Sponer J, Riley KE, Hobza P. Nature and magnitude of aromatic stacking of nucleic acid bases. Phys Chem Chem Phys 2008; 10:2595-610. [PMID: 18464974 DOI: 10.1039/b719370j] [Citation(s) in RCA: 270] [Impact Index Per Article: 16.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
This review summarises recent advances in quantum chemical calculations of base-stacking forces in nucleic acids. We explain in detail the very complex relationship between the gas-phase base-stacking energies, as revealed by quantum chemical (QM) calculations, and the highly variable roles of these interactions in nucleic acids. This issue is rarely discussed in quantum chemical and physical chemistry literature. We further extensively discuss methods that are available for base-stacking studies, complexity of comparison of stacking calculations with gas phase experiments, balance of forces in stacked complexes of nucleic acid bases, and the relation between QM and force field descriptions. We also review all recent calculations on base-stacking systems, including details analysis of the B-DNA stacking. Specific attention is paid to the highest accuracy QM calculations, to the decomposition of the interactions, and development of dispersion-balanced DFT methods. Future prospects of computational studies of base stacking are discussed.
Collapse
Affiliation(s)
- Jirí Sponer
- Institute of Organic Chemistry and Biochemistry, vvi, Academy of Sciences of the Czech Republic and Center for Biomolecules and Complex Molecular Systems, Prague 6, Czech Republic
| | | | | |
Collapse
|
21
|
Wang W, Wang N, Li P, Bu Y, Xie X, Song R. Theoretical studies on the properties of uracil and its dimer upon thioketo substitution. Theor Chem Acc 2008. [DOI: 10.1007/s00214-008-0442-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
22
|
Czyżnikowska Ż, Zaleśny R, Ziółkowski M, Gora RW, Cysewski P. The nature of interactions in uracil dimer: An ab initio study. Chem Phys Lett 2007. [DOI: 10.1016/j.cplett.2007.10.106] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
23
|
Rutledge LR, Campbell-Verduyn LS, Hunter KC, Wetmore SD. Characterization of nucleobase-amino acid stacking interactions utilized by a DNA repair enzyme. J Phys Chem B 2007; 110:19652-63. [PMID: 17004834 DOI: 10.1021/jp061939v] [Citation(s) in RCA: 67] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
The present work characterizes the gas-phase stacking interactions between four aromatic amino acid residues (histidine, phenylalanine, tyrosine, and tryptophan) and adenine or 3-methyladenine due to the proposed utilization of these interactions by enzymes that repair DNA alkylation damage. The MP2 potential energy surfaces of the stacked dimers are considered as a function of four variables (vertical displacement, angle of rotation, horizontal displacement, and tilt angle) using a variety of basis sets. It is found that the maximum stacking interaction energy decreases with the amino acid according to TRP > TYR approximately HIS > PHE for both nucleobases. However, the magnitude of the stacking interaction significantly increases upon alkylation (by 50-115%). Comparison of the stacking energies calculated using our surface scans to those estimated from experimental crystal structures indicates that the stacking interactions within the active site of 3-methyladenine DNA glycosylase can account for 65-75% of the maximum possible stacking interaction between the relevant molecules. The decrease in stacking in the crystal structure arises due to significant differences in the relative orientations of the nucleobase and amino acid. Nevertheless, alkylation is found to significantly increase the stacking energy when the crystal structure geometries are considered. Our calculations provide computational support for suggestions that alkylation enhances the stacking interactions within the active site of DNA repair enzymes, and they give a measure of the magnitude of this enhancement. Our results suggest that alkylation likely plays a more important role in substrate identification and removal than the nature of the aromatic amino acid that interacts with the substrate via stacking interactions.
Collapse
Affiliation(s)
- Lesley R Rutledge
- Department of Chemistry, Mount Allison University, 63C York Street, Sackville, New Brunswick, Canada E4L 1G8
| | | | | | | |
Collapse
|
24
|
Robertazzi A, Platts JA. Gas-phase DNA oligonucleotide structures. A QM/MM and atoms in molecules study. J Phys Chem A 2007; 110:3992-4000. [PMID: 16539422 DOI: 10.1021/jp056626z] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
QM/MM calculations have been employed to investigate the role of hydrogen bonding and pi-stacking in single- and double-stranded DNA oligonucleotides. DFT calculations and Atoms in Molecules analysis on QM/MM-optimized structures allow characterization and estimation of the energies of pi-stacking and hydrogen-bond interactions. This shows that pi-stacking interactions depend on the number and the nature of the DNA bases for single-stranded nucleotides; for instance, guanines are found to be involved in strong hydrogen bonds, whereas adenines interact mainly via stacking interactions. The role of interbase hydrogen bonding was explored: the -NH2 groups of guanine, adenine, and cytosine participate in N-H...O and N-H...N interactions. These are much stronger in single-strand oligonucleotides, where the -NH2 groups are highly nonplanar. In double-stranded DNA, the strong base-pairing hydrogen bonds of complementary bases lead to more planar -NH2 groups, which tend to be involved in pi-stacking interactions rather than H-bonds. The use of AIM also allows us to evaluate the interplay of pi-stacking and H-bonding, suggesting that cooperativity does occur, but is generally limited to about 1-2 kcal/mol.
Collapse
Affiliation(s)
- Arturo Robertazzi
- School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, UK
| | | |
Collapse
|
25
|
McConnell TL, Wetmore SD. How Do Size-Expanded DNA Nucleobases Enhance Duplex Stability? Computational Analysis of the Hydrogen-Bonding and Stacking Ability of xDNA Bases. J Phys Chem B 2007; 111:2999-3009. [PMID: 17388411 DOI: 10.1021/jp0670079] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Computational chemistry (B3LYP, MP2) is used to study the properties of size-expanded DNA nucleobases generated by inserting a benzene spacer into the natural nucleobases. Although the addition of the spacer does not significantly affect the hydrogen-bonding properties of natural nucleobases, the orientation of the base about the glycosidic bond necessary for Watson-Crick binding is destabilized, which could have implications for the selectivity of expanded bases, as well as the stability of expanded duplexes. Consideration of the (stacked) binding energies in the preferred relative orientation of natural and expanded nucleobases aligned according to their centers of mass reveals that the stacking within natural dimers can be increased by up to 50% upon expansion of one nucleobase and up to 90% upon expansion of two nucleobases. The implications of these findings to the stability of expanded duplexes were revealed by considering simplified models of natural and mixed duplexes composed of four nucleobases. Although intra- and interstrand interactions within double helices are typically less than those predicted when nucleobases are stacked according to their centers of mass, some nucleobases utilize their full stacking potential within double helices, where both intra- and interstrand interactions can be significant. Most importantly, increasing the size of nucleobases within the duplex significantly increases both intra- and interstrand stacking interactions. Specifically, some interactions are double the magnitude of the corresponding intrastrand interactions in natural helices, and even greater increases in interstrand interactions are sometimes found. Thus, our work suggests that mixed duplexes composed of natural bases hydrogen bound to expanded bases may exploit the increase in the inherent stacking ability of the expanded bases in more than one way and thereby afford duplexes with greater stability than natural DNA.
Collapse
Affiliation(s)
- Tom L McConnell
- Department of Chemistry, Mount Allison University, 63C York Street, Sackville, New Brunswick, E4L 1G8, Canada
| | | |
Collapse
|
26
|
Hunter KC, Millen AL, Wetmore SD. Effects of Hydrogen-Bonding and Stacking Interactions with Amino Acids on the Acidity of Uracil. J Phys Chem B 2007; 111:1858-71. [PMID: 17256895 DOI: 10.1021/jp066902p] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The effects of hydrogen-bonding interactions with amino acids on the (N1) acidity of uracil are evaluated using (B3LYP) density functional theory. Many different binding arrangements of each amino acid to three uracil binding sites are considered. The effects on the uracil acidity are found to significantly depend upon the nature of the amino acid and the binding orientation, but weakly depend on the binding site. Our results reveal that in some instances small models for the amino acids can be used, while for other amino acids larger models are required to properly describe the binding to uracil. The gas-phase acidity of uracil is found to increase by up to approximately 60 kJ mol(-1) due to discrete hydrogen-bonding interactions. Although (MP2) stacking interactions with aromatic amino acids decrease the acidity of uracil, unexpected increases in the acidity are found when any of the aromatic amino acids, or the backbone, hydrogen bond to uracil. Consideration of enzymatic and aqueous environments leads to decreases in the effects of the amino acids on the acidity of uracil. However, we find that the magnitude of the decrease varies with the nature of the molecule bound, as well as the (gas-phase) binding orientations and strengths, and therefore solvation effects should be considered on a case-by-case basis in future work. Nevertheless, the effects of amino acid interactions within enzymatic environments are as much as approximately 35 kJ mol(-1). The present study has general implications for understanding the nature of active site amino acids in enzymes, such as DNA repair enzymes, that catalyze reactions involving anionic nucleobase intermediates.
Collapse
Affiliation(s)
- Ken C Hunter
- Department of Chemistry, Mount Allison University, 63C York Street, Sackville, New Brunswick E4L 1G8, Canada
| | | | | |
Collapse
|
27
|
Gutlé C, Salpin JY, Cartailler T, Tortajada J, Gaigeot MP. Proton Transfers Induced by Lead(II) in a Uracil Nucleobase: A Study Based on Quantum Chemistry Calculations. J Phys Chem A 2006; 110:11684-94. [PMID: 17034162 DOI: 10.1021/jp0621528] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Within the context of metal biotoxicity, electrospray ionization mass spectrometry experiments (ESIMS) have recently been performed by us on the pyrimidine nucleobases (B) uracil and thymine complexed with lead(II) [Int. J. Mass. Spectrom. 2005, 243, 279]. Among the ions detected, [Pb(B)-H]+ complexes, where the base has been deprotonated, have been identified as producing intense signals. In the same study, quantum calculations based on density functional theory (DFT) have assessed the complexation sites and energies of [Pb(B)-H]+ ions. The present DFT investigations aim at giving an understanding on the energetics and mechanisms associated with uracil's loss of a proton. We specifically assess and quantify the role of lead binding in this process. For that purpose, intra- and intermolecular proton transfers have been considered. We have found that uracil (U) 1,3-tautomerization can be exergonic when uracil is complexed with Pb2+, in opposition to the situation without lead. The corresponding intramolecular processes were nonetheless found to occur at geological time scales. In contrast, the addition of a second body to [Pb(U)]2+ complexes, namely OH- or H2O (as found in the initial water droplet of ESIMS experiments), gives exergonic and fast uracil 1,3-proton transfers. Finally, we have shown that intermolecular proton transfers in uracil-H2O, uracil-OH-, or uracil-uracil complexes are able to explain the experimentally detected [Pb(U)-H]+ ions.
Collapse
Affiliation(s)
- C Gutlé
- Université d'Evry val d'Essonne, Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, LAMBE UMR-8587, Rue Père A. Jarland, F-91025 Evry, France
| | | | | | | | | |
Collapse
|
28
|
Nitrosubstituted aromatic molecules as universal nucleobases: Computational analysis of stacking interactions. Chem Phys Lett 2006. [DOI: 10.1016/j.cplett.2006.07.051] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
29
|
Galindo JF, Bermúdez CI, Daza EE. tRNA structure from a graph and quantum theoretical perspective. J Theor Biol 2006; 240:574-82. [PMID: 16337238 DOI: 10.1016/j.jtbi.2005.10.017] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2005] [Revised: 09/23/2005] [Accepted: 10/25/2005] [Indexed: 10/25/2022]
Abstract
One of the objectives of theoretical biochemistry is to find a suitable representation of molecules allowing us to encode what we know about their structures, interactions and reactivity. Particularly, tRNA structure is involved in some processes like aminoacylation and genetic code translation, and for this reason these molecules represent a biochemical object of the utmost importance requiring characterization. We propose here two fundamental aspects for characterizing and modeling them. The first takes into consideration the connectivity patterns, i.e. the set of linkages between atoms or molecular fragments (a key tool for this purpose is the use of graph theory), and the second one requires the knowledge of some properties related to the interactions taking place within the molecule, at least in an approximate way, and perhaps of its reactivity in certain means. We used quantum mechanics to achieve this goal; specifically, we have used partial charges as a manifestation of the reply to structural changes. These charges were appropriately modified to be used as weighted factors for elements constituting the molecular graph. This new graph-tRNA context allow us to detect some structure-function relationships.
Collapse
Affiliation(s)
- Johan F Galindo
- Grupo de Química Teórica-Universidad Nacional de Colombia, Bogotá, D.C., Colombia
| | | | | |
Collapse
|
30
|
Waller MP, Robertazzi A, Platts JA, Hibbs DE, Williams PA. Hybrid density functional theory for π-stacking interactions: Application to benzenes, pyridines, and DNA bases. J Comput Chem 2006; 27:491-504. [PMID: 16444702 DOI: 10.1002/jcc.20363] [Citation(s) in RCA: 226] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The suitability of a hybrid density functional to qualitatively reproduce geometric and energetic details of parallel pi-stacked aromatic complexes is presented. The hybrid functional includes an ad hoc mixture of half the exact (HF) exchange with half of the uniform electron gas exchange, plus Lee, Yang, and Parr's expression for correlation energy. This functional, in combination with polarized, diffuse basis sets, gives a binding energy for the parallel-displaced benzene dimer in good agreement with the best available high-level calculations reported in the literature, and qualitatively reproduces the local MP2 potential energy surface of the parallel-displaced benzene dimer. This method was further critically compared to high-level calculations recently reported in the literature for a range of pi-stacked complexes, including monosubstituted benzene-benzene dimers, along with DNA and RNA bases, and generally agrees with MP2 and/or CCSD(T) results to within +/-2 kJ mol(-1). We also show that the resulting BH&H binding energy is closely related to the electron density in the intermolecular region. The net result is that the BH&H functional, presumably due to fortuitous cancellation of errors, provides a pragmatic, computationally efficient quantum mechanical tool for the study of large pi-stacked systems such as DNA.
Collapse
Affiliation(s)
- Mark P Waller
- Faculty of Pharmacy, University of Sydney, NSW 2006, Australia
| | | | | | | | | |
Collapse
|
31
|
The MP2 quantum chemistry study on the local minima of guanine stacked with all four nucleic acid bases in conformations corresponding to mean B-DNA. ACTA ACUST UNITED AC 2005. [DOI: 10.1016/j.theochem.2005.06.014] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
32
|
Bera PP, Schaefer HF. (G-H)*-C and G-(C-H)* radicals derived from the guanine.cytosine base pair cause DNA subunit lesions. Proc Natl Acad Sci U S A 2005; 102:6698-703. [PMID: 15814617 PMCID: PMC1100756 DOI: 10.1073/pnas.0408644102] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The radicals generated by the homolytic cleavage of an X-H bond from the guanine.cytosine (G.C) base pair were studied by using carefully calibrated theoretical methods. The gradient-corrected density functional B3LYP was applied in conjunction with double-zeta plus polarization and diffuse function basis sets. Optimized geometries, energies, and vibrational frequencies were obtained for all of the radicals considered. Structural perturbations along with energy relaxation due to radical formation were investigated. Dissociation energies of the G.C base pair and all of the radicals are predicted and compared with the dissociation energy of neutral G.C. The three lowest-energy base pair radicals all involve removal of an H atom from one of the N atoms in G.C. The lowest-energy base pair radical has the hydrogen atom removed from the guanine nitrogen atom used for the sugar phosphate linkage in DNA. This (G-H)(*)-C radical has a dissociation energy (to G-H(*) + C) of 30 kcal/mol, which may be compared with 27 kcal/mol for G.C. All of the radicals that are possible outcomes of direct ionizing radiation or oxidizing species were investigated for the presence of local minima with significant structural changes. Major structural deformations cause strain in the interstrand hydrogen bonding in the DNA double helix. Severe geometry changes were observed when the hydrogen was abstracted from interstrand hydrogen bonding sites, along with sizeable energy changes, indicating the potentially serious consequences to the G.C base pair.
Collapse
Affiliation(s)
- Partha Pratim Bera
- Center for Computational Chemistry, University of Georgia, Athens, GA 30602, USA
| | | |
Collapse
|
33
|
Shishkin OV, Gorb L, Zhikol OA, Leszczynski J. Conformational analysis of canonical 2-deoxyribonucleotides. 2. Purine nucleotides. J Biomol Struct Dyn 2005; 22:227-44. [PMID: 15317483 DOI: 10.1080/07391102.2004.10506998] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
The molecular structure of different conformers of isolated canonical purine 2'-deoxyribonucleotides 2-deoxyadenosine-5'-phosphate (pdA) and 2'-deoxyguanosine-5'-phosphate (pdG) was optimized using the B3LYP/6-31G(d) method. The results of the calculations reveal that the geometrical parameters and relative stability of the conformers significantly depend on the nature of the nucleobase, its orientation, the conformation of the furanose ring, the charge of the phosphate group and the character of the intramolecular hydrogen bonds. Analysis of the electron density distribution in purine nucleotides reveals the existence of a number of intramolecular hydrogen bonds. In general, the south conformer has a lower energy at the anti orientation of the base, and both conformers occur as the most stable for the syn orientation of the nucleobases. The results of the calculations reveal that the geometry and relative energy of the conformers of purine DNTs may be easily tuned by the charge of the phosphate group.
Collapse
Affiliation(s)
- Oleg V Shishkin
- Institute for Scintillation Materials, National Academy of Science of Ukraine, 60 Lenina Ave., Kharkiv 61001, Ukraine.
| | | | | | | |
Collapse
|
34
|
Walsh TR. Exact exchange and Wilson–Levy correlation: a pragmatic device for studying complex weakly-bonded systems. Phys Chem Chem Phys 2005; 7:443-51. [DOI: 10.1039/b415563g] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
35
|
Mignon P, Loverix S, De Proft F, Geerlings P. Influence of Stacking on Hydrogen Bonding: Quantum Chemical Study on Pyridine−Benzene Model Complexes. J Phys Chem A 2004. [DOI: 10.1021/jp049240h] [Citation(s) in RCA: 118] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pierre Mignon
- Eenheid Algemene Chemie (ALGC) and Eenheid van Moleculaire en Cellulaire Interacties, Vlaams Interuniversitair Instituut Biotechnologie (VIB), Faculteit Wetenschappen, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
| | - Stefan Loverix
- Eenheid Algemene Chemie (ALGC) and Eenheid van Moleculaire en Cellulaire Interacties, Vlaams Interuniversitair Instituut Biotechnologie (VIB), Faculteit Wetenschappen, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
| | - Frank De Proft
- Eenheid Algemene Chemie (ALGC) and Eenheid van Moleculaire en Cellulaire Interacties, Vlaams Interuniversitair Instituut Biotechnologie (VIB), Faculteit Wetenschappen, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
| | - Paul Geerlings
- Eenheid Algemene Chemie (ALGC) and Eenheid van Moleculaire en Cellulaire Interacties, Vlaams Interuniversitair Instituut Biotechnologie (VIB), Faculteit Wetenschappen, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
| |
Collapse
|
36
|
Jurečka P, Šponer J, Hobza P. Potential Energy Surface of the Cytosine Dimer: MP2 Complete Basis Set Limit Interaction Energies, CCSD(T) Correction Term, and Comparison with the AMBER Force Field. J Phys Chem B 2004. [DOI: 10.1021/jp049956c] [Citation(s) in RCA: 82] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
37
|
Müller A, Losada M, Leutwyler S. Ab Initio Benchmark Study of (2-Pyridone)2, a Strongly Bound Doubly Hydrogen-Bonded Dimer. J Phys Chem A 2003. [DOI: 10.1021/jp0361024] [Citation(s) in RCA: 82] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Andreas Müller
- Departement für Chemie und Biochemie, Universität Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Martin Losada
- Departement für Chemie und Biochemie, Universität Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Samuel Leutwyler
- Departement für Chemie und Biochemie, Universität Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| |
Collapse
|
38
|
Structure of Stacked Dimers of N-Methylated Watson–Crick Adenine–Thymine Base Pairs. Int J Mol Sci 2003. [DOI: 10.3390/i4100537] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
|
39
|
|
40
|
Šponer J, Hobza P. Molecular Interactions of Nucleic Acid Bases. A Review of Quantum-Chemical Studies. ACTA ACUST UNITED AC 2003. [DOI: 10.1135/cccc20032231] [Citation(s) in RCA: 72] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Ab initio quantum-chemical calculations with inclusion of electron correlation significantly contributed to our understanding of molecular interactions of DNA and RNA bases. Some of the most important findings are introduced in the present overview: structures and energies of hydrogen bonded base pairs, nature of base stacking, interactions between metal cations and nucleobases, nonplanarity of isolated nucleobases and other monomer properties, tautomeric equilibria of nucleobases, out-of-plane hydrogen bonds and amino acceptor interactions. The role of selected molecular interactions in nucleic acids is discussed and representative examples where these interactions occur are given. Also, accuracy of density functional theory, semiempirical methods, distributed multipole analysis and empirical potentials is commented on. Special attention is given to our very recent reference calculations on base stacking and H-bonding. Finally, we briefly comment on the relationship between advanced ab initio quantum-chemical methods and large-scale explicit solvent molecular dynamics simulations of nucleic acids.
Collapse
|
41
|
Hobza P, Sponer J. Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations. J Am Chem Soc 2002; 124:11802-8. [PMID: 12296748 DOI: 10.1021/ja026759n] [Citation(s) in RCA: 328] [Impact Index Per Article: 14.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Stacking energies in low-energy geometries of pyrimidine, uracil, cytosine, and guanine homodimers were determined by the MP2 and CCSD(T) calculations utilizing a wide range of split-valence, correlation-consistent, and bond-functions basis sets. Complete basis set MP2 (CBS MP2) stacking energies extrapolated using aug-cc-pVXZ (X = D, T, and for pyrimidine dimer Q) basis sets equal to -5.3, -12.3, and -11.2 kcal/mol for the first three dimers, respectively. Higher-order correlation corrections estimated as the difference between MP2 and CCSD(T) stacking energies amount to 2.0, 0.7, and 0.9 kcal/mol and lead to final estimates of the genuine stacking energies for the three dimers of -3.4, -11.6, and -10.4 kcal/mol. The CBS MP2 stacking-energy estimate for guanine dimer (-14.8 kcal/mol) was based on the 6-31G(0.25) and aug-cc-pVDZ calculations. This simplified extrapolation can be routinely used with a meaningful accuracy around 1 kcal/mol for large aromatic stacking clusters. The final estimate of the guanine stacking energy after the CCSD(T) correction amounts to -12.9 kcal/mol. The MP2/6-31G(0.25) method previously used as the standard level to calculate aromatic stacking in hundreds of geometries of nucleobase dimers systematically underestimates the base stacking by ca. 1.0-2.5 kcal/mol per stacked dimer, covering 75-90% of the intermolecular correlation stabilization. We suggest that this correction is to be considered in calibration of force fields and other cheaper computational methods. The quality of the MP2/6-31G(0.25) predictions is nevertheless considerably better than suggested on the basis of monomer polarizability calculations. Fast and very accurate estimates of the MP2 aromatic stacking energies can be achieved using the RI-MP2 method. The CBS MP2 calculations and the CCSD(T) correction, when taken together, bring only marginal changes to the relative stability of H-bonded and stacked base pairs, with a slight shift of ca. 1 kcal/mol in favor of H-bonding. We suggest that the present values are very close to ultimate predictions of the strength of aromatic base stacking of DNA and RNA bases.
Collapse
Affiliation(s)
- Pavel Hobza
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Center for Complex Molecular Systems and Biomolecules, Dolejskova 3, 182 23 Prague 8, Czech Republic.
| | | |
Collapse
|
42
|
|
43
|
Li Q, Hu HY, Lam CK, Mak TC. Novel layer- and channel-type anionic host lattices built of thiourea, 2-thiouracil anion and water. ACTA ACUST UNITED AC 2002. [DOI: 10.1016/s1472-7862(03)00069-8] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
44
|
Leininger ML, Nielsen IMB, Colvin ME, Janssen CL. Accurate Structures and Binding Energies for Stacked Uracil Dimers. J Phys Chem A 2002. [DOI: 10.1021/jp013866f] [Citation(s) in RCA: 69] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Matthew L. Leininger
- Sandia National Laboratories, P.O. Box 969, Livermore, California 94551, and Computational Biology Group, Lawrence Livermore National Laboratory L-452, Livermore, California 94550
| | - Ida M. B. Nielsen
- Sandia National Laboratories, P.O. Box 969, Livermore, California 94551, and Computational Biology Group, Lawrence Livermore National Laboratory L-452, Livermore, California 94550
| | - Michael E. Colvin
- Sandia National Laboratories, P.O. Box 969, Livermore, California 94551, and Computational Biology Group, Lawrence Livermore National Laboratory L-452, Livermore, California 94550
| | - Curtis L. Janssen
- Sandia National Laboratories, P.O. Box 969, Livermore, California 94551, and Computational Biology Group, Lawrence Livermore National Laboratory L-452, Livermore, California 94550
| |
Collapse
|
45
|
Ab initio computational approaches to weakly interacting systems in the framework of the valence bond theory: from small to large van der Waals molecules. ACTA ACUST UNITED AC 2002. [DOI: 10.1016/s1380-7323(02)80012-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
|
46
|
|
47
|
Rauhut G. Recent Advances in Computing Heteroatom-Rich Five- and Six-Membered Ring Systems. ADVANCES IN HETEROCYCLIC CHEMISTRY 2001. [DOI: 10.1016/s0065-2725(01)81010-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/04/2023]
|
48
|
Shukla M, Kumar A, Mishra P. An ab initio study of excited state molecular electrostatic potential maps and other related properties of 5-fluorouracil: a comparative study with uracil. ACTA ACUST UNITED AC 2001. [DOI: 10.1016/s0166-1280(00)00584-4] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
49
|
Orozco M, Luque FJ. Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. Chem Rev 2000; 100:4187-4226. [PMID: 11749344 DOI: 10.1021/cr990052a] [Citation(s) in RCA: 454] [Impact Index Per Article: 18.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Modesto Orozco
- Departament de Bioquímica i Biologia Molecular, Facultat de Química, Universitat de Barcelona, Martí i Franqués 1, E-08028 Barcelona, Spain, and Departament de Fisicoquímica, Facultat de Farmàcia, Universitat de Barcelona, Avgda. Diagonal s/n, E-08028 Barcelona, Spain
| | | |
Collapse
|
50
|
Shishkin OV, Gorb L, Leszczynski J. Conformational flexibility of pyrimidine ring in adenine and related compounds. Chem Phys Lett 2000. [DOI: 10.1016/s0009-2614(00)01127-1] [Citation(s) in RCA: 41] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|