High-resolution in situ structures of mammalian respiratory supercomplexes

Mitochondria play a pivotal part in ATP energy production through oxidative phosphorylation, which occurs within the inner membrane through a series of respiratory complexes1–4. Despite extensive in vitro structural studies, determining the atomic details of their molecular mechanisms in physiological states remains a major challenge, primarily because of loss of the native environment during purification. Here we directly image porcine mitochondria using an in situ cryo-electron microscopy approach. This enables us to determine the structures of various high-order assemblies of respiratory supercomplexes in their native states. We identify four main supercomplex organizations: I1III2IV1, I1III2IV2, I2III2IV2 and I2III4IV2, which potentially expand into higher-order arrays on the inner membranes. These diverse supercomplexes are largely formed by ‘protein–lipids–protein’ interactions, which in turn have a substantial impact on the local geometry of the surrounding membranes. Our in situ structures also capture numerous reactive intermediates within these respiratory supercomplexes, shedding light on the dynamic processes of the ubiquinone/ubiquinol exchange mechanism in complex I and the Q-cycle in complex III. Structural comparison of supercomplexes from mitochondria treated under different conditions indicates a possible correlation between conformational states of complexes I and III, probably in response to environmental changes. By preserving the native membrane environment, our approach enables structural studies of mitochondrial respiratory supercomplexes in reaction at high resolution across multiple scales, from atomic-level details to the broader subcellular context.

states in respiratory complexes without compromising the resolution of the cryo-EM maps.Overall, our in situ approach enables investigation of the impacts of diverse mitochondrial diseases and pharmacological treatments by determining reactive protein structures under physiological conditions within mitochondria.

Four main SC types in situ
Three-dimensional (3D) reconstruction unveiled diverse forms of SCs with distinct compositions, including four dominant types: I 1 III 2 IV 1 (type A), I 1 III 2 IV 2 (type B), I 2 III 2 IV 2 (type O) and I 2 III 4 IV 2 (type X) (Fig. 2a, Extended Data Fig. 2 and Supplementary Figs.2-5).Types A and O are similar to previously reported in vitro structures [17][18][19][20] , whereas types B and X represent two new forms not observed in vitro.Our cross-classification results indicate that the type-A SC (Fig. 1h-j) is the most abundant form, determined at resolutions of 1.8-2.4Å in most CI and CIII 2 regions (Supplementary Figs. 2 and 6) and at approximately 2.75 Å in CIV (Supplementary Fig. 2).The high-quality maps, as demonstrated by discernible 'holes' in the side chains of aromatic residues, the precise positions of pyrrole rings and the Fe atom in haem, and S/Fe atoms in SF 4 , enabled us to build accurate atomic models of this form (Fig. 1k, Supplementary Fig. 6 and Supplementary Video 2).For the type-A SC, we modelled 6,327 water molecules into the structure (Extended Data Fig. 3), including those that were likely to have central roles in proton transfer across the membrane.Furthermore, we built 197 structured and associated lipids (Fig. 2c, Extended Data Fig. 4 and Supplementary Video 3) that contribute to stabilization of the protein structure, enhancing the stability of the SCs by facilitating essential protein-lipid-protein interactions (Fig. 2c,d), creating a hydrophobic environment in the Q-binding sites (Extended Data Fig. 4b and Supplementary Video 3) and participating in the hydrogen-bond network through their polar heads (Extended Data Fig. 5 and Supplementary Video 4).We directly visualized mitochondrial inner membranes composed of more dynamic lipids surrounding the proteins and built atomic models (Fig. 2c, Extended Data Fig. 4 and Supplementary Video 3).Although the overall architecture resembles a previously reported structure using in vitro purified protein 17 , our high-resolution in situ structure shows substantial differences at the interaction interfaces among complexes I, III 2 and IV (Supplementary Fig. 7).Previous work has demonstrated that the lipid bilayer can bend to adapt to CI conformations 22 .Consistent with these findings, the membrane enveloping the SC exhibited noticeable curvature, varying across different regions (Fig. 2b).Specifically, the membrane at the CI heel bends towards the mitochondrial matrix, whereas the CIII 2 region displays the opposite curvature pattern (Fig. 2b, Extended Data Fig. 6 and Supplementary Video 5).
Our high-resolution density maps of native respiratory SCs allow for accurate analysis of the interaction interfaces.In type-A SC (Fig. 2a), CI-CIII 2 interactions are mediated through specific contacts between the NDUFA11 I and UQCRQ III subunits, as well as NDUFB9 I , NDUFB4 I and UQCRC1 III .For CI-CIV interactions, the key participating subunits are NDUFB3 I , ND5 I and COX7A IV (Supplementary Fig. 7).One striking feature of the in situ structure is that the interstitial space among CI, CIII 2 and CIV is populated by lipids that mediate the complex interactions.

Article
In particular, no direct protein-protein interactions between CIII 2 and CIV were observed in our structures (Fig. 2d).
In contrast to type A, the type-B SC (Fig. 2a) incorporates a second CIV (CIV′) situated between CIII 2 and the hydrophilic arms of CI.Unlike CIII 2 , which exhibits two-fold symmetry, this extra CIV in the type-B SC is not two-fold symmetric to the first CIV relative to the symmetry axis of CIII 2 .Instead, it displays an approximately 60° rotation (Supplementary Fig. 8a,b), partially enclosing the lipid bilayer between the Q-binding pockets in CI and CIII 2 (Supplementary Fig. 9 and Supplementary Video 6) and forming an architecture that could confine free diffusion of Q 10 and facilitate Q-channelling.In type-B SC, CIV′ establishes new interactions with CI through the NDUFA1 I and COX5A IV′ subunits and with CIII 2 through the UQCR10 III and COX6A IV′ subunits (Supplementary Fig. 7).
In addition, two higher-order SCs were determined through a multilevel cross-classification approach we developed.The I 2 III 2 IV 2 (Fig. 2a) complex at 2.6-3.3Å resolution (Supplementary Fig. 4), designated type O owing to its overall shape, resembles the structure of the human mitochondrial megacomplex 20 (Supplementary Fig. 8e) and can be considered to be a pseudo-C 2 symmetric expansion of type A, sharing the CIII 2 dimer (Supplementary Fig. 8c).In contrast to the detergent-purified type-O SC, we observed that the native type-O SC exhibits pseudo-C 2 rather than strict C 2 symmetry, with the surrounding membrane bending towards the mitochondrial matrix.This observation is not surprising given the highly curved membrane structures of the cristae.
The new SC I 2 III 4 IV 2 (Fig. 2a), termed type X owing to its chromosomelike appearance, is formed by two type-A SCs aligned in a head-to-head manner along a pseudo-two-fold axis (Supplementary Fig. 8d).
No canonical strong protein-protein interactions were evident at the dimerization interfaces between the two type-A SCs, except for potential weak contacts between helices NDUFB3 I and COX4 IV .The dimerization is primarily mediated by 'protein-lipids-protein' interactions in our structure, as evidenced by the numerous lipid molecules filling the interfacial spaces (Supplementary Video 3).
Notably, in all the four main SC forms, CIV shows the most freedom in terms of its binding positions and modes of interaction with other complexes.Different isoforms of CIV subunits have been reported and may affect SC formation.However, using the porcine heart mitochondria, we only identified one specific set of CIV subunit isoforms (Supplementary Table 7), probably because their expression levels are highest in this tissue 23 .Moreover, our structural analyses consistently identified the NDUFA4 subunit as an integral part of CIV across the four types of SC, distinctly separated from CI (Supplementary Fig. 10); this clarifies the positional association of NDUFA4 and is consistent with the findings of a recent study 24 .
Furthermore, the formation of these SCs substantially influences the local curvature of the surrounding membranes (Fig. 2b, Extended Data Fig. 6 and Supplementary Video 5).Vice versa, we speculate that the membrane geometry in turn affects the overall arrangement, distribution and conformation of these SCs; this is worth further investigation.
In addition, we detected potential extensions of respiratory SCs into high-order arrays in the mitochondrial inner membranes.

