A theoretical reinvestigation of the nucleic bases adenine, guanine, cytosine, thymine and uracil using AM1

The Enzymatic Decarboxylation Mechanism of 5-Carboxy Uracil: A Comprehensive Quantum Chemical Study

Dynamic regulation of DNA methylation is an important process for the control of gene expression in mammals. It is believed that in the demethylation pathway of 5-methyl cytosine, the intermediate 5-carboxy cytosine (5caC) can be actively decarboxylated alongside the substitution in the base excision repair. For the active decarboxylation of 5caC, a decarboxylase has not been identified so far. Due to the similar chemistry of the decarboxylation of 5-carboxy uracil (5caU) to uracil (U) in the pyrimidine salvage pathway catalyzed by the iso-orotate decarboxylase (IDCase), the study of this reaction might give valuable insights into the active 5caC decarboxylation process. In this work, we employ quantum chemical and molecular mechanic calculations and find that the catalytic mechanism of IDCase proceeds via a direct decarboxylation mechanism. Detailed investigations on the reaction coordinate reveal that it is a one-step mechanism with concerted proton transfer and C-C bond opening.

Exploring the antioxidant activity of thiaflavan compounds: a quantum chemical study

Conceptual DFT tools (HAT, SPLET, SET-PT, aromaticity index,…) have been used to explore the antioxidant activity of thiaflavan compounds, and predict which derivate should be the best one.

Tautomerism of uracil and related compounds: A mass spectrometry study

It has been demonstrated that uracil has a preponderant tautomeric form, but it is also known that different tautomers co-exist in this equilibrium. In this work, mass spectrometry is used as a helpful tool to analyse the equilibria, using derivative compounds to forbid the presence of some tautomers and ion trap mass spectrometry to follow relevant fragmentation pathways. Theoretical calculations were performed to confirm tautomers abundance by energy minimization in gas phase. Analysis of mass spectra of uracil, three methyl-substituted uracils, 2-thiouracil and three benzouracils suggest that uracil exists mainly as three tautomers in gas phase: one major structure that corresponds to the classical structure of uracil (pyrimidine-2,4(1H,3H)-dione) bearing two carbonyls and two NH moieties, and two minor enolic forms (4-hydroxypyrimidin-2(1H)-one and 2-hydroxypyrimidin-4(1H)-one). Such tautomeric distribution is supported by theoretical calculations, which show that they are the three most stable tautomers.

Vibrational spectra, tautomerism and thermodynamics of anticarcinogenic drug: 5-fluorouracil

The FT-IR and FT-Raman spectra of 5-Fluorouracil were recorded in the solid phase in the regions 400-4000 cm(-1) and 50-4000 cm(-1), respectively. The vibrational spectra were analysed and the observed fundamentals were assigned to different normal modes of vibration. The experimental wavenumbers were compared with the scaled vibrational values using DFT methods: the Ar matrix data were related to gas phase calculations, while the values of the solid state spectra were compared to those with dimer simulations. The study indicates that some features that are characteristic of vibrational spectra of uracil and its derivatives are retained in the spectrum of 5-fluorouracil and it exists in ketonic form in the solid phase. The tautomerism was also studied and the spectra of the two most stable forms were simulated. The calculated wavenumbers have been employed to yield thermodynamic properties.

Encapsulation of ruthenium nitrosylnitrate and DNA purines in nanostructured sol-gel silica matrices

The interactions between DNA purines (guanine and adenine) and the ruthenium complex Ru(NO)(NO(3))(3) were studied within nanostructured silica matrices prepared by a two-step sol-gel process. By infrared analysis in diffuse reflectance mode, it was proved that encapsulation induces a profound modification on the complex, whereas guanine and adenine preserve their structural integrity. The complex undergoes nitrate ligand exchange and co-condenses with the silica oligomers, but the nitrosyl groups remain stable, which is an unusual behavior in Ru nitrosyl complexes. In turn, the doping molecules affect the sol-gel reactions and eventually the silica structure as it forms: the complex yields a microporous structure, and the purine bases are responsible for the creation of macropores due to hydrogen bonding with the silanol groups of the matrix. In a confined environment, the interactions are much stronger for the coencapsulated pair guanine complex. While adenine only establishes hydrogen bonds or van der Waals interactions with the complex, guanine bonds covalently to Ru by one N atom of the imidazole ring, which becomes strongly perturbed, resulting in a deformation of the complex geometry.

