Excited-State Properties and Relaxation Pathways of Selenium-Substituted Guanine Nucleobase in Aqueous Solution and DNA Duplex

Thionated Coumarins: Study of the Intersystem Crossing and the Zero-field Splitting of the Triplet State Using Time-Resolved Transient Optical and Electron Paramagnetic Resonance Spectroscopies

To study the effect of thionation of the carbonyl groups in a chromophore, i. e. replacing the O atom with S atom, on the photophysics, we studied two thionated coumarin derivatives (Cou-S and Cou-6-S) with various steady state and transient spectroscopic methods. Both compounds exhibit red-shifted absorption (up to 4900 cm-1) and strong fluorescence quenching as compared to the unthionated analogues. Femtosecond transient absorption spectra show fast ISC (ca. 10 ps) in the thionated coumarin derivatives, while negligible ISC was observed in the unthionated coumarin. Interestingly, triplet excited state lifetimes of the thionated coumarin (0.14 μs) is much shorter than the unthionated analogues (53.4 μs). Time-resolved electron paramagnetic resonance (TREPR) spectra indicate much larger zero field splitting (ZFS) D parameters (up to 0.287 cm-1) for the T1 state of the thionated coumarins than the unthionated analogues (D=0.1001 cm-1). This large D value is attributed to the strong spin orbital coupling effect. These results demonstrate the advantage and the drawback of thionation-enhanced ISC, i. e. the ISC is efficient, but triplet state lifetimes become substantially shorter. This information is useful for the future design of heavy atom-free triplet photosensitizers for photodynamic therapy, photon upconversion, photocatalytic organic synthesis and photopolymerization, etc.

Analytic First-Order Derivatives of CASPT2 Combined with the Polarizable Continuum Model

The complete active space second-order perturbation theory (CASPT2) is valuable for accurately predicting electronic structures and transition energies. However, optimizing molecular geometries in the solution phase has proven challenging. In this study, we develop analytic first-order derivatives of CASPT2 using an implicit solvation model, specifically the polarizable continuum model, within the open-source package OpenMolcas. Analytic gradients and nonadiabatic coupling vectors are computed by solving a modified Z-vector equation. Comparisons with existing theoretical and experimental results demonstrate that the solvent effects can be qualitatively captured using the developed method.

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.

Mechanistic photophysics of tellurium-substituted cytosine: Electronic structure calculations and nonadiabatic dynamics simulations

Previously, the MS-CASPT2 method was performed to study the static and qualitative photophysics of tellurium-substituted cytosine (TeC). To get quantitative information, we used our recently developed QTMF-FSSH dynamics method to simulate the excited-state decay of TeC. The CASSCF method was adopted to reduce the calculation costs, which was confirmed to provide reliable structures and energies as those of MS-CASPT2. A detailed structural analysis showed that only 5% trajectories will hop to the lower triplet or singlet state via the twisted (S2 /S1 /T2 )T intersection, while 67% trajectories will choose the planar intersections of (S2 /S1 /T3 /T2 /T1 )P and (S2 /S1 /T2 /T1 )P but subsequently become twisted in other electronic states. By contrast, ~28% trajectories will maintain in a plane throughout dynamics. Electronic population revealed that the S2 population will ultrafast transfer to the lower triplet or singlet state. Later, the TeC system will populate in the spin-mixed electronic states composed of S1 , T1 and T2 . At the end of 300 fs, most trajectories (~74%) will decay to the ground state and only 17.4% will survive in the triplet states. Our dynamics simulation verified that tellurium substitution will enhance the intersystem crossings, but the very short triplet lifetime (ca. 125 fs) will make TeC a less effective photosensitizer.

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.

Photophysics of tz Adenine and tz Guanine fluorescent nucleobases embedded into DNA and RNA

UV-VIS photoinduced events of tz A and tz G embedded into DNA and RNA are described by combining the Extended Multi-State Second-Order Perturbation Theory (XMS-CASPT2) and electrostatic embedding molecular mechanics methods (QM/MM). Our results point out that the S1 1 (ππ* La ) state is the bright state in both environments. After the photoexcitation to the S1 1 (ππ* La ) state, the electronic population evolves barrierless towards its minimum, from where the excess of energy can be dissipated by fluorescence. As the minimum energy crossing point structure between the ground and first bright states lies in a high-energy region, the direct internal conversion to the ground state is an unviable mechanism. Other spectroscopic properties (for instance, absorption and Stokes shifts) and comparisons with photochemical properties of canonical nucleobases are also provided.

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.

