Electron tunneling in respiratory complex I

Searching for proton transfer channels in respiratory complex I

We have explored a strategy to identify potential proton transfer channels using computational analysis of a protein structure based on Voronoi partitioning and applied it for the analysis of proton transfer pathways in redox-driven proton-pumping respiratory complex I. The analysis results in a network of connected voids/channels, which represent the dual structure of the protein; we then hydrated the identified channels using our water placement program Dowser++. Many theoretical water molecules found in the channels perfectly match the observed experimental water molecules in the structure; some other predicted water molecules have not been resolved in the experiments. The channels are of varying cross sections. Some channels are big enough to accommodate water molecules that are suitable to conduct protons; others are too narrow to hold water but require only minor conformational changes to accommodate proton transfer. We provide a preliminary analysis of the proton conductivity of the network channels, classifying the proton transfer channels as open, closed, and partially open, and discuss possible conformational changes that can modulate, i.e., open and close, the channels.

Properties of the iron-sulfur cluster electron transfer relay in an [FeFe]-hydrogenase that is tuned for H2 oxidation catalysis

[FeFe]-hydrogenases catalyze the reversible oxidation of H2 from electrons and protons at an organometallic active site cofactor named the H-cluster. In addition to the H-cluster, most [FeFe]-hydrogenases possess accessory FeS cluster (F-cluster) relays that function in mediating electron transfer with catalysis. There is significant variation in the structural properties of F-cluster relays among the [FeFe]-hydrogenases; however, it is unknown how this variation relates to the electronic and thermodynamic properties, and thus the electron transfer properties, of enzymes. Clostridium pasteurianum [FeFe]-hydrogenase II (CpII) exhibits a large catalytic bias for H2 oxidation (compared to H2 production), making it a notable system for examining if F-cluster properties contribute to the overall function and efficiency of the enzyme. By applying a combination of multifrequency and potentiometric electron paramagnetic resonance, we resolved two electron paramagnetic resonance signals with distinct power- and temperature-dependent properties at g = 2.058 1.931 1.891 (F2.058) and g = 2.061 1.920 1.887 (F2.061), with assigned midpoint potentials of -140 ± 18 mV and -406 ± 12 mV versus normal hydrogen electrode, respectively. Spectral analysis revealed features consistent with spin-spin coupling between the two [4Fe-4S] F-clusters, and possible functional models are discussed that account for the contribution of coupling to the electron transfer landscape. The results signify the interplay of electronic coupling and free energy properties and parameters of the FeS clusters to the electron transfer mechanism through the relay and provide new insight as to how relays functionally complement the catalytic directionality of active sites to achieve highly efficient catalysis.

How Mitochondrial Signaling Games May Shape and Stabilize the Nuclear-Mitochondrial Symbiosis

The eukaryotic lineage has enjoyed a long-term "stable" mutualism between nucleus and mitochondrion, since mitochondrial endosymbiosis began about 2 billion years ago. This mostly cooperative interaction has provided the basis for eukaryotic expansion and diversification, which has profoundly altered the forms of life on Earth. While we ignore the exact biochemical details of how the alpha-proteobacterial ancestor of mitochondria entered into endosymbiosis with a proto-eukaryote, in more general terms, we present a signaling games perspective of how the cooperative relationship became established, and has been maintained. While games are used to understand organismal evolution, information-asymmetric games at the molecular level promise novel insights into endosymbiosis. Using a previously devised biomolecular signaling games approach, we model a sender-receiver information asymmetric game, in which the informed mitochondrial sender signals and the uninformed nuclear receiver may take actions (involving for example apoptosis, senescence, regeneration and autophagy/mitophagy). The simulation shows that cellularization is a stabilizing mechanism for Pareto efficient sender/receiver strategic interaction. In stark contrast, the extracellular environment struggles to maintain efficient outcomes, as senders are indifferent to the effects of their signals upon the receiver. Our hypothesis has translational implications, such as in cellular therapy, as mitochondrial medicine matures. It also inspires speculative conjectures about how an analogous human-AI endosymbiosis may be engineered.

H2S Protects from Rotenone-Induced Ferroptosis by Stabilizing Fe-S Clusters in Rat Cardiac Cells

Hydrogen sulfide (H2S) has been recently recognized as an important gasotransmitter with cardioprotections, and iron is vital for various cellular activities. This study explored the regulatory role of H2S on iron metabolism and mitochondrial functions in cultured rat cardiac cells. Rotenone, a mitochondrial complex I inhibitor, was used for establishing an in vitro model of ischemic cell damage. It was first found that rotenone induced oxidative stress and lipid peroxidation and decreased mitochondrial membrane potential and ATP generation, eventually causing cell death. The supplement of H2S at a physiologically relevant concentration protected from rotenone-induced ferroptotic cell death by reducing oxidative stress and mitochondrial damage, maintaining GPx4 expression and intracellular iron level. Deferiprone, an iron chelator, would also protect from rotenone-induced ferroptosis. Further studies demonstrated that H2S inhibited ABCB8-mediated iron efflux from mitochondria to cytosol and promoted NFS1-mediated Fe-S cluster biogenesis. It is also found that rotenone stimulated iron-dependent H2S generation. These results indicate that H2S would protect cardiac cells from ischemic damage through preserving mitochondrial functions and intracellular Fe-S cluster homeostasis.

Respiratory Complex III: A Bioengine with a Ligand-Triggered Electron-Tunneling Gating Mechanism

Respiratory complex III (a.k.a., the bc1 complex) plays a key role in the electron transport chain in aerobic cells. The bc1 complex exhibits multiple unique electron tunneling (ET) processes, such as ET-bifurcation at the Qo site and movement of the Rieske domain. Moreover, we previously discovered that electron tunneling in the low potential arm of the bc1 complex is regulated by a key phenylalanine residue (Phe90). The main goal of the current work is to study the dynamics of the key Phe90 residue in the electron tunneling reaction between heme bL and heme bH as a function of the occupancy of the Qo and Qi binding sites in the bc1 complex. We simulated the molecular dynamics of four model systems of respiratory complex III with different ligands bound at the Qo and Qi binding sites. In addition, we calculated the electron tunneling rate constants between heme bL and heme bH along the simulated molecular dynamics trajectories. The binding of aromatic ligands at the Qo site induces a conformational cascade that properly positions the Phe90 residue, reducing the through-space ET distance from ∼7 to ∼5.5 Å and thus enhancing the electron transfer rate between the heme bL and the heme bH redox pair. Also, the binding of aromatic ligands at the Qi site induces conformational changes that stabilize the Phe90 conformational variation from ∼1.5 to ∼0.5 Å. Hence, our molecular dynamics simulation results show an on-demand two-step conformational connection between the occupancy of the Qo and Qi binding sites and the conformational dynamics of the Phe90 residue. Additionally, our dynamic electron tunneling results confirm our previously reported findings that the Phe90 residue acts as an electron-tunneling gate or switch, controlling the electron transfer rate between the heme bL and heme bH redox systems.

