Fentanyl and the Fluorinated Fentanyl Derivative NFEPP Elicit Distinct Hydrogen-Bond Dynamics of the Opioid Receptor

The development of safe therapeutics to manage pain is of central interest for biomedical applications. The fluorinated fentanyl derivative N-(3-fluoro-1-phenethylpiperidin-4-yl)-N-phenylpropionamide (NFEPP) is potentially a safer alternative to fentanyl because unlike fentanyl─which binds to the μ-opioid receptor (MOR) at both physiological and acidic pH─NFEPP might bind to the MOR only at acidic pH typical of inflamed tissue. Knowledge of the protonation-coupled dynamics of the receptor-drug interactions is thus required to understand the molecular mechanism by which receptor activation initiates cell signaling to silence pain. To this end, here we have carried out extensive atomistic simulations of the MOR in different protonation states, in the absence of opioid drugs, and in the presence of fentanyl vs NFEPP. We used graph-based analyses to characterize internal hydrogen-bond networks that could contribute to the activation of the MOR. We find that fentanyl and NFEPP prefer distinct binding poses and that, in their binding poses, fentanyl and NFEPP partake in distinct internal hydrogen-bond networks, leading to the cytoplasmic G-protein-binding region. Moreover, the protonation state of functionally important aspartic and histidine side chains impacts hydrogen-bond networks that extend throughout the receptor, such that the ligand-bound MOR presents at its cytoplasmic G-protein-binding side, a hydrogen-bonding environment where dynamics depend on whether fentanyl or NFEPP is bound, and on the protonation state of specific MOR groups. The exquisite sensitivity of the internal protein-water hydrogen-bond network to the protonation state and to details of the drug binding could enable the MOR to elicit distinct pH- and opioid-dependent responses at its cytoplasmic G-protein-binding site.

QM/MM Modeling of the Flavin Functionalization in the RutA Monooxygenase

Oxygenase activity of the flavin-dependent enzyme RutA is commonly associated with the formation of flavin-oxygen adducts in the enzyme active site. We report the results of quantum mechanics/molecular mechanics (QM/MM) modeling of possible reaction pathways initiated by various triplet state complexes of the molecular oxygen with the reduced flavin mononucleotide (FMN) formed in the protein cavities. According to the calculation results, these triplet-state flavin-oxygen complexes can be located at both re-side and si-side of the isoalloxazine ring of flavin. In both cases, the dioxygen moiety is activated by electron transfer from FMN, stimulating the attack of the arising reactive oxygen species at the C4a, N5, C6, and C8 positions in the isoalloxazine ring after the switch to the singlet state potential energy surface. The reaction pathways lead to the C(4a)-peroxide, N(5)-oxide, or C(6)-hydroperoxide covalent adducts or directly to the oxidized flavin, depending on the initial position of the oxygen molecule in the protein cavities.

Quantum-mechanical insights into the anisotropic response of the cryptochrome radical pair to a weak magnetic field

Cryptochrome photoreceptors contain a photochemically generated radical pair, which is thought to mediate sensing of the geomagnetic field direction in many living organisms. To gain insight into the response of the cryptochrome to a weak magnetic field, we have studied the quantum-mechanical hyperfine spin states of the radical pair. We identify quantum states responsible for the precise detection of the magnetic field direction, taking into account the strongly axial hyperfine interactions of each radical in the radical pair. The contribution of these states to the formation of the cryptochrome signaling state sharply increases when the magnetic field becomes orthogonal to the hyperfine axis of either radical. Due to such a response, the radical pair may be able to detect the particular field direction normal to the plane containing the hyperfine axes of the radicals.

Spectroscopic and Computational Observation of Glutamine Tautomerization in the Blue Light Sensing Using Flavin Domain Photoreaction

