Next Article in Journal
Structural Characterization of a Rhamnogalacturonan I Domain from Ginseng and Its Inhibitory Effect on Galectin-3
Next Article in Special Issue
Photorelaxation and Photorepair Processes in Nucleic and Amino Acid Derivatives
Previous Article in Journal
Repositioning FDA Drugs as Potential Cruzain Inhibitors from Trypanosoma cruzi: Virtual Screening, In Vitro and In Vivo Studies
Previous Article in Special Issue
Nonradiative Relaxation Mechanisms of UV Excited Phenylalanine Residues: A Comparative Computational Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Perspective

Photophysics and Photochemistry of Canonical Nucleobases’ Thioanalogs: From Quantum Mechanical Studies to Time Resolved Experiments

by
Serra Arslancan
1,
Lara Martínez-Fernández
2,* and
Inés Corral
1,3,*
1
Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, Madrid 28049, Spain
2
Istituto Biostrutture e Bioimmagini-Consiglio Nazionale delle Ricerche, Via Mezzocannone 16, Napoli I-80134, Italy
3
Institute for Advanced Research in Chemical Sciences (IADCHEM), Universidad Autónoma de Madrid, Madrid 28049, Spain
*
Authors to whom correspondence should be addressed.
Molecules 2017, 22(6), 998; https://doi.org/10.3390/molecules22060998
Submission received: 7 May 2017 / Revised: 10 June 2017 / Accepted: 12 June 2017 / Published: 18 June 2017
(This article belongs to the Special Issue Experimental and Computational Photochemistry of Bioorganic Molecules)

Abstract

:
Interest in understanding the photophysics and photochemistry of thiated nucleobases has been awakened because of their possible involvement in primordial RNA or their potential use as photosensitizers in medicinal chemistry. The interpretation of the photodynamics of these systems, conditioned by their intricate potential energy surfaces, requires the powerful interplay between experimental measurements and state of the art molecular simulations. In this review, we provide an overview on the photophysics of natural nucleobases’ thioanalogs, which covers the last 30 years and both experimental and computational contributions. For all the canonical nucleobase’s thioanalogs, we have compiled the main steady state absorption and emission features and their interpretation in terms of theoretical calculations. Then, we revise the main topographical features, including stationary points and interstate crossings, of their potential energy surfaces based on quantum mechanical calculations and we conclude, by combining the outcome of different spectroscopic techniques and molecular dynamics simulations, with the mechanism by which these nucleobase analogs populate their triplet excited states, which are at the origin of their photosensitizing properties.

1. Introduction

Thiobases are the result of substituting an exocyclic carbonyl oxygen atom by a sulfur atom in nucleobases. Despite their structural similarity to their canonical analogs, the photophysics of thiated nucleobases is dramatically different. In contrast to canonical nucleobases, which are able to repopulate the ground state (GS) in ultrafast time scales [1,2,3], the main relaxation pathways in thiated nucleobases drive the system to the most stable triplet excited state in a very efficient manner, that would eventually decay to the GS in longer time scales [4].
Aside from fundamental interest that would seek to understand how the nature of the substituents or their position along the purine and pyrimidine heterocycles [5,6] imprint the topography of the excited state potential energy surfaces (PESs) directing the photophysics of nucleobases, it has been recently recognized that specific thiated nucleobases could have been part of primordial RNA macromolecules [7,8]. Moreover, these systems are currently considered as very valuable prototypes for photosensitizers to be used in pharmacological applications, for their photosensitizing properties and structural similarity with canonical nucleobases, which facilitates their inclusion into parent DNA molecules.
In fact, 2-thiouracil (2t-Ura), 4-thiouridine (4t-Urd), 2-thiocytidine (2t-Cyd) and 2-selenouracil were found in transfer ribonucleic acids (tRNAS) of some bacteria such as Escherichia coli or Bacillus subtiles, yeasts and other higher organisms, and have been related to the improvement of the accuracy and efficiency of the translation process, fostered by the electronic properties and the size of S and Se elements [9,10]. Recent studies on non-enzymatic replication [7,8] have discovered that the exchange of Uracil (U) by 2t-Ura mends important problems along the replication process, such as the slow rate when copying RNA templates or the occurrence of G:T and A:C mispairs [7]. These interesting findings have been interpreted as an indication that thiated compounds might have played a role in primordial RNA replication processes, before the first RNA polymerase enzymes emerged.
As for their pharmacological properties, thiopurines, such as azathioprine, 6-mercaptopurine and 6-thioguanine (6t-Gua), have stood out for their clinical effectiveness as anticancer, anti-inflammatory or immunosuppressant prodrugs for decades [11,12]. Following enzymatic metabolization, a small fraction of these compounds is incorporated as 6t-Gua nucleotides into patients’ DNA undergoing these treatments [13,14,15]. This very small guanine by 6t-Gua substitution (below 0.02%) suffices, however, to increase the sensitivity of 6t-Gua containing DNA to UVA light. The absorbance of UVA radiation by 6t-Gua DNA would precipitate Type I and Type II mechanisms, leading to important photosensitized lesions such as canonical nucleobases’, 6t-Gua’s and protein’s oxidation products, DNA and DNA–protein crosslinks and DNA strand breaks [16].
Nucleobase thioderivatives, such as 4t-Urd, 4-thiothyimine (4t-Thy), and 6t-Gua, have also shown a great utility as photolabels to gain detailed knowledge on the structure of nucleic acids, as well as to recognize interaction points within nucleic acids or in nucleic acids protein aggregates, profiting their photocrosslinking properties upon UVA light absorption [17,18,19,20,21,22,23].
Finally, thiation has also been intensively employed in structural biology studies to understand the chemical, structural and functional alterations introduced by these substituents in the nucleic acids where they were incorporated [24,25,26,27,28].
In this scenario, it is not surprising that the literature on the photophysics and photochemistry of thionucleobases has experienced major growth, especially in the last decade, from both the experimental and computational sides, with the aim of unraveling the differences at the electronic structure level that distinguish these compounds from canonical nucleobases and deciphering the relaxation pathways underlying the very appealing photosensitizing properties of these systems for pharmacology and medicinal chemistry. Considering the abundant literature focused on particular aspects of the photophysics and photochemistry of thiated nucleobases, it is our aim here to compose a concise but complete picture, sewing together contributions from different fields, which summarizes not only the main optical properties of these systems, but also their main deactivation pathways.
As for canonical nucleobases, tautomer equilibria in thiobases have been intensively studied both from computational and experimental points of view. Thiouracils and their N1-methylated derivatives were found to be mainly present in their keto-thione forms, both in gas phase and in solution [29,30,31,32]. On the contrary, theoretical investigations on 2-thiocytosine (2t-Cyt) concluded that in vacuum its thiol form is more stable than the canonical tautomer, although the latter becomes strongly predominant in polar solvents [33,34]. Similar trends have been described also for the tautomeric equilibrium of 6t-Gua [35,36]. Since keto-thione is the predominant form for all thiobases in a biologically relevant environment such as water, most of the theoretical works have been performed on these tautomers and thus it will be on them where we will put the focus of our review.
We will also overlook recent works on the photophysical properties of Se- and Te-nucleobase derivatives [37,38], in line with those of thiated nucleobases, and the characterization of the electronic excited states of ionized thionucleobases, key to the interpretation of UV photoelectron spectra recently measured/simulated for some thiouracil derivatives [39,40].
The review is organized as follows. For all the thionucleobases considered, namely the thiopurine 6t-Gua, the thiopyrimidines 2-thiothymine (2t-Thy) and 2t-Ura, 4t-Thy and 4-thiouracil (4t-Ura), 2t-Cyt and 2,4-dithiothymine (2,4-dtThy), and 2,4-dithiouracil (2,4-dtUra), we will start by compiling the description of the experimental absorption and luminescence spectra and their interpretation with the help of Time Dependent Density Functional Theory (TD-DFT) and high level ab initio calculations. When discussed in the literature, we will briefly comment on the effect of the solvent and glycosylation on the main absorption and emission bands, by comparing with the absorption and emission spectra of the corresponding (deoxy)ribosides. Next, we will summarize the main features of the PES explored with the above-mentioned theoretical methods, sketching the most plausible relaxation mechanisms. Then, the photodynamics of these systems will be discussed based mainly on the results of transient absorption experiments, and in specific cases of time-resolved photoelectron and Resonance Raman spectroscopies. We will also comment on the fate of the triplet excited states, responsible for the photosensitizing properties of these systems. Detailed overviews on the fundamentals and performance of the main theoretical methods and spectroscopies discussed along this review can be found in Refs. [41,42,43,44,45,46,47,48,49].
Finally, the last section of the review concludes with some general remarks on the effect that the degree and position of thiosubstitution has on the absorption spectra, PES and photoexcited dynamics of thionucleobases and a very brief description of the main challenges in the field.

2. Canonical Nucleobases’ Thiated Analogs Photophysics

2.1. Thiopurines

2.1.1. Steady State Absorption and Emission properties

Table 1 collects the experimental absorption and fluorescence emission maxima of 6t-Gua and its nucleosides 6-thio(2′-deoxy)guanosine (6t-Guo) in different solvents. Table 2 also summarizes in vacuum and in solution theoretical vertical absorptions and adiabatic emission energies calculated for 6t-Gua and 6t-Guo. According to these experiments, both the nucleobase and the nucleoside exhibit absorption maxima in the UVA region between 348–316 nm, (3.56–3.92 eV) for all the solvents examined, alkaline media providing the lowest wavelength values, as already noticed by Rubin et al. and Santhosh et al. [50,51]. Other peaks were found around 255, 220 and 205 nm in neutral and acidic conditions, whereas alkaline pHs induce a bathochromic shift in these absorptions which appear at 265, 245 and 210 nm [51]. DFT [38,52,53,54] and multiconfigurational approaches [55] predict the spectroscopic state in 6t-Gua, to absorb between 311–299 nm (3.99–4.15 eV) in gas phase (GP), whilst the computationally more restrictive semi-empirical CNDO (Complete Neglect of Differential Overlap) method [50,51] calculates this transition at significantly lower energies ca. 340 nm. All these methods, however, agree in assigning this transition a ππ* (Highest Occupied Molecular Orbital (HOMO) to Lowest Unocccupied Molecular Orbital (LUMO)) character and locate a very weakly absorbing lower-lying nπ* electronic state in the region 400–369 nm (3.10–3.36 eV). These two transitions are localized within the C-S bond (see Figure 1). TD-DFT predicts a shift of these two absorptions after including solvent-chromophore interactions, i.e., whilst the nπ* absorption is displaced to higher energies (ca. 0.4 eV) both in water and acetonitrile (ACN), the ππ* transition blue shifts by 0.1 eV in the same solvents.
Discrepancies between the different theoretical approaches employed to predict the absorption spectrum of 6t-Gua are larger when predicting higher-lying excited states. Semi-empirical and TD-DFT methods locate a second ππ* excited state centered at ca. 295 nm (4.20 eV), absorbing with different intensities (0.0076–0.48), whilst the more reliable MS-CASPT2/CASSCF (Multi State Complete Active Space Perturbation Theory/Complete Active Space Self-Consistent Field) predicts a quite bright state at higher energies 253 nm (4.90 eV).
The comparison between the absorption features of 6t-Gua and its nucleosides (Table 1) reveals that glysosylation has very little effect in the position of the absorption maxima and only slightly increases the intensity of the UVA and UVC bands.
Fluorescence spectra of 6t-Gua recorded at different pH conditions and excitation wavelengths showed a band with a maximum peaking at 400–468 nm, with a shoulder around 500 nm [51,52]. Fluorescence decay was best fitted with two lifetimes in the range of ns, i.e., τ1 = 1.4 and τ2 = 24.1 ns [52].
Finally, phosphorescence signals centered in the range of 430–550 nm, depending on the pH conditions and on the experimental setup, were also registered [50,52]. Phosphorescence emission was found to decay with a lifetime of 45 ms [52].

2.1.2. Static Description of the PES

A scheme of the potential energy profiles for 6t-Gua based on minimum energy path (MEP) calculations can be found in Figure 2.
The photophysics of 6t-Gua (6t-Guo) is dictated by the existence of two minima in the S1 potential, showing ππCS* and nπCS* characters, and lying 3.80 (3.5) and 3.14 (3.3) eV over the GS equilibrium structure, according to the MS-CASPT2//CASSCF [55] (TD-B3LYP [38]) calculations. These two minima are accessed via the S2(ππCS*)/S1(nπCS*) conical intersection (CoIn) in the vicinity of the Franck Condon (FC) region. Reaching these minima requires in 6t-Gua the lengthening of the C-S bond, besides the out-of-plane movement of the C-S bond in the nπCS* minimum or the pyramidalization of the C2 (see Scheme 1 for atom labeling) coupled to the rotation of the amino group in the ππCS* minimum [55]. The S1 and S2 (corresponding to S1(ππCS*) in Ref. [55]) minima of the deoxynucleoside were found to experience parallel geometrical C-S distortions as those reported for the thionucleobase. No significant out-of-plane deviation of the NH2 group was, however, registered for 6t-Guo [38]. Similar geometries were obtained for equivalent minima of triplet multiplicity located 3.41 (3ππCS*) eV and 3.07 (3CS*) eV relative to the GS minimum in 6t-Gua [55]. A lower energy (2.5 eV), however, was registered for the most stable T1 minimum of 6t-Guo [38]. Interestingly, a third triplet minimum with ππCS* character, identical geometry and thus energy to the nπCS* singlet and triplet minima was also optimized for the thionucleobase [55].
Although a one to one comparison between the calculated adiabatic emission energies and the fluorescence and phosphorescence (Table 1 and Table 2) signals is not possible, since the vibrational structure of the ground and S1 excited state was not considered by the available calculations, the energies calculated for the most stable singlet and triplet minima are consistent with the luminescence signals experimentally obtained.
The same multiconfigurational study [55] locates for 6t-Gua two S1/S0 internal conversion (IC) funnels for the return of population to the S0 from the ππCS* and nπCS* minima. Reaching these funnels requires the perpendicular arrangement of the C-S bond to the purine ring, and the pyramidalization of several C and N centers in the case of the S1(nπCS*)/S0, importantly affecting the energy of this interstate crossing point which lies almost 2 eV above the nπ* minimum. The S1(ππCS*)/S0 crossing, with almost the same energy as the ππCS* minimum, is in contrast easily accessible for the S1 potential. A third S1(ππCC*)/S0 IC funnel, lying at intermediate energies (3.9 eV), was identified with the help of molecular dynamics (MD) simulations on photoexcited 6t-Gua [58]. This interstate crossing is geometrically characterized by the puckering of the C7 and the folding of the NH2 group, and is accessible after surmounting a 0.72 eV energy barrier from the 1CS* minimum.
The two S1 minima, lying both structurally and energetically very close to the minimum energy crossing points (MECP) S1(ππCS*)/T2(nπCS*) and S1(nπCS*)/T(ππ*), were identified as efficient funnels for transferring the population to the triplet manifold, due to the proximity of strongly coupled triplet states (calculated spin-orbit couplings (SOC) of 200 cm−1). Once in the triplet manifold, the system needs to climb energy barriers comprised between 0.1 and 1 eV to access intersystem crossing (ISC) funnels for further decaying to the S0. The SOCs calculated at these regions are small (3–12 cm−1), except for the T1(nπCS*)/S0 crossing, which was predicted to amount to 270 cm−1. It is important to note here that all the stationary and interstate crossing points were connected with MEP calculations.

