Intramolecular Interactions in Derivatives of Uracil Tautomers

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).

Structure-activity features of purines and their receptors: implications in cell physiopathology

The purine molecular structure consists of fused pyrimidine and imidazole rings. Purines are main pieces that conform the structure of nucleic acids which rule the inheritance processes. Purines also work as metabolic intermediates in different cell functions and as messengers in the signaling pathways throughout cellular communication. Purines, mainly ATP and adenosine (ADO), perform their functional and pharmacological properties because of their structural/chemical characteristics that make them either targets of mutagenesis, mother frameworks for designing molecules with controlled effects (e.g. anti-cancer), or chemical donors (e.g., of methyl groups, which represent a potential chemoprotective action against cancer). Purines functions also come from their effect on specific receptors, channel-linked and G-protein coupled for ATP, and exclusively G-coupled receptors for ADO (also known as ADORAs), which are involved in cell signaling pathways, there, purines work as chemical messengers with autocrine, paracrine, and endocrine actions that regulate cell metabolism and immune response in tumor progression which depends on the receptor types involved in these signals. Purines also have antioxidant and anti-inflammatory properties and participate in the cell energy homeostasis. Therefore, purine physiology is important for a variety of functions relevant to cellular health; thus, when these molecules present a homeostatic imbalance, the stability and survival of the cellular systems become compromised.

Effect of the Solvent and Substituent on Tautomeric Preferences of Amine-Adenine Tautomers

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.

Solvent Effect on the Stability and Reverse Substituent Effect in Nitropurine Tautomers

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.

Impact of the Substituents on the Electronic Structure of the Four Most Stable Tautomers of Purine and Their Adenine Analogues

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.

Purine tautomeric preferences and bond-length alternation in relation with protonation-deprotonation and alkali metal cationization

