Curious Case of 2-Selenouracil: Efficient Population of Triplet States and Yet Photostable

Thio and Seleno Derivatives of Angelicin as Efficient Triplet Harvesting Photosensitizers: Implications in Photodynamic Therapy

Photodynamic therapy (PDT) is widely accepted in medical practice for its targeted induction of apoptosis in cancerous cells. Angelicin (Ang) has traditionally been known for its efficacy in cancer treatment and its capability to enter a photoexcited triplet state. This study has comprehensively assessed the effects of substituting individual chalcogen atoms at three specific positions in Angelicin, with the objective of facilitating access to this elusive triplet state to enhance its role as a photosensitizer in PDT. The study scrutinizes various enhancements and factors that are crucial for efficient triplet harvesting. The decrease in singlet-triplet energy gap (ΔEST) and increased spin-orbit coupling (SOC) values present numerous viable pathways for intersystem crossing (ISC), leading to the triplet manifold. The lifetime of ISC, thus, decreases from 10-5 s-1 in Ang to 10-8 s-1 in thioangelicin (TAng) and finally to 10-9 s-1 in selenoangelicin (SeAng). Additionally, this study investigates the two-photon absorption properties of thio and seleno-substituted Angelicin for their potentialities as non-UV photosensitizers. The interplay between electron-withdrawing and electron-donating substitutions in these derivatives significantly enhances the two-photon absorption cross-sections (σ) to as high as 49.3 GM while shifting the absorption wavelengths towards the infrared region enabling them as efficient PDT photosensitizers.

Thio and Seleno-Psoralens as Efficient Triplet Harvesting Photosensitizers for Photodynamic Therapy

The Psoralen (Pso) molecule finds extensive applications in photo-chemotherapy, courtesy of its triplet state forming ability. Sulfur and selenium replacement of exocyclic carbonyl oxygen of organic chromophores foster efficient triplet harvesting with near unity triplet quantum yield. These triplet-forming photosensitizers are useful in Photodynamic Therapy (PDT) applications for selective apoptosis of cancer cells. In this work, we have critically assessed the effect of the sulfur and selenium substitution at the exocyclic carbonyl (TPso and SePso, respectively) and endocyclic oxygen positions of Psoralen. It resulted in a significant redshifted absorption spectrum to access the PDT therapeutic window with increased oscillator strength. The reduction in singlet-triplet energy gap and enhancement in the spin-orbit coupling values increase the number of intersystem crossing (ISC) pathways to the triplet manifold, which shortens the ISC lifetime from 10-5 s for Pso to 10-8 s for TPso and 10-9 s for SePso. The intramolecular photo-induced electron transfer process, a competitive pathway to ISC, is also considerably curbed by exocyclic functionalizations. In addition, a maximum of 115 GM of two-photon absorption (2PA) with IR absorption (660-1050 nm) confirms that the Psoralen skeleton can be effectively tweaked via single chalcogen atom replacement to design a suitable PDT photosensitizer.

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.

Vibrational Coherences of the Photoinduced Mixed-Valent Creutz-Taube Ion Revealed by Excited State Dynamics

A recent study of photoinduced mixed-valency in the one-electron reduced form (μ-pz)[RuII(NH3)5]24+ of the Creutz-Taube ion used transient absorption spectroscopy with vis-NIR broadband detection to uncover a mixed-valent excited state with a typical intervalence charge transfer band and a nanosecond lifetime [Pieslinger et al. Angew. Chem., Int. Ed. 2022, 61, e202211747]. Herein, we use excited state dynamics simulations with implicit solvation to elucidate the electronic and vibrational evolution in the first 10 ps after the optical excitation. A manifold of excited states with weak interaction between the metal centers is populated already at time zero due to the breakdown of the Condon approximation and dominates the population of electronic states at short time scales (<0.5 ps). A long-lived vibrational wave packet mostly confined to oscillations of the metal center-bridge distances is observed. The oscillations are traced to the electronic structure properties of states with weak metal-metal coupling. The long-lived mixed-valent excited state of the Creutz-Taube ion analogue is formed vibrationally cold and has a more compact geometry. While experimentally, intersystem crossing and vibrational relaxation were deduced to be completed within 1 ps, our analysis indicates that both processes might persist at longer times.

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.

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.