2.1.3. Photodynamics

The excited-state dynamics of the nucleoside 6t-Guo and the nucleobase 6t-Gua were respectively studied by means of femtosecond broadband transient absorption spectroscopy (TAS) [52,54] and singlet/triplet surface hopping dynamics using a parameterized semi-empirical Hamiltonian with CASPT2 calculations [58]. The excited state dynamics of the nucleoside (nucleobase) is governed by a single excited state species showing two transient absorption signals at 380 (375) and 620 (520) nm, which evolve in parallel at early time delays and remain stationary after 2 ps. This transient is assigned to a triplet excited state, based on the comparison of the decay lifetimes of the absorption signals in N2 and air-saturated conditions, which are 1.5 times increased in the absence of oxygen [54]. The triplet species in the nucleoside is populated in ultrafast timescales τ1 = 0.31–0.36 ps, both in aqueous solution and in ACN (See Figure 3, and Table 3). For the nucleobase, however, this process was measured to occur in longer time scales τ1 = 0.56 ps [52].
A second lifetime for the nucleoside (nucleobase) τ2 = 32–80 (26) ps was also fitted from the dynamics and assigned to the minor repopulation of the GS from the S1(nπCS*) excited state combined with solvation dynamics of the T1 (See Figure 3, Table 3) [54]. Finally, the longer lifetime in 6t-Guo (6t-Gua) τ3 = 720 (830–1420) ns is ascribed to the ISC process from the triplet manifold to the GS [52,54].
From the above data, it can be concluded that the effect of N9-glycosylation is not homogenous along the dynamics of the chromophore. While N9-glycosylation accelarates the ISC processes for the population of the triplet manifold and the GS, it noticeably slows down the IC process S1→S0.
The complex topography of the excited state PES of thiobases invalidates the interpretation of the photodynamics of these systems based on the FC picture and, moreover, urges the use of MD simulations due to the existence of several competing relaxation pathways [54,58]. The results from the singlet/triplet MD simulations [58] are in line with the experimental results [52,54] and consistent with the topography of the PES predicted by static quantum chemical methods [55]. In the dynamics simulations, the S2 internally converts in ultrafast timescales to the S1. The depopulation of the S1 in favor of the T1 and T2 states is also ultrafast. Only very few direct hops were observed from the S2 to the triplet manifold. The very small population (0.15%) that returns back to the GS, during the 10 ps that the simulations last, proceeds via the S1(ππCC*)/S0 CoIn which is accessed from the S1CS* minimum after overcoming a small energy barrier through which the S1 electronic state evolves and adopts a ππCC* character (Recall Figure 2). Interestingly, the dynamics is able to capture critical features of the topography of the PES, such as the degeneracy between the S1, T1 and T2 electronic states around the region of the most stable S1 minimum, where the system gets trapped during some time. The energetic proximity of these three electronic states translates into multiple recrossings among them along the trajectories examined. At the final time of the simulations, the population is fundamentally confined in the T1 electronic state, with some excited 6t-Gua molecules in the T2 and residual population distributed between the S1 and S0 [58]. The fitted lifetimes extracted from the MD simulations and which amount to τ1 = 122 fs and τ2 ≈ 200 ps are of the same order of magnitude as the experimental results obtained for the 6t-Gua and 6t-Guo.
Finally, a very high triplet yield of 0.9 [58] was inferred from the MD simulations, which is slightly higher than the experimental value of 0.6–0.8 [52,54].

2.1.4. The Fate of the Triplet Excited States of 6t-Gua and 6t-Guo

The very extended use of thiopurine prodrugs, which are enzymatically converted to 6t-Gua nucleotides before being incorporated into DNA, as immunosuppressants, anti-inflammatory and anticancer drugs greatly contributed to intensifying the research on the fate of the triplet excited state of this nucleobase and its nucleoside to better understand their mechanism of action. The chemical evolution of these triplet intermediates importantly depends on the environment where the nucleobase derivative is immersed and where it is expected to undergo Type I and Type II photosensitization reactions, leading either to cellular lesions, if 6t-Gua happens to be incorporated in a DNA macromolecule or to potential toxic products if photosensitization takes place in solution.
In fact, Type I mechanisms rely on the reaction of the photosensitizer with a substrate to generate substrate’s and photosensitizer’s radical species, via hydrogen abstraction or electron transfer processes, that would eventually evolve to oxygenated products or would produce the superoxide radical O2 following an electron transfer process in the presence of molecular oxygen.
Type II photosensitization, in turn, involves an energy transfer process, where triplet excited 6t-Gua transfers its excess of energy to molecular oxygen from the environment, which is excited to singlet oxygen, with the concomitant return of the thiobase back to the GS.
Zhang et al. determined in 2011 the quantum yield of 1O2 production (ΦΔ) for 6t-Gua and 6t-Guo, by following the oxidation of the water soluble phosphine, [2-(dicyclohexylphosphino)ethyl]trimethylammonium chloride by 31P NMR [56]. These authors obtained, depending on the pH conditions, ΦΔ values which oscillated within 0.56–0.58 for 6t-Gua and 0.49–0.55 for 6t-Guo (see Table 3). Only recently, Pollum et al. estimated the 1O2 quantum yield for 6t-Guo using 1O2 phosphorescence and obtained significantly lower values compared to previous studies, ΦΔ = 0.24–0.29 [73]. These unexpectedly low ΦΔ, considering the high triplet quantum yields measured for these systems, were ascribed to the low values determined for the fraction of quenching events (SΔ), which for 6t-Guo only amounts to 0.37, suggesting the concurrence of other competing events [73].
The 1O2 produced from Type II photosensitization reactions has been postulated to be involved in the oxidation reactions of 6t-Gua, canonical nucleobases and proteins that can be considered as major routes for DNA and/or cellular damage [16]. The main oxidation products of 6t-Gua are guanine-6-sulfinate (GSO2) and guanine-6-sulfonate (GSO3) [15,56,74]. Interestingly, guanine is the canonical nucleobase more prone to oxidation, producing 8-oxoguanine (8-oxo-7,8-dihydroguanine, See Scheme 2) [75]. The production of 1O2 has also been found to induce histidine-lysine crosslinking processes in subunits of the PCNA replication and repair complex [76].
In 2013, Zou and coworkers investigated, by combining DFT calculations and experimental measurements of reaction stoichiometry and rate constants, the mechanisms by which the mutagenic GSO2 and GSO3 products are generated [77]. Specifically, these authors invalidate the mechanism, which postulates guanine sulfenate (GSO) as an intermediate in the formation of GSO2 and GSO3 products [15,74,78], and propose instead the more energetically favored formation of GSO2 through a peroxy intermediate (GSOOH), accessed after H transfer from the NH group of the pyrimidine ring towards the terminal O atom of the 6t-Gua-1O2 complex. This peroxy intermediate would then evolve in a second stage to GSO2 product, after the migration of the OH group from the O to the S center (see Scheme 2). Further oxidation of GSO2 by 1O2 would lead to GSO4 and to the final oxidation product GSO3 with the assistance of a water molecule. The very low barriers theoretically predicted for the oxidation 6t-Gua→GSO3 are supported by the fast rate constant measured, which amounts to 4.9 × 109 M−1 s−1. These authors additionally mapped at molecular level a secondary pathway competing with the formation of GSO2 which explains the residual formation of Guanine + sulfur monoxide, SO, in low yields, due to the larger energetic barrier involved.
The complexity of other UVA-photoinduced 6t-Gua photoproducts, such as intra- and interstrand [79] DNA crosslinking, DNA breakage [80] or DNA–protein crosslinking [17,81,82], has prevented the in detail study of the mechanism behind the formation of these photolesions, remaining in some cases the role played by Type I and Type II photosensitization reactions still unknown. In fact, some recent works on 6t-Gua and 6t-Guo have demonstrated, based on theoretical grounds, that electron transfer (Type I mechanisms) reactions between the triplet excited states of 6t-Gua or 6t-Guo and ground state molecular oxygen are energetically favorable, and might therefore offer a plausible explanation for the reduced singlet oxygen quantum yield of this chromophore [73,83]. In fact, Pollum et al. estimate the value of ΔG for the electron transfer reaction in 6t-Guo to amount to −0.72 eV, using the Rehm-Weller equation and the vertical ionization energy of 6t-Guo predicted by DFT calculations in solution [73]. A moderately lower value (−0.26 eV) was calculated in Ref. [83] considering the complex formed between the triplet excited state of 6t-Gua and 3O2 and using the CASPT2//CASSCF protocol, including solvent solute interactions via the Polarizable Continuum Model (PCM) model.

2.2. Thiopyrimidines

2.2.1. 2-Thiopyrimidines

2t-Cyt/Cyd

• Steady State Absorption and Emission properties
The absorption spectra of 2t-Cyt and 2t-Cyd in water were recorded in the early 1980s, showing both a maximum at ca. 270 nm (~4.59 eV, see Table 4) [84]. More recently, the absorption spectrum of 2t-Cyt has been revisited confirming the presence of a maximum at 270 nm (Band I) and ascertaining the existence of a second peak at ~240 nm (Band II, 5.17 eV) and a shoulder ~220 nm (5.64 eV) [34,85]. The two bands are affected in different ways by the change of the polarity and the nature of the solvent from polar-protic (water) to polar-aprotic and less polar-aprotic solvents (DMSO) [85]. Whereas Band I is red-shifted by ~15 nm, but maintains its intensity, when moving from water to DMSO, the position of Band II is less affected but significantly decreases its intensity, vanishing in DMSO [85].
GP MS-CASPT2 calculations (Table 5) for 2t-Cyt predict the first bright state S2Sπ*) at ~330 nm (3.76 eV), which has been connected with the low-energy tail of Band I at 270 nm. The absorption maximum of Band I has been ascribed to the second bright state S4Sπ*) (~4.30–4.40 eV) [34,85]. The oscillator strengths computed for both states confirm this assignment, since the S2Sπ*) is at least five times less intense than S4Sπ*) [34,85]. The S1 (nSπ*) corresponds instead to a dark excited state [34,85]. Higher in energy bright states, S6 and S8, have been recently proposed by Mai et al. to contribute to Band II [85]. These authors have paid especial attention to the performance of MS-CASPT2 in the prediction of solvated absorption spectra and their comparison with the experiment. The computed dipole moments of the electronic excited states together with results from explicit solvent-solute calculations were used to explain the shift of the transitions induced by the solvents: S4Sπ*) with a smaller dipole moment than the GS was found to blue-shift in polar solvents, whereas the opposite was found for S6 and S8, in agreement then with the experimental observations [85]. This blue-shift of the first absorption band (arising from a 1ππ* excited state) in water was previously reported for other purine/pyrimidine bases [5,87,88]. However, the accurate reproduction of the experimental changes in the band intensities was found to be more challenging [85].
• Static Description of the PES
After excitation at 310–320 nm wavelengths, which are the ones used in femtosecond transient absorption experiments, the S2Sπ*) is mainly populated. This excited state is expected to rapidly decay in favor of two lower-lying S1 minima with 1πSπ* and 1nSπ* characters. Indeed, GP MS-CASPT2 MEP calculations reveal that the IC funnel allowing the S2→S1 decay is reached barrierlessly from the FC region, see Figure 2 [34]. Both singlet minima maintain the pyrimidine ring planarity but show some differences for particular bond lengths, the most significant concerns the C5–C6 bond that is larger for the 1πSπ* minimum. Besides being geometrically similar, MS-CASPT2 also predicts excited state features common to both minima. First, they are stable against GS repopulation, since the IC funnels with the S0 were located significantly high in energy (~0.8 eV) over the corresponding minima [34]. By contrast, the presence of close lying singlet/triplet MECP (1πSπ*/3nSπ* or 1nSπ*/3πSπ*), where large SOC were calculated (160 cm−1), identifies these singlet minima as ideal regions of the PES for singlet→triplet population transfer [34]. Thereby, 2t-Cyt does not present any long-lived singlet minima and, although not having yet been reported, the experimental fluorescence quantum yields are expected to be negligible.
Once in the triplet manifold, a CoIn between the two triplet states, 3πSπ*/3nSπ*, was located close in energy to both singlet-triplet MECPs (1πSπ*/3nSπ* and 1nSπ*/3πSπ*). Thus, according to MEP calculations both triplet minima (3nSπ* or 3πSπ*) could be a priori populated (see Figure 2). Both triplet minima structurally resemble their singlet counterparts. Non-radiative GS repopulation from these triplet minima seems also unlikely due to the presence of large energy barriers to access the ISC funnels [34]. Unfortunately, no experimental phosphorescence measurements for 2t-Cyt are so far available.
Two different 3πSπ* minima have been recently characterized by Bai et al. at CASPT2 level for 2t-Cyt (Figure 2). Similarly to 6t-Gua, the two minima differ in their geometries and in their decay mechanisms. Although the relative stability of the two 2t-Cyt triplet minima has not been discussed, the first shows a clear out of plane deviation of the S atom and is expected to easily decay to the GS, whereas the other presents a ring-distorted structure and is expected to live longer [59].
• Photodynamics
2t-Cyt has also been investigated by means of femtosecond TAS using a 308 nm excitation wavelength [34]. For 120 fs time delay onwards, the TAS shows two different absorption maxima, at 350–400 nm and at ~600 nm. For a time delay of 320 fs, the former reaches its maximum whereas the latter blue-shifts (to ~525 nm) and continues to grow also in the ps time regime, being its decay not observed within 20 ps. Two different time constants in the fs regime were found to correctly fit these TAS, τ1 = 210 fs and τ2 = 480 fs and a close-to-unity triplet yield was reported (see Table 3, Figure 3) [34].
GP Multi Reference Configuration Interaction Singles (MRCIS) dynamics simulations showed a fast decay of the 2t-Cyt S2 population, in 160 fs, in agreement with the absence of stable minima along the S2 PES. An average ISC time constant of 250 fs was calculated, which is in good agreement with the experimental value of τ1 = 210 fs (See Figure 3) [34]. Although it is not possible to discern between the two different ISC routes (1πSπ*/3nSπ* or 1nSπ*/3πSπ*), MRCIS dynamics predicts that depopulation from the S1 to the T2 is preferred. Thus, the first time constant would account for the first deactivation step involving S2→S1→T2 relaxation [34]. Once in the triplet manifold, subsequent T2→T1 IC takes place, holding T1 most of the population from 400 fs onwards. Thus, the second experimental time constant τ2 = 480 fs was ascribed to IC in the triplet manifold [34].
Despite their structural similarities, the 1πSπ* and 1nSπ* minima present very different absorption patterns. According to MS-CASPT2 calculations, they respectively absorb at 360 and 600 nm, with different intensities, the 1πSπ* minimum presenting much larger absorption [34]. Thereby, the comparison of experimental TAS at representative time delays with computed vertical absorptions at relevant points of the excited PES has been used to provide further mechanistic information on 2t-Cyt photoinitiated dynamics. As described above, for a 320 fs time delay, the experimental TAS reaches a maximum at ~350 nm and presents some absorption also at ~600 nm, which Mai et al. attribute to a ~10% and ~25% population distribution between the 1πSπ* and 1nSπ* minima, respectively [34]. This residual population at the singlet minima at 320 fs is in agreement with the ISC lifetime (200 fs); in other words for these timings most of the molecules have reached the triplet manifold [34]. For longer time delays, ~4 ps, the recorded experimental TAS is blue-shifted in agreement with the MS-CASPT2 vertical absorption computed for the 3πSπ* minimum centered in the 500–550 nm region [34].
2t-Cyt presents stable long-lived triplet minima, which are both efficiently and rapidly populated from short-lived singlet excited states.