Blue light sensing using flavin (BLUF) domains constitute a family of flavin-binding photoreceptors of bacteria and eukaryotic algae. BLUF photoactivation proceeds via a light-driven hydrogen-bond switch among flavin adenine dinucleotide (FAD) and glutamine and tyrosine side chains, whereby FAD undergoes electron and proton transfer with tyrosine and is subsequently re-oxidized by a hydrogen back-shuttle in picoseconds, constituting an important model system to understand proton-coupled electron transfer in biology. The specific structure of the hydrogen-bond patterns and the prevalence of glutamine tautomeric states in dark-adapted (DA) and light-activated (LA) states have remained controversial. Here, we present a combined femtosecond stimulated Raman spectroscopy (FSRS), computational chemistry, and site-selective isotope labeling Fourier-transform infrared spectroscopy (FTIR) study of the Slr1694 BLUF domain. FSRS showed distinct vibrational bands from the FADS₁ singlet excited state. We observed small but significant shifts in the excited-state vibrational frequency patterns of the DA and LA states, indicating that these frequencies constitute a sensitive probe for the hydrogen-bond arrangement around FAD. Excited-state model calculations utilizing four different realizations of hydrogen bond patterns and glutamine tautomeric states were consistent with a BLUF reaction model that involved glutamine tautomerization to imidic acid, accompanied by a rotation of its side chain. A combined FTIR and double-isotope labeling study, with ¹³C labeling of FAD and ¹⁵N labeling of glutamine, identified the glutamine imidic acid C═N stretch vibration in the LA state and the Gln C═O in the DA state. Hence, our study provides support for glutamine tautomerization and side-chain rotation in the BLUF photoreaction.

Insight into the structural dynamics of light sensitive proteins from time-resolved crystallography and quantum chemical calculations

The structural dynamics underlying molecular mechanisms of light-sensitive proteins can be studied by a variety of experimental and computational biophysical techniques. Here we review recent progress in combining time-resolved crystallography at X-ray free electron lasers and quantum chemical calculations to study structural changes in photoenzymes, photosynthetic proteins, photoreceptors, and photoswitchable fluorescent proteins following photoexcitation.

Rational Control of Off-State Heterogeneity in a Photoswitchable Fluorescent Protein Provides Switching Contrast Enhancement

Reversibly photoswitchable fluorescent proteins are essential markers for advanced biological imaging, and optimization of their photophysical properties underlies improved performance and novel applications. Here we establish a link between photoswitching contrast, one of the key parameters that dictate the achievable resolution in nanoscopy applications, and chromophore conformation in the non-fluorescent state of rsEGFP2, a widely employed label in REversible Saturable OpticaL Fluorescence Transitions (RESOLFT) microscopy. Upon illumination, the cis chromophore of rsEGFP2 isomerizes to two distinct off-state conformations, trans1 and trans2, located on either side of the V151 side chain. Reducing or enlarging the side chain at this position (V151A and V151L variants) leads to single off-state conformations that exhibit higher and lower switching contrast, respectively, compared to the rsEGFP2 parent. The combination of structural information obtained by serial femtosecond crystallography with high-level quantum chemical calculations and with spectroscopic and photophysical data determined in vitro suggests that the changes in switching contrast arise from blue- and red-shifts of the absorption bands associated to trans1 and trans2, respectively. Thus, due to elimination of trans2, the V151A variants of rsEGFP2 and its superfolding variant rsFolder2 display a more than two-fold higher switching contrast than their respective parent proteins, both in vitro and in E. coli cells. The application of the rsFolder2-V151A variant is demonstrated in RESOLFT nanoscopy. Our study rationalizes the connection between structural and photophysical chromophore properties and suggests a means to rationally improve fluorescent proteins for nanoscopy applications.

Protonation States of Molecular Groups in the Chromophore-Binding Site Modulate Properties of the Reversibly Switchable Fluorescent Protein rsEGFP2

