The photochemistry of the light-, oxygen-, and voltage-sensitive domains in the algal blue light receptor phot

Electronic Structure Methods for Simulating Flavin's Spectroscopy and Photophysics: Comparison of Multi-reference, TD-DFT, and Single-Reference Wave Function Methods

The use of flavins and flavoproteins in photocatalytic, sensing, and biotechnological applications has led to a growing interest in computationally modeling the excited-state electronic structure and photophysics of flavin. However, there is limited consensus regarding which computational methods are appropriate for modeling flavin's photophysics. We compare the energies of low-lying excited states of flavin computed with time-dependent density functional theory (TD-DFT), equation-of-motion coupled cluster (EOM-EE-CCSD), scaled opposite-spin configuration interaction [SOS-CIS(D)], multiconfiguration pair-density functional theory (MC-PDFT), and several multireference perturbation theory (MR-PT2) methods. In the first part, we focus on excitation energies of the first singlet excited state (S1) of five different redox and protonation states of flavin, with the goal of finding a suitable active space for MR-PT2 calculations. In the second part, we construct two sets of one-dimensional potential energy surfaces connecting the S0 and S1 equilibrium geometries (S0-S1 path) and the S1 (π,π*) and S2 (n,π*) equilibrium geometries (S1-S2 path). The first path therefore follows a Franck-Condon active mode of flavin while the second path maps crossings points between low-lying singlet and triplet states in flavin. We discuss the similarities and differences in the TD-DFT, EOM-EE-CCSD, SOS-CIS(D), MC-PDFT and MR-PT2 energy profiles along these paths. We find that (TD-)DFT methods are suitable for applications such as simulating the spectra of flavins but are inconsistent with several other methods when used for some geometry optimizations and when describing the energetics of dark (n,π*) states. MR-PT2 methods show promise for the simulation of flavin's low-lying excited states, but the selection of orbitals for the active space and the number of roots used for state averaging must be done carefully to avoid artifacts. Some properties, such as the intersystem crossing geometry and energy between the S1 (π,π*) and T2 (n,π*) states, may require additional benchmarking before they can be determined quantitatively.

VMD as a Platform for Interactive Small Molecule Preparation and Visualization in Quantum and Classical Simulations

Modeling and simulation of small molecules such as drugs and biological cofactors have been both a major focus of computational chemistry for decades and a growing need among computational biophysicists who seek to investigate the interaction of different types of ligands with biomolecules. Of particular interest in this regard are quantum mechanical (QM) calculations that are used to more accurately describe such small molecules, which can be of heterogeneous structures and chemistry, either in purely QM calculations or in hybrid QM/molecular mechanics (MM) simulations. QM programs are also used to develop MM force field parameters for small molecules to be used along with established force fields for biomolecules in classical simulations. With this growing need in mind, here we report a set of software tools developed and closely integrated within the broadly used molecular visualization/analysis program, VMD, that allow the user to construct, modify, and parametrize small molecules and prepare them for QM, hybrid QM/MM, or classical simulations. The tools also provide interactive analysis and visualization capabilities in an easy-to-use and integrated environment. In this paper, we briefly report on these tools and their major features and capabilities, along with examples of how they can facilitate molecular research in computational biophysics that might be otherwise prohibitively complex.

Excited State Resonance Raman of Flavin Mononucleotide: Comparison of Theory and Experiment

Blue light absorbing flavoproteins play important roles in a variety of photobiological processes. Consequently, there have been numerous investigations of their excited state structure and dynamics, in particular by time-resolved vibrational spectroscopy. The isoalloxazine chromophore of the flavoprotein cofactors has been studied in detail by time-resolved Raman, lending it a benchmark status for mode assignments in excited electronic states of large molecules. However, detailed comparisons of calculated and measured spectra have proven challenging, as there are many more modes calculated than are observed, and the role of resonance enhancement is difficult to characterize in excited electronic states. Here we employ a recently developed approach due to Elles and co-workers ( J. Phys. Chem. A 2018, 122, 8308-8319) for the calculation of resonance-enhanced Raman spectra of excited states and apply it to the lowest singlet and triplet excited states of the isoalloxazine chromophore. There is generally good agreement between calculated and observed enhancements, which allows assignment of vibrational bands of the flavoprotein cofactors to be refined. However, some prominently enhanced bands are found to be absent from the calculations, suggesting the need for further development of the theory.