2t-Thy/Thd and 2t-Ura/Urd

• Steady State Absorption and Emission properties
In aqueous solution, the experimental absorption spectrum of the nucleobase 2t-Ura covers the range 360–240 nm (3.44–5.17 eV) with a maximum of absorption around 265 [61]–274 [84] nm (4.68–4.52 eV, Table 4). Interestingly, the absorption band slightly modifies its shape when recorded in ethanol, methanol and ACN. In ethanol or methanol, it presents an additional shoulder at 291 nm (4.26 eV) that increases its intensity when switching the solvent to ACN [84,86]. According to experiments, glycosilation was found to have no effect on the transition coinciding with the absorption maximum S2 (ππ*), however it shifts by 0.05–0.1 eV to higher energies the S1 state and redshifts by 0.2 eV the S3 [84].
Methylation at position 5 has also very little effect on the absorption properties. In fact, 2t-Thy shows maximal absorption around 275 nm in phosphate buffer solution (PBS), whilst two absorption maxima at 275 and 290 nm, with different relative intensities depending on the solvent used, were registered when recording the spectrum in ACN or ethanol [62,63,86]. As for 2t-Ura, Glycosilation only shifts by 2–10 nm the absorption maximum to lower energies, whilst significantly increasing the value of the molar absorptivity coefficient [63,64]. Different quantum chemical methods have been employed to predict the absorption spectrum of these systems, which cover from monoconfigurational DFT [31,37,62,86,90,92] to multiconfigurational [89,90,91] or multireference [89,93] approaches, such as CASPT2/CASSCF and MRCI (See Table 5). For 2t-Ura, these methods predict in GP the lowest lying brightest excited state S2Sπ456*), (the subscripts refer to the pyrimidine centers where the π* orbital is localized), which promotes an electron from the HOMO mainly localized on the S atom to the LUMO extended over the C4-C5-C6 atoms, in the range between 304–215 nm (4.08–5.77 eV) (recall Scheme 1 and Figure 1). This transition is ascribed to the shoulder at 4.26 eV of the experimental spectrum, which becomes apparent in solution [84]. From the wide variety of methods employed in the simulation of the absorption spectrum of 2t-Ura, MRCIS is the one that provides the largest deviation compared to the others. The poor performance of this method has been ascribed to the importance of higher excitations or the use of a larger reference space in the correlation treatment for the correct description of the electronic transitions in this system [89].
Solvent solute interactions included using the PCM model shift this transition by 0.25 eV to higher energies, both for water [31] and ACN [31,92]. All these electronic structure methods also calculate two very weak nπ* transitions flanking the spectroscopic state. The first transition, localized at the C-S bond, was found to be centered ca. 347–322 nm (3.57–3.85 eV). The second, however localized on the C-O bond, was found to peak ca. 300–260 nm (4.13–4.77 eV). These two transitions are also displaced to higher energies when chromophore-water or chromophore-ACN interactions are considered [31,92]. Finally, a second quite intense ππ* transition (πSπ2*), transferring electron density from the C-S bond to the C-O bond and importantly contributing to the absorption band, was calculated in the region 278–243 nm (4.46–5.10 eV). The effect of introducing a methyl group in position 5 (2t-Thy) leads to a blue shift of the calculated electronic transitions, by 0.1–0.4 eV in the case of the first nSπ2* transition and ca. 0.1 eV or less in the case of the lowest lying πSπ456* [62,86,93]. In a water continuum, the S1 and S2 transitions in 2t-Thy were predicted at 3.81 and 4.30 eV by the TD-M06 formalism [37]. MS-CASPT2, however, predicts the second πSπ2* to absorb at the same position for 2t-Ura and 2t-Thy, although shifts to higher energies the second nOπ6* in 2t-Thy, inverting the ordering of these two consecutive states [93].
Kuramochi et al. [62], Taras-Góslińska et al. [63] and Vendrell-Criado et al. [86] studied emission in 2t-Thy nucleobase and its riboside and deoxyriboside. Very small or no fluorescence quantum yields (<10−4) were detected for 2t-Thy [62,86]. Phosphorescence of the nucleobase and the nucleoside in 2-methyltetrahydrofuran (THF) and dichloromethane–methanol (DCM:MeOH) glassy matrices at 77 K and ethanol (EtOH) or deoxygenated ACN solution at room temperature was observed, with maxima at 454 [62], 448 [63], 451 [86], and 480 [63] nm, respectively. A similar value, at 427 nm [86], was obtained for 2t-Ura in ethanol. Phosphoresce quantum yields were also estimated to be very low (5 × 10−4) [63].
• Static description of the PES
2t-Ura and 2t-Thy present a complex excited state PES whose main features are summarized in Figure 4. As a first approach to the problem, Gobbo et al. [91] and Cui et al. [90] locate, by means of CASSCF methods, two minima in the S2 and S1 potential preserving the ππ* and nπ* characters of the FC region (See Figure 4). At this point, it is important to highlight that all the works on 2t-Ura seem to provide a different FC state ordering [89,90,91]. According to Refs. [90] and [91], both minima are characterized by a significant elongation of the C2-S bond. However, whilst Ref. [91] predicts moderate out-of-plane distorted structures for both minima, caused by the pyramidalization of the C6 center, the Cs symmetry constraints imposed in the work of Cui prevent accessing the real minima [90]. A second additional minimum in the S2 potential, with πSπ2* character, but corresponding to the FC S4 electronic state, is latter identified by the group of González [89], using CASPT2 gradients. This minimum maintains the stretched C-S bond distance but shows a strong pyramidalization at the C2 position. Interestingly, González [89] (2t-Ura) and Barbatti (2t-Thy) [93] predict a very similar strongly C2 pyramidalized structure for the S1 nSπ2* minimum, in contrast with the quasi-planar structure calculated by Gobbo and Cui [90,91]. In the triplet manifold, altogether three triplet minima have been identified. Gobbo et al. and Cui et al. locate two minima with nSπ2* and πSπ456* characters in the T2 and T1 potentials (See Figure 4). Similar to what happened for the singlet manifold, CASPT2 predicts for the triplet nSπ2* minimum a very important out-of-plane deviation of the C-S bond relative to the pyrimidine ring. Two minima with πSπ2* and πSπ456* characters along the T1 potential, the first pyramidalized at the C2 position and the second showing a boat like conformation were also found [89,93]. The adiabatic relative energies of the most stable triplet states of 2t-Ura, which were calculated in the range 2.85–3.46 eV, compare very reasonably with the experimental 0-0 emission phosphorescence band from Ref. [94] at 3.17 eV. This is also the case of 2t-Thy where the computed energies for the two ππ* T1 minima amount to 3.20 and 3.23 eV [59], which compare reasonably well with the maxima of the experimental phosphorescence spectra at 2.73–2.76 eV [62,63].
The former stationary points were connected by means of MEP calculations. Two S2/S1 crossings were optimized for each of the two ππ* S2 states to locate the IC funnels for the transfer of population from the spectroscopic state S2 to the lower lying S1. The S2Sπ456*)/S1 CoIn was found to lie close to the position of the S2Sπ456*) minimum, therefore being energetically accessible from this stationary point [89,90,91]. The S2Sπ2*)/S1 CoIn instead lies 0.35 eV above the S2Sπ2*) minimum becoming less accessible compared to the former. These theoretical works also considered the IC funnels for the return of population to the GS. Authors in refs. [90,91] compute the S2Sπ456*)/S0 CoIn respectively ca. 0.05 and 2 eV over the S2 minimum. Moreover, Gobbo et al. [91] predict an energetic barrier to access the S2/S0 CoIn, adding further difficulty to accessing this funnel. Mai et al. fail to optimize the S1(nSπ2*)/S0 degeneracy point, and assume it should not play an important role in the photodynamics of 2t-Ura on the basis of the very large experimental ISC yields. Transfer of population to the triplet manifold has been postulated to take place both from the S2 and from the S1 potentials. Gobbo et al. [91] and Cui et al. [90] locate close lying triplet excited states, showing nπ* characters, at the position of the S2Sπ456*) minimum, where they calculate moderate to large SOCs (20–190 cm−1), in line with El Sayed propensity rules [97]. More probable is, however, the leak of population to triplet potentials at the position (in the vicinity) of the S1 minimum where the calculated SOCs amount to 100–275 cm−1 [89,90,91,93]. Finally, the groups of Fang [90] and Barbatti [59] also explored the paths for the return of the population from the T1 potential to the S0. Cui and Fang (2t-Ura) compute a very high T1/S0 degeneracy points ca. 1.9 eV above the triplet minimum [90]. Bai and Barbatti (2t-Thy) instead calculate this ISC point only 0.3 eV above the C2-pyramidalized T1 minimum.
The complex topography of the excited state PES with multiple possible relaxation pathways called for a dynamics study on this system in order to learn more about the real deactivation mechanism.
• Photodynamics
The dynamics of 2t-Ura and 2t-Thy have been studied using nanosecond to femtosecond transient absorption experiments [60,61,62,63,64], Resonance Raman spectroscopy [92], time-resolved photoelectron spectroscopy [39] and ab initio singlet/triplet non-adiabatic MD simulations [98]. As for 6t-Gua, the dynamics of the 2-thiosubstituted derivatives of Thymine and Uracil is also dominated by a single excited state species absorbing in the visible and UVA regions, which rises in ultrafast time scales, and which shows absorption maxima around 350, 420 and 700 nm. [60,62,63,86] These values, that can slightly vary depending on the specific nucleobase and the solvent considered, are in qualitative agreement with the vertical transient spectra calculated by Bai and Barbatti at DFT/MRCI level of theory, confirming the population of the C2-pyramidalized and boat ππ* triplet minima along the dynamics [59]. The TAS of the ribosides were found to be very similar to that of the nucleobase [63]. This transient absorption signal is preceded by another UV signal (345 nm) populated within few hundreds of femtoseconds and which has been attributed to the absorption of the S1(nπ*) transient [60]. The photodynamics of these species can be then summarized in terms of the following lifetimes τ1(2t-Thy, 2t-Ura) ≤ 200 fs, τ2(2t-Thy) = 0.62 (aqueous solution), 0.32 (ACN) ps and τ2(2t-Ura) = 0.35 (aqueous solution), 0.34 (ACN) ps [60,64], and τ3(2t-Thy) = 2.7 μs, which were respectively assigned to the population of the S1(nπ*), the ISC lifetime and the triplet lifetime of these two nucleobases (See Table 3 and Figure 3) [62,63]. Glycosilation of 2t-Thy was found to reduce by a factor of 1.5 the ISC lifetime [64] and by a factor of 14–22 the lifetimes of the T1 [63], which the authors suggest can be ascribed to the higher density of vibrational states in the nucleosides.
Notably, different time constants for 2t-Ura as compared to previous works in solution were determined by GP time resolved photoionization measurements [39]. According to TR-PES experiments in GP, the S1 would be populated faster (τ1 < 100 fs) than in solution. The time constant extracted from the TR-PES spectra for the S1 ISC process is of the same order as the one calculated for 2t-Thy (τ2 = 775 fs). Quite surprisingly, ISC process for the repopulation of the GS seems to take place remarkably faster (τ3 = 203 ps) compared to other works considering similar systems in solution [60,62,63,86].
MD, using CASPT2 as the electronic structure method for the calculation of the energies, the gradients and the couplings, were undertaken to shed some mechanistic light on the ISC process [98].
In agreement with the experimental results, the dynamics simulations predict S2→S1 IC process and ISC to the triplet manifold in a few hundreds of femtosecond. From all the possible pathways to transfer population to the triplet manifold, the dynamics revealed the S1→T2 route as the main one (40%), followed by the S1→T1 (35%) and S2→T2 (25%). The lifetimes fitted from the dynamics simulations, τ1 = 59 fs, τ2 = 400 fs, and the triplet quantum yield of ΦT = 0.90, are also in line with the experimental findings.
The MD simulations also revealed important features of the relaxation mechanism. For instance, most of the trajectories visit regions of the S2 potential corresponding to the planar πSπ456* minimum, with only few of them accessing C2-pyramidalized regions of the PES, characteristic of the πSπ2* minimum. In contrast, trajectories were found to distribute between the T1 boat conformation and the C2-pyramidalized T1 minima, the second carrying a larger population. The simulations are not able to properly describe the repopulation of the GS from the triplet manifold, due to the very small number of trajectories reaching the GS at the end of the simulation, which is also consistent with the experimental results.
Finally, intersystem crossing and singlet oxygen quantum yields have been also estimated for this family of compounds. The very high triplet quantum yields for 2t-Ura (0.75) [60] and 2t-Thy nucleobase and nucleotides (0.9–1) [62,63,64] contrast with the moderate singlet oxygen yields 0.36 [62,63,64] measured for them (See Table 3). The potential competition between energy transfer and charge transfer complexes has been suggested as one of the most plausible arguments for explaining such modest singlet oxygen yields [63].

2.2.2. 4-Thiopyrimidines

4t-Thy/Thd and 4t-Ura/Urd