Nonadiabatic Molecular Dynamics by Multiconfiguration Pair-Density Functional Theory

We present the first implementation of multiconfiguration pair-density functional theory (MC-PDFT) ab initio molecular dynamics. MC-PDFT is a multireference electronic structure method that in many cases has a similar accuracy (or even better accuracy) the complete active space second-order perturbation theory (CASPT2) at a significantly lower computational cost. In this study, we introduced MC-PDFT analytical gradients into the SHARC molecular dynamics program for ab initio, nonadiabatic molecular dynamics simulations. We verify our implementation by examining the intersystem crossing dynamics of thioformaldehyde, and we observe excellent agreement with recent CASPT2 and experimental findings. Moreover, with MC-PDFT, we could perform dynamics simulations with the 12 electron in 10 orbitals active space that was computationally too expensive for direct dynamics with CASPT2.

One order of magnitude increase of triplet state lifetime observed in deprotonated form selenium substituted uracil

Selenium nucleic acids possess unique properties and have been demonstrated to have a wide range of applications such as in DNA X-ray crystallography and novel medical therapies. However, as a heavy atom, selenium substitution may easily alter the photophysical properties of a nucleic acid by red-shifting the absorption spectra and introducing effective intersystem crossing to triplet excited states. In present work, the excited state dynamics of a naturally occurring selenium substituted uracil (2-selenuracil, 2SeU) is studied by using femtosecond transient absorption spectroscopy as well as quantum chemistry calculations. Ultrafast intersystem crossing to the lowest triplet state (T₁) and effective non-radiative decay of this state to the ground state (S₀) are demonstrated in the neutral form 2SeU. However, the triplet lifetime of the deprotonated form 2SeU is found to be almost one order of magnitude longer than that in the neutral one. Quantum chemistry calculations indicate that the short triplet lifetime in 2SeU is due to excited state population decay through a crossing point between T₁ and S₀. In the deprotonated form, shortening the N1-C2 bond length makes the structural distortion more difficult and brings a larger energy barrier on the pathway to the T₁/S₀ crossing point, resulting in one order of magnitude increase of the triplet state lifetime. Our study reveals one key factor to regulate the triplet lifetime of 2SeU and sets the stage to further investigate the photophysical and photochemical properties of 2SeU-containing DNA/RNA duplexes.

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.

Surface hopping dynamics reveal ultrafast triplet generation promoted by S1-T2-T1 spin-vibronic coupling in 2-mercaptobenzothiazole

We investigate the T₁ formation upon populating the optically "bright" S₂ in 2-mercaptobenzothiazole to interpret the underlying relaxation pathways associated with the experimental decay constants reported by D. Koyama and A. J. Orr-Ewing, Phys. Chem. Chem. Phys., 2016, 18, 26224-26235. Energetics, electronic populations and geometries of various stationary points of low-lying electronic states are computed using the semi-classical ab initio surface hopping dynamics simulations. Estimated decay constants of S₂-S₁ internal conversion (IC) and S₁-T₂ intersystem crossing (ISC) are in excellent agreement with the experiment. The observed ultrafast ISC is analyzed based on the S₁-T₂-T₁ spin-vibronic coupling mechanism. In contrast to the previous assignment of 6 ps to the T₂-T₁ IC, our findings enable us to attribute this decay constant to the combined events of T₂-T₁ IC followed by relaxation of vibrationally hot T₁.

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.

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

