New insights into the mechanism of electron transfer within flavohemoglobins: tunnelling pathways, packing density, thermodynamic and kinetic analyses

Ordered Motions in the Nitric-Oxide Dioxygenase Mechanism of Flavohemoglobin and Assorted Globins with Tightly Coupled Reductases

Nitric-oxide dioxygenases (NODs) activate and combine O₂ with NO to form nitrate. A variety of oxygen-binding hemoglobins with associated partner reductases or electron donors function as enzymatic NODs. Kinetic and structural investigations of the archetypal two-domain microbial flavohemoglobin-NOD have illuminated an allosteric mechanism that employs selective tunnels for O₂ and NO, gates for NO and nitrate, transient O₂ association with ferric heme, and an O₂ and NO-triggered, ferric heme spin crossover-driven, motion-controlled, and dipole-regulated electron-transfer switch. The proposed mechanism facilitates radical-radical coupling of ferric-superoxide with NO to form nitrate while preventing suicidal ferrous-NO formation. Diverse globins display the structural and functional motifs necessary for a similar allosteric NOD mechanism. In silico docking simulations reveal monomeric erythrocyte hemoglobin alpha-chain and beta-chain intrinsically matched and tightly coupled with NADH-cytochrome b₅ oxidoreductase and NADPH-cytochrome P450 oxidoreductase, respectively, forming membrane-bound flavohemoglobin-like mammalian NODs. The neuroprotective neuroglobin manifests a potential NOD role in a close-fitting ternary complex with membrane-bound NADH-cytochrome b₅ oxidoreductase and cytochrome b₅. Cytoglobin interfaces weakly with cytochrome b₅ for O₂ and NO-regulated electron-transfer and coupled NOD activity. The mechanistic model also provides insight into the evolution of O₂ binding cooperativity in hemoglobin and a basis for the discovery of allosteric NOD inhibitors.

Mechanistic Insights on Heme-to-Heme Transmembrane Electron Transfer Within NADPH Oxydases From Atomistic Simulations

NOX5 is a member of the NADPH oxidase family which is dedicated to the production of reactive oxygen species. The molecular mechanisms governing transmembrane electron transfer (ET) that permits to shuttle electrons over the biological membrane have remained elusive for a long time. Using computer simulations, we report conformational dynamics of NOX5 embedded within a realistic membrane environment. We assess the stability of the protein within the membrane and monitor the existence of cavities that could accommodate dioxygen molecules. We investigate the heme-to-heme electron transfer. We find a reaction free energy of a few tenths of eV (ca. −0.3 eV) and a reorganization free energy of around 1.1 eV (0.8 eV after including electrostatic induction corrections). The former indicates thermodynamically favorable ET, while the latter falls in the expected values for transmembrane inter-heme ET. We estimate the electronic coupling to fall in the range of the μeV. We identify electron tunneling pathways showing that not only the W378 residue is playing a central role, but also F348. Finally, we reveal the existence of two connected O₂₋binding pockets near the outer heme with fast exchange between the two sites on the nanosecond timescale. We show that when the terminal heme is reduced, O₂ binds closer to it, affording a more efficient tunneling pathway than when the terminal heme is oxidized, thereby providing an efficient mechanism to catalyze superoxide production in the final step. Overall, our study reveals some key molecular mechanisms permitting reactive oxygen species production by NOX5 and paves the road for further investigation of ET processes in the wide family of NADPH oxidases by computer simulations.

Allostery in the nitric oxide dioxygenase mechanism of flavohemoglobin

The substrates O₂ and NO cooperatively activate the NO dioxygenase function of Escherichia coli flavohemoglobin. Steady-state and transient kinetic measurements support a structure-based mechanistic model in which O₂ and NO movements and conserved amino acids at the E11, G8, E2, E7, B10, and F7 positions within the globin domain control activation. In the cooperative and allosteric mechanism, O₂ migrates to the catalytic heme site via a long hydrophobic tunnel and displaces LeuE11 away from the ferric iron, which forces open a short tunnel to the catalytic site gated by the ValG8/IleE15 pair and LeuE11. NO permeates this tunnel and leverages upon the gating side chains triggering the CD loop to furl, which moves the E and F-helices and switches an electron transfer gate formed by LysF7, GlnE7, and water. This allows FADH₂ to reduce the ferric iron, which forms the stable ferric–superoxide–TyrB10/GlnE7 complex. This complex reacts with internalized NO with a bimolecular rate constant of 10¹⁰ M⁻¹ s⁻¹ forming nitrate, which migrates to the CD loop and unfurls the spring-like structure. To restart the cycle, LeuE11 toggles back to the ferric iron. Actuating electron transfer with O₂ and NO movements averts irreversible NO poisoning and reductive inactivation of the enzyme. Together, structure snapshots and kinetic constants provide glimpses of intermediate conformational states, time scales for motion, and associated energies.