Steric and Electronic Interactions at Gln154 in ZEITLUPE Induce Reorganization of the LOV Domain Dimer Interface

Plants measure light quality, intensity, and duration to coordinate growth and development with daily and seasonal changes in environmental conditions; however, the molecular details linking photochemistry to signal transduction remain incomplete. Two closely related light, oxygen, or voltage (LOV) domain-containing photoreceptor proteins, ZEITLUPE (ZTL) and FLAVIN-BINDING, KELCH REPEAT, F-BOX 1 (FKF1), divergently regulate the protein stability of circadian clock and photoperiodic flowering components to mediate daily and seasonal development. Using structural approaches, we identified that mutations at the Gly46 position led to global rearrangements of the ZTL dimer interface in the isolated ZTL-LOV domain. Specifically, G46S and G46A variants induce a 180° rotation about the ZTL-LOV dimer interface that is coupled to ordering of N- and C-terminal signaling elements. These conformational changes hinge upon rotation of a C-terminal Gln residue (Gln154) analogous to that present in light-state structures of ZTL. In contrast to other LOV proteins, a Q154L variant retains light-state interactions with GIGANTEA (GI), thereby indicating N5 protonation is not required for ZTL signaling. The results presented herein confirm a divergent signaling mechanism within ZTL, whereby steric and electronic effects following adduct formation can be sufficient for signal propagation in LOV proteins containing a Gly residue at position 46. Examination of bacterial LOV structures with Gly residues at the equivalent position suggests that mechanisms of signal transduction in LOV proteins may be fluid across the LOV protein family.

Time-resolved fluctuation during the photochemical reaction of a photoreceptor protein: phototropin1LOV2-linker

Although the relationship between structural fluctuations and reactions is important for elucidating reaction mechanisms, experimental data describing such fluctuations of reaction intermediates are sparse. In order to investigate structural fluctuations during a protein reaction, the compressibilities of intermediate species after photoexcitation of a phot1LOV2-linker, which is a typical LOV domain protein with the C-terminal linker including the J-α helix and used recently for optogenetics, were measured in the time-domain by the transient grating and transient lens methods with a high pressure optical cell. The yield of covalent bond formation between the chromophore and a Cys residue (S state formation) relative to that at 0.1 MPa decreased very slightly with increasing pressure. The fraction of the reactive species that yields the T state (linker-unfolded state) decreased almost proportionally with pressure (0.1-200 MPa) to about 65%. Interestingly, the volume change associated with the reaction was much more pressure sensitive. By combining these data, the compressibility changes for the short lived intermediate (S state) and the final product (T state) formation were determined. The compressibility of the S state was found to increase compared with the dark (D) state, and the compressibility decreased during the transition from the S state to the T state. The compressibility change is discussed in terms of cavities inside the protein. By comparing the crystal structures of the phot1LOV2-linker at dark and light states, we concluded that the cavity volumes between the LOV domain and the linker domain increase in the S state, which explains the enhanced compressibility.

Go in! Go out! Inducible control of nuclear localization

Cells have evolved a variety of mechanisms to regulate the enormous complexity of processes taking place inside them. One mechanism consists in tightly controlling the localization of macromolecules, keeping them away from their place of action until needed. Since a large fraction of the cellular response to external stimuli is mediated by gene expression, it is not surprising that transcriptional regulators are often subject to stimulus-induced nuclear import or export. Here we review recent methods in chemical biology and optogenetics for controlling the nuclear localization of proteins of interest inside living cells. These methods allow researchers to regulate protein activity with exquisite spatiotemporal control, and open up new possibilities for studying the roles of proteins in a broad array of cellular processes and biological functions.

Protein stability: a crystallographer's perspective

Protein stability is a topic of major interest for the biotechnology, pharmaceutical and food industries, in addition to being a daily consideration for academic researchers studying proteins. An understanding of protein stability is essential for optimizing the expression, purification, formulation, storage and structural studies of proteins. In this review, discussion will focus on factors affecting protein stability, on a somewhat practical level, particularly from the view of a protein crystallographer. The differences between protein conformational stability and protein compositional stability will be discussed, along with a brief introduction to key methods useful for analyzing protein stability. Finally, tactics for addressing protein-stability issues during protein expression, purification and crystallization will be discussed.

