1
|
Tarek Ibrahim M, Wait E, Ren P. Quantum Mechanics Characterization of Non-Covalent Interaction in Nucleotide Fragments. Molecules 2024; 29:3258. [PMID: 39064837 PMCID: PMC11279843 DOI: 10.3390/molecules29143258] [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: 05/31/2024] [Revised: 07/03/2024] [Accepted: 07/06/2024] [Indexed: 07/28/2024] Open
Abstract
Accurate calculation of non-covalent interaction energies in nucleotides is crucial for understanding the driving forces governing nucleic acid structure and function, as well as developing advanced molecular mechanics forcefields or machine learning potentials tailored to nucleic acids. Here, we dissect the nucleotides' structure into three main constituents: nucleobases (A, G, C, T, and U), sugar moieties (ribose and deoxyribose), and phosphate group. The interactions among these fragments and between fragments and water were analyzed. Different quantum mechanical methods were compared for their accuracy in capturing the interaction energy. The non-covalent interaction energy was decomposed into electrostatics, exchange-repulsion, dispersion, and induction using two ab initio methods: Symmetry-Adapted Perturbation Theory (SAPT) and Absolutely Localized Molecular Orbitals (ALMO). These calculations provide a benchmark for different QM methods, in addition to providing a valuable understanding of the roles of various intermolecular forces in hydrogen bonding and aromatic stacking. With SAPT, a higher theory level and/or larger basis set did not necessarily give more accuracy. It is hard to know which combination would be best for a given system. In contrast, ALMO EDA2 did not show dependence on theory level or basis set; additionally, it is faster.
Collapse
Affiliation(s)
- Mayar Tarek Ibrahim
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, USA;
| | - Elizabeth Wait
- Interdisciplinary Life Sciences Graduate Program, The University of Texas at Austin, Austin, TX 78712, USA;
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, USA;
- Interdisciplinary Life Sciences Graduate Program, The University of Texas at Austin, Austin, TX 78712, USA;
| |
Collapse
|
2
|
Percot A, Maurel MC, Lambert JF, Zins EL. New insights into the surface Enhanced Raman Scattering (SERS) response of adenine using chemometrics. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 314:124177. [PMID: 38554690 DOI: 10.1016/j.saa.2024.124177] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Revised: 02/15/2024] [Accepted: 03/18/2024] [Indexed: 04/02/2024]
Abstract
The SERS response of adenine is one of the most studied, due to its outstanding exaltation. However, the spectra obtained strongly differ according to the experimental conditions and still remain not well understood. To be able to search for the presence of this molecule in complex environments, it is essential to better understand the SERS response of adenine alone. After a brief presentation of the literature on the subject, we present results suggesting that the experimental spectra would result from the overlap of different spectroscopic signatures, that may probably be due to different non-covalent interactions or different electromagnetic fields experienced by adenine molecules. An independent component analysis is reported. Our results underline the difficulty to precisely analyze the experimental data, the need to continue this research and to constitute data banks that would allow comparing the spectra obtained in different laboratories according to the experimental conditions.
Collapse
Affiliation(s)
- A Percot
- Sorbonne Université, CNRS, MONARIS, UMR8233, F-75005 Paris, France.
| | - M C Maurel
- Sorbonne Université, Muséum National d'Histoire Naturelle, CNRS,MNHN, UMR7205, ISYEB, F-75005 Paris, France
| | - J F Lambert
- Sorbonne Université, CNRS, LAMS, UMR8220, F-75005 Paris, France
| | - E L Zins
- Sorbonne Université, CNRS, MONARIS, UMR8233, F-75005 Paris, France
| |
Collapse
|
3
|
Cataldo PG, Iramain MA, Castillo MV, Manzur ME, Romano E, Brandán SA. Vibrational assignments of cyclic dimers and inter-monomers of adenine relating FT-IR, FT-Raman and UV spectra with SQMFF and DFT calculations. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 306:123540. [PMID: 37898058 DOI: 10.1016/j.saa.2023.123540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 10/13/2023] [Accepted: 10/14/2023] [Indexed: 10/30/2023]
Abstract
In this work, three different cyclic dimers and a tetramer of adenine taken from the experimental structure determined by X-ray diffraction have been studied by combination of experimental FT-IR, FT-Raman and UV-Visible spectra with hybrid B3LYP/6-311++G** and scaled quantum mechanical force field (SQMFF) calculations in order to perform the complete assignments of bands observed in the vibrational spectra. The characteristics of different N-H···N interactions of those three cyclic structures together with the group of IR bands observed between 2865 and 2599 cm-1 have been elucidated considering the tetrameric structure. The cyclic dimers and the tetramer of adenine confirm that the bands observed between 2865 and 2599 cm-1 are not due to N-H···N interactions but to bands of combination, as was previously suggested. The experimental available deuterated IR and terahertz spectra have allowed the complete assignments of regions of higher and lower wavenumbers. Good correlations were acquired comparing the theoretical IR, Raman and UV spectra of three species and the tetramer with the analogous experimental ones, suggesting the presence of all species in both phases. Vibronic bands are observed in the electronic spectra when adenine concentration is increased in aqueous solution evidencing the presence of monomer, tautomers and dimers, as reported by different studies. Similar characteristics of H bonds interactions are predicted for dimers 1 and 2 but different from the dimer 3, as revealed by using NBO and AIM calculations. Different scaled force constants values were found for the cyclic dimers 1 and 2, as compared to the corresponding to dimer 3.
Collapse
Affiliation(s)
- Pablo G Cataldo
- Cátedra de Química General, Instituto de Química Inorgánica, Facultad de Bioquímica. Química y Farmacia, Universidad Nacional de Tucumán, Ayacucho 471, (4000) San Miguel de Tucumán, Tucumán, Argentina
| | - Maximiliano A Iramain
- Cátedra de Química General, Instituto de Química Inorgánica, Facultad de Bioquímica. Química y Farmacia, Universidad Nacional de Tucumán, Ayacucho 471, (4000) San Miguel de Tucumán, Tucumán, Argentina
| | - María V Castillo
- Cátedra de Química General, Instituto de Química Inorgánica, Facultad de Bioquímica. Química y Farmacia, Universidad Nacional de Tucumán, Ayacucho 471, (4000) San Miguel de Tucumán, Tucumán, Argentina
| | - María E Manzur
- Cátedra de Química General, Instituto de Química Inorgánica, Facultad de Bioquímica. Química y Farmacia, Universidad Nacional de Tucumán, Ayacucho 471, (4000) San Miguel de Tucumán, Tucumán, Argentina
| | - Elida Romano
- Cátedra de Química General, Instituto de Química Inorgánica, Facultad de Bioquímica. Química y Farmacia, Universidad Nacional de Tucumán, Ayacucho 471, (4000) San Miguel de Tucumán, Tucumán, Argentina
| | - Silvia Antonia Brandán
- Cátedra de Química General, Instituto de Química Inorgánica, Facultad de Bioquímica. Química y Farmacia, Universidad Nacional de Tucumán, Ayacucho 471, (4000) San Miguel de Tucumán, Tucumán, Argentina.
| |
Collapse
|
4
|
Jezuita A, Wieczorkiewicz PA, Krygowski TM, Szatylowicz H. Influence of the Solvent on the Stability of Aminopurine Tautomers and Properties of the Amino Group. Molecules 2023; 28:molecules28072993. [PMID: 37049758 PMCID: PMC10095612 DOI: 10.3390/molecules28072993] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 03/23/2023] [Accepted: 03/24/2023] [Indexed: 03/30/2023] Open
Abstract
Amino derivatives of purine (2-, 6-, 8-, and N-NH2) have found many applications in biochemistry. This paper presents the results of a systematic computational study of the substituent and solvent effects in these systems. The issues considered are the electron-donating properties of NH2, its geometry, π-electron delocalization in purine rings and tautomeric stability. Calculations were performed in ten environments, with 1 < ε < 109, using the polarizable continuum model of solvation. Electron-donating properties were quantitatively described by cSAR (charge of the substituent active region) parameter and π-electron delocalization by using the HOMA (harmonic oscillator model of aromaticity) index. In aminopurines, NH2 proximity interactions depend on its position and the tautomer. The results show that they are the main factor determining how solvation affects the electron-donating strength and geometry of NH2. Proximity with the NH∙∙∙HN repulsive interaction between the NH2 and endocyclic NH group results in stronger solvent effects than the proximity with two attractive NH∙∙∙N interactions. The effect of amino and nitro (previously studied) substitution on aromaticity was compared; these two groups have, in most cases, the opposite effect, with the largest being in N1H and N3H purine tautomers. The amino group has a smaller effect on the tautomeric preferences of purine than the nitro group. Only in 8-aminopurine do tautomeric preferences change: N7H is more stable than N9H in H2O.
Collapse
Affiliation(s)
- Anna Jezuita
- Faculty of Science and Technology, Jan Dlugosz University in Czestochowa, Al. Armii Krajowej 13/15, 42-200 Czestochowa, Poland
| | - Paweł A. Wieczorkiewicz
- Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664 Warsaw, Poland
- Correspondence: (P.A.W.); (H.S.)
| | | | - Halina Szatylowicz
- Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664 Warsaw, Poland
- Correspondence: (P.A.W.); (H.S.)
| |
Collapse
|
5
|
Stone S, Ray D, Andricioaei I. Force-Field-Dependent DNA Breathing Dynamics: A Case Study of Hoogsteen Base Pairing in A6-DNA. J Chem Inf Model 2022; 62:6749-6761. [PMID: 36049242 PMCID: PMC9795553 DOI: 10.1021/acs.jcim.2c00519] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The Hoogsteen (HG) base pairing conformation, commonly observed in damaged and mutated DNA helices, facilitates DNA repair and DNA recognition. The free energy difference between HG and Watson-Crick (WC) base pairs has been computed in previous studies. However, the mechanism of the conformational transition is not well understood. A detailed understanding of the process of WC to HG base pair transition can provide a deeper understanding of DNA repair and recognition. In an earlier study, we explored the free energy landscape for this process using extensive computer simulation with the CHARMM36 force field. In this work, we study the impact of force field models in describing the WC to HG base pairing transition using meta-eABF enhanced sampling, quasi-harmonic entropy calculation, and nonbonded energy analysis. The secondary structures of both base pairing forms and the topology of the free energy landscapes were consistent over different force field models, although the relative free energy, entropy, and the interaction energies tend to vary. The relative stability of the WC and HG conformations is dictated by a delicate balance between the enthalpic stabilization and the reduced entropy of the structurally rigid HG structure. These findings highlight the impact that subtleties in force field models can have on accurately modeling DNA base pair dynamics and should stimulate further computational investigations into other dynamically important motions in DNA.
Collapse
Affiliation(s)
- Sharon
Emily Stone
- Department
of Chemistry, University of California Irvine, Irvine, California 92697, United States
| | - Dhiman Ray
- Department
of Chemistry, University of California Irvine, Irvine, California 92697, United States
| | - Ioan Andricioaei
- Department
of Chemistry, University of California Irvine, Irvine, California 92697, United States,Department
of Physics and Astronomy, University of
California Irvine, Irvine, California 92697, United States,
| |
Collapse
|
6
|
Intramolecular Interactions in Derivatives of Uracil Tautomers. Molecules 2022; 27:molecules27217240. [PMID: 36364066 PMCID: PMC9656941 DOI: 10.3390/molecules27217240] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Revised: 10/17/2022] [Accepted: 10/21/2022] [Indexed: 11/17/2022] Open
Abstract
The influence of solvents on intramolecular interactions in 5- or 6-substituted nitro and amino derivatives of six tautomeric forms of uracil was investigated. For this purpose, the density functional theory (B97-D3/aug-cc-pVDZ) calculations were performed in ten environments (1 > ε > 109) using the polarizable continuum model (PCM) of solvation. The substituents were characterized by electronic (charge of the substituent active region, cSAR) and geometric parameters. Intramolecular interactions between non-covalently bonded atoms were investigated using the theory of atoms in molecules (AIM) and the non-covalent interaction index (NCI) method, which allowed discussion of possible interactions between the substituents and N/NH endocyclic as well as =O/−OH exocyclic groups. The nitro group was more electron-withdrawing in the 5 than in the 6 position, while the opposite effect was observed in the case of electron donation of the amino group. These properties of both groups were enhanced in polar solvents; the enhancement depended on the ortho interactions. Substitution or solvation did not change tautomeric preferences of uracil significantly. However, the formation of a strong NO∙∙∙HO intramolecular hydrogen bond in the 5-NO2 derivative stabilized the dienol tautomer from +17.9 (unsubstituted) to +5.4 kcal/mol (substituted, energy relative to the most stable diketo tautomer).
Collapse
|
7
|
Yuan Y, Yan H, Cui Z, Liu Z, Su W, Zhang R. Quantum Chemical Calculations with Machine Learning for Multipolar Electrostatics Prediction in RNA: An Application to Pentose. J Chem Inf Model 2022; 62:4122-4133. [PMID: 36036609 DOI: 10.1021/acs.jcim.2c00747] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To develop a realistic electrostatic model that allows for the anisotropy of the atomic electron density, high-rank atomic multipole moments computed by quantum chemical calculations have been studied extensively. However, it is hard to process huge RNA systems only relying on quantum chemical calculations due to its highly computational cost. In this study, we employ five machine learning methods of Gaussian process regression with automatic relevance determination (ARDGPR), Kriging, radial basis function neural networks, Bagging, and generalized regression neural network to predict atomic multipole moments. Atom-atom electrostatic interaction energies are subsequently computed using the predicted atomic multipole moments in the pilot system pentose of RNA. Here, the performance of the five methods is compared in terms of both the multipole moment prediction errors and the electrostatic energy prediction errors. For the predicted high-rank multipole moments of the four elements (O, C, N, and H) in capped pentose, ARDGPR and Kriging consistently outperform the other three methods. Therefore, the multipole moments predicted by the two best methods of ARDGPR and Kriging are then used to predict electrostatic interaction energy of each pentose. Finally, the absolute average energy errors of ARDGPR and Kriging are 1.83 and 4.33 kJ mol-1, respectively. Compared to Kriging, the ARDGPR method achieves a 58% decrease in the absolute average energy error. These satisfactory results demonstrated that the ARDGPR method with the strong feature extraction ability can predict the electrostatic interaction energy of pentose in RNA correctly and reliably.
Collapse
Affiliation(s)
- Yongna Yuan
- School of Information Science & Engineering, Lanzhou University, Lanzhou, China, 730000
| | - Haoqiu Yan
- School of Information Science & Engineering, Lanzhou University, Lanzhou, China, 730000
| | - Zeyang Cui
- School of Information Science & Engineering, Lanzhou University, Lanzhou, China, 730000
| | - Zhenyu Liu
- School of Cyberspace Security, Gansu University of Political Science and Law, Lanzhou, China, 730070
| | - Wei Su
- School of Information Science & Engineering, Lanzhou University, Lanzhou, China, 730000
| | - Ruisheng Zhang
- School of Information Science & Engineering, Lanzhou University, Lanzhou, China, 730000
| |
Collapse
|
8
|
Bhai S, Ganguly B. Role of pH in the stability of cytosine-cytosine mismatch and canonical AT and GC base pairs mediated with silver ion: a DFT study. Struct Chem 2021. [DOI: 10.1007/s11224-021-01814-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
9
|
Jezuita A, Wieczorkiewicz PA, Szatylowicz H, Krygowski TM. Effect of the Solvent and Substituent on Tautomeric Preferences of Amine-Adenine Tautomers. ACS OMEGA 2021; 6:18890-18903. [PMID: 34337229 PMCID: PMC8320138 DOI: 10.1021/acsomega.1c02118] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/21/2021] [Accepted: 06/23/2021] [Indexed: 06/13/2023]
Abstract
Adenine is one of the basic molecules of life; it is also an important building block in the synthesis of new pharmaceuticals, electrochemical (bio)sensors, or self-assembling molecular materials. Therefore, it is important to know the effects of the solvent and substituent on the electronic structure of adenine tautomers and their stability. The four most stable adenine amino tautomers (9H, 7H, 3H, and 1H), modified by substitution (C2- or C8-) of electron-withdrawing NO2 and electron-donating NH2 groups, are studied theoretically in the gas phase and in solvents of different polarities (1 ≤ ε < 109). Solvents have been modeled using the polarizable continuum model. Comparison of the stability of substituted adenine tautomers in various solvents shows that substitution can change tautomeric preferences with respect to the unsubstituted adenine. Moreover, C8 substitution results in slight energy differences between tautomers in polar solvents (<1 kcal/mol), which suggests that in aqueous solution, C8-X-substituted adenine systems may consist of a considerable amount of two tautomers-9H and 7H for X = NH2 and 3H and 9H for X = NO2. Furthermore, solvation enhances the effect of the nitro group; however, the enhancement strongly depends on the proximity effects. This enhancement for the NO2 group with two repulsive N···ON contacts can be threefold higher than that for the NO2 with one attractive NH···ON contact. The proximity effects are even more significant for the NH2 group, as the solvation may increase or decrease its electron-donating ability, depending on the type of proximity.
Collapse
Affiliation(s)
- Anna Jezuita
- Faculty of Science and Technology, Jan Dlugosz University in Czestochowa, Al. Armii Krajowej 113/15, 42-200 Czestochowa, Poland
| | | | - Halina Szatylowicz
- Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664 Warsaw, Poland
| | | |
Collapse
|
10
|
Solvent Effect on the Stability and Reverse Substituent Effect in Nitropurine Tautomers. Symmetry (Basel) 2021. [DOI: 10.3390/sym13071223] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
The solvent effect on the stability and electron-accepting properties (EA) of the nitro group attached to the C2, C6, or C8 position of nitropurine NH tautomers is investigated. For this purpose, the density functional theory (DFT) and the polarizable continuum model (PCM) of solvation in a wide range of solvents (1 < ε < 109) are used. We show that the EA properties of the NO2 group, described by the charge of the substituent active region (cSAR) model, are linearly dependent on the reciprocal of the solvent dielectric constant; in all cases, solvation enhances the EA properties of this group. Furthermore, the sensitivity of EA properties of the nitro group to the solvent effect depends on the proximity effects. It has been shown that the proximity of two endocyclic N atoms (two repulsive interactions) results in higher sensitivity than the asymmetric proximity of the endocyclic N atom and NH group (one repulsive and one attractive interaction). To explain this phenomenon, the geometry of the nitro group in coplanar form and after forcing its rotation around the CN bond is discussed. Relative stabilities of nitropurine tautomers in different solvents are also presented. Differences in the stabilities and solvation energies are explained by aromaticity, electronic structure, and intramolecular interactions of the nitropurine tautomers.
Collapse
|
11
|
Otep S, Ogita K, Yomogita N, Motai K, Wang Y, Tseng YC, Chueh CC, Hayamizu Y, Matsumoto H, Ishikawa K, Mori T, Michinobu T. Cross-Linking of Poly(arylenebutadiynylene)s and Its Effect on Charge Carrier Mobilities in Thin-Film Transistors. Macromolecules 2021. [DOI: 10.1021/acs.macromol.1c00008] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sultan Otep
- Department of Materials Science and Engineering, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| | - Kosuke Ogita
- Department of Materials Science and Engineering, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| | - Naomasa Yomogita
- Department of Materials Science and Engineering, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| | - Kazunori Motai
- Department of Materials Science and Engineering, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| | - Yang Wang
- Department of Materials Science, Fudan University, 2005, Songhu Road, Shanghai 200433, China
| | - Yu-Cheng Tseng
- Department of Chemical Engineering and Advanced Research Center for Green Materials Science and Technology, National Taiwan University, Taipei 10617, Taiwan
| | - Chu-Chen Chueh
- Department of Chemical Engineering and Advanced Research Center for Green Materials Science and Technology, National Taiwan University, Taipei 10617, Taiwan
| | - Yuhei Hayamizu
- Department of Materials Science and Engineering, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| | - Hidetoshi Matsumoto
- Department of Materials Science and Engineering, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| | - Ken Ishikawa
- Department of Materials Science and Engineering, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| | - Takehiko Mori
- Department of Materials Science and Engineering, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| | - Tsuyoshi Michinobu
- Department of Materials Science and Engineering, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| |
Collapse
|
12
|
Jhunjhunwala A, Ali Z, Bhattacharya S, Halder A, Mitra A, Sharma P. On the Nature of Nucleobase Stacking in RNA: A Comprehensive Survey of Its Structural Variability and a Systematic Classification of Associated Interactions. J Chem Inf Model 2021; 61:1470-1480. [PMID: 33570947 DOI: 10.1021/acs.jcim.0c01225] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The astonishing diversity in folding patterns of RNA three-dimensional (3D) structures is crafted by myriads of noncovalent contacts, of which base pairing and stacking are the most prominent. A systematic and comprehensive classification and annotation of these interactions is necessary for a molecular-level understanding of their roles. However, unlike in the case of base pairing, where a widely accepted nomenclature and classification scheme exists in the public domain, currently available classification schemes for base-base stacking need major enhancements to comprehensively capture the necessary features underlying the rich stacking diversity in RNA. Here, we extend the previous stacking classification based on nucleobase interacting faces by introducing a structurally intuitive geometry-cum topology-based scheme. Specifically, a stack is first classified in terms of the geometry described by the relative orientation of the glycosidic bonds, which generates eight basic stacking geometric families for heterodimeric stacks and six of those for homodimeric stacks. Further annotation in terms of the identity of the bases and the region of involvement of purines (five-membered, six-membered, or both rings) leads to the enumeration of 384 distinct RNA base stacks. Based on our classification scheme, we present an algorithm for automated identification of stacks in RNA crystal structures and analyze the stacking context in selected RNA structures. Overall, the work described here is expected to greatly facilitate the structure-based RNA research.
Collapse
Affiliation(s)
- Ayush Jhunjhunwala
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Gachibowli, Hyderabad, Telangana 500032, India
| | - Zakir Ali
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Sohini Bhattacharya
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Gachibowli, Hyderabad, Telangana 500032, India
| | - Antarip Halder
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Gachibowli, Hyderabad, Telangana 500032, India
| | - Abhijit Mitra
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Gachibowli, Hyderabad, Telangana 500032, India
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| |
Collapse
|
13
|
Wieczorkiewicz PA, Szatylowicz H, Krygowski TM. Mutual Relations between Substituent Effect, Hydrogen Bonding, and Aromaticity in Adenine-Uracil and Adenine-Adenine Base Pairs. Molecules 2020; 25:E3688. [PMID: 32823565 PMCID: PMC7464026 DOI: 10.3390/molecules25163688] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2020] [Revised: 07/30/2020] [Accepted: 08/11/2020] [Indexed: 11/16/2022] Open
Abstract
The electronic structure of substituted molecules is governed, to a significant extent, by the substituent effect (SE). In this paper, SEs in selected nucleic acid base pairs (Watson-Crick, Hoogsteen, adenine-adenine) are analyzed, with special emphasis on their influence on intramolecular interactions, aromaticity, and base pair hydrogen bonding. Quantum chemistry methods-DFT calculations, the natural bond orbital (NBO) approach, the Harmonic Oscillator Model of Aromaticity (HOMA) index, the charge of the substituent active region (cSAR) model, and the quantum theory of atoms in molecules (QTAIM)-are used to compare SEs acting on adenine moiety and H-bonds from various substitution positions. Comparisons of classical SEs in adenine with those observed in para- and meta-substituted benzenes allow for the better interpretation of the obtained results. Hydrogen bond stability and its other characteristics (e.g., covalency) can be significantly changed as a result of the SE, and its consequences are dependent on the substitution position. These changes allow us to investigate specific relations between H-bond parameters, leading to conclusions concerning the nature of hydrogen bonding in adenine dimers-e.g., H-bonds formed by five-membered ring nitrogen acceptor atoms have an inferior, less pronounced covalent nature as compared to those formed by six-membered ring nitrogen. The energies of individual H-bonds (obtained by the NBO method) are analyzed and compared to those predicted by the Espinosa-Molins-Lecomte (EML) model. Moreover, both SE and H-bonds can significantly affect the aromaticity of adenine rings; long-distance SEs on π-electron delocalization are also documented.
Collapse
Affiliation(s)
- Paweł A. Wieczorkiewicz
- Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664 Warsaw, Poland;
| | - Halina Szatylowicz
- Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664 Warsaw, Poland;
| | | |
Collapse
|
14
|
Jezuita A, Szatylowicz H, Krygowski TM. Impact of the Substituents on the Electronic Structure of the Four Most Stable Tautomers of Purine and Their Adenine Analogues. ACS OMEGA 2020; 5:11570-11577. [PMID: 32478247 PMCID: PMC7254788 DOI: 10.1021/acsomega.0c00820] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/24/2020] [Accepted: 04/21/2020] [Indexed: 06/11/2023]
Abstract
Substituent effects at the C2-, C8-, and N-positions of adenine and purine on the structural and π-electronic changes in their four tautomers were studied using the B97D3/aug-cc-pvdz computational level. The effect of various substituents (NO2, CN, CHO, Cl, F, H, Me, OMe, OH, and NH2) was characterized by the charge of the substituent active region (cSAR) approach and Hammett substituent constants σ. It has been found that for both adenine and purine derivatives, substituents from the C8-X position have a stronger influence on their electronic structure than from the C2-X and N-X positions. The presence of the amino group in adenine enhances the substituent effect compared to that which occurs in purine. In addition, its electronic structure is more sensitive to the effect of the substituent in 3H and 1H than in the 9H and 7H adenine tautomers. For a given substituent, a large variation in cSAR(X) values is observed, strongly dependent on the substitution position. For 7H and 9H adenine tautomers for C8-X systems, substituents reduce the aromaticity of the five-membered rings but increase the aromaticity of the six-membered rings.
Collapse
Affiliation(s)
- Anna Jezuita
- Faculty
of Chemistry, Opole University, Oleska 48, 45-052 Opole, Poland
| | - Halina Szatylowicz
- Faculty
of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664 Warsaw, Poland
| | | |
Collapse
|
15
|
Ponikvar-Svet M, Zeiger DN, Liebman JF. Interplay of thermochemistry and Structural Chemistry, the journal (volume 30, 2019, issues 1–2) and the discipline. Struct Chem 2020. [DOI: 10.1007/s11224-020-01494-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
16
|
Gop S, Sutradhar R, Chakraborty S, Sinha TP. Tautomeric effect of guanine on stability, spectroscopic and absorbance properties in cytosine–guanine base pairs: a DFT and TD-DFT perspective. Theor Chem Acc 2020. [DOI: 10.1007/s00214-020-2551-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
17
|
Szatylowicz H, Jezuita A, Marek PH, Krygowski TM. Substituent effects on the stability of the four most stable tautomers of adenine and purine. RSC Adv 2019; 9:31343-31356. [PMID: 35527924 PMCID: PMC9072591 DOI: 10.1039/c9ra04615a] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2019] [Accepted: 09/20/2019] [Indexed: 11/21/2022] Open
Abstract
Substituent effects at the C2-, C8- and N-positions of adenine and purine in their four the most stable tautomers are studied by means of B97D3/aug-cc-pvdz computation applying substituents of varying electronic properties: NO2, CN, CHO, Cl, F, H, Me, OMe, OH and NH2. The substituent effect is characterized by the substituent effect stabilization energy (SESE) and substituent Hammett constant σ. For adenine, SESE is obtained with purine as the reference system. Additionally, for both adenine and purine, SESE characteristics are estimated with benzene, imidazole and amino-pyrimidine as reference systems, when possible, taking into account substitution in topologically equivalent positions. The role of a C6-NH2 group in adenine in modifying the substitution effect is observed and discussed. Additionally, the proximity effect for some asymmetric substituents (e.g. CHO, OMe) is recognized and meticulously analyzed.
Collapse
Affiliation(s)
- Halina Szatylowicz
- Faculty of Chemistry, Warsaw University of Technology Noakowskiego 3 00-664 Warsaw Poland
| | - Anna Jezuita
- Faculty of Chemistry, Opole University Oleska 48 45-052 Opole Poland
| | - Paulina H Marek
- Faculty of Chemistry, Warsaw University of Technology Noakowskiego 3 00-664 Warsaw Poland .,Faculty of Chemistry, University of Warsaw Pasteura 1 02-093 Warsaw Poland
| | | |
Collapse
|
18
|
Jezuita A, Szatylowicz H, Marek PH, Krygowski TM. Aromaticity of the most stable adenine and purine tautomers in terms of Hückel's 4N+2 principle. Tetrahedron 2019. [DOI: 10.1016/j.tet.2019.130474] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|