A modeling and simulation perspective on the mechanism and function of respiratory complex I

The mechanisms of action of mitochondrial targeting agents in cancer: inhibiting oxidative phosphorylation and inducing apoptosis

The tumor microenvironment affects the structure and metabolic function of mitochondria in tumor cells. This process involves changes in metabolic activity, an increase in the amount of reactive oxygen species (ROS) in tumor cells compared to normal cells, the production of more intracellular free radicals, and the activation of oxidative pathways. From a practical perspective, it is advantageous to develop drugs that target mitochondria for the treatment of malignant tumors. Such drugs can enhance the selectivity of treatments for specific cell groups, minimize toxic effects on normal tissues, and improve combinational treatments. Mitochondrial targeting agents typically rely on small molecule medications (such as synthetic small molecules agents, active ingredients of plants, mitochondrial inhibitors or autophagy inhibitors, and others), modified mitochondrial delivery system agents (such as lipophilic cation modification or combining other molecules to form targeted mitochondrial agents), and a few mitochondrial complex inhibitors. This article will review these compounds in three main areas: oxidative phosphorylation (OXPHOS), changes in ROS levels, and endogenous oxidative and apoptotic processes.

Horizontal proton transfer across the antiporter-like subunits in mitochondrial respiratory complex I

Respiratory complex I is a redox-driven proton pump contributing to about 40% of total proton motive force required for mitochondrial ATP generation. Recent high-resolution cryo-EM structural data revealed the positions of several water molecules in the membrane domain of the large enzyme complex. However, it remains unclear how protons flow in the membrane-bound antiporter-like subunits of complex I. Here, we performed multiscale computer simulations on high-resolution structural data to model explicit proton transfer processes in the ND2 subunit of complex I. Our results show protons can travel the entire width of antiporter-like subunits, including at the subunit-subunit interface, parallel to the membrane. We identify a previously unrecognized role of conserved tyrosine residues in catalyzing horizontal proton transfer, and that long-range electrostatic effects assist in reducing energetic barriers of proton transfer dynamics. Results from our simulations warrant a revision in several prevailing proton pumping models of respiratory complex I.

From the 'black box' to 'domino effect' mechanism: what have we learned from the structures of respiratory complex I

My group and myself have studied respiratory complex I for almost 30 years, starting in 1994 when it was known as a L-shaped giant 'black box' of bioenergetics. First breakthrough was the X-ray structure of the peripheral arm, followed by structures of the membrane arm and finally the entire complex from Thermus thermophilus. The developments in cryo-EM technology allowed us to solve the first complete structure of the twice larger, ∼1 MDa mammalian enzyme in 2016. However, the mechanism coupling, over large distances, the transfer of two electrons to pumping of four protons across the membrane remained an enigma. Recently we have solved high-resolution structures of mammalian and bacterial complex I under a range of redox conditions, including catalytic turnover. This allowed us to propose a robust and universal mechanism for complex I and related protein families. Redox reactions initially drive conformational changes around the quinone cavity and a long-distance transfer of substrate protons. These set up a stage for a series of electrostatically driven proton transfers along the membrane arm ('domino effect'), eventually resulting in proton expulsion from the distal antiporter-like subunit. The mechanism radically differs from previous suggestions, however, it naturally explains all the unusual structural features of complex I. In this review I discuss the state of knowledge on complex I, including the current most controversial issues.

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.

Docking and molecular simulations reveal a quinone binding site on the surface of respiratory complex I

The first component of the mitochondrial electron transport chain is respiratory complex I. Several high-resolution structures of complex I from different species have been resolved. However, despite these significant achievements, the mechanism of redox-coupled proton pumping remains elusive. Here, we combined atomistic docking, molecular dynamics simulations and site-directed mutagenesis on respiratory complex I from Yarrowia lipolytica to identify a quinone (Q) binding site on its surface near the horizontal amphipathic helices of ND1 and NDUFS7 subunits. The surface-bound Q makes stable interactions with conserved charged and polar residues, including the highly conserved Arg72 from the NDUFS7 subunit. The binding and dynamics of a Q molecule at the surface-binding site raises interesting possibilities about the mechanism of complex I, which are discussed.

