Effect of the microhydration on the tautomerism in the anticarcinogenic drug 5-fluorouracil and relationships with other 5-haloderivatives

On the aromaticity of uracil and its 5-halogeno derivatives as revealed by theoretically derived geometric and magnetic indexes

The problem of aromaticity in heterocyclic rings of uracil and its 5-halogenoderivatives (5XU) was analyzed theoretically by calculating modified harmonic oscillator model of aromaticity (HOMA) for Heterocycle Electron Delocalization (HOMHED), nucleus-independent chemical shift parameters (NICS) and the so-called scan experiments, using helium-3 atom as a magnetic probe. The impact of halogen electronegativity on C5 atom’s NBO charges was also investigated. Water, as a polar environment, has a negligible impact on 5XU aromaticity. The most stable diketo tautomer shows a very low aromaticity while the “rare” dihydroxy form (tautomer No 6) is aromatic and resembles benzene. This is in agreement with traditional drawing of chemical formula of uracil’s six-membered ring, directly showing three alternating single and double bonds in its tautomer No 6. No good correlation between magnetic and geometric indexes of aromaticity for the studied 5XU tautomers was found. Linear correlation between the magnitude of NICS minimum, as well as the distance of the minimum above uracil ring plane center from 3He NMR chemical shift scan plot with respect to halogen electronegativity were observed. A strong linear dependence of magnetic index of aromaticity and the electronegativity of 5X substituent was observed.

Factors Governing the Chemical Stability and NMR Parameters of Uracil Tautomers and Its 5-Halogen Derivatives

We report on the density functional theory (DFT) modelling of structural, energetic and NMR parameters of uracil and its derivatives (5-halogenouracil (5XU), X = F, Cl, Br and I) in vacuum and in water using the polarizable continuum model (PCM) and the solvent model density (SMD) approach. On the basis of the obtained results, we conclude that the intramolecular electrostatic interactions are the main factors governing the stability of the six tautomeric forms of uracil and 5XU. Two indices of aromaticity, the harmonic oscillator model of aromaticity (HOMA), satisfying the geometric criterion, and the nuclear independent chemical shift (NICS), were applied to evaluate the aromaticity of uracil and its derivatives in the gas phase and water. The values of these parameters showed that the most stable tautomer is the least aromatic. A good performance of newly designed xOPBE density functional in combination with both large aug-cc-pVQZ and small STO(1M)-3G basis sets for predicting chemical shifts of uracil and 5-fluorouracil in vacuum and water was observed. As a practical alternative for calculating the chemical shifts of challenging heterocyclic compounds, we also propose B3LYP calculations with small STO(1M)-3G basis set. The indirect spin-spin coupling constants predicted by B3LYP/aug-cc-pVQZ(mixed) method reproduce the experimental data for uracil and 5-fluorouracil well.

Biomolecules of 2-Thiouracil, 4-Thiouracil and 2,4-Dithiouracil: A DFT Study of the Hydration, Molecular Docking and Effect in DNA:RNAMicrohelixes

The molecular structure of 2-thiouracil, 4-thiouracil and 2,4-dithiouracil was analyzed under the effect of the first and second hydration shell by using the B3LYP density functional (DFT) method, and the results were compared to those obtained for the uracil molecule. A slight difference in the water distribution appears in these molecules. On the hydration of these molecules several trends in bond lengths and atomic charges were established. The ring in uracil molecule appears easier to be deformed and adapted to different environments as compared to that when it is thio-substituted. Molecular docking calculations of 2-thiouracil against three different pathogens: Bacillus subtilis, Escherichia coli and Candida albicans were carried out. Docking calculations of 2,4-dithiouracil ligand with various targeted proteins were also performed. Different DNA: RNA hybrid microhelixes with uridine, 2-thiouridine, 4-thiouridine and 2,4-dithiouridine nucleosides were optimized in a simple model with three nucleotide base pairs. Two main types of microhelixes were analyzed in detail depending on the intramolecular H-bond of the 2'-OH group. The weaker Watson-Crick (WC) base pair formed with thio-substituted uracil than with unsubstituted ones slightly deforms the helical and backbone parameters, especially with 2,4-dithiouridine. However, the thio-substitution significantly increases the dipole moment of the A-type microhelixes, as well as the rise and propeller twist parameters.