Flavohaemoglobin: the pre-eminent nitric oxide-detoxifying machine of microorganisms

Flavohaemoglobins were first described in yeast as early as the 1970s but their functions were unclear. The surge in interest in nitric oxide biology and both serendipitous and hypothesis-driven discoveries in bacterial systems have transformed our understanding of this unusual two-domain globin into a comprehensive, yet undoubtedly incomplete, appreciation of its pre-eminent role in nitric oxide detoxification. Here, I focus on research on the flavohaemoglobins of microorganisms, especially of bacteria, and update several earlier and more comprehensive reviews, emphasising advances over the past 5 to 10 years and some controversies that have arisen. Inevitably, in light of space restrictions, details of nitric oxide metabolism and globins in higher organisms are brief.

Impact of azole drugs on energetics, kinetics, and ligand migration pathways of CO photo-dissociation in bacterial flavohemoglobins

Flavohemoglobins (fHbs) are heme proteins found in prokaryotic and eukaryotic microbes. They are involved in NO detoxification through an NO˙ dioxygenase mechanism. The N-terminal heme globin domain allows for binding of gaseous ligands whereas a C-terminal NADH/FADH binding domain facilitates association of redox cofactors necessary for ligand reduction. The NO˙ dioxygenase function is important in facilitating immune resistance by protecting the cell from nitrosative stress brought about by a host organism; as a result, bacterial flavoHbs have recently been considered as targets for the development of new antibiotics. Here, photoacoustic calorimetry and transient absorption spectroscopy have been used to characterize energetics, structural dynamics, and kinetics of CO migration within bacterial flavoHbs from Ralstonia eutropha (FHP) and Staphylococcus aureus (HMP_Sa) in the presence and absence of antibiotic azole compounds. In FHP, the ligand photo-release is associated with ΔH = 26.2 ± 7.0 kcal mol⁻¹ and ΔV = 25.0 ± 1.5 mL mol⁻¹ while in HMP_Sa, ΔH = 34.7 ± 8.0 kcal mol⁻¹ and ΔV = 28.6 ± 17 mL mol⁻¹ were observed, suggesting distinct structural changes associated with ligand escape from FHP and HMP_Sa. In the presence of ketoconazole, the CO escape leads to a more negative enthalpy change and volume change whereas association of miconazole to FHP or HMP_Sa does not impact the reaction volume. These data are in agreement with the computational results that propose distinct binding sites for ketoconazole and miconazole on CO bound FHP. Miconazole or ketoconazole binding to either protein has only a negligible impact on the CO association rates, indicating that azole drugs do not impact flavoHbs interactions with gaseous ligands but may inhibit the NOD activity through preventing the electron transfer between FAD and heme cofactors.

Impact of ketoconalzole and miconazole on structural dynamics of flavohemoglobin.

Quinones and nitroaromatic compounds as subversive substrates of Staphylococcus aureus flavohemoglobin

In microorganisms, flavohemoglobins (FHbs) containing FAD and heme (Fe³⁺, metHb) convert NO. into nitrate at the expense of NADH and O₂. FHbs contribute to bacterial resistance to nitrosative stress. Therefore, inhibition of FHbs functions may decrease the pathogen virulence. We report here a kinetic study of the reduction of quinones and nitroaromatic compounds by S. aureus FHb. We show that this enzyme rapidly reduces quinones and nitroaromatic compounds in a mixed single- and two-electron pathway. The reactivity of nitroaromatics increased upon an increase in their single-electron reduction potential (E¹₇), whereas the reactivity of quinones poorly depended on their E¹₇ with a strong preference for a 2-hydroxy-1,4-naphthoquinone structure. The reaction followed a 'ping-pong' mechanism. In general, the maximal reaction rates were found lower than the maximal presteady-state rate of FAD reduction by NADH and/or of oxyhemoglobin (HbFe²⁺O₂) formation (~130 s^-1, pH 7.0, 25 °C), indicating that the enzyme turnover is limited by the oxidative half-reaction. The turnover studies showed that quinones prefreqently accept electrons from reduced FAD, and not from HbFe²⁺O₂. These results suggest that quinones and nitroaromatics act as 'subversive substrates' for FHb, and may enhance the cytotoxicity of NO. by formation of superoxide and by diverting the electron flux coming from reduced FAD. Because quinone reduction rate was increased by FHb inhibitors such as econazole, ketoconazole, and miconazole, their combined use may represent a novel chemotherapeutical approach.

