1
|
Athmani H, Direm A, Osório FAP, Valverde C. Exploring the nonlinear optical properties of hypoxanthinium salts: a structural and computational analysis. J Mol Model 2024; 30:280. [PMID: 39046583 DOI: 10.1007/s00894-024-06077-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2024] [Accepted: 07/13/2024] [Indexed: 07/25/2024]
Abstract
CONTEXT In this study, we detail the synthesis and crystallographic characterization of an unprecedented structure, specifically hypoxanthinium chloride monohydrate ((I) hereafter), which crystallizes in the monoclinic P21/c space group. A comparative analysis was conducted with four related hypoxanthinium salts: hypoxanthinium bromide monohydrate (II), 9-methylhypoxanthinium chloride monohydrate (III), hypoxanthinium nitrate monohydrate (IV), and hypoxanthinium perchlorate monohydrate (V). This analysis has focused mainly on their crystal packing, hydrogen-bonding networks, and non-classical intermolecular interactions, as elucidated by comprehensive Hirshfeld surface and topological analyses. Theoretical investigation of the nonlinear optical (NLO) properties of the hypoxanthinium derivatives (I-V) was performed using the Density Functional Theory (DFT). METHODS The crystalline environment was simulated using the iterative Supermolecule method (SM), and the static and dynamics linear refractive index, linear polarizability, second-order hyperpolarizability, and the third-order nonlinear susceptibility at the DFT/CAM-B3LYP/6-311++G(d,p) level were computed. The results for the macroscopic third-order nonlinear susceptibility of (II) was found to equalχ 3 = 0.81 × 10 - 20 m 2 / V 2 . By replacing the bromine atom in (II) with a chlorine atom as in (III), theχ 3 value will be multiplied by 2.16, and therefore these results are large enough to suggest the potential application of these crystals as NLO materials.
Collapse
Affiliation(s)
- Hamza Athmani
- Laboratory of Structures, Properties and Interatomic Interactions LASPI2A, Faculty of Sciences and Technology, Abbes Laghrour University, 40.000, Khenchela, Algeria
- Department of Matter Sciences, Faculty of Sciences and Technology, Abbes Laghrour University, 40.000, Khenchela, Algeria
| | - Amani Direm
- Laboratory of Structures, Properties and Interatomic Interactions LASPI2A, Faculty of Sciences and Technology, Abbes Laghrour University, 40.000, Khenchela, Algeria.
- Department of Matter Sciences, Faculty of Sciences and Technology, Abbes Laghrour University, 40.000, Khenchela, Algeria.
| | - Francisco A P Osório
- Instituto de Física, Universidade Federal de Goiás, 74.690-900, Goiânia, GO, Brazil
| | - Clodoaldo Valverde
- Laboratório de Modelagem Molecular Aplicada E Simulação (LaMMAS), Universidade Estadual de Goiás, Campus de CiênciasExatas E Tecnológicas, Anápolis, GO, 75001-970, Brazil.
- Universidade Paulista, Goiânia, GO, 74845-090, Brazil.
| |
Collapse
|
2
|
Alcolea Palafox M, Belskaya NP, Todorov LT, Hristova-Avakoumova NG, Kostova IP. Effect of Lanthanide Ions and Triazole Ligands on the Molecular Properties, Spectroscopy and Pharmacological Activity. Int J Mol Sci 2024; 25:7964. [PMID: 39063204 PMCID: PMC11276792 DOI: 10.3390/ijms25147964] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2024] [Revised: 07/12/2024] [Accepted: 07/17/2024] [Indexed: 07/28/2024] Open
Abstract
The effect of La, Ce, Pr and Nd ions on four Ln(ligand)3 complexes and at three DFT levels of calculation was analyzed. Four ligands were chosen, three of which were based on the 1,2,3-triazole ring. The DFT methods used were B3LYP, CAM-B3LYP and M06-2X. The relationships established were between the geometric parameters, atomic charges, HOMO-LUMO energies and other molecular properties. These comparisons and trends will facilitate the synthesis of new complexes by selecting the ligand and lanthanide ion best suited to the desired property of the complex. The experimental IR and Raman spectra of Ln(2b')3 complexes where Ln = La, Ce, Pr, Nd, Sm, Gd, Dy, Ho and Er ions have been recorded and compared to know the effect of the lanthanide ion on the complex. The hydration in these complexes was also analyzed. Additionally, the effect of the type of coordination center on the ability of an Ln(ligand)3 complex to participate in electron exchange and hydrogen transfer was investigated using two in vitro model systems-DPPH and ABTS.
Collapse
Affiliation(s)
- Mauricio Alcolea Palafox
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense, 28040 Madrid, Spain
| | - Nataliya P. Belskaya
- Department of Technology for Organic Synthesis, Ural Federal University, 19 Mira Str., Yekaterinburg 620012, Russia;
| | - Lozan T. Todorov
- Department of Chemistry, Faculty of Pharmacy, Medical University—Sofia, 2 Dunav Str., 1000 Sofia, Bulgaria;
| | - Nadya G. Hristova-Avakoumova
- Department of Medical Physics and Biophysics, Faculty of Medicine, Medical University—Sofia, 2 Zdrave Str., 1431 Sofia, Bulgaria;
| | - Irena P. Kostova
- Department of Chemistry, Faculty of Pharmacy, Medical University—Sofia, 2 Dunav Str., 1000 Sofia, Bulgaria;
| |
Collapse
|
3
|
Szymańska M, Majerz I. Prototropy, Intramolecular Interactions, Electron Delocalization, and Physicochemical Properties of 1,8-dihydroxy-9-anthrone-DFT-D3 Study of Substituent Effects. MOLECULES (BASEL, SWITZERLAND) 2023; 28:molecules28010344. [PMID: 36615539 PMCID: PMC9822037 DOI: 10.3390/molecules28010344] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 12/27/2022] [Accepted: 12/28/2022] [Indexed: 01/04/2023]
Abstract
1,8-dihydroxy-9-anthrone are tricyclic compounds with a ketone group in the middle ring and two hydroxyl groups substituted in the side-aromatic rings what results in formation of two intramolecular hydrogen bonds in which the oxygen atom from the ketone group is the proton acceptor. 1,8-dihydroxy-9-anthrones in which intramolecular proton transfer between C10 and CO in the middle ring occurs, can exist in a tautomeric keto-enol equilibrium. For anthralin, the most important representative of this group, this equilibrium has been studied previously, but it has not been studied for its derivatives. Substituents in the middle ring change the geometry of 1,8-dihydroxy-9-anthrones so they are also expected to affect the keto-enol equilibrium. It is also important to study the effect of intramolecular hydrogen bonds on the structure of both tautomeric forms. It was found that the nature of the substituent in the middle ring could affect the antioxidant properties of the investigated compound.
Collapse
|
4
|
Palafox MA, Kattan D, de Pedraza Velasco ML, Isasi J, Posada-Moreno P, Rani K, Singh SP, Rastogi VK. Base pairs with 4-amino-3-nitrobenzonitrile: comparison with the natural WC pairs. Dimer and tetramer forms, Infrared and Raman spectra, and several proposed antiviral modified nucleosides. J Biomol Struct Dyn 2022:1-23. [PMID: 35583120 DOI: 10.1080/07391102.2022.2069864] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
Base pairs of 4-amino-3-nitrobenzonitrile (4A-3NBN) molecule with uracil, thymine and cytosine nucleobases were optimized and compared to natural Watson-Crick (WC) pairs. The slightly greater flexibility of the -NO2 group of 4A-3NBN than the N3-H group of the natural nucleobases together with a noticeable higher dipole moment of its pairs can facilitate disruption of the DNA/RNA helix formation. Several new mutagenic modified nucleosides with 4A-3NBN and 3-amino-2-nitrobenzonitrile (3A-2NBN) were proposed as antiviral prodrugs and their base pairs optimized. The special characteristics of these prodrugs appear appropriated for their clinical use. The counterpoise (CP) corrected interaction energies of the base pairs were calculated and compared to the natural ones. The M06-2X DFT method was used for this purpose. The molecular structure of 4A-3NBN was analyzed in detail and the crystal unit cell was simulated by a tetramer form and eight dimer forms. The performance of the B3LYP, X3LYP and M06-2X methods was tested on the vibrational wavenumbers in the monomer, dimer and tetramer forms of 4A-3NBN. The observed IR and Raman bands were assigned according to the optimum dimer II form determined by B3LYP and by the tetramer form calculated by M06-2X, which is the expected unit cell that forms the crystal net. The two best scaling procedures were used.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- M Alcolea Palafox
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense, Madrid, Spain
| | - D Kattan
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense, Madrid, Spain
| | - M L de Pedraza Velasco
- Dpto. de Enfermería, Facultad de Enfermería, Fisioterapia y Podología, UCM, Madrid, Spain
| | - J Isasi
- Dpto. de Química Inorgánica, Facultad de Ciencias Químicas, UCM, Madrid, Spain
| | - P Posada-Moreno
- Dpto. de Enfermería, Facultad de Enfermería, Fisioterapia y Podología, UCM, Madrid, Spain.,Instituto de Investigación Sanitaria del Hospital Clínico San Carlos (IdISSC), Madrid, Spain
| | - Kaushal Rani
- Department of Physics, Meerut College, Meerut, India
| | - S P Singh
- Department of Physics, Dr B R Ambedkar College, Mainpuri, India
| | - V K Rastogi
- Indian Spectroscopy Society, Ghaziabad, India
| |
Collapse
|
5
|
Dutta N, Deb I, Sarzynska J, Lahiri A. Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2022; 169-170:21-52. [PMID: 35065168 DOI: 10.1016/j.pbiomolbio.2022.01.001] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/11/2021] [Revised: 12/11/2021] [Accepted: 01/11/2022] [Indexed: 05/21/2023]
Abstract
Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.
Collapse
Affiliation(s)
- Nivedita Dutta
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, 92, Acharya Prafulla Chandra Road, Kolkata, 700009, West Bengal, India
| | - Indrajit Deb
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, 92, Acharya Prafulla Chandra Road, Kolkata, 700009, West Bengal, India
| | - Joanna Sarzynska
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61-704, Poznan, Poland
| | - Ansuman Lahiri
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, 92, Acharya Prafulla Chandra Road, Kolkata, 700009, West Bengal, India.
| |
Collapse
|
6
|
Brovarets’ OO, Muradova A, Hovorun DM. Novel horizons of the conformationally-tautomeric transformations of the G·T base pairs: quantum-mechanical investigation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2026510] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Ol’ha O. Brovarets’
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
| | - Alona Muradova
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| |
Collapse
|
7
|
Tolosa S, Sansón J, Hidalgo A. A procedure to understanding the C-G to A-T transversion. SMD simulations from guanine oxidation pathways assisted by one H2O2 molecule in the C-G basis pair. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.114123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
8
|
Tolosa S, Sansón J, Hidalgo A. Mechanisms of the T-A to C-G transition studied by SMD simulations: Deamination vs tautomerisation. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.113036] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
Marrero-Carballo R, Tun-Rosado F, Mena-Rejón GJ, Cáceres-Castillo D, Barroso J, Murillo F, Merino G, Quijano-Quiñones RF. The base-catalyzed keto-enol tautomerism of chrysophanol anthrone. A DFT investigation of the base-catalyzed reaction. MOLECULAR SIMULATION 2019. [DOI: 10.1080/08927022.2019.1581937] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Rubén Marrero-Carballo
- Laboratorio de Química Farmacéutica, Facultad de Química, Universidad Autónoma de Yucatán, Mérida, México
| | - Fernando Tun-Rosado
- Laboratorio de Química Teórica, Facultad de Química, Universidad Autónoma de Yucatán, Mérida, México
| | - Gonzalo J. Mena-Rejón
- Laboratorio de Química Farmacéutica, Facultad de Química, Universidad Autónoma de Yucatán, Mérida, México
| | - David Cáceres-Castillo
- Laboratorio de Química Farmacéutica, Facultad de Química, Universidad Autónoma de Yucatán, Mérida, México
| | - Jorge Barroso
- Departamento de Física Aplicada, Centro de Investigación y de Estudios Avanzados, Mérida, México
| | - Fernando Murillo
- Departamento de Física Aplicada, Centro de Investigación y de Estudios Avanzados, Mérida, México
| | - Gabriel Merino
- Departamento de Física Aplicada, Centro de Investigación y de Estudios Avanzados, Mérida, México
| | | |
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
|
Affiliation(s)
- Nadezhda Markova
- Institute of Organic Chemistry with Center of Phytochemistry, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria
| | - Venelin Enchev
- Institute of Organic Chemistry with Center of Phytochemistry, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria
| |
Collapse
|
15
|
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
|
16
|
Masoodi HR, Bagheri S, Ghaderi Z. The influence of Cu + binding to hypoxanthine on stabilization of mismatches involving hypoxanthine and DNA bases: a DFT study. J Biomol Struct Dyn 2018; 37:1923-1934. [PMID: 29757083 DOI: 10.1080/07391102.2018.1475256] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
In the present work, the influence of Cu+ binding to N3- and N7-positions of hypoxanthine on energetic, geometrical and topological properties of hypoxanthine-guanine, hypoxanthine-adenine, hypoxanthine-cytosine, hypoxanthine-thymine and hypoxanthine-hypoxanthine mismatches is theoretically investigated. The calculations, in gas phase, are performed at B3LYP/6-311++G(3df,3pd) level of theory. Unlike the other mispairs, Cu+ binding to N3-position of hypoxanthine causes the proton transfer process from enol form of hypoxanthine to imino forms of adenine and cytosine. This process also occurs in all mismatches having enol form of hypoxanthine when Cu+ binds to N7-position of hypoxanthine. The mismatches are stabilized by hydrogen bonds. The influence of Cu+ on hydrogen bonds is also examined by atoms in molecules (AIM) and natural bond orbital (NBO) analyses. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Hamid Reza Masoodi
- a Faculty of Science, Department of Chemistry , Vali-e-Asr University of Rafsanjan , Rafsanjan , Iran
| | - Sotoodeh Bagheri
- a Faculty of Science, Department of Chemistry , Vali-e-Asr University of Rafsanjan , Rafsanjan , Iran
| | - Zahra Ghaderi
- a Faculty of Science, Department of Chemistry , Vali-e-Asr University of Rafsanjan , Rafsanjan , Iran
| |
Collapse
|
17
|
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
|
18
|
Radhika R, Shankar R, Vijayakumar S, Kolandaivel P. Role of 6-Mercaptopurine in the potential therapeutic targets DNA base pairs and G-quadruplex DNA: insights from quantum chemical and molecular dynamics simulations. J Biomol Struct Dyn 2017; 36:1369-1401. [PMID: 28436311 DOI: 10.1080/07391102.2017.1323013] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The theoretical studies on DNA with the anticancer drug 6-Mercaptopurine (6-MP) are investigated using theoretical methods to shed light on drug designing. Among the DNA base pairs considered, 6-MP is stacked with GC with the highest interaction energy of -46.19 kcal/mol. Structural parameters revealed that structure of the DNA base pairs is deviated from the planarity of the equilibrium position due to the formation of hydrogen bonds and stacking interactions with 6-MP. These deviations are verified through the systematic comparison between X-H bond contraction and elongation and the associated blue shift and red shift values by both NBO analysis and vibrational analysis. Bent's rule is verified for the C-H bond contraction in the 6-MP interacted base pairs. The AIM results disclose that the higher values of electron density (ρ) and Laplacian of electron density (∇2ρ) indicate the increased overlap between the orbitals that represent the strong interaction and positive values of the total electron density show the closed-shell interaction. The relative sensitivity of the chemical shift values for the DNA base pairs with 6-MP is investigated to confirm the hydrogen bond strength. Molecular dynamics simulation studies of G-quadruplex DNA d(TGGGGT)4 with 6-MP revealed that the incorporation of 6-MP appears to cause local distortions and destabilize the G-quadruplex DNA.
Collapse
Affiliation(s)
- R Radhika
- a Department of Physics , Bharathiar University , Coimbatore , India
| | - R Shankar
- a Department of Physics , Bharathiar University , Coimbatore , India
| | - S Vijayakumar
- b Department of Medical Physics , Bharathiar University , Coimbatore , India
| | - P Kolandaivel
- a Department of Physics , Bharathiar University , Coimbatore , India
| |
Collapse
|
19
|
How OH and O– groups affect electronic structure of meta-substituted and para-substituted phenols and phenolates. Struct Chem 2017. [DOI: 10.1007/s11224-017-0965-4] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
20
|
Chiacchio MA, Legnani L, Caramella P, Tejero T, Merino P. Pivotal Neighboring-Group Participation in Substitution versus Elimination Reactions - Computational Evidence for Ion Pairs in the Thionation of Alcohols with Lawesson's Reagent. European J Org Chem 2017. [DOI: 10.1002/ejoc.201700127] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Maria A. Chiacchio
- Dipartimento di Chimica; Università di Pavia; Via Taramelli 12 27100 Pavia Italy
- Dipartimento di Scienze del Farmaco; Università di Catania; V.le A. Doria 6 95125 Catania Italy
| | - Laura Legnani
- Dipartimento di Chimica; Università di Pavia; Via Taramelli 12 27100 Pavia Italy
| | - Pierluigi Caramella
- Dipartimento di Chimica; Università di Pavia; Via Taramelli 12 27100 Pavia Italy
| | - Tomas Tejero
- Instituto de Biocomputación y Física de Sistemas Complejos (BIFI); Universidad de Zaragoza; 50009 Zaragoza Spain
| | - Pedro Merino
- Instituto de Biocomputación y Física de Sistemas Complejos (BIFI); Universidad de Zaragoza; 50009 Zaragoza Spain
| |
Collapse
|
21
|
Szymanski ES, Kimsey IJ, Al-Hashimi HM. Direct NMR Evidence that Transient Tautomeric and Anionic States in dG·dT Form Watson-Crick-like Base Pairs. J Am Chem Soc 2017; 139:4326-4329. [PMID: 28290687 PMCID: PMC5581979 DOI: 10.1021/jacs.7b01156] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The replicative and translational machinery utilizes the unique geometry of canonical G·C and A·T/U Watson-Crick base pairs to discriminate against DNA and RNA mismatches in order to ensure high fidelity replication, transcription, and translation. There is growing evidence that spontaneous errors occur when mismatches adopt a Watson-Crick-like geometry through tautomerization and/or ionization of the bases. Studies employing NMR relaxation dispersion recently showed that wobble dG·dT and rG·rU mismatches in DNA and RNA duplexes transiently form tautomeric and anionic species with probabilities (≈0.01-0.40%) that are in concordance with replicative and translational errors. Although computational studies indicate that these exceptionally short-lived and low-abundance species form Watson-Crick-like base pairs, their conformation could not be directly deduced from the experimental data, and alternative pairing geometries could not be ruled out. Here, we report direct NMR evidence that the transient tautomeric and anionic species form hydrogen-bonded Watson-Crick-like base pairs. A guanine-to-inosine substitution, which selectively knocks out a Watson-Crick-type (G)N2H2···O2(T) hydrogen bond, significantly destabilized the transient tautomeric and anionic species, as assessed by lack of any detectable chemical exchange by imino nitrogen rotating frame spin relaxation (R1ρ) experiments. An 15N R1ρ NMR experiment targeting the amino nitrogen of guanine (dG-N2) provides direct evidence for Watson-Crick (G)N2H2···O2(T) hydrogen bonding in the transient tautomeric state. The strategy presented in this work can be generally applied to examine hydrogen-bonding patterns in nucleic acid transient states including in other tautomeric and anionic species that are postulated to play roles in replication and translational errors.
Collapse
Affiliation(s)
- Eric S. Szymanski
- Department of Biochemistry, Duke University Medical Center, Durham, North Carolina 27710, United States
| | - Isaac J. Kimsey
- Department of Biochemistry, Duke University Medical Center, Durham, North Carolina 27710, United States
| | - Hashim M. Al-Hashimi
- Department of Biochemistry, Duke University Medical Center, Durham, North Carolina 27710, United States
- Department of Chemistry, Duke University, Durham, North Carolina 27710, United States
| |
Collapse
|
22
|
Ö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]
|
23
|
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
|
24
|
Dawley MM, Tanzer K, Carmichael I, Denifl S, Ptasińska S. Dissociative electron attachment to the gas-phase nucleobase hypoxanthine. J Chem Phys 2015; 142:215101. [DOI: 10.1063/1.4921388] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- M. Michele Dawley
- Radiation Laboratory, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Katrin Tanzer
- Institut für Ionenphysik und Angewandte Physik, and Center for Molecular Biosciences Innsbruck (CMBI), Leopold-Franzens Universität Innsbruck, Technikerstr. 25, A-6020 Innsbruck, Austria
| | - Ian Carmichael
- Radiation Laboratory, University of Notre Dame, Notre Dame, Indiana 46556, USA
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Stephan Denifl
- Institut für Ionenphysik und Angewandte Physik, and Center for Molecular Biosciences Innsbruck (CMBI), Leopold-Franzens Universität Innsbruck, Technikerstr. 25, A-6020 Innsbruck, Austria
| | - Sylwia Ptasińska
- Radiation Laboratory, University of Notre Dame, Notre Dame, Indiana 46556, USA
- Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556, USA
| |
Collapse
|
25
|
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
|
26
|
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
|
27
|
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
|
28
|
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
|
29
|
Allehyani BH, Elroby SA, Aziz SG, Hilal RH. Electronic structure of alloxan and its dimers: QM/QD simulations and quantum chemical topology analysis. J Biomol Struct Dyn 2015; 33:2121-32. [DOI: 10.1080/07391102.2014.997291] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- Basmah H. Allehyani
- Faculty of Science, Chemistry Department, King Abdul-Aziz University , Jeddah, Saudi Arabia
| | - Shaaban A. Elroby
- Faculty of Science, Chemistry Department, King Abdul-Aziz University , Jeddah, Saudi Arabia
- Faculty of Science, Chemistry Department, Beni Suef University , Beni Suef, Egypt
| | - Saadalluh G. Aziz
- Faculty of Science, Chemistry Department, King Abdul-Aziz University , Jeddah, Saudi Arabia
| | - Rifaat H. Hilal
- Faculty of Science, Chemistry Department, King Abdul-Aziz University , Jeddah, Saudi Arabia
- Faculty of Science, Chemistry Department, Cairo University , Cairo, Egypt
| |
Collapse
|
30
|
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
|
31
|
Raczyńska ED, Makowski M, Hallmann M, Kamińska B. Geometric and energetic consequences of prototropy for adenine and its structural models – a review. RSC Adv 2015. [DOI: 10.1039/c4ra17280a] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Prototropy for adenine and its convenient models causes parallel changes of geometric (HOMED) and energetic (ΔE) parameters for neutral tautomers.
Collapse
Affiliation(s)
- Ewa D. Raczyńska
- Department of Chemistry
- Warsaw University of Life Sciences (SGGW)
- 02-776 Warszawa
- Poland
| | | | - Małgorzata Hallmann
- Department of Chemistry
- Warsaw University of Life Sciences (SGGW)
- 02-776 Warszawa
- Poland
| | - Beata Kamińska
- Department of Chemistry
- Warsaw University of Life Sciences (SGGW)
- 02-776 Warszawa
- Poland
| |
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
|
Dawley MM, Tanzer K, Cantrell WA, Plattner P, Brinkmann NR, Scheier P, Denifl S, Ptasińska S. Electron ionization of the nucleobases adenine and hypoxanthine near the threshold: a combined experimental and theoretical study. Phys Chem Chem Phys 2014; 16:25039-53. [DOI: 10.1039/c4cp03452j] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
34
|
Ortiz S, Alvarez-Ros MC, Palafox MA, Rastogi VK, Balachandran V, Rathor SK. FT-IR and FT-Raman spectra of 6-chlorouracil: molecular structure, tautomerism and solid state simulation. A comparison between 5-chlorouracil and 6-chlorouracil. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2014; 130:653-668. [PMID: 24856263 DOI: 10.1016/j.saa.2014.04.009] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/04/2014] [Revised: 03/27/2014] [Accepted: 04/01/2014] [Indexed: 06/03/2023]
Abstract
A Raman and IR study of the biomolecule 6-chlorouracil was carried out in the solid state. The unit cell found in the crystal was simulated as a tetramer form by density functional calculations. Specific scale factors and scaling equations deduced from uracil molecule were employed in the predicted wavenumbers of 6-chlorouracil. The scaled wavenumbers were used in the reassignment of the IR and Raman experimental bands. Good reproduction of the experimental wavenumbers is obtained and the % error is very small in the majority of cases. A comparison between the molecular structure and charge distribution of 6-chlorouracil and 5-chlorouracil molecules was presented. The effect of the hydration with the PCM model in the molecular structure and charges was discussed. The optimum tautomers of 6-chlorouracil were optimized and analyzed. Six of them were related to those of uracil molecule. The effect of the halogen substitution in the sixth position of the pyrimidine ring in the stability of the different tautomers was evaluated. HOMO and LUMO orbital energy analysis were carried out.
Collapse
Affiliation(s)
- S Ortiz
- Departamento de Química-Física-I, Facultad de Ciencias Químicas, UCM, Madrid 28040, Spain
| | - M C Alvarez-Ros
- Departamento de Química-Física-I, Facultad de Ciencias Químicas, UCM, Madrid 28040, Spain
| | - M Alcolea Palafox
- Departamento de Química-Física-I, Facultad de Ciencias Químicas, UCM, Madrid 28040, Spain.
| | - V K Rastogi
- R.D. Foundation Group of Institutions, NH-58, Kadrabad, Modinagar, Ghaziabad, India; Indian Spectroscopy Society, KC-68/1, Old Kavinagar, Ghaziabad 201 002, India.
| | - V Balachandran
- Department of Physics, Arignar Anna Government Arts College, Musiri, Tiruchirappalli 621211, India
| | - S K Rathor
- Indian Spectroscopy Society, KC-68/1, Old Kavinagar, Ghaziabad 201 002, India
| |
Collapse
|
35
|
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]
|
36
|
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
|
37
|
Brovarets' OO, Hovorun DM. Atomistic nature of the DPT tautomerisation of the biologically important C·C* DNA base mispair containing amino and imino tautomers of cytosine: a QM and QTAIM approach. Phys Chem Chem Phys 2014; 15:20091-104. [PMID: 24154739 DOI: 10.1039/c3cp52644e] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
A theoretical study of tautomerisation of the biologically important cytosine·cytosine* (C·C*) DNA mismatch with a propeller-like structure (|C4N3N3C4| = 32.4°; C1 symmetry) and cis-oriented N1H glycosidic bonds, formed by the amino and imino tautomers of the C nucleobase, via the asynchronous concerted double proton transfer (DPT) along two H-bonds through the transition state (TSC·C*↔C*·C) (|C4N3N3C4| = 48.5°; C1 symmetry) into the C*·C mispair was carried out for the first time. It was established that the C·C*/C*·C DNA base mispair is associated by the antiparallel N4H···N4 (6.66 kcal mol(-1)), N3H···N3 (6.47 kcal mol(-1)) H-bonds and the O2···O2 van der Waals (vdW) contact (0.33 kcal mol(-1)), while the zwitterionic TSC·C*↔C*·C is stabilized by the parallel N4(+)H···N4(-) (13.55 kcal mol(-1)), N3(+)H···N3(-) (13.20 kcal mol(-1)) H-bonds and the O2(+)···O2(-) vdW contact (0.60 kcal mol(-1)). It was shown that the C·C* ↔ C*·C tautomerisation via the DPT is assisted by the O2···O2 vdW contact, that in contrast to the two others N4H···N4 and N3H···N3 H-bonds exists along the entire intrinsic reaction coordinate (IRC) range. The positive values of the Grunenberg's compliance constants (30.919 and 21.384 Å mdyn(-1) for C·C*/C*·C and TSC·C*↔C*·C, respectively) indicate that the O2···O2 vdW contact is a stabilizing closed-shell interaction. It was found that the middle N3H···N3 H-bond is anti-cooperative with the upper N4H···N4 H-bond and cooperative with the lower O2···O2 vdW contact. The 9 key points, which can be considered as electron-topological "fingerprints" of the asynchronous concerted C·C* ↔ C*·C tautomerisation process via the DPT were revealed along the IRC and examined in detail. It was shown that the C·C*/C*·C base mispair is a thermodynamically and dynamically stable structure. Its lifetime is equal to 1.53 × 10(-7) s at the MP2/cc-pVQZ//B3LYP/6-311++G(d,p) level of theory in vacuum. All 6 low-frequency intermolecular vibrations are able to develop during this time span.
Collapse
Affiliation(s)
- Ol'ha O Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, 150 Akademika Zabolotnoho Str., 03680 Kyiv, Ukraine.
| | | |
Collapse
|
38
|
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
|
39
|
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
|
40
|
Qiu ZM, Wang GL, Wang HL, Xi HP, Hou D. MP2 study on the hydrogen-bonding interaction between 5-fluorouracil and DNA bases: A,C,G,T. Struct Chem 2014. [DOI: 10.1007/s11224-014-0427-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
41
|
Brovarets' OO, Zhurakivsky RO, Hovorun DM. DPT tautomerisation of the wobble guanine·thymine DNA base mispair is not mutagenic: QM and QTAIM arguments. J Biomol Struct Dyn 2014; 33:674-89. [PMID: 24650179 DOI: 10.1080/07391102.2014.897259] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
We have shown for the first time, connecting QM methods with QTAIM analysis and using the methodology of the sweeps of the energetical, electron-topological and geometrical parameters, that the tautomerisation of the wobble guanine·thymine (wG·T) DNA base mispair into the wG(*)·T(*) base mispair induced by the double proton transfer (DPT), which undergoes a concerted asynchronous pathway, is not mutagenic. The wG·T → wG(*)·T(*) DPT tautomerisation does not result in the transition of the G base into its mutagenic tautomeric form G(*) able to mispair with the T base within the Watson-Crick base pairing scheme. This observation is explained by the so-called quantum protection of the wG·T DNA base mispair from its mutagenic tautomerisation - the dynamical non-stability of the tautomerised wG(*)·T(*) base mispair and significantly negative value of the Gibbs free energy of activation for the reverse reaction of the wG·T → wG(*)·T(*) DPT tautomerisation.
Collapse
Affiliation(s)
- Ol'ha O Brovarets'
- a Department of Molecular and Quantum Biophysics , Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine , 150 Akademika Zabolotnoho Str., 03680 Kyiv , Ukraine
| | | | | |
Collapse
|
42
|
Brovarets’ OO, Zhurakivsky RO, Hovorun DM. A QM/QTAIM microstructural analysis of the tautomerisationviathe DPT of the hypoxanthine·adenine nucleobase pair. Mol Phys 2014. [DOI: 10.1080/00268976.2013.877170] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
43
|
Brovarets' OO, Zhurakivsky RO, Hovorun DM. Does the tautomeric status of the adenine bases change upon the dissociation of the A*·Asyn Topal–Fresco DNA mismatch? A combined QM and QTAIM atomistic insight. Phys Chem Chem Phys 2014; 16:3715-25. [DOI: 10.1039/c3cp54708f] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
44
|
Brovarets’ OO, Zhurakivsky RO, Hovorun DM. Structural, energetic and tautomeric properties of the T·T∗/T∗·T DNA mismatch involving mutagenic tautomer of thymine: A QM and QTAIM insight. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2013.12.034] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
45
|
Brovarets' OO. [Structural and energetic properties of the four configurations of the A.T and G.C DNA base pairs]. UKRAINIAN BIOCHEMICAL JOURNAL 2013; 85:104-10. [PMID: 24319980 DOI: 10.15407/ubj85.04.104] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Using the methods of non-empirical quantum chemistry at the MP2/6-311++G(2df,pd)// B3LYP/6-311++G(d,p) level of theory it was established for the first time, that Hoogsteen, reverse Hoogsteen, Watson-Crick and reverse Watson-Crick configurations of the A.T and G.C DNA base pairs are isoelectronic and isomorphic structures with similar dynamic properties. Based on these results, non-ionisation mechanism of the Hoogsteen <"breathing" of the G*.C* DNA base pair, namely transformation of the tautomerised (Lowdin's) G-C base pair with Watson-Crick geometry into the Hoogsteen electroneutral G*.C* H base pair stabilized by the three O6H...N4, N3H...N7 and C8H...02 H-bonds, was postulated. It is suggested that such scenario activates only in those cases, when DNA is not located in aqueous solution, but works together with proteins and cytosine protonation at the N3 atom is precluded.
Collapse
|
46
|
Brovarets' OO. [Under what conditions does G.C Watson-Crick DNA base pair acquire all four configurations characteristic for A.T Watson-Crick DNA base pair?]. UKRAINIAN BIOCHEMICAL JOURNAL 2013; 85:98-103. [PMID: 24319979 DOI: 10.15407/ubj85.04.098] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
At the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory it was established for the first time, that the Löwdin's G*.C* DNA base pair formed by the mutagenic tautomers can acquire, as the A-T Watson-Crick DNA base pair, four biologically important configurations, namely: Watson-Crick, reverse Watson-Crick, Hoogsteen and reverse Hoogsteen. This fact demonstrates rather unexpected role of the tautomerisation of the one of the Watson-Crick DNA base pairs, in particular, via double proton transfer: exactly the G.C-->G*.C* tautomerisation allows to overcome steric hindrances for the implementation of the above mentioned configurations. Geometric, electron-topological and energetic properties of the H-bonds that stabilise the studied pairs, as well as the energetic characteristics of the latters are presented.
Collapse
|
47
|
Brovarets' OO, Zhurakivsky RO, Hovorun DM. Is the DPT tautomerization of the long A·G Watson-Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question. J Comput Chem 2013; 35:451-66. [PMID: 24382756 DOI: 10.1002/jcc.23515] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2013] [Revised: 11/12/2013] [Accepted: 11/30/2013] [Indexed: 02/04/2023]
Abstract
Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum-mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ↔ A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H-bonds) and the N2H···HC2 dihydrogen bond (DH-bond) (0.68 kcal·mol(-1) ), whereas the A*·G* base mispair-by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H-bonds (0.42 kcal·mol(-1) ). The N2H···HC2 DH-bond smoothly and without bifurcation transforms into the C2H···N2 H-bond at the IRC = -10.07 Bohr in the course of the A·G ↔ A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H-bonds, it was observed that the N6H···O6 H-bond is anticooperative to the two others-N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H-bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10(-14) s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low-frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form.
Collapse
Affiliation(s)
- Ol'ha O Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, 150 Akademika Zabolotnoho Str., 03680, Kyiv, Ukraine ; Research and Educational Center "State Key Laboratory of Molecular and Cell Biology", 150 Akademika Zabolotnoho Str., 03680, Kyiv, Ukraine; Department of Molecular Biology, Biotechnology and Biophysics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, 2-h Akademika Hlushkova Ave., 03022, Kyiv, Ukraine
| | | | | |
Collapse
|
48
|
Brovarets’ OO, Hovorun DM. The physicochemical essence of the purine·pyrimidine transition mismatches with Watson-Crick geometry in DNA: A·C*versaA*·C. A QM and QTAIM atomistic understanding. J Biomol Struct Dyn 2013; 33:28-55. [DOI: 10.1080/07391102.2013.852133] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
49
|
Sutradhar D, Zeegers-Huyskens T, Chandra AK. A theoretical investigation on the conformation and the interaction of CHF₂OCF₂CHF₂ (desflurane II) with one water molecule. J Mol Model 2013; 19:5045-52. [PMID: 24077838 DOI: 10.1007/s00894-013-2016-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2013] [Accepted: 09/15/2013] [Indexed: 10/26/2022]
Abstract
The conformation and the interaction of CHF₂OCF₂CHF₂ (desflurane II) with one water molecule is investigated theoretically using the ab initio MP2/aug-cc-pvdz and DFT-based M062X/6-311++G(d,p) methods. The calculations include the optimized geometries, the harmonic frequencies of relevant vibrational modes along with a natural bond orbital (NBO) analysis including the NBO charges, the hybridization of the C atom and the intra- and intermolecular hyperconjugation energies. In the two most stable conformers, the CH bond of the F2HCO- group occupies the gauche position. The hyperconjugation energies are about the same for both conformers and the conformational preference depends on the interaction between the non-bonded F and H atoms. The deprotonation enthalpies of the CH bonds are about the same for both conformers, the proton affinity of the less stable conformer being 3 kcal mol−1 higher. Both conformers of desflurane II interact with water forming cyclic complexes characterized by CH…O and OH…F hydrogen bonds. The binding energies are moderate, ranging from −2.4 to −3.2 kcal mol−1 at the MP2 level. The origin of the blue shifts of the ν(CH) vibrations is analyzed. In three of the complexes, the water molecule acts as an electron donor. Interestingly, in these cases a charge transfer is also directed to the non bonded OH group of the water molecule. This effect seems to be a property of polyfluorinated ethers.
Collapse
|
50
|
Das G, Lyngdoh RHD. Configuration of wobble base pairs having pyrimidines as anticodon wobble bases: significance for codon degeneracy. J Biomol Struct Dyn 2013; 32:1500-20. [PMID: 23968386 DOI: 10.1080/07391102.2013.824822] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Degeneracy of the genetic code was attributed by Crick to imprecise hydrogen-bonded base-pairing at the wobble position during codon-anticodon pairing. The Crick wobble rules define but do not explain the RNA base pair combinations allowed at this position. We select six pyrimidine bases functioning as anticodon wobble bases (AWBs) to study their H-bonded pairing properties with the four major RNA bases using density functional theory at the B3LYP/6-31G(d,p) level. This is done to assess the extent to which the configuration of a solitary RNA wobble base pair may in itself determine specificity and degeneracy of the genetic code by allowing or disallowing the given base pair during codon-anticodon pairing. Calculated values of select configuration markers for the base pairs screen well between allowed and disallowed base pairs for most cases examined here, where the base pair width emerges as an important factor. A few allowed wobble pairs invoke the involvement of RNA nucleoside conformation, as well as involvement of the exocyclic substituent in H-bonding. This study, however, cannot explain the disallowed status of the Ura⋯Gua wobble pair on the basis of configuration alone. Explanation of the allowed status of the V⋯Ura pair requires further study on the mediatory role of water molecules. Apart from these two cases, these computational results are sufficient, on the basis of base pair configuration alone, to account for the specificity and degeneracy of the genetic code for all known cases of codon-anticodon pairing which involve the pyrimidine AWBs studied here.
Collapse
Affiliation(s)
- Gunajyoti Das
- a Department of Chemistry , North-Eastern Hill University , Shillong , 793022 , India
| | | |
Collapse
|