Mitochondrial complex I ROS production and redox signaling in hypoxia

Mitochondria are a main source of cellular energy. Oxidative phosphorylation (OXPHOS) is the major process of aerobic respiration. Enzyme complexes of the electron transport chain (ETC) pump protons to generate a protonmotive force (Δp) that drives OXPHOS. Complex I is an electron entry point into the ETC. Complex I oxidizes nicotinamide adenine dinucleotide (NADH) and transfers electrons to ubiquinone in a reaction coupled with proton pumping. Complex I also produces reactive oxygen species (ROS) under various conditions. The enzymatic activities of complex I can be regulated by metabolic conditions and serves as a regulatory node of the ETC. Complex I ROS plays diverse roles in cell metabolism ranging from physiologic to pathologic conditions. Progress in our understanding indicates that ROS release from complex I serves important signaling functions. Increasing evidence suggests that complex I ROS is important in signaling a mismatch in energy production and demand. In this article, we review the role of ROS from complex I in sensing acute hypoxia.

Impact of weak radiofrequency and static magnetic fields on key signaling molecules, intracellular pH, membrane potential, and cell growth in HT-1080 fibrosarcoma cells

There are substantial concerns that extended exposures to weak radiofrequency (RF) fields can lead to adverse health effects. In this study, HT-1080 fibrosarcoma cells were simultaneously exposed to a static magnetic flux density between 10 [Formula: see text] and 300 [Formula: see text] and RF magnetic fields with amplitudes ranging from 1 nT to 1.5 μT in the frequency range from 1.8 to 7.2 MHz for four days. Cell growth rates, intracellular pH, hydrogen peroxide, peroxynitrite, membrane potential and mitochondrial calcium were measured. Results were dependent on carrier frequency and the magnitude of the RF magnetic field, modulation frequencies and the background static magnetic field (SMF). Iron sulphur (Fe-S) clusters are essential for the generation of reactive oxygen species and reactive nitrogen species (ROS and RNS). We believe the observed changes are associated with hyperfine couplings between the chemically active electrons and nuclear spins. Controlling external magnetic fields may have important clinical implications on aging, cancer, arthritis, and Alzheimer's.

Exploiting Molecular Symmetry to Quantitatively Map the Excited-State Landscape of Iron-Sulfur Clusters

Cuboidal [Fe4S4] clusters are ubiquitous cofactors in biological redox chemistry. In the [Fe4S4]1+ state, pairwise spin coupling gives rise to six arrangements of the Fe valences ("valence isomers") among the four Fe centers. Because of the magnetic complexity of these systems, it has been challenging to understand how a protein's active site dictates both the arrangement of the valences in the ground state as well as the population of excited-state valence isomers. Here, we show that the ground-state valence isomer landscape can be simplified from a six-level system in an asymmetric protein environment to a two-level system by studying the problem in synthetic [Fe4S4]1+ clusters with solution C3v symmetry. This simplification allows for the energy differences between valence isomers to be quantified (in some cases with a resolution of <0.1 kcal/mol) by simultaneously fitting the VT NMR and solution magnetic moment data. Using this fitting protocol, we map the excited-state landscape for a range of clusters of the form [(SIMes)3Fe4S4-X/L]n, (SIMes = 1,3-dimesityl-imidazol-4,5-dihydro-2-ylidene; n = 0 for anionic, X-type ligands and n = +1 for neutral, L-type ligands) and find that a single ligand substitution can alter the relative ground-state energies of valence isomers by at least 103 cm-1. On this basis, we suggest that one result of "non-canonical" amino acid ligation in Fe-S proteins is the redistribution of the valence electrons in the manifold of thermally populated excited states.

Single protonation of the reduced quinone in respiratory complex I drives four-proton pumping

Complex I is a key proton-pumping enzyme in bacterial and mitochondrial respiratory electron transport chains. Using quantum chemistry and electrostatic calculations, we have examined the pKa of the reduced quinone QH-/QH2 in the catalytic cavity of complex I. We find that pKa (QH-/QH2) is very high, above 20. This means that the energy of a single protonation reaction of the doubly reduced quinone (i.e. the reduced semiquinone QH-) is sufficient to drive four protons across the membrane with a potential of 180 mV. Based on these calculations, we propose a possible scheme of redox-linked proton pumping by complex I. The model explains how the energy of the protonation reaction can be divided equally among four pumping units of the pump, and how a single proton can drive translocation of four additional protons in multiple pumping blocks.

Chronobiology Meets Quantum Biology: A New Paradigm Overlooking the Horizon?

Biological processes and physiological functions in living beings are featured by oscillations with a period of about 24 h (circadian) or cycle at the second and third harmonic (ultradian) of the basic frequency, driven by the biological clock. This molecular mechanism, common to all kingdoms of life, comprising animals, plants, fungi, bacteria, and protists, represents an undoubted adaptive advantage allowing anticipation of predictable changes in the environmental niche or of the interior milieu. Biological rhythms are the field of study of Chronobiology. In the last decade, growing evidence hints that molecular platforms holding up non-trivial quantum phenomena, including entanglement, coherence, superposition and tunnelling, bona fide evolved in biosystems. Quantum effects have been mainly implicated in processes related to electromagnetic radiation in the spectrum of visible light and ultraviolet rays, such as photosynthesis, photoreception, magnetoreception, DNA mutation, and not light related such as mitochondrial respiration and enzymatic activity. Quantum effects in biological systems are the field of study of Quantum Biology. Rhythmic changes at the level of gene expression, as well as protein quantity and subcellular distribution, confer temporal features to the molecular platform hosting electrochemical processes and non-trivial quantum phenomena. Precisely, a huge amount of molecules plying scaffold to quantum effects show rhythmic level fluctuations and this biophysical model implies that timescales of biomolecular dynamics could impinge on quantum mechanics biofunctional role. The study of quantum phenomena in biological cycles proposes a profitable “entanglement” between the areas of interest of these seemingly distant scientific disciplines to enlighten functional roles for quantum effects in rhythmic biosystems.

