1
|
Rozenberg M, Fausto R, Reva I. Variable temperature FTIR spectra of polycrystalline purine nucleobases and estimating strengths of individual hydrogen bonds. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 251:119323. [PMID: 33508682 DOI: 10.1016/j.saa.2020.119323] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2020] [Revised: 11/19/2020] [Accepted: 12/05/2020] [Indexed: 06/12/2023]
Abstract
In the first part of this work, we report the FTIR spectra of pure NH and isotopically substituted ND (10-15% D and 80-90% D) polycrystalline hypoxanthine, xanthine, adenine and guanine recorded in the 400-4000 cm-1 range, as a function of temperature (10-300 K). We provide assignments of the stretching and out-of-plane bending amine (NH2) and imine (NH) bands to the distinct H-bonds present in the crystal, based on the temperature sensitivity and isotopic exchange behavior. Empirical correlations between spectral and thermodynamic or structural parameters enabled us to estimate the energies and lengths of H-bonds in the studied nucleobase crystals and to correlate them with literature data. The empirical H-bonding energies are compared with H-bonding and stacking energies computed for hypoxanthine. In the second part, strategies for using the empirical correlations together with information extracted from quantum mechanical data (in particular from the Bader's quantum theory of atoms in molecules, QTAIM) for the evaluation of hydrogen bonding properties are discussed, and their advantages and drawbacks pointed out. The justification for a cooperative use of quantum-mechanical calculations with empirical spectra-energy correlations is discussed.
Collapse
Affiliation(s)
- M Rozenberg
- The Hebrew University of Jerusalem, Department of Inorganic and Analytical Chemistry, Jerusalem, Givat Ram 91904, Israel.
| | - R Fausto
- University of Coimbra, CQC, Department of Chemistry, 3004-535 Coimbra, Portugal.
| | - I Reva
- University of Coimbra, CQC, Department of Chemistry, 3004-535 Coimbra, Portugal; University of Coimbra, CIEPQPF, Department of Chemical Engineering, 3030-790 Coimbra, Portugal.
| |
Collapse
|
2
|
Cuyacot BJR, Durník I, Foroutan-Nejad C, Marek R. Anatomy of Base Pairing in DNA by Interacting Quantum Atoms. J Chem Inf Model 2020; 61:211-222. [PMID: 33112145 DOI: 10.1021/acs.jcim.0c00642] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The formation of purine and pyrimidine base pairs (BPs), which contributes to shaping of the canonical and noncanonical 3D structures of nucleic acids, is one the most investigated phenomena in chemistry and life sciences. In this contribution, the anatomy of the bond energy (BDE) of the base-pairing interaction in 39 different arrangements found experimentally or predicted for DNA structures containing the four common nucleobases (A, C, G, T) in their neutral or protonated forms is described in light of the theory of interacting quantum atoms within the context of the quantum theory of atoms in molecules. The interplay of individual energy components involved in the three stages of the bond formation process (structural deformation, electron-density promotion, and intermolecular interaction) is studied. We recognized that for the neutral BPs, variations in the kinetic and electrostatic contributions to the BDE are rather negligible, leaving the exchange-correlation energy as the main stabilizing component. It is shown that the contribution of the exchange-correlation term can be recovered by including atoms that are formally assumed to be hydrogen bonded (primary interaction). In contrast, to recover the electrostatic component of interaction, one must consider both the primary and secondary (formally nonbonded atoms) interatomic interactions. The results of our study were employed to design new types of BPs with altered bonding anatomy. We demonstrate that improving the electrostatic characteristics of the BPs does not necessarily result in greater interaction energies if weak secondary hydrogen bonding is destroyed. However, the main tuning factor for systems with conserved interacting faces (primary interactions) is the electrostatic component of the interaction energy resulting from the secondary atom-atom electrostatics.
Collapse
Affiliation(s)
- Ben Joseph R Cuyacot
- CEITEC-Central European Institute of Technology, Masaryk University, Kamenice 5/A4, CZ-62500 Brno, Czechia.,Department of Chemistry, Faculty of Science, Masaryk University, Kamenice 5, CZ-62500 Brno, Czechia
| | - Ivo Durník
- CEITEC-Central European Institute of Technology, Masaryk University, Kamenice 5/A4, CZ-62500 Brno, Czechia.,National Center for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, CZ-62500 Brno, Czechia
| | - Cina Foroutan-Nejad
- CEITEC-Central European Institute of Technology, Masaryk University, Kamenice 5/A4, CZ-62500 Brno, Czechia
| | - Radek Marek
- CEITEC-Central European Institute of Technology, Masaryk University, Kamenice 5/A4, CZ-62500 Brno, Czechia.,Department of Chemistry, Faculty of Science, Masaryk University, Kamenice 5, CZ-62500 Brno, Czechia.,National Center for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, CZ-62500 Brno, Czechia
| |
Collapse
|
3
|
Milovanović B, Stanojević A, Etinski M, Petković M. Intriguing Intermolecular Interplay in Guanine Quartet Complexes with Alkali and Alkaline Earth Cations. J Phys Chem B 2020; 124:3002-3014. [DOI: 10.1021/acs.jpcb.0c01165] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Branislav Milovanović
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| | - Ana Stanojević
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| | - Milena Petković
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| |
Collapse
|
4
|
Interaction of (G4)2 and (X4)2 DNA quadruplexes with Cu+, Ag+ and Au+ metal cations: a quantum chemical calculation on structural, energetic and electronic properties. Struct Chem 2019. [DOI: 10.1007/s11224-019-01421-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
5
|
Paragi G, Kupihár Z, Endre G, Fonseca Guerra C, Kovács L. The evaluation of 5-amino- and 5-hydroxyuracil derivatives as potential quadruplex-forming agents. Org Biomol Chem 2018; 15:2174-2184. [PMID: 28054065 DOI: 10.1039/c6ob02574a] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
5-Substituted uracils (NH2 or OH groups in position 5) have been examined theoretically and experimentally as potential building blocks in quadruplex structures. Our high level Density Functional Theory (DFT) calculations showed that the tetramer formation and stacking energies for 5-substituted uracils are similar to the energies of purine-based xanthine (X) or guanine (G) structures. As tetrads of 5-substituted uracils cover almost exactly the same area as purine tetrads, mixed tetrads or quadruplex structures based on X or G and 5-substituted uracil motifs are possible. According to the calculations, 5-hydroxyuracil-based structures are the best candidates for experimental implementation which was corroborated by the existence of higher complexes in the mass spectra of 1-benzyl-5-hydroxyuracil. These pyrimidine-based molecules can be used as efficient building blocks in different applications including aptamers, bio-sensors or - taking into account the larger cavity in the central region of 5-hydroxyuracil structures - as an artificial ion channel.
Collapse
Affiliation(s)
- Gábor Paragi
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Vrije Universiteit Amsterdam (VU), De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands. and MTA-SZTE Supramolecular and Nanostructured Materials Research Group, University of Szeged, Dóm tér 8, 6720 Szeged, Hungary
| | - Zoltán Kupihár
- Department of Medicinal Chemistry, University of Szeged, Dóm tér 8, 6720 Szeged, Hungary.
| | - Gábor Endre
- Department of Medicinal Chemistry, University of Szeged, Dóm tér 8, 6720 Szeged, Hungary.
| | - Célia Fonseca Guerra
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling, Vrije Universiteit Amsterdam (VU), De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands.
| | - Lajos Kovács
- Department of Medicinal Chemistry, University of Szeged, Dóm tér 8, 6720 Szeged, Hungary.
| |
Collapse
|
6
|
Sagi J. In What Ways Do Synthetic Nucleotides and Natural Base Lesions Alter the Structural Stability of G-Quadruplex Nucleic Acids? J Nucleic Acids 2017; 2017:1641845. [PMID: 29181193 PMCID: PMC5664352 DOI: 10.1155/2017/1641845] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2017] [Accepted: 08/15/2017] [Indexed: 01/03/2023] Open
Abstract
Synthetic analogs of natural nucleotides have long been utilized for structural studies of canonical and noncanonical nucleic acids, including the extensively investigated polymorphic G-quadruplexes (GQs). Dependence on the sequence and nucleotide modifications of the folding landscape of GQs has been reviewed by several recent studies. Here, an overview is compiled on the thermodynamic stability of the modified GQ folds and on how the stereochemical preferences of more than 70 synthetic and natural derivatives of nucleotides substituting for natural ones determine the stability as well as the conformation. Groups of nucleotide analogs only stabilize or only destabilize the GQ, while the majority of analogs alter the GQ stability in both ways. This depends on the preferred syn or anti N-glycosidic linkage of the modified building blocks, the position of substitution, and the folding architecture of the native GQ. Natural base lesions and epigenetic modifications of GQs explored so far also stabilize or destabilize the GQ assemblies. Learning the effect of synthetic nucleotide analogs on the stability of GQs can assist in engineering a required stable GQ topology, and exploring the in vitro action of the single and clustered natural base damage on GQ architectures may provide indications for the cellular events.
Collapse
Affiliation(s)
- Janos Sagi
- Rimstone Laboratory, RLI, Carlsbad, CA 92010, USA
| |
Collapse
|
7
|
Havrila M, Stadlbauer P, Islam B, Otyepka M, Šponer J. Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G-Quadruplexes. J Chem Theory Comput 2017; 13:3911-3926. [PMID: 28657760 DOI: 10.1021/acs.jctc.7b00257] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
G-quadruplexes (GQs) are key noncanonical DNA and RNA architectures stabilized by desolvated monovalent cations present in their central channels. We analyze extended atomistic molecular dynamics simulations (∼580 μs in total) of GQs with 11 monovalent cation parametrizations, assessing GQ overall structural stability, dynamics of internal cations, and distortions of the G-tetrad geometries. Majority of simulations were executed with the SPC/E water model; however, test simulations with TIP3P and OPC water models are also reported. The identity and parametrization of ions strongly affect behavior of a tetramolecular d[GGG]4 GQ, which is unstable with several ion parametrizations. The remaining studied RNA and DNA GQs are structurally stable, though the G-tetrad geometries are always deformed by bifurcated H-bonding in a parametrization-specific manner. Thus, basic 10-μs-scale simulations of fully folded GQs can be safely done with a number of cation parametrizations. However, there are parametrization-specific differences and basic force-field errors affecting the quantitative description of ion-tetrad interactions, which may significantly affect studies of the ion-binding processes and description of the GQ folding landscape. Our d[GGG]4 simulations indirectly suggest that such studies will also be sensitive to the water models. During exchanges with bulk water, the Na+ ions move inside the GQs in a concerted manner, while larger relocations of the K+ ions are typically separated. We suggest that the Joung-Cheatham SPC/E K+ parameters represent a safe choice in simulation studies of GQs, though variation of ion parameters can be used for specific simulation goals.
Collapse
Affiliation(s)
- Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Barira Islam
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
8
|
Gresh N, Naseem-Khan S, Lagardère L, Piquemal JP, Sponer JE, Sponer J. Channeling through Two Stacked Guanine Quartets of One and Two Alkali Cations in the Li +, Na +, K +, and Rb + Series. Assessment of the Accuracy of the SIBFA Anisotropic Polarizable Molecular Mechanics Potential. J Phys Chem B 2017; 121:3997-4014. [PMID: 28363025 DOI: 10.1021/acs.jpcb.7b01836] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Stacking of guanine quartets (GQs) can trigger the formation of DNA or RNA quadruple helices, which play numerous biochemical roles. The GQs are stabilized by alkali cations, mainly K+ and Na+, which can reside in, or channel through, the central axis of the GQ stems. Further, ion conduction through GQ wires can be leveraged for nanochemistry applications. G-quadruplex systems have been extensively studied by classical molecular dynamics (MD) simulations using pair-additive force fields or by quantum-chemical (QC) calculations. However, the non-polarizable force fields are very approximate, while QC calculations lack the necessary sampling. Thus, ultimate description of GQ systems would require long-enough simulations using advanced polarizable molecular mechanics (MM). However, to perform such calculations, it is first mandatory to evaluate the method's accuracy using benchmark QC. We report such an evaluation for SIBFA polarizable MM, bearing on the channeling (movement) of an alkali cation (Li+, Na+, K+, or Rb+) along the axis of two stacked G quartets interacting with either one or two ions. The QC energy profiles display markedly different features depending upon the cation but can be retrieved in the majority of cases by the SIBFA profiles. An appropriate balance of first-order (electrostatic and short-range repulsion) and second-order (polarization, charge-transfer, and dispersion) contributions within ΔE is mandatory. With two cations in the channel, the relative weights of the second-order contributions increase steadily upon increasing the ion size. In the G8 complexes with two K+ or two Rb+ cations, the sum of polarization and charge-transfer exceeds the first order terms for all ion positions.
Collapse
Affiliation(s)
- Nohad Gresh
- Laboratoire de Chimie Théorique, Sorbonne Universités , UPMC, UMR7616 CNRS, 75006Paris, France
| | - Sehr Naseem-Khan
- Laboratoire de Chimie Théorique, Sorbonne Universités , UPMC, UMR7616 CNRS, 75006Paris, France
| | - Louis Lagardère
- Laboratoire de Chimie Théorique, Sorbonne Universités , UPMC, UMR7616 CNRS, 75006Paris, France
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, Sorbonne Universités , UPMC, UMR7616 CNRS, 75006Paris, France.,Institut Universitaire de France, Paris Cedex 05, 75231, France.,Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas, 78712, United States
| | - Judit E Sponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Kralovpolska 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jiri Sponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Kralovpolska 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
9
|
Yurenko YP, Novotný J, Marek R. Weak Supramolecular Interactions Governing Parallel and Antiparallel DNA Quadruplexes: Insights from Large-Scale Quantum Mechanics Analysis of Experimentally Derived Models. Chemistry 2017; 23:5573-5584. [PMID: 28225208 DOI: 10.1002/chem.201700236] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2017] [Indexed: 12/30/2022]
Abstract
The topology and energetics of guanine (G) quadruplexes is governed by supramolecular interactions within their strands. In this work, an extensive quantum mechanical (QM) study has been performed to analyze supramolecular interactions that shape the stems of (4+0) parallel (P) and (2+2) antiparallel (AP) quadruplex systems. The large-scale (≈400 atoms) models of P and AP were constructed from high-quality experimental structures. The results provide evidence that each of the P and AP structures is shaped by a distinct network of supramolecular interactions. Analysis of electron topological characteristics of hydrogen bonds in P and AP systems indicates that the P model benefits from stronger intratetrad hydrogen bonding. For intertetrad stacking interactions, both noncovalent interaction plot and energy decomposition analysis approaches suggest that the stem of the P quadruplex benefits more from stacking than that of the AP stem; the difference in energetic stabilization for the two topologies is about 10 %. Stronger hydrogen-bonding and stacking interactions in the stem of the P quadruplex, relative to those in the AP system, can be an important indicator to explain the experimental observations that guanine-rich oligonucleotides tend to form all-parallel stems with an all-anti orientation of nucleobases. However, in addition to intrinsic stabilization, partial desolvation effects, which affect the energetics and dynamics of the G-quadruplex folding process, call for further investigations.
Collapse
Affiliation(s)
- Yevgen P Yurenko
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5/A4, 625 00, Brno, Czech Republic
| | - Jan Novotný
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5/A4, 625 00, Brno, Czech Republic
| | - Radek Marek
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5/A4, 625 00, Brno, Czech Republic
| |
Collapse
|
10
|
Hamlin TA, Poater J, Fonseca Guerra C, Bickelhaupt FM. B-DNA model systems in non-terran bio-solvents: implications for structure, stability and replication. Phys Chem Chem Phys 2017; 19:16969-16978. [DOI: 10.1039/c7cp01908d] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
We have computationally analyzed a comprehensive series of Watson–Crick and mismatched B-DNA base pairs, in the gas phase and in several solvents, including toluene, chloroform, ammonia, methanol and water, using dispersion-corrected density functional theory and implicit solvation.
Collapse
Affiliation(s)
- Trevor A. Hamlin
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM)
- Vrije Universiteit Amsterdam
- NL-1081 HV Amsterdam
- The Netherlands
| | - Jordi Poater
- Departament de Química Inorgànica i Orgánica & Institut de Química Teòrica i Computacional (IQTCUB)
- Universitat de Barcelona
- 08028 Barcelona
- Spain
- ICREA
| | - Célia Fonseca Guerra
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM)
- Vrije Universiteit Amsterdam
- NL-1081 HV Amsterdam
- The Netherlands
- Leiden Institute of Chemistry
| | - F. Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM)
- Vrije Universiteit Amsterdam
- NL-1081 HV Amsterdam
- The Netherlands
- Institute of Molecules and Materials
| |
Collapse
|
11
|
Durec M, Zaccaria F, Fonseca Guerra C, Marek R. Modified Guanines as Constituents of Smart Ligands for Nucleic Acid Quadruplexes. Chemistry 2016; 22:10912-22. [PMID: 27385491 DOI: 10.1002/chem.201601608] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2016] [Indexed: 01/25/2023]
Abstract
Repetitive guanine-rich nucleic acid sequences play a crucial role in maintaining genome stability and the cell life cycle and represent potential targets for regulatory drugs. Recently, it has been demonstrated that guanine-based ligands with a porphyrin core can be used as markers of G-quadruplex assemblies in cell tissues. Herein, model systems of guanine-based ligands are explored by DFT methods. The energies of formation of modified guanine tetrads and those of modified tetrads stacked on the top of natural guanine tetrads have been calculated. The interaction energy has been decomposed into contributions from hydrogen bonding, stacking, and ion coordination and a twist-rise potential energy scan has been performed to find the individual local minima. Energy decomposition analysis reveals the impact of various substituents (F, Cl, Br, I, Me, NMe2 ) on individual energy terms. In addition, cooperative reinforcement in forming the modified and stacked tetrads, as well as the frontier orbitals participating in the hydrogen-bonding framework involving the HOMO-LUMO gap between the occupied σHOMO on the proton-accepting C=O and =N- groups and unoccupied σLUMO on the N-H groups, has been studied. The investigated systems are demonstrated to have a potential in ligand development, mainly due to stacking enhancement compared with natural guanine, which is used as a reference.
Collapse
Affiliation(s)
- Matúš Durec
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5A4, 62500, Brno, Czech Republic.,National Center for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5A4, 62500, Brno, Czech Republic
| | - Francesco Zaccaria
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV, Amsterdam, The Netherlands
| | - Célia Fonseca Guerra
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV, Amsterdam, The Netherlands.
| | - Radek Marek
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5A4, 62500, Brno, Czech Republic. .,National Center for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5A4, 62500, Brno, Czech Republic.
| |
Collapse
|
12
|
Yurenko YP, Novotný J, Nikolaienko TY, Marek R. Nucleotides containing variously modified sugars: energetics, structure, and mechanical properties. Phys Chem Chem Phys 2015; 18:1615-28. [PMID: 26672740 DOI: 10.1039/c5cp05478h] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The influence of various sugar residue modifications on intrinsic energetic, conformational, and mechanical properties of 2'-deoxyribonucleotide-5'-monophosphates (dNs) was comprehensively investigated using modern quantum chemical approaches. In total, fourteen sugar modifications, including double bonds and heteroatoms (S and N) inside the sugar ring, as well as fluorination in various positions, were analyzed. Among hundreds of possible conformational states of dNs, only two - AI and BI, corresponding to the most biologically significant forms of a double-helical DNA, were considered for each dN. It was established that the most of the studied modifications tend to strongly stabilize either AI or BI conformation of dNs both in the gas phase and in aqueous solution (modelled by implicit solvent models). Therefore, some of these modifications can be used as a tool for reducing structural polymorphism of nucleic acids in solution as well as for designing oligonucleotides with specific structural features. The evaluation of relaxed force constants (RFC) for glycosidic bonds suggests that the majority of the studied modifications of the sugar residue yield increased strengths of glycosidic bonds in dNs, and can therefore be used for designing modified nucleic acids with an increased resistance to abasic lesions. The most significant reinforcement of the glycosidic bond occurs in dNs containing the CF2 group instead of the O4' oxygen and the fluorine atom at the 2'-α-position. The calculation of the RFC and vibrational root-mean-square (VRMS) deviations for conformational degrees of freedom revealed a strong dependence between mechanical properties of dNs and their energetic characteristics. In particular, electronic energies of AI and BI conformers of dNs calculated in vacuo are closely connected with the values of relaxed force constants (RFC) for the δ angle: the higher RFC(δ) values correspond to more energetically favorable conformers.
Collapse
Affiliation(s)
- Yevgen P Yurenko
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5, CZ-62500 Brno, Czech Republic.
| | | | | | | |
Collapse
|
13
|
Yurenko YP, Novotný J, Sklenář V, Marek R. Substituting CF2 for O4' in Components of Nucleic Acids: Towards Systems with Reduced Propensity to Form Abasic Lesions. Chemistry 2015; 21:17933-43. [PMID: 26493955 DOI: 10.1002/chem.201502977] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2015] [Indexed: 01/22/2023]
Abstract
Intrinsic structural features and energetics of nucleotides containing variously fluorinated sugars as potential building blocks of DNA duplexes and quadruplexes are explored systematically using the modern methods of density functional theory (DFT) and quantum chemical topology (QCT). Our results suggest that fluorination at the 2'-β or 2'-α,β positions somewhat stabilizes in vacuo the AI relative to the BI conformations. In contrast, substitution of the CF2 group for the O4' atom (O4'-CF2 modification) leads to a preference of the BI relative to AI DNA-like conformers. All the studied modifications result in a noticeable increase in the stability of the glycosidic bond [estimated by the relaxed force constants (RFC) approach], with particularly encouraging results for the O4'-CF2 derivative. Consequently, the O4'-CF2 modified systems are suggested and explored as promising scaffolds for the development of duplex and quadruplex structures with reduced propensity to form abasic lesions and to undergo DNA damage.
Collapse
Affiliation(s)
- Yevgen P Yurenko
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5, 62500 Brno (Czech Republic).
| | - Jan Novotný
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5, 62500 Brno (Czech Republic)
| | - Vladimir Sklenář
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5, 62500 Brno (Czech Republic).,National Center for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 62500 Brno (Czech Republic).,Department of Chemistry, Faculty of Science, Masaryk University, Kamenice 5, 62500 Brno (Czech Republic)
| | - Radek Marek
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5, 62500 Brno (Czech Republic). .,National Center for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 62500 Brno (Czech Republic). .,Department of Chemistry, Faculty of Science, Masaryk University, Kamenice 5, 62500 Brno (Czech Republic).
| |
Collapse
|