Quantum chemical calculations were carried out for deprotonated (P-) and protonated purine (PH+) and for adducts with one alkali metal cation (P-M+ and PM+, where M+ is Li+ or Na+) in the gas phase {B3LYP/6-311+G(d,p)}, a model of perfectly apolar environment, and for selected structures in aqueous solution {PCM(water)//B3LYP/6-311+G(d,p)}, a reference polar medium for biological studies. All potential isomers of purine derivatives were considered, the favored structures indicated, and the preferred sites for protonation/deprotonation and cationization reactions determined. Proton and metal cation basicities of purine in the gas phase were discussed and compared with those of imidazole and pyrimidine. Bond-length alternations in the P, PH+, P-M+, and PM+ forms were quantitatively measured using the harmonic oscillator model of electron delocalization (HOMED) indices and compared with those for P. Variations of the HOMED values when proceeding from the purine structural building blocks, pyrimidine and imidazole, to the bicyclic purine system were also examined. Generally, the isolated NH isomers exhibit a strongly delocalized π-system (HOMED > 0.8). Deprotonation slightly increases the HOMED values, whereas protonation and cationization change the HOMED indices in different way. For bidentate M+-adducts, the HOMED values are larger than 0.9 like for the largely delocalized P-. The HOMED values correlate well in a comprehensive relationship with the relative Gibbs energies (ΔG) calculated for individual isomers whatever the purine form is, neutral, protonated, or cationized. When PCM-DFT model was utilized for P-, PH+, PM+, and P-M+ (M+ = Li+) both electron delocalization and relative stability are different from those for the molecules in vacuo. The solvation effects cause a slight increase in HOMEDs, whereas the ΔEs decrease, but in different ways. Hence, contribution of particular isomers in the isomeric mixtures of PH+, PM+, and P-M+ also varies. HOMED variations for the favored neutral, deprotonated, protonated, and lithiated forms of purine in the gas phase and aqueous solution.

Substituent effects on the stability of the four most stable tautomers of adenine and purine

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: NO₂, CN, CHO, Cl, F, H, Me, OMe, OH and NH₂. 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-NH₂ 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.

Proposal for a chemically consistent way to annotate ions arising from the analysis of reference compounds under ESI conditions: A prerequisite to proper mass spectral database constitution in metabolomics

Nowadays, high-resolution mass spectrometry is widely used for metabolomic studies. Thanks to its high sensitivity, it enables the detection of a large range of metabolites. In metabolomics, the continuous quest for a metabolite identification as complete and accurate as possible has led during the last decade to an ever increasing development of public MS databases (including LC-MS data) concomitantly with bioinformatic tool expansion. To facilitate the annotation process of MS profiles obtained from biological samples, but also to ease data sharing, exchange, and exploitation, the standardization and harmonization of the way to describe and annotate mass spectra seemed crucial to us. Indeed, under electrospray (ESI) conditions, a single metabolite does not produce a unique ion corresponding to its protonated or deprotonated form but could lead to a complex mixture of signals. These MS signals result from the existence of different natural isotopologues of the same compound and also to the potential formation of adduct ions, homomultimeric and heteromultimeric ions, fragment ions resulting from "prompt" in-source dissociations. As a joint reflection process within the French Infrastructure for Metabolomics and Fluxomics (MetaboHUB) and with the purpose of developing a robust and exchangeable annotated MS database made from pure reference compounds (chemical standards) analysis, it appeared to us that giving the metabolomics community some clues to standardize and unambiguously annotate each MS feature was a prerequisite to data entry and further efficient querying of the mass spectral database. The use of a harmonized notation is also mandatory for interlaboratory MS data exchange. Additionally, thorough description of the variety of MS signals arising from the analysis of a unique metabolite might provide greater confidence on its annotation.

Exploring the Role of Consecutive Addition of Nitrogen Atoms on Stability and Reactivity of Hydrogen-Bonded Azine-Water Complexes

The second-order Møller-Plesset perturbation theory (MP2) and density functional theory with dispersion function calculations have been applied to investigate the hydrogen-bonding interaction between azines and water. The study suggests that the ability of nitrogen present in azine to act as a hydrogen-bond acceptor decreases in the order of pyridine (PY) > diazine (DZ) > triazine (TZ) > tetrazine (TTZ) > pentazine (PZ) > hexazine (HZ). Natural bond orbital (NBO) analysis, atoms in molecules, symmetry-adapted perturbation theory (SAPT), and molecular electrostatic potential studies reflect the factors important for hydrogen-bond strength as well as for the structural, electronic, and vibrational changes occurring during complexation. NBO analysis reflects that upon gradual addition of nitrogen atoms, hyperconjugation leads to an increase in the population of antibonding O-H bond, thus causing elongation and weakening of O-H bond in complexes incorporating N···H-O_W interaction, whereas rehybridization leads to an increase in the s character of the carbon hybrid orbital in C-H bond, thus causing contraction and shortening of C-H bond in complexes having C-H···O_W interactions. From the topological analysis, an excellent linear correlation is found to exist between stabilization energy (ΔE _BSSE), electron density (ρ_c), and its Laplacian (∇²ρ_c) at the bond critical points.

Application of the Extended HOMED (Harmonic Oscillator Model of Aromaticity) Index to Simple and Tautomeric Five-Membered Heteroaromatic Cycles with C, N, O, P, and S Atoms

The geometry-based HOMA (Harmonic Oscillator Model of Aromaticity) descriptor, based on the reference compounds of different delocalizations of n- and π-electrons, can be applied to molecules possessing analogous bonds, e.g., only CC, only CN, only CO, etc. For compounds with different heteroatoms and a different number of CC, CX, XX, and XY bonds, its application leads to some discrepancies. For this reason, the structural descriptor was modified and the HOMED (Harmonic Oscillator Model of Electron Delocalization) index defined. In 2010, the HOMED index was parameterized for compounds with C, N and O atoms. For parametrization, the reference molecules of similar delocalizations of n- and π-electrons were employed. In this paper, the HOMED index was extended to compounds containing the CP, CS, NN, NP, PP, NO, NS, PO, and PS bonds. For geometrical optimization of all reference molecules and of all investigated heterocompounds, the same quantum–chemical method {B3LYP/6-311+G(d,p)} was used to eliminate errors of the HOMED estimation. For some tautomeric systems, the Gn methods were also employed to confirm tautomeric preferences. The extended HOMED index was applied to five-membered heterocycles, simple furan and thiophene, and their N and P derivatives as well as for tautomeric pyrrole and phosphole and their N and P derivatives. The effects of additional heteroatom(s) in the ring on the HOMED values for furan are parallel to those for thiophene. For pyrroles, aromaticity dictates the tautomeric preferences. An additional N atom in the ring only slightly affects the HOMED values for the favored and well delocalized NH tautomers. Significant changes take place for their rare CH forms. When intramolecular proton-transfer is considered for phosphole and its P derivatives, the PH tautomers seem to be favored only for 1,2,3-triphosphole/1,2,5-triphosphole and for 1,2,3,5-tetraphosphole. For other phospholes, the CH forms have smaller Gibbs energies than the PH isomers. For phosphazoles, the labile proton in the favored form is linked to the N atom. The PH forms have smaller HOMED indices than the NH tautomers but higher than the CH ones.

An effective and facile synthesis of new blue fluorophores on the basis of an 8-azapurine core

A convenient synthesis of 2-aryl-2,4-dihydro-5H-[1,2,3]triazolo[4,5-d]pyrimidin-5-ones (DTPs) from 3,3-diamino-2-(arylazo)acrylonitriles through a versatile and readily accessible two-step procedure is described. Density functional theory (DFT) calculations were performed to explain the selectivity of the heterocyclization step, which predominantly afforded 6-amino-5-(arylazo)pyrimidin-2(1H)-thiones in chloroform or ethanol, and 2,3-dihydro-1,2,4-triazines in toluene or DMF. Novel 2-aryl-2,4-dihydro-5H-[1,2,3]triazolo[4,5-d]pyrimidin-5-ones were obtained in good yields and showed absorption in the ultraviolet region and good emission in the blue region. The photophysical properties of DTPs were better than those cited in select literature examples of 8-azapurines. Owing to the facile synthesis and good photophysical characteristics in an aqueous medium, the new DTPs should have potential applications as organic fluorophores in fluorescence imaging and materials science.

Effects of Positive and Negative Ionization on Prototropy in Pyrimidine Bases: An Unusual Case of Isocytosine

Intramolecular proton-transfers (prototropic conversions) have been studied for the guanine building block isocytosine (iC), and effects of positive ionization, called one-electron oxidation (iC - e → iC^+•), and negative ionization, called one-electron reduction (iC + e → iC^-•), on tautomeric conversions when proceeding from neutral to ionized isocytosine have been discussed. Although radical cations and radical anions are very short-lived species, the ionization effects could be investigated by quantum-chemical methods. Such kind of studies gives some information about the labile protons and the most basic positions in the neutral and radical forms of the tautomeric system. For investigations, the complete isomeric mixture of isocytosine has been considered and calculations performed in two extreme environments, apolar {DFT(B3LYP)/6-311+G(d,p)} and polar {PCM(water)//DFT(B3LYP)/6-311+G(d,p)}. For selected isomers, the G4 theory has also been applied. There are no good relations for energetic parameters of neutral and ionized forms. Ionization energies depend on localization of labile protons. Tautomeric equilibria for neutral and ionized isocytosine, favored sites of protonation and deprotonation, and favored structures of protonated and deprotonated forms strongly depend on environment. Acidity of iC^+• is close to that of the iC conjugate acid, and basicity of iC^-• is close to that of the iC conjugate base. This increase of acid-base properties of charged radicals explains the proton-transfer in ionized pairs of nucleobases. When compared to other pyrimidine bases such as uracil (U) and cytosine (C), which exhibit analogous tautomeric equilibria between nine prototropic tautomers as isocytosine, the tautomeric preferences for iC, iC^+•, iC^-•, U, U^+•, U^-•, C, C^+•, and C^-• are completely different. The differences suggest that acid-base properties of functional groups, their stabilities, and ionization energies play a principal role in proton-transfers for pyrimidine bases and influence compositions of tautomeric mixtures.

A QM/QTAIM detailed look at the Watson-Crick↔wobble tautomeric transformations of the 2-aminopurine·pyrimidine mispairs

This work is devoted to the careful QM/QTAIM analysis of the evolution of the basic physico-chemical parameters along the intrinsic reaction coordinate (IRC) of the biologically important 2AP·T(WC)↔2AP·T*(w) and 2AP·C*(WC)↔2AP·C(w) Watson-Crick(WC)↔wobble(w) tautomeric transformations obtained at each point of the IRC using original authors' methodology. Established profiles reflect the high similarity between the courses of these processes. Basing on the scrupulous analysis of the profiles of their geometric and electron-topological parameters, it was established that the dipole-active WC↔w tautomerizations of the Watson-Crick-like 2AP·T(WC)/2AP·C*(WC) mispairs, stabilized by the two classical N3H⋯N1, N2H⋯O2 and one weak C6H⋯O4/N4 H-bonds, into the wobble 2AP·T*(w)/2AP·C(w) base pairs, respectively, joined by the two classical N2H⋯N3 and O4/N4H⋯N1 H-bonds, proceed via the concerted stepwise mechanism through the sequential intrapair proton transfer and subsequent large-scale shifting of the bases relative each other, through the planar, highly stable, zwitterionic transition states stabilized by the participation of the four H-bonds - N1⁺H⋯O4^-/N4^-, N1⁺H⋯N3^-, N2⁺H⋯N3^-, and N2⁺H⋯O2^-. Moreover, it was found out that the 2AP·T(WC)↔2AP·T*(w)/2AP·C*(WC)↔2AP·C(w) tautomerization reactions occur non-dissociatively and are accompanied by the consequent replacement of the 10 unique patterns of the specific intermolecular interactions along the IRC. Obtained data are of paramount importance in view of their possible application for the control and management of the proton transfer, e.g. by external electric or laser fields.

A QM/QTAIM research under the magnifying glass of the DPT tautomerisation of the wobble mispairs involving 2-aminopurine

In this study, a comprehensive survey of the changes of the physico-chemical parameters at each point of the IRC for the biologically important T·2AP*(w) ↔ T*·2AP(w) and G·2AP*(w) ↔ G*·2AP(w) DPT tautomerisation reactions involved in the point mutations (transitions and transversions) induced by 2-aminopurine (2AP) in DNA is provided.

Whether 2-aminopurine induces incorporation errors at the DNA replication? A quantum-mechanical answer on the actual biological issue

In this paper, we consider the mutagenic properties of the 2-aminopurine (2AP), which has intrigued molecular biologists, biophysicists and physical chemists for a long time and been widely studied by both experimentalists and theorists. We have shown for the first time using QM calculations, that 2AP very effectively produces incorporation errors binding with cytosine (C) into the wobble (w) C·2AP(w) mispair, which is supported by the N4H⋯N1 and N2H⋯N3 H-bonds and is tautomerized into the Watson-Crick (WC)-like base mispair C*·2AP(WC) (asterisk denotes the mutagenic tautomer of the base), that quite easily in the process of the thermal fluctuations acquires enzymatically competent conformation. 2AP less effectively produces transversions forming the wobble mispair with A base - A·2AP(w), stabilized by the participation of the N6H⋯N1 and N2H⋯N1 H-bonds, followed by further tautomerization A·2AP(w) → A*·2AP(WC) and subsequent conformational transition A*·2AP(WC) → A*·2AP_syn thus acquiring enzymatically competent structure. In this case, incorporation errors occur only in those case, when 2AP belongs to the incoming nucleotide. Thus, answering the question posed in the title of the article, we affirm for certain that 2AP induces incorporation errors at the DNA replication. Obtained results are consistent well with numerous experimental data.

Enhanced Basicity of Push-Pull Nitrogen Bases in the Gas Phase

Nitrogen bases containing one or more pushing amino-group(s) directly linked to a pulling cyano, imino, or phosphoimino group, as well as those in which the pushing and pulling moieties are separated by a conjugated spacer (C═X)_n, where X is CH or N, display an exceptionally strong basicity. The n-π conjugation between the pushing and pulling groups in such systems lowers the basicity of the pushing amino-group(s) and increases the basicity of the pulling cyano, imino, or phosphoimino group. In the gas phase, most of the so-called push-pull nitrogen bases exhibit a very high basicity. This paper presents an analysis of the exceptional gas-phase basicity, mostly in terms of experimental data, in relation with structure and conjugation of various subfamilies of push-pull nitrogen bases: nitriles, azoles, azines, amidines, guanidines, vinamidines, biguanides, and phosphazenes. The strong basicity of biomolecules containing a push-pull nitrogen substructure, such as bioamines, amino acids, and peptides containing push-pull side chains, nucleobases, and their nucleosides and nucleotides, is also analyzed. Progress and perspectives of experimental determinations of GBs and PAs of highly basic compounds, termed as "superbases", are presented and benchmarked on the basis of theoretical calculations on existing or hypothetical molecules.

Effects of ionization on stability of 1-methylcytosine - DFT and PCM studies

Consequences of ionization were studied by quantum-chemical methods (DFT and PCM) for 1-methylcytosine (MC)—a model of the nucleobase cytosine (C) connected with sugar in DNA. For calculations, three prototropic tautomers (one amino and two imino forms) and two imino zwitterions were considered, including conformational or configurational isomerism of exo heterogroups. Ionization and interactions between neighboring groups affect intramolecular proton-transfers, geometric and thermodynamic parameters, and electron delocalization for individual isomers. We discovered that an imino isomer is present in the isomeric mixture in the highest amount for positively ionized MC. Its contribution in neutral and negatively ionized MC is considerably smaller. Acid-base parameters for selected radical ions were estimated in the gas phase and compared to those of neutral MC. Gas-phase acidity of radical cations is close to that of the conjugate acid of MC, and gas-phase basicity of radical anions is close to that of the conjugate base of MC. Various routes of amino-imino conversion between neutral and ionized isomers were considered. Energetic-barrier for intramolecular proton-transfer in MC is close to that in the parent system—formamidine.

Whether the amino–imino tautomerism of 2-aminopurine is involved into its mutagenicity? Results of a thorough QM investigation

2AP* mutagenic tautomer is able to induce only one incorporation error – transversion – by pairing through the H-bonds into the G·2AP* mispair.

Structural grounds for the 2-aminopurine mutagenicity: a novel insight into the old problem of the replication errors

Mutagenic pressure of the 2AP molecule on DNA during its replication is realized via the more intensive generation of the T* mutagenic tautomers through the reaction 2AP·T(WC) → 2AP·T*(w).