QM/MM Studies on Photoisomerization Dynamics of Azobenzene Chromophore Tethered to a DNA Duplex: Local Unpaired Nucleobase Plays a Crucial Role

Quantum Mechanics/Molecular Mechanics Studies on the Photophysical Mechanism of Methyl Salicylate

Methyl salicylate (MS) as a subunit of larger salicylates found in commercial sunscreens has been shown to exhibit keto-enol tautomerization and dual fluorescence emission via excited-state intramolecular proton transfer (ESIPT) after the absorption of ultraviolet (UV) radiation. However, its excited-state relaxation mechanism is unclear. Herein, we have employed the quantum mechanics(CASPT2//CASSCF)/molecular mechanics method to explore the ESIPT and excited-state relaxation mechanism of MS in the lowest three electronic states, that is, S₀, S₁, and T₁ states, in a methanol solution. Based on the optimized geometric and electronic structures, conical intersections and crossing points, and minimum-energy paths combined with the computed linearly interpolated Cartesian coordinate paths, the photophysical mechanism of MS has been proposed. The S₁ state is a spectroscopically bright ¹ππ* state in the Franck-Condon region. From the initially populated S₁ state, there exist three nonradiative relaxation paths to repopulate the S₀ state. In the first one, the S₁ system (i.e., ketoB form) first undergoes an ESIPT path to generate an S₁ tautomer (i.e., enol form) that exhibits a large Stokes shift in experiments. The generated S₁ enol tautomer further evolves toward the nearby S₁/S₀ conical intersection and then hops to the S₀ state, followed by the backward ground-state intramolecular proton transfer (GSIPT) to the initial ketoB form S₀ state. In the second one, the S₁ system first hops through the S₁ → T₁ intersystem crossing (ISC) to the T₁ state, which then further decays to the S₀ state via T₁ → S₀ ISC at the T₁/S₀ crossing point. In the third path, the T₁ system that stems from the S₁ → T₁ ISC process via the S₁/T₁ crossing point first takes place a T₁ ESIPT to generate a T₁ enol tautomer, which can further decay to the S₀ state via T₁-to-S₀ ISC. Finally, the GSIPT occurs to back the system to the initial ketoB form S₀ state. Our present work could contribute to understanding the photophysics of MS and its derivatives.

Theoretical insights into the effect of size and substitution patterns of azobenzene derivatives on the DNA G-quadruplex

Introducing photoswitches into the DNA G-quadruplex provides excellent opportunities to control folding and unfolding of these assemblies, demonstrating their potential in the development of novel nanodevices with medical and nanotechnology applications. Using a quantum mechanics/molecular mechanics (QM/MM) scheme, we carried out a series of simulations to identify the effect of the size and substitution patterns of three azobenzene derivatives (AZ1, AZ2 and AZ3) on the excitation energies of the two lowest excited states of the smallest photoswitchable G-quadruplex reported to date. We demonstrated that the size and the substitution pattern do not affect the ultrafast cis-trans photoiomerization mechanism of the azobenzene derivatives significantly, in agreement with the experiment. However, molecular dynamics simulations revealed that while AZ2 and AZ3 G-quadruplexes are structurally stable during the simulations, the AZ1 G-quadruplex undergoes larger structural changes and shows two ground state populations that differ in the azobenzene backbone adopting two different conformations. AZ1, with para-para substitution pattern, provides more flexibility to the whole G-quadruplex structure compared to AZ2 and AZ3, and can thus facilitate the photoisomerization reaction between a nonpolymorphic, stacked, tetramolecular G-quadruplex and an unstructured state after trans-cis isomerization occurring in a longer time dynamics, in agreement with the experimental findings. The QM/MM simulations of the absorption spectra indicated that the thermal fluctuation plays a more crucial role in the main absorption band of the azobenzene derivatives than the inclusion of the G-quadruplex, implying that the influence of the G-quadruplex environment is minimal. We propose that the latter is attributed to the position of the azobenzene linkers in the G-quadruplexes, i.e. the edgewise loops containing the azobenzene moieties that are located above the G-quartets, not being fully embedded inside or involved in the stacked structure. Our theoretical findings provide support to a recent study of the photoresponsive formation of photoswitchable G-quadruplex motifs.

Ultrafast Nonadiabatic Photoisomerization Dynamics Mechanism for the UV Photoprotection of Stilbenoids in Grape Skin

Natural UV photoprotection plays a vital role in physiological protection. It has been reported that UVC radiation can make resveratrol (RSV) and piceatannol (PIC) accumulate in grape skin. In this work, we demonstrated that RSV and PIC could significantly absorb UVA and UVB, and confirmed their satisfactory photostability. Furthermore, we clarified the UV photoprotection mechanism of typical stilbenoids of RSV and PIC for the first time by using combined femtosecond transient absorption (FTA) spectroscopy and time-dependent density functional theory (TD-DFT) calculations. RSV and PIC can be photoexcited to the excited state after UVA and UVB absorption. Subsequently, the photoisomerized RSV and PIC quickly relax to the ground state via nonadiabatic transition from the S₁ state at a conical intersection (CI) position between potential energy surfaces (PESs) of S₁ and S₀ states. This ultrafast trans-cis photoisomerization will take place within a few tens of picoseconds. As a result, the UV energy absorbed by RSV and PIC could be dissipated by an ultrafast nonadiabatic photoisomerization process.

QM/MM Studies on the Photophysical Mechanism of a Truncated Octocrylene Model

Methyl 2-cyano-3,3-diphenylacrylate (MCDPA) shares the same molecular skeleton with octocrylene (OCR) that is one of the most common molecules used in commercially available sunscreens. However, its excited-state relaxation mechanism is unclear. Herein, we have used the QM(CASPT2//CASSCF)/MM method to explore spectroscopic properties, geometric and electronic structures, relevant conical intersections and crossing points, and excited-state relaxation paths of MCDPA in methanol solution. We found that in the Franck-Condon (FC) region, the V(¹ππ*) state is energetically lower than the V'(¹ππ*) state only by 2.8 kcal/mol and is assigned to experimentally observed maximum absorption band. From these two initially populated singlet states, there exist three nonradiative relaxation paths to repopulate the S₀ state. In the first one, when the V(¹ππ*) state is populated in the FC region, the system diabatically evolves along the V(¹ππ*) state into its minimum where the internal conversion to S₀ occurs. In the second one, the V'(¹ππ*) state is populated in the FC region and the system adiabatically overcomes a barrier of ca. 3.0 kcal/mol to approach the V(¹ππ*) minimum eventually leading to a V(¹ππ*)-to-S₀ internal conversion. In the third one, the V'(¹ππ*) state first hops via the intersystem crossing to the T₂ state, which then decays through the internal conversion to the T₁ state. The T₁ state is finally converted to the S₀ state via the T₁/S₀ crossing point. Our present work contributes to understanding the photophysics of OCR and its variants.