Charge-Transfer Knowledge Graph among Amino Acids Derived from High-Throughput Electronic Structure Calculations for Protein Database

The charge-transfer coupling is an important component in tight-binding methods. Because of the highly complex chemical structure of biomolecules, the anisotropic feature of charge-transfer couplings in realistic proteins cannot be ignored. In this work, we have performed the first large-scale quantitative assessment of charge-transfer preference by calculating the charge-transfer couplings in all 20 × 20 possible amino acid side-chain combinations, which are extracted from available high-quality structures of thousands of protein complexes. The charge-transfer database quantitatively shows distinct features of charge-transfer couplings among millions of amino acid side-chain combinations. The overall distribution of charge-transfer couplings reveals that only one average or representative structure cannot be regarded as the typical charge-transfer preference in realistic proteins. This work provides us an alternative route to comprehensively understand the charge-transfer couplings for the overall distribution of realistic proteins in the foreseen big data scenario.

AMOEBA Polarizable Force Field Parameters of the Heme Cofactor in Its Ferrous and Ferric Forms

We report the first parameters of the heme redox cofactors for the polarizable AMOEBA force field in both the ferric and ferrous forms. We consider two types of complexes, one with two histidine side chains as axial ligands and one with a histidine and a methionine side chain as ligands. We have derived permanent multipoles from second-order Møller-Plesset perturbation theory (MP2). The sets of parameters have been validated in a first step by comparison of AMOEBA interaction energies of heme and a collection of biologically relevant molecules with MP2 and Density Functional Theory (DFT) calculations. In a second validation step, we consider interaction energies with large aggregates comprising around 80 H₂O molecules. These calculations are repeated for 30 structures extracted from semiempirical PM7 DM simulations. Very encouraging agreement is found between DFT and the AMOEBA force field, which results from an accurate treatment of electrostatic interactions. We finally report long (10 ns) MD simulations of cytochromes in two redox states with AMOEBA testing both the 2003 and 2014 AMOEBA water models. These simulations have been carried out with the TINKER-HP (High Performance) program. In conclusion, owing to their ubiquity in biology, we think the present work opens a wide array of applications of the polarizable AMOEBA force field on hemeproteins.

Quantum chemical calculations of tryptophan → heme electron and excitation energy transfer rates in myoglobin

The development of optical multidimensional spectroscopic techniques has opened up new possibilities for the study of biological processes. Recently, ultrafast two-dimensional ultraviolet spectroscopy experiments have determined the rates of tryptophan → heme electron transfer and excitation energy transfer for the two tryptophan residues in myoglobin (Consani et al., Science, 2013, 339, 1586). Here, we show that accurate prediction of these rates can be achieved using Marcus theory in conjunction with time-dependent density functional theory. Key intermediate residues between the donor and acceptor are identified, and in particular the residues Val68 and Ile75 play a critical role in calculations of the electron coupling matrix elements. Our calculations demonstrate how small changes in structure can have a large effect on the rates, and show that the different rates of electron transfer are dictated by the distance between the heme and tryptophan residues, while for excitation energy transfer the orientation of the tryptophan residues relative to the heme is important. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Progress and challenges in simulating and understanding electron transfer in proteins

This Review presents an overview of the most common numerical simulation approaches for the investigation of electron transfer (ET) in proteins. We try to highlight the merits of the different approaches but also the current limitations and challenges. The article is organized into three sections. Section 2 deals with direct simulation algorithms of charge migration in proteins. Section 3 summarizes the methods for testing the applicability of the Marcus theory for ET in proteins and for evaluating key thermodynamic quantities entering the reaction rates (reorganization energies and driving force). Recent studies interrogating the validity of the theory due to the presence of non-ergodic effects or of non-linear responses are also described. Section 4 focuses on the tunneling aspects of electron transfer. How can the electronic coupling between charge transfer states be evaluated by quantum chemistry approaches and rationalized? What interesting physics regarding the impact of protein dynamics on tunneling can be addressed? We will illustrate the different sections with examples taken from the literature to show what types of system are currently manageable with current methodologies. We also take care to recall what has been learned on the biophysics of ET within proteins thanks to the advent of atomistic simulations.