The role of protonation states of the chromophore and its neighboring amino acid side chains of the reversibly switching fluorescent protein rsEGFP2 upon photoswitching is characterized by molecular modeling methods. Numerous conformations of the chromophore-binding site in computationally derived model systems are obtained using the quantum chemistry and QM/MM approaches. Excitation energies are computed using the extended multiconfigurational quasidegenerate perturbation theory (XMCQDPT2). The obtained structures and absorption spectra allow us to provide an interpretation of the observed structural and spectral properties of rsEGFP2 in the active ON and inactive OFF states. The results demonstrate that in addition to the dominating anionic and neutral forms of the chromophore, the cationic and zwitterionic forms may participate in the photoswitching of rsEGFP2. Conformations and protonation forms of the Glu223 and His149 side chains in the chromophore-binding site play an essential role in stabilizing specific protonation forms of the chromophore.

CHARMM Force-Field Parameters for Morphine, Heroin, and Oliceridine, and Conformational Dynamics of Opioid Drugs

Opioid drug binding to specialized G protein-coupled receptors (GPCRs) can lead to analgesia upon activation via downstream G_i protein signaling and to severe side effects via activation of the β-arrestin signaling pathway. Knowledge of how different opioid drugs interact with receptors is essential, as it can inform and guide the design of safer therapeutics. We performed quantum and classical mechanical computations to explore the potential energy landscape of four opioid drugs: morphine and its derivatives heroin and fentanyl and for the unrelated oliceridine. From potential energy profiles for bond twists and from interactions between opioids and water, we derived a set of force-field parameters that allow a good description of structural properties and intermolecular interactions of the opioids. Potential of mean force profiles computed from molecular dynamics simulations indicate that fentanyl and oliceridine have complex energy landscapes with relatively small energy penalties, suggesting that interactions with the receptor could select different binding poses of the drugs.

Mechanism and dynamics of fatty acid photodecarboxylase

Fatty acid photodecarboxylase (FAP) is a photoenzyme with potential green chemistry applications. By combining static, time-resolved, and cryotrapping spectroscopy and crystallography as well as computation, we characterized Chlorella variabilis FAP reaction intermediates on time scales from subpicoseconds to milliseconds. High-resolution crystal structures from synchrotron and free electron laser x-ray sources highlighted an unusual bent shape of the oxidized flavin chromophore. We demonstrate that decarboxylation occurs directly upon reduction of the excited flavin by the fatty acid substrate. Along with flavin reoxidation by the alkyl radical intermediate, a major fraction of the cleaved carbon dioxide unexpectedly transformed in 100 nanoseconds, most likely into bicarbonate. This reaction is orders of magnitude faster than in solution. Two strictly conserved residues, R451 and C432, are essential for substrate stabilization and functional charge transfer.

Frontiers in Multiscale Modeling of Photoreceptor Proteins

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

Four resonance structures elucidate double-bond isomerisation of a biological chromophore

Four resonance structures determining the electronic structure of the chromophore’s ground and first excited states. Changing the relative energies of the structures by hydrogen-bonding interactions tunes all chromophore’s photochemical properties.

Three-dimensional view of ultrafast dynamics in photoexcited bacteriorhodopsin

Bacteriorhodopsin (bR) is a light-driven proton pump. The primary photochemical event upon light absorption is isomerization of the retinal chromophore. Here we used time-resolved crystallography at an X-ray free-electron laser to follow the structural changes in multiphoton-excited bR from 250 femtoseconds to 10 picoseconds. Quantum chemistry and ultrafast spectroscopy were used to identify a sequential two-photon absorption process, leading to excitation of a tryptophan residue flanking the retinal chromophore, as a first manifestation of multiphoton effects. We resolve distinct stages in the structural dynamics of the all-trans retinal in photoexcited bR to a highly twisted 13-cis conformation. Other active site sub-picosecond rearrangements include correlated vibrational motions of the electronically excited retinal chromophore, the surrounding amino acids and water molecules as well as their hydrogen bonding network. These results show that this extended photo-active network forms an electronically and vibrationally coupled system in bR, and most likely in all retinal proteins.

Bacteriorhodopsin (bR) is a light-driven proton pump. Here the authors combine time-resolved crystallography at a free-electron laser, ultrafast spectroscopy and quantum chemistry to study the structural changes following multiphoton photoexcitation of bR and find that they occur within 300 fs not only in the light-absorbing chromophore but also in the surrounding protein.