Substrate profiling of mitochondrial caseinolytic protease P via a site-specific photocrosslinking approach

Approaches for profiling protease substrates are critical for defining protease functions, but remain challenging tasks. We combine genetic code expansion, photocrosslinking and proteomics to identify substrates of the mitochondrial (mt) human caseinolytic protease P (hClpP). Site‐specific incorporation of the diazirine‐bearing amino acid DiazK into the inner proteolytic chamber of hClpP, followed by UV‐irradiation of cells, allows to covalently trap substrate proteins of hClpP and to substantiate hClpP's major involvement in maintaining overall mt homeostasis. In addition to confirming many of the previously annotated hClpP substrates, our approach adds a diverse set of new proteins to the hClpP interactome. Importantly, our workflow allows identifying substrate dynamics upon application of external cues in an unbiased manner. Identification of unique hClpP‐substrate proteins upon induction of mt oxidative stress, suggests that hClpP counteracts oxidative stress by processing of proteins that are involved in respiratory chain complex synthesis and maturation as well as in catabolic pathways.

High-resolution structure and dynamics of mitochondrial complex I-Insights into the proton pumping mechanism

Mitochondrial NADH:ubiquinone oxidoreductase (complex I) is a 1-MDa membrane protein complex with a central role in energy metabolism. Redox-driven proton translocation by complex I contributes substantially to the proton motive force that drives ATP synthase. Several structures of complex I from bacteria and mitochondria have been determined, but its catalytic mechanism has remained controversial. We here present the cryo-EM structure of complex I from Yarrowia lipolytica at 2.1-Å resolution, which reveals the positions of more than 1600 protein-bound water molecules, of which ~100 are located in putative proton translocation pathways. Another structure of the same complex under steady-state activity conditions at 3.4-Å resolution indicates conformational transitions that we associate with proton injection into the central hydrophilic axis. By combining high-resolution structural data with site-directed mutagenesis and large-scale molecular dynamic simulations, we define details of the proton translocation pathways and offer insights into the redox-coupled proton pumping mechanism of 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.

Iron-sulfur flavoenzymes: the added value of making the most ancient redox cofactors and the versatile flavins work together

Iron-sulfur (Fe-S) flavoproteins form a broad and growing class of complex, multi-domain and often multi-subunit proteins coupling the most ancient cofactors (the Fe-S clusters) and the most versatile coenzymes (the flavin coenzymes, FMN and FAD). These enzymes catalyse oxidoreduction reactions usually acting as switches between donors of electron pairs and acceptors of single electrons, and vice versa. Through selected examples, the enzymes' structure−function relationships with respect to rate and directionality of the electron transfer steps, the role of the apoprotein and its dynamics in modulating the electron transfer process will be discussed.

Essential role of accessory subunit LYRM6 in the mechanism of mitochondrial complex I

Respiratory complex I catalyzes electron transfer from NADH to ubiquinone (Q) coupled to vectorial proton translocation across the inner mitochondrial membrane. Despite recent progress in structure determination of this very large membrane protein complex, the coupling mechanism is a matter of ongoing debate and the function of accessory subunits surrounding the canonical core subunits is essentially unknown. Concerted rearrangements within a cluster of conserved loops of central subunits NDUFS2 (β1-β2^S2 loop), ND1 (TMH5-6^ND1 loop) and ND3 (TMH1-2^ND3 loop) were suggested to be critical for its proton pumping mechanism. Here, we show that stabilization of the TMH1-2^ND3 loop by accessory subunit LYRM6 (NDUFA6) is pivotal for energy conversion by mitochondrial complex I. We determined the high-resolution structure of inactive mutant F89A^LYRM6 of eukaryotic complex I from the yeast Yarrowia lipolytica and found long-range structural changes affecting the entire loop cluster. In atomistic molecular dynamics simulations of the mutant, we observed conformational transitions in the loop cluster that disrupted a putative pathway for delivery of substrate protons required in Q redox chemistry. Our results elucidate in detail the essential role of accessory subunit LYRM6 for the function of eukaryotic complex I and offer clues on its redox-linked proton pumping mechanism.