Electron transfer, decoherence, and protein dynamics: insights from atomistic simulations

Electron transfer in biological systems drives the processes of life. From cellular respiration to photosynthesis and enzymatic catalysis, electron transfers (ET) are chemical processes on which essential biological functions rely. Over the last 40 years, scientists have sought understanding of how these essential processes function in biology. One important breakthrough was the discovery that Marcus theory (MT) of electron transfer is applicable to biological systems. Chemists have experimentally collected both the reorganization energies (λ) and the driving forces (ΔG°), two parameters of Marcus theory, for a large variety of ET processes in proteins. At the same time, theoretical chemists have developed computational approaches that rely on molecular dynamics and quantum chemistry calculations to access numerical estimates of λ and ΔG°. Yet another crucial piece in determining the rate of an electron transfer is the electronic coupling between the initial and final electronic wave functions. This is an important prefactor in the nonadiabatic rate expression, since it reflects the probability that an electron tunnels from the electron donor to the acceptor through the intervening medium. The fact that a protein matrix supports electron tunneling much more efficiently than vacuum is now well documented, both experimentally and theoretically. Meanwhile, many chemists have provided examples of the rich physical chemistry that can be induced by protein dynamics. This Account describes our studies of the dynamical effects on electron tunneling. We present our analysis of two examples of natural biological systems through MD simulations and tunneling pathway analyses. Through these examples, we show that protein dynamics sustain efficient tunneling. Second, we introduce two time scales: τcoh and τFC. The former characterizes how fast the electronic coupling varies with nuclear vibrations (which cause dephasing). The latter reflects the time taken by the system to leave the crossing region. In the framework of open quantum systems, τFC is a short time approximation of the characteristic decoherence time of the electronic subsystem in interaction with its nuclear environment. The comparison of the respective values of τcoh and τFC allows us to probe the occurrence of non-Condon effects. We use ab initio MD simulations to analyze how decoherence appears in several biological cofactors. We conclude that we cannot account for its order of magnitude by considering only the atoms or bonds directly concerned with the transfer. Decoherence results from contributions from all atoms of the system appearing with a time delay that increases with the distance from the primarily concerned atoms or bonds. The delay and magnitude of the contributions depend on the chemical nature of the system. Finally, we present recent developments based on constrained DFT for efficient and accurate evaluations of the electronic coupling in ab initio MD simulations. These are promising methods to study the subtle fluctuations of the electronic coupling and the mechanisms of electronic decoherence in biological systems.

ATP binding and aspartate protonation enhance photoinduced electron transfer in plant cryptochrome

Cryptochromes are flavoproteins encountered in most vegetal and animal species. They play a role of blue-light receptors in plants and in invertebrates. The putative resting state of the FAD cofactor in these proteins is its fully oxidized form, FADox. Upon blue-light excitation, the isoalloxazine ring (ISO) may undergo an ultrafast reduction by a nearby tryptophan residue W400. This primary reduction triggers a cascade of electron and proton transfers, ultimately leading to the formation of the FADH° radical. A recent experimental study has shown that the yield of FADH° formation in Arabidopsis cryptochrome can be strongly modulated by ATP binding and by pH, affecting the protonation state of D396 (proton donor to FAD°(-)). Here we provide a detailed molecular analysis of these effects by means of combined classical molecular dynamics simulations and time-dependent density functional theory calculations. When ATP is present and D396 protonated, FAD remains in close contact with W400, thereby enhancing electron transfer (ET) from W400 to ISO*. In contrast, deprotonation of D396 and absence of ATP introduce flexibility to the photoactive site prior to FAD excitation, with the consequence of increased ISO-W400 distance and diminished tunneling rate by almost two orders of magnitude. We show that under these conditions, ET from the adenine moiety of FAD becomes a competitive relaxation pathway. Overall, our data suggest that the observed effects of ATP and pH on the FAD photoreduction find their roots in the earliest stage of the photoreduction process; i.e., ATP binding and the protonation state of D396 determine the preferred pathway of ISO* relaxation.

Protein electron transfer: Dynamics and statistics