• Steady State Absorption and Emission properties
The experimental absorption spectrum of 4t-Ura shows a main band in the 300–400 nm region with a maximum at ca. 328 nm (3.78 eV) in water and in ACN (see Table 6) [61,84,99]. The comparison with experiments performed on 4t-Urd in water solution reveals that the presence of the sugar moiety slightly red-shifts the maximum to ~330 nm, being again not affected by the change of the solvent (ACN) [65,72,84,100,101]. Similar absorption spectra were recorded for thiouracil methyl-derivatives (dm-4tUra) [65,102,103]. The same trends have been reported for the absorption spectra of 4t-Thd, i.e., with a maximum at 335 nm (3.70 eV) [64,67,103] in PBS, which is almost unaffected by the change of the solvent (ACN, 337 nm) [66,68].
Different computational approaches have been applied to predict the absorption spectrum of 4-thioderivatives (Table 7). Dynamically correlated methods (MCQDPT2) deliver for 4t-Ura an energy for the first ππ* state of 3.90 eV (318 nm), which is on top of the experimental value [104]. In contrast, higher energies 298–294 nm (4.16–4.22eV) are provided by monoconfigurational methods, such as TD-B3LYP and EOM-CC2 calculations [31,99,103]. Notably, other more restrictive semi-empirical methods such as INDO/S also find energies for the bright states of 4t-Ura and dm-4tUra close to experimental maximum [65].
The absorption spectra of 4t-Thd and 4t-Thy were also investigated from a computational point of view, both in GP and in solution. MS-CASPT2 calculates the spectroscopic state of 4t-Thy (S2) at 3.87–4.00 eV (320–310 nm) with a ππ* character [83,93]. Including explicit water molecules in the MS-CASPT2 calculation leads to a red-shift of the S2 energy by −0.07 eV (326 nm) [83]. QM(CASPT2)/MM calculations on 4t-Thd provide the same assignment of the bright state, also absorbing at 4.11 eV (302 nm) [95].
The spectroscopic state of 4t-Thy/4t-Thd is still S2 at the TD-B3LYP level of theory and involves an excitation from the HOMO to the LUMO orbitals (Figure 1). The energies of this state range from ~4.3 eV (~290 nm) in the GP to 4.15 eV (~300 nm) in ACN (PCM calculations) and decrease to 312 nm (3.97 eV) with DFT/MRCI or when the M06 functional is used [31,37,66,67,69,93,99].
All theses computational studies characterize the first excited state (S1) of 4t-Thy as a dark state which promotes an electron from the LP of the S atom to the LUMO orbital (Table 7).
• Static description of the PES
A summary of the main deactivation pathways described for 4-thioderivatives is depicted in Figure 4 (top right). Upon excitation of 4t-Thy/4t-Thd at ~330 nm, the spectroscopic state S2 (ππ*), as described above, is mainly populated and its deactivation mechanism has inspired several theoretical studies. Two different deactivation routes have been proposed based on MEPs performed on 4t-Thy coordinated to several explicit water molecules at the MS-CASPT2/CASSCF level of theory [83]. Both pathways converge to singlet ππ* minima, characterized by markedly different properties in terms of energy and geometry. The first S2P minimum exhibits a planar conformation of the pyrimidine ring and lies lower in energy (3.49 eV), whereas the second S2T is distorted along the C5=C6 double bond (recall labeling in Scheme 1) and lies much higher in energy (4.34 eV) [83]. A third singlet minimum was found in this study along the S1 PES presenting nπ* character (2.53 eV) [83]. Accessing the nπ* minimum requires, however, reaching the ππ*/nπ* IC funnel (4.47 eV), the efficiency of this pathway being dictated by the distribution of the trajectories between the two ππ* minima: S2P is well-separated from the CoIn and should be stable against its depopulation; in contrast, the proximity of S2T to the IC funnel indicates that this minimum is prone to decay to lower-lying singlet minimum [83]. Excited state optimizations of 4t-Thd at the QM(CASSCF)/MM level of theory localize, instead, two singlet excited minima, the S2-MIN and S1-MIN placed 3.11 eV and 2.54 eV relative to the GS equilibrium minimum, respectively [95]. The latter resembles the nπ* minimum, but the former seems not to correspond to any of the (ππ*) minima described above, i.e., it is geometrically close to the S2/S1 IC funnel in agreement with the CASPT2 S2T minimum, but its low energy is in line with that of the S2P minimum [95].
According to the most accurate quantum mechanical calculations, the S2 PES of 4t-Thy/4t-Thd presents a complicated topography governed by several local minima, which are expected to decay to lower-lying states within different time scales.
Some experimental hints support these theoretical findings. First, although small (<1 × 10−4), 4t-Thd shows a measurable fluorescence quantum yield, which is consistent with the presence of a quite stable singlet minimum that is partially populated [4,67,68,83]. The nature of this emissive state has been intensively debated. Reichardt et al. [67] first assigned it to either S1 or S2 states, whereas later Wenska et al. [103] suggested that fluorescence should come from the S2 state, consistently with a non-Kasha behavior. Recently, Martínez-Fernández et al. pointed the S2P minimum as the responsible for this fluorescence, based on the large energy required for the deactivation from this state [83]. Importantly, up-conversion experiments reported a bi-exponential decay for the fluorescence, compatible with the intricate topography of the S2 PES governed by different excited minima [4,83].
The fluorescence yield for 4-thiopyrimidines is very small and most of the population is expected to non-radiatively decay from the S2 ππ* to the dark S1 nπ* minimum. Once in this minimum, deactivation to the GS seems improbable since all computational studies find a large activation energy to reach the S1/S0 funnel (1.5–2.0 eV) [83,95]. As an alternative, efficient population of the triplet manifold from the nπ* minimum has been proposed due to the presence of singlet/triplet MECP near the nπ* minimum. MS-CASPT2/CASSCF MEPs locate the S1(nπ*)/T1(ππ*) crossing at the position of the S1 (nπ*) minimum and the QM(CASPT2/CASSCF)/MM results within less than 1 kcal·mol−1 from it [83,95]. Furthermore, the calculated SOCs at this region of the PES are large 60–150 cm−1 reinforcing the hypothesis that considers the nπ* minimum as a doorway for triplet population [69,83,93,95]. Additional singlet/triplet funnels, between the S2 (ππ*) and T2 (nπ*) states, have been also localized close in energy to the S1/S2 IC funnel [69,83,95]. The SOC at this point was calculated to be of the same order of magnitude (90–120 cm−1) as for the S1(nπ*)/T1(ππ*) crossing, which a priori prevents discarding deactivation from this singlet-triplet funnel, at least from a static point of view (see below) [69,83,95]. Independently of the proposed route for triplet population, all the studies conclude that the most stable T1 (ππ*) state (2.28–2.34 eV) will be finally populated and that its decay to the GS is not energetically favored [83,95]. Consistently with the existence of these stable triplet minima, 4t-Thd does present some phosphorescence, although the yield is small 3 × 10−4 [66,67,101,103]. The phosphorescence spectrum has a maximum at ~540 nm [66,69,101,102,103], matching the calculated energies for the T1 (ππ*) minimum both at TD-DFT (538–557 nm) [66,67,69] and CASPT2 (570 nm) levels of theory [83].
• Photodynamics
Time resolved experiments on 4t-Thd have essentially focused on predicting the ISC rates and yields in different solvents. The experimental ISC rates greatly vary depending on the solvent and experimental conditions (i.e., excitation wavelength, laser resolution, etc.). In buffer solution, experimental ISC lifetimes range from 250 fs [67,83] to 10 ps [69], whereas under the same experimental conditions they are longer in ACN, 540 fs (see Table 3) [68,83]. Interestingly, however, all the experiments uniformly determine a quantum yield for T1 population close to unity and a very long lifetime, from tens of ps to μs depending on the solvent, for this triplet state [66,67,68,83,101,102,103].
As far as the TAS are concerned, Reichardt et al. [67,68] discuss two peaks growing simultaneously at ~380 and 550 nm both in water and in ionic liquids, the latter also observed by Harada et al. [69]. The TAS signal peaking at ~550 nm was ascribed to the final triplet state, due to its resemblance with the triplet-triplet absorption band. Importantly, just recently an additional spectral feature in the 650–700 nm region, which was before overlooked and which decays in few ps, has been assigned to one of the two S2 minima, revealing important details of the relaxation mechanism [83].
From a computational point of view, the experimental τISC = 250 fs in water is well reproduced by a simple qualitative approach, which by considering the SOC terms for the S1(nπ*)/T1(ππ*) crossing calculates a τISC for 4t-Thy of 275 fs. [93]. However, the very different τISC values experimentally registered suggest a more complex ISC mechanism, also in agreement with the already mentioned bi-exponential behavior of the fluorescence, (τ1 = 1 ps and τ2 = 3.6 ps) [83]. In this respect, QM/MM MD performed on 4t-Thd in solution shows that upon excitation the system actually explores regions of the PES corresponding to the planar and the twisted S2 minima [83]. It is actually not possible to distinguish the trajectories following the paths related to either one of these two minima, i.e., all of them seem to explore both regions of the S2 PES. These dynamics simulations are, however, able to discern the most likely ISC channel, the S1(nπ*)/T1(ππ*), followed by 87% of the trajectories.
This mechanism is actually supported by the results from the experimental TAS, which were interpreted with the help of MS-CASPT2 vertical spectra at the position of excited state minima. The S2T minimum was found to absorb at 600 nm, whereas the S2P minimum shows a maximum above 650 nm, coinciding with the spectral feature in the experimental TAS recently rediscovered at long wavelengths. Indeed, Martínez-Fernández et al. [83] showed that in order to reproduce the experimental TAS at 0.8 ps a significant population (~40%) must be still present at these minima. This assignment confirms the results of the MD simulations, which situate the photoexcited molecules in planar and twisted regions of the S2 PES along their relaxation to the triplet manifold (see Figure 4) [83]. Experimentally two different τISC rates; one in the ultrafast timescale (τ1 = 0.54 ps) and another one in the picosecond regime (τ2 = 1.8 ps) were registered (see Figure 3) [83]. The long-lived signal at 550 nm was ascribed to a superposition of the two triplet minima.
In short, the deactivation mechanism leading to triplet population in 4t-Thd is complex. Indeed, the intricacy of the PES of 4-thionucleobases is responsible for the modulation of the population rate of the T1 minimum and thus of the different ISC lifetimes experimentally registered.

2.2.3. 2,4-Dithiopyrimidines

Steady State Absorption and Emission properties

The experimental absorption spectra of doubly-thionated bases (Table 8) are characterized by maxima peaking at ~360–340 nm (3.44–3.65 eV) independently of the solvent used, which supposes a red-shift compared to 4-thiopyrimidines (~320 nm) [61,64,84]. Theoretically (Table 9), the first bright state has been predicted in GP at ~330 nm with MS-CASPT2 [93], ~340 nm with DFT/MRCI [93], and 349 nm with TD-B3LYP [31], whereas in water solution the computed wavelengths are ~330 nm at MS-CASPT2 [96], 335 nm at TD-B3LYP [31] and 340 nm at the TD-M06 level of theory [37]. Equivalent calculations for 4t-Thy provided the following values 310, 312, 300 and 311 nm, all the methods being able to reproduce the red-shift of the first absorption band when both carbonyl groups are exchanged by thiocarbonyl groups [31,37,93,95]. The analysis of the molecular orbitals in the series 2t-Thy, 4t-Thy, and 2,4-dtThy performed by Bai et al. concludes that the number and position of sulfur functionalization in the pyrimidine ring leads to different electron density delocalization conditioning the orbital energies and thus the spectra [93].

Static description of the PES

The enhanced photosensitizing properties of 2,4-dtUra/2,4-dtThy [61,64] compared to their monosubstituted counterparts stimulated other theoretical studies, mostly focused on how triplet population occurs in these nucleobases. Following the population of the spectroscopic S2S2π4*) state (3.62–3.72 eV) in GP, MS-CASPT2 studies (based on CASSCF geometry optimizations) predict 2,4-dtThy to undergo either adiabatic relaxation to its minimum or a barrierless IC towards low-lying S1 minimum (nS2π2*), (See Figure 4, bottom) [96]. The population following the first route and reaching the ππ* minimum is expected to eventually relax to the 1nS2π2* minimum at some stage [96]. Due to the proximity of the 1nS4π4*/1nS2π2* crossing to the 1nS2π2* minimum, the authors suggest that another dark minimum, presenting 1nS4π4* character, can also be populated [96]. Access to the GS via IC along the singlet manifold is unfortunately not explored in this work.
At the position of both 1nS2π2* and 1nS4π4* minima, MS-CASPT2//CASSCF calculates very close ISC funnels to the triplets 3πS2π2* and 3πS4π4*, the latter being strongly coupled (SOC amounting to 154 cm−1) to the 1nS4π4* state at its minimum position [96]. This value is in agreement with other SOC values obtained in GP MS-CASPT2 and MRCI studies, 152 cm−1 [93].
Xie et al. [96] paid special attention to the effect that microsolvation and aqueous solution environments have in the computed deactivation pathways of 2,4-dtThy. To this aim, they performed QM(CASSCF)/MM geometry optimizations followed by MS-CASPT2 single points in order to describe the potential energy profiles. At the FC region, water was found to destabilize the (1nS2π2*) dark state, decreasing the energy gap with the bright (1πS2π4*) spectroscopic state, thus, facilitating IC between both states. Then the authors suggest that both the 1nS2π2* and 1πS2π4* minima are populated in solution from the FC region [96]. From either singlet minima, Xie et al. [96] set out an intricate network of possible mechanisms on the way to the most stable 3πS4π4* minimum. However, in contrast to their GP results, in solution the population reaching the 1πS2π4* minimum is not assumed to later decay to the 1nS2π2* singlet minimum, but only to the triplet states. SOCs are not much affected by solvent [96], in agreement with Ref. [93]. Thus, water solvation was found to “entangle” the mechanism leading to the population of the 3πS4π4* minimum, but it is not expected to affect the final outcome of the dynamics.

Photodynamics

From an experimental point of view, doubly-thionated pyrimidines show very rich photophysical properties such as faster ISC rates (180 fs, 2,4-dtThy) and higher singlet oxygen quantum yield (0.49, 2,4-dtUra) compared to their monosubstituted counterparts (see Table 3) [61,64].
To the best of our knowledge, no molecular dynamics simulations are available in the literature for 2,4-dtThy/2,4-dtUra, except for a single study calculating ISC rates for 2,4-dtThy based on SOC terms [93]. The predicted ISC rate is 295 fs, which is significantly faster than the one calculated for 2t-Thy (574 fs), in agreement with the experimental findings [93].
Theoretical works have reported energetically accessible ISC funnels, large SOCs and stable triplet minima for 2,4-dtThy [93,96]. These features, however, are also fingerprints of the mono-substituted thiobases’ PES, what leaves open the question as to why double-thionation leads to shorter ISC lifetimes and larger 1O2 yields. A complete understanding of the photophysics of doubly-thionated pyrimidines would benefit from comprehensive molecular dynamics studies considering both singlet and triplet manifolds, and including TAS interpretation by means of theoretical calculations.

3. General Remarks and Perspectives