Spiers Memorial Lecture. Introductory lecture: the impact of structure on photoinduced processes in nucleic acids and proteins

Light is an important environmental variable and most organisms have evolved means to sense, exploit or avoid it and to repair detrimental effects on their genome. In general, light absorption is the task of specific chromophores, however other biomolecules such as oligonucleotides also do so which can result in undesired outcomes such as mutations and cancer. Given the biological importance of light-induced processes and applications for imaging, optogenetics, photodynamic therapy or photovoltaics, there is a great interest in understanding the detailed molecular mechanisms of photoinduced processes in proteins and nucleic acids. The processes are typically characterized by time-resolved spectroscopic approaches or computation, inferring structural information on transient species from stable ground state structures. Recently, however, structure determination of excited states or other short-lived species has become possible with the advent of X-ray free-electron lasers. This review gives an overview of the impact of structure on the understanding of photoinduced processes in macromolecules, focusing on systems presented at this Faraday Discussion meeting.

Indications of 5' to 3' Interbase Electron Transfer as the First Step of Pyrimidine Dimer Formation Probed by a Dinucleotide Analog

Pyrimidine dimers are the most common DNA lesions generated under UV radiation. To reveal the molecular mechanisms behind their formation, it is of significance to reveal the roles of each pyrimidine residue. We thus replaced the 5'-pyrimidine residue with a photochemically inert xylene moiety (X). The electron-rich X can be readily oxidized but not reduced, defining the direction of interbase electron transfer (ET). Irradiation of the XpT dinucleotide under 254 nm UV light generates two major photoproducts: a pyrimidine (6-4) pyrimidone analog (6-4PP) and an analog of the so-called spore photoproduct (SP). Both products are formed by reaction at C4=O of the photo-excited 3'-thymidine (T), which indicates that excitation of a single "driver" residue is sufficient to trigger pyrimidine dimerization. Our quantum-chemical calculations demonstrated that photo-excited 3'-T accepts an electron from 5'-X. The resulting charge-separated radical pair lowers its energy upon formation of interbase covalent bonds, eventually yielding 6-4PP and SP.

Role of the Molecular Environment in Flavoprotein Color and Redox Tuning: QM Cluster versus QM/MM Modeling

We investigate the origin of the excitation energy shifts induced by the apoprotein in the active site of the bacterial photoreceptor BLUF (Blue Light sensor Using Flavin adenine dinucleotide). In order to compute the vertical excitation energies of three low-lying electronic states, including two π-π* states of flavin (S1 and S2) and a π-π* tyrosine-flavin electron-transfer state (ET), with respect to the energy of the closed-shell ground state (S0), we prepared alternative quantum mechanical (QM) cluster and quantum mechanics/molecular mechanics (QM/MM) models. We found that the excitation energies computed with both types of models correlate with the magnitude of the charge transfer character of the excitation. Accordingly, we conclude that the small charge transfer character of the light absorbing S0-S1 transition and the substantial charge transfer character of the nonabsorbing but redox active S0-ET transition explain the small color changes but substantial redox tuning in BLUF and also in other flavoproteins. Further analysis showed that redox tuning is governed by the electrostatic interaction in the QM/MM model and transfer of charge between the active site and its environment in the QM cluster. Moreover, the wave function polarization of the QM subsystem by the MM subsystem influences the magnitude of the charge transfer, resulting in the QM/MM and QM excitation energies that are not entirely consistent.

Onset of the Electronic Absorption Spectra of Isolated and π-Stacked Oligomers of 5,6-Dihydroxyindole: An Ab Initio Study of the Building Blocks of Eumelanin

