Tautomerism in cytidine

Experimental and theoretical study of the cytosine tautomerism through excited states

CONTEXT

The irradiation of water solution of cytosine with UV light (λmax = 254 nm) shows oxo-hydroxy tautomerism with a rate constant of 6.297 × 10-3 min-1. The order of the reaction implies a tautomeric conversion. After removing the UV light source, we observed a dark reaction with a rate constant of 1.473 × 10-3 min-1 which leads to a restoration of the initial tautomer as before the irradiation. The mechanism of oxo-hydroxy tautomerism of cytosine in water solution was studied in the excited state. It was found that the transformations occur along the 1πσ* excited-state reaction paths which link the Franck-Condon geometries of the tautomers and the conical intersections S0/S1 connected with the H-detachment processes of the corresponding bonds. Furthermore, we established that the conical intersections S0/S1 are also mutually accessible along the 1πσ* excited-state reaction paths.

METHODS

The ground-state equilibrium geometries were optimized at the B3LYP/aug-cc-pVDZ level of theory in water environment according to PCM as well as at the CC2/aug-cc-pVDZ level in the gas phase. The TD B3LYP and CC2 methods were applied for the study of the excited states. The tautomerization mechanisms were studied with the use of the linear interpolation in internal coordinates approach using the optimized geometries of tautomers minima and conical intersections S0/S1 at the CASSCF(6,6)/6-31G* level. All calculations were performed with the GAUSSIAN 16 commercial software.

Electronic Spectra, Excited State Structures and Interactions of Nucleic Acid Bases and Base Assemblies: A Review

A comprehensive review of recent theoretical and experimental advances in the singlet electronic transitions, excited state structures and dynamics of nucleic acid bases (NABs) and base assemblies are presented. It is well known that NABs absorb ultraviolet radiation, but the absorbed energy is efficiently dissipated in the form of ultrafast internal conversion processes believed to occur in the subpicosecond time scale and, therefore, enabling NABs highly photostable. It is not known how much evolutionary role was played in evolving these molecules and the ultimate selection by nature as genetic materials, but it is well accepted that survival-of-fittest prevails. Recently, significant efforts have been continuously paid to understand the mechanism of electronic excitation deactivation, but universally acceptable mechanism is still elusive. However, recent investigations reveal that electronic excited state geometries of DNA bases are usually nonplanar and this structural nonplanarity may facilitate nonradiative deactivation. Investigation of excited state structures is challenging and, therefore, it is not surprising that despite the impressive theoretical and computational advances, this research area is still hampered by the methodological and computational limitations. Further, stacking has significant influence on the emission properties of molecules. The 2-aminopurine, a fluorescent adenine derivative frequently used in studying DNA dynamics, shows significant attenuations in fluorescence quantum yield when incorporated in the DNA. Theoretical and computational bottlenecks limit a thorough theoretical understanding of effect of stacking interactions on the excited state dynamics of NABs. Despite these limitations the investigations of excited state properties are progressing in the right direction and our better understanding of excited state structure and dynamics of NABs and nucleic acids may help to design preventive strategy for radiation induced illness and photostable materials.

Singlet excited-state lifetimes of cytosine derivatives measured by femtosecond transient absorption

Lifetimes of the lowest excited singlet (S1) electronic states of various derivatives of the pyrimidine nucleobase cytosine (Cyt) were measured by the femtosecond transient absorption technique. The bases were excited in room-temperature aqueous solution at 265 nm using approximately 200 fs pump pulses from a titanium-sapphire laser system. The decay of excited-state absorption (ESA) at visible probe wavelengths was used to determine the S1 lifetimes of a variety of modified Cyt compounds at different pH values by global fitting. Identical lifetimes were observed for Cyt and cytidine (Cyd) within experimental uncertainty, but ESA by the ribonucleoside was considerably stronger, suggesting that the ribose group increases the oscillator strength of the S1 --> SN transition. The S1 lifetime of the important minor base 5-methylcytosine (m5Cyt) is 7.2 +/- 0.4 ps at pH 6.8. The same lifetime was measured for the ribonucleoside 5-methylcytidine, but sugar substitution again increased the strength of the ESA signal. Protonation of Cyd and m5Cyt at low pH led to a modest decrease in their S1 lifetimes. On the other hand, deprotonation of Cyt and m5Cyt significantly increased the lifetime of their respective S1 states. These trends support the intermediacy of the n,pi* state localized on the carbonyl oxygen in the nonradiative decay mechanism of Cyt. Longer S1 lifetimes were observed for 5-fluorocytosine and N4-acetylcytosine. Collectively, these results illustrate the great potential of femtosecond laser spectroscopy for investigating excited-state dynamics in DNA and DNA components.