A site-differentiated [4Fe-4S] cluster controls electron transfer reactivity of Clostridium acetobutylicum [FeFe]-hydrogenase I

One of the many functions of reduction–oxidation (redox) cofactors is to mediate electron transfer in biological enzymes catalyzing redox-based chemical transformation reactions. There are numerous examples of enzymes that utilize redox cofactors to form electron transfer relays to connect catalytic sites to external electron donors and acceptors. The compositions of relays are diverse and tune transfer thermodynamics and kinetics towards the chemical reactivity of the enzyme. Diversity in relay design is exemplified among different members of hydrogenases, enzymes which catalyze reversible H₂ activation, which also couple to diverse types of donor and acceptor molecules. The [FeFe]-hydrogenase I from Clostridium acetobutylicum (CaI) is a member of a large family of structurally related enzymes where interfacial electron transfer is mediated by a terminal, non-canonical, His-coordinated, [4Fe–4S] cluster. The function of His coordination was examined by comparing the biophysical properties and reactivity to a Cys substituted variant of CaI. This demonstrated that His coordination strongly affected the distal [4Fe–4S] cluster spin state, spin pairing, and spatial orientations of molecular orbitals, with a minor effect on reduction potential. The deviations in these properties by substituting His for Cys in CaI, correlated with pronounced changes in electron transfer and reactivity with the native electron donor–acceptor ferredoxin. The results demonstrate that differential coordination of the surface localized [4Fe–4S]His cluster in CaI is utilized to control intermolecular and intramolecular electron transfer where His coordination creates a physical and electronic environment that enables facile electron exchange between electron carrier molecules and the iron–sulfur cluster relay for coupling to reversible H₂ activation at the catalytic site.

Histidine coordination of the distal [4Fe–4S] cluster in [FeFe]-hydrogenase was demonstrated to tune the cluster spin-states, spin-pairing and surrounding molecular orbitals to enable more facile electron transfer compared to cysteine coordination.

Resolving Chemical Dynamics in Biological Energy Conversion: Long-Range Proton-Coupled Electron Transfer in Respiratory Complex I

Biological energy conversion is catalyzed by membrane-bound proteins that transduce chemical or light energy into energy forms that power endergonic processes in the cell. At a molecular level, these catalytic processes involve elementary electron-, proton-, charge-, and energy-transfer reactions that take place in the intricate molecular machineries of cell respiration and photosynthesis. Recent developments in structural biology, particularly cryo-electron microscopy (cryoEM), have resolved the molecular architecture of several energy transducing proteins, but detailed mechanistic principles of their charge transfer reactions still remain poorly understood and a major challenge for modern biochemical research. To this end, multiscale molecular simulations provide a powerful approach to probe mechanistic principles on a broad range of time scales (femtoseconds to milliseconds) and spatial resolutions (10¹–10⁶ atoms), although technical challenges also require balancing between the computational accuracy, cost, and approximations introduced within the model. Here we discuss how the combination of atomistic (aMD) and hybrid quantum/classical molecular dynamics (QM/MM MD) simulations with free energy (FE) sampling methods can be used to probe mechanistic principles of enzymes responsible for biological energy conversion. We present mechanistic explorations of long-range proton-coupled electron transfer (PCET) dynamics in the highly intricate respiratory chain enzyme Complex I, which functions as a redox-driven proton pump in bacterial and mitochondrial respiratory chains by catalyzing a 300 Å fully reversible PCET process. This process is initiated by a hydride (H^–) transfer between NADH and FMN, followed by long-range (>100 Å) electron transfer along a wire of 8 FeS centers leading to a quinone biding site. The reduction of the quinone to quinol initiates dissociation of the latter to a second membrane-bound binding site, and triggers proton pumping across the membrane domain of complex I, in subunits up to 200 Å away from the active site. Our simulations across different size and time scales suggest that transient charge transfer reactions lead to changes in the internal hydration state of key regions, local electric fields, and the conformation of conserved ion pairs, which in turn modulate the dynamics of functional steps along the reaction cycle. Similar functional principles, which operate on much shorter length scales, are also found in some unrelated proteins, suggesting that enzymes may employ conserved principles in the catalysis of biological energy transduction processes.

Dehydrogenases of acetic acid bacteria

Acetic acid bacteria (AAB) are a group of bacteria that can oxidize many substrates such as alcohols and sugar alcohols and play important roles in industrial biotechnology. A majority of industrial processes that involve AAB are related to their dehydrogenases, including PQQ/FAD-dependent membrane-bound dehydrogenases and NAD(P)⁺-dependent cytoplasmic dehydrogenases. These cofactor-dependent dehydrogenases must effectively regenerate their cofactors in order to function continuously. For PQQ, FAD and NAD(P)⁺ alike, regeneration is directly or indirectly related to the electron transport chain (ETC) of AAB, which plays an important role in energy generation for aerobic cell growth. Furthermore, in changeable natural habitats, ETC components of AAB can be regulated so that the bacteria survive in different environments. Herein, the progressive cascade in an application of AAB, including key dehydrogenases involved in the application, regeneration of dehydrogenase cofactors, ETC coupling with cofactor regeneration and ETC regulation, is systematically reviewed and discussed. As they have great application value, a deep understanding of the mechanisms through which AAB function will not only promote their utilization and development but also provide a reference for engineering of other industrial strains.

Photoinduced hole hopping through tryptophans in proteins

