Computational characterization of reaction intermediates in the photocycle of the sensory domain of the AppA blue light photoreceptor

The AppA protein with the BLUF (blue light using flavin adenine dinucleotide) domain is a blue light photoreceptor that cycle between dark-adapted and light-induced functional states. We characterized possible reaction intermediates in the photocycle of AppA BLUF. Molecular dynamics (MD), quantum chemical and quantum mechanical-molecular mechanical (QM/MM) calculations were carried out to describe several stable structures of a molecular system modeling the protein. The coordinates of heavy atoms from the crystal structure (PDB code 2IYG) of the protein in the dark state served as starting point for 10 ns MD simulations. Representative MD frames were used in QM(B3LYP/cc-pVDZ)/MM(AMBER) calculations to locate minimum energy configurations of the model system. Vertical electronic excitation energies were estimated for the molecular clusters comprising the quantum subsystems of the QM/MM optimized structures using the SOS-CIS(D) quantum chemistry method. Computational results support the occurrence of photoreaction intermediates that are characterized by spectral absorption bands between those of the dark and light states. They agree with crystal structures of reaction intermediates (PDB code 2IYI) observed in the AppA BLUF domain. Transformations of the Gln63 side chain stimulated by photo-excitation and performed with the assistance of the chromophore and the Met106 side chain are responsible for these intermediates.     

Structures of the chromophore binding sites in BLUF domains as studied by molecular dynamics and quantum chemical calculations

BLUF (blue-light sensing using FAD) domains constitute a new family of flavin-based blue light photoreceptors. The photocycle of BLUF is unique in the sense that a few hydrogen bond rearrangements are accompanied by only slight structural changes in the bound chromophore. The hydrogen bond rearrangements upon illumination have been inferred from spectral changes in the chromophore: approximately 10 nm redshift of the absorption maximum and approximately 16 cm(-1) downshift of the C4=O stretching frequency. However, the exact features of the hydrogen bond network around the active site are still the subject of some controversy. In particular, the orientation of a conserved Gln (Gln63 in AppA) is presently one of the most questioned topics in the field. Here we perform molecular dynamics simulations for the wild-type AppA, AppA1-124C20S, BlrB and T110078 and furthermore quantum chemical calculations to investigate their spectroscopic properties in the dark and signaling states. On the basis of these results, we reveal the dynamic aspect of hydrogen bonding networks at the active site and propose theoretically reasonable models for the dark and signaling states of the BLUF domains.     

Purification and initial characterization of a putative blue light-regulated phosphodiesterase from Escherichia coli

The Escherichia coli protein YcgF contains a photosensory flavin adenine dinucleotide (FAD)-binding BLUF domain covalently linked to an EAL domain, which is predicted to have cyclic-di-guanosine monophosphate (GMP) phosphodiesterase activity. We have cloned, overexpressed and purified this protein, which we refer to as blue light-regulated phosphodiesterase (Blrp) for its putative activity. Blrp undergoes a reversible photocycle after exposure to light in which the spectrum of its photostationary state and kinetics of recovery of the dark state are similar to those of the isolated BLUF domain of the AppA protein. Unlike the AppA BLUF domain, the chromophore environment in the context of full-length Blrp is asymmetric, and the protein does not undergo any detectable global changes on exposure to blue light. When overexpressed in E. coli, Blrp copurifies with certain proteins which suggests that it plays a protective role in response to oxidative stress. Predicted proteins from Klebsiella pneumoniae and from a bacterium in the Sargasso Sea are similar to E. coli Blrp in both their BLUF and EAL domains, which suggests that blue light sensing in these bacteria may follow similar pathways.     

Orientation of a key glutamine residue in the BLUF domain from AppA revealed by mutagenesis, spectroscopy, and quantum chemical calculations

The flavin-adenine-dinucleotide-binding BLUF domain constitutes a new class of blue-light receptors, and the N-terminal domain of AppA is a representative of this family. A crystal structure of the BLUF domain from AppA suggested that a conserved Gln63 forms a hydrogen bond with the flavin N5 atom. Upon light excitation, this residue is proposed to undergo a approximately 180 degrees rotation that leads to a rearrangement of a hydrogen bonding network. However, crystallographic studies on the other BLUF proteins claimed an opposite orientation for the glutamine residue. In this communication, we have revealed the presence of a Gln63-to-N5 hydrogen bond in the dark state of AppA by a combined approach of mutagenesis, spectroscopy, and quantum chemical calculations. The present finding supports the view that the reorientation of the Gln63 side chain is a key event in the signaling state formation of BLUF proteins.     

LOV-like flavin-Cys adduct formation by introducing a Cys residue in the BLUF domain of TePixD