Multiple Q-binding states in CI
CI orchestrates electron transfer from NADH to ubiquinone and concurrently translocates protons across the inner mitochondrial membrane 25 .The Q-binding site in CI features an unusually long, heterogeneous channel structure [26][27][28] .The entry, exit and interactions of Q 10 within this channel are central yet not fully resolved questions in our understanding of the molecular mechanisms of CI.Recent structural studies and molecular dynamic simulations have resolved distinct intermediates of Q-binding states 22,29 .
We identified five main Q-binding states (Fig. 3), including the apo state and a previously described Q-bound state 28,29 (state α), along with three new Q-bound states (states β-δ) (Fig. 3b and Extended Data Fig. 7), which are different from the three Q-binding sites previously reported 29 (Extended Data Fig. 8).The apo state displays only noisy densities in the Q channel.State α resembles the active form of bovine and porcine CI (PDB: 7QSK 28 , 7V2C 29 ) (Extended Data Fig. 7b), with a Q 10 molecule fully occupying the Q-binding site (Fig. 3c).In this state, the 1-carbonyl group of Q 10 is situated approximately 20 Å (Fig. 3b) above the membrane surface (M-distance) and 41.1 Å (Fig. 3d) from the channel entrance (E-distance).The three other states show not only conformational differences in the Q 10 headgroup and the key residue H59 NDUFS2 (Extended Data Fig. 7a) but also long-range structural changes away from the Q 10 headgroup and along the whole Q 10 .In state β, the Q headgroup is angled away from the Fe-S cluster N 2 relative to state α, although the tail region largely remains similar to that in state α (Fig. 3d).These conformational changes result in an M-distance of 16 Å (Fig. 3b) and an E-distance of 38.6 Å (Fig. 3d).
State γ is distinguished by housing a Q 10 molecule with a headgroup aligned more parallel with the membrane and a contorted midsection of the Q tail.Although the M-distance remains consistent with that in state β, the E-distance is reduced by 3 Å (Fig. 3d), a change attributable to its silkworm-like undulatory motion (Fig. 3e).In state δ, a Q 10 molecule only partially fills the Q channel (E-distance = 24.5 Å).Its headgroup aligns flush with the membrane surface (M-distance = 0), while the tail spans the entire lipid bilayer (Fig. 3b).Furthermore, the headgroup in state δ is surrounded by a cluster of acidic, basic and polar residues inside the Q channel (Extended Data Fig. 7).These residues seem to create a highly polarized environment connected to the aqueous matrix, facilitating Q protonation.Taken together, our findings support a model in which Q 10 progresses through the Q channel by peristaltic motion (Fig. 3e).

Active-deactive transitions in CI
Under physiological conditions and without substrates, mammalian CI transitions from an active, ready-to-catalyse state to a substantially deactive resting state 28,30,31 .During ischaemia, a condition characterized by limited oxygen supply, the deactive state emerges owing to cessation of the electron transport chain 28 .To investigate the effects of varying levels of ischaemia on the atomic-level structure of CI in native mitochondria, we exposed porcine hearts to room temperature for  Article durations of 0 min ('fresh'), 40 min ('mild') and more than 4 h ('harsh') before mitochondrial isolation (Fig. 1b).The active and deactive states of CI are delineated by distinct structural hallmarks, including domain movements and conformational alterations around the Q-binding site and in proximal membrane-domain subunits 28 .Using focused 3D classification, we analysed these hallmarks in the SC structures determined from mitochondria under the three conditions (Supplementary Figs.11 and 12).In the fresh sample, approximately 75% of CI in SCs adopted an active state, compared with around 30% and 18% under mild and harsh conditions, respectively (Extended Data Fig. 9).
Previous in vitro structural studies have suggested that the hallmarks of CI can be collectively categorized into either the active or deactive form 26 .To further investigate whether there existed intermediate states involved in the transition from active to deactive (A-D transition), we performed further focused 3D classification targeting the regions around the Q site using all datasets combined, resulting in seven main intermediate classes.We built the models in these states and compared them with the canonical active and deactive structures (Supplementary Fig. 12 and Supplementary Table 8).The comparison unveiled two distinct classes (classes 0 and 7 in Supplementary Table 8) that fully conformed to the hallmarks of the active and deactive states, respectively.By contrast, other classes exhibited only a subset of the hallmarks associated with the deactive state, indicating that there may exist various stages in the transition between the two extreme states in the native mitochondrial environment.