The excited-state properties and relaxation mechanisms after light irradiation of 6-selenoguanine (6SeG) in water and in DNA have been investigated using a quantum mechanics/molecular mechanics (QM/MM) approach with the multistate complete active space second-order perturbation theory (MS-CASPT2) method. In both environments, the S₁ ¹(n_Seπ₅^*) and S₂ ¹(π_Seπ₅^*) states are predicted to be the spectroscopically dark and bright states, respectively. Two triplet states, T₁ ³(π_Seπ₅^*) and T₂ ³(n_Seπ₅^*), are found energetically below the S₂ state. Extending the QM region to include the 6SeG-Cyt base pair slightly stabilizes the S₂ state and destabilizes the S₁, due to hydrogen-bonding interactions, but it does not affect the order of the states. The optimized minima, conical intersections, and singlet-triplet crossings are very similar in water and in DNA, so that the same general mechanism is found. Additionally, for each excited state geometry optimization in DNA, three kind of structures ("up", "down", and "central") are optimized which differ from each other by the orientation of the C═Se group with respect to the surrounding guanine and thymine nucleobases. After irradiation to the S₂ state, 6SeG evolves to the S₂ minimum, near to a S₂/S₁ conical intersection that allows for internal conversion to the S₁ state. Linear interpolation in internal coordinates indicate that the "central" orientation is less favorable since extra energy is needed to surmount the high barrier in order to reach the S₂/S₁ conical intersection. From the S₁ state, 6SeG can further decay to the T₁ ³(π_Seπ₅^*) state via intersystem crossing, where it will be trapped due to the existence of a sizable energy barrier between the T₁ minimum and the T₁/S₀ crossing point. Although this general S₂ → T₁ mechanism takes place in both media, the presence of DNA induces a steeper S₂ potential energy surface, that it is expected to accelerate the S₂ → S₁ internal conversion.

Can modified DNA base pairs with chalcogen bonding expand the genetic alphabet? A combined quantum chemical and molecular dynamics simulation study

A comprehensive (DFT and MD) computational study is presented with the goal to design and analyze model chalcogen-bonded modified nucleobase pairs that replace one (i.e., A_XY:T, G:C_XY, G_XY:C) or two (G_XY:C_X'Y', X/X' = S, Se and Y/Y' = F, Cl, Br) Watson-Crick (WC) hydrogen bonds of the canonical A:T or G:C pair with chalcogen bond(s). DFT calculations on 18 base pair combinations that replace one WC hydrogen bond with a chalcogen bond reveal that the bases favorably interact in the gas phase (binding strengths up to -140 kJ mol^-1) and water (up to -85 kJ mol^-1). Although the remaining hydrogen bond(s) exhibits similar characteristics to those in the canonical base pairs, the structural features of the (Y-XO) chalcogen bond(s) change significantly with the identity of X and Y. The 36 doubly-substituted (G_XY:C_X'Y') base pairs have structural deviations from canonical G:C similar to those of the singly-substituted modifications (G:C_XY or G_XY:C). Furthermore, despite the replacement of two strong hydrogen bonds with chalcogen bonds, some G_XY:C_X'Y' pairs possess comparable binding energies (up to -132 kJ mol^-1 in the gas phase and up to -92 kJ mol^-1 in water) to the most stable G:C_XY or G_XY:C pairs, as well as canonical G:C. More importantly, G:C-modified pairs containing X = Se (high polarizability) and Y = F (high electronegativity) are the most stable, with comparable or slightly larger (by up to 13 kJ mol^-1) binding energies than G:C. Further characterization of the chalcogen bonding in all modified base pairs (AIM, NBO and NCI analyses) reveals that the differences in the binding energies of modified base pairs are mainly dictated by the differences in the strengths of their chalcogen bonds. Finally, MD simulations on DNA oligonucleotides containing the most stable chalcogen-bonded base pair from each of the four classifications (A_XY:T, G:C_XY, G_XY:C and G_XY:C_X'Y') reveal that the singly-modified G:C pairs best retain the local helical structure and pairing stability to a greater extent than the modified A:T pair. Overall, our study identifies two (G:C_SeF and G_SeF:C) promising pairs that retain chalcogen bonding in DNA and should be synthesized and further explored in terms of their potential to expand the genetic alphabet.

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.

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.