Charge Transfer and Chemo-Mechanical Coupling in Respiratory Complex I

A computational study to assess the pathogenicity of single or combinations of missense variants on respiratory complex I

Variants found in the respiratory complex I (CI) subunit genes encoded by mitochondrial DNA can cause severe genetic diseases. However, it is difficult to establish a priori whether a single or a combination of CI variants may impact oxidative phosphorylation. Here we propose a computational approach based on coarse-grained molecular dynamics simulations aimed at investigating new CI variants. One of the primary CI variants associated with the Leber hereditary optic neuropathy (m.14484T>C/MT-ND6) was used as a test case and was investigated alone or in combination with two additional rare CI variants whose role remains uncertain. We found that the primary variant positioned in the E-channel region, which is fundamental for CI function, stiffens the enzyme dynamics. Moreover, a new mechanism for the transition between π- and α-conformation in the helix carrying the primary variant is proposed. This may have implications for the E-channel opening/closing mechanism. Finally, our findings show that one of the rare variants, located next to the primary one, further worsens the stiffening, while the other rare variant does not affect CI function. This approach may be extended to other variants candidate to exert a pathogenic impact on CI dynamics, or to investigate the interaction of multiple variants.

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

Long-range electron proton coupling in respiratory complex I - insights from molecular simulations of the quinone chamber and antiporter-like subunits

Respiratory complex I is a redox-driven proton pump. Several high-resolution structures of complex I have been determined providing important information about the putative proton transfer paths and conformational transitions that may occur during catalysis. However, how redox energy is coupled to the pumping of protons remains unclear. In this article, we review biochemical, structural and molecular simulation data on complex I and discuss several coupling models, including the key unresolved mechanistic questions. Focusing both on the quinone-reductase domain as well as the proton-pumping membrane-bound domain of complex I, we discuss a molecular mechanism of proton pumping that satisfies most experimental and theoretical constraints. We suggest that protonation reactions play an important role not only in catalysis, but also in the physiologically-relevant active/deactive transition of complex I.

Quinone Catalysis Modulates Proton Transfer Reactions in the Membrane Domain of Respiratory Complex I

Complex I is a redox-driven proton pump that drives electron transport chains and powers oxidative phosphorylation across all domains of life. Yet, despite recently resolved structures from multiple organisms, it still remains unclear how the redox reactions in Complex I trigger proton pumping up to 200 Å away from the active site. Here, we show that the proton-coupled electron transfer reactions during quinone reduction drive long-range conformational changes of conserved loops and trans-membrane (TM) helices in the membrane domain of Complex I from Yarrowia lipolytica. We find that the conformational switching triggers a π → α transition in a TM helix (TM3ND6) and establishes a proton pathway between the quinone chamber and the antiporter-like subunits, responsible for proton pumping. Our large-scale (>20 μs) atomistic molecular dynamics (MD) simulations in combination with quantum/classical (QM/MM) free energy calculations show that the helix transition controls the barrier for proton transfer reactions by wetting transitions and electrostatic effects. The conformational switching is enabled by re-arrangements of ion pairs that propagate from the quinone binding site to the membrane domain via an extended network of conserved residues. We find that these redox-driven changes create a conserved coupling network within the Complex I superfamily, with point mutations leading to drastic activity changes and mitochondrial disorders. On a general level, our findings illustrate how catalysis controls large-scale protein conformational changes and enables ion transport across biological membranes.

Fostering discoveries in the era of exascale computing: How the next generation of supercomputers empowers computational and experimental biophysics alike

Over a century ago, physicists started broadly relying on theoretical models to guide new experiments. Soon thereafter, chemists began doing the same. Now, biological research enters a new era when experiment and theory walk hand in hand. Novel software and specialized hardware became essential to understand experimental data and propose new models. In fact, current petascale computing resources already allow researchers to reach unprecedented levels of simulation throughput to connect in silico and in vitro experiments. The reduction in cost and improved access allowed a large number of research groups to adopt supercomputing resources and techniques. Here, we outline how large-scale computing has evolved to expand decades-old research, spark new research efforts, and continuously connect simulation and observation. For instance, multiple publicly and privately funded groups have dedicated extensive resources to develop artificial intelligence tools for computational biophysics, from accelerating quantum chemistry calculations to proposing protein structure models. Moreover, advances in computer hardware have accelerated data processing from single-molecule experimental observations and simulations of chemical reactions occurring throughout entire cells. The combination of software and hardware has opened the way for exascale computing and the production of the first public exascale supercomputer, Frontier, inaugurated by the Oak Ridge National Laboratory in 2022. Ultimately, the popularization and development of computational techniques and the training of researchers to use them will only accelerate the diversification of tools and learning resources for future generations.