Catalytic states of CIII 2 in reaction
CIII 2 transfers electrons from ubiquinol to cytochrome c (cyt c) and contributes to the proton gradient for ATP synthesis 32,33 .Comprising three core subunits (cyt b, cyt c 1 and the Rieske iron-sulfur protein (ISP)) along with eight auxiliary subunits, each CIII hosts four metal centres (haems b H and b L in cyt b, the [2Fe-2S] cluster in ISP and haem c 1 in cyt c 1 ) and two Q sites: Q i for Q reduction and Q o for QH 2 oxidation [34][35][36] .From our high-resolution maps, with local resolutions between 1.8 and 2.4 Å (Supplementary Figs. 2 and 6) except for the dynamic ISP, we unambiguously identified all reactive centres and built atomic models of all endogenous ligands (Fig. 4a,b and Extended Data Fig. 4).Focused 3D classification showed that endogenous Q 10 ligands bound to both Q sites adopt multiple conformations representing different reaction stages.Further well-defined densities in the Q-binding pockets were identified as structured lipids with their phosphate headgroups tightly bound to the surrounding protein regions.Tails of these lipids create a hydrophobic, dynamic environment for Q-binding and release.
During the Q-cycle, Q 10 at the Q i site undergoes a two-step reduction, acquiring two protons from the mitochondrial matrix to form QH 2 (ref.37).Our structure showed a hydrogen-bond network near the Q i site, comprising water molecules (Fig. 4c,d, Extended Data Fig. 5 and Supplementary Video 4), polar amino acid residues and the phosphate heads of three lipids (two cardiolipin and one phosphatidylethanolamine) (Fig. 4d and Extended Data Fig. 5), which form a Grotthuss-competent system for proton transfer 25,38 .The network could be divided into two hydrogen-bond chains primarily constituted of water molecules that link Q 10 to the matrix (Supplementary Video 4).One of the two chains further branches above the three lipids that is directly involved in facilitating the proton transfer through the hydrogen network.This bifurcated chain engages two charged amino acids, K227 MTCYB and D228 MTCYB , with D228 MTCYB linked to the 1-carbonyl of Q 10 (Fig. 4c and Extended Data Fig. 5).The second chain, comprising 5 waters that surround H201 MTCYB (Fig. 4c), is directly connected to 4-carbonyl of Q 10 .These chains provide an explanation for how the rapid proton transfer process is achieved in the Q i site.
At the Q o site, QH 2 releases two protons into the intermembrane space as part of its reaction cycle 39 .Our high-resolution map clearly resolved a prevalent Q-binding conformation at this site, characterized by π-π stacking interactions with F274 MTCYB (Fig. 4e).Similar to the Q i site, there was also a hydrogen-bonded water chain near the Q o site for proton transfer.Notably, residue E271 MTCYB , engaged in this water chain, exhibited clear dual conformations (Conf. 1 and Conf.2) and is likely to play a critical part in proton transfer (Fig. 4e).The two conformations resembled those of two previously reported crystal structures in the apo state 40 (Conf.1) and with the Q o inhibitor stigmatellin bound 41 (Conf.2), respectively (Fig. 4f).Further structural comparison indicated that the headgroup of stigmatellin lies closer to the ISP than the endogenous QH 2 captured in our in situ structure.Focused 3D classification uncovered a less-populated (approximately 13%) Q-binding state (Fig. 4g) at the Q o site with its headgroup shifted around 6 Å closer towards the ISP compared with the dominant form (Fig. 4h,i).We propose that this state, rather than the dominant class, represents a transient electron-hopping phase from QH 2 to [2Fe-2S] in ISP.
The transfer of electrons from QH 2 at the Q o site to cyt c 1 is mediated by movement of the Rieske domain 33,37 .Although early structural studies have hinted at its dynamic role in electron transfer, questions remain about how this motion is initiated, regulated and coupled with Q-binding states during the catalytic cycle.In addition to the apo form, we resolved four distinct Q-bound states (Fig. 4h,i).In state I, the Q 10 headgroup is situated near the [2Fe-2S] cluster at a distance of 10.8 Å, with the Rieske domain exclusively adopting the b position.In states II, III and IV, Q 10 takes on the prevalent conformation revealed in our cryo-EM maps, while the Rieske domain undergoes a stepwise movement from the b position towards the c position, eventually making contact with cyt c 1 for electron transfer in state IV.To further analyse the detailed conformational changes among these states, we measured two characteristic distances, one from [2Fe-2S] to Q 10 and the other to haem c 1 (Fig. 4h).Our results indicate a possible sequential progression from state I to state IV, characterized by increasing [2Fe-2S]-Q 10 distances and decreasing [2Fe-2S]-haem c 1 distances.On the basis of these observations, we propose a model that outlines the coupling between ISP movement and Q-binding states.

CI and CIII 2 conformational correlation
To examine the potential correlations between the conformational states of CI and CIII 2 in the same mitochondrial SC, we conducted focused conformational analyses of paired CI and CIII 2 from samples treated differently (Fig. 1b).In the fresh sample, CI predominantly existed in its active state; concurrently, the ISP domain of CIII 2 adopted the b position (Extended Data Fig. 9b).Under harsh conditions, most CI particles were classified into the deactive form, whereas the ISP domain of CIII 2 preferentially adopted the c position (Extended Data Fig. 9c).Notably, for the mild treatment, despite CI predominantly being in the deactive state, the ISP domain of CIII 2 primarily remained in the b position, indicating that it may maintains catalytic functionality for a short duration after deactivation of CI (Extended Data Fig. 9c).Previous research indicates that the deactive state of CI arises as a result of disrupted electron transfer through the respiratory chain owing to limited oxygen supply under ischaemic conditions 28 .High-resolution structural analyses using X-ray crystallography and single-particle cryo-EM also demonstrate that CIII 2 adopts the c position in the absence of substrate 40 .This behaviour of CI/CIII 2 conformational distributions could be attributed to direct effects of the treatments or to environmental changes induced by different treatments, potentially affecting the turnover of respiratory chain substrates.

Possible functional roles of SCs
The existence of mammalian respiratory SCs has been established by both in situ cryo-ET 42 and in vitro approaches 17,18 , including single-particle cryo-EM and blue-native gel analysis 7,43,44 .However, details of the actual forms and conformations of these SCs within native membranes have remained elusive.The question of whether SCs offer enzymatic or functional advantages continues to be a subject of debate.Studies have not conclusively shown that SCs facilitate enhanced catalytic activity by directly channelling intermediate substrates 45,46 .A recent study suggested that SCs are non-essential for mouse bioenergetics and physiology 44 , yet the authors maintained an open stance on their potential alternative roles, including even distribution of complexes throughout the membrane, regulation of protein stability 47 , minimization of reactive oxygen species production 48 and prevention of age-associated protein aggregation 49 .
Our in situ structures indicate a distance of approximately 100 Å between the Q-entrances of CI and CIII 2 , consistent with previous evidence 17,18 .The space between the two Q-entrances is filled with diffusive lipid molecules and seems to lack a confined space for direct substrate channelling in most of the SC structures we determined; the exception is type B, which seems to possess a semiclosed channel between CI and CIII 2 (Supplementary Fig. 9).