Ultrafast decay of electronically excited singlet cytosine via a pi,pi* to n(O),pi* state switch

Singlet fluorescence lifetimes of adenosine, cytidine, guanosine, and thymidine, determined by femtosecond pump-probe spectroscopy (Pecourt, J.-M. L.; Peon, J.; Kohler, B. J. Am. Chem. Soc. 2000, 122, 9348. Pecourt, J.-M. L.; Peon, J.; Kohler, B. J. Am. Chem. Soc. 2001, 123, 10370), show that the excited states produced by 263 nm light in these nucleosides decay in the subpicosecond range (290-720 fs). Ultrafast radiationless decay to the ground state greatly reduces the probability of photochemical damage. In this work we present a theoretical study of isolated cytosine, the chromophore of cytidine. The experimental lifetime of 720 fs indicates that there must be an ultrafast decay channel for this species. We have documented the possible decay channels and approximate energetics, using a valence-bond derived analysis to rationalize the structural details of the paths. The mechanism favored by our calculations and the experimental data involves (1) a two-mode decay coordinate composed of initial bond length inversion followed by internal vibrational energy redistribution (IVR) to populate a carbon pyramidalization mode, (2) a state switch between the pi,pi* and nO,pi* (excitation from oxygen lone pair) excited states, and (3) decay to the ground state through a conical intersention. A second decay path through the nN,pi* state (excitation from the nitrogen lone pair), with a higher barrier, involves out-of-plane bending of the amino substituent.

Study of polynucleotide conformation by resolution-enhanced ultraviolet spectroscopy poly(rC) and poly(dC)

Self-deconvolution and the fourth derivative of ultraviolet absorption spectra have been used to study stacked single-stranded and double-helix structures of different cytosine-containing polynucleotides for the first time. These compounds were studied under different solution conditions (pH and organic solvents) and at low temperatures. The red shift of the lower band (B2u band plus possibly some n-->pi* transition) of the absorption spectra in the cytosine-containing polynucleotides and the appearance of new peaks in the deconvoluted and derivative spectra in the 280-310 nm region are attributed mainly to cytosine-cytosine stacking interactions. In particular, the fourth-derivative peaks at wavelengths higher than 290 nm can be associated to coupling of electronic transitions of cytosine bases. The nature of the electronic transitions producing the absorption bands which are resolved in the aforementioned fourth-derivative peaks is discussed. It is concluded that the resolution-enhancement techniques used in this work, i.e. self-deconvolution and fourth derivative, complement each other and are useful methods to study structural changes of single-stranded and double-stranded polynucleotides allowing, at the same time, more information to be obtained about specific stacking interactions than classical absorption spectrophotometry.

Reactivity of propylene oxides towards deoxycytidine and identification of reaction products

Deoxycytidine was reacted with four epoxides of varying alkylating rates: propylene oxide, glycidol, epichlorohydrin and trichloropropylene oxide. Deoxycytidine was chosen to compare the reactivities of these epoxides as all sites of possible alkylation, the oxygen and both nitrogens, are involved in base pairing in DNA. Reaction products were separated on HPLC. Products of the least and most reactive epoxides, propylene oxide and trichloropropylene oxide, were characterized by UV, IR, 360 MHz NMR and MS analysis. For the epichlorohydrin and glycidol reactions and all analytical separations, products were characterized by their HPLC retention times and UV spectra. While differences in reactivity among the epoxides towards specific nitrogen and oxygen sites were found, total reactivity correlated with Taft sigma-values of the substituent groups.

Solvent effects on the kinetics and thermodynamics of stacking in poly(cytidylic acid)

