Interplay between Conjugation and Size-Driven Delocalization Leads to Characteristic Properties of Substituted Thymines

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.

Quantum effects in photosensitization: the case of singlet oxygen generation by thiothymines

Photosensitization is the indirect electronic excitation of a molecule with the aid of a photosensitizer and is a bimolecular nonradiative energy transfer. In this study, we have attempted to elucidate its mechanism, and we do this by calculating rate constants of photosensitization of oxygen by thiothymines (2-thiothymine, 4-thiothymine and 2,4 dithiothymine). The rate constants have been calculated using two approaches: (a) a classical limit of Fermi's Golden Rule (FGR), and (b) a time-dependent variant of FGR, where the treatment is purely quantum mechanical. The former approach has previously been developed for bimolecular systems and has been applied to the photosensitization reactions studied here. The latter approach, however, has so far only been used for unimolecular reactions, and in this work, we describe how it can be adapted for bimolecular reactions. Experimentally, all three thiothymines are known to have significant singlet oxygen yields, which are indicative of similar rates. Rate constants calculated using the time-dependent variant of FGR are similar across all three thiothymines. While the classical approximation gives reasonable rate constants for 2-thiothymine, it severely underestimates them for 4-thiothymine and 2,4 dithiothymine, by several orders of magnitude. This work indicates the importance of quantum effects in driving photosensitization.

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.

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.

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.

Triplet Decay Dynamics in Sulfur-Substituted Thymine: How Position of Substitution Matters

Sulfur-substituted analogues of thymine are of three types depending on the position of sulfur substitution: 2-thiothymine (2tThy), 4-thothymine (4tThy), and 2,4-dithiothymine (dtThy). These molecules, on photoexcitation, are known to form in their triplet state with near unity yield. Consequently, they are able to photosensitize ground state molecular oxygen to singlet oxygen, a property which makes them potential drugs for photodynamic therapy (PDT). The singlet oxygen yield is directly correlated with the triplet lifetime of the thiothymine, which in turn is governed by its triplet decay dynamics. In this work, the dependence of the triplet decay dynamics on the position of sulfur substitution is investigated by comparatively studying all three thiothymines. The topology of the triplet potential energy surface and decay mechanism of 2tThy is found to be distinctly different from 4tThy and dtThy. The fundamental reason for this is the different electronic natures of the two C═X (X = O, S) moieties in each molecule, one of which is conjugated with a C═C bond, while the other is not. Further, it is shown that the triplet lifetime of 2tThy can be increased by manipulating the energetic ordering of its molecular orbitals with unobtrusive substitutions, thus making it a better candidate for a PDT drug.

Photochemical and Photodynamical Properties of Sulfur-Substituted Nucleic Acid Bases

Sulfur-substituted nucleobases (a.k.a., thiobases) are among the world's leading prescriptions for chemotherapy and immunosuppression. Long-term treatment with azathioprine, 6-mercaptopurine and 6-thioguanine has been correlated with the photoinduced formation of carcinomas. Establishing an in-depth understanding of the photochemical properties of these prodrugs may provide a route to overcoming these carcinogenic side effects, or, alternatively, a basis for developing effective compounds for targeted phototherapy. In this review, a broad examination is undertaken, surveying the basic photochemical properties and excited-state dynamics of sulfur-substituted analogs of the canonical DNA and RNA nucleobases. A molecular-level understanding of how sulfur substitution so remarkably perturbs the photochemical properties of the nucleobases is presented by combining experimental results with quantum-chemical calculations. Structure-property relationships demonstrate the impact of site-specific sulfur substitution on the photochemical properties, particularly on the population of the reactive triplet state. The value of fundamental photochemical investigations for driving the development of ultraviolet-A chemotherapeutics is showcased. The most promising photodynamic agents identified thus far have been investigated in various carcinoma cell lines and shown to decrease cell proliferation upon exposure to ultraviolet-A radiation. Overarching principles have been elucidated for the impact that sulfur substitution of the carbonyl oxygen has on the photochemical properties of the nucleobases.