Article
However, this is not the dominant form, indicating that the majority of Q/QH 2 is not necessarily shuttled between the CI and CIII 2 within the same SC.
Our study shows that the local geometry of the slightly curved membrane surrounding various forms of SCs fits the overall shape of the planar regions of mitochondrial cristae, consistent with previous findings that SCs predominantly localize in these regions [50][51][52] .The optimal spatial localization and organization of CI, CIII 2 and CIV can ensure a homogenous distribution of SCs within the inner membrane 53 .Notably, our current cryo-EM classification attempts did not detect noticeable cyt c density, indicating that binding of cyt c may occur instantly, whereas its diffusion within the intermembrane space may be the rate-limiting step.The homogenous arrangement of CI, CIII 2 and CIV helps to minimize the diffusion distance for Q/QH 2 and cyt c, potentially enhancing catalytic efficiency 54 .
Despite a lack of direct Q-channelling structure between CI and CIII 2 , we speculate that the variable local membrane curvature could regulate Q/QH 2 diffusion.Our structures show that the membrane near CI exhibits notable convex curvature towards the matrix, whereas it is distinctly concave near CIII 2 , with a more planar region in between.Studies have demonstrated that membrane curvature can affect the lateral mobility and sorting of lipids, peptides and proteins, thereby influencing the behaviour of small molecules within the membrane 55,56 .Similar rules can be applied to the diffusion and sorting of Q/QH 2 , potentially enriching the required substrates at the Q sites of CI and CIII 2 .
Moreover, conventional approaches such as digitonin solubilization, blue-native gel analysis, cross-linking mass spectrometry and, particularly, mutagenesis on the basis of in vitro structures, may have limitations in their ability to accurately reflect the actual forms and physiological states of SCs in mitochondria.In the future, studies integrating functional and physiological analyses with the in situ imaging approach used in our current work will help to overcome these limitations and enhance our understanding of SC functions.

A-D states and catalytic intermediates
In vitro purified mammalian CI consists of a mixture of the active and deactive resting states 45 .The closed structure, characterized by a reduced angle and tight interactions between the transmembrane and hydrophilic arms, corresponds to the active resting state, whereas the open structure is indicative of the deactive resting state.Consistently, our study also identified two main classes of CI, namely the active and deactive states, characterized by the closed and open conformations.At present, there are two opposing opinions about opening and closing of the two arms in the A-D transition: (1) that the open state is not linked to the catalytic cycle 45 ; and (2) that the catalysis involves a succession of open and closed states, with the enzyme settling in these states upon cessation of catalysis 12 .
In the fresh-condition batch, the active state predominated, with an increase in deactive states upon incubation at room temperature, indicating that the deactive state may arise with gradual depletion of substrates.In addition, active-apo and all different Q/QH 2 -binding states were dominantly in the closed conformation (except state β) (Supplementary Fig. 11), implying that the closed form may be sufficient for the catalytic cycle.The global conformational changes between the A-D forms within the native mitochondrial membrane were considerably smaller than those observed in detergent-purified 12 and membrane scaffold protein nanodisc reconstituted respiratory CI 28 , indicating that the energetic cost required for the conformational changes in the native lipid bilayer environment may be higher than those under in vitro conditions (Extended Data Fig. 10).Notably, our extensive 3D classification of in situ images did not yield a class that resembled the reported relocation of the helix ND6-THM4 (Supplementary Fig. 13) 12 .We speculate that such a relocation may be energetically unfavourable, as this would require substantial reorganization of surrounding lipids associated with CI.

In situ imaging beyond mitochondria
Our study used in situ cryo-EM imaging techniques to directly visualize SC structures within mitochondria.This approach bypasses the limitations associated with traditional in vitro purification methods, which could result in loss of native ligands and physiological states, leading to compositional and conformational artefacts.Our approach preserves both the native membrane environment and the electrochemical proton gradient.This enabled concurrent analysis of SC structures in reaction at both atomic detail and large-scale organization, achieving sub-2 Å local resolution.Optimized in silico classification of the in situ cryo-EM data enabled us to capture a series of dynamic ligand-binding states in the respiratory chain, eliminating the requirement for analogues or reaction inhibitors commonly used in conventional methods and thus providing new insight into the mechanisms of mitochondrial SCs that reflect physiological behaviour.Notably, the preservation of SC structures in native membranes enabled direct visualization of numerous surrounding lipids, which have pivotal roles at the interaction interfaces among the complexes as well as within the Q-binding pockets.These findings have potential extended applications in studies to elucidate the effects of many analogues, drugs and inhibitors on mammalian mitochondria, shedding light on SC behaviours under various physiological and pathological conditions, including cardiopathy, diabetes, neurodegeneration and cancer.Furthermore, our in situ approach could be extended to structural studies of native membrane protein complexes within various organelles, such as cilia, chloroplasts, Golgi and lysosomes.

Online content
Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41586-024-07488-9.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Preparation of porcine mitochondria and cryo-EM grids
Mitochondria were isolated from porcine hearts following a modified version of a protocol originally described by A. L. Smith 57 .Before mitochondrial extraction, the pig hearts were subjected to three distinct treatment conditions: (1) fresh-immediately placed on ice for all subsequent procedures; (2) mild-incubated at room temperature for 40 min and put on ice to quickly cool down for isolation; and (3) harsh-incubated at room temperature for more than 4 h before being cooled on ice.The isolated mitochondria were then resuspended in a solution containing 0.25 M sucrose, 10 mM Tris-buffered with H 2 SO 4 and 0.2 mM EDTA at pH 7.8.The suspension was adjusted to achieve a final optical density at 600 nm of 1.3 absorbance units.
For cryo-EM grid preparation, 3.3 μl of the mitochondrial suspension was applied to each Quantifoil holey carbon grid (R2/1, 300 mesh gold).Grids were incubated for 5 s in a Vitrobot Mark IV (Thermo Fisher Scientific) chamber maintained at 8 °C and 95% relative humidity.Excess solution was blotted using standard Vitrobot filter paper before the grids were rapidly plunged into liquid ethane at a temperature of approximately −170 °C.