Finding Excited-State Minimum Energy Crossing Points on a Budget: Non-Self-Consistent Tight-Binding Methods

The automated exploration and identification of minimum energy conical intersections (MECIs) is a valuable computational strategy for the study of photochemical processes. Due to the immense computational effort involved in calculating non-adiabatic derivative coupling vectors, simplifications have been introduced focusing instead on minimum energy crossing points (MECPs), where promising attempts were made with semiempirical quantum mechanical methods. A simplified treatment for describing crossing points between almost arbitrary diabatic states based on a non-self-consistent extended tight-binding method, GFN0-xTB, is presented. Involving only a single diagonalization of the Hamiltonian, the method can provide energies and gradients for multiple electronic states, which can be used in a derivative coupling-vector-free scheme to calculate MECPs. By comparison with high-lying MECIs of benchmark systems, it is demonstrated that the identified geometries are good starting points for further MECI refinement with ab initio methods.

Sulphur- and Selenium-for-Oxygen Replacement as a Strategy to Obtain Dual Type I/Type II Photosensitizers for Photodynamic Therapy

The effect on the photophysical properties of sulfur- and selenium-for-oxygen replacement in the skeleton of the oxo-4-dimethylaminonaphthalimide molecule (DMNP) has been explored at the density functional (DFT) level of theory. Structural parameters, excitation energies, singlet–triplet energy gaps (ΔES-T), and spin–orbit coupling constants (SOC) have been computed. The determined SOCs indicate an enhanced probability of intersystem crossing (ISC) in both the thio- and seleno-derivatives (SDMNP and SeDMNP, respectively) and, consequently, an enhancement of the singlet oxygen quantum yields. Inspection of Type I reactions reveals that the electron transfer mechanisms leading to the generation of superoxide is feasible for all the compounds, suggesting a dual Type I/Type II activity.

Resolving Ultrafast Photoinitiated Dynamics of the Hachimoji 5-Aza-7-Deazaguanine Nucleobase: Impact of Synthetically Expanding the Genetic Alphabet

The guanine derivative, 5-aza-7-deazaguanine (^5N7C G) has recently been proposed as one of four unnatural bases, termed Hachimoji (8-letter) to expand the genetic code. We apply steady-state and time-resolved spectroscopy to investigate its electronic relaxation mechanism and probe the effect of atom substitution on the relaxation mechanism in polar protic and polar aprotic solvents. Mapping of the excited state potential energy surfaces is performed, from which the critical points are optimized by using the state-of-art Extended Multi-State Complete Active Space Second-Order Perturbation Theory. It is demonstrated that excitation to the lowest energy ¹ ππ* state of ^5N7C G results in complex dynamics leading to ca. 10 to 30-fold slower relaxation (depending on solvent) compared to guanine. A significant conformational change occurs at the S₁ minimum, resulting in a 10-fold greater fluorescence quantum yield compared to guanine. The fluorescence quantum yield and S₁ decay lifetime increase going from water to acetonitrile to propanol. The solvent-dependent results are supported by the quantum chemical calculations showing an increase in the energy barrier between the S₁ minimum and the S₁ /S₀ conical intersection going from water to propanol. The longer lifetimes might make ^5N7C G more photochemical active to adjacent nucleobases than guanine or other nucleobases within DNA.

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.

Quantum mechanics/molecular mechanics studies on the mechanistic photophysics of sunscreen oxybenzone in methanol solution