Eumelanin is a naturally occurring skin pigment which is responsible for developing a suntan. The complex structure of eumelanin consists of π-stacked oligomers of various indole derivatives, such as the monomeric building block 5,6-dihydroxyindole (DHI). In this work, we present an ab initio wave-function study of the absorption behavior of DHI oligomers and of doubly and triply π-stacked species of these oligomers. We have simulated the onset of the electronic absorption spectra by employing the MP2 and the linear-response CC2 methods. Our results demonstrate the effect of an increasing degree of oligomerization of DHI and of an increasing degree of π-stacking of DHI oligomers on the onset of the absorption spectra and on the degree of red-shift toward the visible region of the spectrum. We find that π-stacking of DHI and its oligomers substantially red-shifts the onset of the absorption spectra. Our results also suggest that the optical properties of biological eumelanin cannot be simulated by considering the DHI building blocks alone, but instead the building blocks indole-semiquinone and indole-quinone have to be considered as well. This study contributes to advancing the understanding of the complex photophysics of the eumelanin biopolymer.

Single Hydrogen Bond Donation from Flavin N5 to Proximal Asparagine Ensures FAD Reduction in DNA Photolyase

The spread of the absorbance of the stable FADH(•) radical (300-700 nm) allows CPD photolyase to highly efficiently form FADH(-), making it functional for DNA repair. In this study, FTIR spectroscopy detected a strong hydrogen bond, from FAD N5-H to the carbonyl group of the Asn378 side chain, that is modulated by the redox state of FAD. The observed characteristic frequency shifts were reproduced in quantum-mechanical models of the flavin binding site, which were then employed to elucidate redox tuning governed by Asn378. We demonstrate that enhanced hydrogen bonding of the Asn378 side chain with the FADH(•) radical increases thermodynamic stabilization of the radical state, and further ensures kinetic stabilization and accumulation of the fully reduced FADH(-) state.

Evidence for Tautomerisation of Glutamine in BLUF Blue Light Receptors by Vibrational Spectroscopy and Computational Chemistry

BLUF (blue light sensor using flavin) domains regulate the activity of various enzymatic effector domains in bacteria and euglenids. BLUF features a unique photoactivation through restructuring of the hydrogen-bonding network as opposed to a redox reaction or an isomerization of the chromophore. A conserved glutamine residue close to the flavin chromophore plays a central role in the light response, but the underlying modification is still unclear. We labelled this glutamine with (15)N in two representative BLUF domains and performed time-resolved infrared double difference spectroscopy. The assignment of the signals was conducted by extensive quantum chemical calculations on large models with 187 atoms reproducing the UV-vis and infrared signatures of BLUF photoactivation. In the dark state, the comparatively low frequency of 1,667 cm(-1) is assigned to the glutamine C=O accepting a hydrogen bond from tyrosine. In the light state, the signature of a tautomerised glutamine was extracted with the C=N stretch at ~1,691 cm(-1) exhibiting the characteristic strong downshift by (15)N labelling. Moreover, an indirect isotope effect on the flavin C4=O stretch was found. We conclude that photoactivation of the BLUF receptor does not only involve a rearrangement of hydrogen bonds but includes a change in covalent bonds of the protein.

Analysis of the Primary Photocycle Reactions Occurring in the Light, Oxygen, and Voltage Blue-Light Receptor by Multiconfigurational Quantum-Chemical Methods

The photocycle reactions occurring between the flavin mononucleotide cofactor and the reactive cysteine residue in the blue-light photoreceptor domain light, oxygen, and voltage (LOV) were modeled for a system consisting of lumiflavin and thiomethanol. The electronic structure and energies of the reactive species were estimated using the CASSCF and MCQDPT2 quantum-chemical methods. The reaction pathway for the S-C4a covalent adduct formation in the triplet state was determined. Concerted electron and proton transfer from the thiol to the flavin in the triplet electronic state results in a biradical complex that is, however, unstable because its structure corresponds to a triplet-singlet crossing. The covalent adduct dissociation in the ground electronic state is a reverse of the photoreaction proceeding via a single energy barrier for hydrogen transfer. Thus, both photo- and dark reactions were found to be single-step chemical transformations occurring without stable intermediates. The photoreaction yielding the S-C4a covalent adduct is an intrinsic property of the isoalloxazine-thiol complex in the specific geometry arranged by the protein in LOV. The S-C4a covalent adduct between lumiflavin and thiomethanol is rather stable implying that in LOV its dissociation is facilitated by the protein.