The main effects induced by thionation in the natural DNA and RNA nucleobases’ photophysics are nowadays well established. First, exocyclic oxygen by sulfur substitution significantly red-shifts the maximum of the absorption spectrum. In contrast to natural DNA and RNA nucleobases, which are potent UVB chromophores, thiated nucleobases absorb maximally in the UVA region of the electromagnetic spectrum, (see Figure 5) being this an important side effect in pharmacological applications due to the large amount of UVA photons that reach the Earth’s surface.
Interestingly, experiments and theoretical simulations have shown that for pyrimidines the magnitude of this shift depends on the position and degree substitution, being modest in case of 2-thiopyrimidines but meaningful for 4- and 2,4-dithioderivatives [61,64]. Although calculations were found to overestimate the shifts in the absorption bands for 2-thiopyrimidines, they have provided useful hints to understand the changes in the absorption spectra from an electronic point of view. Bai et al. [93] actually succeeded to explain the observed spectral shifts on the basis of the different delocalization of the electron density provoked by the degree and pattern of substitution and the changes in the energies of the molecular orbitals involved in the absorption bands of 2t-Thy, 4t-Thy and 2,4-dtThy. It is also important to note that these shifts do not affect to the same extent all the electronic states and can therefore lead to a change of the electronic state ordering compared to canonical nucleobases. In fact, the lowest lying electronic state in all the thioderivatives presented in this review corresponds to a dark (S1 nπ*), something not common in canonical nucleobases [42].
Importantly, the exchange of one or several carbonyl groups by thiocarbonyl groups in natural nucleobases undoubtedly conditions the topography of their PES and thus of their photoinitiated dynamics. As a consequence of the important shift of the first nπ* state towards lower energies and below the spectroscopic ππ* electronic state, the first step in the relaxation mechanisms of all canonical nucleobases’ thiated analogs is an IC process towards the first excited PES. This deactivation channel has been shown to be preferred to direct ISC from the S2. The main differences between thiated and natural nucleobases lie on the subsequent relaxation from S1. Even though recent studies in solution have located small energy barriers associated to GS deactivation from the S1, it has been generally assumed that this is the most favorable mechanism in canonical nucleobases [42,106]. By contrast, this deactivation route is strongly prevented in sulfur-substituted nucleobases since the S1/S0 funnels are energetically inaccessible due to the significant stabilization of the S1 minima (Figure 6). Furthermore, another secondary effect of sulfur substitution is that other S1/S0 CoIn different to the preferred ones for canonical nucleobases may come into play [34]. Besides, oxygen by sulfur substitution also provokes a decrease in the singlet-triplet energy gaps and the enhancement of the SOC (thiobases 100–200 cm−1; canonical bases < 50 cm−1) due to the heavy atom effect of S. Both effects are translated into energetically accessible and very efficient singlet-triplet crossings at the vicinity of the S1 minima of thioderivatives, which correspond to secondary routes in canonical nucleobases (Figure 6). Indeed, ultrafast ISCs and close to unity triplet yields were experimentally observed and inferred from MD simulations [4,34,58,83,98]. Although Pollum et al. [61,64] have registered enhanced photosensitizing properties in doubly thionated nucleobases, the reasons for their faster ISC and closer to unity triplet yields are still unclear since singly- and doubly-substituted thiobases share the main features of the PES. The long life of their most stable triplet state, characterized in the case of 2t-Thy, 2t-Ura, 6t-Gua, 2t-Cyt by double well potential energy profile, is also common to all the thiobases [55,59,89].
Thiobases are known to be important precursors of 1O2 in solution via Type II photosensitizing mechanisms, with yields comprised between 0.14 and 0.58 (see Table 3) [4,73]. The significant differences between the triplet and singlet oxygen quantum yields suggest that other mechanisms, such as charge transfer reactions with O2, might compete in solution [62,63,65,73,83]. In cellular or biological environments, the situation is much more complex since Type I and Type II mechanisms can involve a broad spectrum of biomolecules leading to UVA induced lesions including canonical nucleobases’, 6t-Gua’s and protein’s oxidation products, DNA and DNA–protein crosslinks and DNA strand breaks [16,79,82,107].
Although the basis for the first stages of the photosensitizing reactions of thiated nucleobases that lead to triplet excited states’ population have been clearly established by experimental and computational works, a great piece of work still remains to be done for the characterization of the mechanisms at molecular level leading to photoproducts’ formation, which are ultimately responsible for the mechanism of action of these photosensitizers in medicinal chemistry.

Acknowledgments