Dynamic conformational switching underlies TFIIH function in transcription and DNA repair and impacts genetic diseases

Transcription factor IIH (TFIIH) is a protein assembly essential for transcription initiation and nucleotide excision repair (NER). Yet, understanding of the conformational switching underpinning these diverse TFIIH functions remains fragmentary. TFIIH mechanisms critically depend on two translocase subunits, XPB and XPD. To unravel their functions and regulation, we build cryo-EM based TFIIH models in transcription- and NER-competent states. Using simulations and graph-theoretical analysis methods, we reveal TFIIH's global motions, define TFIIH partitioning into dynamic communities and show how TFIIH reshapes itself and self-regulates depending on functional context. Our study uncovers an internal regulatory mechanism that switches XPB and XPD activities making them mutually exclusive between NER and transcription initiation. By sequentially coordinating the XPB and XPD DNA-unwinding activities, the switch ensures precise DNA incision in NER. Mapping TFIIH disease mutations onto network models reveals clustering into distinct mechanistic classes, affecting translocase functions, protein interactions and interface dynamics.

On the role of ubiquinone in the proton translocation mechanism of respiratory complex I

Complex I converts oxidoreduction energy into a proton electrochemical gradient across the inner mitochondrial or bacterial cell membrane. This gradient is the primary source of energy for aerobic synthesis of ATP. Oxidation of reduced nicotinamide adenine dinucleotide (NADH) by ubiquinone (Q) yields NAD⁺ and ubiquinol (QH₂ ), which is tightly coupled to translocation of four protons from the negatively to the positively charged side of the membrane. Electrons from NADH oxidation reach the iron-sulfur centre N2 positioned near the bottom of a tunnel that extends circa 30 Å from the membrane domain into the hydrophilic domain of the complex. The tunnel is occupied by ubiquinone, which can take a distal position near the N2 centre or proximal positions closer to the membrane. Here, we review important structural, kinetic and thermodynamic properties of ubiquinone that define its role in complex I function. We suggest that this function exceeds that of a mere substrate or electron acceptor and propose that ubiquinone may be the redox element of complex I coupling electron transfer to proton translocation.

Structure of respiratory complex I – An emerging blueprint for the mechanism

Complex I is one of the major respiratory complexes, conserved from bacteria to mammals. It oxidises NADH, reduces quinone and pumps protons across the membrane, thus playing a central role in the oxidative energy metabolism. In this review we discuss our current state of understanding the structure of complex I from various species of mammals, plants, fungi, and bacteria, as well as of several complex I-related proteins. By comparing the structural evidence from these systems in different redox states and data from mutagenesis and molecular simulations, we formulate the mechanisms of electron transfer and proton pumping and explain how they are conformationally and electrostatically coupled. Finally, we discuss the structural basis of the deactivation phenomenon in mammalian complex I.

Electron Transport Lipids Fold Within Membrane-Like Interfaces

