Quantum Mechanics/Molecular Mechanics Study on the Photoreactions of Dark- and Light-Adapted States of a Blue-Light YtvA LOV Photoreceptor

Deciphering Photoreceptors Through Atomistic Modeling from Light Absorption to Conformational Response

In this review, we discuss the successes and challenges of the atomistic modeling of photoreceptors. Throughout our presentation, we integrate explanations of the primary methodological approaches, ranging from quantum mechanical descriptions to classical enhanced sampling methods, all while providing illustrative examples of their practical application to specific systems. To enhance the effectiveness of our analysis, our primary focus has been directed towards the examination of applications across three distinct photoreceptors. These include an example of Blue Light-Using Flavin (BLUF) domains, a bacteriophytochrome, and the orange carotenoid protein (OCP) employed by cyanobacteria for photoprotection. Particular emphasis will be placed on the pivotal role played by the protein matrix in fine-tuning the initial photochemical event within the embedded chromophore. Furthermore, we will investigate how this localized perturbation initiates a cascade of events propagating from the binding pocket throughout the entire protein structure, thanks to the intricate network of interactions between the chromophore and the protein.

Quantum mechanics/molecular mechanics studies on mechanistic photophysics of cytosine aza-analogues: 2,4-diamino-1,3,5-triazine and 2-amino-1,3,5-triazine in aqueous solution