Respiratory complex I plays a key role in energy metabolism. Cryo-EM structure of a mutant accessory subunit LYRM6 from the yeast Yarrowia lipolytica and molecular dynamics simulations reveal conformational changes at the interface between LYRM6 and subunit ND3, propagated further into the complex. These findings offer insight into the mechanism of proton pumping by respiratory complex I.

Respiratory complex I: Bottleneck at the entrance of quinone site requires conformational change for its opening

The structure of the entire respiratory complex I is now known at reasonably high resolution for many species - bacteria, yeast, and several mammals, including human. The structure reveals an almost 30 angstrom tunnel-like chamber for ubiquinone binding in the core part of the enzyme, at the joint between the membrane and hydrophilic arms of the enzyme. Here we characterize the geometric bottleneck forming the entrance of the quinone reaction chamber. Computer simulations of quinone/quinol passage through the bottleneck suggest that in all structures available, from bacterial to human, this bottleneck is too narrow for the quinone or quinol to pass and that a conformational change is required to open the channel. Moreover, the bottleneck is too narrow even for isoprenoid tail free passage. The closed structure can be an artifact of the crystallization packing forces, low temperature, or other unnatural conditions occurring in the structural data acquisition procedure that affect this flexible part of the enzyme. Two of the helices forming the bottleneck are in direct contact with the subunit (ND3) that was recently demonstrated to be involved in conformational changes during the redox proton pumping cycle, which indicates flexibility of that part of the enzyme. We conclude that the published structures are all locked in the unfunctional states and do not represent correctly the functional enzyme; we discuss possible ways to open the structure in the context of possible mechanisms of the enzyme.

The coupling mechanism of mammalian respiratory complex I

Mitochondrial complex I couples NADH:ubiquinone oxidoreduction to proton pumping by an unknown mechanism. Here, we present cryo-electron microscopy structures of ovine complex I in five different conditions, including turnover, at resolutions up to 2.3 to 2.5 angstroms. Resolved water molecules allowed us to experimentally define the proton translocation pathways. Quinone binds at three positions along the quinone cavity, as does the inhibitor rotenone that also binds within subunit ND4. Dramatic conformational changes around the quinone cavity couple the redox reaction to proton translocation during open-to-closed state transitions of the enzyme. In the induced deactive state, the open conformation is arrested by the ND6 subunit. We propose a detailed molecular coupling mechanism of complex I, which is an unexpected combination of conformational changes and electrostatic interactions.

Key role of quinone in the mechanism of respiratory complex I

Complex I is the first and the largest enzyme of respiratory chains in bacteria and mitochondria. The mechanism which couples spatially separated transfer of electrons to proton translocation in complex I is not known. Here we report five crystal structures of T. thermophilus enzyme in complex with NADH or quinone-like compounds. We also determined cryo-EM structures of major and minor native states of the complex, differing in the position of the peripheral arm. Crystal structures show that binding of quinone-like compounds (but not of NADH) leads to a related global conformational change, accompanied by local re-arrangements propagating from the quinone site to the nearest proton channel. Normal mode and molecular dynamics analyses indicate that these are likely to represent the first steps in the proton translocation mechanism. Our results suggest that quinone binding and chemistry play a key role in the coupling mechanism of complex I.

Complex I (NADH:ubiquinone oxidoreductase) is the first enzyme of the respiratory chain in bacteria and mitochondria. Here, the authors present cryo-EM and crystal structures of T. thermophilus complex I in different conformational states and further analyse them by Normal Mode Analysis and molecular dynamics simulations and conclude that quinone redox reactions are important for the coupling mechanism of complex I.

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.

A novel m.11406 T > A mutation in mitochondrial ND4 gene causes MELAS syndrome