Herein, we have employed the QM(CASPT2//CASSCF)/MM method to explore the photophysical and photochemical mechanism of oxybenzone (OB) in methanol solution. Based on the optimized minima, conical intersections and crossing points, and minimum-energy reaction paths related to excited-state intramolecular proton transfer (ESIPT) and excited-state decay paths in the ¹ππ*, ¹nπ*, ³ππ*, ³nπ*, and S₀ states, we have identified several feasible excited-state relaxation pathways for the initially populated S₂(¹ππ*) state to decay to the initial enol isomer' S₀ state. The major one is the singlet-mediated and stretch-torsion coupled ESIPT pathway, in which the system first undergoes an essentially barrierless ¹ππ* ESIPT process to generate the ¹ππ* keto species, and finally realizes its ground state recovery through the subsequent carbonyl stretch-torsion facilitating S₁ → S₀ internal conversion (IC) and the reverse ground-state intramolecular proton transfer (GSIPT) process. The minor ones are related to intersystem crossing (ISC) processes. At the S₂(¹ππ*) minimum, an S₂(¹ππ*)/S₁(¹nπ*)/T₂(³nπ*) three-state intersection region helps the S₂ system branch into the T₁ state through a S₂ → S₁ → T₁ or S₂ → T₂ → T₁ process. Once it has reached the T₁ state, the system may relax to the S₀ state via direct ISC or via subsequent nearly barrierless ³ππ* ESIPT to yield the T₁ keto tautomer and ISC. The resultant S₀ keto species significantly undergoes reverse GSIPT and only a small fraction yields the trans-keto form that relaxes back more slowly. However, due to small spin-orbit couplings at T₁/S₀ crossing points, the ISC to S₀ state occurs very slowly. The present work rationalizes not only the ultrafast excited-state decay dynamics of OB but also its phosphorescence emission at low temperature.

Theoretical Studies on the Excited-State Decay Mechanism of Homomenthyl Salicylate in a Gas Phase and an Acetonitrile Solution

Here, we employ the CASPT2//CASSCF and QM(CASPT2//CASSCF)/MM approaches to explore the photochemical mechanism of homomenthyl salicylate (HMS) in vacuum and an acetonitrile solution. The results show that in both cases, the excited-state relaxation mainly involves a spectroscopically "bright" S₁(¹ππ*) state and the lower-lying T₁ and T₂ states. In the major relaxation pathway, the photoexcited S₁ keto system first undergoes an essentially barrierless excited-state intramolecular proton transfer (ESIPT) to generate the S₁ enol minimum, near which a favorable S₁/S₀ conical intersection decays the system to the S₀ state followed by a reverse ground-state intramolecular proton transfer (GSIPT) to repopulate the initial S₀ keto species. In the minor one, an S₁/T₂/T₁ three-state intersection in the keto region makes the T₁ state populated via direct and T₂-mediated intersystem crossing (ISC) processes. In the T₁ state, an ESIPT occurs, which is followed by ISC near a T₁/S₀ crossing point in the enol region to the S₀ state and finally back to the S₀ keto species. In addition, a T₁/S₀ crossing point near the T₁ keto minimum can also help the system decay to the S₀ keto species. However, small spin-orbit couplings between T₁ and S₀ at these T₁/S₀ crossing points make ISC to the S₀ state very slow and make the system trapped in the T₁ state for a while. The present work rationalizes not only the ultrafast excited-state decay dynamics of HMS but also its low quantum yield of phosphorescence at 77 K.

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.

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.

Photophysical Deactivation Mechanisms of the Pyrimidine Analogue 1-Cyclohexyluracil

The photophysical relaxation mechanisms of 1-cyclohexyluracil, in vacuum and water, were investigated by employing the Multi-State CASPT2 (MS-CASPT2, Multi-State Complete Active-Space Second-Order Perturbation Theory) quantum chemical method and Dunning's cc-pVDZ basis sets. In both environments, our results suggest that the primary photophysical event is the population of the S11(ππ*) bright state. Afterwards, two likely deactivation pathways can take place, which is sustained by linear interpolation in internal coordinates defined via Z-Matrix scans connecting the most important characteristic points. The first one (Route 1) is the same relaxation mechanism observed for uracil, its canonical analogue, i.e., internal conversion to the ground state through an ethylenic-like conical intersection. The other route (Route 2) is the direct population transfer from the S11(ππ*) bright state to the T23(nπ*) triplet state via an intersystem crossing process involving the (S11(ππ*)/T23(nπ*))STCP singlet-triplet crossing point. As the spin-orbit coupling is not too large in either environment, we propose that most of the electronic population initially on the S11(ππ*) state returns to the ground following the same ultrafast deactivation mechanism observed in uracil (Route 1), while a smaller percentage goes to the triplet manifold. The presence of a minimum on the S11(ππ*) potential energy hypersurface in water can help to understand why experimentally it is noticed suppression of the triplet states population in polar protic solvent.

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.

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.