The authors acknowledge financial support from the MINECO (FPI grant (SA), and projects FIS2013-42002-R and CTQ2015-63997-C2), the Comunidad de Madrid (project FOTOCARBON-CM S2013/MIT-2841), and the European COST Action XLIC CM1204.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Crespo-Hernández, C.E.; Cohen, B.; Hare, P.M.; Kohler, B. Ultrafast excited-state dynamics in nucleic acids. Chem. Rev. 2004, 104, 1977–2020. [Google Scholar] [CrossRef] [PubMed]
  2. Barbatti, M.; Borin, A.C.; Ullrich, S. Photoinduced phenomena in nucleic acids I, II. In Topics in Current Chemistry; Springer: Berlin, Germany, 2014. [Google Scholar]
  3. Gustavsson, T.; Improta, R.; Markovitsi, D. DNA/RNA: Building blocks of life under UV irradiation. J. Phys. Chem. Lett. 2010, 1, 2025–2030. [Google Scholar] [CrossRef]
  4. Pollum, M.; Martínez-Fernández, L.; Crespo-Hernández, C.E. Photochemistry of nucleic acid bases and their thio- and aza-analogues in solution. In Photoinduced Phenomena in Nucleic Acids I: Nucleobases in the Gas Phase and in Solvents; Barbatti, M., Borin, A.C., Ullrich, S., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 245–327. [Google Scholar]
  5. Crespo-Hernández, C.E.; Martínez-Fernández, L.; Rauer, C.; Reichardt, C.; Mai, S.; Pollum, M.; Marquetand, P.; González, L.; Corral, I. Electronic and structural elements that regulate the excited-state dynamics in purine nucleobase derivatives. J. Am. Chem. Soc. 2015, 137, 4368–4381. [Google Scholar] [CrossRef] [PubMed]
  6. Matsika, S. Modified nucleobases. In Photoinduced Phenomena in Nucleic Acids I: Nucleobases in the Gas Phase and in Solvents; Barbatti, M., Borin, A.C., Ullrich, S., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 209–243. [Google Scholar]
  7. Szostak, J.W. The eightfold path to non-enzymatic RNA replication. J. Syst. Chem. 2012, 3, 2. [Google Scholar] [CrossRef]
  8. Zhang, S.; Blain, J.C.; Zielinska, D.; Gryaznov, S.M.; Szostak, J.W. Fast and accurate nonenzymatic copying of an RNA-like synthetic genetic polymer. Proc. Natl. Acad. Sci. USA 2013, 110, 17732–17737. [Google Scholar] [CrossRef] [PubMed]
  9. Caton-Williams, J.; Huang, Z. Biochemistry of selenium-derivatized naturally occurring and unnatural nucleic acids. Chem. Biodivers. 2008, 5, 396–407. [Google Scholar] [CrossRef] [PubMed]
  10. Ajitkumar, P.; Cherayil, J.D. Thionucleosides in transfer ribonucleic acid: Diversity, structure, biosynthesis, and function. Microbiol. Rev. 1988, 52, 103–113. [Google Scholar] [PubMed]
  11. Aarbakke, J.; Janka-Schaub, G.; Elion, G.B. Thiopurine biology and pharmacology. Trends Pharmacol. Sci. 1997, 18, 3–7. [Google Scholar] [CrossRef]
  12. Relling, M.V.; Dervieux, T. Pharmacogenetics and cancer therapy. Nat. Rev. Cancer 2001, 1, 99–108. [Google Scholar] [CrossRef] [PubMed]
  13. Warren, D.J.; Andersen, A.; Slørdal, L. Quantitation of 6-thioguanine residues in peripheral blood leukocyte DNA obtained from patients receiving 6-mercaptopurine-based maintenance therapy. Cancer Res. 1995, 55, 1670–1674. [Google Scholar] [PubMed]
  14. Cuffari, C.; Li, D.Y.; Mahoney, J.; Barnes, Y.; Bayless, T.M. Peripheral blood mononuclear cell DNA 6-thioguanine metabolite levels correlate with decreased interferon-γ production in patients with crohn’s disease on aza therapy. Dig. Dis. Sci. 2004, 49, 133–137. [Google Scholar] [CrossRef] [PubMed]
  15. O’Donovan, P.; Perrett, C.M.; Zhang, X.; Montaner, B.; Xu, Y.-Z.; Harwood, C.A.; McGregor, J.M.; Walker, S.L.; Hanaoka, F.; Karran, P. Azathioprine and UVA light generate mutagenic oxidative DNA damage. Science 2005, 309, 1871–1874. [Google Scholar] [CrossRef] [PubMed]
  16. Brem, R.; Karran, P. Multiple forms of DNA damage caused by UVA photoactivation of DNA 6-thioguanine. J. Photochem. Photobiol. 2012, 88, 5–13. [Google Scholar] [CrossRef] [PubMed]
  17. Favre, A.; Saintomé, C.; Fourrey, J.-L.; Clivio, P.; Laugâa, P. Thionucleobases as intrinsic photoaffinity probes of nucleic acid structure and nucleic acid-protein interactions. J. Photochem. Photobiol. B 1998, 42, 109–124. [Google Scholar] [CrossRef]
  18. Meisenheimer, K.M.; Koch, T.H. Photocross-linking of nucleic acids to associated proteins. Crit. Rev. Biochem. Mol. Biol. 1997, 32, 101–140. [Google Scholar] [CrossRef] [PubMed]
  19. Favre, A.; Yaniv, M.; Michels, A.M. The photochemistry of 4-thiouridine in escherichia coli t-RNAİVal. Biochem. Biophys. Res. Commun. 1969, 37, 266–271. [Google Scholar] [CrossRef]
  20. Favre, A.; Moreno, G.; Blondel, M.O.; Kliber, J.; Vinzens, F.; Salet, C. 4-thiouridine photosensitized RNA-protein crosslinking in mammalian cells. Biochem. Biophys. Res. Commun. 1986, 141, 847–854. [Google Scholar] [CrossRef]
  21. Favre, A.; Bezerra, R.; Hajnsdorf, E.; Dubreuil, Y.L.; Expert-Bezançon, A. Substitution of uridine in vivo by the intrinsic photoactivable probe 4-thiouridine in escherichia coli RNA. Eur. J. Biochem. 1986, 160, 441–449. [Google Scholar] [CrossRef] [PubMed]
  22. Sontheimer, E.J. Site-specific RNA crosslinking with 4-thiouridine. Mol. Biol. Rep. 1994, 20, 35–44. [Google Scholar] [CrossRef] [PubMed]
  23. Harris, M.E.; Christian, E.L. Chapter 7—RNA crosslinking methods. In Methods Enzymol.; Academic Press: Amsterdam, The Netherlands, 2009; Volume 468, pp. 127–146. [Google Scholar]
  24. Kumar, R.K.; Davis, D.R. Synthesis and studies on the effect of 2-thiouridine and 4-thiouridine on sugar conformation and RNA duplex stability. Nucleic Acids Res. 1997, 25, 1272–1280. [Google Scholar] [CrossRef] [PubMed]
  25. Smith, W.S.; Sierzputowska-Gracz, H.; Sochacka, E.; Malkiewicz, A.; Agris, P.F. Chemistry and structure of modified uridine dinucleosides are determined by thiolation. J. Am. Chem. Soc. 1992, 114, 7989–7997. [Google Scholar] [CrossRef]
  26. Sierzputowska-Gracz, H.; Sochacka, E.; Malkiewicz, A.; Kuo, K.; Gehrke, C.W.; Agris, P.F. Chemistry and structure of modified uridines in the anticodon, wobble position of transfer RNA are determined by thiolation. J. Am. Chem. Soc. 1987, 109, 7171–7177. [Google Scholar] [CrossRef]
  27. Vormbrock, R.; Morawietz, R.; Gassen, H.G. Codon—Anticodon interaction studied with trinucleoside diphosphates containing 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, or 2-thiocytidine. Biochim. Biophys. Acta 1974, 340, 348–358. [Google Scholar] [CrossRef]
  28. Park, E.; Baron, R.; Landgraf, R. Higher-order association states of cellular ERBB3 probed with photo-cross-linkable aptamers. Biochemistry 2008, 47, 11992–12005. [Google Scholar] [CrossRef] [PubMed]
  29. Les, A.; Adamowicz, L. Tautomerism of 2- and 4-thiouracil. Ab initio theoretical study. J. Am. Chem. Soc. 1990, 112, 1504–1509. [Google Scholar] [CrossRef]
  30. Rostkowska, H.; Szczepaniak, K.; Nowak, M.J.; Leszczynski, J.; KuBulat, K.; Person, W.B. Thiouracils. 2. Tautomerism and infrared spectra of thiouracils. Matrix-isolation and ab initio studies. J. Am. Chem. Soc. 1990, 112, 2147–2160. [Google Scholar] [CrossRef]
  31. Shukla, M.K.; Leszczynski, J. Electronic transitions of thiouracils in the gas phase and in solutions:  Time-dependent density functional theory (TD-DFT) study. J. Phys. Chem. A 2004, 108, 10367–10375. [Google Scholar] [CrossRef]
  32. Puzzarini, C.; Biczysko, M.; Barone, V.; Peña, I.; Cabezas, C.; Alonso, J.L. Accurate molecular structure and spectroscopic properties for nucleobases: A combined computational—Microwave investigation of 2-thiouracil as a case study. Phys. Chem. Chem. Phys. 2013, 15, 16965–16975. [Google Scholar] [CrossRef] [PubMed]
  33. Podolyan, Y.; Gorb, L.; Blue, A.; Leszczynski, J. A theoretical investigation of tautomeric equilibria and proton transfer in isolated and hydrated thiocytosine. J. Mol. Struct. THEOCHEM 2001, 549, 101–109. [Google Scholar] [CrossRef]
  34. Mai, S.; Pollum, M.; Martínez-Fernández, L.; Dunn, N.; Marquetand, P.; Corral, I.; Crespo-Hernández, C.E.; González, L. The origin of efficient triplet state population in sulfur-substituted nucleobases. Nat. Commun. 2016, 7, 13077. [Google Scholar] [CrossRef] [PubMed]
  35. Leszczynski, J. Tautomers of 6-thioguanine: Structures and properties. J. Phys. Chem. 1993, 97, 3520–3524. [Google Scholar] [CrossRef]
  36. Civcir, P.U. Tautomerism of 6-thioguanine in the gas and aqueous phases using AM1 and PM3 methods. J. Mol. Struct. THEOCHEM 2001, 536, 161–171. [Google Scholar] [CrossRef]
  37. Pirillo, J.; De Simone, B.C.; Russo, N. Photophysical properties prediction of selenium- and tellurium-substituted thymidine as potential UVA chemotherapeutic agents. Theor. Chem. Acc. 2016, 135, 8–12. [Google Scholar] [CrossRef]
  38. Pirillo, J.; Mazzone, G.; Russo, N.; Bertini, L. Photophysical properties of S, Se and Te-substituted deoxyguanosines: Insight into their ability to act as chemotherapeutic agents. J. Chem. Inf. Model. 2017, 57, 234–242. [Google Scholar] [CrossRef] [PubMed]
  39. Yu, H.; Sanchez-Rodriguez, J.A.; Pollum, M.; Crespo-Hernández, C.E.; Mai, S.; Marquetand, P.; González, L.; Ullrich, S. Internal conversion and intersystem crossing pathways in UV excited, isolated uracils and their implications in prebiotic chemistry. Phys. Chem. Chem. Phys. 2016, 18, 20168–20176. [Google Scholar] [CrossRef] [PubMed]
  40. Ruckenbauer, M.; Mai, S.; Marquetand, P.; González, L. Photoelectron spectra of 2-thiouracil, 4-thiouracil, and 2,4-dithiouracil. J. Chem. Phys. 2016, 144, 074303–074310. [Google Scholar] [CrossRef] [PubMed]
  41. González, L.; Escudero, D.; Serrano-Andrés, L. Progress and challenges in the calculation of electronic excited states. Chem. Phys. Chem. 2012, 13, 28–51. [Google Scholar] [CrossRef] [PubMed]
  42. Improta, R.; Santoro, F.; Blancafort, L. Quantum mechanical studies on the photophysics and the photochemistry of nucleic acids and nucleobases. Chem. Rev. 2016, 116, 3540–3593. [Google Scholar] [CrossRef] [PubMed]
  43. Ruckebusch, C.; Sliwa, M.; Pernot, P.; de Juan, A.; Tauler, R. Comprehensive data analysis of femtosecond transient absorption spectra: A review. J. Photochem. Photobiol. C 2012, 13, 1–27. [Google Scholar] [CrossRef]
  44. Stolow, A.; Bragg, A.E.; Neumark, D.M. Femtosecond time-resolved photoelectron spectroscopy. Chem. Rev. 2004, 104, 1719–1758. [Google Scholar] [CrossRef] [PubMed]
  45. Efremov, E.V.; Ariese, F.; Gooijer, C. Achievements in Resonance Raman spectroscopy: Review of a technique with a distinct analytical chemistry potential. Anal. Chim. Acta 2008, 606, 119–134. [Google Scholar] [CrossRef] [PubMed]
  46. Thiel, W. Semiempirical quantum–chemical methods. WIREs Comput. Mol. Sci. 2014, 4, 145–157. [Google Scholar] [CrossRef]
  47. Persico, M.; Granucci, G. An overview of nonadiabatic dynamics simulations methods, with focus on the direct approach versus the fitting of potential energy surfaces. Theor. Chem. Acc. 2014, 133, 1–28. [Google Scholar] [CrossRef]
  48. Lasorne, B.; Worth, G.A.; Robb, M.A. Excited-state dynamics. WIREs Comput. Mol. Sci. 2011, 1, 460–475. [Google Scholar] [CrossRef]
  49. Barbatti, M. Nonadiabatic dynamics with trajectory surface hopping method. WIREs Comput. Mol. Sci. 2011, 1, 620–633. [Google Scholar] [CrossRef]
  50. Rubin, Y.V.; Blagoi, Y.P.; Bokovoy, V.A. 6-thioguanine luminescence probe to study DNA and low-molecular-weight systems. J. Fluoresc. 1995, 5, 263–272. [Google Scholar] [CrossRef] [PubMed]
  51. Santhosh, C.; Mishra, P.C. Electronic structures and spectra of 6-mercaptopurine and 6-thioguanine. Spectrochim. Acta Part A 1993, 49, 985–993. [Google Scholar] [CrossRef]
  52. Ashwood, B.; Jockusch, S.; Crespo-Hernández, C.E. Excited-state dynamics of the thiopurine prodrug 6-thioguanine: Can N9-glycosylation affect its phototoxic activity? Molecules 2017, 22, 379–393. [Google Scholar] [CrossRef] [PubMed]
  53. Gomzi, V. TDDFT study of nucleobase thioanalogues and oxo-derivatives excited states. J. Theor. Comput. Chem. 2009, 08, 71–83. [Google Scholar] [CrossRef]
  54. Reichardt, C.; Guo, C.; Crespo-Hernández, C.E. Excited-state dynamics in 6-thioguanosine from the femtosecond to microsecond time scale. J. Phys. Chem. B 2011, 115, 3263–3270. [Google Scholar] [CrossRef] [PubMed]
  55. Martínez-Fernández, L.; González, L.; Corral, I. An ab initio mechanism for efficient population of triplet states in cytotoxic sulfur substituted DNA bases: The case of 6-thioguanine. Chem. Commun. 2012, 48, 2134–2136. [Google Scholar] [CrossRef] [PubMed]
  56. Zhang, Y.; Zhu, X.; Smith, J.; Haygood, M.T.; Gao, R. Direct observation and quantitative characterization of singlet oxygen in aqueous solution upon UVA excitation of 6-thioguanines. J. Phys. Chem. B 2011, 115, 1889–1894. [Google Scholar] [CrossRef] [PubMed]
  57. Chuan, D.; Wen, Y.; Shaomin, S.; Pin, Y. Determination of thioguanine in pharmaceutical preparations by paper substrate room temperature phosphorimetry. Analyst 2000, 125, 1327–1330. [Google Scholar] [CrossRef] [PubMed]
  58. Martínez-Fernández, L.; Corral, I.; Granucci, G.; Persico, M. Competing ultrafast intersystem crossing and internal conversion: A time resolved picture for the deactivation of 6-thioguanine. Chem. Sci. 2014, 5, 1336–1347. [Google Scholar] [CrossRef]
  59. Bai, S.; Barbatti, M. On the decay of the triplet state of thionucleobases. Phys. Chem. Chem. Phys. 2017, 19, 12674–12682. [Google Scholar] [CrossRef] [PubMed]
  60. Pollum, M.; Crespo-Hernández, C.E. Communication: The dark singlet state as a doorway state in the ultrafast and efficient intersystem crossing dynamics in 2-thiothymine and 2-thiouracil. J. Chem. Phys. 2014, 140, 071101–071105. [Google Scholar] [CrossRef] [PubMed]
  61. Pollum, M.; Jockusch, S.; Crespo-Hernández, C.E. Increase in the photoreactivity of uracil derivatives by doubling thionation. Phys. Chem. Chem. Phys. 2015, 17, 27851–27861. [Google Scholar] [CrossRef] [PubMed]
  62. Kuramochi, H.; Kobayashi, T.; Suzuki, T.; Ichimura, T. Excited-state dynamics of 6-aza-2-thiothymine and 2-thiothymine: Highly efficient intersystem crossing and singlet oxygen photosensitization. J. Phys. Chem. B 2010, 114, 8782–8789. [Google Scholar] [CrossRef] [PubMed]
  63. Taras-Goślińska, K.; Burdziński, G.; Wenska, G. Relaxation of the T1 excited state of 2-thiothymine, its riboside and deoxyriboside-enhanced nonradiative decay rate induced by sugar substituent. J. Photochem. Photobiol. A 2014, 275, 89–95. [Google Scholar] [CrossRef]
  64. Pollum, M.; Jockusch, S.; Crespo-Hernández, C.E. 2,4-dithiothymine as a potent UVA chemotherapeutic agent. J. Am. Chem. Soc. 2014, 136, 17930–17933. [Google Scholar] [CrossRef] [PubMed]
  65. Milder, S.J.; Kliger, D.S. Spectroscopy and photochemistry of thiouracils: Implications for the mechanism of photocrosslinking in tRNA. J. Am. Chem. Soc. 1985, 107, 7365–7373. [Google Scholar] [CrossRef]
  66. Harada, Y.; Suzuki, T.; Ichimura, T.; Xu, Y.-Z. Triplet formation of 4-thiothymidine and its photosensitization to oxygen studied by time-resolved thermal lensing technique. J. Phys. Chem. B 2007, 111, 5518–5524. [Google Scholar] [CrossRef] [PubMed]
  67. Reichardt, C.; Crespo-Hernández, C.E. Room-temperature phosphorescence of the DNA monomer analogue 4-thiothymidine in aqueous solutions after UVA excitation. J. Phys. Chem. Lett. 2010, 1, 2239–2243. [Google Scholar] [CrossRef]
  68. Reichardt, C.; Crespo-Hernández, C.E. Ultrafast spin crossover in 4-thiothymidine in an ionic liquid. Chem. Commun. 2010, 46, 5963–5965. [Google Scholar] [CrossRef] [PubMed]
  69. Harada, Y.; Okabe, C.; Kobayashi, T.; Suzuki, T.; Ichimura, T.; Nishi, N.; Xu, Y.-Z. Ultrafast intersystem crossing of 4-thiothymidine in aqueous solution. J. Phys. Chem. Lett. 2010, 1, 480–484. [Google Scholar] [CrossRef]
  70. Foote, C.S.; Dobrowolski, D.C. Oxygen Radicals in Chemistry and Biology; Bors, W., Saran, M., Tait, D., Eds.; De Gruyter: Berlin, Germany, 1984; pp. 465–470. [Google Scholar]
  71. Heihoff, K.; Redmond, R.W.; Braslavsky, S.E.; Rougee, M.; Salet, C.; Favre, A.; Bensasson, R.V. Quantum yields of triplet and O2(1Δg) formation of 4-thiouridine in water and acetonitrile. Phochem. Photobiol. 1990, 51, 635–641. [Google Scholar]
  72. Salet, C.; Bensasson, R.V.; Favre, A. Studies on the triplet excited state of 4-thiouridine. Photochem. Photobiol. 1983, 38, 521–525. [Google Scholar] [CrossRef]
  73. Pollum, M.; Ortiz-Rodríguez, L.A.; Jockusch, S.; Crespo-Hernández, C.E. The triplet state of 6-thio-2′-deoxyguanosine: Intrinsic properties and reactivity toward molecular oxygen. Photochem. Photobiol. 2016, 92, 286–292. [Google Scholar] [CrossRef] [PubMed]
  74. Ren, X.; Li, F.; Jeffs, G.; Zhang, X.; Xu, Y.-Z.; Karran, P. Guanine sulphinate is a major stable product of photochemical oxidation of DNA 6-thioguanine by UVA irradiation. Nucleic Acids Res. 2010, 38, 1832–1840. [Google Scholar] [CrossRef] [PubMed]
  75. Cooke, M.S.; Duarte, T.L.; Cooper, D.; Chen, J.; Nandagopal, S.; Evans, M.D. Combination of azathioprine and UVA irradiation is a major source of cellular 8-oxo-7,8-dihydro-2’-deoxyguanosine. DNA Repair 2008, 7, 1982–1989. [Google Scholar] [CrossRef] [PubMed]
  76. Montaner, B.; O’Donovan, P.; Reelfs, O.; Perrett, C.M.; Zhang, X.; Xu, Y.Z.; Ren, X.; Macpherson, P.; Frith, D.; Karran, P. Reactive oxygen-mediated damage to a human DNA replication and repair protein. EMBO Rep. 2007, 8, 1074–1079. [Google Scholar] [CrossRef] [PubMed]
  77. Zou, X.; Zhao, H.; Yu, Y.; Su, H. Formation of guanine-6-sulfonate from 6-thioguanine and singlet oxygen: A combined theoretical and experimental study. J. Am. Chem. Soc. 2013, 135, 4509–4515. [Google Scholar] [CrossRef] [PubMed]
  78. Zhang, X.; Jeffs, G.; Ren, X.; O’Donovan, P.; Montaner, B.; Perrett, C.M.; Karran, P.; Xu, Y.-Z. Novel DNA lesions generated by the interaction between therapeutic thiopurines and UVA light. DNA Repair 2007, 6, 344–354. [Google Scholar] [CrossRef] [PubMed]
  79. Brem, R.; Daehn, I.; Karran, P. Efficient DNA interstrand crosslinking by 6-thioguanine and UVA radiation. DNA Repair 2011, 10, 869–876. [Google Scholar] [CrossRef] [PubMed]
  80. Brem, R.; Li, F.; Montaner, B.; Reelfs, O.; Karran, P. DNA breakage and cell cycle checkpoint abrogation induced by a therapeutic thiopurine and UVA radiation. Oncogene 2010, 29, 3953–3963. [Google Scholar] [CrossRef] [PubMed]
  81. Wang, Z.; Rana, T.M. RNA−protein interactions in the Tat-trans-activation response element complex determined by site-specific photo-cross-linking. Biochemistry 1998, 37, 4235–4243. [Google Scholar] [CrossRef] [PubMed]
  82. Gueranger, Q.; Kia, A.; Frith, D.; Karran, P. Crosslinking of DNA repair and replication proteins to DNA in cells treated with 6-thioguanine and UVA. Nucleic Acids Res. 2011, 39, 5057–5066. [Google Scholar] [CrossRef] [PubMed]
  83. Martínez-Fernández, L.; Granucci, G.; Pollum, M.; Crespo-Hernández, C.E.; Persico, M.; Corral, I. Decoding the molecular basis for the population mechanism of the triplet phototoxic precursors in UVA light-activated pyrimidine anticancer drugs. Chem. Eur. J. 2017, 23, 2619–2627. [Google Scholar] [CrossRef] [PubMed]
  84. Igarashi-Yamamoto, N.; Tajiri, A.; Hatano, M.; Shibuya, S.; Ueda, T. Ultraviolet absorption, circular dichroism and magnetic circular dichroism studies of sulfur-containing nucleic acid bases and their nucleosides. Biochim. Biophys. Acta 1981, 656, 1–15. [Google Scholar] [CrossRef]
  85. Mai, S.; Ashwood, B.; Marquetand, P.; Crespo-Hernández, C.E.; González, L. Solvatochromic effects on the absorption spectrum of 2-thiocytosine. J. Phys. Chem. B 2017, 121, 5187–5196. [Google Scholar] [CrossRef] [PubMed]
  86. Vendrell-Criado, V.; Sáez, J.A.; Lhiaubet-Vallet, V.; Cuquerella, M.C.; Miranda, M.A. Photophysical properties of 5-substituted 2-thiopyrimidines. Photochem. Photobiol. Sci. 2013, 12, 1460–1465. [Google Scholar] [CrossRef] [PubMed]
  87. Martínez-Fernández, L.; Pepino, A.J.; Segarra-Martí, J.; Banyasz, A.; Garavelli, M.; Improta, R. Computing the absorption and emission spectra of 5-methylcytidine in different solvents: A test-case for different solvation models. J. Chem. Theor. Comput. 2016, 12, 4430–4439. [Google Scholar] [CrossRef] [PubMed]
  88. Improta, R.; Barone, V. Excited states behavior of nucleobases in solution: Insights from computational studies. In Photoinduced Phenomena in Nucleic Acids I: Nucleobases in the Gas Phase and in Solvents; Barbatti, M., Borin, A.C., Ullrich, S., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 329–357. [Google Scholar]
  89. Mai, S.; Marquetand, P.; González, L. A static picture of the relaxation and intersystem crossing mechanisms of photoexcited 2-thiouracil. J. Phys. Chem. A 2015, 119, 9524–9533. [Google Scholar] [CrossRef] [PubMed]
  90. Cui, G.; Fang, W.-H. State-specific heavy-atom effect on intersystem crossing processes in 2-thiothymine: A potential photodynamic therapy photosensitizer. J. Chem. Phys. 2013, 138, 044315–044323. [Google Scholar] [CrossRef] [PubMed]
  91. Gobbo, J.P.; Borin, A.C. 2-thiouracil deactivation pathways and triplet states population. Comput. Theor. Chem. 2014, 1040–1041, 195–201. [Google Scholar] [CrossRef]
  92. Jiang, J.; Zhang, T.-S.; Xue, J.-D.; Zheng, X.; Cui, G.; Fang, W.-H. Short-time dynamics of 2-thiouracil in the light absorbing S2(ππ) state. J. Chem. Phys. 2015, 143, 175103–175111. [Google Scholar] [CrossRef] [PubMed]
  93. Bai, S.; Barbatti, M. Why replacing different oxygens of thymine with sulfur causes distinct absorption and intersystem crossing. J. Phys. Chem. A 2016, 120, 6342–6350. [Google Scholar] [CrossRef] [PubMed]
  94. Taherian, M.R.; Maki, A.H. Optically detected magnetic resonance study of the phosphorescent states of thiouracils. Chem. Phys. 1981, 55, 85–96. [Google Scholar] [CrossRef]
  95. Cui, G.; Thiel, W. Intersystem crossing enables 4-thiothymidine to act as a photosensitizer in photodynamic therapy: An ab initio QM/MM study. J. Phys. Chem. Lett. 2014, 5, 2682–2687. [Google Scholar] [CrossRef] [PubMed]
  96. Xie, B.-B.; Wang, Q.; Guo, W.-W.; Cui, G. The excited-state decay mechanism of 2,4-dithiothymine in the gas phase, microsolvated surroundings, and aqueous solution. Phys. Chem. Chem. Phys. 2017, 19, 7689–7698. [Google Scholar] [CrossRef] [PubMed]
  97. El-Sayed, M.A. Spin—Orbit coupling and the radiationless processes in nitrogen heterocyclics. J. Chem. Phys. 1963, 38, 2834–2838. [Google Scholar] [CrossRef]
  98. Mai, S.; Marquetand, P.; González, L. Intersystem crossing pathways in the noncanonical nucleobase 2-thiouracil: A time-dependent picture. J. Phys. Chem. Lett. 2016, 7, 1978–1983. [Google Scholar] [CrossRef] [PubMed]
  99. Zou, X.; Dai, X.; Liu, K.; Zhao, H.; Song, D.; Su, H. Photophysical and photochemical properties of 4-thiouracil: Time-resolved IR spectroscopy and DFT studies. J. Phys. Chem. B 2014, 118, 5864–5872. [Google Scholar] [CrossRef] [PubMed]
  100. Wenska, G.; Taras-Goślińska, K.; Skalski, B.; Maciejewski, A.; Burdziński, G.; Karolczak, J. Putative phototautomerization of 4-thiouridine in the S2 excited state revealed by fluorescence study using picosecond laser spectroscopy. J. Photochem. Photobiol. A 2006, 181, 12–18. [Google Scholar] [CrossRef]
  101. Taras-Goślińska, K.; Wenska, G.; Skalski, B.; Maciejewski, A.; Burdziński, G.; Karolczak, J. Spectral and photophysical properties of the lowest excited triplet state of 4-thiouridine and its 5-halogeno derivatives. J. Photochem. Photobiol. A 2004, 168, 227–233. [Google Scholar] [CrossRef]
  102. Taras-Goślińska, K.; Wenska, G.; Skalski, B.; Maciejewski, A.; Burdziński, G.; Karolczak, J. Intra- and intermolecular electronic relaxation of the second excited singlet and the lowest excited triplet states of 1,3-dimethyl-4-thiouracil in solution. Photochem. Photobiol. 2002, 75, 448–456. [Google Scholar] [CrossRef]
  103. Wenska, G.; Taras-Goślińska, K.; Łukaszewicz, A.; Burdzinski, G.; Koput, J.; Maciejewski, A. Mechanism and dynamics of intramolecular triplet state decay of 1-propyl-4-thiouracil and its α-methyl-substituted derivatives studied in perfluoro-1,3-dimethylcyclohexane. Photochem. Photobiol. Sci. 2011, 10, 1294–1302. [Google Scholar] [CrossRef] [PubMed]
  104. Shukla, M.K.; Leszczynski, J. Multiconfigurational self-consistent field study of the excited state properties of 4-thiouracil in the gas phase. J. Phys. Chem. A 2004, 108, 7241–7246. [Google Scholar] [CrossRef]
  105. Wenska, G.; Koput, J.; Burdziński, G.; Taras-Goślińska, K.; Maciejewski, A. Photophysical and photochemical properties of the T1 excited state of thioinosine. J. Photochem. Photobiol. A 2009, 206, 93–101. [Google Scholar] [CrossRef]
  106. Pepino, A.J.; Segarra-Martí, J.; Nenov, A.; Improta, R.; Garavelli, M. Resolving ultrafast photoinduced deactivations in water-solvated pyrimidine nucleosides. J. Phys. Chem. Lett. 2017, 8, 1777–1783. [Google Scholar] [CrossRef] [PubMed]
  107. Wenska, G.; Filipiak, P.; Taras-Goślińska, K.; Sobierajska, A.; Gdaniec, Z. Orientation-dependent quenching of the triplet excited 6-thiopurine by nucleobases. J. Photochem. Photobiol. A 2011, 217, 55–61. [Google Scholar] [CrossRef]