BLUF and LOV are blue-light sensor domains that possess flavin as a common chromophore but exhibit distinct photoreactions. Ile66 located in the BLUF domain of a cyanobacterial photosensor protein, TePixD, was replaced with Cys to mimic the LOV domain. Light-induced Fourier transform infrared spectra of the I66C TePixD showed that a flavin-Cys adduct, typical of the photoinduced intermediates of LOV domains, was formed in the I66C BLUF domain. This result demonstrates that different types of flavin photoreactions can be realized in the same domain if key amino acids are properly arranged near the flavin and the domain structure itself is not a crucial factor to determine the photoreaction type.     

Light-induced flipping of a conserved glutamine sidechain and its orientation in the AppA BLUF domain

The AppA BLUF domain is a blue light photoreceptor containing flavin. Conserved glutamine 63 is necessary for the photocycle of the protein, and its side chain has been proposed to flip in response to blue light illumination. Recently published crystal structures of AppA WT and the AppA mutant C20S describe contradictory conclusions regarding the orientation of the conserved glutamine 63 side chain in the dark. Here, we present evidence from NMR spectroscopy confirming light-induced flipping of the glutamine side chain to form a strong hydrogen bond between the glutamine 63 side chain carbonyl group and the tyrosine 21 side chain hydroxyl proton in the light-induced state. Our conclusions are consistent with published data from UV/vis absorbance and FTIR spectroscopy, as well as the crystal structure of AppA WT.     

BLUF Hydrogen network dynamics and UV/Vis spectra: A combined molecular dynamics and quantum chemical study

Blue light sensing using flavin (BLUF) protein photoreceptor domains change their hydrogen bond network after photoexcitation. To explore this phenomenon, BLUF domains from R. sphaeroides were simulated using Amber99 molecular dynamics (MD). Five starting configurations were considered, to study different BLUF proteins (AppA/BlrB), Trp conformations (“W_in”/“W_out”), structure determination (X‐ray/NMR), and finally, His protonation states. We found dependencies of the hydrogen bonds on almost all parameters. Our data show an especially strong correlation of the Trp position and hydrogen bonds involving Gln63. The latter is in some contradiction to earlier results (Obanayama et al., Photochem. Photobiol. 2008, 84 10031010). Possible origins and implications are discussed. Our calculations support conjectures that Gln63 is more flexible with Trp104 in W_in position. Using snapshots from MD and time‐dependent density functional theory, UV/vis spectra for the chromophore were determined, which account for molecular motion of the protein under ambient conditions. In accord with experiment, it is found that the UV/vis spectra of BLUF bound flavin are red‐shifted and thermally broadened for all calculated π → π* transitions, relative to gas phase flavin at T = 0 K. However, differences in the spectra between the various BLUF configurations cannot be resolved with the present approach. © 2012 Wiley Periodicals, Inc.     

Initial characterization of the primary photochemistry of AppA, a blue-light-using flavin adenine dinucleotide-domain containing transcriptional antirepressor protein from Rhodobacter sphaeroides: a key role for reversible intramolecular proton transfer from the flavin adenine dinucleotide chromophore to a conserved tyrosine?

The flavin adenine dinucleotide (FAD)-containing photoreceptor protein AppA (in which the FAD is bound to a novel so-called BLUF domain) from the purple nonsulfur bacterium Rhodobacter sphaeroides was previously shown to be photoactive by the formation of a slightly redshifted long-lived intermediate that is thought to be the signaling state. In this study, we provide further characterization of the primary photochemistry of this photoreceptor protein using UV-Vis and Fourier-transform infrared spectroscopy, pH measurements and site-directed mutagenesis. Available evidence indicates that the FAD chromophore of AppA may be protonated in the receptor state, and that it becomes exposed to solvent in the signaling state. Furthermore, experimental data lead to the suggestion that intramolecular proton transfer (that may involve [anionic] Tyr-17) forms the basis for the stabilization of the signaling state.     

Light-generated paramagnetic intermediates in BLUF domains

Blue-light sensitive photoreceptory BLUF domains are flavoproteins, which regulate various, mostly stress-related processes in bacteria and eukaryotes. The photoreactivity of the flavin adenine dinucleotide (FAD) cofactor in three BLUF domains from Rhodobacter sphaeroides, Synechocystis sp. PCC 6803 and Escherichia coli have been studied at low temperature using time-resolved electron paramagnetic resonance. Photoinduced flavin triplet states and radical-pair species have been detected on a microsecond time scale. Differences in the electronic structures of the FAD cofactors as reflected by altered zero-field splitting parameters of the triplet states could be correlated with changes in the amino-acid composition of the various BLUF domains' cofactor binding pockets. For the generation of the light-induced, spin-correlated radical-pair species in the BLUF domain from Synechocystis sp. PCC 6803, a tyrosine residue near the flavin's isoalloxazine moiety plays a critical role.     

