Quantum mechanical/molecular mechanical trajectory surface hopping molecular dynamics simulation by spin-flip time-dependent density functional theory

Temperature Controlled Decay and Pendulum Dynamics of Green Fluorescent Protein (GFP) Chromophore

The excited-state dynamics of the GFP chromophore, HBDI- (anionic p-hydroxybenzylidene-2,3-dimethylimidazolinone), were investigated through a combination of theoretical nonadiabatic molecular dynamics (NAMD) simulations and femtosecond transient absorption spectroscopy (fs-TA). The NAMD simulations revealed that the primary dynamics in excited states involve the formation of a P-twisted intermediate (S1min,P), which undergoes pendulum-like oscillations with respect to ϕ = 90°. This motion serves as a reservoir for the excited-state population and the primary source of fluorescence. Rather than a direct channel from the major S1min,P, a coordinated pathway of S1min,P → S1min → S1min,I → S0 is responsible for the decay to the ground state, emphasizing the importance of planar intermediate (S1min) formation. The experimental fs-TA spectra confirmed these dynamics, revealing three distinct time scales (340-470 fs, 1.4 ps, and 8.3 ps), corresponding to the formation of S1min,P and its decay governed by the coordinated pathway. At low temperatures, the coordinated decay pathway is suppressed, leading to prolonged fluorescence lifetimes, consistent with low-temperature experimental results. This study presents a new model for the excited-state dynamics of GFP chromophore, suggesting that pendulum motion and the coordinated decay pathway play a crucial role in regulating fluorescence intensity.

A density functional study of the photocatalytic degradation of polycaprolactone by the decatungstate anion in acetonitrile solution

A recent experimental study has reported that decatungstate [W10O32]4- can degrade various polyesters in the presence of light and molecular oxygen [Li et al., Nanoscale, 2023, 15, 15038]. We apply density functional theory to the photocatalyst-polycaprolactone model complex in acetonitrile solution and elucidate the degradation mechanisms and catalytic cycle. We consider hydrogen atom transfer (HAT) and single electron transfer (SET) mechanisms. The potential energy profiles show that the former proceeds exergonically in a single step but that the latter involves a subsequent proton transfer and finally yields HAT products as well. Oxygenated polymer species can regain the transferred hydrogen and regenerate the reduced photocatalyst. We propose a photocatalytic cycle that realizes both the photocatalyst regeneration and the polymer degradation.

Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges

Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin-orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm-Dancoff approximation to TDDFT, spin-flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine-Ag20 complex), we emphasize that caution is needed to properly account for system-environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.

Temperature-Dependent Fluorescence of mPlum Fluorescent Protein from 295 to 20 K

The development of bright fluorescent proteins (FPs) emitting beyond 600 nm continues to be of interest both from a fundamental perspective in understanding protein-chromophore interactions and from a practical perspective as these FPs would be valuable for cellular imaging. We previously reported ultrafast spectral observations of the excited-state dynamics in mPlum resulting from interconversion between direct hydrogen bonding and water-mediated hydrogen bonding between the chromophore acylimine carbonyl and the Glu16 side chain. Here, we report temperature-dependent steady-state and time-resolved fluorescence measurements of mPlum and its E16H variant, which does not contain a side-chain permitting hydrogen bonding with the acylimine carbonyl. Lowering the temperature of the system freezes interconversion between the hydrogen-bonding states, thus revealing the spectral signatures of the two states. Analysis of the temperature-dependent spectra assuming Boltzmann populations of the two states yields a 205 cm^-1 energy difference. This value agrees with the predictions from a quantum mechanics/molecular mechanics study of mPlum (198 cm^-1). This study demonstrates the first use of cryogenic spectroscopy to quantify the energetics and timescales of FP chromophore structural states that were only previously obtained from computational methods and further confirms the importance of acylimine hydrogen-bonding dynamics to the fluorescence spectral shifts of red FPs.

Fluorescence Modulation by Ultrafast Chromophore Twisting Events: Developing a Powerful Toolset for Fluorescent-Protein-Based Imaging

The advancement of modern life sciences has benefited tremendously from the discovery and development of fluorescent proteins (FPs), widely expressed in live cells to track a myriad of cellular events. The chromophores of various FPs can undergo many ultrafast photophysical and/or photochemical processes in the electronic excited state and emit fluorescence with different colors. However, the chromophore becomes essentially nonfluorescent in solution environment due to its intrinsic twisting capability upon photoexcitation. To study "microscopic" torsional events and their effects on "macroscopic" fluorescence, we have developed an integrated ultrafast characterization platform involving femtosecond transient absorption (fs-TA) and wavelength-tunable femtosecond stimulated Raman spectroscopy (FSRS). A wide range of naturally occurring, circularly permuted, non-canonical amino-acid-decorated FPs and FP-based optical highlighters with photochromicity, photoconversion, and/or photoswitching capabilities have been recently investigated in great detail. Twisting conformational motions were elucidated to exist in all of these systems but to various extents. The associated different ultrafast pathways can be monitored via frequency changes of characteristic Raman bands during primary events and functional processes. The mapped electronic and structural dynamics information is crucial and has shown great potential and initial success for the rational design of proteins and other photoreceptors with novel functions and fluorescence properties.