Figure 1. HOMO (H) and LUMO (L) molecular orbitals involved in the ππ* lowest spectroscopic state in 2t-Cyt, 2t-Thy, 6t-Gua, 4t-Thy and 2,4-dtThy.
Figure 1. HOMO (H) and LUMO (L) molecular orbitals involved in the ππ* lowest spectroscopic state in 2t-Cyt, 2t-Thy, 6t-Gua, 4t-Thy and 2,4-dtThy.
Molecules 22 00998 g001
Figure 2. Potential energy profiles for 6t-Gua (Left) and 2t-Cyt (Right) inferred from quantum mechanical studies. Energies in eV relative to the ground state. For 6t-Gua, underlined energies were taken from Ref. [58]. The rest were extracted from Ref. [55]. For 2t-Cyt, energies were taken from Ref. [34]. The double well profile of T1 is based on Ref. [59].
Figure 2. Potential energy profiles for 6t-Gua (Left) and 2t-Cyt (Right) inferred from quantum mechanical studies. Energies in eV relative to the ground state. For 6t-Gua, underlined energies were taken from Ref. [58]. The rest were extracted from Ref. [55]. For 2t-Cyt, energies were taken from Ref. [34]. The double well profile of T1 is based on Ref. [59].
Molecules 22 00998 g002
Scheme 1. Structures of the thiosubstituted analogs of the nucleobases considered in this work. Atom numbering is given for the heterocyclic rings of the nucleobases.
Scheme 1. Structures of the thiosubstituted analogs of the nucleobases considered in this work. Atom numbering is given for the heterocyclic rings of the nucleobases.
Molecules 22 00998 sch001
Figure 3. Proposed kinetic models for nucleobase thioderivatives inferred from the literature.
Figure 3. Proposed kinetic models for nucleobase thioderivatives inferred from the literature.
Molecules 22 00998 g003
Scheme 2. Oxidation products of 6-thioguanine and guanine (Top); and their generation mechanism according to Ref. [77] (Bottom).
Scheme 2. Oxidation products of 6-thioguanine and guanine (Top); and their generation mechanism according to Ref. [77] (Bottom).
Molecules 22 00998 sch002
Figure 4. Potential Energy Profiles for 2t-Ura (Top left) and 4t-Thy (Top right) and 2,4-dtThy (Bottom) inferred from quantum mechanical studies. Energies in eV relative to the ground state. For 2t-Ura, underlined energies were extracted from Ref. [90], italic from Ref. [91] and the rest from Ref. [89]. For 4t-Thy, underlined energies were taken from Ref. [95] and the rest were extracted from Ref. [83]. The energies signaled with an asterisk denote geometries optimized including seven water molecules. For 2,4-dtThy, energies were taken from Ref. [96].
Figure 4. Potential Energy Profiles for 2t-Ura (Top left) and 4t-Thy (Top right) and 2,4-dtThy (Bottom) inferred from quantum mechanical studies. Energies in eV relative to the ground state. For 2t-Ura, underlined energies were extracted from Ref. [90], italic from Ref. [91] and the rest from Ref. [89]. For 4t-Thy, underlined energies were taken from Ref. [95] and the rest were extracted from Ref. [83]. The energies signaled with an asterisk denote geometries optimized including seven water molecules. For 2,4-dtThy, energies were taken from Ref. [96].
Molecules 22 00998 g004
Figure 5. Experimental absorption spectra recorded in PBS, pH = 7.4 for Thy, Guo and their thioderivatives. The data for Thy and it derivatives (Top) was taken from Ref. [64]. The data for Guo and 6t-Guo (Bottom) was taken from Ref. [73] and the spectrum for 6t-Gua was extracted from Ref. [52].
Figure 5. Experimental absorption spectra recorded in PBS, pH = 7.4 for Thy, Guo and their thioderivatives. The data for Thy and it derivatives (Top) was taken from Ref. [64]. The data for Guo and 6t-Guo (Bottom) was taken from Ref. [73] and the spectrum for 6t-Gua was extracted from Ref. [52].
Molecules 22 00998 g005
Figure 6. General scheme summarizing the preferred route towards the population of triplet states upon UVA excitation in thiated nucleobases.
Figure 6. General scheme summarizing the preferred route towards the population of triplet states upon UVA excitation in thiated nucleobases.
Molecules 22 00998 g006
Table 1. Steady-state absorption and emission properties of 6-thioguanine nucleobase and nucleosides at room temperature unless otherwise indicated. F: Fluorescence; P: Phosphorescence; N.D.: Not detected.
Table 1. Steady-state absorption and emission properties of 6-thioguanine nucleobase and nucleosides at room temperature unless otherwise indicated. F: Fluorescence; P: Phosphorescence; N.D.: Not detected.
λmax (nm) εmax (M−1 cm−1 × 104)λemission (nm)Solvent
FP
6t-Gua3371.80NaOH, pH = 10 [56]
3402.10Tris buffer, pH = 7.4 [56]
340, 255, 220, 205400H2O, pH = 2.3 [51]
340, 255, 220, 205400, 500 1H2O, pH = 7.8 [51]
320, 265, 245, 210400H2O, pH = 12.8 [51]
341, 254, 220 1, 2041.80, 0.60, 1.40, 1.90PBS, pH = 7.4 [52]
347 22.08478pH = 1.0, 77 K [50]
343 22.22459pH = 5.0, 77 K [50]
323 21.85442pH = 10.0, 77 K [50]
321 22.124421 N NaOH, 77 K [50]
468, 437 1, 494 1430–550Tris buffer, pH = 7.4, 77 K [52]
478NaOH, pH = 12.1 [57]
6t-Guo3201.80NaOH, pH = 10 [56]
3422.00Tris buffer, pH = 7.4 [56]
3422.30N.D.N.D.PBS, pH = 7 [54]
3462.20N.D.N.D.ACN [54]
341, 258, 227 1, 2082.30, 0.80, 1.20, 2.10PBS, pH = 7.4 [52]
348 2481pH = 1.0 ,77 K [50]
344 2459pH = 5.0, 77 K [50]
316 2431pH = 12.0, 77 K [50]
1 Band shoulder. 2 Absorption maxima at room temperature.
Table 2. Computed absorption and emission properties of 6-thioguanine nucleobase and nucleosides. Only information relative to the lower lying spectroscopic state and preceding electronic states is included. In parenthesis the oscillator strengths are shown.
Table 2. Computed absorption and emission properties of 6-thioguanine nucleobase and nucleosides. Only information relative to the lower lying spectroscopic state and preceding electronic states is included. In parenthesis the oscillator strengths are shown.
λmax (nm)SolventExcitationλemission (nm)Method
FP
6t-Gua [55]369 (0.000)S1 1nπ*395MS-CASPT2//
CASSCF(14,12)/ANO-L
306 (0.535)S2 1ππ*
400T1 3ππ*404
375T2 3nπ*
6t-Gua [53]397 (0.000)S1 TD-B3LYP/6-311++G(d,p)
311 (0.203)S2
6t-Gua [52]381 (0.000)S1 1nπ*TD-PBE0/6-311++G(d,p)/
IEFPCM
346 (0.000)H2OS1 1nπ*
299 (0.260)S2 1ππ*
310 (0.400)H2OS2 1ππ*
6t-Gua [50]337 (0.470)S2 1ππ*CNDO/S
6t-Gua [51]400 (0.000)S1 1nπ*CNDO/s-CI
340 (0.290)S2 1ππ*
6t-Guo [54]384 (0.000)S1 1nπ*TD-PBE0/6-311++G(d,p)/
IEFPCM
341 (0.000)H2OS1 1nπ*
342 (0.000)ACNS1 1nπ*
305 (0.280)S2 1ππ*
313 (0.450)H2OS2 1ππ*
313 (0.450)ACNS2 1ππ*
6t-Guo [38]359 (0.000)S1 1nπ*376TD-B3LYP/6-31+G*
319 (0.434)S2 1ππ*
457T1 3ππ*496
357 (0.000)S1 1nπ*TD-M06/6-31+G*
314 (0.455)S2 1ππ*
Table 3. Photophysical parameters (lifetimes in ps unless otherwise indicated and triplet and singlet oxygen quantum yields) for thiated nucleobase analogs.
Table 3. Photophysical parameters (lifetimes in ps unless otherwise indicated and triplet and singlet oxygen quantum yields) for thiated nucleobase analogs.
τ ISC (ps)Φ ISCLifetimesΦΔ
2t-Ura0.35 ± 0.06 8 [60]0.75 ± 0.20 8 [60] τ1 < 0.20, τ2 = 0.35 ± 0.06 8, 0.34 ± 0.09 1 [60]
0.34 ± 0.09 1 [60]
0.36 ± 0.03 8 [61]
2t-Thy0.62 ± 0.07 8 [60]1.00 ± 0.05 3 [62]τ1 < 0.20, τ2 = 0.62 ± 0.07 8, 0.32 ± 0.09 1 [60], τ3 = 2.7 ± 0.5 (μs) [62] 0.36 ± 0.02 3 [62]
0.32 ± 0.09 1 [60]0.9 ± 0.1 2 [63]τ1 < 0.10, τ2 = 0.78, τ3 = 203 [39]
2t-Thd0.41 ± 0.06 8 [64]0.9 ± 0.1 2 [63]
2t-Cyt0.21 ± 0.05 8 [34]>0.9 8 [34]τ1 = 0.21 ± 0.05, τ2 = 0.48 ± 0.06 [34]
dm-4tUra1.0 ± 0.1 9 [65]
4t-Ura0.24 ± 0.02 8 [61]0.9 ± 0.15 8 [61]0.49 ± 0.02 3 [61]
4t-Thd0.24 ± 0.02 8 [64] 1.0 ± 0.1 2 [66]0.50 ± 0.10 3 [66]
0.24 ± 0.02 8 [67]0.85 ± 0.15 8 [67]τ1 = 0.24 ± 0.02, τ2 = 84 ± 2 [67]0.42 ± 0.02 3 [64]
0.54 ± 0.01 1 [68]τ1 = 0.54 ± 0.01, τ2 = 45 ± 5 [68], τ1 = 0.54, τ2 = 1.8 [57]
~10 1 [69]
4t-Urd 0.7 [70]
0.9 ± 0.1 9 [65]0.18 ± 0.04 10 [4,71]
0.02 2 [72]0.50 ±0.02 3 [4,71]
2,4-dtUra0.22 ± 0.04 8 [61]0.9 ± 0.15 8 [61]0.49 ± 0.02 3 [61]
2,4-dtThy0.18 ± 0.02 8 [64] >0.9 8 [64]0.46 ± 0.02 3 [64]
6t-Gua0.58 ± 0.08 6 [56]
0.56 ± 0.18 5 [56]
0.56 ± 0.06 8 [52]≥0.6± 0.2 8 [52]τ1 = 0.56 ± 0.06, τ2 = 26 ± 3, τ3 = 830 ± 70 (ns), τ3 = 1420 ± 180 (ns) (N2) [52] 0.21 ± 0.02 4 [52]
0.23 ± 0.02 5 [52]
6t-Guo0.31 ± 0.05 7 [54]0.8 ± 0.2 7 [54]τ1 = 0.31 ± 0.05, τ2 = 80 ± 15, τ3 = 460 ± 15 (ns), τ3 = 720 ± 10 (ns) (N2) [54]0.49 ± 0.09 6 [56]
0.24 ± 0.02 5 [73]
0.14 ± 0.02 4 [73]
0.29 ± 0.02 3 [73]
0.36 ± 0.04 1 [54]τ1 = 0.36 ± 0.04, τ2 = 32 ± 5 [54] 0.55 ± 0.08 5 [56]
1 ACN; 2 ACN, Ar; 3 ACN, O2; 4 D2O, Tris buffer, pH = 7.4; 5 D2O, Tris buffer, pH = 7.4, O2; 6 NaOH/D2O, pH = 10; 7 PBS, pH = 7.0; 8 PBS, pH = 7.4; 9 H2O, Ar, 10 H2O, O2.
Table 4. Steady-state absorption and emission properties of 2-thiopyrimidine derivatives at room temperature unless otherwise indicated. F: Fluorescence; P: Phosphorescence; N.D.: Not detected.
Table 4. Steady-state absorption and emission properties of 2-thiopyrimidine derivatives at room temperature unless otherwise indicated. F: Fluorescence; P: Phosphorescence; N.D.: Not detected.
λmax (nm)εmax (M−1 cm−1 × 104)λemission (nm)Φ (×10−4)Solvent
FPFP
2t-Ura2741.45H2O, pH = 7 [84]
2681.19ACN [84]
274, 291 1_MeOH [84]
271, 290 1N.D.405, 427, 4557.0EtOH, 77 K [86]
265PBS, pH = 7.4 [61]
2t-Urd2751.67H2O, pH = 7 [84]
2681.04ACN [84]
2t-Thy275, 290 11.29, 1.33N.D.N.D.ACN [62]
454THF, 77 K [62]
2901.24N.D.4805.0, 9.8 2ACN, Ar [63]
448DCM:MeOH, 77 K [63]
275, 290 1N.D.425, 451, 4837.0EtOH, 77 K [86]
275PBS, pH = 7.4 [64]
2t-Thd2851.41N.D.N.D.<1 3ACN, Ar [63]
448DCM:MeOH, 77 K [63]
277PBS, pH = 7.4 [64]
2t-Cyt270, 242, 220 11.82H2O, pH = 7 [84]
269, 241, 219 1PBS, pH = 7.4 [34]
283, 2331.67, 1.19ACN [85]
2861.69EtOAc [85]
2851.80DMSO [85]
280, 242, 210 11.79, 1.42, 0.65EtOH [85]
278, 242, 212 11.78, 1.48, 0.71MeOH [85]
269, 242, 218 11.80, 1.80, 0.90PBS, pH = 7.4 [85]
5m-2tCyt2662.10H2O, pH = 7 [84]
2t-Cyd274, 270 12.30, 1.40H2O, pH = 7 [84]
1 Band shoulder. 2 Calculated to infinite dilution. 3 Estimated due to the sensitivity limit of the instrument.
Table 5. Computed absorption and emission properties of 2-thiopyrimidine derivatives. Only information concerning the lowest lying spectroscopic state and preceding states is included. In parenthesis the oscillator strengths are shown. F: Fluorescence; P: Phosphorescence.
Table 5. Computed absorption and emission properties of 2-thiopyrimidine derivatives. Only information concerning the lowest lying spectroscopic state and preceding states is included. In parenthesis the oscillator strengths are shown. F: Fluorescence; P: Phosphorescence.
λmax (nm)SolventExcitationλemission (nm)Method
FP
2t-Ura 1 [89]329 (0.000)S1 1nπ*358MS(6/4)-CASPT2/ANO-RCC-VQZQ//
SA(6/4)-CASSCF(16,12)
292 (0.350)S2 1ππ*
353T1 3ππ*386
322 (0.000)S1 1nπ*MRCIS(6,5)/cc-pVDZ//
SA(4/2)-CASSCF(12,9)
215 (0.110)S2 1ππ*
2t-Ura 1 [90]323S1 1nπ*325CASPT2/cc-pVDZ//
SA(4/3)-CASSCF(16,11)
278S2 1ππ*
T1 3ππ*358
347S1 1nπ*TD-B3LYP/aug-cc-pVDZ
304S2 1ππ*
2t-Ura [91]340 (0.000)S1 1nπ*369SS-CASPT2/ANO-L-vDZP//
SA(6/6)-CASSCF(14,10)
303 (0.236)S2 1ππ*349
395T1 3ππ*435
2t-Ura [31]346 (0.000)S1 1nπ*TD-B3LYP/6-311++G(d,p)
304 (0.051)S2 1ππ*
317 (0.000)H2OS1 1nπ*
286 (0.145)H2OS2 1ππ*
318 (0.000)ACNS1 1nπ*
287 (0.144)ACNS2 1ππ*
2t-Ura [92]320 (0.000)ACNS1 1nπ*TD-B3LYP/6-31+G(d)/PCM
283 (0.108)ACNS2 1ππ*
2t-Ura [86]332 (0.000)EtOHS1 1nπ*TD-B3LYP/aug-cc-pVDZ/PCM
298 (0.110)EtOHS2 1ππ*
2t-Thy [37]325H2OS1 1nπ*TD-M06/6-31+G*/PCM
288H2OS2 1ππ*
2t-Thy 1,2 [93]286 (0.578)1ππ*MS-CASPT2(14,10)/ANO-RCC-VTZP
295 (0.243)1ππ*DFT/MRCI
2t-Thy [62]315 (0.000)ACNS1 1nπ*TD-B3LYP/6-31+G(d,p)/PCM
295 (0.145)ACNS2 1ππ*
2t-Thy [86]323 (0.000) EtOHS1 1nπ*TD-B3LYP/aug-cc-pVDZ/PCM
301 (0.125)EtOHS2 1ππ*
2t-Thy [53]339 (0.000)S1TD-B3LYP/6-311++G(d,p)
308 (0.054)S2
2t-Cyt 1 [34]340 (0.001)S1 1nπ*420MS-CASPT2/CASSCF(14,10)/
ANO-RCC-VQZP//RI-MP2/cc-pVQZ
332 (0.097)S2 1ππ*410
368T1 3ππ*435
2t-Cyt [85]358 (0.000)S1 1nπ*MS-CASPT2(14,10)/cc-pVDZ/PCM//BP86/aug-cc-pVDZ/COSMO
315 (0.000)EtOAcS1 1nπ*
309 (0.000)ACNS1 1nπ*
309 (0.000)DMSOS1 1nπ*
286 (0.000)EtOHS1 1nπ*
284 (0.000)MeOHS1 1nπ*
276 (0.020)H2OS1 1nπ*
357 (0.020)S2 1ππ*
315 (0.010)EtOAcS2 1ππ*
311 (0.020)ACNS2 1ππ*
307 (0.020)DMSOS2 1ππ*
294 (0.030)EtOHS2 1ππ*
292 (0.030)MeOHS2 1ππ*
293 (0.010)H2OS2 1ππ*
2t-Cyt [53]403 (0.000)S1TD-B3LYP/6-311++G(d,p)
371 (0.010)S2
1 This work reports additional values obtained at other levels of theory. 2 This work reports values computed considering other solvents.
Table 6. Steady-state absorption and emission properties of 4-thiopyrimidine derivatives at room temperature unless otherwise indicated. F: Fluorescence; P: Phosphorescence; N.D.: Not detected, PFDMCH: perfluoro-1,3-dimethylcyclohexane.
Table 6. Steady-state absorption and emission properties of 4-thiopyrimidine derivatives at room temperature unless otherwise indicated. F: Fluorescence; P: Phosphorescence; N.D.: Not detected, PFDMCH: perfluoro-1,3-dimethylcyclohexane.
λmax (nm)εmax (M−1 cm−1 × 104)λemission (nm)Φ (×10−4)Solvent
FPFP
4t-Ura3271.28ACN [84]
3281.67H2O, pH = 7 [84]
3271.94ACN [99]
330550 3.0 1H2O [17]
dm-4tUra3281.98415 528 0.6 114.0 1CCl4, He [102]
317416 535 1.5 125.0 1PFDMCH, He [102]
3291.88420 5350.4H2O [102]
4t-Urd 3312.30H2O, pH = 7 [84]
3311.19ACN [84]
Ac-4tUrd 3282.06420 550 1.02.0ACN [100]
3302.06420540123 1CCl4, He [101]
25 1ACN, Ar [105]
4t-Thd 337 410550ACN [68]
335 400542 PBS, pH = 7.4 [67]
3371.94 N.D.N.D.ACN [66]
494, 512, 560 2EtOH, 77 K [66]
1 Calculated to infinite dilution. 2 Band shoulder.
Table 7. Computed absorption and emission properties of 4-thiopyrimidine derivatives. In parenthesis the oscillator strengths are shown. F: Fluorescence; P: Phosphorescence.
Table 7. Computed absorption and emission properties of 4-thiopyrimidine derivatives. In parenthesis the oscillator strengths are shown. F: Fluorescence; P: Phosphorescence.
λmax (nm)SolventExcitationλemission (nm)Method
FP
4t-Ura [31]446 (0.000)S1 1nπ*TD-B3LYP/6-311++G(d,p)
295 (0.298)S2 1ππ*
387 (0.000)H2OS1 1nπ*
298 (0.437)H2OS2 1ππ*
389 (0.000)ACNS1 1nπ*
299(0.440)ACNS2 1ππ*
4t-Ura [104]401S1 1nπ*CASSCF
259S2 1ππ*
440S1 1nπ*MCQDPT2
318S2 1ππ*
4t-Ura [99] 444 (0.000) S1 1nπ*TD-B3LYP/6-311++G(d,p)
294 (0.299)S2 1ππ*
dm-4tUra [103]408S1 1nπ*EOM-CC2/aug-cc-pVDZ
296S2 1ππ*
4t-Thy 1 [93]310 (0.662)1ππ*MS-CASPT2 (14,10)/ANO-RCC/VTZP
312 (0.427)1ππ*DFT/MRCI
4t-Thy [53]438 (0.000)S1TD-B3LYP/6-311++G(d,p)
298 (0.282)S2
4t-Thy [37]403H2OS1 1nπ*TD-M06/6-31+G*/PCM
312H2OS2 1ππ*
4t-Thy [83]419 (0.000)H2O (PCM + 7 H2O )S1 1nπ*MS-CASPT2 (12,9)/ANO-L
326 (0.546)H2O (PCM + 7 H2O )S2 1ππ*430
448T1 3ππ*570
4t-Thd 1 [67]385 (0.000)H2O (PCM + 2 H2O )S1 1nπ*TD-PBE0/IEFPCM/6-311++G(d,p)
310 (0.509)H2O (PCM + 2 H2O )S2 1ππ*
478H2O (PCM + 2 H2O )T1 3ππ*557 2
4t-Thd [66]292 (0.379)S2 1ππ*TD-B3LYP/6-31G(d,p)
4t-Thd [69]~435S1 1nπ*TD-B3LYP/6-31G(d,p)
~294S2 1ππ*
4t-Thd [95]413H2O (MM)S1 1nπ*QM(CASPT2//CASSCF(10,8))/6-31G*/MM
302H2O (MM)S2 1ππ*
1 This work reports additional values obtained at other levels of theory and computed considering other solvents. 2 Calculated at the same level of theory without the explicit water molecules for 4t-Thd.
Table 8. Steady-state absorption and emission properties of 2,4-dithiopyrimidine derivatives at room temperature. F: Fluorescence; P: Phosphorescence.
Table 8. Steady-state absorption and emission properties of 2,4-dithiopyrimidine derivatives at room temperature. F: Fluorescence; P: Phosphorescence.
λmax (nm)εmax (M−1 cm−1 × 104)λemission (nm)Solvent
FP
2,4-dtUra 3510.46H2O, pH = 7 [84]
3510.10ACN [84]
2,4-dtUrd 3451.14 H2O, pH = 7 [84]
3511.00 ACN [84]
2,4-dtThy3630.97PBS [64]
Table 9. Computed absorption and emission properties of 2,4-dithiopyrimidine derivatives. In parenthesis the oscillator strengths are shown. F: Fluorescence; P: Phosphorescence.
Table 9. Computed absorption and emission properties of 2,4-dithiopyrimidine derivatives. In parenthesis the oscillator strengths are shown. F: Fluorescence; P: Phosphorescence.
λmax (nm)SolventExcitationλemission (nm)Method
FP
2,4-dtUra [31]449 (0.000)S1 1nπ*TD-B3LYP/6-311++G(d,p)
349 (0.014)S2 1ππ*
396 (0.000)H2OS1 1nπ*
335 (0.063)H2OS2 1ππ*
399 (0.000)ACNS1 1nπ*
336 (0.063)ACNS2 1ππ*
2,4-dtThy 1 [93]333 (0.216)1ππ*MS-CASPT2(14,10)/ANO-RCC-VTZP
343 (0.078)1ππ*DFT/MRCI
2,4-dtThy 1 [96]400H2OS1 1nπ*446QM(MS-CASPT2)/MM
334H2OS2 1ππ*
445H2OT1 1ππ*474
2,4-dtThy [53]443 (0.000)S1TD-B3LYP/6-311++G(d,p)
361 (0.014)S2
2,4-dtThy [37]409H2O S1TD-M06/6-31+G*/PCM
340H2OS2
1 This work reports additional values obtained at other levels of theory and computed considering other solvents.

Share and Cite

MDPI and ACS Style

Arslancan, S.; Martínez-Fernández, L.; Corral, I. Photophysics and Photochemistry of Canonical Nucleobases’ Thioanalogs: From Quantum Mechanical Studies to Time Resolved Experiments. Molecules 2017, 22, 998. https://doi.org/10.3390/molecules22060998

AMA Style

Arslancan S, Martínez-Fernández L, Corral I. Photophysics and Photochemistry of Canonical Nucleobases’ Thioanalogs: From Quantum Mechanical Studies to Time Resolved Experiments. Molecules. 2017; 22(6):998. https://doi.org/10.3390/molecules22060998

Chicago/Turabian Style

Arslancan, Serra, Lara Martínez-Fernández, and Inés Corral. 2017. "Photophysics and Photochemistry of Canonical Nucleobases’ Thioanalogs: From Quantum Mechanical Studies to Time Resolved Experiments" Molecules 22, no. 6: 998. https://doi.org/10.3390/molecules22060998

APA Style

Arslancan, S., Martínez-Fernández, L., & Corral, I. (2017). Photophysics and Photochemistry of Canonical Nucleobases’ Thioanalogs: From Quantum Mechanical Studies to Time Resolved Experiments. Molecules, 22(6), 998. https://doi.org/10.3390/molecules22060998

Article Metrics

Back to TopTop