The electronic structure of the four nucleotide bases in DNA, of their stacks, and of their homopolynucleotides in the absence and presence of water

Using the ab initio Hartree-Fock crystal orbital method in its linear combination of atomic orbital form, the energy band structure of the four homo-DNA-base stacks and those of poly(adenilic acid), polythymidine, and polycytidine were calculated both in the absence and presence of their surrounding water molecules. For these computations Clementi's double zeta basis set was applied. To facilitate the interpretation of the results, the calculations were supplemented by the calculations of the six narrow bands above the conduction band of poly(guanilic acid) with water. Further, the sugar-phosphate chain as well as the water structures around poly(adenilic acid) and polythymidine, respectively, were computed. Three important features have emerged from these calculations. (1) The nonbase-type or water-type bands in the fundamental gap are all close to the corresponding conduction bands. (2) The very broad conduction band (1.70 eV) of the guanine stack is split off to seven narrow bands in the case of poly(guanilic acid) (both without and with water) showing that in the energy range of the originally guanine-stack-type conduction band, states belonging to the sugar, to PO(4)(-), to Na(+), and to water mix with the guanine-type states. (3) It is apparent that at the homopolynucleotides with water in three cases the valence bands are very similar (polycytidine, because it has a very narrow valence band, does not fall into this category). We have supplemented these calculations by the computation of correlation effects on the band structures of the base stacks by solving the inverse Dyson equation in its diagonal approximation taken for the self-energy the MP2 many body perturbation theory expression. In all cases the too large fundamental gap decreased by 2-3 eV. In most cases the widths of the valence and conduction bands, respectively, decreased (but not in all cases). This unusual behavior is most probably due to the rather large complexity of the systems. From all this emerges the following picture for the charge transport in DNA: There is a possibility in short segments of the DNA helix of a Bloch-type conduction of holes through the nucleotide base stacks of DNA combined with hopping (and in a lesser degree with tunneling). The motivation of this large scale computation was that recently in Zurich (ETH) they have performed high resolution x-ray diffraction experiments on the structure of the nucleosomes. The 8 nucleohistones in them are wrapped around by a DNA superhelix of 147 base pairs in the DNA B form. The most recent investigations have shown that between the DNA superhelix (mostly from its PO(4) (-) groups) there is a charge transfer to the positively charged side chains (first of all arginines and lysines) of the histones at 120 sites of the superhelix. This would cause a hole conduction in DNA and an electronic one in the proteins.

((15)N(5))-labeled adenine derivatives: synthesis and studies of tautomerism by (15)N NMR spectroscopy and theoretical calculations

Since the nitrogens of nucleosides and nucleotides play an important role in the molecular recognition of these compounds, (15)N NMR became a method of choice in this field. Fully (15)N-labeled adenine, required in the latter studies, was obtained in four synthetic steps, in a good yield. Likewise, ((15)N(5))-2-hexylthioether-adenine and ((15)N(5))-8-Br-adenine were obtained in five synthetic steps from the relatively inexpensive (15)N sources: (15)N-NH(4)Cl, (15)N-NH(4)OH, (15)N-NaNO(2). Full (15)N labeling of these adenine prototypes enabled to obtain high-resolution (15)N NMR spectra of these bases at 60.8 MHz. Furthermore, the spectra suggested the existence of the N3-H species in the tautomeric mixtures of these compounds in solution, in addition to the well-reported N9-H (major) and N7-H (minor) tautomers. These observations were also supported by quantum mechanical calculations of the tautomeric equilibria in the gas phase and in solution of the above-mentioned adenine compounds. The gas-phase tautomeric equilibria were estimated using density functional theory and second-order perturbation theory methods. Solvent effects were included by means of both continuum and discrete solvation models. The observation of the existence of the N3-H tautomer has a clear impact on the possible H-bonding patterns of these adenine prototypes and on their molecular recognition by various biological macromolecules. The above(15)N-labeled analogues are expected to find use as (15)N NMR probes for numerous biochemical studies.