Separation of photo-induced radical pair in cryptochrome to a functionally critical distance

Cryptochrome is a blue light receptor that acts as a sensor for the geomagnetic field and assists many animals in long-range navigation. The magnetoreceptor function arises from light-induced formation of a radical pair through electron transfer between a flavin cofactor (FAD) and a triad of tryptophan residues. Here, this electron transfer is investigated by quantum chemical and classical molecular dynamics calculations. The results reveal how sequential electron transfer, assisted by rearrangement of polar side groups in the cryptochrome interior, can yield a FAD-Trp radical pair state with the FAD and Trp partners separated beyond a critical distance. The large radical pair separation reached establishes cryptochrome's sensitivity to the geomagnetic field through weakening of distance-dependent exchange and dipole-dipole interactions. It is estimated that the key secondary electron transfer step can overcome in speed both recombination (electron back-transfer) and proton transfer involving the radical pair reached after primary electron transfer.

Computational spectroscopy, dynamics, and photochemistry of photosensory flavoproteins

Extensive interest in photosensory proteins stimulated computational studies of flavins and flavoproteins in the past decade. This review is dedicated to the three central topics of these studies: calculations of flavin UV-visible and IR spectra, simulated dynamics of photoreceptor proteins, and flavin photochemistry. Accordingly, this chapter is divided into three parts; each part describes corresponding computational protocols, summarizes computational results, and discusses the emerging mechanistic picture.

Photoinduced electron transfer facilitates tautomerization of the conserved signaling glutamine side chain in BLUF protein light sensors

The BLUF domain (sensor of blue light using flavin adenine dinucleotide) from a bacterial photoreceptor protein AppA undergoes a cascade of chemical transformations, including hydrogen bond rearrangements around the flavin adenine dinucleotide (FAD) chromophore, in response to light illumination. These transformations are initiated by photoinduced electron and proton transfer from a tyrosine residue to the photoexcited flavin which is assisted by a glutamine residue. According to the recent studies, the proton-coupled electron transfer leads to formation of a radical-pair intermediate Tyr•···FADH• and a tautomeric EE form of glutamine in the ground electronic state. This intermediate is a precursor of the light-induced state of the BLUF photoreceptor implicated in biological signaling. In order to describe evolution of the radical pair, we computed reaction pathways on the ground state potential energy surface employing quantum-chemical calculations in the DFT PBE0/cc-pVDZ approximation for a molecular cluster mimicking the chromophore containing pocket of the AppA BLUF protein. We found a minimum-energy pathway comprised of the following consecutive reaction steps: (1) rotation of the imidic group of the EE glutamine side chain around the Cγ-Cδ bond; (2) flip of the OεH group and formation of the ZE form of the glutamine side chain; and (3) biradical recombination via coupled proton and electron transfer, leading to the ZZ form of the glutamine side chain. The potential-energy barriers for stages 1-3 do not exceed 9 kcal/mol. Energy barrier 3 describing the ZE to ZZ glutamine tautomerization is significantly smaller in the BLUF model than in isolated glutamine, since tautomerization in BLUF is facilitated by electron transfer and radical recombination. Thus, our study shows that tautomerization of the conserved glutamine is coupled to the light-induced electron transfer process in BLUF and, thus, is a viable candidate for the photoactivation mechanism which at present is very much debated.

Decrypting cryptochrome: revealing the molecular identity of the photoactivation reaction