Hole hopping through tryptophan/tyrosine chains enables rapid unidirectional charge transport over long distances. We have elucidated structural and dynamical factors controlling hopping speed and efficiency in two modified azurin constructs that include a rhenium(I) sensitizer, Re(His)(CO)₃(dmp)⁺, and one or two tryptophans (W₁, W₂). Experimental kinetics investigations showed that the two closely spaced (3 to 4 Å) intervening tryptophans dramatically accelerated long-range electron transfer (ET) from Cu^I to the photoexcited sensitizer. In our theoretical work, we found that time-dependent density-functional theory (TDDFT) quantum mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) trajectories of low-lying triplet excited states of Re^I(His)(CO)₃(dmp)⁺-W₁(-W₂) exhibited crossings between sensitizer-localized (*Re) and charge-separated [Re^I(His)(CO)₃(dmp^•-)/(W₁ ^•+ or W₂ ^•+)] (CS1 or CS2) states. Our analysis revealed that the distances, angles, and mutual orientations of ET-active cofactors fluctuate in a relatively narrow range in which the cofactors are strongly coupled, enabling adiabatic ET. Water-dominated electrostatic field fluctuations bring *Re and CS1 states to a crossing where *Re(CO)₃(dmp)⁺←W₁ ET occurs, and CS1 becomes the lowest triplet state. ET is promoted by solvation dynamics around *Re(CO)₃(dmp)⁺(W₁); and CS1 is stabilized by Re(dmp^•-)/W₁ ^•+ electron/hole interaction and enhanced W₁ ^•+ solvation. The second hop, W₁ ^•+←W₂, is facilitated by water fluctuations near the W₁/W₂ unit, taking place when the electrostatic potential at W₂ drops well below that at W₁ ^•+ Insufficient solvation and reorganization around W₂ make W₁ ^•+←W₂ ET endergonic, shifting the equilibrium toward W₁ ^•+ and decreasing the charge-separation yield. We suggest that multiscale TDDFT/MM/MD is a suitable technique to model the simultaneous evolution of photogenerated excited-state manifolds.

Electron transfer pathways from quantum dynamics simulations

This work explores the possibility of simulating an electron transfer process between a donor and an acceptor in real time using time-dependent density functional theory electron dynamics. To achieve this objective, a central issue to resolve is the definition of the initial state. This must be a non-equilibrium electronic state able to trigger the charge transfer dynamics; here, two schemes are proposed to prepare such states. One is based on the combination of the density matrices of the donor and acceptor converged separately with appropriate charges (for example, -1 for the donor and +1 for the acceptor). The second approach relied on constrained DFT to localize the charge on each fragment. With these schemes, electron transfer processes are simulated in different model systems of increasing complexity: an atomic hydrogen dimer, a polyacetylene chain, and the active site of the T. cruzi hybrid type A heme peroxidase, for which two possible electron transfer paths have been postulated. For the latter system, the present methodology applied in a hybrid Quantum Mechanics - Molecular Mechanics framework allows us to establish the relative probabilities of each path and provides insight into the inhibition of the electron transfer provoked by the substitution of tryptophan by phenylalanine in the W233F mutant.

How the Donor/Acceptor Spin States Affect the Electronic Couplings in Molecular Charge-Transfer Processes?

The electronic coupling matrix element H_AB is an essential ingredient of most electron-transfer theories. H_AB depends on the overlap between donor and acceptor wave functions and is affected by the involved states’ spin. We classify the spin-state effects into three categories: orbital occupation, spin-dependent electron density, and density delocalization. The orbital occupancy reflects the diverse chemical nature and reactivity of the spin states of interest. The effect of spin-dependent density is related to a more compact electron density cloud at lower spin states due to decreased exchange interactions between electrons. Density delocalization is strongly connected with the covalency concept that increases the spatial extent of the diabatic state’s electron density in specific directions. We illustrate these effects with high-level ab initio calculations on model direct donor–acceptor systems relevant to metal oxide materials and biological electron transfer. Obtained results can be used to benchmark existing methods for H_AB calculations in complicated cases such as spin-crossover materials or antiferromagnetically coupled systems.

Functional Water Wires Catalyze Long-Range Proton Pumping in the Mammalian Respiratory Complex I

The respiratory complex I is a gigantic (1 MDa) redox-driven proton pump that reduces the ubiquinone pool and generates proton motive force to power ATP synthesis in mitochondria. Despite resolved molecular structures and biochemical characterization of the enzyme from multiple organisms, its long-range (∼300 Å) proton-coupled electron transfer (PCET) mechanism remains unsolved. We employ here microsecond molecular dynamics simulations to probe the dynamics of the mammalian complex I in combination with hybrid quantum/classical (QM/MM) free energy calculations to explore how proton pumping reactions are triggered within its 200 Å wide membrane domain. Our simulations predict extensive hydration dynamics of the antiporter-like subunits in complex I that enable lateral proton transfer reactions on a microsecond time scale. We further show how the coupling between conserved ion pairs and charged residues modulate the proton transfer dynamics, and how transmembrane helices and gating residues control the hydration process. Our findings suggest that the mammalian complex I pumps protons by tightly linked conformational and electrostatic coupling principles.

Mimicking Natural Photosynthesis: Designing Ultrafast Photosensitized Electron Transfer into Multiheme Cytochrome Protein Nanowires

Efficient nanomaterials for artificial photosynthesis require fast and robust unidirectional electron transfer (ET) from photosensitizers through charge-separation and accumulation units to redox-active catalytic sites. We explored the ultrafast time-scale limits of photo-induced charge transfer between a Ru(II)tris(bipyridine) derivative photosensitizer and PpcA, a 3-heme c-type cytochrome serving as a nanoscale biological wire. Four covalent attachment sites (K28C, K29C, K52C, and G53C) were engineered in PpcA enabling site-specific covalent labeling with expected donor-acceptor (DA) distances of 4-8 Å. X-ray scattering results demonstrated that mutations and chemical labeling did not disrupt the structure of the proteins. Time-resolved spectroscopy revealed three orders of magnitude difference in charge transfer rates for the systems with otherwise similar DA distances and the same number of covalent bonds separating donors and acceptors. All-atom molecular dynamics simulations provided additional insight into the structure-function requirements for ultrafast charge transfer and the requirement of van der Waals contact between aromatic atoms of photosensitizers and hemes in order to observe sub-nanosecond ET. This work demonstrates opportunities to utilize multi-heme c-cytochromes as frameworks for designing ultrafast light-driven ET into charge-accumulating biohybrid model systems, and ultimately for mimicking the photosynthetic paradigm of efficiently coupling ultrafast, light-driven electron transfer chemistry to multi-step catalysis within small, experimentally versatile photosynthetic biohybrid assemblies.

