Hydration structure of uracil as studied by 1D and 2D heteronuclear Overhauser spectroscopy: evidence for the formation of a trihydrate in the first solvation shell

How proton transfer affects the helical parameters in DNA:DNA microhelices

Proton transfer reactions are a widespread phenomenon in many areas of the life sciences and it is one of the origins of the spontaneous point mutations during DNA replication. Because of its importance, many studies have been reported on these reactions. However, the present work is the first one focused on the structural geometrical changes by double proton transfer (DPT). Thus, different Watson-Crick (WC) pairs were optimized first in a simple model with one nucleoside base pair, and in a microhelix form with three nucleoside base pairs. The canonical and few tautomeric forms were considered in DNA:DNA microhelices with A-type and B-type helical forms. The stability of these structures and how the DPT process affects the main geometrical parameters was analyzed, in particular the deformation of the helical parameters. The M06-2X DFT method was used for this purpose. The purine/pyrimidine ring in the keto form appears easier to be deformed than when it is in the enol form. The weaker WC base pair formed with mixed microhelices than with nucleobases alone and the significant deformation of the helical and backbone parameters with the DPT appears to complicate this process in microhelices.Communicated by Ramaswamy H. Sarma.

Density Functional Method Study on the Cooperativity of Intermolecular H-bonding and π-π+ Stacking Interactions in Thymine-[Cnmim]Br (n = 2, 4, 6, 8, 10) Microhydrates

The exploration of the ionic liquids’ mechanism of action on nucleobase’s structure and properties is still limited. In this work, the binding model of the 1-alkyl-3-methylimidazolium bromide ([Cnmim]Br, n = 2, 4, 6, 8, 10) ionic liquids to the thymine (T) was studied in a water environment (PCM) and a microhydrated surroundings (PCM + wH2O). Geometries of the mono-, di-, tri-, and tetra-ionic thymine (T-wH2O-y[Cnmim]+-xBr−, w = 5~1 and x + y = 0~4) complexes were optimized at the M06-2X/6-311++G(2d, p) level. The IR and UV-Vis spectra, QTAIM, and NBO analysis for the most stable T-4H2O-Br−-1, T-3H2O-[Cnmim]+-Br−-1, T-2H2O-[Cnmim]+-2Br−-1, and T-1H2O-2[Cnmim]+-2Br−-1 hydrates were presented in great detail. The results show that the order of the arrangement stability of thymine with the cations (T-[Cnmim]+) by PCM is stacking > perpendicular > coplanar, and with the anion (T-Br−) is front > top. The stability order for the different microhydrates is following T-5H2O-1 < T-4H2O-Br−-1 < T-3H2O-[Cnmim]+-Br−-1 < T-2H2O-[Cnmim]+-2Br−-1 < T-1H2O-2[Cnmim]+-2Br−-1. A good linear relationship between binding EB values and the increasing number (x + y) of ions has been found, which indicates that the cooperativity of interactions for the H-bonding and π-π+ stacking is varying incrementally in the growing ionic clusters. The stacking model between thymine and [Cnmim]+ cations is accompanied by weaker hydrogen bonds which are always much less favorable than those in T-xBr− complexes; the same trend holds when the clusters in size grow and the length of alkyl chains in the imidazolium cations increase. QTAIM and NBO analytical methods support the existence of mutually reinforcing hydrogen bonds and π-π cooperativity in the systems.

Solvation of Uracil and Its Derivatives by DMSO: A DFT-Supported 1H NMR and 13C NMR Study

¹H NMR and ¹³C NMR spectra of uracil, thymine, 5-hydroxymethyluracil, 5,6-dihydrouracil, and 5,6-dihydrothymine in DMSO-d₆ solutions have been measured. Additionally, molecular structures as well as NMR parameters of these compounds and their various solvates have been calculated using DFT B3LYP/6-311++G(2d,p) PCM(DMSO) method. The analysis of the chemical shift data for these compounds has shown that, indeed, in DMSO solutions they occur as equilibrium mixtures of free molecules and solvates in which solute and solvent molecules are joined by NH···O or OH···O hydrogen bonds. The populations of particular species present in the solutions have been estimated. Moreover, it has been found that 5,6-dihydrothymine exists in DMSO solution preferentially in conformation with the methyl group occupying the pseudoequatorial position. This finding is based on the molecular energy calculations and remains in full agreement with the interpretation of NMR data and theoretical calculations of NMR parameters.

