Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water

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.

Synthesis of 6-Selanyl-2-triazolylpurine Derivatives Using 2,6-Bistriazolylpurines as Starting Materials

Two pathways toward 6-selanyl-2-triazolylpurine derivatives were designed. The first method involved the synthesis of 2-chloro-6-selanylpurine derivatives, further SNAr reaction with NaN3, and following CuAAC using different alkynes. The second method was based on the synthesis of 2,6-bistriazolylpurine derivatives as starting materials followed by SNAr reaction with commercial or in situ generated selenols as nucleophiles. A series of 2-chloro-6-selanylpurine derivatives were obtained in yields up to 84%. It was found that in the latter compounds, 6-selanyl moiety was the better leaving group compared to 2-chlorosubstituent in SNAr reactions. On the other hand, the SNAr reaction between 2,6-bistriazolylpurines and selenols or diselenides was successful, and 13 examples of 6-selanyl-2-triazolylpurine derivatives were obtained in yields up to 87%. This direct approach for the Se-C bond formation proved the ability of the 1,2,3-triazolyl ring at the C6 position of purine to act as a good leaving group.

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.

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.

Role of Ultrafast Internal Conversion and Intersystem Crossing in the Nonadiabatic Relaxation Dynamics of ortho-Nitrobenzaldehyde

ortho-Nitrobenzaldehyde (oNBA) is a well-known photoactivated acid and a prototypical photolabile nitro-aromatic compound. Despite extensive investigations, the ultrafast relaxation dynamics of oNBA is still not properly understood, especially concerning the role of the triplet states. In this work, we provide an in-depth picture of this dynamics by combining single- and multireference electronic structure methods with potential energy surface exploration and nonadiabatic dynamics simulations using the Surface Hopping including ARbitary Couplings (SHARC) approach. Our results reveal that the initial decay from the bright ππ* state to the S₁ minimum is barrierless. It involves three changes in electronic structure from ππ* (ring) to nπ* (nitro group), to nπ* (aldehyde group), and then to another nπ* (nitro group). The decay of the ππ* takes 60-80 fs and can be tracked with time-resolved luminescence spectroscopy, where we predict for the first time a short-lived coherence of the luminescence energy with a 25 fs period. Intersystem crossing can occur already during the S₄ → S₁ deactivation cascade but also from S₁, with a time constant of about 2.4 ps and such that first a triplet ππ* state localized on the nitro group is populated. The triplet population first evolves into an nπ* and then quickly undergoes hydrogen transfer to form a biradical intermediate, from where the ketene is eventually produced. The majority of the excited population decays from S₁ through two conical intersections of equal utilization, a previously unreported one involving a scissoring motion of the nitro group that leads back to the oNBA ground state and the one involving hydrogen transfer that leads to the ketene intermediate.

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.