Theoretical study of proton coupled electron transfer reaction in the light state of the AppA BLUF photoreceptor

The BLUF (blue light sensor using flavin adenine dinucleotide) domain is widely studied as a prototype for proton coupled electron transfer (PCET) reactions in biological systems. In this work, the photo-induced concerted PCET reaction from the light state of the AppA BLUF domain is investigated. To model the simultaneous transfer of two protons in the reaction, two-dimensional potential energy surfaces for the double proton transfer are first calculated for the locally excited and charge transfer states, which are then used to obtain the vibrational wave function overlaps and the vibrational energy levels. Contributions to the PCET rate constant from each pair of vibronic states are then analyzed using the theory based on the Fermi's golden rule. We show that, the recently proposed light state structure of the BLUF domain with a tautomerized Gln63 residue is consistent with the concerted transfer of one electron and two protons. It is also found that, thermal fluctuations of the protein structure, especially the proton donor-acceptor distances, play an important role in determining the PCET reaction rate. © 2018 Wiley Periodicals, Inc.     

Simultaneous analysis of photoinduced electron transfer in wild type and mutated AppAs

Photoinduced electron transfer (PET) from Tyr21 to isoalloxazine (Iso) in the excited state (Iso*) is considered to be an initial step of the photosensing function of the blue-light sensing using flavin adenine dinucleotide (BLUF) component of the anti-repressor of the photosynthetic regulation (AppA). The PET mechanism was investigated via fluorescence dynamics of AppA and Kakitani and Mataga (KM) theories as well as by molecular dynamic (MD) simulation. The local structures of both the Y21F and W104F mutant AppAs around the Iso binding sites were quite different from those of the wild type (WT) AppA. The distances between Iso and Trp104 in Y21F, and between Iso and Tyr21 in W104F were shorter by 0.06 nm and 0.02 nm, respectively, compared to the WT. The frequency factor, ν₀, in Tyr21 was 1.15-fold greater than that in Trp104. The critical distance between adiabatic and non-adiabatic PET processes, R₀, was found to be very long in the AppA Tyr21. The large values of ν₀ and R₀ for Tyr21 of AppA compared to those in a non photosensing flavoprotein, FMN binding protein (FBP), were elucidated by hydrogen bond (H bond) chain between Tyr21 and Iso through Gln63. Interaction energies among Iso*, Trp104, Tyr21 and Gln63 in WT were calculated using the semi-empirical PM3 method. The amount of the transferred charge from Trp104 to Iso* in the WT exhibited a maximum at an interaction energy of around ?20 kcal/mol, but decreased as the interaction energy (absolute value) increased.     

Ultrafast vibrational spectroscopy of the flavin chromophore

Ultrafast time-resolved infrared (TRIR) spectra of flavin adenine dinucleotide (FAD) and the anion of lumiflavin (Lf-) are described. Ground-state recovery and excited-state decay of FAD reveal a common dominant ultrafast relaxation and a minor slower component. The Lf- transient lacks a fast component. No intermediate species are observed, suggesting that the quenching mechanism is internal conversion promoted by interaction of the adenine and isoalloxazine rings in FAD. Modes are assigned, and the potential for extension of the TRIR method to photoactive proteins is discussed.     

From crystallographic data to the solution structure of photoreceptors: the case of the AppA BLUF domain

Photoreceptor proteins bind a chromophore, which, upon light absorption, modifies its geometry or its interactions with the protein, finally inducing the structural change needed to switch the protein from an inactive to an active or signaling state. In the Blue Light-Using Flavin (BLUF) family of photoreceptors, the chromophore is a flavin and the changes have been connected with a rearrangement of the hydrogen bond network around it on the basis of spectroscopic changes measured for the dark-to-light conversion. However, the exact conformational change triggered by the photoexcitation is still elusive mainly because a clear consensus on the identity not only of the light activated state but also of the dark one has not been achieved. Here, we present an integrated investigation that combines microsecond MD simulations starting from the two conflicting crystal structures available for the AppA BLUF domain with calculations of NMR, IR and UV-Vis spectra using a polarizable QM/MM approach. Thanks to such a combined analysis of the three different spectroscopic responses, a robust characterization of the structure of the dark state in solution is given together with the uncovering of important flaws of the most popular molecular mechanisms present in the literature for the dark-to-light activation.

With an integrated molecular dynamics and QM/MM strategy we characterize the dark-state structure of a BLUF photoreceptor and ration alize the discrepancy between published crystal structures.     