Pathogenic point mutations of mitochondrial DNA (mtDNA) are associated with a large number of heterogeneous diseases involving multiple systems with which patients may present with a wide range of clinical phenotypes. In this study, we describe a novel heteroplasmic missense mutation, m.11406 T > A, of the ND4 gene encoding the subunit 4 of mitochondrial complex I in a 32-year-old woman with recurrent epileptic seizure, headache and bilateral hearing loss. Skeletal muscle histochemistry demonstrated that approximately 20% of fibers were cytochrome C oxidase (COX) deficient with increased activity of succinate dehydrogenase (SDH). Further investigations in muscle specimens showed significantly reduced level of ND4 protein. It is interesting that the subunits of complex I (ND1 and NDFUB8) and complex IV(CO1) were also remarkably decreased. These findings indicate that ND1, NDFUB8 and CO1 are more susceptible than other subunits to mutations in the mitochondrial ND4 gene.

Respiratory complex I - structure, mechanism and evolution

Respiratory complex I is an intricate multi-subunit membrane protein with a central function in aerobic energy metabolism. During the last years, structures of mitochondrial complex I and respiratory supercomplexes were determined by cryo-EM at increasing resolution. Structural and computational studies have shed light on the dynamics of proton translocation pathways, the interaction of complex I with lipids and the unusual access pathway of ubiquinone to the active site. Recent advances in understanding complex I function include characterization of specific conformational changes that are critical for proton pumping. Cryo-EM structures of the NADH dehydrogenase-like (NDH) complex of photosynthesis and a bacterial membrane bound hydrogenase (MBH) have provided a broader perspective on the complex I superfamily.

Respiratory complex I - Mechanistic insights and advances in structure determination

Complex I is the largest and most intricate redox-driven proton pump of the respiratory chain. The structure of bacterial and mitochondrial complex I has been determined by X-ray crystallography and cryo-EM at increasing resolution. The recent cryo-EM structures of the complex I-like NDH complex and membrane bound hydrogenase open a new and more comprehensive perspective on the complex I superfamily. Functional studies and molecular modeling approaches have greatly advanced our understanding of the catalytic cycle of complex I. However, the molecular mechanism by which energy is extracted from the redox reaction and utilized to drive proton translocation is unresolved and a matter of ongoing debate. Here, we review progress in structure determination and functional characterization of complex I and discuss current mechanistic models.

Kinetics and Efficiency of Energy-Transducing Enzymes

Complexes I to IV, with the exception of Complex II, are redox-driven proton pumps that convert redox energy of oxygen reduction to proton gradient across the mitochondrial or bacterial membrane; in turn, the created electrochemical gradient drives the adenosine triphosphate synthesis in the cells by utilizing complex V of the chain. Here we address a general question of the efficiency of such enzymes, considering them as molecular machines that couple endergonic and exergonic reactions and converting one form of free energy into another. One well-known example of the efficiency is given by Carnot's theorem for heat engines. Here we extend the concept to respiratory enzymes and specifically focus on the proton pumping by Complex I of the respiratory chain, nicotinamide adenine dinucleotide dehydrogenase. To discuss the efficiency issues, we develop a model of enzyme kinetics, which generalizes the Michaelis-Menten model. Our model includes several substrates and products and, in general, can be considered as Generalized Michaelis-Menten Kinetic model. The model might be useful for describing complex enzyme kinetics, regardless of the efficiency issues that are addressed in this paper.

Mutations in a conserved loop in the PSST subunit of respiratory complex I affect ubiquinone binding and dynamics

Respiratory complex I catalyses the reduction of ubiquinone (Q) from NADH coupled to proton pumping across the inner membrane of mitochondria. The electrical charging of the inner mitochondrial membrane drives the synthesis of ATP, which is used to power biochemical reactions of the cell. The recent surge in structural data on complex I from bacteria and mitochondria have contributed to significant understanding of its molecular architecture. However, despite these accomplishments, the role of various subdomains in redox-coupled proton pumping remains entirely unclear. In this work, we have mutated conserved residues in the loop of the PSST subunit that faces the ~30 Å long unique Q-binding tunnel of respiratory complex I. The data show a drastic decrease in Q reductase activity upon mutating several residues despite full assembly of the complex. In-silico modeling and multiple microsecond long molecular dynamics simulations of wild-type and enzyme variants with exchanges of conserved arginine residues revealed remarkable ejection of the bound Q from the site near terminal electron donor N2. Based on experiments and long-time scale molecular simulations, we identify microscopic elements that dynamically control the diffusion of Q and are central to redox-coupled proton pumping in respiratory complex I.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Role of Second Quinone Binding Site in Proton Pumping by Respiratory Complex I