Structural role of two histidines in the (6-4) photolyase reaction

Photolyases (PHRs) are DNA repair enzymes that revert UV-induced photoproducts, either cyclobutane pyrimidine dimers (CPD) or (6-4) photoproducts (PPs), into normal bases to maintain genetic integrity. (6-4) PHR must catalyze not only covalent bond cleavage, but also hydroxyl or amino group transfer, yielding a more complex mechanism than that postulated for CPD PHR. Previous mutation analysis revealed the importance of two histidines in the active center, H354 and H358 for Xenopus (6-4) PHR, whose mutations significantly lowered the enzymatic activity. Based upon highly sensitive FTIR analysis of the repair function, here we report that both H354A and H358A mutants of Xenopus (6-4) PHR still maintain their repair activity, although the efficiency is much lower than that of the wild type. Similar difference FTIR spectra between the wild type and mutant proteins suggest a common mechanism of repair in which (6-4) PP binds to the active center of each mutant, and is released after repair, as occurs in the wild type. Similar FTIR spectra also suggest that a decrease in volume by the H-to-A mutation is possibly compensated by the addition of water molecule( s). Such a modified environment is sufficient for the repair function that is probably controlled by proton-coupled electron transfer between the enzyme and substrate. On the other hand, two histidines must work in a concerted manner in the active center of the wild-type enzyme, which significantly raises the repair efficiency.

LOV-based optogenetic devices: light-driven modules to impart photoregulated control of cellular signaling

The Light-Oxygen-Voltage domain family of proteins is widespread in biology where they impart sensory responses to signal transduction domains. The small, light responsive LOV modules offer a novel platform for the construction of optogenetic tools. Currently, the design and implementation of these devices is partially hindered by a lack of understanding of how light drives allosteric changes in protein conformation to activate diverse signal transduction domains. Further, divergent photocycle properties amongst LOV family members complicate construction of highly sensitive devices with fast on/off kinetics. In the present review we discuss the history of LOV domain research with primary emphasis on tuning LOV domain chemistry and signal transduction to allow for improved optogenetic tools.

Exploring the multiscale signaling behavior of phototropin1 from Chlamydomonas reinhardtii using a full-residue space kinetic Monte Carlo molecular dynamics technique

Devising analysis tools for elucidating the regulatory mechanism of complex enzymes has been a challenging task for many decades. It generally requires the determination of the structural-dynamical information of protein solvent systems far from equilibrium over multiple length and time scales, which is still difficult both theoretically and experimentally. To cope with the problem, we introduce a full-residue space multiscale simulation method based on a combination of the kinetic Monte Carlo and molecular dynamics techniques, in which the rates of the rate-determining processes are evaluated from a biomolecular forcefield on the fly during the simulation run by taking into account the full space of residues. To demonstrate its reliability and efficiency, we explore the light-induced functional behavior of the full-length phototropin1 from Chlamydomonas reinhardtii (Cr-phot1) and its various subdomains. Our results demonstrate that in the dark state the light oxygen voltage-2-Jα (LOV2-Jα) photoswitch inhibits the enzymatic activity of the kinase, whereas the LOV1-Jα photoswitch controls the dimerization with the LOV2 domain. This leads to the repulsion of the LOV1-LOV2 linker out of the interface region between both LOV domains, which results in a positively charged surface suitable for cell-membrane interaction. By contrast, in the light state, we observe that the distance between both LOV domains is increased and the LOV1-LOV2 linker forms a helix-turn-helix (HTH) motif, which enables gene control through nucleotide binding. Finally, we find that the kinase is activated through the disruption of the Jα-helix from the LOV2 domain, which is followed by a stretching of the activation loop (A-loop) and broadening of the catalytic cleft of the kinase.

FTIR spectroscopy of flavin-binding photoreceptors

