Structural basis of mitochondrial membrane bending by the I–II–III2–IV2 supercomplex

1Mitochondrial energy conversion requires an electron transport chain (ETC) that generates a proton motive force across the inner mitochondrial membrane to drive the essential adenosine triphosphate (ATP) formation by F1Fo-ATP synthase. The ETC consists of four multisubunit membrane complexes: complex I (CI; NADH:ubiquinone (UQ) oxidoreductase), complex II (CII; succinate:UQ oxidoreductase), complex III (CIII; cytochrome bc1) and complex IV (CIV; cytochrome c oxidase). Biochemical and structural analyses have shown that these components can organize into supercomplexes containing CI, CIII dimer (CIII2) and CIV,. CII transfers electrons from succinate via its covalently bound flavin adenine dinucleotide (FAD) and iron–sulfur clusters to UQ and is also a component of the tricarboxylic acid cycle, making a functional link between the two central metabolic pathways. Although CII has been suggested to interact with mammalian ETC complexes–, it was not experimentally found as part of any characterized supercomplex. In addition, for the bioenergetic process to occur, a specific topology of the crista membranes that form functionally distinct high-potential compartments is critical. An established mechanism for maintenance of such a topology relies on oligomerization of ATP synthase and its specific interplay with lipids–. In ciliates, the inner mitochondrial membrane is organized as tubular cristae, which was previously explained by the helical row assembly of ATP synthase,,.

2We purified the intact respiratory supercomplex from mitochondria of the ciliate protist Tetrahymena thermophila and determined its structure by single-particle cryo-electron microscopy (cryo-EM) (Extended Data Fig. 1 and Supplementary Table 1). At an overall resolution of 2.9 Å, the structure revealed that CI, CII, CIII2 and a CIV dimer (CIV2) associated into a 5.8-MDa supercomplex (Fig. 1). When viewed along the membrane plane, the assembly of more than 300 transmembrane helices displays a bent shape, indicating that the accommodating membrane adopts a local curvature with a radius of approximately 20 nm (Fig. 1). Masked refinements resolved individual structures that together form an assembly of 150 different protein subunits and 311 bound lipids (Extended Data Figs. 1 and 2a–f and Supplementary Tables 1 and 2). CIV2 is associated with the long side of the membrane region of CI, opposite CIII2. This arrangement is markedly different compared with known mammalian supercomplexes, and correlates with the acquisition of four ciliate-specific CI subunits that would clash with the position of CIV as seen in mammals (Fig. 1d and Extended Data Fig. 3). CII is anchored between CI and CIV, highlighting the unique architecture and composition of the native supercomplex (Fig. 1). In-gel activity assays confirmed co-migration of functional CI, CII and CIV in a high-molecular-mass band, which we assigned to the intact supercomplex (Extended Data Fig. 2g; and see source data for Extended Data Fig. 2).

3CIV2 is the most divergent of the four ETC complexes (Extended Data Figs. 4 and 5). We modelled 105 lipids, four UQs and 53 protein chains per monomer, of which four are mitochondrially encoded. We found that two of those subunits, previously annotated as ciliate-specific Ymf67 and Ymf68 (also known as COX3) represent complementary protein fragments, with coding genes split in the mitochondrial genome by a tRNATrp gene insertion (Extended Data Fig. 6). Each fragment is extended by over 400 and 200 residues, respectively. Together, they form a functional COX3 (COX3a and COX3b), including the conserved seven transmembrane helix fold (Extended Data Fig. 6e). In our structure, COX3a and COX3b extend throughout the CIV membrane region, and COX3b has evolved interactions with CI subunits on the matrix side, thereby mediating the supercomplex assembly (Fig. 2a,b). In particular, COX3b forms a contact with a peripheral amphipathic helix of NDUCA1, which is part of a zinc-free γ-carbonic anhydrase (γ-CA) heterotrimer (Fig. 2c). The γ-CA heterotrimer was previously reported in viridiplantae and ciliates (both diaphoretickes), and our structure demonstrates that it acts as a structural scaffold within the supercomplex architecture. Another CI–CIV contact at the same site involves COXTT2 with a N-terminal globin-like domain that interacts with NDUFA3 and was not resolved in the individual CIV2 structure, suggesting that it becomes ordered to mediate supercomplex formation (Fig. 2a,c). The second major interaction site at the CI–CIV interface, which is on the lumenal side of the membrane, also involves a fragmented protein subunit, this time from CI. Consistent with the observation with respect to the protein splitting in CIV, we modelled the N-terminal extension of the ND5 fragment (ND5a), as well as the newly identified protein subunit NDUTT16 (from CI) (Fig. 2a,b,d). NDUTT16 engages in interactions with at least four subunits of CIV, as well as an interfacial CIV haem group (Fig. 2d and Extended Data Fig. 5f).

