1
|
Poddar S, Levitus M. Buffer-Dependent Photophysics of 2-Aminopurine: Insights into Fluorescence Quenching and Excited-State Interactions. J Phys Chem B 2024; 128:2640-2651. [PMID: 38452253 DOI: 10.1021/acs.jpcb.3c07269] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/09/2024]
Abstract
2-Aminopurine (2AP) is the most widely used fluorescent nucleobase analogue in DNA and RNA research. Its unique photophysical properties and sensitivity to environmental changes make it a useful tool for understanding nucleic acid dynamics and DNA-protein interactions. We studied the effect of ions present in commonly used buffer solutions on the excited-state photophysical properties of 2AP. Fluorescence quenching was negligible for tris(hydroxymethyl)aminomethane (TRIS), but significant for phosphate, carbonate, 3-(N-morpholino) propanesulfonic acid (MOPS), and 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES) buffers. Results indicate that the two tautomers of 2AP (7H, 9H) are quenched by phosphate ions to different extents. Quenching by the H2PO4- ion is more pronounced for the 7H tautomer, while the opposite is true for the HPO42- ion. For phosphate ions, the results of the time-resolved fluorescence study cannot be explained using a simple collisional quenching mechanism. Instead, results are consistent with transient interactions between 2AP and the phosphate ions. We postulate that excited-state interactions between the 2AP tautomers and an H-bond acceptor (phosphate and carbonate) result in significant quenching of the singlet-excited state of 2AP. Such interactions manifest in biexponential fluorescence intensity decays with pre-exponential factors that vary with quencher concentration, and downward curvatures of the Stern-Volmer plots.
Collapse
Affiliation(s)
- Souvik Poddar
- School of Molecular Sciences, Arizona State University, 551 E. University Drive, Tempe, Arizona 85287, United States
- The Biodesign Institute Center for Single Molecule Biophysics, Arizona State University, Tempe, Arizona 85287, United States
| | - Marcia Levitus
- School of Molecular Sciences, Arizona State University, 551 E. University Drive, Tempe, Arizona 85287, United States
- The Biodesign Institute Center for Single Molecule Biophysics, Arizona State University, Tempe, Arizona 85287, United States
| |
Collapse
|
2
|
Brovarets’ OO, Muradova A, Hovorun DM. Novel horizons of the conformationally-tautomeric transformations of the G·T base pairs: quantum-mechanical investigation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2026510] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Ol’ha O. Brovarets’
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
| | - Alona Muradova
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| |
Collapse
|
3
|
Fan Q, Zhu L, Ren H, Lin H, Wu G. A new metal-organic framework of 3,9-diazatetraasterane-1,5,7,11-tetracarboxylic acid-3,6,9,12-tetraphenyl with sodium ion: Synthesis, characterization and DFT calculations. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138469] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
4
|
Abstract
Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.
Collapse
|
5
|
Brovarets’ OO, Hovorun DM. Key microstructural mechanisms of the 2-aminopurine mutagenicity: Results of extensive quantum-chemical research. J Biomol Struct Dyn 2019; 37:2716-2732. [DOI: 10.1080/07391102.2018.1495577] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Ol’ha O. Brovarets’
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, 2-h Akademika Hlushkova Ave, Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, 2-h Akademika Hlushkova Ave, Kyiv, Ukraine
| |
Collapse
|
6
|
Brovarets’ OO, Tsiupa KS, Hovorun DM. Novel pathway for mutagenic tautomerization of classical А∙Т DNA base pairs via sequential proton transfer through quasi-orthogonal transition states: A QM/QTAIM investigation. PLoS One 2018; 13:e0199044. [PMID: 29949602 PMCID: PMC6021055 DOI: 10.1371/journal.pone.0199044] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2018] [Accepted: 05/30/2018] [Indexed: 12/13/2022] Open
Abstract
In this paper we have theoretically predicted a novel pathway for the mutagenic tautomerization of the classical A∙T DNA base pairs in the free state, the Watson-Crick A·Т(WC), reverse Watson-Crick A·Т(rWC), Hoogsteen A·Т(H) and reverse Hoogsteen A·Т(rH) pairs, via sequential proton transfer accompanied by a significant change in the mutual orientation of the bases. Quantum-mechanical (QM) calculations were performed at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level in vacuum phase, along with Bader's quantum theory of Atoms in Molecules (QTAIM). These processes involve transition states (TSs) with quasi-orthogonal structures (symmetry C1), which are highly polar, tight ion pairs (A-, N6H2-deprotonated)∙(T+, O4/O2-protonated). Gibbs free energies of activation for the A∙T(WC) / A∙T(rWC) ↔ A*∙Т(rwWC) / A*∙Т(wWC) tautomeric transitions (~43.5 kcal∙mol-1) are lower than for the A∙T(H) / A∙T(rH) ↔ A*N7∙Т(rwH) / A*N7∙Т(wH) tautomerisations (~53.0 kcal∙mol-1) (rare tautomers are marked by an asterisk; w-wobble configured tautomerisation products). The (T)N3+H⋯N1-(A), (T)O4+H⋯N1-(A) / (T)N3+H⋯N1-(A) and (T)O2+H⋯N1-(A) H-bonds are found in the transition states TSA-·T+A·T(WC)↔A*·T(rwWC) / TSA-·T+A·T(rWC)↔A*·T(wWC). However, in the transition state TSA-·T+A·Т(H)↔A*N7·T(rwH) / TSA-·T+A·Т(rH)↔A*N7·T(wH), the (T)N3+H⋯N7-(A), (T)O4+H⋯N7-(A) / (T)N3+H⋯N7-(A) and (T)O2+H⋯N7-(A) H-bonds are supplemented by the attractive (T)O4+/O2+⋯N6-(A) van der Waals contacts. It was demonstrated that the products of the tautomerization of the classical A∙T DNA base pairs-A*∙Т(rwWC), A*N7∙Т(rwH) and A*N7∙Т(wH) (symmetry Cs)-further transform via double proton transfer into the energetically favorable wobble A∙T*(rwWC), A∙T*(rwH) and A∙T*O2(wH) base mispairs (symmetry Cs).
Collapse
Affiliation(s)
- Ol’ha O. Brovarets’
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| | - Kostiantyn S. Tsiupa
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| |
Collapse
|
7
|
Brovarets’ OO, Hovorun DM. Atomistic mechanisms of the double proton transfer in the H-bonded nucleobase pairs: QM/QTAIM computational lessons. J Biomol Struct Dyn 2018; 37:1880-1907. [DOI: 10.1080/07391102.2018.1467795] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- Ol’ha O. Brovarets’
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine , Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv , Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine , Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv , Kyiv, Ukraine
| |
Collapse
|
8
|
Brovarets' OO, Voiteshenko IS, Hovorun DM. Physico-chemical profiles of the wobble ↔ Watson-Crick G*·2AP(w) ↔ G·2AP(WC) and A·2AP(w) ↔ A*·2AP(WC) tautomerisations: a QM/QTAIM comprehensive survey. Phys Chem Chem Phys 2018; 20:623-636. [PMID: 29227488 DOI: 10.1039/c7cp05139e] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
This study is intended to clarify in detail the tautomeric transformations of the wobble (w) G*·2AP(w) and A·2AP(w) nucleobase mispairs involving 2-aminopurine (2AP) into the Watson-Crick (WC) G·2AP(WC) and A*·2AP(WC) base mispairs (asterisks denote mutagenic tautomers of the DNA bases), respectively, by quantum-mechanical methods and Bader's Quantum Theory of Atoms in Molecules. Our previously reported methodology has been used, which allows the evolution of the physico-chemical parameters to be tracked along the entire internal reaction coordinate (IRC), not exclusively in the stationary states of these reactions. These biologically important G*·2AP(w) ↔ G·2AP(WC) and A·2AP(w) ↔ A*·2AP(WC) w ↔ WC tautomerisations, which are involved in mutagenic tautomerically-conformational pathways, determine the origin of the transitions and transversions induced by 2AP. In addition, it is established that they proceed through planar, highly stable, zwitterionic transition states and they exhibit similar physico-chemical profiles and stages of sequential intrapair proton transfer, followed by spatial rearrangement of the nucleobases relative to each other within the base pairs. These w ↔ WC tautomerisations occur non-dissociatively and are accompanied by a significant alteration in geometry (from wobble to Watson-Crick and vice versa) and redistribution of the specific intermolecular interactions, which can be divided into 10 patterns including AHB H-bonds and loosened A-H-B covalent bridges along the IRC of tautomerisation. Based on the redistribution of the geometrical and electron-topological parameters of the intrapair hydrogen bonds, exactly 9 key points have been allocated to characterize the evolution of these reactions.
Collapse
Affiliation(s)
- Ol'ha O Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, 150 Akademika Zabolotnoho Str., 03680 Kyiv, Ukraine.
| | | | | |
Collapse
|
9
|
Brovarets' OO, Tsiupa KS, Hovorun DM. The A·T(rWC)/A·T(H)/A·T(rH) ↔ A·T*(rwWC)/A·T*(wH)/A·T*(rwH) mutagenic tautomerization via sequential proton transfer: a QM/QTAIM study. RSC Adv 2018; 8:13433-13445. [PMID: 35542561 PMCID: PMC9079753 DOI: 10.1039/c8ra01446a] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2018] [Accepted: 03/23/2018] [Indexed: 12/14/2022] Open
Abstract
In this study for the first time we have revealed by QM and QTAIM calculations at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of QM theory the novel routes of the mutagenic tautomerization of three biologically important A·T DNA base pairs – reverse Watson–Crick A·T(rWC), Hoogsteen A·T(H) and reverse Hoogsteen A·T(rH) – followed by their rebuilding into the wobble (w) A·T*(rwWC), A·T*(wH) and A·T*(rwH) base mispairs by the participation of the mutagenic tautomers of the DNA bases (denoted by asterisk) and vice versa, thus complementing the physico-chemical property of the canonical A·T(WC) Watson–Crick DNA base pair reported earlier (Brovarets' et al., RSC Adv., 2015, 5, 99594–99605). These non-dissociative tautomeric transformations in the classical A·T(rWC), A·T(H) and A·T(rH) DNA base pairs proceed similarly to the canonical A·T(WC) DNA base pair via the intrapair sequential proton transfer with shifting towards major or minor grooves of DNA followed by further double proton transfer along the intermolecular H-bonds and are controlled by the plane symmetric and highly stable transition states – tight ion pairs formed by the A+ nucleobase, protonated by the N1/N7 nitrogen atoms, and T− nucleobase, deprotonated by the N3H imino group. Comparison of the estimated populations of the tautomerised states (10−21 to 10−14) with similar characteristics for the canonical A·T(WC) DNA base pair (10−8 to 10−7) leads authors to the conclusion, that only a base pair with WC architecture can be a building block of the DNA macromolecule as a genetic material, which is able for the evolutionary self-development. Among all four classical DNA base pairs, only A·T(WC) DNA base pair can ensure the proper rate of the spontaneous point errors of replication in DNA. We discovered tautomeric wobbling of the classical A·T DNA base pairs. This data evidence, that only a base pair with Watson–Crick architecture can be a building block of the DNA macromolecule as a genetic material, which is able for the evolutionary self-development.![]()
Collapse
Affiliation(s)
- Ol'ha O. Brovarets'
- Department of Molecular and Quantum Biophysics
- Institute of Molecular Biology and Genetics
- National Academy of Sciences of Ukraine
- 03680 Kyiv
- Ukraine
| | - Kostiantyn S. Tsiupa
- Department of Molecular and Quantum Biophysics
- Institute of Molecular Biology and Genetics
- National Academy of Sciences of Ukraine
- 03680 Kyiv
- Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics
- Institute of Molecular Biology and Genetics
- National Academy of Sciences of Ukraine
- 03680 Kyiv
- Ukraine
| |
Collapse
|
10
|
Radhika R, Shankar R, Vijayakumar S, Kolandaivel P. Role of 6-Mercaptopurine in the potential therapeutic targets DNA base pairs and G-quadruplex DNA: insights from quantum chemical and molecular dynamics simulations. J Biomol Struct Dyn 2017; 36:1369-1401. [PMID: 28436311 DOI: 10.1080/07391102.2017.1323013] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The theoretical studies on DNA with the anticancer drug 6-Mercaptopurine (6-MP) are investigated using theoretical methods to shed light on drug designing. Among the DNA base pairs considered, 6-MP is stacked with GC with the highest interaction energy of -46.19 kcal/mol. Structural parameters revealed that structure of the DNA base pairs is deviated from the planarity of the equilibrium position due to the formation of hydrogen bonds and stacking interactions with 6-MP. These deviations are verified through the systematic comparison between X-H bond contraction and elongation and the associated blue shift and red shift values by both NBO analysis and vibrational analysis. Bent's rule is verified for the C-H bond contraction in the 6-MP interacted base pairs. The AIM results disclose that the higher values of electron density (ρ) and Laplacian of electron density (∇2ρ) indicate the increased overlap between the orbitals that represent the strong interaction and positive values of the total electron density show the closed-shell interaction. The relative sensitivity of the chemical shift values for the DNA base pairs with 6-MP is investigated to confirm the hydrogen bond strength. Molecular dynamics simulation studies of G-quadruplex DNA d(TGGGGT)4 with 6-MP revealed that the incorporation of 6-MP appears to cause local distortions and destabilize the G-quadruplex DNA.
Collapse
Affiliation(s)
- R Radhika
- a Department of Physics , Bharathiar University , Coimbatore , India
| | - R Shankar
- a Department of Physics , Bharathiar University , Coimbatore , India
| | - S Vijayakumar
- b Department of Medical Physics , Bharathiar University , Coimbatore , India
| | - P Kolandaivel
- a Department of Physics , Bharathiar University , Coimbatore , India
| |
Collapse
|
11
|
Srivastava R. Complexes of DNA bases and Watson–Crick base pairs interaction with neutral silver Agn (n = 8, 10, 12) clusters: a DFT and TDDFT study. J Biomol Struct Dyn 2017; 36:1050-1062. [DOI: 10.1080/07391102.2017.1310059] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Ruby Srivastava
- Center for Molecular Modeling, CSIR-Indian Institute of Chemical Technology, Hyderabad 500607, India
| |
Collapse
|
12
|
Tolosa S, Hidalgo A, Sansón JA. Theoretical study of enzymatically catalyzed tautomerization of carbon acids in aqueous solution: quantum calculations and steered molecular dynamics simulations. J Mol Model 2016; 22:44. [PMID: 26815031 DOI: 10.1007/s00894-016-2914-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2015] [Accepted: 01/11/2016] [Indexed: 10/22/2022]
Abstract
The thermodynamics and kinetics of enzymatically assisted reactions of carbon acids were studied theoretically in this work. Quantum electronic (QE) structure calculations and steered molecular dynamics (SMD) simulations were carried out. Three 3-butenal tautomerization reactions that proceed from the β,γ-unsaturated reactant (R) to the α,β-unsaturated carbon acid product (P) and occur in two elementary steps through an intermediate (I) were studied, ignoring or including the surrounding aqueous medium in the calculations. The Gibbs free energies of activation of the R ⇆ I enolization and I ⇆ P ketonization steps were found to decrease considerably when residues simulating enzymes were introduced into these processes. Although the processes became slightly more favorable thermodynamically when the solution was included in the simulations, they became less favorable kinetically. The results from SMD simulations of these reactions were qualitatively consistent with the values we obtained using QE as well as those found by other authors in similar studies. Our simulations also allowed us to perform a detailed study of these reactions in solution.
Collapse
Affiliation(s)
- Santiago Tolosa
- Departamento de Ingeniería Química y Química Física, Universidad de Extremadura, Avda. Elvas s/n, 06071, Badajoz, Spain.
| | - Antonio Hidalgo
- Departamento de Ingeniería Química y Química Física, Universidad de Extremadura, Avda. Elvas s/n, 06071, Badajoz, Spain
| | - Jorge A Sansón
- Departamento de Ingeniería Química y Química Física, Universidad de Extremadura, Avda. Elvas s/n, 06071, Badajoz, Spain
| |
Collapse
|
13
|
Brovarets' OO, Pérez-Sánchez H. Whether the amino–imino tautomerism of 2-aminopurine is involved into its mutagenicity? Results of a thorough QM investigation. RSC Adv 2016. [DOI: 10.1039/c6ra24277d] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
2AP* mutagenic tautomer is able to induce only one incorporation error – transversion – by pairing through the H-bonds into the G·2AP* mispair.
Collapse
Affiliation(s)
- Ol'ha O. Brovarets'
- Department of Molecular and Quantum Biophysics
- Institute of Molecular Biology and Genetics
- National Academy of Sciences of Ukraine
- 03680 Kyiv
- Ukraine
| | - Horacio Pérez-Sánchez
- Computer Science Department
- Bioinformatics and High Performance Computing (BIO-HPC) Research Group
- Universidad Católica San Antonio de Murcia (UCAM)
- Murcia
- Spain
| |
Collapse
|
14
|
Izadi Z, Farrokhpour H, Tabrizchi M. Sodium affinity of caffeine and adenine: the effect of microsolvation and electrostatic field of solvent on the sodium affinity. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1115132] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
15
|
Zhang JK, Ma XF, Wang XW, Cai JQ, Wu XL, Niu YP, Zhou N, Hua RM. X-ray structure and conformational study of tris(2-(dimethylamino)phenyl)methanol salt derivatives: Roles of anions and hydrogen bond. J Mol Struct 2015. [DOI: 10.1016/j.molstruc.2015.07.039] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
16
|
Masoodi HR, Bagheri S, Abareghi M. The effects of tautomerization and protonation on the adenine-cytosine mismatches: a density functional theory study. J Biomol Struct Dyn 2015. [PMID: 26198186 DOI: 10.1080/07391102.2015.1072734] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
In the present work, we demonstrate the results of a theoretical study concerned with the question how tautomerization and protonation of adenine affect the various properties of adenine-cytosine mismatches. The calculations, in gas phase and in water, are performed at B3LYP/6-311++G(d,p) level. In gas phase, it is observed that any tautomeric form of investigated mismatches is more stabilized when adenine is protonated. As for the neutral mismatches, the mismatches containing amino form of cytosine and imino form of protonated adenine are more stable. The role of aromaticity on the stability of tautomeric forms of mismatches is investigated by NICS(1)ZZ index. The stability of mispairs decreases by going from gas phase to water. It can be explained using dipole moment parameter. The influence of hydrogen bonds on the stability of mismatches is examined by atoms in molecules and natural bond orbital analyses. In addition to geometrical parameters and binding energies, the study of the topological properties of electron charge density aids in better understanding of these mispairs.
Collapse
Affiliation(s)
- Hamid Reza Masoodi
- a Faculty of Science, Department of Chemistry , Vali-e-Asr University of Rafsanjan , P.O. Box 77176, Rafsanjan , Iran
| | - Sotoodeh Bagheri
- a Faculty of Science, Department of Chemistry , Vali-e-Asr University of Rafsanjan , P.O. Box 77176, Rafsanjan , Iran
| | - Mahsa Abareghi
- a Faculty of Science, Department of Chemistry , Vali-e-Asr University of Rafsanjan , P.O. Box 77176, Rafsanjan , Iran
| |
Collapse
|
17
|
Brovarets' OO, Hovorun DM. DPT tautomerisation of the G·A(syn) and A*·G*(syn) DNA mismatches: a QM/QTAIM combined atomistic investigation. Phys Chem Chem Phys 2015; 16:9074-85. [PMID: 24695821 DOI: 10.1039/c4cp00488d] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
By applying a combined QM and QTAIM atomistic computational approach we have established for the first time that the G·A(syn) and A*·G*(syn) DNA mismatches (rare tautomers are marked with an asterisk), causing spontaneous transversions with substantially various probabilities, radically differ from each other in their ability to tautomerise through the double proton transfer (DPT). The A*·G*(syn) mismatch tautomerises quite easily (ΔΔG(TS) ≈ 4·kT at room temperature) into the A·G*(syn) mismatch through the asynchronous concerted mechanism, whereas the G·A(syn) base mispair does not tautomerise via the DPT at all, since there is no local minimum corresponding to the tautomerised G*·A*(syn) mismatch on the potential energy surface. It was established that the A·G*(syn) base mispair is a dynamically unstable H-bonded complex with an extremely short lifetime of 2.17 × 10(-13) s. Consequently, the obtained results allow us to believe that spontaneous or forced dissociation of both the G·A(syn) and A*·G*(syn) DNA mismatches by the DNA-polymerase occurs with the preservation of the tautomeric status of the bases.
Collapse
Affiliation(s)
- Ol'ha O Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, 150 Akademika Zabolotnoho Str., 03680 Kyiv, Ukraine.
| | | |
Collapse
|
18
|
Siddiqui SA, Rasheed T, Bouarissa N, Al-Hajry A. Possible use of BN-modified fullerene as a nano-biosensor to detect adenine–thymine Watson–Crick base pair in mutagenic tautomeric form: Theoretical approach. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2015. [DOI: 10.1142/s0219633615500030] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/07/2022]
Abstract
The present work deals with the theoretical investigation of electronic structure features and stability of adenine–thymine (AT) and rare tautomer of adenine–thymine (rAT) base pairs along with their complexes with Cu 2+ cation and their interactions with BN doped fullerene ( C 58 BN ). All the calculations have been performed with density functional theory using B3LYP functional. Electronic structures of the two base pairs are almost identical. Hence, it is rather difficult to distinguish between the two base pairs on the basis of their electronic properties. As per our theoretical calculations, we have observed that, BN modified fullerene could act as a nano-biosensor for detection of mispairing between these two complementary bases as well as their Cu 2+ complexes.
Collapse
Affiliation(s)
- Shamoon Ahmad Siddiqui
- Promising Centre for Sensors and Electronic Devices, Najran University, Najran, KSA
- Department of Physics, College of Arts and Science, Najran University, Najran, KSA
| | - Tabish Rasheed
- Department of Applied Sciences, School of Engineering and Technology, Sharda University, Plot No. 32–34, Knowledge Park III, Greater Noida, N.C.R., U.P., India-201306, India
| | - Nadir Bouarissa
- Department of Physics, Faculty of Science, University of M'sila, 28000 M'sila, Algeria
| | - A. Al-Hajry
- Promising Centre for Sensors and Electronic Devices, Najran University, Najran, KSA
- Department of Physics, College of Arts and Science, Najran University, Najran, KSA
| |
Collapse
|
19
|
Brovarets’ OO, Yurenko YP, Hovorun DM. The significant role of the intermolecular CH⋯O/N hydrogen bonds in governing the biologically important pairs of the DNA and RNA modified bases: a comprehensive theoretical investigation. J Biomol Struct Dyn 2014; 33:1624-52. [DOI: 10.1080/07391102.2014.968623] [Citation(s) in RCA: 62] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
20
|
Brovarets' OO, Hovorun DM. Atomistic nature of the DPT tautomerisation of the biologically important C·C* DNA base mispair containing amino and imino tautomers of cytosine: a QM and QTAIM approach. Phys Chem Chem Phys 2014; 15:20091-104. [PMID: 24154739 DOI: 10.1039/c3cp52644e] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
A theoretical study of tautomerisation of the biologically important cytosine·cytosine* (C·C*) DNA mismatch with a propeller-like structure (|C4N3N3C4| = 32.4°; C1 symmetry) and cis-oriented N1H glycosidic bonds, formed by the amino and imino tautomers of the C nucleobase, via the asynchronous concerted double proton transfer (DPT) along two H-bonds through the transition state (TSC·C*↔C*·C) (|C4N3N3C4| = 48.5°; C1 symmetry) into the C*·C mispair was carried out for the first time. It was established that the C·C*/C*·C DNA base mispair is associated by the antiparallel N4H···N4 (6.66 kcal mol(-1)), N3H···N3 (6.47 kcal mol(-1)) H-bonds and the O2···O2 van der Waals (vdW) contact (0.33 kcal mol(-1)), while the zwitterionic TSC·C*↔C*·C is stabilized by the parallel N4(+)H···N4(-) (13.55 kcal mol(-1)), N3(+)H···N3(-) (13.20 kcal mol(-1)) H-bonds and the O2(+)···O2(-) vdW contact (0.60 kcal mol(-1)). It was shown that the C·C* ↔ C*·C tautomerisation via the DPT is assisted by the O2···O2 vdW contact, that in contrast to the two others N4H···N4 and N3H···N3 H-bonds exists along the entire intrinsic reaction coordinate (IRC) range. The positive values of the Grunenberg's compliance constants (30.919 and 21.384 Å mdyn(-1) for C·C*/C*·C and TSC·C*↔C*·C, respectively) indicate that the O2···O2 vdW contact is a stabilizing closed-shell interaction. It was found that the middle N3H···N3 H-bond is anti-cooperative with the upper N4H···N4 H-bond and cooperative with the lower O2···O2 vdW contact. The 9 key points, which can be considered as electron-topological "fingerprints" of the asynchronous concerted C·C* ↔ C*·C tautomerisation process via the DPT were revealed along the IRC and examined in detail. It was shown that the C·C*/C*·C base mispair is a thermodynamically and dynamically stable structure. Its lifetime is equal to 1.53 × 10(-7) s at the MP2/cc-pVQZ//B3LYP/6-311++G(d,p) level of theory in vacuum. All 6 low-frequency intermolecular vibrations are able to develop during this time span.
Collapse
Affiliation(s)
- Ol'ha O Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, 150 Akademika Zabolotnoho Str., 03680 Kyiv, Ukraine.
| | | |
Collapse
|
21
|
Brovarets' OO, Zhurakivsky RO, Hovorun DM. DPT tautomerisation of the wobble guanine·thymine DNA base mispair is not mutagenic: QM and QTAIM arguments. J Biomol Struct Dyn 2014; 33:674-89. [PMID: 24650179 DOI: 10.1080/07391102.2014.897259] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
We have shown for the first time, connecting QM methods with QTAIM analysis and using the methodology of the sweeps of the energetical, electron-topological and geometrical parameters, that the tautomerisation of the wobble guanine·thymine (wG·T) DNA base mispair into the wG(*)·T(*) base mispair induced by the double proton transfer (DPT), which undergoes a concerted asynchronous pathway, is not mutagenic. The wG·T → wG(*)·T(*) DPT tautomerisation does not result in the transition of the G base into its mutagenic tautomeric form G(*) able to mispair with the T base within the Watson-Crick base pairing scheme. This observation is explained by the so-called quantum protection of the wG·T DNA base mispair from its mutagenic tautomerisation - the dynamical non-stability of the tautomerised wG(*)·T(*) base mispair and significantly negative value of the Gibbs free energy of activation for the reverse reaction of the wG·T → wG(*)·T(*) DPT tautomerisation.
Collapse
Affiliation(s)
- Ol'ha O Brovarets'
- a Department of Molecular and Quantum Biophysics , Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine , 150 Akademika Zabolotnoho Str., 03680 Kyiv , Ukraine
| | | | | |
Collapse
|
22
|
Brovarets’ OO, Zhurakivsky RO, Hovorun DM. A QM/QTAIM microstructural analysis of the tautomerisationviathe DPT of the hypoxanthine·adenine nucleobase pair. Mol Phys 2014. [DOI: 10.1080/00268976.2013.877170] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
23
|
Brovarets' OO, Zhurakivsky RO, Hovorun DM. Does the tautomeric status of the adenine bases change upon the dissociation of the A*·Asyn Topal–Fresco DNA mismatch? A combined QM and QTAIM atomistic insight. Phys Chem Chem Phys 2014; 16:3715-25. [DOI: 10.1039/c3cp54708f] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
24
|
Brovarets’ OO, Zhurakivsky RO, Hovorun DM. Structural, energetic and tautomeric properties of the T·T∗/T∗·T DNA mismatch involving mutagenic tautomer of thymine: A QM and QTAIM insight. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2013.12.034] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
25
|
Brovarets' OO. [Structural and energetic properties of the four configurations of the A.T and G.C DNA base pairs]. UKRAINIAN BIOCHEMICAL JOURNAL 2013; 85:104-10. [PMID: 24319980 DOI: 10.15407/ubj85.04.104] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Using the methods of non-empirical quantum chemistry at the MP2/6-311++G(2df,pd)// B3LYP/6-311++G(d,p) level of theory it was established for the first time, that Hoogsteen, reverse Hoogsteen, Watson-Crick and reverse Watson-Crick configurations of the A.T and G.C DNA base pairs are isoelectronic and isomorphic structures with similar dynamic properties. Based on these results, non-ionisation mechanism of the Hoogsteen <"breathing" of the G*.C* DNA base pair, namely transformation of the tautomerised (Lowdin's) G-C base pair with Watson-Crick geometry into the Hoogsteen electroneutral G*.C* H base pair stabilized by the three O6H...N4, N3H...N7 and C8H...02 H-bonds, was postulated. It is suggested that such scenario activates only in those cases, when DNA is not located in aqueous solution, but works together with proteins and cytosine protonation at the N3 atom is precluded.
Collapse
|
26
|
Brovarets' OO. [Under what conditions does G.C Watson-Crick DNA base pair acquire all four configurations characteristic for A.T Watson-Crick DNA base pair?]. UKRAINIAN BIOCHEMICAL JOURNAL 2013; 85:98-103. [PMID: 24319979 DOI: 10.15407/ubj85.04.098] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
At the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory it was established for the first time, that the Löwdin's G*.C* DNA base pair formed by the mutagenic tautomers can acquire, as the A-T Watson-Crick DNA base pair, four biologically important configurations, namely: Watson-Crick, reverse Watson-Crick, Hoogsteen and reverse Hoogsteen. This fact demonstrates rather unexpected role of the tautomerisation of the one of the Watson-Crick DNA base pairs, in particular, via double proton transfer: exactly the G.C-->G*.C* tautomerisation allows to overcome steric hindrances for the implementation of the above mentioned configurations. Geometric, electron-topological and energetic properties of the H-bonds that stabilise the studied pairs, as well as the energetic characteristics of the latters are presented.
Collapse
|
27
|
Brovarets' OO, Zhurakivsky RO, Hovorun DM. Is the DPT tautomerization of the long A·G Watson-Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question. J Comput Chem 2013; 35:451-66. [PMID: 24382756 DOI: 10.1002/jcc.23515] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2013] [Revised: 11/12/2013] [Accepted: 11/30/2013] [Indexed: 02/04/2023]
Abstract
Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum-mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ↔ A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H-bonds) and the N2H···HC2 dihydrogen bond (DH-bond) (0.68 kcal·mol(-1) ), whereas the A*·G* base mispair-by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H-bonds (0.42 kcal·mol(-1) ). The N2H···HC2 DH-bond smoothly and without bifurcation transforms into the C2H···N2 H-bond at the IRC = -10.07 Bohr in the course of the A·G ↔ A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H-bonds, it was observed that the N6H···O6 H-bond is anticooperative to the two others-N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H-bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10(-14) s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low-frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form.
Collapse
Affiliation(s)
- Ol'ha O Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, 150 Akademika Zabolotnoho Str., 03680, Kyiv, Ukraine ; Research and Educational Center "State Key Laboratory of Molecular and Cell Biology", 150 Akademika Zabolotnoho Str., 03680, Kyiv, Ukraine; Department of Molecular Biology, Biotechnology and Biophysics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, 2-h Akademika Hlushkova Ave., 03022, Kyiv, Ukraine
| | | | | |
Collapse
|
28
|
Brovarets' OO, Hovorun DM. Atomistic understanding of the C·T mismatched DNA base pair tautomerization via the DPT: QM and QTAIM computational approaches. J Comput Chem 2013; 34:2577-90. [PMID: 23955922 DOI: 10.1002/jcc.23412] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2013] [Accepted: 07/25/2013] [Indexed: 11/09/2022]
Abstract
It was established that the cytosine·thymine (C·T) mismatched DNA base pair with cis-oriented N1H glycosidic bonds has propeller-like structure (|N3C4C4N3| = 38.4°), which is stabilized by three specific intermolecular interactions-two antiparallel N4H…O4 (5.19 kcal mol(-1)) and N3H…N3 (6.33 kcal mol(-1)) H-bonds and a van der Waals (vdW) contact O2…O2 (0.32 kcal mol(-1)). The C·T base mispair is thermodynamically stable structure (ΔG(int) = -1.54 kcal mol(-1) ) and even slightly more stable than the A·T Watson-Crick DNA base pair (ΔG(int) = -1.43 kcal mol(-1)) at the room temperature. It was shown that the C·T ↔ C*·T* tautomerization via the double proton transfer (DPT) is assisted by the O2…O2 vdW contact along the entire range of the intrinsic reaction coordinate (IRC). The positive value of the Grunenberg's compliance constants (31.186, 30.265, and 22.166 Å/mdyn for the C·T, C*·T*, and TS(C·T ↔ C*·T*), respectively) proves that the O2…O2 vdW contact is a stabilizing interaction. Based on the sweeps of the H-bond energies, it was found that the N4H…O4/O4H…N4, and N3H…N3 H-bonds in the C·T and C*·T* base pairs are anticooperative and weaken each other, whereas the middle N3H…N3 H-bond and the O2…O2 vdW contact are cooperative and mutually reinforce each other. It was found that the tautomerization of the C·T base mispair through the DPT is concerted and asynchronous reaction that proceeds via the TS(C·T ↔ C*·T*) stabilized by the loosened N4-H-O4 covalent bridge, N3H…N3 H-bond (9.67 kcal mol(-1) ) and O2…O2 vdW contact (0.41 kcal mol(-1) ). The nine key points, describing the evolution of the C·T ↔ C*·T* tautomerization via the DPT, were detected and completely investigated along the IRC. The C*·T* mispair was revealed to be the dynamically unstable structure with a lifetime 2.13·× 10(-13) s. In this case, as for the A·T Watson-Crick DNA base pair, activates the mechanism of the quantum protection of the C·T DNA base mispair from its spontaneous mutagenic tautomerization through the DPT.
Collapse
Affiliation(s)
- Ol'ha O Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, 150 Zabolotnoho Str., Kyiv, 03680, Ukraine; Research and Educational Center, State Key Laboratory of Molecular and Cell Biology, 150 Zabolotnoho Str., Kyiv, 03680, Ukraine; Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, 2 Hlushkova Ave., Kyiv, 03022, Ukraine
| | | |
Collapse
|
29
|
Brovarets’ OO, Hovorun DM. Why the tautomerization of the G·C Watson–Crick base pairviathe DPT does not cause point mutations during DNA replication? QM and QTAIM comprehensive analysis. J Biomol Struct Dyn 2013; 32:1474-99. [DOI: 10.1080/07391102.2013.822829] [Citation(s) in RCA: 71] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
30
|
Brovarets’ OO, Hovorun DM. Prototropic tautomerism and basic molecular principles of hypoxanthine mutagenicity: an exhaustive quantum-chemical analysis. J Biomol Struct Dyn 2013; 31:913-36. [PMID: 22962845 DOI: 10.1080/07391102.2012.715041] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
|
31
|
Brovarets’ OO, Zhurakivsky RO, Hovorun DM. The physico-chemical mechanism of the tautomerisation via the DPT of the long Hyp∗·Hyp Watson–Crick base pair containing rare tautomer: A QM and QTAIM detailed look. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2013.05.067] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
32
|
Brovarets’ OO, Yurenko YP, Hovorun DM. Intermolecular CH···O/N H-bonds in the biologically important pairs of natural nucleobases: a thorough quantum-chemical study. J Biomol Struct Dyn 2013; 32:993-1022. [PMID: 23730732 DOI: 10.1080/07391102.2013.799439] [Citation(s) in RCA: 109] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
33
|
Brovarets’ OO, Zhurakivsky RO, Hovorun DM. DPT tautomerization of the long A∙A* Watson-Crick base pair formed by the amino and imino tautomers of adenine: combined QM and QTAIM investigation. J Mol Model 2013; 19:4223-37. [DOI: 10.1007/s00894-013-1880-2] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2013] [Accepted: 05/02/2013] [Indexed: 12/01/2022]
|
34
|
Brovarets OO, Hovorun DM. Can tautomerization of the A·T Watson–Crick base pairviadouble proton transfer provoke point mutations during DNA replication? A comprehensive QM and QTAIM analysis. J Biomol Struct Dyn 2013; 32:127-54. [DOI: 10.1080/07391102.2012.755795] [Citation(s) in RCA: 112] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
35
|
Brovarets' OO, Zhurakivsky RO, Hovorun DM. The physico-chemical "anatomy" of the tautomerization through the DPT of the biologically important pairs of hypoxanthine with DNA bases: QM and QTAIM perspectives. J Mol Model 2013; 19:4119-37. [PMID: 23292249 DOI: 10.1007/s00894-012-1720-9] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2012] [Accepted: 12/03/2012] [Indexed: 10/27/2022]
Abstract
The biologically important tautomerization of the Hyp·Cyt, Hyp·Thy and Hyp·Hyp base pairs to the Hyp·Cyt, Hyp·Thy and Hyp·Hyp base pairs, respectively, by the double proton transfer (DPT) was comprehensively studied in vacuo and in the continuum with a low dielectric constant (ε = 4) corresponding to hydrophobic interfaces of protein-nucleic acid interactions by combining theoretical investigations at the B3LYP/6-311++G(d,p) level of QM theory with QTAIM topological analysis. Based on the sweeps of the energetic, electron-topological, geometric and polar parameters, which describe the course of the tautomerization along the intrinsic reaction coordinate (IRC), it was proved that the tautomerization through the DPT is concerted and asynchronous process for the Hyp·Cyt and Hyp·Thy base pairs, while concerted and synchronous for the Hyp·Hyp homodimer. The continuum with ε = 4 does not affect qualitatively the course of the tautomerization reaction for all studied complexes. The nine key points along the IRC of the Hyp·Cyt↔Hyp·Cyt and Hyp·Thy↔Hyp·Thy tautomerizations and the six key points of the Hyp·Hyp↔Hyp·Hyp tautomerization have been identified and fully characterized. These key points could be considered as electron-topological "fingerprints" of concerted asynchronous (for Hyp·Cyt and Hyp·Thy) or synchronous (for Hyp·Hyp) tautomerization process via the DPT. It was found, that in the Hyp·Cyt, Hyp·Thy, Hyp·Hyp and Hyp·Hyp base pairs all H-bonds are significantly cooperative and mutually reinforce each other, while the C2H…O2 H-bond in the Hyp·Cyt base pair and the O6H…O4 H-bond in the Hyp·Thy base pair behave anti-cooperatively, i.e., they become weakened, while two others become strengthened.
Collapse
Affiliation(s)
- Ol'ha O Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, 150 Zabolotnoho Str, 03680, Kyiv, Ukraine
| | | | | |
Collapse
|
36
|
Nikolaienko TY, Bulavin LA, Hovorun DM. Bridging QTAIM with vibrational spectroscopy: the energy of intramolecular hydrogen bonds in DNA-related biomolecules. Phys Chem Chem Phys 2012; 14:7441-7. [PMID: 22514024 DOI: 10.1039/c2cp40176b] [Citation(s) in RCA: 111] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Physical properties of over 8000 intramolecular hydrogen bonds (iHBs), including 2901 ones of the types OH···O, OH···N, NH···O and OH···C, in 4244 conformers of the DNA-related molecules (four canonical 2'-deoxyribonucleotides, 1,2-dideoxyribose-5-phosphate, and 2-deoxy-D-ribose in its furanose, pyranose and linear forms) have been investigated using quantum theory of atoms in molecules (QTAIM) and vibrational analysis. It has been found that for all iHBs with positive red-shift of the proton donating group stretching frequency the shift value correlates with ρ(cp)-the electron charge density at the (3,-1)-type bond critical point. Combining QTAIM and spectroscopic data new relationships for estimation of OH···O, OH···N, NH···O and OH···C iHB enthalpy of formation (kcal mol(-1)) with RMS error below 0.8 kcal mol(-1) have been established: E(OH···O) = -3.09 + 239·ρ(cp), E(OH···N) = 1.72 + 142·ρ(cp), E(NH···O) = -2.03 + 225·ρ(cp), E(OH···C) = -0.29 + 288·ρ(cp), where ρ(cp) is in e a(0)(-3) (a(0)- the Bohr radius). It has been shown that XHY iHBs with red-shift values over 40 cm(-1) are characterized by the following minimal values of the XHY angle, ρ(cp) and nubla(2)ρ(cp): 112°, 0.005 e a(0)(-3) and 0.016 e a(0)(-5), respectively. New relationships have been used to reveal the strongest iHBs in canonical 2'-deoxy- and ribonucleosides and the O(5')H···N(3) H-bond in ribonucleoside guanosine was found to have the maximum energy (8.1 kcal mol(-1)).
Collapse
Affiliation(s)
- Tymofii Yu Nikolaienko
- Molecular Physics Department, Faculty of Physics, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine.
| | | | | |
Collapse
|
37
|
Goh GB, Knight JL, Brooks CL. Constant pH Molecular Dynamics Simulations of Nucleic Acids in Explicit Solvent. J Chem Theory Comput 2011; 8:36-46. [PMID: 22337595 DOI: 10.1021/ct2006314] [Citation(s) in RCA: 74] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The nucleosides of adenine and cytosine have pKa values of 3.50 and 4.08, respectively, and are assumed to be unprotonated under physiological conditions. However, evidence from recent NMR and X-Ray crystallography studies has revealed the prevalence of protonated adenine and cytosine in RNA macromolecules. Such nucleotides with elevated pKa values may play a role in stabilizing RNA structure and participate in the mechanism of ribozyme catalysis. With the work presented here, we establish the framework and demonstrate the first constant pH MD simulations (CPHMD) for nucleic acids in explicit solvent in which the protonation state is coupled to the dynamical evolution of the RNA system via λ-dynamics. We adopt the new functional form λ(Nexp) for λ that was recently developed for Multi-Site λ-Dynamics (MSλD) and demonstrate good sampling characteristics in which rapid and frequent transitions between the protonated and unprotonated states at pH = pKa are achieved. Our calculated pKa values of simple nucleotides are in a good agreement with experimentally measured values, with a mean absolute error of 0.24 pKa units. This work demonstrates that CPHMD can be used as a powerful tool to investigate pH-dependent biological properties of RNA macromolecules.
Collapse
Affiliation(s)
- Garrett B Goh
- Department of Chemistry, University of Michigan, 930 N. University, Ann Arbor, Michigan 48109, United States
| | | | | |
Collapse
|
38
|
Lu P, Shi Z, Walker AV. Selective electroless deposition of copper on organic thin films with improved morphology. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2011; 27:13022-13028. [PMID: 21936558 DOI: 10.1021/la202839z] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
We have investigated the selective electroless deposition (ELD) of Cu on functionalized self-assembled monolayers (SAMs). Previous studies have demonstrated that Cu deposits on -COOH and -CH(3) terminated SAMs using ELD. However, the deposited films were rough and contained irregular crystallites. Further, the copper penetrated through the film. In this Article, we demonstrate that copper can be selectively deposited on -COOH terminated SAMs with improved morphology and without penetration of copper through the organic layer. The method employs a Cu(II) seed layer and an additive, adenine or guanine. We demonstrate the efficacy of the technique on photopatterned -CH(3)/-COOH SAMs. Copper is observed to deposit only atop the -COOH terminated SAM area and not on the -CH(3) terminated SAM. The use of a Cu(II) seed layer increased the Cu ELD rate on both -COOH and -CH(3) terminated SAMs. The deposited copper layer strongly adheres to the -COOH terminated SAMs because the copper layer nucleates at Cu(2+)-carboxylate complexes. In contrast, the deposited copper layer can easily be removed from the -CH(3) terminated SAM surface because there is no specific copper-surface interaction. The additives adenine and guanine mediate the interaction of Cu(2+) and the deprotonated -COOH terminated SAMs via the formation of additive-carboxylate complexes. These complexes lead to significantly reduced copper penetration through the SAM. In the case of adenine, the diffusion of copper through the organic film was eliminated. This new technique for copper deposition will facilitate the development of inexpensive molecular electronics, sensors, and other nanotechological devices.
Collapse
Affiliation(s)
- Peng Lu
- Department of Materials Science and Engineering, University of Texas at Dallas, 800 West Campbell Road, RL 10, Richardson, Texas 75080, United States
| | | | | |
Collapse
|
39
|
Samijlenko SP, Yurenko YP, Stepanyugin AV, Hovorun DM. Tautomeric equilibrium of uracil and thymine in model protein-nucleic acid contacts. Spectroscopic and quantum chemical approach. J Phys Chem B 2010; 114:1454-61. [PMID: 20047292 DOI: 10.1021/jp909099a] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
This work deals with tautomeric transformations of uracil (Ura) and thymine (Thy) in their model complexes with the deprotonated carboxylic group. Essential changes in the UV spectra of the bases upon their interaction with NaAc, vanishing signals of both imino protons in (1)H NMR spectra, and a perceptible decrease in intensity of both IR bands, related to the stretching vibrations nu(C=O) of the carbonyl groups, imply involvement of enolic tautomers. Results of quantum chemical calculations of the double complexes of the Ura(Thy) tautomers with CH(3)COO(-) at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory proved to be incompatible with the spectral features: despite the fact that the complexes of the enolic tautomers are much closer in energy to the diketo ones as compared to isolated tautomers, the energy gap between them is such that in tautomeric equilibrium dominate diketo forms. Calculations of triple complexes of the type CH(3)COO(-):Ura(Thy) tautomer:Na(+), taking into account the effect of the Na(+) coordination with tautomers, show that three triple complexes formed by enolic tautomers appeared more stable than those formed by diketo ones. This makes the UV and (1)H NMR data understandable, but the high residual intensity of the nu(C=O) bands in the IR spectra remains unclear. At that ion, Na(+) itself was not able to disturb the tautomeric equilibrium in the coordination complexes of the type Ura(Thy) tautomer:Na(+). To evaluate the DMSO effect, the CPCM solvation model was applied to triple complexes of the Ura tautomers. It appeared that in the solution there is coexistence between the diketo and enolic tautomers in a ratio of 53%:47%. This makes possible reconciliation of our experimental data. The biological significance of high-energy tautomers of nucleotide bases is discussed.
Collapse
Affiliation(s)
- Svitlana P Samijlenko
- Department of molecular and quantum biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, vul. Zabolotnoho150, Kyiv, 03143, Ukraine
| | | | | | | |
Collapse
|
40
|
Vogt N, Dorofeeva OV, Sipachev VA, Rykov AN. Molecular Structure of 9H-Adenine Tautomer: Gas-Phase Electron Diffraction and Quantum-Chemical Studies. J Phys Chem A 2009; 113:13816-23. [DOI: 10.1021/jp905755u] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Natalja Vogt
- Chemieinformationssysteme, Universität Ulm, 89069 Ulm, Germany, and Department of Chemistry, Moscow State University, Moscow 119992, Russia
| | - Olga V. Dorofeeva
- Chemieinformationssysteme, Universität Ulm, 89069 Ulm, Germany, and Department of Chemistry, Moscow State University, Moscow 119992, Russia
| | - Victor A. Sipachev
- Chemieinformationssysteme, Universität Ulm, 89069 Ulm, Germany, and Department of Chemistry, Moscow State University, Moscow 119992, Russia
| | - Anatolii N. Rykov
- Chemieinformationssysteme, Universität Ulm, 89069 Ulm, Germany, and Department of Chemistry, Moscow State University, Moscow 119992, Russia
| |
Collapse
|
41
|
Li D, Ai H. Catalysis Effects of Water Molecules and of Charge on Intramolecular Proton Transfer of Uracil. J Phys Chem B 2009; 113:11732-42. [PMID: 19645458 DOI: 10.1021/jp9031833] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Dejie Li
- School of Chemistry and Chemical Engineering, University of Jinan, 250022 P.R. China
| | - Hongqi Ai
- School of Chemistry and Chemical Engineering, University of Jinan, 250022 P.R. China
| |
Collapse
|
42
|
|
43
|
Platonov MO, Samijlenko SP, Sudakov OO, Kondratyuk IV, Hovorun DM. To what extent can methyl derivatives be regarded as stabilized tautomers of xanthine? SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2005; 62:112-4. [PMID: 16257701 DOI: 10.1016/j.saa.2004.12.012] [Citation(s) in RCA: 36] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/07/2004] [Revised: 12/01/2004] [Accepted: 12/02/2004] [Indexed: 05/05/2023]
Abstract
Taking xanthine (Xan) as an example, validity of an approach to experimental investigations of nucleotide bases' tautomeric equilibrium, based on the use of methyl derivatives corresponding to their prototropic tautomers, was studied by (1)H NMR in dimethylsulfoxide (DMSO) and by quantum chemical calculations at the B3LYP/6-311++G(d,p) level of theory. From (1)H NMR spectra of m(7)Xan, m(9)Xan, Xan and m(3)Xan conclusion was made that the N7H tautomeric forms of the last two compounds dominate in solution, which was supported by quantum chemical data. Calculated relative energies of the N9H tautomers of Xan and m(3)Xan (8.74 and 9.57 kcal/mol, accordingly) are rather close to the m(9)Xan one (9.11 kcal/mol). Nonspecific influence of DMSO modelled by the COSMO algorithm therewith reduces these values by approximately 2-3 kcal/mol. The data obtained imply that methyl derivatives are rather good models of high-energy tautomers of nucleotide bases, if their relative energies are not less than a few kcal/mol.
Collapse
Affiliation(s)
- M O Platonov
- Institute of Molecular Biology and Genetics, Ukrainian National Academy of Sciences, 150 vul. Zabolotnoho, Kyiv 03143, Ukraine
| | | | | | | | | |
Collapse
|