Cryo-ET data collection
Grids were initially screened for optimal ice conditions using a 200 kV Glacios microscope (Thermo Fisher Scientific) at the Yale Science Hill Electron Microscopy Facility.Selected grids were subsequently transferred to a 300 kV Titan Krios microscope (Thermo Fisher Scientific), equipped with a Bioquantum Energy Filter and a K3 direct electron detector (Gatan), for high-resolution data acquisition at the Yale West Campus Electron Microscopy Facility.Automated data collection was facilitated using SerialEM software 58 and Gatan DigitalMicrograph.All images were captured in superresolution mode, with a physical pixel size of 6.1 Å (effectively 3.05 Å in superresolution).A total of eight tilt series were collected, targeting a relatively high defocus range, from −6 μm to −10 μm, for better contrast to guarantee a more reliable initial reconstruction.A grouped dose-symmetric scheme, spanning from −60° to 60° at 2° increments, was used for tilt series acquisition, with an accumulated dose of 100 e − /Å 2 .

Cryo-ET reconstruction and subtomogram averaging
Tomogram reconstruction was streamlined using custom scripts.Initial frame alignment was performed using MotionCorr2 (ref.59), followed by micrograph binning at a factor of two.Tilt series stacks were generated using in-house scripts.All tilt series were aligned and reconstructed using AreTomo 1.2. 5 (ref.60).Initial contrast transfer function (CTF) parameters were estimated with GCTF 61 and cryoSPARC 62 .Raw micrographs and reconstructed results were visualized and diagnosed using IMOD 63 and ChimeraX 64 .

Single-particle cryo-EM data collection
Automated data acquisition was performed using either a Glacios or a Titan Krios electron microscope (Thermo Fisher Scientific).The Glacios was equipped with a K3 direct electron detector (Gatan) and operated at 200 kV at a pixel size of 0.434 Å in superresolution mode, with an objective aperture of 100 μm.The Titan Krios, also equipped with a K3 direct electron detector, was operated at 300 kV at a pixel size of 0.416 Å in superresolution mode with a Gatan energy filter.Automatic data collection was facilitated using the SerialEM software package 58 .Multishot acquisition parameters were set at 3 × 3 holes per imaging location, with four exposures per hole at 200 kV and five exposures per hole at 300 kV.The total electron dose was fractionated to 42 e − /Å 2 for the Glacios and 50 e − /Å 2 for the Titan Krios, distributed across 45 frames at 40 ms per frame.Defocus parameters ranged from −1.0 μm to −3.0 μm for the 200 kV dataset and from −1.3 μm to −3.0 μm for the 300 kV datasets.Details of the data collection are summarized in Supplementary Tables 1-6.

Preprocessing
For all datasets, motion correction was performed using MotionCor2 (ref.59) or cryoSPARC 62 .The CTF of each motion-corrected micrograph was estimated using Gctf 61 or cryoSPARC 62 .Particles were picked with Gautomatch or cryoSPARC using an iterative sorting strategy as described below.Cryo-EM scripts used for real-time data transfer and on-the-fly preprocessing can be downloaded from https://github.com/JackZhang-Lab.

Overall particle selection and sorting strategy
Owing to the challenges posed by low signal-to-noise ratios and a highly congested macromolecular environment (Extended Data Fig. 1a), traditional particle selection methodologies were insufficient for generating datasets amenable to reliable two-dimensional (2D) classification, ab initio three-dimensional (3D) reconstruction and subsequent local refinement.To address this issue, we implemented an iterative strategy to optimize particle selection and sorting.The approach involved several rounds of iterative 2D particle picking, 2D classification and 3D analyses including ab initio 3D reconstruction, 3D classification and multilevel local refinement.Unlike the conventional particle selection approach, our strategy used Gautomatch and cryoSPARC 62 for template matching to gradually increase the resolution of 3D projections as the reconstructions were progressively improved over cycles.We used several independent sources of references to cross-validate the final results.To maximize the yield of high-quality particles, particles from the classes that show clear features of SCs in all cycles were merged for subsequent 3D cross-classification.More details of the strategy are explained in the following sections.

Initial 3D reconstruction with surrounding membranes (type A)
Conventional 2D classification failed to generate meaningful class averages using images selected from our in situ cryo-EM micrographs of mitochondria for three main reasons: (1) thick samples that led to low signal-to-noise ratios and large defocus variations, (2) a crowded environment that affected particle detection and alignment, and (3) strong membrane signals that dominated the alignment, leading to blurred averages of protein regions (Extended Data Fig. 1b).
To address this, we initially used the strong membrane signals and focused on the side views surrounded by membranes using 2D classification.These side views in principle contained sufficient orientational information for a complete 3D reconstruction.At the outset, protein signals were completely averaged out in the 2D classification, whereas the membranes were well aligned owing to the strong side-view signals (Extended Data Fig. 1c).We then conducted several cycles of 2D classification to focus only on particles exhibiting clear membrane signals.
Through comprehensive 2D analyses, we found that regions potentially harbouring mitochondrial SCs exhibited special features of local curvature.Specifically, these regions were characterized by membrane signals that seemed to be concave towards the matrix direction, indicative of the presence of CIII 2 (Extended Data Fig. 1c-e).By merging particles from classes with characteristic concave membranes surrounding CIII 2 and conducting further 2D classification, we achieved improved 2D averages showing clear membrane features around CIII 2 (Extended Data Fig. 1d).Notably, extra protein densities adjacent to CIII 2 were obvious, probably representing CI or CIV densities.However, it was unclear how many types of respiratory SC exist in native mitochondria and whether CI, CIII 2 and CIV always appear in the form of SCs or just partially.
To further address these observations and obtain unbiased density maps, we used four independent methods to generate initial references: (1) cryo-ET subvolume averaging (Supplementary Fig. 1), (2) ab initio reconstruction using particles assigned to the 2D averages with visible protein densities (Extended Data Fig. 1f) and characteristic CIII 2 membrane features (Extended Data Fig. 1g), (3) ab initio 3D reconstruction using particles after membrane signal subtraction (Extended Data Fig. 1h), and (4) models generated from random selection of unsorted particles or random noise (Extended Data Fig. 1i).All these references were combined for 3D classification and subsequently used for local refinement and focused classification (Extended Data Fig. 1i).Given that in situ cryo-EM datasets are more heterogeneous than conventional single-particle datasets, we included 'false references' generated from approach (4) for better classification.Finally, particles corresponding to classes showing clear features of type-A SC were re-extracted and merged for further classification and refinement (Extended Data Fig. 1i).