Light-induced difference Fourier transform infrared (FTIR) spectroscopy is a powerful, sensitive, and informative method to study structure-function relationships in photoreceptive proteins. Strong absorption of water in the IR region is always problematic in this method, but if water content in the sample is controlled during measurements, this method can provide useful information on a single protein-bound water molecule. We established three kinds of sample preparations: hydrated film, redissolved sample, and concentrated solution. Hydrated films were used for the measurements of LOV and BLUF domains, where accurate difference FTIR spectra were obtained in the whole mid-IR region (4,000-800 cm(-1)). Vibrations of S-H stretch of cysteine, O-H stretch of water, and O-H stretch of tyrosine provide important information on hydrogen bonds in these proteins. Redissolved samples were used for the measurements of (6-4) photolyase, in which enzymatic turnover takes place. From the illumination time-dependence of excess amount of substrate, it is possible to isolate the signal originating from the binding of enzyme to substrate. If proteins are less tolerant to drying, as for example cryptochromes of the DASH type, concentrated solution is used. Detailed methodological aspects in light-induced difference FTIR spectroscopy are reviewed by mainly focusing on our results.

Fourier-transform infrared study of the photoactivation process of Xenopus (6-4) photolyase

Photolyases (PHRs) are blue light-activated DNA repair enzymes that maintain genetic integrity by reverting UV-induced photoproducts into normal bases. The flavin adenine dinucleotide (FAD) chromophore of PHRs has four different redox states: oxidized (FAD(ox)), anion radical (FAD(•-)), neutral radical (FADH(•)), and fully reduced (FADH(-)). We combined difference Fourier-transform infrared (FTIR) spectroscopy with UV-visible spectroscopy to study the detailed photoactivation process of Xenopus (6-4) PHR. Two photons produce the enzymatically active, fully reduced PHR from oxidized FAD: FAD(ox) is converted to semiquinone via light-induced one-electron and one-proton transfers and then to FADH(-) by light-induced one-electron transfer. We successfully trapped FAD(•-) at 200 K, where electron transfer occurs but proton transfer does not. UV-visible spectroscopy following 450 nm illumination of FAD(ox) at 277 K defined the FADH(•)/FADH(-) mixture and allowed calculation of difference FTIR spectra among the four redox states. The absence of a characteristic C=O stretching vibration indicated that the proton donor is not a protonated carboxylic acid. Structural changes in Trp and Tyr are suggested by UV-visible and FTIR analysis of FAD(•-) at 200 K. Spectral analysis of amide I vibrations revealed structural perturbation of the protein's β-sheet during initial electron transfer (FAD(•-) formation), a transient increase in α-helicity during proton transfer (FADH(•) formation), and reversion to the initial amide I signal following subsequent electron transfer (FADH(-) formation). Consequently, in (6-4) PHR, unlike cryptochrome-DASH, formation of enzymatically active FADH(-) did not perturb α-helicity. Protein structural changes in the photoactivation of (6-4) PHR are discussed on the basis of these FTIR observations.

Signals of LOV1: a computer simulation study on the wildtype LOV1-domain of Chlamydomonas reinhardtii and its mutants

Phototropins are photoreceptors regulating the blue-light response in plants and bacteria. They consist of two LOV (light oxygen voltage sensitive) domains each containing a non-covalently bound flavin-mononucleotide (FMN) chromophore, which are connected to a serine/threonine-kinase. Upon illumination, the LOV-domains undergo conformational changes, triggering a signal cascade in the organism through kinase activation. Here, we present results from molecular dynamics simulations in which we investigate the signal transduction pathway of the wildtype LOV1-domain of Chlamydomonas reinhardtii and a methyl-mercaptan (MM) adduct of its Cys57Gly-mutant at the molecular level. In particular, we analyzed the effect of covalent-bond formation between the reactive cysteine Cys57 and the FMN-reaction center, as well as the subsequent charge redistribution, on the spatio-dynamical behavior of the LOV1-domain. We compare the calculation results with experimental data and demonstrate that these adduct state characteristics have an important influence on the response of this photosensor. The light-induced changes implicate primarily an alteration of the surface charge distribution through rearrangement of the highly flexible Cα-, Dα- and Eα-helices including the Glu51-Lys91-salt bridge on the hydrophilic side of the protein domain and a β-sheet tightening process via coupling of the Aβ- and Bβ-strands. Our findings confirm the aptitude of the LOV1-domain to function as a dimerization partner, allowing the green alga to adapt its reproduction and growth speed to the environmental conditions.