The excited-state properties and photophysics of cytosine aza-analogues, i.e., 2,4-diamino-1,3,5-triazine (2,4-DT) and 2-amino-1,3,5-triazine (2-AT) in solution have been systematically explored using the QM(MS-CASPT2//CASSCF)/MM approach. The excited-state nonradiative relaxation mechanisms for the initially photoexcited S₁(ππ*) state decay back to the S₀ state are proposed in terms of the present computed minima, surface crossings (conical intersections and singlet-triplet crossings), and excited-state decay paths in the S₁, S₂, T₁, T₂, and S₀ states. Upon photoexcitation to the bright S₁(ππ*) state, 2,4-DT quickly relaxes to its S₁ minimum and then overcomes a small energy barrier of 5.1 kcal mol^-1 to approach a S₁/S₀ conical intersection, where the S₁ system hops to the S₀ state through S₁ → S₀ internal conversion (IC). In addition, at the S₁ minimum, the system could partially undergo intersystem crossing (ISC) to the T₁ state, followed by further ISC to the S₀ state via the T₁/S₀ crossing point. In the T₁ state, an energy barrier of 7.9 kcal mol^-1 will trap 2,4-DT for a while. In parallel, for 2-AT, the system first relaxes to the S₁ minimum and then S₁ → S₀ IC or S₁ → T₁ → S₀ ISCs take place to the S₀ state by surmounting a large barrier of 15.3 kcal mol^-1 or 11.9 kcal mol^-1, respectively, which heavily suppress electronic transition to the S₀ state. Different from 2,4-DT, upon photoexcitation in the Franck-Condon region, 2-AT can quickly evolve in an essentially barrierless manner to nearby S₂/S₁ conical intersection, where the S₂ and T₁ states can be populated. Once it hops to the S₂ state, the system will overcome a relatively small barrier (6.6 kcal mol^-1vs. 15.3 kcal mol^-1) through IC to the S₀ state. Similarly, an energy barrier of 11.9 kcal mol^-1 heavily suppresses the T₁ state transformation to the S₀ state. The present work manifests that the amination/deamination of the triazine rings can affect some degree of different vertical and adiabatic excitation energies and nonradiative decay pathways in solution. It not only rationalizes excited-state decay dynamics of 2,4-DT and 2-AT in aqueous solution but could also provide insights into the understanding of the photophysics of aza-nucleobases.

Quantum mechanics/molecular mechanics studies on the mechanistic photophysics of sunscreen oxybenzone in methanol solution

Herein, we have employed the QM(CASPT2//CASSCF)/MM method to explore the photophysical and photochemical mechanism of oxybenzone (OB) in methanol solution. Based on the optimized minima, conical intersections and crossing points, and minimum-energy reaction paths related to excited-state intramolecular proton transfer (ESIPT) and excited-state decay paths in the ¹ππ*, ¹nπ*, ³ππ*, ³nπ*, and S₀ states, we have identified several feasible excited-state relaxation pathways for the initially populated S₂(¹ππ*) state to decay to the initial enol isomer' S₀ state. The major one is the singlet-mediated and stretch-torsion coupled ESIPT pathway, in which the system first undergoes an essentially barrierless ¹ππ* ESIPT process to generate the ¹ππ* keto species, and finally realizes its ground state recovery through the subsequent carbonyl stretch-torsion facilitating S₁ → S₀ internal conversion (IC) and the reverse ground-state intramolecular proton transfer (GSIPT) process. The minor ones are related to intersystem crossing (ISC) processes. At the S₂(¹ππ*) minimum, an S₂(¹ππ*)/S₁(¹nπ*)/T₂(³nπ*) three-state intersection region helps the S₂ system branch into the T₁ state through a S₂ → S₁ → T₁ or S₂ → T₂ → T₁ process. Once it has reached the T₁ state, the system may relax to the S₀ state via direct ISC or via subsequent nearly barrierless ³ππ* ESIPT to yield the T₁ keto tautomer and ISC. The resultant S₀ keto species significantly undergoes reverse GSIPT and only a small fraction yields the trans-keto form that relaxes back more slowly. However, due to small spin-orbit couplings at T₁/S₀ crossing points, the ISC to S₀ state occurs very slowly. The present work rationalizes not only the ultrafast excited-state decay dynamics of OB but also its phosphorescence emission at low temperature.

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.

Cysteinyl radicals in chemical synthesis and in nature

Nature harnesses the unique properties of cysteinyl radical intermediates for a diverse range of essential biological transformations including DNA biosynthesis and repair, metabolism, and biological photochemistry. In parallel, the synthetic accessibility and redox chemistry of cysteinyl radicals renders them versatile reactive intermediates for use in a vast array of synthetic applications such as lipidation, glycosylation and fluorescent labelling of proteins, peptide macrocyclization and stapling, desulfurisation of peptides and proteins, and development of novel therapeutics. This review provides the reader with an overview of the role of cysteinyl radical intermediates in both chemical synthesis and biological systems, with a critical focus on mechanistic details. Direct insights from biological systems, where applied to chemical synthesis, are highlighted and potential avenues from nature which are yet to be explored synthetically are presented.

Hydrogen-Bond Network Determines the Early Photoisomerization Processes of Cph1 and AnPixJ Phytochromes

Phytochrome proteins are light receptors that play a pivotal role in regulating the life cycles of plants and microorganisms. Intriguingly, while cyanobacterial phytochrome Cph1 and cyanobacteriochrome AnPixJ use the same phycocyanobilin (PCB) chromophore to absorb light, their excited‐state behavior is very different. We employ multiscale calculations to rationalize the different early photoisomerization mechanisms of PCB in Cph1 and AnPixJ. We found that their electronic S₁, T₁, and S₀ potential minima exhibit distinct geometric and electronic structures due to different hydrogen bond networks with the protein environment. These specific interactions influence the S₁ electronic structures along the photoisomerization paths, ultimately leading to internal conversion in Cph1 but intersystem crossing in AnPixJ. This explains why the excited‐state relaxation in AnPixJ is much slower (ca. 100 ns) than in Cph1 (ca. 30 ps). Further, we predict that efficient internal conversion in AnPixJ can be achieved upon protonating the carboxylic group that interacts with PCB.

QM calculations predict the energetics and infrared spectra of transient glutamine isomers in LOV photoreceptors

Photosensory receptors containing the flavin-binding light-oxygen-voltage (LOV) domain are modular proteins that fulfil a variety of biological functions ranging from gene expression to phototropism. The LOV photocycle is initiated by blue-light and involves a cascade of intermediate species, including an electronically excited triplet state, that leads to covalent bond formation between the flavin mononucleotide (FMN) chromophore and a nearby cysteine residue. Subsequent conformational changes in the polypeptide chain arise due to the remodelling of the hydrogen bond network in the cofactor binding pocket, whereby a conserved glutamine residue plays a key role in coupling FMN photochemistry with LOV photobiology. Although the dark-to-light transition of LOV photosensors has been previously addressed by spectroscopy and computational approaches, the mechanistic basis of the underlying reactions is still not well understood. Here we present a detailed computational study of three distinct LOV domains: EL222 from Erythrobacter litoralis, AsLOV2 from the second LOV domain of Avena sativa phototropin 1, and RsLOV from Rhodobacter sphaeroides LOV protein. Extended protein-chromophore models containing all known crucial residues involved in the initial steps (femtosecond-to-microsecond) of the photocycle were employed. Energies and rotational barriers were calculated for possible rotamers and tautomers of the critical glutamine side chain, which allowed us to postulate the most energetically favoured glutamine orientation for each LOV domain along the assumed reaction path. In turn, for each evolving species, infrared difference spectra were constructed and compared to experimental EL222 and AsLOV2 transient infrared spectra, the former from original work presented here and the latter from the literature. The good agreement between theory and experiment permitted the assignment of the majority of observed bands, notably the ∼1635 cm-1 transient of the adduct state to the carbonyl of the glutamine side chain after rotation. Moreover, both the energetic and spectroscopic approaches converge in suggesting a facile glutamine flip at the adduct intermediate for EL222 and more so for AsLOV2, while for RsLOV the glutamine keeps its initial configuration. Additionally, the computed infrared shifts of the glutamine and interacting residues could guide experimental research addressing early events of signal transduction in LOV proteins.

Peripheral Methionine Residues Impact Flavin Photoreduction and Protonation in an Engineered LOV Domain Light Sensor

Proton-coupled electron transfer reactions play critical roles in many aspects of sensory phototransduction. In the case of flavoprotein light sensors, reductive quenching of flavin excited states initiates chemical and conformational changes that ultimately transmit light signals to downstream targets. These reactions generally require neighboring aromatic residues and proton-donating side chains for rapid and coordinated electron and proton transfer to flavin. Although photoreduction of flavoproteins can produce either the anionic (ASQ) or neutral semiquinone (NSQ), the factors that favor one over the other are not well understood. Here we employ a biologically active variant of the light-oxygen-voltage (LOV) domain protein VVD devoid of the adduct-forming Cys residue (VVD-III) to probe the mechanism of flavin photoreduction and protonation. A series of isosteric and conservative residue replacements studied by rate measurements, fluorescence quantum yields, FTIR difference spectroscopy, and molecular dynamics simulations indicate that tyrosine residues facilitate charge recombination reactions that limit sustained flavin reduction, whereas methionine residues facilitate radical propagation and quenching and also gate solvent access for flavin protonation. Replacement of a single surface Met residue with Leu favors formation of the ASQ over the NSQ and desensitizes photoreduction to oxidants. In contrast, increasing site hydrophilicity by Gln substitution promotes rapid NSQ formation and weakens the influence of the redox environment. Overall, the photoreactivity of VVD-III can be understood in terms of redundant electron donors, internal hole quenching, and coupled proton transfer reactions that all depend upon protein conformation, dynamics, and solvent penetration.

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.

Phototropin-Inspired Chemoselective Synthesis of Unsymmetrical Disulfides: Aerobic Oxidative Heterocoupling of Thiols Using Flavin Photocatalysis

Inspired by the photochemical mechanism of a plant blue-light receptor, a unique flavin-based photocatalytic system was developed for the chemoselective heterocoupling of two different thiols, which enabled the facile synthesis of unsymmetrical disulfides. Owing to the redox- and photo-organocatalysis of flavin, the coupling reaction took place under mild metal-free conditions and visible light irradiation with the use of air, which is recognized as the ideal green oxidant.

Unraveling the Mechanism of a LOV Domain Optogenetic Sensor: A Glutamine Lever Induces Unfolding of the Jα Helix

Light-activated protein domains provide a convenient, modular, and genetically encodable sensor for optogenetics and optobiology. Although these domains have now been deployed in numerous systems, the precise mechanism of photoactivation and the accompanying structural dynamics that modulate output domain activity remain to be fully elucidated. In the C-terminal light-oxygen-voltage (LOV) domain of plant phototropins (LOV2), blue light activation leads to formation of an adduct between a conserved Cys residue and the embedded FMN chromophore, rotation of a conserved Gln (Q513), and unfolding of a helix (Jα-helix) which is coupled to the output domain. In the present work, we focus on the allosteric pathways leading to Jα helix unfolding in Avena sativa LOV2 (AsLOV2) using an interdisciplinary approach involving molecular dynamics simulations extending to 7 μs, time-resolved infrared spectroscopy, solution NMR spectroscopy, and in-cell optogenetic experiments. In the dark state, the side chain of N414 is hydrogen bonded to the backbone N-H of Q513. The simulations predict a lever-like motion of Q513 after Cys adduct formation resulting in a loss of the interaction between the side chain of N414 and the backbone C═O of Q513, and formation of a transient hydrogen bond between the Q513 and N414 side chains. The central role of N414 in signal transduction was evaluated by site-directed mutagenesis supporting a direct link between Jα helix unfolding dynamics and the cellular function of the Zdk2-AsLOV2 optogenetic construct. Through this multifaceted approach, we show that Q513 and N414 are critical mediators of protein structural dynamics, linking the ultrafast (sub-ps) excitation of the FMN chromophore to the microsecond conformational changes that result in photoreceptor activation and biological function.

Computer-Assisted Design of Environmentally Friendly and Light-Stable Fluorescent Dyes for Textile Applications

Five potentially environmentally friendly and light-stable hemicyanine dyes were designed based on integrated consideration of photo, environmental, and computational chemistry as well as textile applications. Two of them were synthesized and applied in dyeing polyacrylonitrile (PAN), cotton, and nylon fabrics, and demonstrated the desired properties speculated by the programs. The computer-assisted analytical processes includes estimation of the maximum absorption and emission wavelengths, aquatic environmental toxicity, affinity to fibers, and photo-stability. This procedure could effectively narrow down discovery of new potential dye structures, greatly reduce and prevent complex and expensive preparation processes, and significantly improve the development efficiency of novel environmentally friendly dyes.

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.

Open-Source, Python-Based Redevelopment of the ChemShell Multiscale QM/MM Environment

ChemShell is a scriptable computational chemistry environment with an emphasis on multiscale simulation of complex systems using combined quantum mechanical and molecular mechanical (QM/MM) methods. Motivated by a scientific need to efficiently and accurately model chemical reactions on surfaces and within microporous solids on massively parallel computing systems, we present a major redevelopment of the ChemShell code, which provides a modern platform for advanced QM/MM embedding models. The new version of ChemShell has been re-engineered from the ground up with a new QM/MM driver module, an improved parallelization framework, new interfaces to high performance QM and MM programs, and a user interface written in the Python programming language. The redeveloped package is capable of performing QM/MM calculations on systems of significantly increased size, which we illustrate with benchmarks on zirconium dioxide nanoparticles of over 160000 atoms.

Electronic spectra of flavin in different redox and protonation states: a computational perspective on the effect of the electrostatic environment

This study discusses how UV/vis absorption spectra of flavin in different redox and protonation states are shifted by the nearby electrostatic microenvironment.

Blue-Light Receptors for Optogenetics

Sensory photoreceptors underpin light-dependent adaptations of organismal physiology, development, and behavior in nature. Adapted for optogenetics, sensory photoreceptors become genetically encoded actuators and reporters to enable the noninvasive, spatiotemporally accurate and reversible control by light of cellular processes. Rooted in a mechanistic understanding of natural photoreceptors, artificial photoreceptors with customized light-gated function have been engineered that greatly expand the scope of optogenetics beyond the original application of light-controlled ion flow. As we survey presently, UV/blue-light-sensitive photoreceptors have particularly allowed optogenetics to transcend its initial neuroscience applications by unlocking numerous additional cellular processes and parameters for optogenetic intervention, including gene expression, DNA recombination, subcellular localization, cytoskeleton dynamics, intracellular protein stability, signal transduction cascades, apoptosis, and enzyme activity. The engineering of novel photoreceptors benefits from powerful and reusable design strategies, most importantly light-dependent protein association and (un)folding reactions. Additionally, modified versions of these same sensory photoreceptors serve as fluorescent proteins and generators of singlet oxygen, thereby further enriching the optogenetic toolkit. The available and upcoming UV/blue-light-sensitive actuators and reporters enable the detailed and quantitative interrogation of cellular signal networks and processes in increasingly more precise and illuminating manners.

Electron transfer pathways in a light, oxygen, voltage (LOV) protein devoid of the photoactive cysteine

Blue-light absorption by the flavin chromophore in light, oxygen, voltage (LOV) photoreceptors triggers photochemical reactions that lead to the formation of a flavin-cysteine adduct. While it has long been assumed that adduct formation is essential for signaling, it was recently shown that LOV photoreceptor variants devoid of the photoactive cysteine can elicit a functional response and that flavin photoreduction to the neutral semiquinone radical is sufficient for signal transduction. Currently, the mechanistic basis of the underlying electron- (eT) and proton-transfer (pT) reactions is not well understood. We here reengineered pT into the naturally not photoreducible iLOV protein, a fluorescent reporter protein derived from the Arabidopsis thaliana phototropin-2 LOV2 domain. A single amino-acid substitution (Q489D) enabled efficient photoreduction, suggesting that an eT pathway is naturally present in the protein. By using a combination of site-directed mutagenesis, steady-state UV/Vis, transient absorption and electron paramagnetic resonance spectroscopy, we investigate the underlying eT and pT reactions. Our study provides strong evidence that several Tyr and Trp residues, highly conserved in all LOV proteins, constitute the eT pathway for flavin photoreduction, suggesting that the propensity for photoreduction is evolutionary imprinted in all LOV domains, while efficient pT is needed to stabilize the neutral semiquinone radical.

Dual Photochemical Reaction Pathway in Flavin-Based Photoreceptor LOV Domain: A Combined Quantum-Mechanics/Molecular-Mechanics Investigation

The primary photochemical reaction of the light, oxygen, and voltage (LOV) domain of the blue-light photosensor YtvA of Bacillus subtilis were investigated using high-level QM(DFT/MRCI)/MM methods. After blue-light excitation, the S_γ atom of the reactive cysteine forms a covalent bond with the C_4a of the flavin mononucleotide (FMN) ring. Two conformations for the side chain of reactive cysteine with occupancies of 70% (conf A) and 30% (conf B) are observed in the X-ray crystallographic structures of the YtvA-LOV ( Möglich , A. ; Moffat , K. J. Mol. Biol. 2007 , 373 , 112 - 126 ). In conf A, the thiol group is directed toward the dimethylbenzene moiety of the FMN ring whereas it is placed directly above the N₅ atom of the FMN ring in conf B. Starting from both conformations, the singlet and triplet excited pathways were evaluated. The singlet states excited from conf A decay nonradiatively to the triplet states by intersystem crossing (ISC). After the formation of a neutral biradical, the triplet states cross over to the electronic ground state by a second ISC and the adducts are efficiently formed. The singlet states excited from conf B are located near the S₁/S₀ conical intersection (CIn). A major fraction returns to the initial states through the CIn. The rest may directly reach the adduct state. Thus, the photoexcitation has a dual reaction pathway. In YtvA-LOV, it is inferred that the efficient triplet excitation from conf A was chosen by bypassing the less efficient singlet excitation from conf B.