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.     

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.     

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.     

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.     

Photoexcitation of the blue light using FAD photoreceptor AppA results in ultrafast changes to the protein matrix

Photoexcitation of the flavin chromophore in the BLUF photosensor AppA results in a conformational change that leads to photosensor activation. This conformational change is mediated by a hydrogen-bonding network that surrounds the flavin, and photoexcitation is known to result in changes in the network that include a strengthening of hydrogen bonding to the flavin C4═O carbonyl group. Q63 is a key residue in the hydrogen-bonding network, and replacement of this residue with a glutamate results in a photoinactive mutant. While the ultrafast time-resolved infrared (TRIR) spectrum of Q63E AppA(BLUF) is characterized by flavin carbonyl modes at 1680 and 1650 cm(-1), which are similar in frequency to the analogous modes from the light activated state of the wild-type protein, a band is also observed in the TRIR spectrum at 1724 cm(-1) that is unambiguously assigned to the Q63E carboxylic acid based on U-(13)C labeling of the protein. Light absorption instantaneously (<100 fs) bleaches the 1724 cm(-1) band leading to a transient absorption at 1707 cm(-1). Because Q63E is not part of the isoalloxazine electronic transition, the shift in frequency must arise from a sub picosecond perturbation to the flavin binding pocket. The light-induced change in the frequency of the Q63E side chain is assigned to an increase in hydrogen-bond strength of 3 kcal mol(-1) caused by electronic reorganization of the isoalloxazine ring in the excited state, providing direct evidence that the protein matrix of AppA responds instantaneously to changes in the electronic structure of the chromophore and supporting a model for photoactivation of the wild-type protein that involves initial tautomerization of the Q63 side chain.     

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.     

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.     

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.     

In Vivo Effects on Photosynthesis Gene Expression of Base Pair Exchanges in the Gene Encoding the Light-responsive BLUF Domain of AppA in Rhodobacter Sphaeroides

The Rhodobacter sphaeroides protein AppA has the unique quality of sensing and transmitting light and redox signals. By acting as antirepressor to the PpsR protein, it acts as a major regulator in photosynthesis gene expression. In this study, we show that by introducing amino acid exchanges into the AppA protein, the in vivo activity as an antirepressor can be greatly altered. The tryptophan 104 to phenylalanine (W104F) base exchange greatly diminished blue-light sensitivity of the BLUF domain. From the obtained in vivo data, the difference in thermal recovery rate of the signaling state of the BLUF domain between the wild type and mutated protein was calculated, predicting an about 10-fold faster recovery in the mutant, which is consistent with in vitro data. Introduction of a tyrosine 21 to phenylalanine (Y21F) or to cysteine (Y21C) mutation led to a complete loss of AppA antirepressor activity, while additionally leading to an increase of photosynthesis gene expression after illumination with high blue-light quantities. Interestingly, this effect is not visible in a W104F/Y21F double mutant that again shows a wild-type–like behavior of the BLUF domain after blue-light illumination, thus restoring the activity of AppA.     

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.     

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.     

The photoreaction mechanism in the bacterial blue light receptor BLUF according to metadynamics modeling

The mechanism of chemical transformations in the blue light photoreceptor domains (BLUF) implies the isomerization of the glutamine side chain. The Helmholtz energy profiles for the side-chain isomerization of the tautomeric form of glutamine in the BLUF domain of the bacterial protein AppA were calculated using metadynamics and the potentials that were obtained using quantum mechanics-molecular mechanics approximation (QM/MM).     

Inherent Conformational Preferences of Ac‐Gln‐Gln‐NHBn: Sidechain Hydrogen Bonding Supports a β‐Turn in the Gas Phase

Gas‐phase single‐conformation spectroscopy is used to study Ac‐Gln‐Gln‐NHBn in order to probe the interplay between sidechain hydrogen bonding and backbone conformational preferences. This small, amide‐rich peptide offers many possibilities for backbone–backbone, sidechain–backbone, and sidechain–sidechain interactions. The major conformer observed experimentally features a type‐I β‐turn with a canonical 10‐membered ring C=O—H?N hydrogen bond between backbone amide groups. In addition, the C=O group of each Gln sidechain participates in a seven‐membered ring hydrogen bond with the backbone NH of the same residue. Thus, sidechain hydrogen‐bonding potential is satisfied in a manner that is consistent with and stabilizes the β‐turn secondary structure. This turn‐forming propensity may be relevant to pathogenic amyloid formation by polyglutamine segments in human proteins.     

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.     