Structure and conformational analysis of the anti-HIV reverse transcriptase inhibitor AZT using MP2 and DFT methods. Differences with the natural nucleoside thymidine. Simulation of the 1st phosphorylation step with ATP

A comprehensive quantum-chemical investigation of the conformational landscape of the HIV-1 reverse transcriptase inhibitor AZT (3'-azido-3'-deoxythymidine) nucleoside analogue was carried out. The whole conformational parameters (χ, γ, β, δ, ϕ, P, νmax) were analysed as well as the NBO charges. The search located at least 55 stable structures, 9 of which were by MP2 within a 1 kcal mol(-1) electronic energy range of the global minimum. Most conformers were anti or high-anti around the glycoside bond and with North sugar ring puckering angles. The distribution of all the conformers according to the ranges of stability of the characteristic torsional angles was established. The results obtained were in accordance with those found in related anti-HIV nucleoside analogues. The best conformer in the anti form corresponded to the calculated values by MP2 of χ = -126.9°, β = 176.4° and γ = 49.1°. An analysis of the lowest vibrations in conformer C1 was carried out. The first hydration shell was simulated and the structural differences with the natural nucleoside deoxythymidine (dT) were determined. The first phosphorylation step was simulated by interacting ATP with the best hydrated clusters of AZT and dT. The Na cations act as a bridge between the phosphate moieties of ATP making it easy for -P3O3 to receive the H5' proton from AZT or dT. A proton-transfer mechanism is proposed through the water molecules. When the number of the water molecules surrounding AZT is lower than 8, the first phosphorylation step of AZT can be carried out. However, the appropriate orientation of the O5'-H in dT avoids this limitation and it can be performed with large numbers of water molecules.

FT-IR and FT-Raman spectra of 5-fluoroorotic acid with solid state simulation by DFT methods

FT-Raman and FT-IR studies of the biomolecule 5-fluoroorotic acid in the solid state were carried out. The unit cell found in the crystal was simulated as a tetramer form by density functional calculations. They were performed to clarify wavenumber assignments of the experimental observed bands in the spectra. Correlations with the molecule of uracil were made, and specific scale equations were employed to scale the wavenumbers of 5-fluoroorotic acid. Good reproduction of the experimental wavenumbers is obtained and the % error is very small in the majority of the bands. This fact confirms our simplified solid state model. The molecular structure was fully optimized using DFT and MP2 methods. The relative stability of both the syn and anti conformations was investigated, and the anti-form was found to be slightly more stable, by 7.49 kJ/mol at the MP2 level. The structures of all possible tautomeric forms were determined. The keto-form appeared as the most stable one. The NBO atomic charges and several thermodynamic parameters were also calculated.

FT-IR and FT-Raman spectra of 6-chlorouracil: molecular structure, tautomerism and solid state simulation. A comparison between 5-chlorouracil and 6-chlorouracil

A Raman and IR study of the biomolecule 6-chlorouracil was carried out in the solid state. The unit cell found in the crystal was simulated as a tetramer form by density functional calculations. Specific scale factors and scaling equations deduced from uracil molecule were employed in the predicted wavenumbers of 6-chlorouracil. The scaled wavenumbers were used in the reassignment of the IR and Raman experimental bands. Good reproduction of the experimental wavenumbers is obtained and the % error is very small in the majority of cases. A comparison between the molecular structure and charge distribution of 6-chlorouracil and 5-chlorouracil molecules was presented. The effect of the hydration with the PCM model in the molecular structure and charges was discussed. The optimum tautomers of 6-chlorouracil were optimized and analyzed. Six of them were related to those of uracil molecule. The effect of the halogen substitution in the sixth position of the pyrimidine ring in the stability of the different tautomers was evaluated. HOMO and LUMO orbital energy analysis were carried out.