Optical absorption and magnetic circular dichroism spectra of thiouracils: a quantum mechanical study in solution

The excited electronic states of thiouracils, the analogues of uracil where the carbonyl oxygens are substituted by sulphur atoms, have been investigated by computing the magnetic circular dichroism (MCD) and one-photon absorption (OPA) spectra at the TD-DFT level of theory.

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.

Predicting the NMR spectra of nucleotides by DFT calculations: cyclic uridine monophosphate

We present an experimental and quantum chemical NMR study of the mononucleotide cyclic uridine monophosphate in water. Spectral parameters ((1)H and (13)C chemical shifts and (1)H--(1)H, (13)C--(1)H, (31)P--(13)C and (31)P--(1)H spin-spin coupling constants) have been carefully obtained experimentally and calculated using DFT methods including the solvent effect and the conformational flexibility of the solute. This study confirms that the (1)H and (13)C spectra of polar, flexible molecules in aqueous solution can be predicted with a high level of accuracy, comparable to that obtained for less complex systems.

Tautomerism of Uracil and 5-Bromouracil in a Microcosmic Environment with Water and Metal Ions. What Roles Do Metal Ions Play?

The base tautomerization processes of uracil/5-bromouracil were investigated in a microcosmic environment with both H2O and Na+ (W-M environment). It was found that uracil was more stable in the W-M environment than in the microcosmic environment with only water, which suggested that the metal ions and water work cooperated to maintain the classical nucleic acid bases. However, 5-bromouracil, a chemical mutagen, was found to be less stable than uracil in the W-M environment. Why the 5-bromouracil is easier to tautomerize and therefore induce gene mutation was explained to some extent. Further research revealed that the water molecule would assist the tautomerization in the W-M environment. However, the metal ions in different regions play absolutely opposite roles: in one region, the metal ions can prevent the base from tautomerizing, whereas in another region, the metal ion can assist the tautomerization process. Furthermore, from the viewpoint of ionization of the base, it seems BrU has a stronger tendency to lose the proton at N3, which is an intrinsic consequence of the bromine atom and is not affected by the metal cation.

Microsolvation Pattern of the Hydrated Radical Anion of Uracil: U-(H2O)n (n = 3−5)

The microsolvation patterns of the uracil radical anion in water clusters U-(H2O)n with n ranging from 3 to 5 were investigated by the density functional theory approach. The electron detachment energies (VDE) of the stable anionic complexes with different numbers of hydration water are predicted. The linear dependence of the VDE value of the most stable anionic complexes with respect to the hydration number suggests the importance of the clustered waters in the microsolvation of the radical anion of the nucleobases. The formation of the water clusters is found to be necessary in the most stable conformers of the tri-, tetra-, and pentahydrated radical anion of uracil. The microsolvation pattern with three or more well-separated hydration water molecules in the first hydration layer is less stable than the arrangement with the waters in tight clusters. The charge transfer between the anionic uracil and the hydration water is high. Good agreement between the experimental and the theoretical vertical detachment energy yield in this study further demonstrates the practicability of the B3LYP/DZP++ approach in the study of radical anions of the DNA subunits.

Combined use of Overhauser spectroscopy and NMR diffusion experiments for mapping the hydration structure of nucleosides: structure and dynamics of uridine in water

Complementary results from 13C intermolecular nuclear Overhauser effects (NOE), 1H-13C heteronuclear Overhauser spectroscopy (HOSEY) and 1H-NMR diffusion measurements were used for probing the structure of the first solvation shell of uridine in water. It is demonstrated that a cyclic dihydrate is formed. The two water molecules produce two hydrogen bonds with the two oxygen atoms from the pyrimidine ring and accept only one hydrogen bond from the amide proton. The dihydrate has only a short lifetime as compared with the rotational correlation time of the free nucleoside. The chemical exchange constant of the amide proton with water is then estimated by diffusion experiments. The results are consistent with previous data obtained for uracil in water and provide interesting information about water accessibility in nucleic acid bases.

Solvent Effect on the Singlet Excited-State Lifetimes of Nucleic Acid Bases: A Computational Study of 5-Fluorouracil and Uracil in Acetonitrile and Water