Cross-classification of multiple SCs
Around the reconstructed SC I 1 III 2 IV 1 (type A) map, we observed extra densities, clearly indicating that more proteins bound to the type-A SC to form larger SCs.We suspected that more types of SC existed in native mitochondrial membranes.Preliminary results from both large single-particle 3D classification at low resolutions and cryo-ET subvolume averaging confirmed this speculation.To further improve the accuracy of 3D classification for high-resolution refinement, we deliberately provided extra false references generated from random subdatasets using discarded particles from previous cycles.These false references served to randomly absorb low-quality and falsely picked particles, leading to a relatively clean dataset for the target class.We then accumulated particles classified into good classes, defined by clear secondary structures, over several cycles.Owing to the crowded mitochondrial environment, misclassified and misaligned particles were always present.To address this, we reorganized the particles by merging those that fell into classes generating similar 3D maps.We selected multiple references from different classes, including those considered 'bad' and reperformed 3D classification on each subdataset.Afterwards, we recombined all subsets of different classes that were considered 'good' and reclassified them.On the basis of these results, we then merged all the particles belonging to a specific target from previous cycles and performed a further cycle of 3D classification on the merged dataset.This further classification used high-resolution references generated from previous classification cycles and local refinement to discard low-quality or misclassified particles.
After numerous rounds of cross-classification followed by local refinement, we identified various other types of SC, including the three other main classes: type B (I 1 III 2 IV 2 ), type O (I 2 III 2 IV 2 ) and type X (I 2 III 4 IV 2 ).In addition to the four main classes, other classes such as I 1 III 2 , I 4 III 4 IV 4 and even higher-order assemblies were observed; however, they were not subject to further refinement in this study owing to the low population.Subsequently, we cross-validated our classification results by providing a set of references lacking the correct form of the SC for subclassification of each class.We also performed further reference-free 2D classification after 3D classification and refinement to verify different forms of SC.This allowed us to visualize the distinct features of the four main classes from 2D averages directly, without imposing any references.Only datasets converging to the correct form of supercomplex, regardless of the initial references used, were included in the final multilevel local refinement and focused 3D classification.

Multilevel local refinement and focused 3D classification
A hierarchical masking strategy was used for local refinement on all four main types of SC.Specifically, the mask size was incrementally reduced to focus on distinct regions of each type of respiratory SC, ensuring stable local refinement.We partitioned the type-A SC into five principal domains: (1) CI hydrophilic region, (2) CI hydrophobic region, (3) CIII 2 , (4) CIV and (5) lipid environment.
Before the multilevel local refinement, the type-A SC was refined to 3.39 Å overall using images binned two times (1.664 Å per pixel after binning) with 1,113,902 high-quality particles.This included type B, type O and type X, as they all share the type-A region.We recentred and re-extracted these particles, generating 1,050,463 final particles for subsequent local refinement (particles near the edges were excluded after re-extraction).Initially, the resolutions of CI, CIII 2 and CIV worsened slightly (approximately 3.5 Å) after the first cycle of refinement using the unbinned particles (0.832 Å per pixel).Further improvement was achieved by optimizing several local refinement parameters, including optimization of mask sizes, global CTF, local CTF refinement, local angular refinement and non-uniform refinement 67 .
By iteratively applying these techniques, we refined the maps of the hydrophilic region of CI and the hydrophobic regions of CI, CIII 2 and CIV to average resolutions of 2.46 Å, 2.58 Å, 2.31 Å and 2.66 Å, respectively (Supplementary Fig. 2 and Supplementary Tables 1-6).Even smaller regional masks, focused on CI and CIII 2 , further improved local resolutions.Local resolutions in most of the protein regions of CIII 2 ranged from 1.8 to 2.4 Å (Supplementary Fig. 2c).Focused classification and refinement for specific subdomains, such as the Q/QH 2 binding sites, yielded further improvements that aided in model building.For more complex regions, such as the lipid environment surrounding the transmembrane regions of the SCs and Q/QH 2 binding sites, further levels of focused classification and local refinement were performed.To ensure seamless integration of adjacent regions, all local masks were manually created so that pairs of adjacent masks contained sufficiently large areas for the generation of final composite maps using the smaller regions individually refined.All locally refined segments were integrated into a composite map in ChimeraX 64 .
Similar multilevel refinement approaches were used to determine the structures of other forms of respiratory SC.Detailed parameters and refinement results are summarized in Supplementary Figs.2-5 and Supplementary Tables 1-6.

Membrane signal detection and weakening
One of the critical bottlenecks limiting high-resolution cryo-EM reconstruction of membrane proteins in their native environment is the severe signal interference from surrounding membranes.This interference can significantly affect several steps in cryo-EM data analysis, including ab initio reconstruction, Euler angle determination, and 2D and 3D classification, as well as refinement of alignment parameters.To address this issue, we developed a computational toolkit to detect membrane signals from 2D averages, estimate the local geometry of detected membranes, and suppress or remove these signals to substantially improve the alignment reliability of mitochondrial complexes in native membrane environments.
Initially, we generated a series of 15-30 computationally simulated 2D projections of lipid bilayers, with local curvatures ranging from 0 nm −1 to 0.02 nm −1 .These simulated 2D membranes served as templates for detection of the side-view signals of mitochondrial membranes using Gautomatch.Subsequently, three to five cycles of 2D classification were performed to discard low-quality and non-membrane particles, resulting in a subset of particles showing clear side views of lipid bilayers.We then estimated the approximate orientation and centre of each individual lipid bilayer on the basis of its corresponding 2D average using the Radon transform.Local curvature was determined by maximizing the cross-correlation between each 2D average and a series of simulated lipid bilayers.These curves were rotated and translated using alignment parameters from 2D classification generated by cryoSPARC 62 .Centres of each membrane segment were refined by maximizing the normalized cross-correlation between the raw image and transformed 2D average.Using these estimated parameters, we approximated the principal signals of each membrane segment by locally averaging the image intensities along the membrane curve within a soft mask, which was around 25% larger than the typical lipid bilayer we estimated.Membrane signals that had dominated the alignment in the raw images were weakened to enhance protein signal contributions for subsequent reconstruction, alignment, classification and local refinement.This improved the signal contributions from protein regions for the initial alignment, akin to the critical effects observed in our previously described microtubule signal subtraction method 68,69 .Finally, alignment and classification parameters were applied to the raw images along with membrane signals for subsequent local refinement and focused classification.