Toward assessing the position-dependent contributions of backbone hydrogen bonding to beta-sheet folding thermodynamics employing amide-to-ester perturbations

An amide-to-ester backbone substitution in a protein is accomplished by replacing an alpha-amino acid residue with the corresponding alpha-hydroxy acid, preserving stereochemistry, and conformation of the backbone and the structure of the side chain. This substitution replaces the amide NH (a hydrogen bond donor) with an ester O (which is not a hydrogen bond donor) and the amide carbonyl (a strong hydrogen bond acceptor) with an ester carbonyl (a weaker hydrogen bond acceptor), thus perturbing folding energetics. Amide-to-ester perturbations were used to evaluate the thermodynamic contribution of each hydrogen bond in the PIN WW domain, a three-stranded beta-sheet protein. Our results reveal that removing a hydrogen bond donor destabilizes the native state more than weakening a hydrogen bond acceptor and that the degree of destabilization is strongly dependent on the location of the amide bond replaced. Hydrogen bonds near turns or at the ends of beta-strands are less influential than hydrogen bonds that are protected within a hydrophobic core. Beta-sheet destabilization caused by an amide-to-ester substitution cannot be directly related to hydrogen bond strength because of differences in the solvation and electrostatic interactions of amides and esters. We propose corrections for these differences to obtain approximate hydrogen bond strengths from destabilization energies. These corrections, however, do not alter the trends noted above, indicating that the destabilization energy of an amide-to-ester mutation is a good first-order approximation of the free energy of formation of a backbone amide hydrogen bond.     

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.     

How can a single second sphere amino acid substitution cause reduction midpoint potential changes of hundreds of millivolts?

The active site metal ion of superoxide dismutase (SOD) is reduced and reoxidized as it disproportionates superoxide to dioxygen and hydrogen peroxide. Thus, the reduction midpoint potential (Em) is a critical determinant of catalytic activity. In E. coli Fe-containing SOD (FeSOD), reduction of Fe3+ is accompanied by protonation of a coordinated OH-, to produce Fe2+ coordinated by H2O. The coordinated solvent's only contact with the protein beyond the active site is a conserved Gln residue. Mutation of this Gln to His or Glu resulted in elevation of the Em by 220 mV and more than 660 mV, respectively [Yikilmaz et al., Biochemistry 2006, 45, 1151-1161], despite the fact that overall protein structure was preserved, His is a chemically conservative replacement for Gln, and neutral Glu is isostructural and isoelectronic with Gln. Therefore, we have investigated several possible bases for the elevated Em's, including altered Fe electronic structure, altered active site electrostatics, altered H-bonding and altered redox-coupled proton transfer. Using EPR, MCD, and NMR spectroscopies, we find that the active site electronic structures of the two mutants resemble that of the WT enzyme, for both oxidation states, and Q69E-FeSOD's apparent deviation from WT-like Fe3+ coordination in the oxidized state can be explained by increased affinity for a small anion. Spontaneous coordination of an exogenous anion can only stabilize oxidized Q69E-Fe3+SOD and, therefore, cannot account for the increased Em of Q69E FeSOD. WT-like anion binding affinities and active site pK's indicate that His69 of Q69H-FeSOD is neutral in both oxidation states, like Gln69 of WT-FeSOD, whereas Glu69 appears to be neutral in the oxidized state but ionized in the reduced state of Q69E-FeSOD. A 1.1 A resolution crystal structure of Q69E-Fe2+SOD indicates that Glu69 accepts a strong H-bond from coordinated solvent in the reduced state, in contrast to the case in WT-FeSOD where Gln69 donates an H-bond. These data and DFT calculations lead to the proposal that the elevated Em of Q69E-FeSOD can be substantially explained by (1) relief from enforced H-bond donation in the reduced state, (2) Glu69's capacity to provide a proton for proton-coupled Fe3+ reduction, and (3) strong hydrogen bond acceptance in the reduced state, which stabilizes coordinated H2O. Our results thus support the hypothesis that the protein matrix can apply significant redox tuning via its influence over redox-coupled proton transfer and the energy associated with it.