Mutations in NDUFS1 Cause Metabolic Reprogramming and Disruption of the Electron Transfer

Complex I (CI) is the first enzyme of the mitochondrial respiratory chain and couples the electron transfer with proton pumping. Mutations in genes encoding CI subunits can frequently cause inborn metabolic errors. We applied proteome and metabolome profiling of patient-derived cells harboring pathogenic mutations in two distinct CI genes to elucidate underlying pathomechanisms on the molecular level. Our results indicated that the electron transfer within CI was interrupted in both patients by different mechanisms. We showed that the biallelic mutations in NDUFS1 led to a decreased stability of the entire N-module of CI and disrupted the electron transfer between two iron–sulfur clusters. Strikingly interesting and in contrast to the proteome, metabolome profiling illustrated that the pattern of dysregulated metabolites was almost identical in both patients, such as the inhibitory feedback on the TCA cycle and altered glutathione levels, indicative for reactive oxygen species (ROS) stress. Our findings deciphered pathological mechanisms of CI deficiency to better understand inborn metabolic errors.

Long-range proton-coupled electron transfer in biological energy conversion: towards mechanistic understanding of respiratory complex I

Biological energy conversion is driven by efficient enzymes that capture, store and transfer protons and electrons across large distances. Recent advances in structural biology have provided atomic-scale blueprints of these types of remarkable molecular machinery, which together with biochemical, biophysical and computational experiments allow us to derive detailed energy transduction mechanisms for the first time. Here, I present one of the most intricate and least understood types of biological energy conversion machinery, the respiratory complex I, and how its redox-driven proton-pump catalyses charge transfer across approximately 300 Å distances. After discussing the functional elements of complex I, a putative mechanistic model for its action-at-a-distance effect is presented, and functional parallels are drawn to other redox- and light-driven ion pumps.

Medical hypothesis: Neurodegenerative diseases arise from oxidative damage to electron tunneling proteins in mitochondria

Mitochondria likely arose from serial endosymbiosis by early eukaryotic cells and control electron flow to molecular oxygen to facilitate energy transformation. Mitochondria translate between the quantum and macroscopic worlds and utilize quantum tunneling of electrons to reduce activation energy barriers to electron flow. Electron tunneling has been extensively characterized in Complex I of the electron transport chain. Age-related increases in oxidative damage to these electron tunneling systems may account for decreased energy storage found in aged and neurodegenerative disease tissues, such as those from sufferers of amyotrophic lateral sclerosis (ALS), Alzheimer's disease (AD) and Parkinson's disease (PD). This hypothesis is testable. If correct, this hypothesis supports pre-symptomatic, mitochondrially-directed oxygen free radical scavenging therapies.

A density functional theory study of the hydrogenation and reduction of the thio-spinel Fe3S4{111} surface

The mineral greigite, Fe3S4, shows promising electro-reduction activity, especially towards carbon dioxide conversion to small organic molecules. We have employed density functional theory calculations with correction for the long-range dispersion forces to investigate the behavior of hydrogen on the greigite{111} surface. We have studied the adsorption, diffusion, surface reduction and associative (i.e. Volmer-Tafel mechanism) and molecular desorption of hydrogen as a function of its coverage. We found that (i) the H ad-atoms adsorb on S sites far from metallic centres in the topmost surface layer; (ii) the reduction of greigite by hydrogen is energetically unfavorable at any surface coverage; and (iii) molecular hydrogen evolution has a transition state at ∼0.5 eV above the energy of the reactants on Fe3S4{111}, which is very similar to the barrier found experimentally on Pt{111}. We have also determined the electrode potential under room conditions at which the H2 evolution reaction becomes energetically barrierless.

Two Tryptophans Are Better Than One in Accelerating Electron Flow through a Protein

We have constructed and structurally characterized a Pseudomonas aeruginosa azurin mutant Re126WWCuI , where two adjacent tryptophan residues (W124 and W122, indole separation 3.6-4.1 Å) are inserted between the CuI center and a Re photosensitizer coordinated to the imidazole of H126 (ReI(H126)(CO)3(4,7-dimethyl-1,10-phenanthroline)+). CuI oxidation by the photoexcited Re label (*Re) 22.9 Å away proceeds with a ∼70 ns time constant, similar to that of a single-tryptophan mutant (∼40 ns) with a 19.4 Å Re-Cu distance. Time-resolved spectroscopy (luminescence, visible and IR absorption) revealed two rapid reversible electron transfer steps, W124 → *Re (400-475 ps, K 1 ≅ 3.5-4) and W122 → W124•+ (7-9 ns, K 2 ≅ 0.55-0.75), followed by a rate-determining (70-90 ns) CuI oxidation by W122•+ ca. 11 Å away. The photocycle is completed by 120 μs recombination. No photochemical CuI oxidation was observed in Re126FWCuI , whereas in Re126WFCuI , the photocycle is restricted to the ReH126W124 unit and CuI remains isolated. QM/MM/MD simulations of Re126WWCuI indicate that indole solvation changes through the hopping process and W124 → *Re electron transfer is accompanied by water fluctuations that tighten W124 solvation. Our finding that multistep tunneling (hopping) confers a ∼9000-fold advantage over single-step tunneling in the double-tryptophan protein supports the proposal that hole-hopping through tryptophan/tyrosine chains protects enzymes from oxidative damage.

New Perspectives on Iron Uptake in Eukaryotes

All eukaryotic organisms require iron to function. Malfunctions within iron homeostasis have a range of physiological consequences, and can lead to the development of pathological conditions that can result in an excess of non-transferrin bound iron (NTBI). Despite extensive understanding of iron homeostasis, the links between the “macroscopic” transport of iron across biological barriers (cellular membranes) and the chemistry of redox changes that drive these processes still needs elucidating. This review draws conclusions from the current literature, and describes some of the underlying biophysical and biochemical processes that occur in iron homeostasis. By first taking a broad view of iron uptake within the gut and subsequent delivery to tissues, in addition to describing the transferrin and non-transferrin mediated components of these processes, we provide a base of knowledge from which we further explore NTBI uptake. We provide concise up-to-date information of the transplasma electron transport systems (tPMETSs) involved within NTBI uptake, and highlight how these systems are not only involved within NTBI uptake for detoxification but also may play a role within the reduction of metabolic stress through regeneration of intracellular NAD(P)H/NAD(P)⁺ levels. Furthermore, we illuminate the thermodynamics that governs iron transport, namely the redox potential cascade and electrochemical behavior of key components of the electron transport systems that facilitate the movement of electrons across the plasma membrane to the extracellular compartment. We also take account of kinetic changes that occur to transport iron into the cell, namely membrane dipole change and their consequent effects within membrane structure that act to facilitate transport of ions.