Migrating birds fly thousands of miles or more, often without visual cues and in treacherous winds, yet keep direction. They employ for this purpose, apparently as a powerful navigational tool, the photoreceptor protein cryptochrome to sense the geomagnetic field. The unique biological function of cryptochrome supposedly arises from a photoactivation reaction involving radical pair formation through electron transfer. Radical pairs, indeed, can act as a magnetic compass; however, the cryptochrome photoreaction pathway is not fully resolved yet. To reveal this pathway and underlying photochemical mechanisms, we carried out a combination of quantum chemical calculations and molecular dynamics simulations on plant ( Arabidopsis thaliana ) cryptochrome. The results demonstrate that after photoexcitation a radical pair forms, becomes stabilized through proton transfer, and decays back to the protein's resting state on time scales allowing the protein, in principle, to act as a radical pair-based magnetic sensor. We briefly relate our findings on A. thaliana cryptochrome to photoreaction pathways in animal cryptochromes.

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.

Coupling between the BLUF and EAL domains in the blue light-regulated phosphodiesterase BlrP1

The first biochemical and structural characterization of the full-length active photoreceptor BlrP1 from Klebsiella pneumoniae was recently reported by Barends et al. [Nature 459:1015-1018, (2009)]. The light-regulated catalytic function of its C-terminal c-di-guanosine monophosphate phosphodiesterase, the EAL (Glu-Ala-Leu) domain, is activated by the N-terminal sensor of blue light using the flavin adenine dinucleotide (BLUF) domain. We performed molecular dynamics simulations on the dimeric BlrP1 protein in order to examine the coupling regions that are presumably involved in transmitting light-induced structural changes which occur in the BLUF domain to the EAL domain. According to the results of simulations and an analysis of the hydrogen bonding between the respective polypeptide chains, the region containing the site on the α3α4 loop of BLUF is responsible for communication between the photosensing and catalytic domains in the dimeric BlrP1 protein.

Electronic structure of (6-4) DNA photoproduct repair involving a non-oxetane pathway

Mutagenic pyrimidine-pyrimidone (6-4) photoproducts are one of the main DNA lesions induced by solar UV radiation. These lesions can be photoreversed by (6-4) photolyases. The originally published repair mechanism involves rearrangement of the lesion into an oxetane intermediate upon binding to the (6-4) photolyase, followed by light-induced electron transfer from the reduced flavin cofactor. In a recent crystallographic study on a (6-4) photoproduct complexed with (6-4) photolyase from Drosophila melanogaster no oxetane was observed, raising the possibility of a non-oxetane repair mechanism. Using quantum-chemical calculations we find that in addition to repair via an oxetane, a direct transfer of the hydroxyl group results in reversal of the radical anion (6-4) photoproduct. In both mechanisms, the transition states have high energies and correspond to avoided crossings of the ground and excited electronic states. To study whether the repair can proceed via these state crossings, the excited-state potential energy curves were computed. The radical excitation energies and accessibility of the nonadiabatic repair path were found to depend on hydrogen bonds and the protonation state of the lesion. On the basis of the energy calculations, a nonadiabatic repair of the excited (6-4) lesion radical anion via hydroxyl transfer is probable. This repair mechanism is in line with the recent structural data on the (6-4) photolyase from D. melanogaster .

Primary reactions of the LOV2 domain of phototropin studied with ultrafast mid-infrared spectroscopy and quantum chemistry

