1
|
Dantsu Y, Zhang Y, Zhang W. Insight into the structures of unusual base pairs in RNA complexes containing a primer/template/adenosine ligand. RSC Chem Biol 2023; 4:942-951. [PMID: 37920395 PMCID: PMC10619131 DOI: 10.1039/d3cb00137g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Accepted: 08/29/2023] [Indexed: 11/04/2023] Open
Abstract
In the prebiotic RNA world, the self-replication of RNA without enzymes can be achieved through the utilization of 2-aminoimidazole activated nucleotides as efficient substrates. The mechanism of RNA nonenzymatic polymerization has been extensively investigated biophysically and structurally by using the model of an RNA primer/template complex which is bound by the imidazolium-bridged or triphosphate-bridged diguanosine intermediate. However, beyond the realm of the guanosine substrate, the structural insight into how alternative activated nucleotides bind and interact with the RNA primer/template complex remains unexplored, which is important for understanding the low reactivity of adenosine and uridine substrates in RNA primer extension, as well as its relationship with the structures. Here we use crystallography as a method and determine a series of high-resolution structures of RNA primer/template complexes bound by ApppG, a close analog of the dinucleotide intermediate containing adenosine and guanosine. The structures show that ApppG ligands bind to the RNA template through both Watson-Crick and noncanonical base pairs, with the primer 3'-OH group far from the adjacent phosphorus atom of the incoming substrate. The structures indicate that when adenosine is included in the imidazolium-bridged intermediate, the complexes are likely preorganized in a suboptimal conformation, making it difficult for the primer to in-line attack the substrate. Moreover, by co-crystallizing the RNA primer/template with chemically activated adenosine and guanosine monomers, we successfully observe the slow formation of the imidazolium-bridged intermediate (Ap-AI-pG) and the preorganized structure for RNA primer extension. Overall, our studies offer a structural explanation for the slow rate of RNA primer extension when using adenosine-5'-phosphoro-2-aminoimidazolide as a substrate during nonenzymatic polymerization.
Collapse
Affiliation(s)
- Yuliya Dantsu
- Department of Biochemistry and Molecular Biology, Indiana University School of Medicine Indianapolis IN 46202 USA
| | - Ying Zhang
- Department of Biochemistry and Molecular Biology, Indiana University School of Medicine Indianapolis IN 46202 USA
| | - Wen Zhang
- Department of Biochemistry and Molecular Biology, Indiana University School of Medicine Indianapolis IN 46202 USA
| |
Collapse
|
2
|
Odai K, Umesaki K. Kinetic Study of Transition Mutations from G-C to A-T Base Pairs in Watson-Crick DNA Base Pairs: Double Proton Transfers. J Phys Chem A 2021; 125:8196-8204. [PMID: 34516113 DOI: 10.1021/acs.jpca.1c05604] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
According to the Löwdin model [ Rev. Mod. Phys. 1963, 35, 724-732], the Watson-Crick guanine-cytosine (G-C) base pair is tautomerized (G*-C*) with a small probability and then replication of G*-C* produces G*-thymine (T) and adenine (A)-C* base pairs. On the basis of this model and our previous work [ J. Phys. Chem. B 2020, 124, 1715-1722], we first calculated the intrinsic reaction coordinates from G*-T to G-T* using density functional theory and evaluated the probability of G*-T tautomerization to G-T* by double proton transfer (DPT) on the basis of the transition state theory. Similarly, we calculated the probability of A-C* tautomerization to A*-C by DPT. Then, according to these probabilities, we calculated the probability of transition mutations from G-C to A-T after 2 replications. The calculated probability was 1.31 × 10-8, a value consistent with the mutation rate previously reported by Drake et al. [ Proc. Natl. Acad. Sci. U.S.A. 1991, 88, 7160-7164]. Our results suggest that DPT is one cause of the G-C → A-T transition. To investigate differences in the optical properties between G*-T and G-T* and between A-C* and A*-C, we also evaluated the infrared absorption spectra and Raman intensities for these base pairs.
Collapse
Affiliation(s)
- Kei Odai
- School of Science and Engineering, Kokushikan University, Setagaya-ku, Tokyo 154-8515, Japan
| | - Keisho Umesaki
- School of Science and Engineering, Kokushikan University, Setagaya-ku, Tokyo 154-8515, Japan
| |
Collapse
|
3
|
Brovarets' OO, Oliynyk TA, Hovorun DM. Novel Tautomerisation Mechanisms of the Biologically Important Conformers of the Reverse Löwdin, Hoogsteen, and Reverse Hoogsteen G *·C * DNA Base Pairs via Proton Transfer: A Quantum-Mechanical Survey. Front Chem 2019; 7:597. [PMID: 31620420 PMCID: PMC6759773 DOI: 10.3389/fchem.2019.00597] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2019] [Accepted: 08/12/2019] [Indexed: 11/22/2022] Open
Abstract
For the first time, in this study with the use of QM/QTAIM methods we have exhaustively investigated the tautomerization of the biologically-important conformers of the G*·C* DNA base pair-reverse Löwdin G*·C*(rWC), Hoogsteen G*'·C*(H), and reverse Hoogsteen G*'·C*(rH) DNA base pairs-via the single (SPT) or double (DPT) proton transfer along the neighboring intermolecular H-bonds. These tautomeric reactions finally lead to the formation of the novel G· C O 2 * (rWC), G N 2 * · C(rWC), G*'N2·C(rWC), G N 7 * · C(H), and G*'N7·C(rH) DNA base mispairs. Gibbs free energies of activation for these reactions are within the range 3.64-31.65 kcal·mol-1 in vacuum under normal conditions. All TSs are planar structures (Cs symmetry) with a single exception-the essentially non-planar transition state TSG*·C*(rWC)↔G+·C-(rWC) (C1 symmetry). Analysis of the kinetic parameters of the considered tautomerization reactions indicates that in reality only the reverse Hoogsteen G*'·C*(rH) base pair undergoes tautomerization. However, the population of its tautomerised state G*'N7·C(rH) amounts to an insignificant value-2.3·10-17. So, the G*·C*(rWC), G*'·C*(H), and G*'·C*(rH) base pairs possess a permanent tautomeric status, which does not depend on proton mobility along the neighboring H-bonds. The investigated tautomerization processes were analyzed in details by applying the author's unique methodology-sweeps of the main physical and chemical parameters along the intrinsic reaction coordinate (IRC). In general, the obtained data demonstrate the tautomeric mobility and diversity of the G*·C* DNA base pair.
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 Pharmacology, Bohomolets National Medical University, Kyiv, Ukraine
| | - Timothy A. Oliynyk
- 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
- Department of Pathophysiology, Bohomolets National Medical University, Kyiv, Ukraine
| |
Collapse
|
4
|
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
|
5
|
Brovarets' OO, Tsiupa KS, Dinets A, Hovorun DM. Unexpected Routes of the Mutagenic Tautomerization of the T Nucleobase in the Classical A·T DNA Base Pairs: A QM/QTAIM Comprehensive View. Front Chem 2018; 6:532. [PMID: 30538979 PMCID: PMC6277528 DOI: 10.3389/fchem.2018.00532] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2018] [Accepted: 10/12/2018] [Indexed: 01/24/2023] Open
Abstract
In this paper using quantum-mechanical (QM) calculations in combination with Bader's quantum theory of "Atoms in Molecules" (QTAIM) in the continuum with ε = 1, we have theoretically demonstrated for the first time that revealed recently highly-energetic conformers of the classical A·T DNA base pairs - Watson-Crick [A·T(wWC)], reverse Watson-Crick [A·T(wrWC)], Hoogsteen [A·T(wH)] and reverse Hoogsteen [A·T(wrH)] - act as intermediates of the intrapair mutagenic tautomerization of the T nucleobase owing to the novel tautomerisation pathways: A·T(wWC)↔A·T*(w⊥ WC); A·T(wrWC)↔A·T O 2 * (w⊥ rWC); A·T(wH)↔A·T*(w⊥ H); A·T(wrH)↔A·T O 2 * (w⊥ rH). All of them occur via the transition states as tight ion pairs (A+, protonated by the N6H2 amino group)·(T-, deprotonated by the N3H group) with quasi-orthogonal geometry, which are stabilized by the participation of the strong (A)N6+H···O4-/O2-(T) and (A)N6+H···N3-(T) H-bonds. Established tautomerizations proceed through a two-step mechanism of the protons moving in the opposite directions along the intermolecular H-bonds. Initially, proton moves from the N3H imino group of T to the N6H2 amino group of A and then subsequently from the protonated N6+H3 amino group of A to the O4/O2 oxygen atom of T, leading to the products - A·T*(w⊥ WC), A·T O 2 * (w⊥ rWC), A·T*(w⊥ H), and A·T O 2 * (w⊥ rH), which are substantially non-planar, conformationally-labile complexes. These mispairs are stabilized by the participation of the (A)N6H/N6H'···N3(T) and (T)O2H/O4H···N6(A) H-bonds, for which the pyramidalized amino group of A is their donor and acceptor. The Gibbs free energy of activation of these mutagenic tautomerizations lies in the range of 27.8-29.8 kcal·mol-1 at T = 298.15 K in the continuum with ε = 1.
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 Pharmacology, Bogomolets National Medical University, 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
| | - Andrii Dinets
- Department of Surgery #4, Bogomolets National Medical University, 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 Pathophysiology, Bogomolets National Medical University, Kyiv, Ukraine
| |
Collapse
|
6
|
Tolosa S, Sansón J, Hidalgo A. Theoretical study of mechanisms for double proton transfer in adenine–uracil base pair via steered molecular dynamic simulations. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.06.018] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
7
|
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
|
8
|
Jin L, Lv M, Zhao M, Wang R, Zhao C, Lu J, Wang L, Wang W, Wei Y. Formic acid catalyzed isomerization of protonated cytosine: a lower barrier reaction for tautomer production of potential biological importance. Phys Chem Chem Phys 2018; 19:13515-13523. [PMID: 28497833 DOI: 10.1039/c7cp01008g] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Tautomerism in nucleotide bases is one of the possible mechanisms of DNA mutation. In spite of numerous studies on the structure and energy of protonated cytosine tautomers, little information is available on the process of their intra- and intermolecular tautomerizations. The catalytic ability of H2O, HCOOH, and the HCOOHH2O group to facilitate the tautomerism of the Cyt2t+ to CytN3+ isomer has been studied. It is shown that the activation free energies of tautomerism in the gas phase are 161.17, 58.96, 26.06, and 15.69 kJ mol-1, respectively, when the reaction is carried out in the absence and presence of H2O, HCOOH, or the HCOOHH2O group. The formation of a doubly hydrogen bonded transition state is central to lowering the activation free energy and facilitating the intramolecular hydrogen atom transfer that is required for isomerization. In the aqueous phase, although the solvent effects of water significantly decrease the activation free energy of intramolecular tautomerization, the isomerization of the Cyt2t+ to CytN3+ isomer remains unfavorable, and the HCOOH and HCOOHH2O group mediated mechanisms are still more favorable. Meanwhile, conventional transition state theory (CTST) followed by Wigner tunneling correction is then applied to estimate the rate constants. The rate constant with Wigner tunneling correction for direct tautomerization is obviously smaller than that of HCOOH-mediated tautomerization, which is the most plausible mechanism. Finally, another important finding is that the product complex (CytN3+HCOOH) is in the rapid tautomeric equilibrium with the reaction complex (Cyt2t+HCOOH) (τ99.9% = 3.84 × 10-12 s), which is implemented by the mechanism of the concerted synchronous double proton transfer. Its lifetime of the formed CytN3+HCOOH complex (τ = 8.33 × 10-9 s) is almost one order of magnitude larger than the time required for the replication machinery to forcibly dissociate a base pair into the monomers during DNA replication (several ns), which is further dissociated into the CytN3+ and HCOOH monomers. The results of the present study demonstrate the feasibility of acid catalysis for DNA base isomerization reactions that would otherwise be forbidden.
Collapse
Affiliation(s)
- Lingxia Jin
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi Sci-Tech University, Hanzhong 723001, China.
| | - Mengdan Lv
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi Sci-Tech University, Hanzhong 723001, China.
| | - Mengting Zhao
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi Sci-Tech University, Hanzhong 723001, China.
| | - Rui Wang
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi Sci-Tech University, Hanzhong 723001, China.
| | - Caibin Zhao
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi Sci-Tech University, Hanzhong 723001, China.
| | - Jiufu Lu
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi Sci-Tech University, Hanzhong 723001, China.
| | - Ling Wang
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi Sci-Tech University, Hanzhong 723001, China.
| | - Wenliang Wang
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an 710062, China.
| | - Yawen Wei
- Institute of publication Science, Chang'an University, Xi'an 710064, China
| |
Collapse
|
9
|
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
|
10
|
Jin L, Song X, Cao Z, Luo L, Zhao C, Lu J, Zhang Q. The isomerization of cytosine: Intramolecular hydrogen atom transfer mediated through formic acid. J PHYS ORG CHEM 2018. [DOI: 10.1002/poc.3831] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Lingxia Jin
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science; Shaanxi University of Technology; Hanzhong China
| | - Xiaoling Song
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science; Shaanxi University of Technology; Hanzhong China
| | - Zhe Cao
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science; Shaanxi University of Technology; Hanzhong China
| | - LiYang Luo
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science; Shaanxi University of Technology; Hanzhong China
| | - Caibin Zhao
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science; Shaanxi University of Technology; Hanzhong China
| | - Jiufu Lu
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science; Shaanxi University of Technology; Hanzhong China
| | - Qiang Zhang
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science; Shaanxi University of Technology; Hanzhong China
| |
Collapse
|
11
|
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
|
12
|
Tolosa S, Sansón JA, Hidalgo A. Theoretical study of mechanisms for the hydrolytic deamination of cytosine via steered molecular dynamic simulations. RSC Adv 2018; 8:34867-34876. [PMID: 35547048 PMCID: PMC9087476 DOI: 10.1039/c8ra07390b] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2018] [Accepted: 10/04/2018] [Indexed: 12/31/2022] Open
Abstract
Gibbs free energy profiles of the cytosine deamination assisted by a water molecule in a discrete aqueous medium were obtained by the application of Steered Molecular Dynamic (SMD) simulations. Two pathways were considered to explain the mechanism of this process, where the water molecule attacks the C–N bond to give an intermediate (an amino–hydroxy–oxo structure in the A-path, and a hydroxy–oxo in the B-path) as the determinant step of reaction. Stationary structures along both energy profiles were analyzed at molecular dynamics level, obtaining states with higher free energies than those from electronic calculations in the gas phase and in solution described as a continuous medium. From the results obtained, the more complex A-pathway, with five steps, was kinetically the most favorable (with an endergonic reaction energy of 7.41 kcal mol−1, a high barrier of 67.53 kcal mol−1, and a small velocity constant k2 = 1.80 × 10−37 s−1), concluding that the uracil base can participate in a spontaneous genetic mutation since the uracil–ammonia complex has a long lifetime of 6.10 × 1027 s. This process turns out exergonic and faster when carried out in gas phase simulation or electronic calculation with a continuous medium, due to the disappearance of explicit water molecules that can compete with the assistant molecule. Gibbs free energy profiles of the cytosine deamination assisted by a water molecule in a discrete aqueous medium were obtained by the application of Steered Molecular Dynamic (SMD) simulations.![]()
Collapse
Affiliation(s)
- S. Tolosa
- Departamento de Ingeniería Química y Química Física
- Universidad de Extremadura
- Badajoz
- Spain
| | - J. A. Sansón
- Departamento de Ingeniería Química y Química Física
- Universidad de Extremadura
- Badajoz
- Spain
| | - A. Hidalgo
- Departamento de Ingeniería Química y Química Física
- Universidad de Extremadura
- Badajoz
- Spain
| |
Collapse
|
13
|
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
|
14
|
Abdalla S, Obaid A, Al-Marzouki FM. Effects of Environmental Factors and Metallic Electrodes on AC Electrical Conduction Through DNA Molecule. NANOSCALE RESEARCH LETTERS 2017; 12:316. [PMID: 28454482 PMCID: PMC5407417 DOI: 10.1186/s11671-017-2076-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/23/2016] [Accepted: 04/13/2017] [Indexed: 06/07/2023]
Abstract
BACKGROUND Deoxyribonucleic acid (DNA) is one of the best candidate materials for various device applications such as in electrodes for rechargeable batteries, biosensors, molecular electronics, medical- and biomedical-applications etc. Hence, it is worthwhile to examine the mechanism of charge transport in the DNA molecule, however, still a question without a clear answer is DNA a molecular conducting material (wire), semiconductor, or insulator? The answer, after the published data, is still ambiguous without any confirmed and clear scientific answer. DNA is found to be always surrounded with different electric charges, ions, and dipoles. These surrounding charges and electric barrier(s) due to metallic electrodes (as environmental factors (EFs)) play a substantial role when measuring the electrical conductivity through λ-double helix (DNA) molecule suspended between metallic electrodes. We found that strong frequency dependence of AC-complex conductivity comes from the electrical conduction of EFs. This leads to superimposing serious incorrect experimental data to measured ones. METHODS At 1 MHz, we carried out a first control experiment on electrical conductivity with and without the presence of DNA molecule. If there are possible electrical conduction due to stray ions and contribution of substrate, we will detected them. This control experiment revealed that there is an important role played by the environmental-charges around DNA molecule and any experiment should consider this role. RESULTS AND DISCUSSION We have succeeded to measure both electrical conductivity due to EFs (σ ENV) and electrical conductivity due to DNA molecule (σ DNA) independently by carrying the measurements at different DNA-lengths and subtracting the data. We carried out measurements as a function of frequency (f) and temperature (T) in the ranges 0.1 Hz < f < 1 MHz and 288 K < T < 343 K. The measured conductivity (σ MES) portrays a metal-like behavior at high frequencies near 1 MHz. However, we found that σ DNA was far from this behavior because the conduction due to EFs superimposes σ DNA, in particular at low frequencies. By measuring the electrical conductivity at different lengths: 40, 60, 80, and 100 nm, we have succeeded not only to separate the electrical conduction of the DNA molecule from all EFs effects that surround the molecule, but also to present accurate values of σ DNA and the dielectric constant of the molecule ε'DNA as a function of temperature and frequency. Furthermore, in order to explain these data, we present a model describing the electrical conduction through DNA molecule: DNA is a classical semiconductor with charges, dipoles and ions that result in creation of localized energy-states (LESs) in the extended bands and in the energy gap of the DNA molecule. CONCLUSIONS This model explains clearly the mechanism of charge transfer mechanism in the DNA, and it sheds light on why the charge transfer through the DNA can lead to insulating, semiconducting, or metallic behavior on the same time. The model considers charges on DNA, in the extended bands, either could be free to move under electric field or localized in potential wells/hills. Localization of charges in DNA is an intrinsic structural-property of this solitaire molecule. At all temperatures, the expected increase in thermal-induced charge is attributed to the delocalization of holes (or/and electrons) in potential hills (or/and potential wells) which accurately accounts for the total electric and dielectric behavior through DNA molecule. We succeeded to fit the experimental data to the proposed model with reasonable magnitudes of potential hills/wells that are in the energy range from 0.068 eV.
Collapse
Affiliation(s)
- S. Abdalla
- Department of Physics, Faculty of Science, King Abdulaziz University Jeddah, P.O. Box 80203, Jeddah, 21589 Saudi Arabia
| | - A. Obaid
- Department of Chemistry, Faculty of Science, King Abdulaziz University Jeddah, P.O. Box 80203, Jeddah, 21589 Saudi Arabia
| | - F. M. Al-Marzouki
- Department of Physics, Faculty of Science, King Abdulaziz University Jeddah, P.O. Box 80203, Jeddah, 21589 Saudi Arabia
| |
Collapse
|
15
|
Zhen JP, Wei XC, Shi WJ, Huang ZY, Jin B, Zhou YK. Cooperativity effect involving drug-DNA/RNA intermolecular interaction: A B3LYP-D3 and MP2 theoretical investigation on ketoprofen⋯cytosine⋯H 2O system. J Biomol Struct Dyn 2017; 36:3587-3606. [PMID: 29092677 DOI: 10.1080/07391102.2017.1400469] [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] [Indexed: 10/18/2022]
Abstract
In order to examine the origin of the drug action and design new DNA/RNA-targeted drugs, the cooperativity effect involving drug-DNA/RNA intermolecular interaction in ketoprofen⋯cytosine⋯H2O ternary system were investigated by the B3LYP, B3LYP-D3, and MP2 methods with the 6-311++G(2d,p) basis set. The thermodynamic cooperativity was also evaluated at 310.15 K. The N-H⋯O, O-H⋯O, O-H⋯N, C-H⋯N, and C-H⋯O H bonds coexist in ternary complexes. The intermolecular interactions obtained by B3LYP-D3 are close to those calculated by MP2. The steric effects and van der Waals interactions have little influence on the cooperativity effects. The anti-cooperativity effect in ket⋯cyt⋯H2O is far more notable than the cooperativity effect, and the stability of the cyclic structure with anti-cooperativity effect is higher than that of the linear structure with cooperativity effect, as is confirmed by the AIM (atoms in molecules) and RDG (reduced density gradient) analysis. Thus, it can be inferred that, in the presence of H2O, the anti-cooperativity effect plays a dominant role in the drug-DNA/RNA interaction, and the nature of the hydration in the binding of drugs to DNA/RNA bases is the H-bonding anti-cooperativity effect. Furthermore, the drug always links simultaneously with DNA/RNA base and H2O, and only in this way can the biological activity of drugs play a role. In most cases, the enthalpy change is the major factor driving the cooperativity, as is different from most of biomacromolecule complexes.
Collapse
Affiliation(s)
- Jun-Ping Zhen
- a Molecular Imaging Laboratory, Department of Radiology , The Second Hospital of Shanxi Medical University , Taiyuan 030053 , China
| | - Xiao-Chun Wei
- a Molecular Imaging Laboratory, Department of Radiology , The Second Hospital of Shanxi Medical University , Taiyuan 030053 , China
| | - Wen-Jing Shi
- a Molecular Imaging Laboratory, Department of Radiology , The Second Hospital of Shanxi Medical University , Taiyuan 030053 , China
| | - Zhu-Yuan Huang
- a Molecular Imaging Laboratory, Department of Radiology , The Second Hospital of Shanxi Medical University , Taiyuan 030053 , China
| | - Bo Jin
- a Molecular Imaging Laboratory, Department of Radiology , The Second Hospital of Shanxi Medical University , Taiyuan 030053 , China
| | - Yu-Kun Zhou
- a Molecular Imaging Laboratory, Department of Radiology , The Second Hospital of Shanxi Medical University , Taiyuan 030053 , China
| |
Collapse
|
16
|
Ebrahimi S, Dabbagh HA, Eskandari K. Arrangement and nature of intermolecular hydrogen bonding in complex biomolecular systems: modeling the vitamin C---L-alanine interaction. Struct Chem 2017. [DOI: 10.1007/s11224-017-1046-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
17
|
The role of intramolecular H-bonds predominant effects in myricetin higher antioxidant activity. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.06.020] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
18
|
Brovarets' OO, Voiteshenko IS, Pérez-Sánchez H, Hovorun DM. A QM/QTAIM research under the magnifying glass of the DPT tautomerisation of the wobble mispairs involving 2-aminopurine. NEW J CHEM 2017. [DOI: 10.1039/c7nj00717e] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this study, a comprehensive survey of the changes of the physico-chemical parameters at each point of the IRC for the biologically important T·2AP*(w) ↔ T*·2AP(w) and G·2AP*(w) ↔ G*·2AP(w) DPT tautomerisation reactions involved in the point mutations (transitions and transversions) induced by 2-aminopurine (2AP) in DNA is provided.
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
| | - Ivan S. Voiteshenko
- Department of Molecular Biotechnology and Bioinformatics
- Institute of High Technologies
- Taras Shevchenko National University of Kyiv
- 03022 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)
- 30107 Guadalupe (Murcia)
- Spain
| | - 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
|
19
|
Palanivel U, Lakshmipathi S. Hydrogen bonds in Zif268 proteins – a theoretical perspective. J Biomol Struct Dyn 2016; 34:1607-24. [DOI: 10.1080/07391102.2015.1085903] [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]
Affiliation(s)
- Umadevi Palanivel
- Department of Physics, Bharathiar University, Coimbatore 641 046, Tamil Nadu, India
| | | |
Collapse
|
20
|
Vega-Rodríguez S, Jiménez-Cataño R, Leyva E. Tautomerism in substituted pyridofuroxans: A theoretical study. COMPUT THEOR CHEM 2016. [DOI: 10.1016/j.comptc.2016.06.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
21
|
Özdemir N. Quantum chemical study of tautomerism in 2-[(4-phenylthiazol-2-yl)hydrazonomethyl]phenol. COMPUT THEOR CHEM 2016. [DOI: 10.1016/j.comptc.2016.04.013] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
|
22
|
Markova N, Pejov L, Stoyanova N, Enchev V. Hybrid MC/QC simulations of water-assisted proton transfer in nucleosides. Guanosine and its analog acyclovir. J Biomol Struct Dyn 2016; 35:1168-1188. [PMID: 27092850 DOI: 10.1080/07391102.2016.1179594] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Abstract
To provide an in-depth insight into the molecular basis of spontaneous tautomerism in DNA and RNA base pairs, a hybrid Monte Carlo (MC)-quantum chemical (QC) methodology is implemented to map two-dimensional potential energy surfaces along the reaction coordinates of solvent-assisted proton transfer processes in guanosine and its analog acyclovir in aqueous solution. The solvent effects were simulated by explicit inclusion of water molecules that model the relevant part of the first hydration shell around the solute. The position of these water molecules was estimated by carrying out a classical Metropolis Monte Carlo simulation of dilute water solutions of the guanosine (Gs) and acyclovir (ACV) and subsequently analyzing solute-solvent intermolecular interactions in the statistically-independent MC-generated configurations. The solvent-assisted proton transfer processes were further investigated using two different ab initio MP2 quantum chemical approaches. In the first one, potential energy surfaces of the 'bare' finite solute-solvent clusters containing Gs/ACV and four water molecules (MP2/6-31+G(d,p) level) were explored, while within the second approach, these clusters were embedded in 'bulk' solvent treated as polarizable continuum (C-PCM/MP2/6-31+G(d,p) level of theory). It was found that in the gas phase and in water solution, the most stable tautomer for guanosine and acyclovir is the 1H-2-amino-6-oxo form followed by the 2-amino-6-(sZ)-hydroxy form. The energy barriers of the water-assisted proton transfer reaction in guanosine and in acyclovir are found to be very similar - 11.74 kcal mol-1 for guanosine and 11.16 kcal mol-1 for acyclovir, and the respective rate constants (k = 1.5 × 101 s-1, guanosine and k = 4.09 × 101 s-1, acyclovir), are sufficiently large to generate the 2-amino-6-(sZ)-hydroxy tautomer. The analysis of the reaction profiles in both compounds shows that the proton transfer processes occur through the asynchronous concerted mechanism.
Collapse
Affiliation(s)
- Nadezhda Markova
- a Institute of Organic Chemistry , Bulgarian Academy of Sciences , 1113 Sofia , Bulgaria
| | - Ljupco Pejov
- b Faculty of Natural Sciences and Mathematics , Institute of Chemistry , Skopje , Macedonia.,c Research Centre for Environment and Materials , Macedonian Academy of Sciences and Arts , Krste Misirkov 2, 1000 Skopje , Macedonia
| | - Nina Stoyanova
- a Institute of Organic Chemistry , Bulgarian Academy of Sciences , 1113 Sofia , Bulgaria
| | - Venelin Enchev
- a Institute of Organic Chemistry , Bulgarian Academy of Sciences , 1113 Sofia , Bulgaria
| |
Collapse
|
23
|
Brovarets' OO, Hovorun DM. A novel conception for spontaneous transversions caused by homo-pyrimidine DNA mismatches: a QM/QTAIM highlight. Phys Chem Chem Phys 2016. [PMID: 26219928 DOI: 10.1039/c5cp03211c] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
We have firstly shown that the T·T(w) and C·C(w) DNA mismatches with wobble (w) geometry stay in slow tautomeric equilibrium with short T·T*(WC) and C·C*(WC) Watson-Crick (WC) mispairs. These non-dissociative tautomeric rearrangements are controlled by the plane-symmetric, highly stable, highly polar and zwitterionic transition states. The obtained results allow us to understand in what way the T·T(w) and C·C(w) mismatches acquire enzymatically competent T·T*(WC) and C·C*(WC) conformations directly in the hydrophobic recognition pocket of a high-fidelity DNA-polymerase, thereby producing thermodynamically non-equilibrium spontaneous transversions. The simplest numerical estimation of the frequency ratio of the TT to CC spontaneous transversions satisfactorily agrees with experimental data.
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
|
24
|
Optical absorption and emission properties of benzene-expanded Janus AT nucleobase analogues: A DFT study. Struct Chem 2016. [DOI: 10.1007/s11224-016-0743-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
25
|
Brovarets' OO, Hovorun DM. Tautomeric transition between wobble A·C DNA base mispair and Watson-Crick-like A·C* mismatch: microstructural mechanism and biological significance. Phys Chem Chem Phys 2016; 17:15103-10. [PMID: 25994250 DOI: 10.1039/c5cp01568e] [Citation(s) in RCA: 60] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Here, we use MP2/DFT quantum-chemical methods combined with Quantum Theory of Atoms in Molecules to study the tautomeric transition between wobble A·C(w) mismatch and Watson-Crick-like A·C*(WC) base mispair, proceeding non-dissociatively via sequential proton transfer between bases through the planar, highly stable and zwitterionic TS(A∙C-)(A∙C(W)<-->A∙C&(WC)) transition state joined by the participation of (A)N6(+)H∙∙∙N4(-)(C), (A)N1(+)H∙∙∙N4(-)(C) and (A)C2(+)H∙∙∙N3(-)(C) H-bonds. Notably, the A·C(w) ↔ A·C*(WC) tautomerization reaction is accompanied by 10 unique patterns of the specific intermolecular interactions that consistently replace each other. Our data suggest that biologically significant A·C(w) → A·C*(WC) tautomerization is a kinetically controlled pathway for formation of the enzymatically competent Watson-Crick-like A·C*(WC) DNA base mispair in the essentially hydrophobic recognition pocket of the high-fidelity DNA-polymerase, responsible for the occurrence of spontaneous point AC/CA incorporation errors during DNA biosynthesis.
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
|
26
|
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
|
27
|
Brovarets' OO, Pérez-Sánchez H, Hovorun DM. Structural grounds for the 2-aminopurine mutagenicity: a novel insight into the old problem of the replication errors. RSC Adv 2016. [DOI: 10.1039/c6ra17787e] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Mutagenic pressure of the 2AP molecule on DNA during its replication is realized via the more intensive generation of the T* mutagenic tautomers through the reaction 2AP·T(WC) → 2AP·T*(w).
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
| | - 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
|
28
|
Cardamone S, Caine BA, Blanch E, Lizio MG, Popelier PLA. The computational prediction of Raman and ROA spectra of charged histidine tautomers in aqueous solution. Phys Chem Chem Phys 2016; 18:27377-27389. [DOI: 10.1039/c6cp05744f] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Histidine is a key component of a number of enzymatic mechanisms, and undertakes many functionalities in biochemical systems.
Collapse
Affiliation(s)
- Salvatore Cardamone
- Manchester Institute of Biotechnology (MIB)
- Manchester M1 7DN
- UK
- School of Chemistry
- University of Manchester
| | - Beth A. Caine
- Manchester Institute of Biotechnology (MIB)
- Manchester M1 7DN
- UK
- School of Chemistry
- University of Manchester
| | - Ewan Blanch
- School of Science
- Royal Melbourne Institute of Technology (RMIT)
- Melbourne
- Australia
| | - Maria G. Lizio
- Manchester Institute of Biotechnology (MIB)
- Manchester M1 7DN
- UK
- School of Chemistry
- University of Manchester
| | - Paul L. A. Popelier
- Manchester Institute of Biotechnology (MIB)
- Manchester M1 7DN
- UK
- School of Chemistry
- University of Manchester
| |
Collapse
|
29
|
By how many tautomerisation routes the Watson–Crick-like A·C* DNA base mispair is linked with the wobble mismatches? A QM/QTAIM vision from a biological point of view. Struct Chem 2015. [DOI: 10.1007/s11224-015-0687-4] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
30
|
Das D, Dutta A, Mondal P. Interaction of aquated form of ruthenium(III) anticancer complexes with normal and mismatch base pairs: A density functional theoretical study. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2015.08.020] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
31
|
Sukker GM, Elroby SA, Hilal R. Gas-phase acidity and dynamics of the protonation processes of carbidopa and levodopa. A QM/QD study. J Biomol Struct Dyn 2015; 34:2268-80. [PMID: 26511889 DOI: 10.1080/07391102.2015.1113385] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Abstract
The present work details, our efforts to investigate and compare the acid-base properties of levodopa (LD) and carbidopa (CD), a drug combination being used in the treatment of Parkinson's disease. Protonation and deprotonation were examined for all possible sites in the two drugs. Proton affinity and proton detachment enthalpies were computed at the B3LYP/6-311++G** level of theory. Results of the present work reveal that CD is more basic and can abstract protons in solution much more efficiently than LD. Furthermore, at all deportation sites considered, CD is more acidic than LD. DFT-based ADMP, dynamic simulations have been performed to explore the dynamics of the protonation processes in LD and CD. Thus, while the dynamics of the protonation process of LD is very straightforward leading to the formation of the corresponding cation, the protonation process in CD is very complex involving a major geometry change and rearrangement. Results of the present work reveal that the active species in acid medium is not CD in its normal geometry but a carbonyl hydrazine form instead. The presence of the carbonyl group β to the hydrazine group may very well underlie its enhanced activity which allows it to bind to the active site of the DDC enzyme. The relative stabilities of various water-water-CD complexes have been computed and compared.
Collapse
Affiliation(s)
- Ghader M Sukker
- a Faculty of Science, Chemistry Department , King Abdulaziz University , Jeddah , Saudi Arabia
| | - Shaaban A Elroby
- a Faculty of Science, Chemistry Department , King Abdulaziz University , Jeddah , Saudi Arabia.,c Faculty of Science, Chemistry Department , Beni-Suef University , Beni-Suef , Egypt
| | - Rifaat Hilal
- a Faculty of Science, Chemistry Department , King Abdulaziz University , Jeddah , Saudi Arabia.,b Faculty of Science, Chemistry Department , Cairo University , Giza 12613 , Egypt
| |
Collapse
|
32
|
Brovarets' OO, Hovorun DM. Proton tunneling in the A∙T Watson-Crick DNA base pair: myth or reality? J Biomol Struct Dyn 2015; 33:2716-20. [PMID: 26362836 DOI: 10.1080/07391102.2015.1092886] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The results and conclusions reached by Godbeer et al. in their recent work, that proton tunneling in the A∙T(WC) Watson-Crick (WC) DNA base pair occurs according to the Löwdin's (L) model, but with a small (~10(-9)) probability were critically analyzed. Here, it was shown that this finding overestimates the possibility of the proton tunneling at the A∙T(WC)↔A*∙T*(L) tautomerization, because this process cannot be implemented as a chemical reaction. Furthermore, it was outlined those biologically important nucleobase mispairs (A∙A*↔A*∙A, G∙G*↔G*∙G, T∙T*↔T*∙T, C∙C*↔C*∙C, H∙H*↔H*∙H (H - hypoxanthine)) - the players in the field of the spontaneous point mutagenesis - where the tunneling of protons is expected and for which the application of the model proposed by Godbeer et al. can be productive.
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.,b Department of Molecular Biotechnology and Bioinformatics , Institute of High Technologies, Taras Shevchenko National University of Kyiv , 2-h Akademika Hlushkova Ave., 03022 Kyiv , Ukraine
| | - Dmytro M Hovorun
- 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.,b Department of Molecular Biotechnology and Bioinformatics , Institute of High Technologies, Taras Shevchenko National University of Kyiv , 2-h Akademika Hlushkova Ave., 03022 Kyiv , Ukraine
| |
Collapse
|
33
|
|
34
|
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
|
35
|
Structural significance of modified nucleoside 5-taurinomethyl-2-thiouridine, τm5s2U, found at ‘wobble’ position in anticodon loop of human mitochondrial tRNALys. Struct Chem 2015. [DOI: 10.1007/s11224-015-0642-4] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
36
|
|
37
|
|
38
|
Brovarets' OO, Hovorun DM. How many tautomerization pathways connect Watson-Crick-like G*·T DNA base mispair and wobble mismatches? J Biomol Struct Dyn 2015; 33:2297-315. [PMID: 25932960 DOI: 10.1080/07391102.2015.1046936] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
In this study, we have theoretically demonstrated the intrinsic ability of the wobble G·T(w)/G*·T*(w)/G·T(w1)/G·T(w2) and Watson-Crick-like G*·T(WC) DNA base mispairs to interconvert into each other via the DPT tautomerization. We have established that among all these transitions, only one single G·T(w) ↔ G*·T(WC) pathway is eligible from a biological perspective. It involves short-lived intermediate - the G·T*(WC) base mispair - and is governed by the planar, highly stable, and zwitterionic [Formula: see text] transition state stabilized by the participation of the unique pattern of the five intermolecular O6(+)H⋯O4(-), O6(+)H⋯N3(-), N1(+)H⋯N3(-), N1(+)H⋯O2(-), and N2(+)H⋯O2(-) H-bonds. This non-dissociative G·T(w) ↔ G*·T(WC) tautomerization occurs without opening of the pair: Bases within mispair remain connected by 14 different patterns of the specific intermolecular interactions that successively change each other along the IRC. Novel kinetically controlled mechanism of the thermodynamically non-equilibrium spontaneous point GT/TG incorporation errors has been suggested. The mutagenic effect of the analogues of the nucleotide bases, in particular 5-bromouracil, can be attributed to the decreasing of the barrier of the acquisition by the wobble pair containing these compounds of the enzymatically competent Watson-Crick's geometry via the intrapair mutagenic tautomerization directly in the essentially hydrophobic recognition pocket of the replication DNA-polymerase machinery. Proposed approaches are able to explain experimental data, namely growth of the rate of the spontaneous point incorporation errors during DNA biosynthesis with increasing temperature.
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., Kyiv 03680 , Ukraine.,b Department of Molecular Biotechnology and Bioinformatics , Institute of High Technologies, Taras Shevchenko National University of Kyiv , 2-h Akademika Hlushkova Ave., Kyiv 03022 , Ukraine
| | - Dmytro M Hovorun
- a Department of Molecular and Quantum Biophysics , Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine , 150 Akademika Zabolotnoho Str., Kyiv 03680 , Ukraine.,b Department of Molecular Biotechnology and Bioinformatics , Institute of High Technologies, Taras Shevchenko National University of Kyiv , 2-h Akademika Hlushkova Ave., Kyiv 03022 , Ukraine
| |
Collapse
|
39
|
Brovarets' OO, Hovorun DM. How does the long G·G* Watson-Crick DNA base mispair comprising keto and enol tautomers of the guanine tautomerise? The results of a QM/QTAIM investigation. Phys Chem Chem Phys 2015; 16:15886-99. [PMID: 24964351 DOI: 10.1039/c4cp01241k] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The double proton transfer (DPT) in the long G·G* Watson-Crick base mispair (|C6N1(G*)N1C6(G)| = 36.4°; C1 symmetry), involving keto and enol tautomers of the guanine (G) nucleobase, along two intermolecular neighboring O6H···O6 (8.39) and N1···HN1 (6.14 kcal mol(-1)) H-bonds that were established to be slightly anti-cooperative, leads to its transformation into the G*·G base mispair through a single transition state (|C6N1N1C6| = 37.1°; C1), namely to the interconversion into itself. It was shown that the G·G* ↔ G*·G tautomerisation via the DPT is assisted by the third specific contact, that sequentially switches along the intrinsic reaction coordinate (IRC) in an original way: (G)N2H···N2(G*) H-bond (-25.13 to -10.37) → N2···N2 van der Waals contact (-10.37 to -9.23) → (G)N2···HN2(G*) H-bond (-9.23 to 0.79) → (G*)N2···HN2(G) H-bond (0.79 to 7.35 Bohr). The DPT tautomerisation was found to proceed through the asynchronous concerted mechanism by employing the QM/QTAIM approach and the methodology of the scans of the geometric, electron-topological, energetic, polar and NBO properties along the IRC. Nine key points, that can be considered as part of the tautomerisation repertoire, have been established and analyzed in detail. Furthermore, it was shown that the G·G* or G*·G base mispair is a thermodynamically and dynamically stable structure with a lifetime of 8.22 × 10(-10) s and all 6 low-frequency intermolecular vibrations are able to develop during this time span. Lastly, our results highlight the importance of the G·G* ↔ G*·G DPT tautomerisation, which can have implications for biological and chemical sensing applications.
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
|
40
|
Arba M, Kartasasmita RE, Tjahjono DH. Molecular docking and dynamics simulations on the interaction of cationic porphyrin-anthraquinone hybrids with DNA G-quadruplexes. J Biomol Struct Dyn 2015; 34:427-38. [PMID: 25808513 DOI: 10.1080/07391102.2015.1033015] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
A series of cationic porphyrin-anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.
Collapse
Affiliation(s)
- Muhammad Arba
- a School of Pharmacy , Bandung Institute of Technology , Jalan Ganesha 10, Bandung 40132 , Indonesia.,b Department of Chemistry , Halu Oleo University , Jl. HEA Mokodompit, Kendari 93232 , Indonesia
| | - Rahmana E Kartasasmita
- a School of Pharmacy , Bandung Institute of Technology , Jalan Ganesha 10, Bandung 40132 , Indonesia
| | - Daryono H Tjahjono
- a School of Pharmacy , Bandung Institute of Technology , Jalan Ganesha 10, Bandung 40132 , Indonesia
| |
Collapse
|
41
|
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
|
42
|
Arslan NB, Özdemir N. Direct and solvent-assisted keto-enol tautomerism and hydrogen-bonding interactions in 4-(m-chlorobenzylamino)-3-phenyl-4,5-dihydro-1H-1,2,4-triazol-5-one: a quantum-chemical study. J Mol Model 2015; 21:19. [PMID: 25617212 DOI: 10.1007/s00894-015-2574-8] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2014] [Accepted: 01/04/2015] [Indexed: 10/24/2022]
Abstract
The tautomeric equilibrium of the title triazole compound was computationally analyzed at the B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) levels of theory. The solvent effect was considered for three solvents (chloroform, methanol, and water). Two distinct mechanisms were applied: a direct intramolecular transfer using the polarizable continuum model (PCM) and a solvent-assisted mechanism. The calculations indicated that the keto form is more stable in all cases. It was found that the barrier heights for the tautomerization reaction are very high, indicating a relatively disfavored process. Although the barrier heights for solvent-assisted reactions are significantly lower than those for the unassisted tautomerization reaction, implying the importance of the superior catalytic effect of the solvents, monosolvation was not found to be sufficient for the reaction to occur. Finally, the two intermolecular hydrogen-bonding interactions in the crystal structure were investigated in the gas phase; according to the calculated energies and structural parameters, the order of stability is N3-H3···O1 > N1-H1···O1.
Collapse
Affiliation(s)
- N Burcu Arslan
- Department of Computer Education and Instructional Technology, Faculty of Education, Giresun University, 28100, Giresun, Turkey
| | | |
Collapse
|
43
|
Chi S, Xie W, Zhang J, Xu S. Theoretical insight into the structural mechanism for the binding of vinblastine with tubulin. J Biomol Struct Dyn 2015; 33:2234-54. [PMID: 25588192 DOI: 10.1080/07391102.2014.999256] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Vinblastine (VLB) is one of vinca alkaloids with high cytotoxicity toward cancer cells approved for clinical use. However, because of drug resistance, toxicity, and other side effects caused from the use of VLB, new vinca alkaloids with higher cytotoxicity toward cancer cells and other good qualities need to develop. One strategy is to further study and better understand the essence why VLB possesses the high cytotoxicity toward cancer cells. In present work, by using molecular simulation, molecular docking, density functional calculation, and the crystal structure of α,β-tubulin complex, we find two modes labeled in catharanthine moiety (CM) and vindoline moiety (VM) modes of VLB bound with the interface of α,β-tubulin to probe the essence why VLB has the high cytotoxicity toward cancer cells. In the CM mode, nine key residues B-Ser178, B-Asp179, B-Glu183, B-Tyr210, B-Asp226, C-Lys326, C-Asp327, C-Lys336, and C-Lys352 from the α,β-tubulin complex are determined as the active sites for the interaction of VLB with α,β-tubulin. Some of them such as B-Ser178, B-Glu183, B-Tyr210, B-Asp226, C-Lys326, C-Asp327, and C-Lys336 are newly identified as the active sites in present work. The affinity between VLB and the active pocket within the interface of α,β-tubulin is -60.8 kJ mol(-1) in the CM mode. In the VM mode, that is a new mode established in present paper, nine similar key residues B-Lys176, B-Ser178, B-Asp179, B-Glu183, B-Tyr210, B-Asp226, C-Lys326, C-Asp327, and C-Lys336 from the α,β-tubulin complex are found as the active sites for the interaction with VLB. The difference is from one key residue C-Lys352 in the CM mode changed to the key residue B-Lys176 in the VM mode. The affinity between VLB and the active pocket within the interface of α,β-tubulin is -96.3 kJ mol(-1) in the VM mode. Based on the results obtained in present work, and because VLB looks like two faces, composed of CM and VM both to have similar polar active groups, to interact with the active sites, we suggest double-faces sticking mechanism for the binding of VLB to the interface of α,β-tubulin. The double-faces sticking mechanism can be used to qualitatively explain high cytotoxicity toward cancer cells of vinca alkaloids including vinblastine, vincristine, vindestine, and vinorelbine approved for clinical use and vinflunine still in a phase III clinical trial. Furthermore, this mechanism will be applied to develop novel vinca alkaloids with much higher cytotoxicity toward cancer cells.
Collapse
Affiliation(s)
- Shaoming Chi
- a Key Laboratory of Education Ministry for Medicinal Chemistry of Natural Resource , College of Chemical Science and Technology, Yunnan University , Kunming 650091 , China
| | | | | | | |
Collapse
|
44
|
Kato RB, Silva FT, Pappa GL, Belchior JC. Genetic algorithms coupled with quantum mechanics for refinement of force fields for RNA simulation: a case study of glycosidic torsions in the canonical ribonucleosides. Phys Chem Chem Phys 2015; 17:2703-14. [DOI: 10.1039/c4cp03779k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We report the use of genetic algorithms (GA) as a method to refine force field parameters in order to determine RNA energy.
Collapse
Affiliation(s)
- Rodrigo B. Kato
- Department of Computer Science
- Universidade Federal de Minas Gerais
- Belo Horizonte
- Brazil
| | - Frederico T. Silva
- Department of Chemistry
- Universidade Federal de Minas Gerais
- Belo Horizonte
- Brazil
| | - Gisele L. Pappa
- Department of Computer Science
- Universidade Federal de Minas Gerais
- Belo Horizonte
- Brazil
| | - Jadson C. Belchior
- Department of Chemistry
- Universidade Federal de Minas Gerais
- Belo Horizonte
- Brazil
| |
Collapse
|
45
|
Brovarets’ OO, Hovorun DM. Novel physico-chemical mechanism of the mutagenic tautomerisation of the Watson–Crick-like A·G and C·T DNA base mispairs: a quantum-chemical picture. RSC Adv 2015. [DOI: 10.1039/c5ra11773a] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Novel routes for the mutagenic tautomerisation of the long A·G and short C·T Watson–Crick DNA base mispairs via sequential DPT are reported, pursuing the goal of an estimation of their contribution into spontaneous point replication errors in DNA.
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
| | - 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
|
46
|
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]
|
47
|
Sassi M, Carter DJ, Uberuaga BP, Stanek CR, Mancera RL, Marks NA. Hydrogen Bond Disruption in DNA Base Pairs from 14C Transmutation. J Phys Chem B 2014; 118:10430-5. [DOI: 10.1021/jp508118f] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Michel Sassi
- Nanochemistry Research Institute & Department of Chemistry, Curtin University, GPO Box U1987, Perth, WA 6845, Australia
- Fundamental
and Computational Sciences Directorate, Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Damien J. Carter
- Nanochemistry Research Institute & Department of Chemistry, Curtin University, GPO Box U1987, Perth, WA 6845, Australia
| | - Blas P. Uberuaga
- Material
Science and Technology Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Christopher R. Stanek
- Material
Science and Technology Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Ricardo L. Mancera
- School
of Biomedical Sciences, CHIRI Biosciences, Curtin University, GPO Box U1987, Perth, WA 6845, Australia
| | - Nigel A. Marks
- Discipline of Physics & Astronomy, Curtin University, GPO Box U1987, Perth, WA 6845, Australia
| |
Collapse
|
48
|
Brovarets’ OO, Hovorun DM. Does the G·G*synDNA mismatch containing canonical and rare tautomers of the guanine tautomerise through the DPT? A QM/QTAIM microstructural study. Mol Phys 2014. [DOI: 10.1080/00268976.2014.927079] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
49
|
Brovarets' OO, Hovorun DM. The nature of the transition mismatches with Watson-Crick architecture: the G*·T or G·T* DNA base mispair or both? A QM/QTAIM perspective for the biological problem. J Biomol Struct Dyn 2014; 33:925-45. [PMID: 24842163 DOI: 10.1080/07391102.2014.924879] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
This study provides the first accurate investigation of the tautomerization of the biologically important guanine*·thymine (G*·T) DNA base mispair with Watson-Crick geometry, involving the enol mutagenic tautomer of the G and the keto tautomer of the T, into the G·T* mispair (∆G = .99 kcal mol(-1), population = 15.8% obtained at the MP2 level of quantum-mechanical theory in the continuum with ε = 4), formed by the keto tautomer of the G and the enol mutagenic tautomer of the T base, using DFT and MP2 methods in vacuum and in the weakly polar medium (ε = 4), characteristic for the hydrophobic interfaces of specific protein-nucleic acid interactions. We were first able to show that the G*·T↔G·T* tautomerization occurs through the asynchronous concerted double proton transfer along two antiparallel O6H···O4 and N1···HN3 H-bonds and is assisted by the third N2H···O2 H-bond, that exists along the entire reaction pathway. The obtained results indicate that the G·T* base mispair is stable from the thermodynamic point of view complex, while it is dynamically unstable structure in vacuum and dynamically stable structure in the continuum with ε = 4 with lifetime of 6.4·10(-12) s, that, on the one side, makes it possible to develop all six low-frequency intermolecular vibrations, but, on the other side, it is by three orders less than the time (several ns) required for the replication machinery to forcibly dissociate a base pair into the monomers during DNA replication. One of the more significant findings to emerge from this study is that the short-lived G·T* base mispair, which electronic interaction energy between the bases (-23.76 kcal mol(-1)) exceeds the analogical value for the G·C Watson-Crick nucleobase pair (-20.38 kcal mol(-1)), "escapes from the hands" of the DNA replication machinery by fast transforming into the G*·T mismatch playing an indirect role of its supplier during the DNA replication. So, exactly the G*·T mismatch was established to play the crucial role in the spontaneous point mutagenesis.
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., Kyiv 03680 , Ukraine
| | | |
Collapse
|
50
|
Allehyani BH, Elroby SA, Aziz SG, Hilal RH. Protonation and deprotonation enthalpies of alloxan and implications for the structure and energy of its complexes with water: a computational study. J Biomol Struct Dyn 2014; 33:897-910. [DOI: 10.1080/07391102.2014.919238] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Basmah H. Allehyani
- Faculty of Science, Chemistry Department, King Abdullaziz University, Jeddah, Saudi Arabia
| | - Shaaban A. Elroby
- Faculty of Science, Chemistry Department, King Abdullaziz University, Jeddah, Saudi Arabia
- Faculty of Science, Chemistry Department, Benisuief University, Benisuief, Egypt
| | - Saadullah G. Aziz
- Faculty of Science, Chemistry Department, King Abdullaziz University, Jeddah, Saudi Arabia
| | - Rifaat H. Hilal
- Faculty of Science, Chemistry Department, King Abdullaziz University, Jeddah, Saudi Arabia
- Faculty of Science, Chemistry Department, Cairo University, Cairo, Egypt
| |
Collapse
|