Photoreaction in BLUF receptors: proton-coupled electron transfer in the flavin-Gln-Tyr system

Photoinduced electron transfer from tyrosine to the flavin chromophore is involved in activation of BLUF (sensor of blue light using FAD) photoreceptors. We studied the electron transfer (ET) coupled with proton-transfer (PT) reactions, by means of XMCQDPT2//CASSCF calculations on a molecular cluster model. By defining a minimum active space in the CASSCF calculations, we could compute the entire photoreaction pathway. We find that the crossing of the locally excited and ET states is located along the flavin bond-stretching coordinate. The ET state is stabilized by a proton transfer from the electron donor to the electron acceptor. We mapped two different PT pathways from tyrosine to flavin via the conserved glutamine. These reactions generate a tautomeric form of glutamine. Along the PT coordinates, we find geometries where the ET and the electronic ground states degenerate. At the state crossing structures, either formation of the ground state biradical intermediate or a relaxation back to the Franck-Condon minimum takes places. The computed relaxation pathways reveal that the hydrogen bonds involving glutamine in the chromophore-binding pocket control BLUF photoefficiency.     

Why BLUF photoreceptors with roseoflavin cofactors lose their biological functionality

The photophysics of roseoflavin in three different environments is investigated by using ab initio and quantum mechanics/molecular mechanics methods. Intramolecular charge transfer is shown to be responsible for the quenching of the fluorescence in the gas phase, and in the water environment. However, for the roseoflavin incorporated into the blue light using flavin (BLUF) protein environment (substituting the native flavin) no such deactivation is found. The conical intersection between the locally excited state of the chromophore and the charge transfer state involving the tyrosine residue, which in the native BLUF domain is responsible for initiating the photocycle, is missing for the roseoflavin substituted protein. This explains the experimental observations of the lack of any photocycle, and the loss of the biological function of the BLUF photoreceptor reported earlier.     

Modulation of the photocycle of a LOV domain photoreceptor by the hydrogen-bonding network

An extended hydrogen-bonding (HB) network stabilizes the isoalloxazine ring of the flavin mononucleotide (FMN) chromophore within the photosensing LOV domain of blue-light protein receptors, via interactions between the C(2)═O, N(3)H, C(4)═O, and N(5) groups and conserved glutamine and asparagine residues. In this work we studied the influence of the HB network on the efficiency, kinetics, and energetics of a LOV protein photocycle, involving the reversible formation of a FMN-cysteine covalent adduct. The following results were found for mutations of the conserved amino acids N94, N104, and Q123 in the Bacillus subtilis LOV protein YtvA: (i) Increased (N104D, N94D) or strongly reduced (N94A) rate of adduct formation; this latter mutation extends the lifetime of the flavin triplet state, i.e., adduct formation, more than 60-fold, from 2 μs for the wild-type (WT) protein to 129 μs. (ii) Acceleration of the overall photocycle for N94S, N94A, and Q123N, with recovery lifetimes 20, 45, and 85 times faster than for YtvA-WT, respectively. (iii) Slight modifications of FMN spectral features, correlated with the polarization of low-energy transitions. (iv) Strongly reduced (N94S) or suppressed (Q123N) structural volume changes accompanying adduct formation, as determined by optoacoustic spectroscopy. (v) Minor effects on the quantum yield, with the exception of a considerable reduction for Q123N, i.e., 0.22 vs 0.49 for YtvA-WT. The data stress the importance of the HB network in modulating the photocycle of LOV domains, while at the same time establishing a link with functional responses.     

Structure and photoreaction of photoactive yellow protein, a structural prototype of the PAS domain superfamily

Photoactive yellow protein (PYP) is a water-soluble photosensor protein found in purple photosynthetic bacteria. Unlike bacterial rhodopsins, photosensor proteins composed of seven transmembrane helices and a retinal chromophore in halophilic archaebacteria, PYP is a highly soluble globular protein. The alpha/beta fold structure of PYP is a structural prototype of the PAS domain superfamily, many members of which function as sensors for various kinds of stimuli. To absorb a photon in the visible region, PYP has a p-coumaric acid chromophore binding to the cysteine residue via a thioester bond. It exists in a deprotonated trans form in the dark. The primary photochemical event is photo-isomerization of the chromophore from trans to cis form. The twisted cis chromophore in early intermediates is relaxed and finally protonated. Consequently, the chromophore becomes electrostatically neutral and rearrangement of the hydrogen-bonding network triggers overall structural change of the protein moiety, in which local conformational change around the chromophore is propagated to the N-terminal region. Thus, it is an ideal model for protein conformational changes that result in functional change, responding to stimuli and expressing physiological activity. In this paper, recent progress in investigation of the photoresponse of PYP is reviewed.