The first comprehensive quantum mechanical study of solvent effects on the behavior of the two lowest energy excited states of uracil derivatives is presented. The absorption and emission spectra of uracil and 5-fluorouracil in acetonitrile and aqueous solution have been computed at the time-dependent density-functional theory level, using the polarizable continuum model (PCM) to take into account bulk solvent effects. The computed spectra and the solvent shifts provided by our method are close to their experimental counterpart. The S0/S1 conical intersection, located in the presence of hydrogen-bonded solvent molecules by CASSCF (8/8) calculations, indicates that the mechanism of ground-state recovery, involving out-of-plane motion of the 5 substituent, does not depend on the nature of the solvent. Extensive explorations of the excited-state surfaces in the Franck-Condon (FC) region show that solvent can modulate the accessibility of an additional decay channel, involving a dark n/pi* excited state. This finding provides the first unifying explanation for the experimental trend of 5-fluorouracil excited-state lifetime in different solvents. The microscopic mechanisms underlying solvent effects on the excited-state behavior of nucleobases are discussed.

Solute-solvent contact by intermolecular cross relaxation. I. The nature of the water-hydrophobic interface

Intermolecular cross-relaxation rates between solute and solvent were measured by {1H} 19F nuclear magnetic resonance experiments in aqueous molecular solutions of ammonium perfluoro-octanoate and sodium trifluoroacetate. The experiments performed at three different magnetic fields provide frequency-dependent cross-relaxation rates which demonstrate clearly the lack of extreme narrowing for nuclear spin relaxation by diffusionally modulated intermolecular interactions. Supplemented by suitable intramolecular cross-relaxation, longitudinal relaxation, and self-diffusion data, the obtained cross-relaxation rates are evaluated within the framework of recent relaxation models and provide information about the hydrophobic hydration. In particular, water dynamics around the trifluoromethyl group in ammonium perfluoro-octanoate are more retarded than that in the smaller trifluoroacetate.

Determining chemical exchange rates of the uracil labile protons by NMR diffusion experiments

The exchange rates of the amido-protons of uracil with water were determined by NMR diffusion experiments and the results showed a factor 2 difference in lability between them, which was confirmed by more classical 2D-NMR exchange experiments.

Singlet Excited-State Behavior of Uracil and Thymine in Aqueous Solution: A Combined Experimental and Computational Study of 11 Uracil Derivatives

The excited-state properties of uracil, thymine, and nine other derivatives of uracil have been studied by steady-state and time-resolved spectroscopy. The excited-state lifetimes were measured using femtosecond fluorescence upconversion in the UV. The absorption and emission spectra of five representative compounds have been computed at the TD-DFT level, using the PBE0 exchange-correlation functional for ground- and excited-state geometry optimization and the Polarizable Continuum Model (PCM) to simulate the aqueous solution. The calculated spectra are in good agreement with the experimental ones. Experiments show that the excited-state lifetimes of all the compounds examined are dominated by an ultrafast (<100 fs) component. Only 5-substituted compounds show more complex behavior than uracil, exhibiting longer excited-state lifetimes and biexponential fluorescence decays. The S(0)/S(1) conical intersection, located at CASSCF (8/8) level, is indeed characterized by pyramidalization and out of plane motion of the substituents on the C5 atom. A thorough analysis of the excited-state Potential Energy Surfaces, performed at the PCM/TD-DFT(PBE0) level in aqueous solution, shows that the energy barrier separating the local S(1) minimum from the conical intersection increases going from uracil through thymine to 5-fluorouracil, in agreement with the ordering of the experimental excited-state lifetime.

NMR diffusion and nuclear Overhauser investigation of the hydration properties of thymine: influence of the methyl group

The absence of preferential hydration in thymine and its lowest water accessibility with respect to uracil were evidenced by NMR diffusion and HOESY experiments; the hydration differences observed between these pyrimidine bases were attributed to the electronic rather than steric properties of the methyl group.

The Influence of Microhydration on the Ionization Energy Thresholds of Uracil and Thymine