Flavin redox switching of protein functions

Flavin cofactors impart remarkable catalytic diversity to enzymes, enabling them to participate in a broad array of biological processes. The properties of flavins also provide proteins with a versatile redox sensor that can be utilized for converting physiological signals such as cellular metabolism, light, and redox status into a unique functional output. The control of protein functions by the flavin redox state is important for transcriptional regulation, cell signaling pathways, and environmental adaptation. A significant number of proteins that have flavin redox switches are found in the Per-Arnt-Sim (PAS) domain family and include flavoproteins that act as photosensors and respond to changes in cellular redox conditions. Biochemical and structural studies of PAS domain flavoproteins have revealed key insights into how flavin redox changes are propagated to the surface of the protein and translated into a new functional output such as the binding of a target protein in a signaling pathway. Mechanistic details of proteins unrelated to the PAS domain are also emerging and provide novel examples of how the flavin redox state governs protein-membrane interactions in response to appropriate stimuli. Analysis of different flavin switch proteins reveals shared mechanistic themes for the regulation of protein structure and function by flavins.

Effect of computational methodology on the conformational dynamics of the protein photosensor LOV1 from Chlamydomonas reinhardtii

LOV domains are the light-sensitive protein domains of plant phototropins and bacteria. They photochemically form a covalent bond between a flavin mononucleotide (FMN) chromophore and a cysteine, attached to the apo-protein, upon irradiation with blue light, which triggers a signal in the adjacent kinase. Although their signaling state has been well characterized through experimental means, their signal transduction pathway as well as dark-state activity are generally only poorly understood. Here we show results from molecular dynamics simulations where we investigated the effect of thermostating and long-range electrostatics on the solution structure and dynamical behavior of the wild-type LOV1 domain from the green algae Chlamydomonas reinhardtii in the dark. We demonstrate that these computational issues can dramatically affect the conformational fluctuations of such protein domains by suppressing configurations far from equilibrium or destabilizing local configurations, leading to artificial changes of the protein secondary structure as well as the H-bond network formed by the amino acids and the FMN. By comparing our calculation results with recent experimental data, we show that the non-invasive thermostating strategy, where the protein solute is only indirectly coupled to the thermostat via the solvent, in conjunction with the particle-mesh Ewald technique, provides dark-state conformers, which are in consistency with experimental observations. Moreover, our calculations indicate that the LOV1 domains can alter the intersystem crossing rate and rate of adduct formation by adjusting the population distribution of these dark-state conformers. This might permit them to function as a modulator of the signal intensity under low light conditions.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (doi:10.1007/s12154-011-0060-z) contains supplementary material, which is available to authorized users.

FMN binding and photochemical properties of plant putative photoreceptors containing two LOV domains, LOV/LOV proteins

LOV domains function as blue light-sensing modules in various photoreceptors in plants, fungi, algae, and bacteria. A LOV/LOV protein (LLP) has been found from Arabidopsis thaliana (AtLLP) as a two LOV domain-containing protein. However, its function remains unknown. We isolated cDNA clones coding for an LLP homolog from tomato (Solanum lycopersicum) and two homologs from the moss Physcomitrella patens. The tomato LLP (SlLLP) contains two LOV domains (LOV1 and LOV2 domains), as in AtLLP. Most of the amino acids required for association with chromophore are conserved in both LOV domains, except that the amino acid at the position equivalent to the cysteine essential for cysteinyl adduct formation is glycine in the LOV1 domain as in AtLLP. When expressed in Escherichia coli, SlLLP binds FMN and undergoes a self-contained photocycle upon irradiation of blue light. Analyses using mutant SlLLPs revealed that SlLLP binds FMN in both LOV domains, although the LOV1 domain does not show spectral changes on irradiation. However, when Gly(66) in the LOV1 domain, which is located at the position equivalent to the essential cysteine of LOV domains, is replaced by cysteine, the mutated LOV1 domain shows light-induced spectral changes. In addition, all four LOV domains of P. patens LLPs (PpLLP1 and PpLLP2) show the typical features of LOV domains, including the reactive cysteine in each. This study shows that plants have a new LOV domain-containing protein family with the typical biochemical and photochemical properties of other LOV domain-containing proteins such as the phototropins.