Membrane modelling and geometry analysis
The in situ mitochondrial respiratory chain complexes largely preserved the native state of the membrane architecture, as evidenced by exceptionally clear density maps (Extended Data Fig. 2) compared with previously published in vitro structures.This high fidelity in density was observable in both the final 3D reconstructions and the post-3D-refinement 2D class averages, enabling direct modelling of native membrane structures.
The model building for the inner membrane structures surrounding the mitochondrial SCs involved a four-step procedure.First, discrete points were sampled from the raw signals in a given density map-such as the type-A SC-on the basis of binarized membrane density.A 2D plane was fitted by least-square minimization; the normal vector of each SC was estimated and the coordinate system was rotated so that this vector aligned with the z axis.Second, these sampled discrete points were used to generate two smooth, curved surfaces with a thickness of around 4 nm.Third, planar phospholipid bilayer structures were generated to match the geometry of these estimated surfaces.Finally, the information from the second and third steps was integrated to geometrically deform each planar membrane structure into a smooth, curved surface.
To optimize the initial sampling for membrane model building, we categorized the membrane structures surrounding the protein into three distinct groups: structured lipids, surface-associated lipids and generic bilayer lipids.The first category, structured lipids, included lipids that are closely associated with the transmembrane regions of the protein.This close association enabled identification and direct atomic-level modelling of these specific lipid species, which have also been observed in previously reported structures purified using detergent.The second category, surface-associated lipids, comprised lipids situated around the immediate periphery of the protein, forming a pseudo-lattice structure.Within this lattice, partial phosphatidyl head groups and hydrophobic tails could be discerned.Our in situ density maps allowed us to unambiguously determine the locations of individual lipids in this category; however, the current quality of the density maps does not permit identification of the specific types of lipid present.The third category, generic bilayer lipids, represented a region farther from the protein where only the density features corresponding to the bilayers could be observed.We used a generic phospholipid membrane model to approximate the probable horizontal positions of the phosphatidyl headgroups.Owing to the fluid nature of the lipid bilayer and the high level of noise in the density maps, the central positions of these generic bilayer lipids may still vary among different subclasses even after focused classification.However, the average geometric features and the central locations of the membranes were notably consistent across each of the four main types of SC.Therefore, these generic bilayer lipids were used solely for calibrating the central locations and orientations of the phospholipid bilayer, rather than representing the actual positions of individual phospholipid molecules within the bilayer of each SC.This approach facilitated analysis of the overall geometric changes among the SCs, albeit not at the level of individual phospholipid molecule structures.
To achieve a sufficiently smooth model for the generic bilayer lipids, we performed real-space refinement of the initial structures using the Coot software 70 .The refined structures were subjected to further smoothing using a local Gaussian filter to minimize residual noise in localized membrane regions.This step enabled precise estimation of the contour map and the local curvature at each point (Fig. 2b).We used the CHARM-GUI web service 71 to generate a simulated rectangular planar phospholipid bilayer.This planar structure was then mapped on to the curved surfaces that were obtained after Gaussian smoothing.This mapping process yielded a curved membrane model that optimally fit the density map.From these estimated surfaces, information about the local geometry of the membranes surrounding the mitochondrial SCs could be directly retrieved for subsequent geometry analyses and comparisons.

Model building, refinement and validation
The atomic models were built manually using Coot 72 .First, highresolution structures of bovine CI (PDB: 7QSK), bovine CIII 2 (PDB: 2A06) and bovine CIV (PDB: 5XDQ) were fitted into the corresponding map as a rigid body using ChimeraX 64 .Then, the fitted model was manually mutated, adjusted and real-space refined to correct errors in local regions to best match the density maps using Coot 72 .The final model was refined using phenix.real_space_refine 73with geometric constraints and validated using MolProbity 74 .Figures were generated using UCSF ChimeraX 64 and PyMOL.

Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.Extended Data Fig. 9 | Distribution of the conformational states of CI and CIII under different treatments.a, CIII utilizes QH 2 from either CI or II as its substrate.b, In the 'fresh' sample, approximately 75% of CI within SCs is observed to be in an active state.This contrasts with the stability of the active state under varying conditions, where only about 30% and 18% of CI retain their active configuration under 'mild' and 'harsh' experimental conditions, respectively.Notably, in the 'fresh' sample, CI is primarily found in its active form.Simultaneously, within CIII, the ISP domain exhibits a distribution between the b-position (76%) and c-position (24%).c, In the 'fresh' sample, approximately 25% of CI is found in a deactive state.This proportion increases significantly under experimental conditions, with about 70% of CI in a deactive state under 'mild' conditions and further rising to 82% under 'harsh' conditions.Concurrently, the distribution of the ISP domain within CIII in the b-position changes to 60.4% under 'mild' conditions and decreases to 23.3% under 'harsh' conditions.Conversely, the prevalence of the c-position within CIII shifts to 39.6% under 'mild' conditions and escalates to 65.6% under 'harsh' conditions.

Fig. 1 |
Fig. 1 | In situ single-particle cryo-EM and cryo-ET analysis of mammalian mitochondrial respiratory SCs.a,b, Grid preparation of porcine mitochondria.Mitochondria were extracted from porcine hearts treated under different conditions (fresh, mild and harsh) (a) and directly frozen on to cryo-EM grids (b) for subsequent in situ imaging.c, A representative image of porcine mitochondria under low magnification.Hole diameter, 2 μm.d, Representative cryo-EM micrograph of mitochondria for single-particle analysis (SPA).Scale bar, 20 nm.e, Representative 2D class averages showing different types of SC after 3D classification.f,g, Representative tomographic slices (f) and corresponding subvolume averages (g) of a reconstructed tomogram (Tomo).Scale bar, 100 nm.h,i, Side (h) and top (i) views of a representative high-resolution map of an SC in the native mitochondrial inner membrane.j, A molecular model of the high-resolution SC with the surrounding membrane built.k, High-resolution features shown by representative density of amino acid residues and endogenous ligands.CDL, cardiolipin.

Fig. 2 |
Fig. 2 | Architecture, membrane curvature and interaction interfaces of the four types of respiratory SC. a, Top and side views and cartoon models of SCs I 1 III 2 IV 1 , I 1 III 2 IV 2 , I 2 III 2 IV 2 and I 2 III 4 IV 2 with models of the surrounding membranes.These views highlight the impact of different SC compositions on local membrane curvature.b, Contour maps of the native membrane around the four main types of SC, viewed from the mitochondrial matrix side.Red and blue indicate high and low altitudes, respectively.These gradients clearly demonstrate the common feature in which the membranes surrounding the CI heel and CIII 2 regions are convex and concave, respectively.They also show distinct local curvature differences among the four types.c, Interaction interfaces among CI (pink), CIII 2 (light yellow) and CIV (light blue) in type-A SC.Blue mesh represents the density maps of lipid molecules filling the interstitial space among these complexes.d, Representative local density maps and atomic modes of lipids built in the interface between CIII 2 and CIV.

