The respiratory complex I transduces redox energy into an electrochemical proton gradient in aerobic respiratory chains, powering energy-requiring processes in the cell. However, despite recently resolved molecular structures, the mechanism of this gigantic enzyme remains poorly understood. By combining large-scale quantum and classical simulations with site-directed mutagenesis and biophysical experiments, we show here how the conformational state of buried ion-pairs and water molecules control the protonation dynamics in the membrane domain of complex I and establish evolutionary conserved long-range coupling elements. We suggest that an electrostatic wave propagates in forward and reverse directions across the 200 Å long membrane domain during enzyme turnover, without significant dissipation of energy. Our findings demonstrate molecular principles that enable efficient long-range proton–electron coupling (PCET) and how perturbation of this PCET machinery may lead to development of mitochondrial disease.

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.

Photosynthetic organisms capture light energy to drive their energy metabolism, and employ the chemical reducing power to convert carbon dioxide (CO 2 ) into organic molecules. Photorespiration, however, significantly reduces the photosynthetic yields. To survive under low CO 2 concentrations, cyanobacteria evolved unique carbon-concentration mechanisms that enhance the efficiency of photosynthetic CO 2 fixation, for which the molecular principles have remained unknown. We show here how modular adaptations enabled the cyanobacterial photosynthetic complex I to concentrate CO 2 using a redox-driven proton-pumping machinery. Our cryo-electron microscopy structure at 3.2 Å resolution shows a catalytic carbonic anhydrase module that harbours a Zn 2+ active site, with connectivity to protonpumping subunits that are activated by electron transfer from photosystem I. Our findings illustrate molecular principles in the photosynthetic complex I machinery that enabled cyanobacteria to survive in drastically changing CO 2 conditions.

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.

The regioselective activation of C–H bonds in complex molecules containing several of them is still an exciting challenge in chemistry. However, many enzymes catalyze these processes, and much can be learned from the way they do it. For example, each mammalian lipoxygenase isoenzyme abstracts a hydrogen atom from essentially a unique carbon position. This paper presents a comprehensive study at the B3LYP­(QM)/CHARMM­(MM) level of the hydrogen abstraction reaction from arachidonic acid (AA) catalyzed by rabbit 15-lipoxygenase (15-rLO). Most of the products of this reaction arise from the initial hydrogen abstraction from the carbon C13 of AA. Nevertheless, we have shown that 15-rLO seems able to catalyze not only the abstraction of H13 but also that of H10. After having studied 20 of these hydrogen transfers initiated from different snapshots of an extensive MD sampling of the Michaelis complex, we have even concluded that the reaction mechanisms for both abstractions are identical (proton-coupled electron transfer processes), with transition state structures matching their geometries around the shifting hydrogen. Despite that similarity, the average potential energy barrier for the H13 abstractions is 4.0 kcal/mol lower than for the H10 abstractions, in good agreement with the experimental C13:C10 ratio of 97:3. We have found that a subtle steric hindrance by Leu597 and Ile663 is the main cause for that difference. Driving the strict regiospecificity exhibited by 15-rLO appears to be the essential function of the bulky side chains of those conserved residues, in this way making possible the vital physiological role of 15-rLO and, probably, of all the mammalian lipoxygenase isoenzymes. The understanding of how nature uses residues with the bulkiest aliphatic side chains to achieve the selective activation of C–H bonds can stimulate the design of efficient biocatalysts to that aim.

ALOX15 (12/15-lipoxygenase) orthologs have been implicated in maturational degradation of intracellular organelles and in the biosynthesis of antiinflammatory and proresolving eicosanoids. Here we hypothesized that lower mammals (mice, rats, pigs) express 12-lipoxygenating ALOX15 orthologs. In contrast, 15-lipoxygenating isoforms are found in higher primates (orangutans, men), and these results suggest an evolution of ALOX15 specificity. To test this hypothesis we first cloned and characterized ALOX15 orthologs of selected Catarrhini representing different stages of late primate evolution and found that higher primates (men, chimpanzees) express 15-lipoxygenating orthologs. In contrast, lower primates (baboons, rhesus monkeys) express 12-lipoxygenating enzymes. Gibbons, which are flanked in evolution by rhesus monkeys (12-lipoxygenating ALOX15) and orangutans (15-lipoxygenating ALOX15), express an ALOX15 ortholog with pronounced dual specificity. To explore the driving force for this evolutionary alterations, we quantified the lipoxin synthase activity of 12-lipoxygenating (rhesus monkey, mouse, rat, pig, humIle418Ala) and 15-lipoxygenating (man, chimpanzee, orangutan, rabbit, ratPhe353Ala) ALOX15 variants and found that, when normalized to their arachidonic acid oxygenase activities, the lipoxin synthase activities of 15-lipoxygenating ALOX15 variants were more than fivefold higher (P < 0.01). Comparative molecular dynamics simulations and quantum mechanics/molecular mechanics calculations indicated that, for the 15-lipoxygenating rabbit ALOX15, the energy barrier for C13-hydrogen abstraction (15-lipoxygenation) was 17 kJ/mol lower than for arachidonic acid 12-lipoxygenation. In contrast, for the 12-lipoxygenating Ile418Ala mutant, the energy barrier for 15-lipoxygenation was 10 kJ/mol higher than for 12-lipoxygenation. Taken together, our data suggest an evolution of ALOX15 specificity, which is aimed at optimizing the biosynthetic capacity for antiinflammatory and proresolving lipoxins.

Complex I functions as an initial electron acceptor in aerobic respiratory chains that reduces quinone and pumps protons across a biological membrane. This remarkable charge transfer process extends ca. 300 Å and it is initiated by a poorly understood proton-coupled electron transfer (PCET) reaction between nicotinamide adenine dinucleotide (NADH) and a protein-bound flavin (FMN) cofactor. We combine here large-scale density functional theory calculations and quantum/classical models with atomistic molecular dynamics simulations to probe the energetics and dynamics of the NADH-driven PCET reaction in complex I. We find that the reaction takes place by concerted hydrogen atom (H•) transfer that couples to an electron transfer (eT) between the aromatic ring systems of the cofactors and further triggers reduction of the nearby FeS centers. In bacterial, Escherichia coli-like complex I isoforms, reduction of the N1a FeS center increases the binding affinity of the oxidized NAD+ that prevents the nucleotide from leaving prematurely. This electrostatic trapping could provide a protective gating mechanism against reactive oxygen species formation. We also find that proton transfer from the transient FMNH• to a nearby conserved glutamate (Glu97) residue favors eT from N1a onward along the FeS chain and modulates the binding of a new NADH molecule. The PCET in complex I isoforms with low-potential N1a centers is also discussed. On the basis of our combined results, we propose a putative mechanistic model for the NADH-driven proton/electron-transfer reaction in complex I.

The respiratory complex I is a redox-driven proton pump that employs the free energy released from quinone reduction to pump protons across its complete ca. 200 Å wide membrane domain. Despite recently resolved structures and molecular simulations, the exact mechanism for the proton transport process remains unclear. Here we combine large-scale molecular simulations with quantum chemical density functional theory (DFT) models to study how contacts between neighboring antiporter-like subunits in the membrane domain of complex I affect the proton transfer energetics. Our combined results suggest that opening of conserved Lys/Glu ion pairs within each antiporter-like subunit modulates the barrier for the lateral proton transfer reactions. Our work provides a mechanistic suggestion for key coupling effects in the long-range force propagation process of complex I.