4Our finding of the split core subunits gaining a capacity of establishing intercomplex contacts to stabilize the supercomplex that curves the membrane suggests a putative mechanism by which gene fragmentation, followed by its expansion, might convey subunit function, which was also observed in mitoribosomes,. Overall, the CI–CIV2 interface involves 25 subunits, forming an extensive buried interface of approximately 2,300 Å2 with a curved membrane region (Fig. 2a,b), and thus a single mutation is unlikely to disrupt the curvature induction.

5T. thermophila CII binds in a wedge-shaped gap formed by CI–CIV in our structure (Fig. 1a–c). In addition to the four canonical subunits (SDHA–D), it is composed of 11 ciliate-specific subunits SDHTT1–11 (Fig. 2e and Extended Data Fig. 7). The matrix module SDHA and SDHB forms a conserved head region, containing both the covalently bound FAD and the three iron–sulfur clusters (Extended Data Fig. 7b). The membrane anchor is formed by two small subunits, SDHC (7 kDa) and mitochondria-encoded SDHD (5 kDa), which could only be assigned by locating topologically conserved transmembrane helices in the map (Extended Data Fig. 7a,c,d). At the lumen, an approximately 70-kDa module (SDHTT2, SDHTT4, SDHTT5 and SDHTT8) anchors CII to CI–CIV (Fig. 2e,f and Extended Data Fig. 7a). SDHTT5 interacts with a helix of NDUTT5 protruding from the membrane arm of CI, and with the Surf1-like protein subunit COXTT8 (CIV) (Fig. 2a,e). At the same position, COXTT8 interacts with CI via the N-terminal ND5a extension, and with Ymf75, COXTT27 and COXTT18 at the CIV dimer interface. Thus, the three complexes—CI, CII and CIV—are connected in the lumen (Fig. 2f).