Charge transfer through a fragment of the respiratory complex I and its regulation: an atomistic simulation approach

We simulate electron transfer within a fragment of the NADH:ubiquinone oxidoreductase (respiratory complex I) of the hyperthermophilic bacterium Aquifex aeolicus. We apply molecular dynamics simulations, thermodynamic integration, and a thermodynamic network least squares analysis to compute two key parameters of Marcus' theory of charge transfer, the thermodynamic driving force and the reorganization energy. Intramolecular contributions to the Gibbs free energy differences of electron and hydrogen transfer processes, ΔG, are accessed by calibrating against experimental redox titration data. This approach permits the computation of the interactions between the species NAD+, FMNH2, N1a-, and N3-, and the construction of a free energy surface for the flow of electrons within the fragment. We find NAD+ to be a strong candidate for the regulation of charge transfer.

Fragment-orbital tunneling currents and electronic couplings for analysis of molecular charge-transfer systems

In theoretical charge-transfer research, calculation of the electronic coupling element is crucial for examining the degree of the electronic donor-acceptor interaction. The tunneling current (TC), representing the magnitudes and directions of electron flow, provides a way of evaluating electronic couplings, along with the ability of visualizing how electrons flow in systems. Here, we applied the TC theory to π-conjugated organic dimer systems, in the form of our fragment-orbital tunneling current (FOTC) method, which uses the frontier molecular-orbitals of system fragments as diabatic states. For a comprehensive test of FOTC, we assessed how reasonable the computed electronic couplings and the corresponding TC densities are for the hole- and electron-transfer databases HAB11 and HAB7. FOTC gave 12.5% mean relative unsigned error with regard to the high-level ab initio reference. The shown performance is comparable with that of fragment-orbital density functional theory, which gave the same error by 20.6% or 13.9% depending on the formulation. In the test of a set of nucleobase π stacks, we showed that the original TC expression is also applicable to nondegenerate cases under the condition that the overlap between the charge distributions of diabatic states is small enough to offset the energy difference. Lastly, we carried out visual analysis on the FOTC densities of thiophene dimers with different intermolecular alignments. The result depicts an intimate topological connection between the system geometry and electron flow. Our work provides quantitative and qualitative grounds for FOTC, showing it to be a versatile tool in characterization of molecular charge-transfer systems.

SIRT3 deregulation is linked to mitochondrial dysfunction in Alzheimer's disease

Alzheimer's disease (AD) is the leading cause of dementia in the elderly. Despite decades of study, effective treatments for AD are lacking. Mitochondrial dysfunction has been closely linked to the pathogenesis of AD, but the relationship between mitochondrial pathology and neuronal damage is poorly understood. Sirtuins (SIRT, silent mating type information regulation 2 homolog in yeast) are NAD-dependent histone deacetylases involved in aging and longevity. The objective of this study was to investigate the relationship between SIRT3 and mitochondrial function and neuronal activity in AD. SIRT3 mRNA and protein levels were significantly decreased in AD cerebral cortex, and Ac-p53 K320 was significantly increased in AD mitochondria. SIRT3 prevented p53-induced mitochondrial dysfunction and neuronal damage in a deacetylase activity-dependent manner. Notably, mitochondrially targeted p53 (mito-p53) directly reduced mitochondria DNA-encoded ND2 and ND4 gene expression resulting in increased reactive oxygen species (ROS) and reduced mitochondrial oxygen consumption. ND2 and ND4 gene expressions were significantly decreased in patients with AD. p53-ChIP analysis verified the presence of p53-binding elements in the human mitochondrial genome and increased p53 occupancy of mitochondrial DNA in AD. SIRT3 overexpression restored the expression of ND2 and ND4 and improved mitochondrial oxygen consumption by repressing mito-p53 activity. Our results indicate that SIRT3 dysfunction leads to p53-mediated mitochondrial and neuronal damage in AD. Therapeutic modulation of SIRT3 activity may ameliorate mitochondrial pathology and neurodegeneration in AD.

Robust Production, Crystallization, Structure Determination, and Analysis of [Fe-S] Proteins: Uncovering Control of Electron Shuttling and Gating in the Respiratory Metabolism of Molybdopterin Guanine Dinucleotide Enzymes

[Fe-S] clusters are essential cofactors in all domains of life. They play many biological roles due to their unique abilities for electron transfer and conformational control. Yet, producing and analyzing Fe-S proteins can be difficult and even misleading if not done anaerobically. Due to unique redox properties of [Fe-S] clusters and their oxygen sensitivity, they pose multiple challenges and can lose enzymatic activity or cause their component proteins to be structurally disordered due to [Fe-S] cluster oxidation and loss in air. Here we highlight tested protocols and strategies enabling efficient and stable [Fe-S] protein production, purification, crystallization, X-ray diffraction data collection, and structure determination. From multiple high-resolution anaerobic crystal structures, we furthermore analyze exemplary data defining [Fe-S] clusters, substrate entry, and product exit for the functional oxidation states of type II molybdo-bis(molybdopterin guanine dinucleotide) (Mo-bisMGD) enzymes. Notably, these enzymes perform electron shuttling between quinone pools and specific substrates to catalyze respiratory metabolism. The identified structure-activity relationships for this enzyme class have broad implications germane to perchlorate environments on Earth and Mars extending to an alternative mechanism underlying metabolic origins for the evolution of the oxygen atmosphere. Integrated structural analyses of type II Mo-bisMGD enzymes unveil novel distinctive shared molecular mechanisms for dynamic control of substrate entry and product release gated by hydrophobic residues. Collective findings support a prototypic model for type II Mo-bisMGD enzymes including insights for a fundamental molecular mechanistic understanding of selectivity and regulation by a conformationally gated channel with general implications for [Fe-S] cluster respiratory enzymes.

