1
|
Pandey J, Prajapati P, Tandon P, Sinha K, Ayala AP, Ellena J. Vibrational and conformational analysis of structural phase transition in Estradiol 17β valerate with temperature. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 263:120219. [PMID: 34332239 DOI: 10.1016/j.saa.2021.120219] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2021] [Revised: 07/03/2021] [Accepted: 07/19/2021] [Indexed: 06/13/2023]
Abstract
Estradiol 17β valerate (E2V) is a hormonal medicine widely used in hormone replacement therapy. E2V undergoes a reversible isosymmetric structural phase transition at low temperature (̴ 250 K) which results from the reorientation of the valerate chain. The reversible isosymmetric structural phase transition follows Ehrenfest's classification when described as first-order and Buerger's classification when classified as order-disorder. The conformational difference also induces changes in molecular torsional angles and on the hydrogen bond pattern. In combination with density functional theory (DFT) calculations, vibrational spectroscopy has been used to correlate the valerate chain modes with the modifications of the dihedral angles on phase transition. We are expecting improvement in our understanding of the phase transition mechanism driven by the temperature. The Conformational analysis reveals the feasible structures corresponding to changes in the dihedral angles associated with the valerate chain. The infrared spectra of calculated conformers are in good agreement with the experimental spectra of E2V structure recorded at room temperature revealing that the changes in valerate chain modes at 1115 cm-1, 1200 cm-1and 1415 cm-1 fingerprint the molecular conformation. An investigation made to determine the ligand-protein interaction of E2V through docking against estrogen receptor (ER) reveals the inhibitive and agonist nature of E2V.
Collapse
Affiliation(s)
- Jaya Pandey
- Department of Physics, University of Lucknow, Lucknow 226007, India
| | - Preeti Prajapati
- Department of Physics, University of Lucknow, Lucknow 226007, India
| | - Poonam Tandon
- Department of Physics, University of Lucknow, Lucknow 226007, India.
| | - Kirti Sinha
- Department of Physics, University of Lucknow, Lucknow 226007, India
| | - Alejandro Pedro Ayala
- Departamento de Fìsica, Universidade Federal do Ceará, C.P. 6030, 60.455-900 Fortaleza, CE Brazil
| | - Javier Ellena
- Instituto de Fìsica de São Carlos, Universidade de São Paulo,C.P 369, 13560-970,São Carlos, S P, Brazil
| |
Collapse
|
2
|
Brovarets' OO, Muradova A, Hovorun DM. A Quantum-Mechanical Looking Behind the Scene of the Classic G·C Nucleobase Pairs Tautomerization. Front Chem 2020; 8:574454. [PMID: 33330362 PMCID: PMC7732530 DOI: 10.3389/fchem.2020.574454] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Accepted: 09/02/2020] [Indexed: 11/13/2022] Open
Abstract
For the first time, at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory, a comprehensive quantum-mechanical investigation of the physico-chemical mechanism of the tautomeric wobblization of the four biologically-important G·C nucleobase pairs by the participation of the monomers in rare, in particular mutagenic, tautomeric forms (marked with an asterisk) was provided. These novel tautomeric transformations (wobblization or shifting of the bases within the pair) are intrinsically inherent properties of the G·C nucleobase pairs. In this study, we have obtained intriguing results, lying far beyond the existing representations. Thus, it was shown that Löwdin's G*·C*(WC) base pair does not tautomerize according to the wobblization mechanism. Tautomeric wobblization of the G*·C*(rWC) (relative Gibbs free energy ΔG = 0.00/relative electronic energy ΔE = 0.00 kcal·mol-1) ("r"-means the configuration of the base pair in reverse position; "WC"-the classic Watson-Crick configuration) and G*t·C*(H) (ΔG = -0.19/ΔE = 0.29 kcal·mol-1) ("H"-Hoogsteen configuration;"t" denotes the O6H hydroxyl group in the trans position) base pairs are preceded by the stages of the base pairs tautomerization by the single proton transfer (SPT). It was established that the G*t·C*(rH) (ΔG = 2.21/ΔE = 2.81 kcal·mol-1) base pair can be wobbled through two different pathways via the traditional one-stage mechanism through the TSs, which are tight G+·C- ion pairs, stabilized by the participation of only two intermolecular H-bonds. It was found out that the G·C base pair is most likely incorporated into the DNA/RNA double helix with parallel strands in the G*·C*(rWC), G·C*(rwwc), and G*·C(rwwc) ("w"-wobble configuration of the pair) tautomeric forms, which are in rapid tautomeric equilibrium with each other. It was proven that the G*·C*(rWC) nucleobase pair is also in rapid tautomeric equilibrium with the eight tautomeric forms of the so-called Levitt base pair. It was revealed that a few cases of tautomerization via the DPT of the nucleobase pairs by the participation of the C8H group of the guanine had occurred. The biological role of the obtained results was also made apparent.
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
| | - 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
|
Energy of the CH⋯O H-bonds and others specific contacts in the quercetin molecule: QM/QTAIM approximation formulas. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.113456] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
4
|
A Never-Ending Conformational Story of the Quercetin Molecule: Quantum-Mechanical Investigation of the O3′H and O4′H Hydroxyl Groups Rotations. APPLIED SCIENCES-BASEL 2020. [DOI: 10.3390/app10031147] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The quercetin molecule is known to be an effective pharmaceutical compound of a plant origin. Its chemical structure represents two aromatic A and B rings linked through the C ring containing oxygen and five OH hydroxyl groups attached to the 3, 3′, 4′, 5, and 7 positions. In this study, a novel conformational mobility of the quercetin molecule was explored due to the turnings of the O3′H and O4′H hydroxyl groups, belonging to the B ring, around the exocyclic C-O bonds. It was established that the presence of only three degrees of freedom of the conformational mobility of the O3′H and O4′H hydroxyl groups is connected with their concerted behavior, which is controlled by the non-planar (in the case of the interconverting planar conformers) or locally non-planar (in other cases) TSsO3′H/O4′H transition states, in which O3′H and O4′H hydroxyl groups are oriented by the hydrogen atoms towards each other. We also explored the number of the physico-chemical and electron-topological characteristics of all intramolecular-specific contacts—hydrogen bonds and attractive van der Waals contacts at the conformers and also at the transition states. Long-terms perspectives for the investigations of the structural bases of the biological activity of this legendary molecule have been shortly described.
Collapse
|
5
|
A Hidden Side of the Conformational Mobility of the Quercetin Molecule Caused by the Rotations of the O3H, O5H and O7H Hydroxyl Groups: In Silico Scrupulous Study. Symmetry (Basel) 2020. [DOI: 10.3390/sym12020230] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
In this study at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of quantum-mechanical theory it was explored conformational variety of the isolated quercetin molecule due to the mirror-symmetrical hindered turnings of the O3H, O5H and O7H hydroxyl groups, belonging to the A and C rings, around the exocyclic C–O bonds. These dipole active conformational transformations proceed through the 72 transition states (TSs; C1 point symmetry) with non-orthogonal orientation of the hydroxyl groups relatively the plane of the A or C rings of the molecule (HO7C7C8/HO7C7C6 = ±(89.9–93.3), HO5C5C10 = ±(108.9–114.4) and HO3C3C4 = ±(113.6–118.8 degrees) (here and below signs ‘±’ corresponds to the enantiomers)) with Gibbs free energy barrier of activation ΔΔGTS in the range 3.51–16.17 kcal·mol−1 under the standard conditions (T = 298.1 K and pressure 1 atm): ΔΔGTSO7H (3.51–4.27) < ΔΔGTSO3H (9.04–11.26) < ΔΔGTSO5H (12.34–16.17 kcal mol−1). Conformational dynamics of the O3H and O5H groups is partially controlled by the intramolecular specific interactions O3H…O4, C2′/C6′H…O3, O3H…C2′/C6′, O5H…O4 and O4…O5, which are flexible and cooperative. Dipole-active interconversions of the enantiomers of the non-planar conformers of the quercetin molecule (C1 point symmetry) is realized via the 24 TSs with C1 point symmetry (HO3C3C2C1 = ±(11.0–19.1), HC2′/C6′C1′C2 = ±(0.6–2.9) and C3C2C1′C2′/C3C2C1′C6′ = ±(1.7–9.1) degree; ΔΔGTS = 1.65–5.59 kcal·mol−1), which are stabilized by the participation of the intramolecular C2′/C6′H…O1 and O3H…HC2′/C6′ H-bonds. Investigated conformational rearrangements are rather quick processes, since the time, which is necessary to acquire thermal equilibrium does not exceed 6.5 ns.
Collapse
|
6
|
Brovarets’ OO, Hovorun DM. A new era of the prototropic tautomerism of the quercetin molecule: A QM/QTAIM computational advances. J Biomol Struct Dyn 2019; 38:4774-4800. [PMID: 31711364 DOI: 10.1080/07391102.2019.1691660] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/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
| | - 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. Intramolecular tautomerization of the quercetin molecule due to the proton transfer: QM computational study. PLoS One 2019; 14:e0224762. [PMID: 31751372 PMCID: PMC6874073 DOI: 10.1371/journal.pone.0224762] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2019] [Accepted: 10/20/2019] [Indexed: 12/21/2022] Open
Abstract
Quercetin molecule (3, 3', 4', 5, 7-pentahydroxyflavone, C15H10O7) is an important flavonoid compound of natural origin, consisting of two aromatic A and B rings linked through the C ring with endocyclic oxygen atom and five hydroxyl groups attached to the 3, 3', 4', 5 and 7 positions. This molecule is found in many foods and plants, and is known to have a wide range of therapeutic properties, like an anti-oxidant, anti-toxic, anti-inflammatory etc. In this study for the first time we have revealed and investigated the pathways of the tautomeric transformations for the most stable conformers of the isolated quercetin molecule (Brovarets' & Hovorun, 2019) via the intramolecular proton transfer. Energetic, structural, dynamical and polar characteristics of these transitions, in particular relative Gibbs free and electronic energies, characteristics of the intramolecular specific interactions-H-bonds and attractive van der Waals contacts, have been analysed in details. It was demonstrated that the most probable process among all investigated is the proton transfer from the O3H hydroxyl group of the C ring to the C2' carbon atom of the C2'H group of the B ring along the intramolecular O3H…C2' H-bond with the further formation of the C2'H2 group. It was established that the proton transfer from the hydroxyl groups to the carbon atoms of the neighboring CH groups is assisted at the transition states by the strong intramolecular HCH…O H-bond (~28.5 kcal∙mol-1). The least probable path of the proton transfer-from the C8H group to the endocyclic O1 oxygen atom-causes the decyclization of the C ring in some cases. It is shortly discussed the biological importance of the obtained results.
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
| | - 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, 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
|
9
|
Brovarets’ OO, Hovorun DM. Conformational diversity of the quercetin molecule: a quantum-chemical view. J Biomol Struct Dyn 2019; 38:2817-2836. [DOI: 10.1080/07391102.2019.1656671] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/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
| | - 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
|
Brovarets’ OO, Hovorun DM. Conformational transitions of the quercetin molecule via the rotations of its rings: a comprehensive theoretical study. J Biomol Struct Dyn 2019; 38:2865-2883. [DOI: 10.1080/07391102.2019.1645734] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/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 Pharmacology, Bohomolets 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 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
|
11
|
Alcolea Palafox M, Franklin Benial AM, K Rastogi V. Biomolecules of 2-Thiouracil, 4-Thiouracil and 2,4-Dithiouracil: A DFT Study of the Hydration, Molecular Docking and Effect in DNA:RNAMicrohelixes. Int J Mol Sci 2019; 20:E3477. [PMID: 31311161 PMCID: PMC6678171 DOI: 10.3390/ijms20143477] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2019] [Revised: 07/10/2019] [Accepted: 07/11/2019] [Indexed: 12/17/2022] Open
Abstract
The molecular structure of 2-thiouracil, 4-thiouracil and 2,4-dithiouracil was analyzed under the effect of the first and second hydration shell by using the B3LYP density functional (DFT) method, and the results were compared to those obtained for the uracil molecule. A slight difference in the water distribution appears in these molecules. On the hydration of these molecules several trends in bond lengths and atomic charges were established. The ring in uracil molecule appears easier to be deformed and adapted to different environments as compared to that when it is thio-substituted. Molecular docking calculations of 2-thiouracil against three different pathogens: Bacillus subtilis, Escherichia coli and Candida albicans were carried out. Docking calculations of 2,4-dithiouracil ligand with various targeted proteins were also performed. Different DNA: RNA hybrid microhelixes with uridine, 2-thiouridine, 4-thiouridine and 2,4-dithiouridine nucleosides were optimized in a simple model with three nucleotide base pairs. Two main types of microhelixes were analyzed in detail depending on the intramolecular H-bond of the 2'-OH group. The weaker Watson-Crick (WC) base pair formed with thio-substituted uracil than with unsubstituted ones slightly deforms the helical and backbone parameters, especially with 2,4-dithiouridine. However, the thio-substitution significantly increases the dipole moment of the A-type microhelixes, as well as the rise and propeller twist parameters.
Collapse
Affiliation(s)
- M Alcolea Palafox
- Departamento de Química-Física, Facultad de CienciasQuímicas, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| | | | - V K Rastogi
- Indian Spectroscopy Society, KC 68/1, Old Kavinagar, Ghaziabad 201002, India
| |
Collapse
|
12
|
Alcolea Palafox M. Effect of the sulfur atom on S2 and S4 positions of the uracil ring in different DNA:RNA hybrid microhelixes with three nucleotide base pairs. Biopolymers 2019; 110:e23247. [PMID: 30676643 DOI: 10.1002/bip.23247] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2018] [Revised: 11/03/2018] [Accepted: 11/07/2018] [Indexed: 02/06/2023]
Abstract
The effect of the sulphur atom on the uracil ring was analyzed in different DNA:RNA microhelixes with three nucleotide base-pairs, including uridine, 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, cytidine, adenosine and guanosine. Distinct backbone and helical parameters were optimized at different density functional (DFT) levels. The Watson-Crick pair with 2-thiouridine appears weaker than with uridine, but its interaction with water molecules appears easier. Two types of microhelixes were found, depending on the H-bond of H2' hydroxyl atom: A-type appears with the ribose ring in 3 E-envelope C3' -endo, and B-type in 2 E-envelope C2' -endo. B-type is less common but it is more stable and with higher dipole-moment. The sulphur atoms significantly increase the dipole-moment of the microhelix, as well as the rise and propeller twist parameters. Simulations with four Na atoms H-bonded to the phosphate groups, and further hydration with explicit water molecules were carried out. A re-definition of the numerical value calculation of several base-pair and base-stacking parameters is suggested.
Collapse
Affiliation(s)
- Mauricio Alcolea Palafox
- Departamento de Química-Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Madrid, Spain
| |
Collapse
|
13
|
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
|
14
|
Non-dissociative structural transitions of the Watson-Crick and reverse Watson-Crick А·Т DNA base pairs into the Hoogsteen and reverse Hoogsteen forms. Sci Rep 2018; 8:10371. [PMID: 29991693 PMCID: PMC6039495 DOI: 10.1038/s41598-018-28636-y] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2018] [Accepted: 06/21/2018] [Indexed: 11/08/2022] Open
Abstract
In this study it was theoretically shown that discovered by us recently (Brovarets' et al., Frontiers in Chemistry, 2018, 6:8; doi: 10.3389/fchem.2018.00008) high-energetical, significantly non-planar (symmetry C1), short-lived wobbled conformers of the classical Watson-Crick А·Т(WC), reverse Watson-Crick А·Т(rWC), Hoogsteen А·Т(Н) and reverse Hoogsteen А·Т(rН) DNA base pairs are the intermediates of their pairwise А∙Т(WC)/А∙Т(rWC) ↔ А∙Т(H)/А∙Т(rH) conformational transformations. These transitions do not require for their realization the energy-consumable anisotropic rotation of the amino group of A around the exocyclic C6-N6 bond. They are controlled by the non-planar transition states with quasi-orthogonal geometry (symmetry C1) joined by the single intermolecular (Т)N3H···N6(А) H-bond (~4 kcal∙mol-1). The Gibbs free energies of activation for these non-dissociative, dipole-active conformational transitions consist 7.33 and 7.81 kcal∙mol-1, accordingly. Quantum-mechanical (QM) calculations in combination with Bader's quantum theory of "Atoms in Molecules" (QTAIM) have been performed at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of QM theory in the continuum with ε = 4 under normal conditions.
Collapse
|
15
|
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
|
16
|
Alcolea Palafox M, Kattan D, Afseth N. FT-IR spectra of the anti-HIV nucleoside analogue d4T (Stavudine). Solid state simulation by DFT methods and scaling by different procedures. J Mol Struct 2018. [DOI: 10.1016/j.molstruc.2017.12.079] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
17
|
Khoee S, Hashemi A, Molavipordanjani S. Synthesis and characterization of IUdR loaded PEG/PCL/PEG polymersome in mixed DCM/DMF solvent: Experimental and molecular dynamics insights into the role of solvent composition and star architecture in drug dispersion and diffusion. Eur J Pharm Sci 2018; 114:1-12. [DOI: 10.1016/j.ejps.2017.11.019] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2017] [Revised: 11/19/2017] [Accepted: 11/20/2017] [Indexed: 12/19/2022]
|
18
|
Brovarets' OO, Tsiupa KS, Hovorun DM. Surprising Conformers of the Biologically Important A·T DNA Base Pairs: QM/QTAIM Proofs. Front Chem 2018; 6:8. [PMID: 29536003 PMCID: PMC5835050 DOI: 10.3389/fchem.2018.00008] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2017] [Accepted: 01/11/2018] [Indexed: 11/22/2022] Open
Abstract
For the first time novel high-energy conformers-A·T(wWC) (5.36), A·T(wrWC) (5.97), A·T(wH) (5.78), and A·T(wrH) (ΔG = 5.82 kcal·mol-1) (See Graphical Abstract) were revealed for each of the four biologically important A·T DNA base pairs - Watson-Crick A·T(WC), reverse Watson-Crick A·T(rWC), Hoogsteen A·T(H) and reverse Hoogsteen A·T(rH) at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of quantum-mechanical theory in the continuum with ε = 4 under normal conditions. Each of these conformers possesses substantially non-planar wobble (w) structure and is stabilized by the participation of the two anti-parallel N6H/N6H'…O4/O2 and N3H…N6 H-bonds, involving the pyramidalized amino group of the A DNA base as an acceptor and a donor of the H-bonding. The transition states - TSA·T(WC)↔A·T(wWC), TSA·T(rWC)↔A·T(wrWC), TSA·T(H)↔A·T(wH), and TSA·T(rH)↔A·T(wrH), controlling the dipole-active transformations of the conformers from the main plane-symmetric state into the high-energy, significantly non-planar state and vice versa, were localized. They also possess wobble structures similarly to the high-energy conformers and are stabilized by the participation of the N6H/N6H'…O4/O2 and N3H…N6 H-bonds. Discovered conformers of the A·T DNA base pairs are dynamically stable short-lived structures [lifetime τ = (1.4-3.9) ps]. Their possible biological significance and future perspectives have been briefly discussed.
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
|
19
|
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
|
20
|
Bednarko J, Stachurski O, Wielińska J, Kozakiewicz K, Liberek B, Nowacki A. Threocytidines: Insight into the Conformational Preferences of Artificial Threose Nucleic Acid (TNA) Building Blocks in B3LYP Studies. J Mol Graph Model 2018; 80:157-172. [PMID: 29366882 DOI: 10.1016/j.jmgm.2018.01.007] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Revised: 01/10/2018] [Accepted: 01/11/2018] [Indexed: 10/18/2022]
Abstract
A systematic DFT conformational studies of four building blocks of TNA with cytosine attached to the C1' atom of the α-L-threofuranose moiety are presented. Structures bearing 2'-OR and 3'-OR substituents, where R represents H, CH3 and phosphate groups, were used in the studies using a B3LYP functional in the gas phase. The χ angle (C2-N1-C1'-O4'), the ν0-ν4 endocyclic torsion angles and the exocyclic torsion angles ε (X-O2'-C2'-C1') and γ (X-O3'-C3'-C2') geometry parameter variations were taken into consideration. Three energy minima, high-anti, anti and syn, were found for the rotation about the C1'-N1 bond. The high-anti orientation of the base with respect to the sugar moiety, turned out to be preferred, regardless of the substituents at the C2' and C3' positions. Other orientations are at least 1.65 kcal/mol higher in Gibbs free energy than the high-anti one. It has been shown that intramolecular H-bonds and the anomeric effect of phosphate groups strongly affect the conformational preferences of the studied compounds. Further, the structure of substituents attached to the sugar moiety influence the pucker of the furanoid ring. The furanoid ring in the global minima of the compound with two OH groups (TC1) in the 2' and 3' positions, and the compound having a 3'-phosphate group (TC2), adopt roughly the same conformation located at the southern range of the pseudorotation wheel, and thus are close to those found in the B type DNA helix. The low-energy high-anti rotamers of the geometry with the phosphate group attached to the sugar ring in the 2' position (TC3) and the geometry with two methoxyl groups (TC4) have their furanoid rings in conformations resembling those found in A DNA and RNA helices (the northern range of the pseudorotation wheel).
Collapse
Affiliation(s)
- Justyna Bednarko
- Faculty of Chemistry, Department of Organic Chemistry, University of Gdańsk, Wita Stwosza 63, PL-80-308 Gdańsk, Poland
| | - Oktawian Stachurski
- Faculty of Chemistry, Department of Organic Chemistry, University of Gdańsk, Wita Stwosza 63, PL-80-308 Gdańsk, Poland
| | - Justyna Wielińska
- Faculty of Chemistry, Department of Organic Chemistry, University of Gdańsk, Wita Stwosza 63, PL-80-308 Gdańsk, Poland
| | - Karol Kozakiewicz
- Faculty of Chemistry, Department of Organic Chemistry, University of Gdańsk, Wita Stwosza 63, PL-80-308 Gdańsk, Poland
| | - Beata Liberek
- Faculty of Chemistry, Department of Organic Chemistry, University of Gdańsk, Wita Stwosza 63, PL-80-308 Gdańsk, Poland
| | - Andrzej Nowacki
- Faculty of Chemistry, Department of Organic Chemistry, University of Gdańsk, Wita Stwosza 63, PL-80-308 Gdańsk, Poland
| |
Collapse
|
21
|
Alcolea Palafox M, Rastogi VK, Singh SP. FT-IR and FT-Raman spectra of 5-chlorocytosine: Solid state simulation and tautomerism. Effect of the chlorine substitution in the Watson-Crick base pair 5-chlorodeoxycytidine-deoxyguanosine. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2018; 188:418-435. [PMID: 28756258 DOI: 10.1016/j.saa.2017.07.018] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/26/2016] [Revised: 04/26/2017] [Accepted: 07/13/2017] [Indexed: 06/07/2023]
Abstract
The laser Raman and IR spectra of 5-chlorocytosine have been recorded and accurately assigned in the solid state using Density functional calculations (DFT) together with the linear scaling equation procedure (LSE) and the solid state simulation of the crystal unit cell through a tetramer form. These results remarkably improve those reported previously by other authors. Several new scaling equations were proposed to be used in related molecules. The six main tautomers of the biomolecule 5-chlorocytosine were determined and optimized at the MP2 and CCSD levels, using different basis sets. The relative stabilities were compared with those obtained in cytosine and their 5-halo derivatives. Several relationships between energies, geometric parameters and NBO atomic charges were established. The effect of the chlorine substitution in the fifth position was evaluated through the stability of the Watson-Crick (WC) base pair of 5-chlorodeoxycytidine with deoxyguanosine, and through their vibrational spectra.
Collapse
Affiliation(s)
- M Alcolea Palafox
- Nofima AS - the Norwegian Institute of Food, Fisheries and Aquaculture Research, PB 210, N-1431 Ås, Norway; Departamento de Química-Física I, Facultad de Ciencias Químicas, Universidad Complutense, Ciudad Universitaria, Madrid 28040, Spain.
| | - V K Rastogi
- Internet Lab, R.D. Foundation Group of Institutions, NH-58, Kadrabad, Modinagar, Ghaziabad, India; Indian Spectroscopy Society, KC 68/1, Old Kavinagar, Ghaziabad 201 002, India.
| | - S P Singh
- Department of Physics, Dr B R Ambedkar Govt Degree College, Mainpuri 205 001, India
| |
Collapse
|
22
|
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
|
23
|
Brovarets' OO, Voiteshenko IS, Pérez-Sánchez H, Hovorun DM. A QM/QTAIM detailed look at the Watson-Crick↔wobble tautomeric transformations of the 2-aminopurine·pyrimidine mispairs. J Biomol Struct Dyn 2017; 36:1649-1665. [PMID: 28514900 DOI: 10.1080/07391102.2017.1331864] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
This work is devoted to the careful QM/QTAIM analysis of the evolution of the basic physico-chemical parameters along the intrinsic reaction coordinate (IRC) of the biologically important 2AP·T(WC)↔2AP·T*(w) and 2AP·C*(WC)↔2AP·C(w) Watson-Crick(WC)↔wobble(w) tautomeric transformations obtained at each point of the IRC using original authors' methodology. Established profiles reflect the high similarity between the courses of these processes. Basing on the scrupulous analysis of the profiles of their geometric and electron-topological parameters, it was established that the dipole-active WC↔w tautomerizations of the Watson-Crick-like 2AP·T(WC)/2AP·C*(WC) mispairs, stabilized by the two classical N3H⋯N1, N2H⋯O2 and one weak C6H⋯O4/N4 H-bonds, into the wobble 2AP·T*(w)/2AP·C(w) base pairs, respectively, joined by the two classical N2H⋯N3 and O4/N4H⋯N1 H-bonds, proceed via the concerted stepwise mechanism through the sequential intrapair proton transfer and subsequent large-scale shifting of the bases relative each other, through the planar, highly stable, zwitterionic transition states stabilized by the participation of the four H-bonds - N1+H⋯O4-/N4-, N1+H⋯N3-, N2+H⋯N3-, and N2+H⋯O2-. Moreover, it was found out that the 2AP·T(WC)↔2AP·T*(w)/2AP·C*(WC)↔2AP·C(w) tautomerization reactions occur non-dissociatively and are accompanied by the consequent replacement of the 10 unique patterns of the specific intermolecular interactions along the IRC. Obtained data are of paramount importance in view of their possible application for the control and management of the proton transfer, e.g. by external electric or laser fields.
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
| | - Ivan S Voiteshenko
- 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
| | - Horacio Pérez-Sánchez
- c Computer Science Department , Bioinformatics and High Performance Computing (BIO-HPC) Research Group, Universidad Católica San Antonio de Murcia (UCAM) , Guadalupe, Murcia 30107 , Spain
| | - 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
|
24
|
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
|
25
|
A DFT analysis of the molecular structure, vibrational spectra and other molecular properties of 5-nitrouracil and comparison with uracil. J Mol Struct 2016. [DOI: 10.1016/j.molstruc.2015.10.096] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
26
|
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
|
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
|
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]
|
29
|
Palafox MA. Structure and conformational analysis of the anti-HIV reverse transcriptase inhibitor AZT using MP2 and DFT methods. Differences with the natural nucleoside thymidine. Simulation of the 1st phosphorylation step with ATP. Phys Chem Chem Phys 2015; 16:24763-83. [PMID: 25317835 DOI: 10.1039/c4cp03695f] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
A comprehensive quantum-chemical investigation of the conformational landscape of the HIV-1 reverse transcriptase inhibitor AZT (3'-azido-3'-deoxythymidine) nucleoside analogue was carried out. The whole conformational parameters (χ, γ, β, δ, ϕ, P, νmax) were analysed as well as the NBO charges. The search located at least 55 stable structures, 9 of which were by MP2 within a 1 kcal mol(-1) electronic energy range of the global minimum. Most conformers were anti or high-anti around the glycoside bond and with North sugar ring puckering angles. The distribution of all the conformers according to the ranges of stability of the characteristic torsional angles was established. The results obtained were in accordance with those found in related anti-HIV nucleoside analogues. The best conformer in the anti form corresponded to the calculated values by MP2 of χ = -126.9°, β = 176.4° and γ = 49.1°. An analysis of the lowest vibrations in conformer C1 was carried out. The first hydration shell was simulated and the structural differences with the natural nucleoside deoxythymidine (dT) were determined. The first phosphorylation step was simulated by interacting ATP with the best hydrated clusters of AZT and dT. The Na cations act as a bridge between the phosphate moieties of ATP making it easy for -P3O3 to receive the H5' proton from AZT or dT. A proton-transfer mechanism is proposed through the water molecules. When the number of the water molecules surrounding AZT is lower than 8, the first phosphorylation step of AZT can be carried out. However, the appropriate orientation of the O5'-H in dT avoids this limitation and it can be performed with large numbers of water molecules.
Collapse
Affiliation(s)
- M Alcolea Palafox
- Chemical Physics Department, Chemistry Faculty, Complutense University, Ciudad Universitaria, Madrid-28040, Spain.
| |
Collapse
|
30
|
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
|
31
|
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
|
32
|
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]
|
33
|
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
|
34
|
Alvarez-Ros MC, Palafox MA. Conformational analysis, molecular structure and solid state simulation of the antiviral drug acyclovir (zovirax) using density functional theory methods. Pharmaceuticals (Basel) 2014; 7:695-722. [PMID: 24915059 PMCID: PMC4078516 DOI: 10.3390/ph7060695] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2014] [Revised: 05/22/2014] [Accepted: 05/26/2014] [Indexed: 11/23/2022] Open
Abstract
The five tautomers of the drug acyclovir (ACV) were determined and optimised at the MP2 and B3LYP quantum chemical levels of theory. The stability of the tautomers was correlated with different parameters. On the most stable tautomer N1 was carried out a comprehensive conformational analysis, and the whole conformational parameters (R, β, Φ, φ1, φ2, φ3, φ4, φ5) were studied as well as the NBO Natural atomic charges. The calculations were carried out with full relaxation of all geometrical parameters. The search located at least 78 stable structures within 8.5 kcal/mol electronic energy range of the global minimum, and classified in two groups according to the positive or negative value of the torsional angle j1. In the nitrogen atoms and in the O2' and O5' oxygen atoms of the most stable conformer appear a higher reactivity than in the natural nucleoside deoxyguanosine. The solid state was simulated through a dimer and tetramer forms and the structural parameters were compared with the X-ray crystal data available. Several general conclusions were emphasized.
Collapse
Affiliation(s)
- Margarita Clara Alvarez-Ros
- Departamento de Química-Física I. Facultad de Ciencias Químicas, Universidad Complutense, Ciudad Universitaria, Madrid 28040, Spain.
| | - Mauricio Alcolea Palafox
- Departamento de Química-Física I. Facultad de Ciencias Químicas, Universidad Complutense, Ciudad Universitaria, Madrid 28040, Spain.
| |
Collapse
|
35
|
El-Sayed AA, Tamara Molina A, Álvarez-Ros M, Alcolea Palafox M. Conformational analysis of the anti-HIV Nikavir prodrug: comparisons with AZT and Thymidine, and establishment of structure–activity relationships/tendencies in other 6′-derivatives. J Biomol Struct Dyn 2014; 33:723-48. [DOI: 10.1080/07391102.2014.909743] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
| | - A. Tamara Molina
- Facultad de Ciencias Químicas, Departamento de Química-Física I, Universidad Complutense, Ciudad Universitaria, Madrid 28040, Spain
| | - M.C. Álvarez-Ros
- Facultad de Ciencias Químicas, Departamento de Química-Física I, Universidad Complutense, Ciudad Universitaria, Madrid 28040, Spain
| | - M. Alcolea Palafox
- Facultad de Ciencias Químicas, Departamento de Química-Física I, Universidad Complutense, Ciudad Universitaria, Madrid 28040, Spain
| |
Collapse
|