Electron transfer between redox proteins participating in energy chains of biology is required to proceed with high energetic efficiency, minimizing losses of redox energy to heat. Within the standard models of electron transfer, this requirement, combined with the need for unidirectional (preferably activationless) transitions, is translated into the need to minimize the reorganization energy of electron transfer. This design program is, however, unrealistic for proteins whose active sites are typically positioned close to the polar and flexible protein-water interface to allow inter-protein electron tunneling. The high flexibility of the interfacial region makes both the hydration water and the surface protein layer act as highly polar solvents. The reorganization energy, as measured by fluctuations, is not minimized, but rather maximized in this region. Natural systems in fact utilize the broad breadth of interfacial electrostatic fluctuations, but in the ways not anticipated by the standard models based on equilibrium thermodynamics. The combination of the broad spectrum of static fluctuations with their dispersive dynamics offers the mechanism of dynamical freezing (ergodicity breaking) of subsets of nuclear modes on the time of reaction/residence of the electron at a redox cofactor. The separation of time-scales of nuclear modes coupled to electron transfer allows dynamical freezing. In particular, the separation between the relaxation time of electro-elastic fluctuations of the interface and the time of conformational transitions of the protein caused by changing redox state results in dynamical freezing of the latter for sufficiently fast electron transfer. The observable consequence of this dynamical freezing is significantly different reorganization energies describing the curvature at the bottom of electron-transfer free energy surfaces (large) and the distance between their minima (Stokes shift, small). The ratio of the two reorganization energies establishes the parameter by which the energetic efficiency of protein electron transfer is increased relative to the standard expectations, thus minimizing losses of energy to heat. Energetically efficient electron transfer occurs in a chain of conformationally quenched cofactors and is characterized by flattened free energy surfaces, reminiscent of the flat and rugged landscape at the stability basin of a folded protein.

Microbial eukaryote globins

A bioinformatics survey of about 120 protist and 240 fungal genomes and transcriptomes revealed a broad array of globins, representing five of the eight subfamilies identified in bacteria. Most conspicuous is the absence of protoglobins and globin-coupled sensors, except for a two-domain globin in Leishmanias, that comprises a nucleotidyl cyclase domain, and the virtual absence of truncated group 3 globins. In contrast to bacteria, co-occurrence of more than two globin subfamilies appears to be rare in protists. Although globins were lacking in the Apicomplexa and the Microsporidia intracellular pathogens, they occurred in the pathogenic Trypanosomatidae, Stramenopiles and certain fungi. Flavohaemoglobins (FHbs) and related single-domain globins occur across the protist groups. Fungi are unique in having FHbs co-occurring with sensor single-domain globins (SSDgbs). Obligately biotrophic fungi covered in our analysis lack globins. Furthermore, SSDgbs occur only in a heterolobosean amoeba, Naegleria and the stramenopile Hyphochytrium. Of the three subfamilies of truncated Mb-fold globins, TrHb1s appear to be the most widespread, occurring as multiple copies in chlorophyte and ciliophora genomes, many as multidomain proteins. Although the ciliates appear to have only TrHb1s, the chlorophytes have Mb-like globins and TrHb2s, both closely related to the corresponding plant globins. The presently available number of protist genomes is inadequate to provide a definitive census of their globins. Bayesian molecular analyses of single-domain 3/3 Mb-fold globins suggest a close relationship of chlorophyte and haptophyte globins, including choanoflagellate and Capsaspora globins to land plant symbiotic and non-symbiotic haemoglobins and to vertebrate neuroglobins.

Hemoglobin: a nitric-oxide dioxygenase

Members of the hemoglobin superfamily efficiently catalyze nitric-oxide dioxygenation, and when paired with native electron donors, function as NO dioxygenases (NODs). Indeed, the NOD function has emerged as a more common and ancient function than the well-known role in O2 transport-storage. Novel hemoglobins possessing a NOD function continue to be discovered in diverse life forms. Unique hemoglobin structures evolved, in part, for catalysis with different electron donors. The mechanism of NOD catalysis by representative single domain hemoglobins and multidomain flavohemoglobin occurs through a multistep mechanism involving O2 migration to the heme pocket, O2 binding-reduction, NO migration, radical-radical coupling, O-atom rearrangement, nitrate release, and heme iron re-reduction. Unraveling the physiological functions of multiple NODs with varying expression in organisms and the complexity of NO as both a poison and signaling molecule remain grand challenges for the NO field. NOD knockout organisms and cells expressing recombinant NODs are helping to advance our understanding of NO actions in microbial infection, plant senescence, cancer, mitochondrial function, iron metabolism, and tissue O2 homeostasis. NOD inhibitors are being pursued for therapeutic applications as antibiotics and antitumor agents. Transgenic NOD-expressing plants, fish, algae, and microbes are being developed for agriculture, aquaculture, and industry.