Observation of dielectric universalities in albumin, cytochrome C and Shewanella oneidensis MR-1 extracellular matrix

The electrodynamics of metals is well understood within the Drude conductivity model; properties of insulators and semiconductors are governed by a gap in the electronic states. But there is a great variety of disordered materials that do not fall in these categories and still respond to external field in an amazingly uniform manner. At radiofrequencies delocalized charges yield a frequency-independent conductivity σ 1(ν) whose magnitude exponentially decreases while cooling. With increasing frequency, dispersionless conductivity starts to reveal a power-law dependence σ 1(ν)∝ν s with s < 1 caused by hopping charge carriers. At low temperatures, such Universal Dielectric Response can cross over to another universal regime with nearly constant loss ε″∝σ1/ν = const. The powerful research potential based on such universalities is widely used in condensed matter physics. Here we study the broad-band (1-1012 Hz) dielectric response of Shewanella oneidensis MR-1 extracellular matrix, cytochrome C and serum albumin. Applying concepts of condensed matter physics, we identify transport mechanisms and a number of energy, time, frequency, spatial and temperature scales in these biological objects, which can provide us with deeper insight into the protein dynamics.

Electron-transfer chain in respiratory complex I

Complex I is a part of the respiration energy chain converting the redox energy into the cross-membrane proton gradient. The electron-transfer chain of iron-sulfur cofactors within the water-soluble peripheral part of the complex is responsible for the delivery of electrons to the proton pumping subunit. The protein is porous to water penetration and the hydration level of the cofactors changes when the electron is transferred along the chain. High reaction barriers and trapping of the electrons at the iron-sulfur cofactors are prevented by the combination of intense electrostatic noise produced by the protein-water interface with the high density of quantum states in the iron-sulfur clusters caused by spin interactions between paramagnetic iron atoms. The combination of these factors substantially lowers the activation barrier for electron transfer compared to the prediction of the Marcus theory, bringing the rate to the experimentally established range. The unique role of iron-sulfur clusters as electron-transfer cofactors is in merging protein-water fluctuations with quantum-state multiplicity to allow low activation barriers and robust operation. Water plays a vital role in electron transport energetics by electrowetting the cofactors in the chain upon arrival of the electron. A general property of a protein is to violate the fluctuation-dissipation relation through nonergodic sampling of its landscape. High functional efficiency of redox enzymes is a direct consequence of nonergodicity.

Self-assembly of an electronically conductive network through microporous scaffolds

Electron transfer spanning significant distances through a microporous structure was established via the self-assembly of an electronically conductive iridium oxide nanowire matrix enveloping the pore walls. Microporous formations were simulated using two scaffold materials of varying physical and chemical properties; paraffin wax beads, and agar gel. Following infiltration into the micropores, iridium nanoparticles self-assembled at the pore wall/ethanol interface. Subsequently, cyclic voltammetry was employed to electrochemically crosslink the metal, erecting an interconnected, and electronically conductive metal oxide nanowire matrix. Electrochemical and spectral characterization techniques confirmed the formation of oxide nanowire matrices encompassing lengths of at least 1.6mm, 400× distances previously achieved using iridium nanoparticles. Nanowire matrices were engaged as biofuel cell anodes, where electrons were donated to the nanowires by a glucose oxidizing enzyme.

Iron-Sulfur Clusters in Mitochondrial Metabolism: Multifaceted Roles of a Simple Cofactor

Iron-sulfur metabolism is essential for cellular function and is a key process in mitochondria. In this review, we focus on the structure and assembly of mitochondrial iron-sulfur clusters and their roles in various metabolic processes that occur in mitochondria. Iron-sulfur clusters are crucial in mitochondrial respiration, in which they are required for the assembly, stability, and function of respiratory complexes I, II, and III. They also serve important functions in the citric acid cycle, DNA metabolism, and apoptosis. Whereas the identification of iron-sulfur containing proteins and their roles in numerous aspects of cellular function has been a long-standing research area, that in mitochondria is comparatively recent, and it is likely that their roles within mitochondria have been only partially revealed. We review the status of the field and provide examples of other cellular iron-sulfur proteins to highlight their multifarious roles.

Temperature-dependent ESR and computational studies on antiferromagnetic electron transfer in the yeast NADH dehydrogenase Ndi1

The electron transfer mechanism of Ndi1.

Polarizability of the active site of cytochrome c reduces the activation barrier for electron transfer

Enzymes in biology’s energy chains operate with low energy input distributed through multiple electron transfer steps between protein active sites. The general challenge of biological design is how to lower the activation barrier without sacrificing a large negative reaction free energy. We show that this goal is achieved through a large polarizability of the active site. It is polarized by allowing a large number of excited states, which are populated quantum mechanically by electrostatic fluctuations of the protein and hydration water shells. This perspective is achieved by extensive mixed quantum mechanical/molecular dynamics simulations of the half reaction of reduction of cytochrome c. The barrier for electron transfer is consistently lowered by increasing the number of excited states included in the Hamiltonian of the active site diagonalized along the classical trajectory. We suggest that molecular polarizability, in addition to much studied electrostatics of permanent charges, is a key parameter to consider in order to understand how enzymes work.

Electron tunneling in proteins program

We developed a unique integrated software package (called Electron Tunneling in Proteins Program or ETP) which provides an environment with different capabilities such as tunneling current calculation, semi-empirical quantum mechanical calculation, and molecular modeling simulation for calculation and analysis of electron transfer reactions in proteins. ETP program is developed as a cross-platform client-server program in which all the different calculations are conducted at the server side while only the client terminal displays the resulting calculation outputs in the different supported representations. ETP program is integrated with a set of well-known computational software packages including Gaussian, BALLVIEW, Dowser, pKip, and APBS. In addition, ETP program supports various visualization methods for the tunneling calculation results that assist in a more comprehensive understanding of the tunneling process. © 2016 Wiley Periodicals, Inc.

Novel Inhibitors for a Novel Binding Site in Respiratory Complex III