6Between SDHB and SDHD, we identified a ligand, which we assign to ubiquinone bound to the conserved proximal Qp site, (Extended Data Fig. 7e). On the matrix side, the 36-kDa soluble subunit SDHTT1 contains a bis-histidine C-type haem covalently bound by a single cysteine residue (Extended Data Fig. 7a). Although it is exposed to the membrane region, at a distance of approximately 60 Å from the Fe3S4 cluster, the non-canonical haem C is located too far to participate in direct CII electron transfer (Extended Data Fig. 7a). To provide further evidence of the presence of additional haem groups in the supercomplex, we recorded visible absorption redox spectra of the purified sample (Extended Data Fig. 8. Deconvolution of the merged absorption bands characteristic of B-type and C-type haems clearly revealed the contribution of at least one additional haem group in the supercomplex, in addition to those present in CIII, with specific absorption at 556 nm.

7The presence of a functional ETC with CII is consistent with previous observations that T. thermophila can utilize succinate to drive cellular respiration. Our native structure with bound CII, which contributes to the ubiquinol pool, demonstrates that assembly of the supercomplex is not limited to the proton-pumping respiratory complexes (CI, CIII and CIV). Beyond decreasing cytochrome c transfer distance, this suggests a potential role of supercomplex formation in mediating increased UQ diffusion, as suggested in analogous membrane systems with high protein to lipid ratios,. Furthermore, the tubular membrane morphology may require anchoring of CII to the curved supercomplex to retain it in the functionally relevant cristae and prevent diffusion to flat membrane regions.

8CIII2 in our structure is tilted with respect to CI by 37° (Fig. 3a). This tilted arrangement offsets the transmembrane region, consistent with its curved membrane environment. The interface involves 20 subunits and 19 bound native lipids interacting within the matrix, transmembrane and lumenal domains (Extended Data Fig. 9). When compared with the mammalian counterpart, CIII2 is rotated by 41° and shifted approximately 14 Å due to acquisition of four CI subunits, as well as the CIII subunit UQCRTT1 (refs.,–) (Fig. 3a and Extended Data Fig. 9a–c). This arrangement results in a specific CI–CIII2 contact with one copy of the Rieske iron–sulfur protein of CIII (UQCRFS1) interacting with the CI membrane arm (Fig. 3b). The interaction site is further augmented by a hitherto unidentified protein UQCRTT3, which interacts with the lumenal head domain of UQCRFS1 and wedges in between the interface to CYC1 (Fig. 3b).

9We traced the membrane-accessible UQ tunnel, lined by UQCR10 and ciliate-specific subunit UQCRTT2 (Fig. 3c), leading to the conserved COB haem bH, where the density for a bound (semi)-UQ was observed in the Qi site (Fig. 3d and Extended Data Fig. 9e,f). Furthermore, we observed map density features close to the two conserved haem bL groups, which probably correspond to ubiquinols bound at the Qo sites (Fig. 3d,e and Extended Data Fig. 9e,f). The distances between the Qo site, haem bL and haem bH within each CIII monomer are consistent with those observed in mammalian CIII2, with the two haem bL in COB being bridged by a non-canonical UQ. We detected density for two copies of the flexible UQCRFS1 head domain (Extended Data Fig. 9g), which contrasts with recent work that found only the head domain proximal to the CI quinone tunnel to display flexibility, whereas the distal domain at the CI interface was proposed to be non-functional in electron transport. Using focused 3D classification for the distal UQCRFS1 head domain, we then identified two classes probably representing the extremes of the head domain movement: from the B state where the Fe2S2 cluster is distanced from haem c1 to the C state where the Fe2S2 cluster is closest, (Fig. 3f,g and Extended Data Fig. 9i). This movement of the UQCRFS1 head domain is coupled to conformational changes in the unidentified UQCRTT3 protein, suggesting a potential role for this subunit in regulation of CIII2 activity. Thus, our observation of the distal UQCRFS1 head domain flexibility, together with the haem bL-wedged UQ, suggests that functional symmetry is maintained in the ciliate CIII2 despite deviation from the structural symmetry.

10To investigate whether the membrane-bending capacity of the supercomplex is biologically relevant, we performed cryo-electron tomography of isolated mitochondrial membranes. Cryo-tomograms revealed approximately 40-nm tubular cristae densely packed with helical ATP synthase rows and supercomplexes, identified by the conspicuous CI matrix arm (Fig. 4a). To elucidate the supercomplex architecture in situ, we performed subtomogram averaging and obtained a map at 28 Å resolution (Extended Data Fig. 1e). The subtomogram average confirmed the presence of the supercomplex, which fits our atomic model (Extended Data Fig. 1f), suggesting that this is the most abundant form. The appearance of a tubular membrane density in the subtomogram average suggests that the supercomplex adopts a preferred orientation, with its CI–CIV2 interface approximately aligned with the long axis of the tube (Fig. 4b). Furthermore, the curved membrane region of the supercomplex subtends an angle of approximately 130°, indicating that it contributes to the tubular shape of the cristae.

11To elucidate the membrane-shaping activity of the supercomplex, we performed coarse-grained molecular dynamics simulations. When placed into a planar lipid bilayer, the supercomplex induces a curved membrane topology, displacing the membrane by 18 nm from the original plane (Fig. 4c and Supplementary Video 1), in contrast to the protein-free lipid bilayer, which remains relatively flat (Supplementary Video 2). Furthermore, the annular lipid shell surrounding the complex in the equilibrated system displays a highly curved architecture, supportive of an active role in membrane curvature induction (Fig. 4d and Supplementary Video 3). In addition, we observed lipid pockets in the transmembrane interfaces between subcomplexes, which suggests that their maintenance is crucial for the integrity of the supercomplex (Extended Data Fig. 2a,b), as was also shown for photosynthetic membranes. Molecular dynamics simulations of CIV2 in a membrane suggest that it can also induce membrane bending in its immediate vicinity, similar to the supercomplex (Extended Data Fig. 10a,b). However, for the formation and stability of a tubular architecture, juxtapositioning of individual subcomplexes is probably necessary. This is also supported by the increased arc length of the detergent belt observed in the cryo-EM structure of the supercomplex, when compared with CIV2 (Extended Data Fig. 10b). Finally, molecular dynamics simulations of the supercomplex lacking CII showed a similar wrapping of the lipid bilayer as for the full supercomplex, with lipids filling the generated void (Extended Data Fig. 10c). This indicates that CII association may contribute to complex stability while retaining the enzyme in the tubular cristae regions, where ubiquinol is required as a substrate.

12Our results indicate that cristae shaping involves both the respiratory supercomplex and the ATP synthase that together generate membrane tubulation. Although the coiled ATP synthase rows fix the helix diameter at 130 nm (refs.,), supercomplexes serve the function of confining a narrow crista diameter of around 40 nm, which allows tight packing of cristae, thereby increasing the surface area of the bioenergetic membrane. This membrane-shaping organization of the respiratory supercomplex is markedly different from the mammalian homologue, which resides in the flat crista regions. Furthermore, because every crista represents an independent functional compartment, restriction of the crista diameter by the respiratory supercomplex probably serves to minimize the volume, thereby potentially contributing to higher local concentration of electron carriers, ensuring that proton translocation results in an increased local proton motive force, ultimately optimizing conditions for ATP synthesis. This is consistent with the observation that mutant yeast strains with large, balloon-like cristae display respiratory defects,. Thus, our findings show how respiratory supercomplexes together with other factors can organize the architecture of the bioenergetic membrane, providing a mechanism for enabling its functional specialization.

13T. thermophila cells were grown at 36 °C and harvested as previously described,, and cell pellets were resuspended in homogenization buffer (20 mM HEPES-KOH pH 7.5, 350 mM d-mannitol, 5 mM EDTA and 1x protease-inhibitor tablet) and lysed in a Dounce homogenizer on ice. Intact mitochondria were isolated first by differential centrifugation of the lysate; following membrane solubilization, the lysate was cleared at 30,000g for 30 min at 4 °C, and finally, on a discontinuous sucrose gradient with 15%, 23%, 32% and 60% w/v sucrose in buffer SEM (20 mM HEPES-KOH pH 7.5, 250 mM sucrose and 1 mM EDTA) at 14,1371g for 60 min at 4 °C in an SW28 rotor. Intact mitochondria sedimented to the interface between 32% and 60% sucrose and were collected from the gradient, snap-frozen in liquid nitrogen and stored at −80 °C. The isolated mitochondria were lysed in buffer A (25 mM HEPES-KOH pH 7.5, 25 mM KCl, 5 mM MgCl2 and 4% w/v digitonin) for 1 h on ice. This procedure was previously confirmed as a gentle solubilization method. Following mitochondrial membrane solubilization, cleared lysate was placed on a sucrose cushion (25 mM HEPES-KOH pH 7.5, 25 mM KCl, 5 mM MgCl2, 0.1% w/v digitonin and 30% w/v sucrose) in Ti70 tubes and centrifuged at 164,685g for 1 h. The pellet was gently washed and finally resuspended in buffer D (25 mM HEPES-KOH pH 7.5, 25 mM KCl, 5 mM MgCl2 and 0.1% w/v digitonin). Before loading sample material on a size-exclusion chromatography column, larger aggregates were pelleted at 30,000g for 20 min at 4 °C. Cleared sample was loaded on a Superose 6 Increase 3.2/300 column equilibrated in buffer D, collecting elution fractions of 100 µl throughout the run. Peak fractions were used immediately for cryo-grid preparation.

14The sample obtained from the sucrose cushion step was analysed for haem content and supercomplex composition using UV–visible difference spectroscopy and clear native-PAGE (CN-PAGE) combined with in-gel activity assays. UV–visible difference spectra were recorded between 390 nm and 675 nm using a home-built spectrophotometer. Protein samples were diluted as necessary in 50 mM HEPES and 0.1% digitonin at pH 8.0. Spectra were measured from sodium dithionite-reduced minus air-oxidized spectra. When multiple absorption bands overlapped, spectra were deconvoluted using the peak analysis function in OriginPro 2015 (OriginLab Corporation).

15NativePAGE 3 to 12%, Bis-Tris, 1.0 mm, Mini Protein Gel (Invitrogen) pre-cast gels were used for CN-PAGE. The gels were loaded with a protein ladder (NativeMark, Invitrogen) and four identical sample lanes where the protein samples were mixed with NativePAGE Sample Buffer (final concentration, 50 mM BisTris, 6N HCl, 50 mM NaCl, 10% w/v glycerol and 0.001% Ponceau S, pH 7.2) as per the manufacturer’s instruction. Electrophoresis was conducted at 4 °C, first at 150 V for 30 min with NativePage Light Blue Cathode buffer (50 mM BisTris, 50 mM Tricine, pH 6.8, and 0.002% Coomassie G-250) and then at 250 V for 150 min with NativePAGE Anode Buffer (50 mM BisTris and 50 mM Tricine, pH 6.8). In-gel activity assays were performed following published protocols. In brief, each sample lane from CN-PAGE was incubated with an aqueous solution to reveal (1) protein bands (0.02% Coomassie G-250, overnight) or the presence of active (2) CI (2 mM Tris pH 7.4, 2.5 mg ml−1 nitrotetrazolium blue chloride (NBT) and 0.1 mg ml−1 NADH, 15 min), (3) CII (5 mM Tris pH 7.4, 2.5 mg ml−1 NBT, 84 mM succinic acid and 0.2 mM phenazine methosulfate, 30–40 min) and (4) CIV (0.05 mM KPi pH 7.4, 0.5 mg ml−1 3,3′-diaminobenzidine (DAB) and 1 mg ml−1 cytochrome c from Saccharomyces cerevisiae, overnight). The reactions were stopped by incubation in 10% (v/v) acetic acid, followed by multiple exchanges of water. The gels shown in this paper are representative of two experiments from two separate supercomplex preparations.

16Supercomplex eluted at a concentration of approximately 10 mg ml−1. Aliquots of the peak fraction were diluted in buffer D to 0.75 mg ml−1 before applied to cryo-grids. Quantifoil R2/2-300 grids floated with a homemade 3-nm amorphous carbon layer were glow-discharged immediately before applying a 3 µl sample. Grids were vitrified using liquid ethane cooled by liquid nitrogen in a Vitrobot Mark IV, with a 30-s wait time before blotting grids for 3 s at blot force of 0. A total of 26,063 movies were collected on a Titan Krios (Thermo Fisher Scientific) operated at 300 kV at a nominal magnification of ×165,000 (0.83 Å per pixel) with a Quantum K2 camera (Gatan) using a slit width of 20 eV. An objective lens aperture of 70 µm was used, with an exposure rate of 4.26 electrons per pixel per second with 5-s exposure fractionated into 20 frames, and the defocus range was 0.6–2.6 μm.

17Motion correction was performed in the internal implementation of RELION-3.1 (ref. ), followed by contrast transfer function estimation by CTFFIND4. Initial rounds of particle picking and 2D classification, followed by ab initio reconstruction, 3D classification and preliminary refinement of the supercomplex. Template-based particle picking in RELION was then used to pick and extract 1,664,103 particles. 2D classification and 3D heterogeneous refinement steps in cryoSPARC v2 (ref. ) were then used to separate supercomplex particles from copurified ATP synthase, resulting in a final 138,746 supercomplex particles used for subsequent refinement. Following a consensus refinement in cryoSPARC, per-particle contrast transfer function refinement and Bayesian polishing were performed in RELION-3.1. For final refinements in cryoSPARC, particles were downsampled from a 724-pixel box to 480 pixels, resulting in a pixel size of 1.25 Å per pixel. Masked refinements of the respective supercomplex subregions resulted in map resolutions of 2.9 Å for the entire supercomplex, 2.8 Å for CI, 3.0 Å for CII, 2.8 Å for CIII and 2.6 Å for CIV2. Reported map resolutions are according to gold-standard Fourier shell correlation using the 0.143-criterion. To assess flexibility of the Rieske subunit wedged in the CI–CIII2 interface, we performed focused 3D classification in RELION-3.1 using pre-aligned particles with a mask on the extended area around the headgroup of the Rieske subunit. Classification into ten classes resulted in maps confirming flexibility of the structural element, with two classes corresponding closely to the previously reported B state and C state (Fig. 3f,g).

18Crude mitochondrial pellets were resuspended in an equal volume of buffer containing 20 mM HEPES-KOH pH 7.4, 2 mM EDTA, 250 mM sucrose and mixed in a 1:1 ratio with 5-nm colloidal gold solution (Sigma Aldrich) and vitrified as described above on glow-discharged Quantifoil R2/2 Au 200 mesh grids. Tilt series were acquired on a Titan Krios operated at 300 kV with a K3 camera (slit width of 20 eV) using SerialEM or the EPU software (Thermo Fisher Scientific). Mitochondrial membranes were imaged at a nominal magnification of ×42,000 (2.11 Å per pixel) and an exposure rate of 19.5 electrons per pixel per second with a 3 e− Å−2 exposure per tilt fractionated into five frames, with tilt series acquired using the exposure-symmetric scheme to ±60° tilt and a 3° tilt increment. Following motion correction in motionCor2, tomographic reconstruction from the tilt series was performed in IMOD using phase-flipping and a binning factor of 2. Tomograms were contrast enhanced using nonlinear anisotropic diffusion filtering to facilitate manual particle picking of supercomplex particles based on the matrix arm of CI. Subtomogram averaging was performed in PEET. Initial references were generated from the data by averaging after rotating subvolumes into a common orientation with respect to the membrane based on manually assigned vectors. Following initial rounds of averaging to generate a suitable reference, data were manually split into half-sets and refined independently, following low-pass filtering to 50 Å. Averaging of 360 particles from 12 tomograms resulted in a 28 Å subtomogram average.

19Manual model building was performed in Coot and new subunits identified directly for the cryo-EM map. For identified canonical subunits, homology models were generated using SWISS-MODEL. Bound cardiolipins were unambiguously identified from their head group density. Other natively bound lipids were tentatively modelled as phosphatidylcholine, phosphatidylethanolamine or phosphatidic acid based on head group densities. Real-space refinement of atomic models was performed in PHENIX using secondary structure restraints. Atomic model statistics were calculated using MolProbity.

20Given the mild solubilization conditions we used, for CIII2, the cryo-EM map showed density located on the pseudo-C2 symmetry axis between the two COB haem bL molecules, displaying planar map features consistent with the quinone moiety of UQ. The density clearly indicates that UQ can bind in two orientations, related by the symmetry rotation of the dimer. In either orientation, the quinone moiety is positioned close to a haem bL, where it potentially could accept electrons for transfer across the dimer axis. In the recent amphipol CIII2 structure, the isoprenoid tail of UQ was modelled in the equivalent position; however, the planar density for the quinone was missing. This orientation-equivalent binding of UQ between the two COB haem bL molecules, together with the B-state and C-state Rieske conformations, suggest a maintained functional symmetry of ciliate CIII2 within the supercomplex.

21In CI, we identified 49 canonical subunits and 21 subunits that we assign as phylum-specific. In each CIV monomer, we identified 11 subunits homologous to mammalian CIV (COX1, COX2, COX3a, COX3b, COX5B, COX6A, COX6B, COX6C, COX7A, COX7C and NDUFA4) and 42 ciliate-specific subunits, most of which are peripherally associated around the mitochondrial protein core. Three of the mammalian subunits missing in T. thermophila CIV (COX4, COX7B and COX8) are at the interface where two mitochondrial carriers are bound. The mitochondrially encoded core subunit COX3 is split into two fragments. Most of the TM helices are contributed by the C-terminal of COX3b, which is encoded by the mitochondrial ymf68 gene. The newly annotated Ymf68 is structurally conserved, apart from the missing helix (H1), which is structurally replaced by Ymf67. We therefore assign ymf67 and ymf68 of the ciliate mitochondrial genome as separately encoding the COX3a and COX3b subunit fragments, respectively. On the T. thermophila mitochondrial genome, ymf67 and ymf68 genes are located on the same strand, but are separated by the gene encoding tRNATrp, suggesting that a transposition event may have led to the fragmentation of the original gene encoding COX3. tRNA genes are known to be among the most motile elements in metazoan mitochondrial DNA. Both COX3a and COX3b fragments have evolved substantial subunit extensions that thread through the augmented CIV monomer unit to recruit lineage-specific subunits and mediate supercomplex assembly.

22In the CIV dimer, the dimer interface of 17,000 Å2 is dominated by 16 species-specific subunits. Furthermore, when aligned on the CIV core, a comparison of the mammalian and ciliate structures reveals that the two dimers display markedly different architectures, dimer axes and distances between COX1 cores. This suggests that the dimerization of the ciliate CIV probably evolved through the acquisition of lineage-specific subunits and reflects the constraints of the unique tubular membrane environment.

23In addition, each CIV monomer complex contains two different Surf1-like proteins, which were reported to complement defects causing Leigh syndrome in humans. In our structure, two Surf1-like proteins are permanently attached to CIV and display similar overall structures, consisting of a lumen-exposed soluble domain and a TM helix hairpin. The two Surf1 proteins are facing each other, bound on opposite sides of each CIV core.

24The presence of subunit extension and accessory subunits in CIV generates a pronounced cavity around the cytochrome c-binding site. However, overlaying of a T. thermophila cytochrome c homology model suggests that the canonical binding site is not obstructed. Cytochrome c binding is known to be driven by electrostatic interactions with the CuA domain of COX2, which forms a negatively charged patch in mammals. This structural feature is positively charged in T. thermophila, interacting mainly with H1 of cytochrome c, which displays a flipped polarity. We conclude that the experimentally observed functional incompatibility of T. thermophila CIV and mammalian cytochrome c is not due to divergent architecture, but to an inverted surface charge of the binding pocket.

25We performed coarse-grained (CG) molecular dynamics simulations on the entire T. thermophila supercomplex structure using Martini3 forcefield to study the rearrangement of the lipid bilayer around the highly bent protein assembly. Using the martinize2 (version 2.6) tool, we transformed the atomistic structure into a CG model (atoms clustered into Martini beads). Using the small molecules database and the existing topologies of phospholipids available in the Martini3 forcefield, we generated the force-field parameters for cardiolipin. Cofactors and resolved lipids in the structure were not included in the simulated model system, and only protein was simulated to study the dynamics of lipid molecules around it. First, the CG model of the protein structure was minimized for 100 steps in vacuum to remove possible steric clashes. Then, the minimized CG supercomplex was embedded in a large (75 nm × 75 nm) hybrid membrane slab (POPE:POPC:CL in 4:2:1 ratio) using the insane.py script. The CG protein–membrane system was solvated using standard Martini3 water beads and 100 mM Na+ and Cl− ions. Starting from this initial position, the simulation system was minimized keeping all beads free, first in double precision to resolve steric clashes between the lipids (maximum of 500 steps) and then in regular single precision (maximum of 10,000 steps). After minimization, with 4,000 kJ mol−1 nm−2 harmonic constraints on the backbone beads, the system was equilibrated using velocity-rescaling thermostat and Berendsen barostat for 10 ns. During production runs, the 4,000 kJ mol−1 nm−2 harmonic constraints on the backbone beads were maintained. Velocity-rescaling thermostat and Parrinello–Rahman barostat were used for temperature (310 K) and pressure (1 bar) control in the production phase. Coulombic interactions were treated with the reaction-field algorithm using εr = 15 (ref. ). The Verlet cut-off scheme was implemented with a Lennard–Jones cut-off of 1.1 nm (ref. ). The time step of the CG molecular dynamics simulations was 20 fs. Initial simulation replicas showed incomplete or unstable wrapping of the membrane around the protein, so we translated the lipid bilayer patch in z-direction and altered the insertion angle of the supercomplex to find an initial position that allowed the membrane to equilibrate and wrap fully around the protein (Systems T1-T7, simulation lengths of 0.9–2.8 µs, total of 9.6 µs). After finding the correct insertion of the protein into the membrane, we initiated three independent simulation replicas (Systems P1-P3, simulation lengths of 10 µs each, total of 30 µs). The simulations were performed using the Gromacs software (version 2021). In addition to molecular dynamics simulations of the entire supercomplex in the lipid bilayer, we also performed molecular dynamics simulations of the pure lipid bilayer (see above; three replicas, 10 µs each) and of CIV2, CI–CIII2 and CI–CIII2–CIV2 subcomplexes (1–3 replicas, 1–5 µs each).

26Images and videos were rendered using ChimeraX and Visual Molecular Dynamics. To analyse the T. thermophila cytochrome c-binding site, the mammalian cytochrome c-bound CIV structure (Protein Data Bank ID: 5IY5) was overlaid to the T. thermophila structure. Using AlphaFold2 (ref. ), a T. thermophila cytochrome c structure was predicted and overlaid to both the mammalian structure and the T. thermophila structures. The composite map of the complete respiratory supercomplex was generated in ChimeraX. This map was only used for visualization, and not for atomic model refinement; instead, a consensus map was used. The buried areas of the CI–CII–CIII2–CIV2 supercomplex and the CIV2 interface were calculated in ChimeraX.

27Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

28Any 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 10.1038/s41586-023-05817-y.