Respiratory complex I performs the reduction of quinone (Q) to quinol (QH₂) and pumps protons across the membrane. Structural data on complex I have provided spectacular insights into the electron and proton transfer paths, as well as into the long (~30 Å) and unique substrate binding channel. However, due to missing structural information on Q binding modes, it remains unclear how Q reduction drives long range (~20 nm) redox-coupled proton pumping in complex I. Here we applied multiscale computational approaches to study the dynamics and redox chemistry of Q and QH₂. Based on tens of microseconds of atomistic molecular dynamics (MD) simulations of bacterial and mitochondrial complex I, we find that the dynamics of Q is remarkably rapid and it diffuses from the N2 binding site to another stable site near the entrance of the Q channel in microseconds. Analysis of simulation trajectories also reveal the presence of yet another Q binding site 25–30 Å from the N2 center, which is in remarkable agreement with the electron density observed in recent cryo electron microscopy structure of complex I from Yarrowia lipolytica. Quantum chemical computations on the two Q binding sites closer to the entrance of the Q tunnel reveal redox-coupled protonation reactions that may be important in driving the proton pump of complex I.

Exploring the quinone/inhibitor-binding pocket in mitochondrial respiratory complex I by chemical biology approaches

NADH-quinone oxidoreductase (respiratory complex I) couples NADH-to-quinone electron transfer to the translocation of protons across the membrane. Even though the architecture of the quinone-access channel in the enzyme has been modeled by X-ray crystallography and cryo-EM, conflicting findings raise the question whether the models fully reflect physiologically relevant states present throughout the catalytic cycle. To gain further insights into the structural features of the binding pocket for quinone/inhibitor, we performed chemical biology experiments using bovine heart sub-mitochondrial particles. We synthesized ubiquinones that are oversized (SF-UQs) or lipid-like (PC-UQs) and are highly unlikely to enter and transit the predicted narrow channel. We found that SF-UQs and PC-UQs can be catalytically reduced by complex I, albeit only at moderate or low rates. Moreover, quinone-site inhibitors completely blocked the catalytic reduction and the membrane potential formation coupled to this reduction. Photoaffinity-labeling experiments revealed that amiloride-type inhibitors bind to the interfacial domain of multiple core subunits (49 kDa, ND1, and PSST) and the 39-kDa supernumerary subunit, although the latter does not make up the channel cavity in the current models. The binding of amilorides to the multiple target subunits was remarkably suppressed by other quinone-site inhibitors and SF-UQs. Taken together, the present results are difficult to reconcile with the current channel models. On the basis of comprehensive interpretations of the present results and of previous findings, we discuss the physiological relevance of these models.

Redox-coupled quinone dynamics in the respiratory complex I

Complex I couples the free energy released from quinone (Q) reduction to pump protons across the biological membrane in the respiratory chains of mitochondria and many bacteria. The Q reduction site is separated by a large distance from the proton-pumping membrane domain. To address the molecular mechanism of this long-range proton-electron coupling, we perform here full atomistic molecular dynamics simulations, free energy calculations, and continuum electrostatics calculations on complex I from Thermus thermophilus We show that the dynamics of Q is redox-state-dependent, and that quinol, QH₂, moves out of its reduction site and into a site in the Q tunnel that is occupied by a Q analog in a crystal structure of Yarrowia lipolytica We also identify a second Q-binding site near the opening of the Q tunnel in the membrane domain, where the Q headgroup forms strong interactions with a cluster of aromatic and charged residues, while the Q tail resides in the lipid membrane. We estimate the effective diffusion coefficient of Q in the tunnel, and in turn the characteristic time for Q to reach the active site and for QH₂ to escape to the membrane. Our simulations show that Q moves along the Q tunnel in a redox-state-dependent manner, with distinct binding sites formed by conserved residue clusters. The motion of Q to these binding sites is proposed to be coupled to the proton-pumping machinery in complex I.