Lipoquinones, such as ubiquinones (UQ) and menaquinones (MK), function as essential lipid components of the electron transport system (ETS) by shuttling electrons and protons to facilitate the production of ATP in eukaryotes and prokaryotes. Lipoquinone function in membrane systems has been widely studied, but the exact location and conformation within membranes remains controversial. Lipoquinones, such as Coenzyme Q (UQ-10), are generally depicted simply as "Q" in life science diagrams or in extended conformations in primary literature even though specific conformations are important for function in the ETS. In this study, our goal was to determine the location, orientation, and conformation of UQ-2, a truncated analog of UQ-10, in model membrane systems and to compare our results to previously studied MK-2. Herein, we first carried out a six-step synthesis to yield UQ-2 and then demonstrated that UQ-2 adopts a folded conformation in organic solvents using 1H-1H 2D NOESY and ROESY NMR spectroscopic studies. Similarly, using 1H-1H 2D NOESY NMR spectroscopic studies, UQ-2 was found to adopt a folded, U-shaped conformation within the interface of an AOT reverse micelle model membrane system. UQ-2 was located slightly closer to the surfactant-water interface compared to the more hydrophobic MK-2. In addition, Langmuir monolayer studies determined UQ-2 resided within the monolayer water-phospholipid interface causing expansion, whereas MK-2 was more likely to be compressed out and reside within the phospholipid tails. All together these results support the model that lipoquinones fold regardless of the headgroup structure but that the polarity of the headgroup influences lipoquinone location within the membrane interface. These results have implications regarding the redox activity near the interface as quinone vs. quinol forms may facilitate locomotion of lipoquinones within the membrane. The location, orientation, and conformation of lipoquinones are critical for their function in generating cellular energy within membrane ETS, and the studies described herein shed light on the behavior of lipoquinones within membrane-like environments.

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.

Structure of the peripheral arm of a minimalistic respiratory complex I

Respiratory complex I drives proton translocation across energy-transducing membranes by NADH oxidation coupled with (ubi)quinone reduction. In humans, its dysfunction is associated with neurodegenerative diseases. The Escherichia coli complex represents the structural minimal form of an energy-converting NADH:ubiquinone oxidoreductase. Here, we report the structure of the peripheral arm of the E. coli complex I consisting of six subunits, the FMN cofactor, and nine iron-sulfur clusters at 2.7 Å resolution obtained by cryo electron microscopy. While the cofactors are in equivalent positions as in the complex from other species, individual subunits are adapted to the absence of supernumerary proteins to guarantee structural stability. The catalytically important subunits NuoC and D are fused resulting in a specific architecture of functional importance. Striking features of the E. coli complex are scrutinized by mutagenesis and biochemical characterization of the variants. Moreover, the arrangement of the subunits sheds light on the unknown assembly of the complex.

Interactions of Truncated Menaquinones in Lipid Monolayers and Bilayers

Menaquinones (MK) are hydrophobic molecules that consist of a naphthoquinone headgroup and a repeating isoprenyl side chain and are cofactors used in bacterial electron transport systems to generate cellular energy. We have previously demonstrated that the folded conformation of truncated MK homologues, MK-1 and MK-2, in both solution and reverse micelle microemulsions depended on environment. There is little information on how MKs associate with phospholipids in a model membrane system and how MKs affect phospholipid organization. In this manuscript, we used a combination of Langmuir monolayer studies and molecular dynamics (MD) simulations to probe these questions on truncated MK homologues, MK-1 through MK-4 within a model membrane. We observed that truncated MKs reside farther away from the interfacial water than ubiquinones are are located closer to the phospholipid tails. We also observed that phospholipid packing does not change at physiological pressure in the presence of truncated MKs, though a difference in phospholipid packing has been observed in the presence of ubiquinones. We found through MD simulations that for truncated MKs, the folded conformation varied, but MKs location and association with the bilayer remained unchanged at physiological conditions regardless of side chain length. Combined, this manuscript provides fundamental information, both experimental and computational, on the location, association, and conformation of truncated MK homologues in model membrane environments relevant to bacterial energy production.

Poor Person's pH Simulation of Membrane Proteins

pH conditions are central to the functioning of all biomolecules. However, implications of pH changes are nontrivial on a molecular scale. Though a rigorous microscopic definition of pH exists, its implementation in classical molecular dynamics (MD) simulations is cumbersome, and more so in large integral membrane systems. In this chapter, an integrative pipeline is described that combines Multi-Conformation Continuum Electrostatics (MCCE) computations with MD simulations to capture the effect of transient protonation states on the coupled conformational changes in transmembrane proteins. The core methodologies are explained, and all the software required to set up this pipeline are outlined with their key parameters. All associated analyses of structure and function are provided using two case studies, namely those of bioenergetic complexes: NADH dehydrogenase (complex I) and V_o domain of V-type ATPase. The hybrid MCCE-MD pipeline has allowed the discovery of hydrogen bond networks, ligand binding pathways, and disease-causing mutations.