A new binding site and potential novel inhibitors of the respiratory complex III are described. The site is located at the opposite side of the enzyme with respect to ubiquinol binding site (Qo site), and distinctly different from both Qo and Qi sites (hence designated as Non-Q binding site, NQ). NQ site binding pocket extends up close to Phe90 residue, an internal switch (LH switch) that regulates electron transfer between heme bL and heme bH of the low potential redox chain. Docking studies and molecular dynamics simulations of different molecules to the NQ site revealed potential ligands which exhibit a novel inhibitory effect for bc1 complex by switching the LH switch to "off" conformation, thereby shutting down electron transfer in the low potential redox chain. Moreover, the novel inhibitors have lower binding affinity for both Qo and Qi sites, and hence do not interfere with binding of the natural ligands to those sites. The inhibitory activity of those novel ligands in bc1 complex is suggested to promote the production of reactive oxygen species (ROS) at the Qo site. Hence those ligands are potential candidates for designing new "mitocan" drugs.

The multitude of iron-sulfur clusters in respiratory complex I

Respiratory complex I couples the electron transfer from NADH to ubiquinone with the translocation of protons across the membrane. Complex I contains one non-covalently bound flavin mononucleotide and, depending on the species, up to ten iron-sulfur (Fe/S) clusters as cofactors. The reason for the presence of the multitude of Fe/S clusters in complex I remained enigmatic for a long time. The question was partly answered by investigations on the evolution of the complex revealing the stepwise construction of the electron transfer domain from several modules. Extension of the ancestral to the modern electron input domain was associated with the acquisition of several Fe/S-proteins. The X-ray structure of the complex showed that the NADH oxidation-site is connected with the quinone-reduction site by a chain of seven Fe/S-clusters. Fast enzyme kinetics revealed that this chain of Fe/S-clusters is used to regulate electron-tunneling rates within the complex. A possible function of the off-pathway cluster N1a is discussed. This article is part of a Special Issue entitled 'EBEC 2016: 19th European Bioenergetics Conference, Riva del Garda, Italy, July 2-6, 2016', edited by Prof. Paolo Bernardi.

Molecular simulation and modeling of complex I

Molecular modeling and molecular dynamics simulations play an important role in the functional characterization of complex I. With its large size and complicated function, linking quinone reduction to proton pumping across a membrane, complex I poses unique modeling challenges. Nonetheless, simulations have already helped in the identification of possible proton transfer pathways. Simulations have also shed light on the coupling between electron and proton transfer, thus pointing the way in the search for the mechanistic principles underlying the proton pump. In addition to reviewing what has already been achieved in complex I modeling, we aim here to identify pressing issues and to provide guidance for future research to harness the power of modeling in the functional characterization of complex I. This article is part of a Special Issue entitled Respiratory complex I, edited by Volker Zickermann and Ulrich Brandt.

Genome-scale metabolic reconstructions and theoretical investigation of methane conversion in Methylomicrobium buryatense strain 5G(B1)

BACKGROUND

Methane-utilizing bacteria (methanotrophs) are capable of growth on methane and are attractive systems for bio-catalysis. However, the application of natural methanotrophic strains to large-scale production of value-added chemicals/biofuels requires a number of physiological and genetic alterations. An accurate metabolic model coupled with flux balance analysis can provide a solid interpretative framework for experimental data analyses and integration.

RESULTS

A stoichiometric flux balance model of Methylomicrobium buryatense strain 5G(B1) was constructed and used for evaluating metabolic engineering strategies for biofuels and chemical production with a methanotrophic bacterium as the catalytic platform. The initial metabolic reconstruction was based on whole-genome predictions. Each metabolic step was manually verified, gapfilled, and modified in accordance with genome-wide expression data. The final model incorporates a total of 841 reactions (in 167 metabolic pathways). Of these, up to 400 reactions were recruited to produce 118 intracellular metabolites. The flux balance simulations suggest that only the transfer of electrons from methanol oxidation to methane oxidation steps can support measured growth and methane/oxygen consumption parameters, while the scenario employing NADH as a possible source of electrons for particulate methane monooxygenase cannot. Direct coupling between methane oxidation and methanol oxidation accounts for most of the membrane-associated methane monooxygenase activity. However the best fit to experimental results is achieved only after assuming that the efficiency of direct coupling depends on growth conditions and additional NADH input (about 0.1-0.2 mol of incremental NADH per one mol of methane oxidized). The additional input is proposed to cover loss of electrons through inefficiency and to sustain methane oxidation at perturbations or support uphill electron transfer. Finally, the model was used for testing the carbon conversion efficiency of different pathways for C1-utilization, including different variants of the ribulose monophosphate pathway and the serine cycle.

CONCLUSION

We demonstrate that the metabolic model can provide an effective tool for predicting metabolic parameters for different nutrients and genetic perturbations, and as such, should be valuable for metabolic engineering of the central metabolism of M. buryatense strains.

A theoretical multiscale treatment of protein-protein electron transfer: The ferredoxin/ferredoxin-NADP(+) reductase and flavodoxin/ferredoxin-NADP(+) reductase systems

In the photosynthetic electron transfer (ET) chain, two electrons transfer from photosystem I to the flavin-dependent ferredoxin-NADP(+) reductase (FNR) via two sequential independent ferredoxin (Fd) electron carriers. In some algae and cyanobacteria (as Anabaena), under low iron conditions, flavodoxin (Fld) replaces Fd as single electron carrier. Extensive mutational studies have characterized the protein-protein interaction in FNR/Fd and FNR/Fld complexes. Interestingly, even though Fd and Fld share the interaction site on FNR, individual residues on FNR do not participate to the same extent in the interaction with each of the protein partners, pointing to different electron transfer mechanisms. Despite of extensive mutational studies, only FNR/Fd X-ray structures from Anabaena and maize have been solved; structural data for FNR/Fld remains elusive. Here, we present a multiscale modelling approach including coarse-grained and all-atom protein-protein docking, the QM/MM e-Pathway analysis and electronic coupling calculations, allowing for a molecular and electronic comprehensive analysis of the ET process in both complexes. Our results, consistent with experimental mutational data, reveal the ET in FNR/Fd proceeding through a bridge-mediated mechanism in a dominant protein-protein complex, where transfer of the electron is facilitated by Fd loop-residues 40-49. In FNR/Fld, however, we observe a direct transfer between redox cofactors and less complex specificity than in Fd; more than one orientation in the encounter complex can be efficient in ET.