The Raman laser temperature-jump technique has been used to measure the kinetics of the coil to helix reaction of poly(cytidylic acid) [poly(C)] in aqueous cosolvent mixtures. The rate of helix formation has a low activation energy and is proportional to reciprocal solvent viscosity. The observations suggest helix formation is rotationally diffusion controlled. The rate of coil formation in poly(C) has an activation energy of approximately 11 kcal/mol, presumably reflecting the electronic stacking interactions which stabilize the helix. Viscous cosolvents, glycerol or sucrose, slow down the rate of coil formation; acetonitrile and formamide at 5 mol % increase the rate relative to that in water. The polar cosolvents may specifically attack a cytosine stack. The absorbance vs. temperature profiles for poly(C) are analyzed with the one-dimensional Ising model. When only optical data are used, the cooperativity parameter, sigma, and the enthalpy, delta H, cannot be uniquely determined. A method is proposed that allows determination of sigma by combining spectroscopic and calorimetric data. The values of sigma derived for poly(C) are between 0.8 and 1.0, and delta H is about -9 kcal/mol of stack. An alternative method using integration of the excess heat capacity curve and extrapolation to fully stacked and random coil species yields a delta H of -7 kcal/mol of stack.

Stacking properties of a highly hydrophobic dinucleotide sequence, N6, N6-dimethyladenylyl(3' leads to 5')N6, N6-dimethyladenosine, occurring in 16--18-S ribosomal RNA

The thermal denaturation ultraviolet absorption spectra of N6,N6-dimethyladenylyl(3' leads to 5')-N6,N6-dimethyladenosine (m2(6)Apm2(6)A), which is a common sequence in 16--18-S ribosomal RNA, in aqueous buffer at pH 7 have been measured over the temperature range 3-90 degrees C. These data have been used to determine the thermodynamic quantitites associated with the intramolecular stacking equilibria. At 25 degrees C in neutral aqueous solution m2(6)Apmw(6)A exists mainly (about 81%) as a stacked form, so that the stacking interactions are stronger than those in the parent unmethylated adenylyl-(3'-5')adenosine (ApA), where about 52% is stacked. From the parameters of delta H and delta S, it is concluded that 'hidden' hydrophobic inteactions are of prime importance in the enhanced stability of m2(6)Apm2(6)A. Transphosphorylation reaction of ApA and m2(6)Apm2(6)A to form the corresponding cyclic 2',3'-phosphates has been studied. First-order rate constants at 25 degrees C for the reactions, which are base-catalyzed, have been obtained. Insertion of two methyl groups at N-6 of ApA reduces the rate of transphosphorylation. Effects of stacking on rates are discussed in the light of reaction mechanisms.

The nature of stacking interations in polynucleotides. Molecular states in Oligo- and polyribocytidylic acids by relaxation analysis

The dynamics of the helix-coil transition of single-stranded poly(C) (polyribocytidylate) and CpC (cytidyly(3'-5')cytosine) was investigated by an improved cable temperature-jump technique. The single-strand relaxation was characterized by following the ultraviolet (uv) absorbance changes at 248 and 280 nm. Poly(C) and CpC showed single relaxation processes with amplitudes corresponding to those expected from equilibrium melting curves. The relaxation time contants in the range of 25-100 ns were independent of the nucleotide concentration, but strongly dependent upon temperature. Using thermodynanic parameters obtained from circular dichroism (CD) and uv absorbance melting curves, the following rate constants k (at 20 degrees C, 1.05 M ionic strength, pH 7) and activation enthalpies EA were calculated for poly (C): helix formation kR = 1.11 X 10(-7) s-1 (EAR = 2.6 kcal); helix dissociation kD = 2.1 X 10(6) s-1 (EAD = 11.9 kcal). The rate constants obtained for CpC were higher by a factor of about 2 in kR and 12 in kD, whereas the activation enthalpies closely corresponded to those found for the polymer. In addition to the single-stranded helix-coil relaxation, poly(C) and CpC exhibit a relaxation process with a time constant below 25 ns and maximum amplitudes at wavelengths lambda greater than or equal to 285 nm. The same process is found in cytidine and is attributed to hydration equilibria. The hydration reaction can be considered to be in equilibrium during the entire time range of the helix-coil transition and thus the data obtained for the helix-coil transition can be described by a simple two-state model. The rate parameters indicate the existence of relatively high energy barriers in the helix-coil transition and provide strong evidence evidence against an oscillating dimer model. If there is an ensemble of substates for one of the states (as may be expected for the coil form), the energy difference between the populated substates is small compared with the energy difference between the major conformational states.