Sulphur-atom positional engineering in perylenimide: structure-property relationships and H-aggregation directed type-I photodynamic therapy

An innovative design strategy of placing sulfur (S)-atoms within the pendant functional groups and at carbonyl positions in conventional perylenimide (PNI-O) has been demonstrated to investigate the condensed state structure-property relationship and potential photodynamic therapy (PDT) application. Incorporation of simply S-atoms at the peri-functionalized perylenimide (RPNI-O) leads to an aggregation-induced enhanced emission luminogen (AIEEgen), 2-hexyl-8-(thianthren-1-yl)-1H-benzo[5,10]anthra[2,1,9-def]isoquinoline-1,3(2H)-dione (API), which achieves a remarkable photoluminescence quantum yield (Φ PL) of 0.85 in aqueous environments and established novel AIE mechanisms. Additionally, substitution of the S-atom at the carbonyl position in RPNI-O leads to thioperylenimides (RPNI-S): 2-hexyl-8-phenyl-1H-benzo[5,10]anthra[2,1,9-def]isoquinoline-1,3(2H)-dithione (PPIS), 8-([2,2'-bithiophen]-5-yl)-2-hexyl-1H-benzo[5,10]anthra[2,1,9-def]isoquinoline-1,3(2H)-dithione (THPIS), and 2-hexyl-8-(thianthren-1-yl)-1H-benzo[5,10]anthra[2,1,9-def]isoquinoline-1,3(2H)-dithion (APIS), with distinct photophysical properties (enlarged spin-orbit coupling (SOC) and Φ PL ≈ 0.00), and developed diverse potent photosensitizers (PSs). The present work provides a novel SOC enhancement mechanism via pronounced H-aggregation. Surprisingly, the lowest singlet oxygen quantum yield (Φ Δ) and theoretical calculation suggest the specific type-I PDT for RPNI-S. Interestingly, RPNI-S efficiently produces superoxide (O2˙-) due to its remarkably lower Gibbs free energy (ΔG) values (THPIS: -40.83 kcal mol-1). The non-toxic and heavy-atom free very specific thio-based PPIS and THPIS PSs showed selective and efficient PDT under normoxia, as a rare example.

Unexpected longer T1 lifetime of 6-sulfur guanine than 6-selenium guanine: the solvent effect of hydrogen bonds to brake the triplet decay

The decay of the T1 state to the ground state is an essential property of photosensitizers because it decides the lifetime of excited states and, thus, the time window for sensitization. The sulfur/selenium substitution of carbonyl groups can red-shift absorption spectra and enhance the triplet yield because of the large spin-orbit coupling, modifying nucleobases to potential photosensitizers for various applications. However, replacing sulfur with selenium will also cause a much shorter T1 lifetime. Experimental studies found that the triplet decay rate of 6-seleno guanine (6SeGua) is 835 times faster than that of 6-thio guanine (6tGua) in aqueous solution. In this work, we reveal the mechanism of the T1 decay difference between 6SeGua and 6tGua by computing the activation energy and spin-orbit coupling for rate calculation. The solvent effect of water is treated with explicit microsolvation and implicit solvent models. We find that the hydrogen bond between the sulfur atom of 6tGua and the water molecule can brake the triplet decay, which is weaker in 6SeGua. This difference is crucial to explain the relatively long T1 lifetime of 6tGua in an aqueous solution. This insight emphasizes the role of solvents in modulating the excited state dynamics and the efficiency of photosensitizers, particularly in aqueous environments.

Disruptive Model That Explains for the Long-Lived Triplet States Observed for 2-Thiocytosine upon UVA Radiation

The possible role of radical species in the formation of the long-lived triplet states observed for 2-thiocytosine upon UV irradiation was theoretically investigated. It is predicted that the radical fragments arising from the homolytic rupture of the NH group of the thiobase can be yielded upon ultraviolet-A radiation. Recombination of the radicals through the most favorable singlet channel yields the lowest-lying tautomer of the 2-thiocytosine (the amino-thiol form) through a barrierless pathway. The rebounding of the radical fragments along the triplet channels that emerge from the attack of the hydrogen to the nitrogen atoms next to the C-S bond leads to stable structures for the amino-thion-N1H and amino-thion-N3H tautomers. These results allow for the rationalization of the near-unity triplet yields observed when this pure light-atom organic molecule is exposed to UV irradiation, without invoking intersystem crossings between the electronic states of different spin-multiplicities. A similar study for cytosine showed that the energy required to induce the homolytic breaking of the N-H bond of the nucleobase is not attainable under UVA radiation. This result is consistent with the experimental fact that no triplet states are observed when this molecule is exposed to that light.

Dual Functionality of 6-Methylthioguanine: Synergistic Effects Enhancing the Photolability of DNA Nucleobases

A small chemical modification of the nucleobase structure can significantly enhance the photoactivity of DNA, which may incur DNA damage, thus holding promising applications in photochemotherapy treatment of cancers or pathogens. However, single substitution confers only limited phototoxicity to DNA. Herein, we combine femtosecond and nanosecond time-resolved spectroscopy with high-level ab initio calculations to disentangle the excited-state dynamics of 6-methylthioguanine (me6-TG) under variable wavelength UVA excitation (310-330 nm). We find that double substitution of nucleobases (thionation and methylation) boosts the photoactivity by introducing more reactive channels. Intriguingly, 1nNπ*, rather than 1nSπ*, acts as the doorway state engendering the formation of the long-lived reactive triplet state in me6-TG. The 1nNπ* induces a low spin-orbit coupling of 8.3 cm-1, which increases the intersystem crossing (ISC) time (2.91 ± 0.14 ns). Despite the slowed ISC, the triplet quantum yield (ΦT) still accounts for a large fraction (0.6 ± 0.1), consistent with the potential energy surface that favors excited-state bifurcation to 1nNπ*min (3.36 ± 0.15 ps) rather than 1ππ*min (5.05 ± 0.26 ps), such that the subsequent ISC to triplet via 1nNπ*min constitutes the main relaxation pathway in me6-TG. Although this ΦT is inferior to its single-substituted predecessor 6-thioguanine (6-TG, 0.8 ± 0.2), the effect of thionation in synergy with methylation opens a unique C-S bond cleavage pathway through crossing to a repulsive 1πσ* state, generating thiyl radicals as highly reactive intermediates that may invoke biological damage. This photodissociation channel is extremely difficult for conventional nucleobases. These findings demonstrate the synergistic effects of double functionality substitution in modulating excited-state dynamics and enhancing the photolabile character of DNA nucleobases, providing inspirations for the rational design of advanced photodynamic and photochemotherapy approaches.

Electronic spectroscopy of gemcitabine and derivatives for possible dual-action photodynamic therapy applications

In this paper, we explore the molecular basis of combining photodynamic therapy (PDT), a light-triggered targeted anticancer therapy, with the traditional chemotherapeutic properties of the well-known cytotoxic agent gemcitabine. A photosensitizer prerequisite is significant absorption of biocompatible light in the visible/near IR range, ideally between 600 and 1000 nm. We use highly accurate multiconfigurational CASSCF/MS-CASPT2/MM and TD-DFT methodologies to determine the absorption properties of a series of gemcitabine derivatives with the goal of red-shifting the UV absorption band toward the visible region and facilitating triplet state population. The choice of the substitutions and, thus, the rational design is based on important biochemical criteria and on derivatives whose synthesis is reported in the literature. The modifications tackled in this paper consist of: (i) substitution of the oxygen atom at O2 position with heavier atoms (O → S and O → Se) to red shift the absorption band and increase the spin-orbit coupling, (ii) addition of a lipophilic chain at the N7 position to enhance transport into cancer cells and slow down gemcitabine metabolism, and (iii) attachment of aromatic systems at C5 position to enhance red shift further. Results indicate that the combination of these three chemical modifications markedly shifts the absorption spectrum toward the 500 nm region and beyond and drastically increases spin-orbit coupling values, two key PDT requirements. The obtained theoretical predictions encourage biological studies to further develop this anticancer approach.

Ligand Dictated Photosensitization of Iridium(III) Dithiocarbamate Complexes for Photodynamic Therapy

Organelle-targeted photosensitizers (PSs) for photodynamic therapy (PDT) are considered as an effective therapeutic strategy for the development of next generation PSs with the least side effects and high therapeutic efficacy. However, multiorganelle targeted PSs eliciting PDT via both type I and type II mechanisms are scarce. Herein, a series of cyclometalated iridium(III) complexes were formulated [Ir(C∧N)2(S∧S)] (C∧N = 2-phenylpyridine (ppy) and 2-(thiophen-2-yl)pyridine (thpy); S∧S = diethyldithiocarbamate (DEDTC), morpholine-N-dithiocarbamate (MORDTC) and methoxycarbonodithioate (MEDTC)) and the newly designed complexes Ir2@DEDTC and Ir1@MEDTC were characterized by single crystal X-ray crystallography. Complexes containing thpy as C∧N ligand exhibit excellent photophysical properties such as red-shifted emission, high singlet oxygen quantum yield (ϕΔ) and longer photoluminescence lifetime when compared with complexes containing ppy ligands. Ir2@DEDTC exhibits the highest ϕΔ and photoluminescence lifetimes among the synthesized complexes. Therefore, Ir2@DEDTC was chosen to evaluate the photosensitizing ability to produce reactive oxygen species (ROS). Upon blue light irradiation (456 nm), it efficiently produces ROS, i.e., hydroxy radical (•OH) and singlet oxygen (1O2), which was confirmed by electron paramagnetic resonance (EPR) spectroscopy. In vitro photocytotoxicity toward HCT116, HeLa, and PC3 cell lines showed that out of all the synthesized complexes, Ir2@DEDTC has the highest photocytotoxic index (PI > 400) value. Ir2@DEDTC is efficiently taken up by the HCT116 cell line and accumulated mainly in the lysosome and mitochondria of the cells, and after PDT treatment, it elicits cell shrinkage, membrane blebbing, and DNA fragmentation. The phototherapeutic efficacy of Ir2@DEDTC has been investigated against 3D spheroids considering its ability to mimic some of the basic features of solid tumors. The morphology was drastically altered in the Ir2@DEDTC treated 3D spheroid after the light irradiation unleashed the potential of the Ir(III) dithiocarbamate complex as a superior PS for PDT. Hence, mitochondria and lysosome targeted photoactive cyclometalated Ir(III) dithiocarbamate complex exerting oxidative stress via both type I and type II PDT can be regarded as a dual-organelle targeted two-pronged approach for enhanced PDT.

The impact of the chalcogen-substitution element and initial spectroscopic state on excited-state relaxation pathways in nucleobase photosensitizers: a combination of static and dynamic studies

The substitution of oxygen with chalcogen in carbonyl group(s) of canonical nucleobases gives an impressive triplet generation, enabling their promising applications in medicine and other emerging techniques. The excited-state relaxation S2(ππ*) → S1(nπ*) → T1(ππ*) has been considered the preferred path for triplet generation in these nucleobase derivatives. Here, we demonstrate enhanced quantum efficiency of direct intersystem crossing from S2 to triplet manifold upon substitution with heavier chalcogen elements. The excited-state relaxation dynamics of sulfur/selenium substituted guanines in a vacuum is investigated using a combination of static quantum chemical calculations and on-the-fly excited-state molecular dynamics simulations. We find that in sulfur-substitution the S2 state predominantly decays to the S1 state, while upon selenium-substitution the S2 state deactivation leads to simultaneous population of the S1 and T2,3 states in the same time scale and multi-state quasi-degeneracy region S2/S1/T2,3. Interestingly, the ultrafast deactivation of the spectroscopic S3 state of both studied molecules to the S1 state occurs through a successive S3 → S2 → S1 path involving a multi-state quasi-degeneracy S3/S2/S1. The populated S1 and T2 states will cross the lowest triplet state, and the S1 → T intersystem crossing happens in a multi-state quasi-degeneracy region S1/T2,3/T1 and is accelerated by selenium-substitution. The present study reveals the influence of both the chalcogen substitution element and initial spectroscopic state on the excited-state relaxation mechanism of nucleobase photosensitizers and also highlights the important role of multi-state quasi-degeneracy in mediating the complex relaxation process. These theoretical results provide additional insights into the intrinsic photophysics of nucleobase-based photosensitizers and are helpful for designing novel photo-sensitizers for real applications.

New Gradient Correction Scheme for Electronically Nonadiabatic Dynamics Involving Multiple Spin States

It has been recommended that the best representation to use for trajectory surface hopping (TSH) calculations is the fully adiabatic basis in which the Hamiltonian is diagonal. Simulations of intersystem crossing processes with conventional TSH methods require an explicit computation of nonadiabatic coupling vectors (NACs) in the molecular-Coulomb-Hamiltonian (MCH) basis, also called the spin-orbit-free basis, in order to compute the gradient in the fully adiabatic basis (also called the diagonal representation). This explicit requirement destroys some of the advantages of the overlap-based algorithms and curvature-driven algorithms that can be used for the most efficient TSH calculations. Therefore, although these algorithms allow one to perform NAC-free simulations for internal conversion processes, one still requires NACs for intersystem crossing. Here, we show that how the NAC requirement is circumvented by a new computation scheme called the time-derivative-matrix scheme.

Photochemical and thermochemical pathways to S2 and polysulfur formation in the atmosphere of Venus

Polysulfur species have been proposed to be the unknown near-UV absorber in the atmosphere of Venus. Recent work argues that photolysis of one of the (SO)₂ isomers, cis-OSSO, directly yields S₂ with a branching ratio of about 10%. If correct, this pathway dominates polysulfur formation by several orders of magnitude, and by addition reactions yields significant quantities of S₃, S₄, and S₈. We report here the results of high-level ab-initio quantum-chemistry computations that demonstrate that S₂ is not a product in cis-OSSO photolysis. Instead, we establish a novel mechanism in which S₂ is formed in a two-step process. Firstly, the intermediate S₂O is produced by the coupling between the S and Cl atmospheric chemistries (in particular, SO reaction with ClS) and in a lesser extension by O-abstraction reactions from cis-OSSO. Secondly, S₂O reacts with SO. This modified chemistry yields S₂ and subsequent polysulfur abundances comparable to the photolytic cis-OSSO mechanism through a more plausible pathway. Ab initio quantification of the photodissociations at play fills a critical data void in current atmospheric models of Venus.

Polysulfur compounds have been ascribed as the unknown near-UV absorbers in Venusian atmosphere and play a key role in the sulfur chemical cycle of this planet. Here, authors establish their production from (SO)₂ on the grounds of quantifications of photochemical and thermal pathways involved in the sulfur chemical cycle of the planet.

Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water

Rationalizing the photochemistry of nucleobases where an oxygen is replaced by a heavier atom is essential for applications that exploit near-unity triplet quantum yields. Herein, we report on the ultrafast excited-state deactivation mechanism of 6-selenoguanine (6SeGua) in water by combining nonadiabatic trajectory surface-hopping dynamics with an electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. We find that the predominant relaxation mechanism after irradiation starts on the bright singlet S₂ state that converts internally to the dark S₁ state, from which the population is transferred to the triplet T₂ state via intersystem crossing and finally to the lowest T₁ state. This S₂ → S₁ → T₂ → T₁ deactivation pathway is similar to that observed for the lighter 6-thioguanine (6tGua) analogue, but counterintuitively, the T₁ lifetime of the heavier 6SeGua is shorter than that of 6tGua. This fact is explained by the smaller activation barrier to reach the T₁/S₀ crossing point and the larger spin-orbit couplings of 6SeGua compared to 6tGua. From the dynamical simulations, we also calculate transient absorption spectra (TAS), which provide two time constants (τ₁ = 131 fs and τ₂ = 191 fs) that are in excellent agreement with the experimentally reported value (τ_exp = 130 ± 50 fs) (Farrel et al. J. Am. Chem. Soc. 2018, 140, 11214). Intersystem crossing itself is calculated to occur with a time scale of 452 ± 38 fs, highlighting that the TAS is the result of a complex average of signals coming from different nonradiative processes and not intersystem crossing alone.

Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems

Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.

Dynamics of Excitons in Conjugated Molecules and Organic Semiconductor Systems

The exciton, an excited electron-hole pair bound by Coulomb attraction, plays a key role in photophysics of organic molecules and drives practically important phenomena such as photoinduced mechanical motions of a molecule, photochemical conversions, energy transfer, generation of free charge carriers, etc. Its behavior in extended π-conjugated molecules and disordered organic films is very different and very rich compared with exciton behavior in inorganic semiconductor crystals. Due to the high degree of variability of organic systems themselves, the exciton not only exerts changes on molecules that carry it but undergoes its own changes during all phases of its lifetime, that is, birth, conversion and transport, and decay. The goal of this review is to give a systematic and comprehensive view on exciton behavior in π-conjugated molecules and molecular assemblies at all phases of exciton evolution with emphasis on rates typical for this dynamic picture and various consequences of the above dynamics. To uncover the rich variety of exciton behavior, details of exciton formation, exciton transport, exciton energy conversion, direct and reverse intersystem crossing, and radiative and nonradiative decay are considered in different systems, where these processes lead to or are influenced by static and dynamic disorder, charge distribution symmetry breaking, photoinduced reactions, electron and proton transfer, structural rearrangements, exciton coupling with vibrations and intermediate particles, and exciton dissociation and annihilation as well.

Ultrafast Excited-State Decay Mechanisms of 6-Thioguanine Followed by Sub-20 fs UV Transient Absorption Spectroscopy

Understanding the primary steps following UV photoexcitation in sulphur-substituted DNA bases (thiobases) is fundamental for developing new phototherapeutic drugs. However, the investigation of the excited-state dynamics in sub-100 fs time scales has been elusive until now due to technical challenges. Here, we track the ultrafast decay mechanisms that lead to the electron trapping in the triplet manifold for 6-thioguanine in an aqueous solution, using broadband transient absorption spectroscopy with a sub-20 fs temporal resolution. We obtain experimental evidence of the fast internal conversion from the S₂(ππ*) to the S₁(nπ*) states, which takes place in about 80 fs and demonstrates that the S₁(nπ*) state acts as a doorway to the triplet population in 522 fs. Our results are supported by MS-CASPT2 calculations, predicting a planar S₂(ππ*) pseudo-minimum in agreement with the stimulated emission signal observed in the experiment.

Mechanistic photophysics and photochemistry of unnatural bases and sunscreen molecules: insights from electronic structure calculations

Photophysics and photochemistry are basic subjects in the study of light-matter interactions and are ubiquitous in diverse fields such as biology, energy, materials, and environment. A full understanding of mechanistic photophysics and photochemistry underpins many recent advances and applications. This contribution first provides a short discussion on the theoretical calculation methods we have used in relevant studies, then we introduce our latest progress on the mechanistic photophysics and photochemistry of two classes of molecular systems, namely unnatural bases and sunscreens. For unnatural bases, we disclose the intrinsic driving forces for the ultrafast population to reactive triplet states, impacts of the position and degree of chalcogen substitutions, and the effects of complex environments. For sunscreen molecules, we reveal the photoprotection mechanisms that dissipate excess photon energy to the surroundings by ultrafast internal conversion to the ground state. Finally, relevant theoretical challenges and outlooks are discussed.

Intersystem crossing and internal conversion dynamics with GAIMS-TeraChem: Excited state relaxation in 2-cyclopentenone

Excited states relaxation in complex molecules often involves two types of nonradiative transitions, internal conversion (IC) and intersystem crossing (ISC). In the situations when the timescales of IC and ISC are comparable, an interplay between these two types of transitions can lead to complex nonadiabatic dynamics on multiple electronic states of different characters and spin multiplicities. We demonstrate that the generalized ab initio multiple spawning (GAIMS) method interfaced with the fast graphics processing unit-based TeraChem electronic structure code can be used to model such nonadiabatic dynamics involving both the IC and ISC transitions in molecules of moderate size. We carried out 1500 fs GAIMS simulations leading to the creation of up to 2500 trajectory basis functions to study the excited states relaxation in 2-cyclopentenone. After a vertical excitation from the ground state to the bright S₂ state, the molecule quickly relaxes to the S₁ state via conical intersection. The following relaxation proceeds along two competing pathways: one involves IC to the ground state, and the other is dominated by ISC to the low-lying triplet states. The time constants describing the population transfer between the six lowest singlet and triplet states predicted by the GAIMS dynamics are in good agreement with the characteristic times of IC and ISC obtained from the analysis of the time-resolved photoelectron spectrum.

Mechanistic Photophysics of Tellurium-Substituted Uracils: Insights from Multistate Complete-Active-Space Second-Order Perturbation Calculations

The photophysical mechanisms of tellurium-substituted uracils were studied at the multistate complete-active-space second-order perturbation level with a particular focus on how the position and number of tellurium substitutions affect their nonadiabatic relaxation processes. Electronic structure analysis reveals that the lowest several excited states are closely concerned with the n and π orbitals at the Te₇-C₂ [Te₈-C₄] moiety of 2-tellurouracil (2TeU) [4TeU and 24TeU]. Both planar and twisted minima were optimized for 2TeU, whereas only planar ones were obtained for 4TeU and 24TeU, except for a twisted T₁ minimum of 4TeU. Based on intersection structures and linearly interpolated internal coordinate paths, we proposed several feasible excited-state deactivation paths. It is found that the relaxation channels for 2TeU are more complicated than those of 4TeU and 24TeU. The electronic population transfer to the T₁ state for 2TeU is easier than that for 4TeU and 24TeU in consideration of the barrier heights from the S₂ Franck-Condon point to the S₂/S₁ or S₂/T₂ intersections. In addition, the recovery of the ground state from the T₁ state for 2TeU will be more efficient than that for the other two systems as well.

An efficient implementation of spin-orbit coupling within the framework of semiempirical orthogonalization-corrected methods for ultrafast intersystem crossing dynamics

We implement spin-orbit coupling (SOC) within the framework of semiempirical orthogonalization-corrected methods (OMx). The excited-state wavefunction is generated from configuration interaction with single excitations (CIS). The SOC Hamiltonian in terms of the one-electron Breit-Pauli operator with effective nuclear charges is adopted in this work. Benchmark calculations show that SOCs evaluated using the OMx/CIS method agree very well with those obtained from time-dependent density functional theory. As a particularly attractive application, we incorporate SOCs between singlet and triplet states into Tully's fewest switches surface hopping algorithm to enable excited-state nonadiabatic dynamics simulations, treating internal conversion and intersystem crossing on an equal footing. This semiempirical dynamics simulation approach is applied to investigate ultrafast intersystem crossing processes in core-substituted naphthalenediimides.

Intramolecular Charge Transfer in the Azathioprine Prodrug Quenches Intersystem Crossing to the Reactive Triplet State in 6-Mercaptopurine†

The thiopurine prodrugs 6-mercaptopurine and azathioprine are among the world's essential medications for acute lymphoblastic leukemia, immunosuppression and several autoimmune conditions. Thiopurine prodrugs are efficient UVA absorbers and singlet oxygen generators and the long-term treatment with these prodrugs correlates with a high incidence of sunlight-induced skin cancer in patients. In this contribution, we show that the electronic relaxation mechanisms and photochemical properties of azathioprine are remarkably different from those of 6-mercaptopurine upon absorption of UVA radiation. UVA excitation of 6-mercaptopurine results in nearly 100% triplet yield and up to 30% singlet oxygen generation, whereas excitation of azathioprine with UVA leads to triplet yields of 15-3% depending on pH of the aqueous solution and <1% singlet oxygen generation. While photoexcitation of 6-mercaptopurine and other thiopurine prodrugs can facilitate oxidatively generated cell damage, azathioprine's poor photosensitization ability reveals the use of interchromophoric charge-transfer interactions for the molecular design of photostable prodrugs exhibiting a remarkable reduction in photocytotoxic side effects before drug metabolization.

Ribose Alters the Photochemical Properties of the Nucleobase in Thionated Nucleosides

Substitution of exocyclic oxygen with sulfur was shown to substantially influence the properties of RNA/DNA bases, which are crucial for prebiotic chemistry and photodynamic therapies. Upon UV irradiation, thionucleobases were shown to efficiently populate triplet excited states and can be involved in characteristic photochemistry or generation of singlet oxygen. Here, we show that the photochemistry of a thionucleobase can be considerably modified in a nucleoside, that is, by the presence of ribose. Our transient absorption spectroscopy experiments demonstrate that thiocytosine exhibits 5 times longer excited-state lifetime and different excited-state absorption features than thiocytidine. On the basis of accurate quantum chemical simulations, we assign these differences to the dominant population of a shorter-lived triplet nπ* state in the nucleoside and longer-lived triplet ππ* states in the nucleobase. This explains the distinctive photoanomerziation of thiocytidine and indicates that the nucleoside will be a less efficient phototherapeutic agent with regard to singlet oxygen generation.

Photochemistry and Non-adiabatic Photodynamics of the HOSO Radical

Hydroxysulfinyl radical (HOSO) is important due to its involvement in climate geoengineering upon SO₂ injection and generation of the highly hygroscopic H₂SO₄. Its photochemical behavior in the upper atmosphere is, however, uncertain. Here we present the ultraviolet-visible photochemistry and photodynamics of this species by simulating the atmospheric conditions with high-level quantum chemistry methods. Photocleavage to HO + SO arises as the major solar-induced channel, with a minor contribution of H + SO₂ photoproducts. The efficient generation of SO is relevant due to its reactivity with O₃ and the consequent depletion of ozone in the stratosphere.

MS-CASPT2 studies on the mechanistic photophysics of tellurium-substituted guanine and cytosine

Sulfur-substituted nucleobases are highly promising photosensitizers that are widely used in photodynamic therapy, and there are numerous studies exploring their unique photophysical behaviors. However, relevant photophysical investigations on selenium and tellurium substitutions are still rare. Herein, the high-level multistate complete-active-space second-order perturbation (MS-CASPT2) method was performed for the first time to explore the excited-state relaxation processes of tellurium-substituted guanine (TeG) and cytosine (TeC). Based on the electronic state properties in the Franck-Condon (FC) region, we found that the lowest five (S0, S1, S2, T1, and T2) and six (S0, S1, S2, T1, T2 and T3) states will participate in the nonadiabatic transition processes of TeG and TeC systems, respectively. In these electronic states, two kinds of minimum and intersection structures (i.e., planar and twisted structures) were obtained for both TeG and TeC systems. The linearly interpolated internal coordinate (LIIC) paths and spin-orbit coupling (SOC) constants revealed several possible planar and twisted excited-state decay channels, which could lead the systems to the lowest reactive triplet state of T1. Small energy barriers in the T1 state will trap the TeG and TeC systems for a while before they finally populate to the ground state. Although tellurium substitution would further redshift the absorption wavelength and enhance the intersystem crossing (ISC) rate to the T1 state compared with sulfur and selenium substitutions, the rapid ISC process of T1 → S0 may make it a less effective photosensitizer to sensitize the molecular oxygen. We believe our present work will provide important mechanistic insights into the photophysics of tellurium-substituted nucleobases.

Modeling Spin-Crossover Dynamics

In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.

On the population of triplet states of 2-seleno-thymine

The population and depopulation mechanisms leading to the lowest-lying triplet states of 2-Se-Thymine were studied at the MS-CASPT2/cc-pVDZ level of theory. Several critical points on different potential energy hypersurfaces were optimized, including minima, conical intersections, and singlet-triplet crossings. The accessibility of all relevant regions on the potential energy hypersurfaces was investigated by means of minimum energy paths and linear interpolation in internal coordinates techniques. Our analysis indicates that, after the population of the bright S2 state in the Franck-Condon region, the first photochemical event is a barrierless evolution towards one of its two minima. After that, three viable photophysical deactivation paths can take place. In one of them, the population in the S2 state is transferred to the T2 state via intersystem crossing and subsequently to the T1 state by internal conversion. Alternatively, the S1 state could be accessed by internal conversion through two distinct conical intersections with S2 state followed by singlet-triplet crossing with the T2 state. The absence of a second minimum on the T1 state and a small energy barrier on pathway along the potential energy surface towards the ground state from the lowest triplet state are attributed as potential reasons to explain why the lifetime of the triplet state of 2-Se-Thymine might be reduced in comparison with its thio-analogue.

Thioguanine restoration through type I photosensitization-superoxide oxidation-glutathione reduction cycles

UVA-induced deleterious effect of thiopurine prodrugs including azathioprine, 6-mercaptopurine and 6-thioguanine (6-TG) increases the risk of cancer development due to the incorporation of 6-TG in patients' DNA. The catalytic mechanism by which thiobases act as a sustained oxidant producer has yet to be explored, especially through the Type I electron transfer pathway that produces superoxide radicals (O2˙-). Under Fenton-like conditions O2˙- radicals convert to extremely reactive hydroxyl radicals (˙OH), thus carrying even higher risk of biological damage than that induced by the well-studied type II reaction. By monitoring 6-TG/UVA-induced photochemistry in mass spectra and superoxide radicals (O2˙-) via nitro blue tetrazolium (NBT) reduction, this work provides two new findings: (1) in the presence of reduced glutathione (GSH), the production of O2˙-via the type I reaction is enhanced 10-fold. 6-TG thiyl radicals are identified as the primary intermediate formed in the reaction of 6-TG with O2˙-. The restoration of 6-TG and concurrent generation of O2˙- occur via a 3-step-cycle: 6-TG type I photosensitization, O2˙- oxidation and GSH reduction. (2) In the absence of GSH, 6-TG thiyl radicals undergo oxygen addition and sulfur dioxide removal to form carbon radicals (C6) which further convert to thioether by reacting with 6-TG molecules. These findings help explain not only thiol-regulation in a biological system but chemoprevention of cancer.

Behind the scenes of spin-forbidden decay pathways in transition metal complexes

The interpretation of the ultrafast photophysics of transition metal complexes following photo-absorption is quite involved as the heavy metal center leads to a complicated and entangled singlet-triplet manifold. This opens up multiple pathways for deactivation, often with competitive rates. As a result, intersystem crossing (ISC) and phosphorescence are commonly observed in transition metal complexes. A detailed understanding of such an excited-state structure and dynamics calls for state-of-the-art experimental and theoretical methodologies. In this review, we delve into the inability of non-relativistic quantum theory to describe spin-forbidden transitions, which can be overcome by taking into account spin-orbit coupling, whose importance grows with increasing atomic number. We present the quantum chemical theory of phosphorescence and ISC together with illustrative examples. Finally, a few applications are highlighted, bridging the gap between theoretical studies and experimental applications, such as photofunctional materials.

Internal Conversion and Intersystem Crossing with the Exact Factorization

We present a detailed derivation of the generalized coupled-trajectory mixed quantum-classical (G-CT-MQC) algorithm based on the exact-factorization equations. The ultimate goal is to propose an algorithm that can be employed for molecular dynamics simulations of nonradiative phenomena, as the spin-allowed internal conversions and the spin-forbidden intersystem crossings. Internal conversions are nonadiabatic processes driven by the kinetic coupling between electronic states, whereas intersystem crossings are mediated by the spin-orbit coupling. In this paper, we discuss computational issues related to the suitable representation for electronic dynamics and the different natures of kinetic and spin-orbit coupling. Numerical studies on model systems allow us to test the performance of the G-CT-MQC algorithm in different situations.

Single-atom replacement as a general approach towards visible-light/near-infrared heavy-atom-free photosensitizers for photodynamic therapy

Photodynamic therapy has become an emerging strategy for the treatment of cancer. This technology relies on the development of photosensitizers (PSs) that convert molecular oxygen to cytotoxic reactive oxygen species upon exposure to light. In this study, we have developed a facile and general strategy for obtaining visible light/near-infrared-absorbing PSs by performing a simple sulfur-for-oxygen replacement within existing fluorophores. Thionation of carbonyl groups within existing fluorophore cores leads to an improvement of the singlet oxygen quantum yield and molar absorption coefficient at longer wavelengths (deep to 600-800 nm). Additionally, these thio-based PSs lack dark cytotoxicity but exhibit significant phototoxicity against monolayer cancer cells and 3D multicellular tumor spheroids with IC50 in the micromolar range. To achieve tumor-specific delivery, we have conjugated these thio-based PSs to an antibody and demonstrated their tumor-specific therapeutic activity.

Selenium substitution effects on excited-state properties and photophysics of uracil: a MS-CASPT2 study

The photophysics of selenium-substituted nucleobases has attracted recent experimental attention because they could serve as potential photosensitizers in photodynamic therapy. Herein, we present a comprehensive MS-CASPT2 study on the spectroscopic and excited-state properties, and photophysics of 2-selenouracil (2SeU), 4-selenouracil (4SeU), and 2,4-selenouracil (24SeU). Relevant minima, conical intersections, crossing points, and excited-state relaxation paths in the lowest five electronic states (i.e., S₀, S₁, S₂, T₂, and T₁) are explored. On the basis of these results, their photophysical mechanisms are proposed. Upon photoirradiation to the bright S₂ state, 2SeU quickly relaxes to its S₂ minimum and then moves in an essentially barrierless way to a nearby S₂/S₁ conical intersection near which the S₁ state is populated. Next, the S₁ system arrives at an S₁/T₂/T₁ intersection where a large S₁/T₁ spin-orbit coupling of 430.8 cm^-1 makes the T₁ state populated. In this state, a barrier of 6.8 kcal mol^-1 will trap 2SeU for a while. In parallel, for 4SeU or 24SeU, the system first relaxes to the S₂ minimum and then overcomes a small barrier to approach an S₂/S₁ conical intersection. Once hopping to the S₁ state, there exists an extended region with very close S₁, T₂, and T₁ energies. Similarly, a large S₁/T₁ spin-orbit coupling of 426.8 cm^-1 drives the S₁→ T₁ intersystem crossing process thereby making the T₁ state populated. Similarly, an energy barrier heavily suppresses electronic transition to the S₀ state. The present work manifests that different selenium substitutions on uracil can lead to a certain extent of different vertical and adiabatic excitation energies, excited-state properties, and relaxation pathways. These insights could help understand the photophysics of selenium-substituted nucleobases.

Simulating the Nonadiabatic Relaxation Dynamics of 4-(N,N-Dimethylamino)benzonitrile (DMABN) in Polar Solution

The compound 4-(N,N-dimethylamino)benzonitrile (DMABN) represents the archetypal system for dual fluorescence, a rare photophysical phenomenon in which a given fluorophore shows two distinct emission bands. Despite extensive studies, the underlying mechanism remains the subject of debate. In the present contribution, we address this issue by simulating the excited-state relaxation process of DMABN as it occurs in polar solution. The potential energy surfaces for the system are constructed with the use of the additive quantum mechanics/molecular mechanics (QM/MM) method, and the coupled dynamics of the electronic wave function and the nuclei is propagated with the semiclassical fewest switches surface hopping method. The DMABN molecule, which comprises the QM subsystem, is treated with the use of the second-order algebraic diagrammatic construction (ADC(2)) method with the imposition of spin-opposite scaling (SOS). It is verified that this level of theory achieves a realistic description of the excited-state potential energy surfaces of DMABN. The simulation results qualitatively reproduce the main features of the experimentally observed fluorescence spectrum, thus allowing the unambiguous assignment of the two fluorescence bands: the normal band is due to the near-planar locally excited (LE) structure of DMABN, while the so-called "anomalous" second band arises from the twisted intramolecular charge transfer (TICT) structure. The transformation of the LE structure into the TICT structure takes place directly via intramolecular rotation, and is not mediated by another excited-state structure. In particular, the oft-discussed rehybridized intramolecular charge transfer (RICT) structure, which is characterized by a bent nitrile group, does not play a role in the relaxation process.

Ultrafast dynamics of photo-excited 2-thiopyridone: Theoretical insights into triplet state population and proton transfer pathways

Ultrafast non-adiabatic dynamics of the small heteroaromatic compound 2-thiopyridone has been studied with surface hopping simulations based on multi-configurational quantum chemistry. Initial excitation of the bright S 2 ( π , π * ) state is found to promptly relax to S 1 ( n ,   π * ) through in-plane motion. The subsequent dynamics are oppositely driven by out-of-plane motion, which results in both complex population transfers among all the available states and intersystem crossing predominantly through the "El-Sayed forbidden" S 1 ( n ,   π * ) to T 2 ( n ,   π * ) channel, through significant mixing of electronic excitation characters. Despite this complexity, the femto- to picosecond triplet population, expected from several spectroscopic measurements, is well described as a simple exponential decay of the singlet state manifold. No proton transfer is found in the reported trajectories, but two mechanisms for its possible mediation in previously reported experiments are proposed based on the observed structural dynamics: (i) ultrafast intra-molecular transfer driven by the initially coherent in-plane motion and (ii) inter-molecular solvent-mediated transfer driven by the out-of-plane modes that dominate the later motion.

Excited State Lifetimes of Sulfur-Substituted DNA and RNA Monomers Probed Using the Femtosecond Fluorescence Up-Conversion Technique

Sulfur-substituted DNA and RNA nucleobase derivatives (a.k.a., thiobases) are an important family of biomolecules. They are used as prodrugs and as chemotherapeutic agents in medical settings, and as photocrosslinker molecules in structural-biology applications. Remarkably, excitation of thiobases with ultraviolet to near-visible light results in the population of long-lived and reactive triplet states on a time scale of hundreds of femtoseconds and with near-unity yields. This efficient nonradiative decay pathway explains the vanishingly small fluorescence yields reported for the thiobases and the scarcity of fluorescence lifetimes in the literature. In this study, we report fluorescence lifetimes for twelve thiobase derivatives, both in aqueous solution at physiological pH and in acetonitrile. Excitation is performed at 267 and 362 nm, while fluorescence emission is detected at 380, 425, 450, 525, or 532 nm. All the investigated thiobases reveal fluorescence lifetimes that decay in a few hundreds of femtoseconds and with magnitudes that depend and are sensitive to the position and degree of sulfur-atom substitution and on the solvent environment. Interestingly, however, three thiopyrimidine derivatives (i.e., 2-thiocytidine, 2-thiouridine, and 4-thiothymidine) also exhibit a small amplitude fluorescence component of a few picoseconds in aqueous solution. Furthermore, the N-glycosylation of thiobases to form DNA or RNA nucleoside analogues is demonstrated as affecting their fluorescence lifetimes. In aqueous solution, the fluorescence decay signals exciting at 267 nm are equal or slower than those collected exciting at 362 nm. In acetonitrile, however, the fluorescence decay signals recorded upon 267 nm excitation are, in all cases, faster than those measured exciting at 362 nm. A comparison to the literature values show that, while both the DNA and RNA nucleobase and thiobase derivatives exhibit sub-picosecond fluorescence lifetimes, the 1ππ* excited-state population in the nucleobase monomers primarily decay back to the ground state, whereas it predominantly populates long-lived and reactive triplet states in thiobase monomers.

Spin-Orbit Interactions in Ultrafast Molecular Processes

We investigate spin-orbit interactions in ultrafast molecular processes employing the exact factorization of the electron-nuclear wave function. We revisit the original derivation by including spin-orbit coupling, and show how the dynamics driven by the time-dependent potential energy surface alleviates inconsistencies arising from different electronic representations. We propose a novel trajectory-based scheme to simulate spin-forbidden non-radiative processes, and we show its performance in the treatment of excited-state dynamics where spin-orbit effects couple different spin multiplets.

Intersystem crossing in tunneling regime: T1 → S0 relaxation in thiophosgene

The T1 excited state relaxation in thiophosgene has attracted much attention as a relatively simple model for the intersystem crossing (ISC) transitions in polyatomic molecules. The very short (20-40 ps) T1 lifetime predicted in several theoretical studies strongly disagrees with the experimental values (∼20 ns) indicating that the kinetics of T1 → S0 ISC is not well understood. We use the nonadiabatic transition state theory (NA-TST) with the Zhu-Nakamura transition probability and the multireference perturbation theory (CASPT2) to show that the T1 → S0 ISC occurs in the quantum tunneling regime. We also introduce a new zero-point vibrational energy correction scheme that improves the accuracy of the predicted ISC rate constants at low internal energies. The predicted lifetimes of the T1 vibrational states are between one and two orders of magnitude larger than the experimental values. This overestimation is attributed to the multidimensional nature of quantum tunneling that facilitates ISC transitions along the non-minimum energy path and is not accounted for in the one-dimensional NA-TST.

MS-CASPT2 Studies on the Photophysics of Selenium-Substituted Guanine Nucleobase

The MS-CASPT2 method has been employed to optimize minimum-energy structures of 6-selenoguanine (6SeGua) and related two- and three-state intersection structures in and between the lowest five electronic states, i.e., S₂(¹ππ*), S₁(¹ nπ*), T₂(³ nπ*), T₁(³ππ*), and S₀. In combination with MS-CASPT2 calculated linearly interpolated internal coordinate paths, the photophysical mechanism of 6SeGua has been proposed. The initially populated S₂(¹ππ*) state decays to either S₁(¹ nπ*) or T₂(³ nπ*) states through a three-state S₂/S₁/T₂ intersection point. The large S₂/T₂ spin-orbit coupling of 435 cm^-1, according to the classical El-Sayed rule, benefits the S₂ → T₂ intersystem crossing process. The S₁(¹ nπ*) state that stems from the S₂ → S₁ internal conversion process at the S₂/S₁/T₂ intersection point can further jump to the T₂(³ nπ*) state through the S₁ → T₂ intersystem crossing process. This process does not comply with the El-Sayed rule, but it is still related to a comparatively large spin-orbit coupling of 39 cm^-1 and is expected to occur relatively fast. Finally, the T₂(³ nπ*) state, which is populated from the above S₂ → T₂ and S₁ → T₂ intersystem crossing processes, decays to the T₁(³ππ*) state via an internal conversion process. Because there is merely a small energy barrier of 0.11 eV separating the T₁(³ππ*) minimum and an energetically allowed two-state T₁/S₀ intersection point, the T₁(³ππ*) state still can decay to the S₀ state quickly, which is also enhanced by a large T₁/S₀ spin-orbit coupling of 252 cm^-1. Our proposed mechanism explains experimentally observed ultrafast intersystem crossing processes in 6SeGua and its 835-fold acceleration of the T₁ state decay to the S₀ state compared with 6tGua. Finally, we have found that the ground-state electronic structure of 6SeGua has more apparent multireference character.

Excited-State Dynamics of Thienoguanosine, an Isomorphic Highly Fluorescent Analogue of Guanosine

Thienoguanosine (^th G) is an isomorphic analogue of guanosine with promising potentialities as fluorescent DNA label. As a free probe in protic solvents, ^th G exists in two tautomeric forms, identified as the H1, being the only one observed in nonprotic solvents, and H3 keto-amino tautomers. We herein investigate the photophysics of ^th G in solvents of different polarity, from water to dioxane, by combining time-resolved fluorescence with PCM/TD-DFT and CASSCF calculations. Fluorescence lifetimes of 14.5-20.5 and 7-13 ns were observed for the H1 and H3 tautomers, respectively, in the tested solvents. In methanol and ethanol, an additional fluorescent decay lifetime (≈3 ns) at the blue emission side (λ≈430 nm) as well as a 0.5 ns component with negative amplitude at the red edge of the spectrum, typical of an excited-state reaction, were observed. Our computational analysis explains the solvent effects observed on the tautomeric equilibrium. The main radiative and nonradiative deactivation routes have been mapped by PCM/TD-DFT calculations in solution and CASSCF in the gas phase. The most easily accessible conical intersection, involving an out-of plane motion of the sulfur atom in the five-membered ring of ^th G, is separated by a sizeable energy barrier (≥0.4 eV) from the minimum of the spectroscopic state, which explains the large experimental fluorescence quantum yield.

Photostability of oxazoline RNA-precursors in UV-rich prebiotic environments

Pentose aminooxazolines and oxazolidinone thiones are considered as the key precursors which could have enabled the formation of RNA nucleotides under the conditions of early Earth. UV-irradiation experiments and quantum-chemical calculations demonstrate that these compounds are remarkably photostable and could accumulate over long periods of time in UV-rich prebiotic environments to undergo stereoisomeric purification.

Nonadiabatic dynamics: The SHARC approach

We review the Surface Hopping including ARbitrary Couplings (SHARC) approach for excited-state nonadiabatic dynamics simulations. As a generalization of the popular surface hopping method, SHARC allows simulating the full-dimensional dynamics of molecules including any type of coupling terms beyond nonadiabatic couplings. Examples of these arbitrary couplings include spin-orbit couplings or dipole moment-laser field couplings, such that SHARC can describe ultrafast internal conversion, intersystem crossing, and radiative processes. The key step of the SHARC approach consists of a diagonalization of the Hamiltonian including these couplings, such that the nuclear dynamics is carried out on potential energy surfaces including the effects of the couplings-this is critical in any applications considering, for example, transition metal complexes or strong laser fields. We also give an overview over the new SHARC2.0 dynamics software package, released under the GNU General Public License, which implements the SHARC approach and several analysis tools. The review closes with a brief survey of applications where SHARC was employed to study the nonadiabatic dynamics of a wide range of molecular systems. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and KineticsSoftware > Simulation MethodsSoftware > Quantum Chemistry.