Fig. 3 |
Fig. 3 | High-resolution in situ structures reveal multiple Q/QH 2 binding states within the Q-binding channel.a, Cartoon models of CI active-apo state, with subunits constituting the Q-binding channel highlighted.b, Variations in distances between the Q 10 headgroup and membrane surfaces across fully occupied, two intermediate, and half-occupied states.The height of the membrane surface was estimated using the average of lipid headgroups surrounding the CI 'heel'.c, Spatial positioning of Q 10 within the Q-binding channel for the four distinct binding states.d, Comparisons of three other binding states with the fully occupied Q 10 (transparent stick).The distance from the Q 10 headgroup to the quinone-binding channel entrance varied among different binding states.e, Schematic depiction of the silkworm-like undulatory motion of Q 10 within the Q-binding channel.

Fig. 4 |
Fig. 4 | Dynamic intermediates of CIII 2 revealed by high-resolution in situ cryo-EM.a, The structures resolve all reactive centres.b, The endogenous Q 10 and lipid molecules in the CIII 2 Q-binding pocket and around the pocket entrance.These lipids are organized in a relatively ordered manner.c, Different views of the hydrogen-bonded network and the water chains for proton transfer.Red and purple spheres denote waters in the two branches of the bifurcated proton-influx path, whereas cyan spheres indicate waters in the single-wired proton-influx path.d, Schematic illustration of (1) bifurcated proton-influx path, (2) single-wired proton-influx path and (3) proton-outflow path.e, Structural details of the Q o -binding site in the most abundant class.E271 MTCYB displays dual conformations in this class.Q 10 and amino acids are coloured purple and green, respectively.Blue spheres denote waters in the proton-outflow path near the Q o site.f, Comparison of E271 conformations between the apo-form CIII (PDB: 1NTK) and its complex with Q o inhibitor stigmatellin (PDB: 2A06) from the same view as shown in e. g, Structural details of the Q-binding site in state I with the closest distance to [2Fe-2S].h, Four main Q-bound states of CIII as revealed by high-resolution in situ structures.The Rieske domain sequentially moves from b position to c position to shuttle electrons.This movement is coupled with the Q-binding state at the Q o site.Density maps for Q 10 headgroups are shown as blue mesh, whereas those for the Rieske domain and cyt c 1 are represented as transparent surfaces.i, Schematic representation of the coupling of ISP movement to Q-binding states.

Extended Data Fig. 1 |
Schematic Overview of Initial Single-Particle Cryo-EM Data Analysis Workflow.a, High-magnification micrograph exemplifying a typical single-particle cryo-EM sample.b, Preliminary results of reference-free 2D classification.c, Representative outcomes of 2D sub-classification focused on particles displaying discernible side views of mitochondrial membranes.d, Illustrative 2D classification of mitochondrial supercomplex particles derived from the prior 2D classes; features reveal distinct concave membrane morphologies in the CIII 2 regions.e, Exemplary 2D class averages after extensive particle sorting and membrane signal weakening.f-h, 3D reconstructions generated using particles selected from steps (c), (d), and (e), respectively.i, Typical results of the iterative cross-3D classification.Reference maps and particle classification/alignment parameters were progressively refined over multiple cycles.Upon reaching this stage, reliable 3D reference models were obtained for subsequent multi-level refinement and focused classification.Particles classified into SCs were merged for all ensuing data processing steps.Extended Data Fig. 2 | Evaluation of the 3D classification by post-3D-refinement 2D classification.a-d, Atomic models (top) from three representative views, the corresponding projections from the same orientations (middle) and representative reference-free 2D class averages of type-A (a), -B (b), -O (c) and -X (d) SCs after 3D classification and refinement.Extended Data Fig. 3 | Water molecules in CI, CIII 2 , and CIV.a-c, The water molecules bound in CI, CIII 2 and CIV.Blue dots represent water molecules in the hydrophilic regions, while red dots indicate water molecules near the hydrophobic regions.The core subunits of CI, CIII 2 and CIV are color-coded as shown in the figure.Extended Data Fig. 4 | Endogenous cofactors and representative structured lipids revealed by the high resolution in-situ cryo-EM structures.a, Representative density maps of endogenous cofactors and atomic models fitted.The cartoon in the center represents the contour of the type-A supercomplex projection with the ligand position assigned.b, Structured and associated lipids identified from the high-resolution cryo-EM density of type-A supercomplex.Extended Data Fig. 5 | Complex hydrogen-bonded networks near the Q i sites for proton transfer in CIII 2 .a, These networks consist mainly of water molecules, supported by their interacting residues, and polar headgroups of lipids, forming the proton uptake path near the Q i site and the proton release path near the Q i site.b, The high-resolution map enables us to confidently build water molecules in Q i site and their surrounding residues and lipids that form these intricate hydrogen-bonded networks.c, The proton-transfer water chain and local density map near the Q o Site.Extended Data Fig. 6 | Influence of supercomplexes on surrounding membrane curvature.a-d, Side views illustrating the impact of SCs on membrane curvature.Extended Data Fig. 7 | Different Q-binding states and conformational changes in the Q-channel.a, Superpositions of the four different Q-binding states with the active-apo state (grey) reveal not only conformational alterations in the Q 10 headgroup and H59 NDUFS2 (left panel) but also significant changes in long-range structures away from Q 10 headgroup (right panel).b-e, Superimposition of the porcine in-situ supercomplex structure with the bovine CI (grey) incorporated into MSP nanodisks indicates that our fully occupied state resembles the active form of bovine CI with a Q 10 fully occupying the Q-site (State-α).Residues interacting with Q 10 in the other three binding states (State-β, -γ and -δ).

Data Fig. 10 |
Comparison of global conformational changes between active/deactive states from different studies.For clarity, the centers of the core subunits are used for studying the conformational changes.The centers are calculated by averaging the coordinates of C α of each subunit.a-c, Displacement distances of hydrophilic and hydrophobic core subunits in CI active-apo (cyan) and deactive-class7 (this work, porcine, pink) (a); CI active-apo (7QSL, cyan) and deactive-apo (7QSN, bovine, pink) (b); and CI Native-closed (6ZKO, cyan) and CI deactive-open1 (6ZKS, ovine, pink) (c).d, Bar chart statistics of the distances in a, b, and c. *Because the ND6-TMH4 has a big conformational change, this distance comparison is calculated without ND6-TMH4.