Flavin-based electron bifurcation allows enzymes to redistribute energy among electrons by coupling endergonic and exergonic electron transfer reactions. Diverse bifurcating enzymes employ a two-flavin electron transfer flavoprotein (ETF) that accepts hydride from NADH at a flavin (the so-called bifurcating FAD, Bf-FAD). The Bf-FAD passes one electron exergonically to a second flavin thereby assuming a reactive semiquinone state able to reduce ferredoxin or flavodoxin semiquinone. The flavin that accepts one electron and passes it on via exergonic electron transfer is known as the electron transfer FAD (ET-FAD) and is believed to correspond to the single FAD present in canonical ETFs, in domain II. The Bf-FAD is believed to be the one that is unique to bifurcating ETFs, bound between domains I and III. This very reasonable model has yet to be challenged experimentally. Herein we used site-directed mutagenesis to disrupt FAD binding to the presumed Bf site between domains I and III, in the Bf-ETF from Rhodopseudomonas palustris ( Rpa ETF). The resulting protein contained only 0.80 ± 0.05 FAD, plus 1.21 ± 0.04 bound AMP as in canonical ETFs. The flavin was not subject to reduction by NADH, confirming absence of Bf-FAD. The retained FAD displayed visible circular dichroism (CD) similar to that of the ET-FAD of Rpa ETF. Likewise, the mutant underwent two sequential one-electron reductions forming and then consuming anionic semiquinone, reproducing the reactivity of the ET-FAD. These data confirm that the retained FAD in domain II corresponds the ET-FAD. Quantum chemical calculations of the absorbance and CD spectra of each of WT Rpa ETF's two flavins reproduced the observed differences between their CD and absorbance signatures. The calculations for the flavin bound in domain II agreed better with the spectra of the ET-flavin, and those calculated based on the flavin between domains I and III agreed better with spectra of the Bf-flavin. Thus calculations independently confirm the locations of each flavin. We conclude that the site in domain II harbours the ET-FAD whereas the mutated site between domains I and III is the Bf-FAD site, confirming the accepted model by two different tests.
more »
« less
Understanding flavin electronic structure and spectra
Flavins have emerged as central to electron bifurcation, signaling, and countless enzymatic reactions. In bifurcation, two electrons acquired as a pair are separated in coupled transfers wherein the energy of both is concentrated on one of the two. This enables organisms to drive demanding reactions based on abundant low‐grade chemical fuel. To enable incorporation of this and other flavin capabilities into designed materials and devices, it is essential to understand fundamental principles of flavin electronic structure that make flavins so reactive and tunable by interactions with protein. Emerging computational tools can now replicate spectra of flavins and are gaining capacity to explain reactivity at atomistic resolution, based on electronic structures. Such fundamental understanding can moreover be transferrable to other chemical systems. A variety of computational innovations have been critical in reproducing experimental properties of flavins including their electronic spectra, vibrational signatures, and nuclear magnetic resonance (NMR) chemical shifts. A computational toolbox for understanding flavin reactivity moreover must be able to treat all five oxidation and protonation states, in addition to excited states that participate in flavoprotein's light‐driven reactions. Therefore, we compare emerging hybrid strategies and their successes in replicating effects of hydrogen bonding, the surrounding dielectric, and local electrostatics. These contribute to the protein's ability to modulate flavin reactivity, so we conclude with a survey of methods for incorporating the effects of the protein residues explicitly, as well as local dynamics. Computation is poised to elucidate the factors that affect a bound flavin's ability to mediate stunningly diverse reactions, and make life possible.
This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Electronic Structure Theory > Combined QM/MM Methods Theoretical and Physical Chemistry > Spectroscopy
- Award ID(s):
- 1808433
- NSF-PAR ID:
- 10363804
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- WIREs Computational Molecular Science
- Volume:
- 12
- Issue:
- 2
- ISSN:
- 1759-0876
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Flavins are central to countless enzymes but display different reactivities depending on their environments. This is understood to reflect modulation of the flavin electronic structure. To understand changes in orbital natures, energies, and correlation over the ring system, we begin by comparing seven flavin variants differing at C8, exploiting their different electronic spectra to validate quantum chemical calculations. Ground state calculations replicate a Hammett trend and reveal the significance of the flavin π-system. Comparison of higher-level theories establishes CC2 and ACD(2) as methods of choice for characterization of electronic transitions. Charge transfer character and electron correlation prove responsive to the identity of the substituent at C8. Indeed, bond length alternation analysis demonstrates extensive conjugation and delocalization from the C8 position throughout the ring system. Moreover, we succeed in replicating a particularly challenging UV/Vis spectrum by implementing hybrid QM/MM in explicit solvents. Our calculations reveal that the presence of nonbonding lone pairs correlates with the change in the UV/Vis spectrum observed when the 8-methyl is replaced by NH2, OH, or SH. Thus, our computations offer routes to understanding the spectra of flavins with different modifications. This is a first step toward understanding how the same is accomplished by different binding environments.more » « less
-
more » « less
Abstract Modern vibrational spectroscopy is more than just an analytical tool. Information about the electronic structure of a molecule, the strength of its bonds, and its conformational flexibility is encoded in the normal vibrational modes. On the other hand, normal vibrational modes are generally delocalized, which hinders the direct access to this information, attainable only via local vibration modes and associated local properties. Konkoli and Cremer provided an ingenious solution to this problem by deriving local vibrational modes from the fundamental normal modes, obtained in the harmonic approximation of the potential, via mass‐decoupled Euler–Lagrange equations. This review gives a general introduction into the local vibrational mode theory of Konkoli and Cremer, elucidating how this theory unifies earlier attempts to obtain easy to interpret chemical information from vibrational spectroscopy: (a) the local mode theory furnishes bond strength descriptors derived from force constant matrices with a physical basis, (b) provides the highly sought after extension of the Badger rule to polyatomic molecules, (c) and offers a simpler way to derive localized vibrations compared to the complex route via overtone spectroscopy. Successful applications are presented, including a new measure of bond strength, a new detailed analysis of infrared/Raman spectra, and the recent extension to periodic systems, opening a new avenue for the characterization of bonding in crystals. At the end of this review the LMODEA software is introduced, which performs the local mode analysis (with minimal computational costs) after a harmonic vibrational frequency calculation optionally using measured frequencies as additional input.
This article is categorized under:
Structure and Mechanism > Molecular Structures
Theoretical and Physical Chemistry > Spectroscopy
Software > Quantum Chemistry
Electronic Structure Theory > Ab Initio Electronic Structure Methods
-
more » « less
Abstract Interfaces, the boundary that separates two or more chemical compositions and/or phases of matter, alters basic chemical and physical properties including the thermodynamics of selectivity, transition states, and pathways of chemical reactions, nucleation events and phase growth, and kinetic barriers and mechanisms for mass transport and heat transport. While progress has been made in advancing more interface‐sensitive experimental approaches, their interpretation requires new theoretical methods and models that in turn can further elaborate on the microscopic physics that make interfacial chemistry so unique compared to the bulk phase. In this review, we describe some of the most recent theoretical efforts in modeling interfaces, and what has been learned about the transport and chemical transformations that occur at the air–liquid and solid–liquid interfaces.
This article is categorized under:
Structure and Mechanism > Reaction Mechanisms and Catalysis
Structure and Mechanism > Computational Materials Science
Software > Quantum Chemistry
Software > Simulation Methods
-
more » « less
Abstract Network analysis of infectious disease in wildlife can reveal traits or individuals critical to pathogen transmission and help inform disease management strategies. However, estimates of contact between animals are notoriously difficult to acquire. Researchers commonly use telemetry technologies to identify animal associations, but such data may have different sampling intervals and often captures a small subset of the population. The objectives of this study were to outline best practices for telemetry sampling in network studies of infectious disease by determining (a) the consequences of telemetry sampling on our ability to estimate network structure, (b) whether contact networks can be approximated using purely spatial contact definitions and (c) how wildlife spatial configurations may influence telemetry sampling requirements.
We simulated individual movement trajectories for wildlife populations using a home range‐like movement model, creating full location datasets and corresponding ‘complete’ networks. To mimic telemetry data, we created ‘sample’ networks by subsampling the population (10%–100% of individuals) with a range of sampling intervals (every minute to every 3 days). We varied the definition of contact for sample networks, using either spatiotemporal or spatial overlap, and varied the spatial configuration of populations (random, lattice or clustered). To compare complete and sample networks, we calculated seven network metrics important for disease transmission and assessed mean ranked correlation coefficients and percent error between complete and sample network metrics.
Telemetry sampling severely reduced our ability to calculate global node‐level network metrics, but had less impact on local and network‐level metrics. Even so, in populations with infrequent associations, high intensity telemetry sampling may still be necessary. Defining contact in terms of spatial overlap generally resulted in overly connected networks, but in some instances, could compensate for otherwise coarse telemetry data.
By synthesizing movement and disease ecology with computational approaches, we characterized trade‐offs important for using wildlife telemetry data beyond ecological studies of individual movement, and found that careful use of telemetry data has the potential to inform network models. Thus, with informed application of telemetry data, we can make significant advances in leveraging its use for a better understanding and management of wildlife infectious disease.