Nonadiabatic dynamics with spin-flip vs linear-response time-dependent density functional theory: A case study for the protonated Schiff base C5H6NH2. J Chem Phys 2021;155:124111. [PMID: 34598550 DOI: 10.1063/5.0062757]

Nonadiabatic trajectory surface hopping simulations are reported for trans-C₅H₆NH₂ ⁺, a model of the rhodopsin chromophore, using the augmented fewest-switches algorithm. Electronic structure calculations were performed using time-dependent density functional theory (TDDFT) in both its conventional linear-response (LR) and its spin-flip (SF) formulations. In the SF-TDDFT case, spin contamination in the low-lying singlet states is removed by projecting out the lowest triplet component during iterative solution of the TDDFT eigenvalue problem. The results show that SF-TDDFT qualitatively describes the photoisomerization of trans-C₅H₆NH₂ ⁺, with favorable comparison to previous studies using multireference electronic structure methods. In contrast, conventional LR-TDDFT affords qualitatively different photodynamics due to an incorrect excited-state potential surface near the Franck-Condon region. In addition, the photochemistry (involving pre-twisting of the central double bond) appears to be different for SF- and LR-TDDFT, which may be a consequence of different conical intersection topographies afforded by these two methods. The present results contrast with previous surface-hopping studies suggesting that the LR-TDDFT method's incorrect topology around S₁/S₀ conical intersections is immaterial to the photodynamics.

Resolving the ultrafast dynamics of the anionic green fluorescent protein chromophore in water

The chromophore of the green fluorescent protein (GFP) is critical for probing environmental influences on fluorescent protein behavior. Using the aqueous system as a bridge between the unconfined vacuum system and a constricting protein scaffold, we investigate the steric and electronic effects of the environment on the photodynamical behavior of the chromophore. Specifically, we apply ab initio multiple spawning to simulate five picoseconds of nonadiabatic dynamics after photoexcitation, resolving the excited-state pathways responsible for internal conversion in the aqueous chromophore. We identify an ultrafast pathway that proceeds through a short-lived (sub-picosecond) imidazolinone-twisted (I-twisted) species and a slower (several picoseconds) channel that proceeds through a long-lived phenolate-twisted (P-twisted) intermediate. The molecule navigates the non-equilibrium energy landscape via an aborted hula-twist-like motion toward the one-bond-flip dominated conical intersection seams, as opposed to following the pure one-bond-flip paths proposed by the excited-state equilibrium picture. We interpret our simulations in the context of time-resolved fluorescence experiments, which use short- and long-time components to describe the fluorescence decay of the aqueous GFP chromophore. Our results suggest that the longer time component is caused by an energetically uphill approach to the P-twisted intersection seam rather than an excited-state barrier to reach the twisted intramolecular charge-transfer species. Irrespective of the location of the nonadiabatic population events, the twisted intersection seams are inefficient at facilitating isomerization in aqueous solution. The disordered and homogeneous nature of the aqueous solvent environment facilitates non-selective stabilization with respect to I- and P-twisted species, offering an important foundation for understanding the consequences of selective stabilization in heterogeneous and rigid protein environments.

Trajectory Surface Hopping Approach to Condensed-Phase Nonradiative Relaxation Dynamics Using Divide-and-Conquer Spin-Flip Time-Dependent Density-Functional Tight Binding

Nonradiative relaxation of excited molecules is central to many crucial issues in photochemistry. Condensed phases are typical contexts in which such problems are considered, and the nonradiative relaxation dynamics are expected to be significantly affected by interactions with the environment, for example, a solvent. We developed a nonadiabatic molecular dynamics simulation technique that can treat the nonradiative relaxation and explicitly include the environment in the calculations without a heavy computational burden. Specifically, we combined trajectory surface hopping with Tully's fewest-switches algorithm, a tight-binding approximated version of spin-flip time-dependent density-functional theory, and divide-and-conquer (DC) spatial fragmentation scheme. Numerical results showed that this method can treat systems with thousands of atoms within reasonable computational resources, and the error arising from DC fragmentation is negligibly small. Using this method, we obtained molecular insights into the solvent dependence of the photoexcited-state dynamics of trans-azobenzene, which demonstrate the importance of the environment for condensed-phase nonradiative relaxation.

Fast Nonadiabatic Molecular Dynamics via Spin-Flip Time-Dependent Density-Functional Tight-Binding Approach: Application to Nonradiative Relaxation of Tetraphenylethylene with Locked Aromatic Rings

Nonadiabatic dynamics around conical intersections between ground and excited states are crucial to understand excited-state phenomena in complex chemical systems. With this background in mind, we present an approach combining fewest-switches trajectory surface hopping and spin-flip (SF) time-dependent (TD) density-functional tight binding (DFTB), which is a simplified version of SF-TD density functional theory (DFT) with semiempirical parametrizations, for computationally efficient nonadiabatic molecular dynamics simulations. The estimated computational time of the SF-TD-DFTB approach is several orders of magnitude lower than that of SF-TD-DFT. In addition, the proposed method reproduces the time scales and quantum yields in photoisomerization reactions of azobenzene at a level comparable with conventional ab initio approaches, demonstrating reasonable accuracy. Finally, we report a practical application of the developed technique to explore the nonradiative relaxation processes of tetraphenylethylene and its derivative with torsionally locked aromatic rings and discuss the effect of locking the rings on the excited-state lifetime.