Proton affinity of uracil. A computational study of protonation sites

Relative stabilities of uracil tautomers and cations formed by gas-phase protonation were studied computationally with the B3LYP, MP2, QCISD, and QCISD(T) methods and with basis sets expanding from 6-31G(d,p) to 6-311+G(3df,2p). In accordance with a previous density functional theory study, the dioxo tautomer 1a was the most stable uracil isomer in the gas phase. Gibbs free energy calculations using effective QCISD(T)/6-311+G(3df,2p) energies suggested >99.9% of 1a at equilibrium at 523 K. The most stable ion isomer corresponded to N-1 protonated 2,4-dihydroxypyrimidine, which however is not formed by direct protonation of 1a. The topical proton affinities in 1a followed the order O-8 > O-7 > C-5 > N-3 > N-1. The thermodynamic proton affinity of 1a was calculated as 858 kJ mol(-1) at 298 K. A revision is suggested for the current estimate included in the ion thermochemistry database.

Density functional theory studies on tautomeric stability and infrared spectra of 2-chloroadenine

The tautomeric stability and vibrational IR spectrum of 2-chloroadenine were studied using density functional theory (DFT) at B3LYP/6-31G* level. The amino N(9)H tautomer of 2-chloroadenine was predicted to be most stable. A scaled quantum mechanical (SQM) force field approach was used to calculate the vibrational frequencies of amino N(9)H form of 2-chloroadenine. The force constant scale factors were transferred from those of purine. The mean deviation between the predicted vibrational frequencies and the observed ones is 8.0 cm(-1). The results made it possible to give complete assignments of the IR spectrum of this molecule.

Molecular recognition of modified adenine nucleotides by the P2Y(1)-receptor. 2. A computational approach

The molecular recognition of C2- or C8-substituted ATP derivatives by the P2Y(1)-receptor (P2Y(1)-R) is analyzed using ab initio quantum mechanical calculations. Parameters that may determine ligand specificity toward P2Y(1)-R were examined on reduced models and correlated with the biochemical data for the parent compounds. These include tautomerism and protonation energy in the gas and aqueous phases, as well as molecular electrostatic potential (MEP) and dipole moment vector. The calculated electronic parameters cannot explain the inactivity of the C8-substituted ATP derivatives, nor the difference in activity among the C2-substituted ATP analogues. These results indicate that neither tautomerism nor changes in the electronic distribution of the adenine ring play a major role in determining binding specificity of adenine nucleotides to the receptor. It is suggested that the higher potency of the C2-substituted ATP derivatives, compared to ATP, might be due to interaction between the C2 side chain heteroatom and the receptor. Furthermore, the interaction of the C2 alkyl side chain with a hydrophobic pocket at the receptor binding site is suggested. In addition, NMR data in the companion paper indicate that the inactivity of the C8-substituted ATP analogues may be due to steric and conformational, rather than electronic, effects.

Elucidation of primary radiation damage in DNA through application of ab initio molecular orbital theory

This review summarizes the results of recent ab initio molecular orbital calculations performed on DNA constituents that attempt to further our understanding of the radiation-induced damage to DNA. The results reviewed include calculations performed on the four individual DNA bases, the base pairs in gas phase and modelled aqueous phase, the deoxyribose moiety, and a portion of the sugar-phosphate backbone. The emphasis is on the electron affinities and ionization potentials of the radical species calculated under various conditions (i.e. gas phase, aqueous phase, proton transfer, base stacking), as it has been shown that the initial ion radical distribution is largely a function of these two properties. Theoretical studies of the electronic excited states of the individual bases and radioprotection of the biomolecule by various thiol compounds are also reviewed. Finally, a summary is provided to allow for further elaboration of the current model for radiation damage to DNA and to show the present advantages and limitations of ab initio theory in the investigation of such processes.