Deactivation blocks proton pathways in the mitochondrial complex I

The electron transport chain of mitochondria is initiated by the respiratory complex I that converts chemical energy into a proton motive force to power synthesis of adenosine triphosphate. On a chemical level, complex I catalyzes elementary electron and proton transfer processes that couple across large molecular distances of >300 Å. However, under low oxygen concentrations, the respiratory chain operates in reverse mode and produces harmful reactive oxygen species. To avoid cell damage, the mitochondrial complex I transitions into a deactive state that inhibits turnover by molecular principles that remain elusive. By combining large-scale molecular simulations with cryo-electron microscopy data, we show here that complex I deactivation blocks the communication between proton pumping and redox modules by conformational and hydration changes.

Cellular respiration is powered by membrane-bound redox enzymes that convert chemical energy into an electrochemical proton gradient and drive the energy metabolism. By combining large-scale classical and quantum mechanical simulations with cryo-electron microscopy data, we resolve here molecular details of conformational changes linked to proton pumping in the mammalian complex I. Our data suggest that complex I deactivation blocks water-mediated proton transfer between a membrane-bound quinone site and proton-pumping modules, decoupling the energy-transduction machinery. We identify a putative gating region at the interface between membrane domain subunits ND1 and ND3/ND4L/ND6 that modulates the proton transfer by conformational changes in transmembrane helices and bulky residues. The region is perturbed by mutations linked to human mitochondrial disorders and is suggested to also undergo conformational changes during catalysis of simpler complex I variants that lack the “active”-to-“deactive” transition. Our findings suggest that conformational changes in transmembrane helices modulate the proton transfer dynamics by wetting/dewetting transitions and provide important functional insight into the mammalian respiratory complex I.

Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient

Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.

A conserved arginine residue is critical for stabilizing the N2 FeS cluster in mitochondrial complex I

Respiratory complex I (NADH:ubiquinone oxidoreductase), the first enzyme of the electron-transport chain, captures the free energy released by NADH oxidation and ubiquinone reduction to translocate protons across an energy-transducing membrane and drive ATP synthesis during oxidative phosphorylation. The cofactor that transfers the electrons directly to ubiquinone is an iron–sulfur cluster (N2) located in the NDUFS2/NUCM subunit. A nearby arginine residue (R121), which forms part of the second coordination sphere of the N2 cluster, is known to be posttranslationally dimethylated but its functional and structural significance are not known. Here, we show that mutations of this arginine residue (R121M/K) abolish the quinone-reductase activity, concomitant with disappearance of the N2 signature from the electron paramagnetic resonance (EPR) spectrum. Analysis of the cryo-EM structure of NDUFS2-R121M complex I at 3.7 Å resolution identified the absence of the cubane N2 cluster as the cause of the dysfunction, within an otherwise intact enzyme. The mutation further induced localized disorder in nearby elements of the quinone-binding site, consistent with the close connections between the cluster and substrate-binding regions. Our results demonstrate that R121 is required for the formation and/or stability of the N2 cluster and highlight the importance of structural analyses for mechanistic interpretation of biochemical and spectroscopic data on complex I variants.

Is the NDUFV2 subunit of the hydrophilic complex I domain a key determinant of animal longevity?

Complex I, a component of the electron transport chain, plays a central functional role in cell bioenergetics and the biology of free radicals. The structural and functional N module of complex I is one of the main sites of the generation of free radicals. The NDUFV2 subunit/N1a cluster is a component of this module. Furthermore, the rate of free radical production is linked to animal longevity. In this review, we explore the hypothesis that NDUFV2 is the only conserved core subunit designed with a regulatory function to ensure correct electron transfer and free radical production, that low gene expression and protein abundance of the NDUFV2 subunit is an evolutionary adaptation needed to achieve a longevity phenotype, and that these features are determinants of the lower free radical generation at the mitochondrial level and a slower rate of aging of long-lived animals.

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.

Coupling of quinone dynamics to proton pumping in respiratory complex I

