An official website of the United States government
Flavoproteins are a ubiquitous class of redox proteins, enzymes, and photoreceptors that derive their versatility from the flavin cofactor—a prosthetic group that serves as the main locus of their spectral, photophysical, and (photo)chemical properties. It is thus common for computational modeling of flavoproteins to employ a hybrid approach that treats the flavin quantum mechanically and the remaining atoms classically. Such quantum mechanical/molecular mechanical (QM/MM) methods have proven powerful for studying flavoproteins so far, but users are often faced with a choice between treating the flavin electronic structure with ab initio wave function methods or using more approximate methods that allow for more extensive sampling of the protein dynamics. Herein, we present APEC-F 2.0, an automated QM/MM workflow that uses several open-source software packages to construct QM/MM models of flavoproteins. Exploiting the rigidity of flavin’s tricyclic isoalloxazine ring, the APEC approach iteratively optimizes flavin’s geometry in a static MM environment that represents a dynamic protein using a superposition of configurations generated from molecular dynamics. The automation of the code enables the systematic construction of QM/MM models using a common protocol and is suitable for comparing flavin’s spectral, electronic, and chemical properties in different redox, protonation, or excited states in a wide range of flavoproteins.
more »
« less
Vibrational spectroscopy of flavoproteins
The flavin cofactor performs many functions in the cell based on the ability of the isoalloxazine ring to undergo one- or two-electron reduction and form covalent adducts with reactants such as amino acids. In addition, the strong visible absorption of the cofactor is also the basis for flavin-dependent photoreceptors. Vibrational spectroscopy is uniquely suited to studying the mechanism of flavoproteins since the frequency of the vibrational modes is very sensitive to the electronic structure and environment of the isoalloxazine ring. This chapter describes the mechanistic information that can be gained using vibrational spectroscopy as well experimental challenges and approaches that are used to obtain and interpret the complex data contained in a vibrational spectrum.
more »
« less
- Award ID(s):
- 2042704
- PAR ID:
- 10315212
- Editor(s):
- Palfey, Bruce A.
- Date Published:
- Journal Name:
- Methods in enzymology
- Volume:
- 620
- ISSN:
- 0076-6879
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
FTIR spectroscopy accompanied by quantum chemical simulations can reveal important information about molecular structure and intermolecular interactions in the condensed phase. Simulations typically account for the solvent either through cluster quantum mechanical (QM) models, polarizable continuum models (PCM), or hybrid quantum mechanical/molecular mechanical (QM/MM) models. Recently, we studied the effect of aqueous solvent interactions on the vibrational frequencies of lumiflavin, a minimal flavin model, using cluster QM and PCM models. Those models successfully reproduced the relative frequencies of four prominent stretching modes of flavin’s isoalloxazine ring in the diagnostic 1450–1750 cm−1 range but poorly reproduced the relative band intensities. Here, we extend our studies on this system and account for solvation through a series of increasingly sophisticated models. Only by combining elements of QM clusters, QM/MM, and PCM approaches do we obtain an improved agreement with the experiment. The study sheds light more generally on factors that can impact the computed frequencies and intensities of IR bands in solution.more » « less
-
Dihydroorotate dehydrogenases (DHODs) catalyze the transfer of electrons between dihydroorotate and specific oxidant substrates. Class 1B DHODs (DHODBs) use NAD+ as the oxidant substrate and have a heterodimeric structure that incorporates two active sites, each with a flavin cofactor. One Fe2S2 center lies roughly equidistant between the flavin isoalloxazine rings. This arrangement allows for simultaneous association of reductant and oxidant substrates. Here we describe a series of experiments designed to reveal sequences and contingencies in DHODB chemistry. From these data it was concluded that the resting state of the enzyme is FAD•Fe2S2•FMN. Reduction by either NADH or DHO results in two electrons residing on the FMN cofactor that has a 47 mV higher reduction potential than the FAD. The FAD•Fe2S2•FMNH2 state accumulates with a bisemiquinone state that is an equilibrium accumulation formed from a partial transfer of one electron to the FAD. Pyrimidine reduction is reliant on the availability of the Cys135 proton, as the C135S variant slows orotate reduction by ∼40-fold. The rate of pyrimidine reduction is modulated by occupancy of the FAD site; NADH•FAD•Fe2S2•FMNH2•orotate complex can reduce the pyrimidine at 16 s–1, while NAD+•FAD•Fe2S2•FMNH2•orotate complex reduces the pyrimidine at 5.4 s–1 and the FAD•Fe2S2•FMNH2•orotate complex at 0.6 s–1. This set of effector states account for the apparent discrepancy in the slowest rate observed in transient state single turnover reactions with limiting NADH and the limiting rate observed in steady state.more » « less
-
null (Ed.)Cyanobacteriochromes (CBCRs) are promising optogenetic tools for their diverse absorption properties with a single compact cofactor-binding domain. We previously uncovered the ultrafast reversible photoswitching dynamics of a red/green photoreceptor AnPixJg2, which binds phycocyanobilin (PCB) that is unavailable in mammalian cells. Biliverdin (BV) is a mammalian cofactor with a similar structure to PCB but exhibits redder absorption. To improve the AnPixJg2 feasibility in mammalian applications, AnPixJg2_BV4 with only four mutations has been engineered to incorporate BV. Herein, we implemented femtosecond transient absorption (fs-TA) and ground state femtosecond stimulated Raman spectroscopy (GS-FSRS) to uncover transient electronic dynamics on molecular time scales and key structural motions responsible for the photoconversion of AnPixJg2_BV4 with PCB (Bpcb) and BV (Bbv) cofactors in comparison with the parent AnPixJg2 (Apcb). Bpcb adopts the same photoconversion scheme as Apcb, while BV4 mutations create a less bulky environment around the cofactor D ring that promotes a faster twist. The engineered Bbv employs a reversible clockwise/counterclockwise photoswitching that requires a two-step twist on ~5 and 35 picosecond (ps) time scales. The primary forward Pfr → Po transition displays equal amplitude weights between the two processes before reaching a conical intersection. In contrast, the primary reverse Po → Pfr transition shows a 2:1 weight ratio of the ~35 ps over 5 ps component, implying notable changes to the D-ring-twisting pathway. Moreover, we performed pre-resonance GS-FSRS and quantum calculations to identify the Bbv vibrational marker bands at ~659,797, and 1225 cm−1. These modes reveal a stronger H-bonding network around the BV cofactor A ring with BV4 mutations, corroborating the D-ring-dominant reversible photoswitching pathway in the excited state. Implementation of BV4 mutations in other PCB-binding GAF domains like AnPixJg4, AM1_1870g3, and NpF2164g5 could promote similar efficient reversible photoswitching for more directional bioimaging and optogenetic applications, and inspire other bioengineering advances.more » « less
-
null (Ed.)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.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1016/bs.mie.2019.03.011