Phototropins, major blue-light receptors in plants, are sensitive to blue light through a pair of flavin mononucleotide (FMN)-binding light oxygen and voltage (LOV) domains, LOV1 and LOV2. LOV2 undergoes a photocycle involving light-driven covalent adduct formation between a conserved cysteine and the FMN C(4a) atom. Here, the primary reactions of Avena sativa phototropin 1 LOV2 (AsLOV2) were studied using ultrafast mid-infrared spectroscopy and quantum chemistry. The singlet excited state (S1) evolves into the triplet state (T1) with a lifetime of 1.5 ns at a yield of approximately 50%. The infrared signature of S1 is characterized by absorption bands at 1657 cm(-1), 1495-1415 cm(-1), and 1375 cm(-1). The T1 state shows infrared bands at 1657 cm(-1), 1645 cm(-1), 1491-1438 cm(-1), and 1390 cm(-1). For both electronic states, these bands are assigned principally to C=O, C=N, C-C, and C-N stretch modes. The overall downshifting of C=O and C=N bond stretch modes is consistent with an overall bond-order decrease of the conjugated isoalloxazine system upon a pi-pi* transition. The configuration interaction singles (CIS) method was used to calculate the vibrational spectra of the S1 and T1 excited pipi* states, as well as respective electronic energies, structural parameters, electronic dipole moments, and intrinsic force constants. The harmonic frequencies of S1 and T1, as calculated by the CIS method, are in satisfactory agreement with the evident band positions and intensities. On the other hand, CIS calculations of a T1 cation that was protonated at the N(5) site did not reproduce the experimental FMN T1 spectrum. We conclude that the FMN T1 state remains nonprotonated on a nanosecond timescale, which rules out an ionic mechanism for covalent adduct formation involving cysteine-N(5) proton transfer on this timescale. Finally, we observed a heterogeneous population of singly and doubly H-bonded FMN C(4)=O conformers in the dark state, with stretch frequencies at 1714 cm(-1) and 1694 cm(-1), respectively.

Crystal Structures of the AppA BLUF Domain Photoreceptor Provide Insights into Blue Light-mediated Signal Transduction

Proteins containing a sensor of blue light using FAD (BLUF) domain control diverse cellular processes, such as gene expression, nucleotide metabolism and motility, by relaying blue light signals to distinct output units. Despite its crucial and widespread functions, the mechanism of BLUF signal transduction has remained elusive. We determined crystal structures of the dark-adapted state and of a photo-excited, red-shifted photocycle intermediate of the BLUF unit of AppA, a purple bacterial photoreceptor involved in the light-dependent regulation of photosynthesis gene expression. In contrast to a recently published crystal structure of the AppA BLUF domain determined in the presence of detergent molecules, our structural model of the dark state corresponds well to those reported for the BLUF domains of Tll0078 and BlrB. This establishes that a highly conserved methionine (Met106 in AppA) is next to the active site glutamine (Gln63 in AppA), which is of relevance for the latter's orientation in the dark state and for the mechanism of the photoreaction. The comparison of the dark-adapted and photointermediate state structures shows light-induced conformational alterations, which suggest a path for signal propagation. In particular, we observe a significant movement of the Met106 side-chain. Met106 thereby changes its mode of interaction with Gln63, which supports a light-dependent rotation of the latter. In view of other BLUF structures available, our data further suggest that the hydrogen bond between Asn45 and the backbone carbonyl of His105 breaks upon illumination. The ensuing extensive structural rearrangement of beta-strand 5 is predicted to involve a flip of Met106 out of the flavin-binding pocket and Trp104 moving in to fill the void. We propose that the blue light signal is transmitted towards the surface of the BLUF domain via His44, which serves as a reporter of active site changes.

Structure of a bacterial BLUF photoreceptor: insights into blue light-mediated signal transduction

Light is an essential environmental factor, and many species have evolved the capability to respond to it. Blue light is perceived through three flavin-containing photoreceptor families: cryptochromes, light-oxygen-voltage, and BLUF (sensor of blue light using flavin adenine dinucleotide, FAD) domain proteins. BLUF domains are present in various proteins from Bacteria and lower Eukarya. They are fully modular and can relay signals to structurally and functionally diverse output units, most of which are implicated in nucleotide metabolism. We present the high resolution crystal structure of the dark resting state of BlrB, a short BLUF domain-containing protein from Rhodobacter sphaeroides. The structure reveals a previously uncharacterized FAD-binding fold. Along with other lines of evidence, it suggests mechanistic aspects for the photocycle that is characterized by a red-shifted absorbance of the flavin. The isoalloxazine ring of FAD binds in a cleft between two helices, whereas the adenine ring points into the solvent. We propose that the adenine ring serves as a hook mediating the interaction with its effector/output domain. The structure suggests a unique photochemical signaling switch in which the absorption of light induces a structural change in the rim surrounding the hook, thereby changing the protein interface between BLUF and the output domain.