Respiratory complex I (NADH:quinone oxidoreductase) plays a central role in generating the proton electrochemical gradient in mitochondrial and bacterial membranes, which is needed to generate ATP. Several high-resolution structures of complex I have been determined, revealing its intricate architecture and complementing the biochemical and biophysical studies. However, the molecular mechanism of long-range coupling between ubiquinone (Q) reduction and proton pumping is not known. Computer simulations have been applied to decipher the dynamics of Q molecule in the ~30 Å long Q tunnel. In this short report, we discuss the binding and dynamics of Q at computationally predicted Q binding sites, many of which are supported by structural data on complex I. We suggest that the binding of Q at these sites is coupled to proton pumping by means of conformational rearrangements in the conserved loops of core subunits.

Molecular Mechanism by which Cobra Venom Cardiotoxins Interact with the Outer Mitochondrial Membrane

Cardiotoxin CTII from Najaoxiana cobra venom translocates to the intermembrane space (IMS) of mitochondria to disrupt the structure and function of the inner mitochondrial membrane. At low concentrations, CTII facilitates ATP-synthase activity, presumably via the formation of non-bilayer, immobilized phospholipids that are critical in modulating ATP-synthase activity. In this study, we investigated the effects of another cardiotoxin CTI from Najaoxiana cobra venom on the structure of mitochondrial membranes and on mitochondrial-derived ATP synthesis. By employing robust biophysical methods including 31P-NMR and 1H-NMR spectroscopy, we analyzed the effects of CTI and CTII on phospholipid packing and dynamics in model phosphatidylcholine (PC) membranes enriched with 2.5 and 5.0 mol% of cardiolipin (CL), a phospholipid composition that mimics that in the outer mitochondrial membrane (OMM). These experiments revealed that CTII converted a higher percentage of bilayer phospholipids to a non-bilayer and immobilized state and both cardiotoxins utilized CL and PC molecules to form non-bilayer structures. Furthermore, in order to gain further understanding on how cardiotoxins bind to mitochondrial membranes, we employed molecular dynamics (MD) and molecular docking simulations to investigate the molecular mechanisms by which CTII and CTI interactively bind with an in silico phospholipid membrane that models the composition similar to the OMM. In brief, MD studies suggest that CTII utilized the N-terminal region to embed the phospholipid bilayer more avidly in a horizontal orientation with respect to the lipid bilayer and thereby penetrate at a faster rate compared with CTI. Molecular dynamics along with the Autodock studies identified critical amino acid residues on the molecular surfaces of CTII and CTI that facilitated the long-range and short-range interactions of cardiotoxins with CL and PC. Based on our compiled data and our published findings, we provide a conceptual model that explains a molecular mechanism by which snake venom cardiotoxins, including CTI and CTII, interact with mitochondrial membranes to alter the mitochondrial membrane structure to either upregulate ATP-synthase activity or disrupt mitochondrial function.

Hydrogen bond network analysis reveals the pathway for the proton transfer in the E-channel of T. thermophilus Complex I

Complex I, NADH-ubiquinone oxidoreductase, is the first enzyme in the mitochondrial and bacterial aerobic respiratory chain. It pumps four protons through four transiently open pathways from the high pH, negative, N-side of the membrane to the positive, P-side driven by the exergonic transfer of electrons from NADH to a quinone. Three protons transfer through subunits descended from antiporters, while the fourth, E-channel is unique. The path through the E-channel is determined by a network analysis of hydrogen bonded pathways obtained by Monte Carlo sampling of protonation states, polar hydrogen orientation and water occupancy. Input coordinates are derived from molecular dynamics trajectories comparing oxidized, reduced (dihydro) and no menaquinone-8 (MQ). A complex proton transfer path from the N- to the P-side is found consisting of six clusters of highly connected hydrogen-bonded residues. The network connectivity depends on the presence of quinone and its redox state, supporting a role for this cofactor in coupling electron and proton transfers. The N-side is more organized with MQ-bound complex I facilitating proton entry, while the P-side is more connected in the apo-protein, facilitating proton exit. Subunit Nqo8 forms the core of the E channel; Nqo4 provides the N-side entry, Nqo7 and then Nqo10 join the pathway in the middle, while Nqo11 contributes to the P-side exit.