How the green alga Chlamydomonas reinhardtii keeps time

The unicellular green alga Chlamydomonas reinhardtii has two flagella and a primitive visual system, the eyespot apparatus, which allows the cell to phototax. About 40 years ago, it was shown that the circadian clock controls its phototactic movement. Since then, several circadian rhythms such as chemotaxis, cell division, UV sensitivity, adherence to glass, or starch metabolism have been characterized. The availability of its entire genome sequence along with homology studies and the analysis of several sub-proteomes render C. reinhardtii as an excellent eukaryotic model organism to study its circadian clock at different levels of organization. Previous studies point to several potential photoreceptors that may be involved in forwarding light information to entrain its clock. However, experimental data are still missing toward this end. In the past years, several components have been functionally characterized that are likely to be part of the oscillatory machinery of C. reinhardtii since alterations in their expression levels or insertional mutagenesis of the genes resulted in defects in phase, period, or amplitude of at least two independent measured rhythms. These include several RHYTHM OF CHLOROPLAST (ROC) proteins, a CONSTANS protein (CrCO) that is involved in parallel in photoperiodic control, as well as the two subunits of the circadian RNA-binding protein CHLAMY1. The latter is also tightly connected to circadian output processes. Several candidates including a significant number of ROCs, CrCO, and CASEIN KINASE1 whose alterations of expression affect the circadian clock have in parallel severe effects on the release of daughter cells, flagellar formation, and/or movement, indicating that these processes are interconnected in C. reinhardtii. The challenging task for the future will be to get insights into the clock network and to find out how the clock-related factors are functionally connected. In this respect, system biology approaches will certainly contribute in the future to improve our understanding of the C. reinhardtii clock machinery.

Photoreaction of mutated LOV photoreceptor domains from Chlamydomonas reinhardtii with aliphatic mercaptans: implications for the mechanism of wild type LOV

Irradiation of the LOV1 domain from the blue-light photoreceptor phototropin of the green alga Chlamydomonas reinhardtii leads to the formation of a covalent adduct of the sulfur atom of cysteine 57 to the carbon C(4a) in the chromophore FMN. This reaction is not possible in the mutant LOV1-C57G in which this cysteine is replaced by glycine. Irradiation of LOV1-C57G in the absence of oxygen but in the presence of aliphatic mercaptans or thioethers leads to the formation of a species with an absorption maximum at 615 nm, which is identified as the neutral radical FMNH . When oxygen is admitted, the reaction is completely reversible. Irradiation of LOV1-C57G in the presence of methylmercaptan CH(3)SH under oxygen-free conditions yields, in addition to FMNH , a third species with a single absorption maximum at 379 nm. This species is stable against oxygen and is also formed when the irradiation is performed in the presence of oxygen. This species is assigned to the adduct between CH(3)SH and FMN. In aqueous solution the photoreaction of CH(3)SH with FMN leads to the fully reduced hydroquinone form FMNH(2) or its anion FMNH(-). Adduct formation apparently requires the protein cage. After formation, the adduct is stable for hours inside the protein, but decomposes immediately upon denaturation. The implications of these observations for the mechanism of adduct formation in wild type LOV domains are discussed.

Phototropin blue-light receptors

Phototropins are blue-light receptors controlling a range of responses that serve to optimize the photosynthetic efficiency of plants. These include phototropism, light-induced stomatal opening, and chloroplast movements in response to changes in light intensity. Since the isolation of the Arabidopsis PHOT1 gene in 1997, phototropins have been identified in ferns and mosses where their physiological functions appear to be conserved. Arabidopsis contains two phototropins, phot1 and phot2, that exhibit overlapping functions in addition to having unique physiological roles. Phototropins are light-activated serine/threonine protein kinases. Light sensing by the phototropins is mediated by a repeated motif at the N-terminal region of the protein known as the LOV domain. Photoexcitation of the LOV domain results in receptor autophosphorylation and an initiation of phototropin signaling. Here we summarize the photochemical and biochemical events underlying phototropin activation in addition to the current knowledge of the molecular mechanisms associated with photoreceptor signaling.