In the present study the ionization energy thresholds (IET's) of uracil and thymine have been calculated (with the B3LYP, PMP2, and P3 levels of theory using the standard 6-31++G(d,p) basis set) with one to three water molecules placed in the first hydration shell. Then (B3LYP) polarizable continuum model (PCM) calculations were performed with one to three waters of the hydration shell included. Calculations show there is a distinct effect of microhydration on uracil and thymine. For uracil, one added water results in a decrease in the IET of about 0.15 eV. The second and third water molecules cause a further decrease by about 0.07 eV each. For thymine, the first water molecule is seen to decrease the IET by about 0.1 eV, while the second and third water molecules cause a further decrease of less than 0.1 eV each. The changes in IET calculated here for thymine with one to three waters of hydration are smaller than the experimental values determined by Kim et al. (Kim, S. K.; Lee, W.; Herschbach, D. R. J. Phys. Chem. 1996, 100, 7933). Preliminary results presented here indicate that the experimental results may involve keto-enol tautomers of thymine. The results of placing the microhydrated structures of uracil and thymine in a PCM cavity was seen to make very little difference in the IET when compared to the IET of ordinary uracil or thymine in a PCM cavity. The implications are that accurate calculations of the IET's of uracil and thymine can be obtained by simply considering long-range solvation effects.

Effect of Hydration on the Lowest Singlet ππ* Excited-State Geometry of Guanine: A Theoretical Study

An ab-initio computational study was performed to investigate the effect of explicit hydration on the ground and lowest singlet PiPi* excited-state geometry and on the selected stretching vibrational frequencies corresponding to the different NH sites of the guanine acting as hydrogen-bond donors. The studied systems consisted of guanine interacting with one, three, five, six, and seven water molecules. Ground-state geometries were optimized at the HF level, while excited-state geometries were optimized at the CIS level. The 6-311G(d,p) basis set was used in all calculations. The nature of potential energy surfaces was ascertained via the harmonic vibrational frequency analysis; all structures were found minima at the respective potential energy surfaces. The changes in the geometry and the stretching vibrational frequencies of hydrogen-bond-donating sites of the guanine in the ground and excited state consequent to the hydration are discussed. It was found that the first solvation shell of the guanine can accommodate up to six water molecules. The addition of the another water molecule distorts the hydrogen-bonding network by displacing other neighboring water molecules away from the guanine plane.

Absorption and Fluorescence Spectra of Uracil in the Gas Phase and in Aqueous Solution: A TD-DFT Quantum Mechanical Study

Here we present the first computations of fluorescence spectra in aqueous solution at an accurate quantum mechanical level. From a methodological point of view, our study shows that by only taking into account both bulk effects and explicit solvent molecules it is possible to reproduce solvent effects on the energy and the intensities of the electronic spectra, especially for what concerns pi/pi* transition. The computed absorption and fluorescence spectra are in a good agreement with the available experimental results. The energy ordering between the lowest energy n-pi* and the pi/pi* transitions in uracil strongly depends on the nature of the embedding medium. The geometry of the first solvation shell is remarkably sensitive to the specific electronic state, suggesting that solvent degrees of freedom can act as S1/S2 coupling modes.

Intermolecular Overhauser effects in fluoroalcohol solutions of cyclo-alanylglycine

Interactions between the diketopiperazine cyclo-alanylglycine and four fluorinated alcohols in water-fluoroalcohol mixtures were examined by (1)H[(19)F] intermolecular nuclear Overhauser effects (NOE) experiments. The alcohols studied were trifluoroethanol, hexafluoroacetone trihydrate, 1,1,1,3,3,3-hexafluoroisopropanol and perfluoro-t-butanol. The experimental methods used permit detection of solvent-solute NOEs of 0.1% or less. Solute and solvent diffusion coefficients were determined and apparent molecular radii of the fluoroalcohols estimated. Using these data, it was shown that observed (1)H[(19)F] intermolecular NOEs are consistent with expectations based on theory. A method for extending conventional theory to take into account the shape of a solute and the exposure of its hydrogens to solvent is described. This approach gives reasonable agreement with experimental results, particularly if it is assumed that solute-solvent interactions take place in such a way that the fluorines of a fluoroalcohol are preferentially oriented toward the solute during solute-solvent encounters. The results support the suggestion that intermolecular (1)H[(19)F] NOEs may become a useful tool for studies of peptide and protein conformations in fluoroalcohol-water solvent mixtures.