Protein-induced membrane strain drives supercomplex formation

  1. Maximilian C Pöverlein
  2. Alexander Jussupow
  3. Hyunho Kim
  4. Ville RI Kaila  Is a corresponding author
  1. Department of Biochemistry and Biophysics, Stockholm University, Sweden
20 figures, 5 videos, 8 tables and 1 additional file

Figures

Structure and function of Complex I, Complex III2, and the SCI/III2.

(A) The structure of the OXPHOS complexes with lipid headgroups from atomistic molecular dynamics (MD) simulations. (B, C) Overview of the structure and function of redox-driven proton pumping in CI and the Q-cycle of CIII2. Molecular models of the (D) SCI/III2, (E) CI, and (F) CIII2 used in the atomistic MD simulations (water and ions were omitted for visual clarity).

SC formation affects the concentration of lipids and the membrane thickness.

(A–C) OXPHOS complexes viewed from the N-side of the membrane. Membrane shift induced by (D) the SC, and (E, F) the isolated CI and CIII2. The membrane shift relative to the average membrane plane is colored by the shift in z-position viewed from the N-side. (G) Local membrane thickness around SC (top), CI and CIII2 (bottom). (H) Local cardiolipin concentration. (I) Local concentration of Q (Q + QH2). (J) Cryo-electron microscopic (cryo-EM) difference map. (K) Membrane shift determined from the cryo-EM map. (L) Membrane thickness determined from cryo-EM map. (M) Density of Q/QH2 (headgroup) from coarse-grained molecular dynamics (cgMD) simulations. Insets: Protein–Q interactions from atomistic molecular dynamics (MD). (N) Density of cardiolipin (headgroup) from cgMD simulations. Insets: Protein–CDL interactions from atomistic MD.

Conformational dynamics of the OXPHOS proteins and the SC.

(A) Subunits comprising the SC interface. Inset: Interaction near the DED-loop. (B) Normal modes derived from essential dynamics analysis for the SC (left), CI (middle), and CIII2 (right). See Videos 15 for normal modes of the SC, CI, and CIII2. (C) Decomposition of interaction energies within the SC. (D) Distribution of the opening/closing mode (left) and the twisting mode (right) for the isolated CI (in orange) and the CI within the SC (in blue). (E) Differences in the dynamics of CIII2 for minimum and maximum values of Mode 1. (F) Distribution of the CI domain angle in the SC (in blue) and CI (in orange). (G, J) Conformational changes in the Rieske subunit between SC (in blue) and CIII2 (in green) affect (G) the QH2 binding in the proximal Qo site and (J) the Rieske FeS-heme c1 distance in the distal monomer. (H) Distribution of the dihedral angle in TM3 of ND6. (I) Overview of the SC showing the locations of ND6 and Rieske subunits, as well as the proximal and distal protomers of CIII2.

Lattice model of SC formation and crowding effects in the IMM.

(A) Possible specific and non-specific interactions between CI (blue/white) and CIII2 (red) and their respective energies in the lattice model. (B) SC population as a function of the specific interaction energy (Especific) and molecular strain (Estrain). (C) Specific and non-specific assemblies as a function of the protein–lipid ratio with varying specific interaction energies (Especific). (D) Representative protein arrangement in crowded IMMs. (E) Average edge-to-edge distance distributions and nearest neighbor distance (red line) for specific protein–Q/protein contacts. (F) Nearest neighbor distance between Q and CI as a function of the Q/QH2 ratio. The dashed line indicates the 7 nm distance between active sites in the SC.

Proposed thermodynamic and functional effects of SCs.

(A) The individual OXPHOS complexes, CI and CIII2, induce prominent molecular strain in the surrounding lipid membrane (dark green area) due to a hydrophobic mismatch between the protein and the membrane. Condensation of locally strained membrane patches entropically drives the SC formation and leads to an overall reduction of the membrane strain (light green area around the SC), favoring the accumulation of cardiolipin around the SC. (B) During normal respiratory conditions, each OXPHOS complex is surrounded by multiple (ca. 6) quinol/quinone molecules that can act as substrates for the proteins. At limiting QH2 concentration, the quinol diffusion between CI and CIII2 becomes rate-limiting and leads to a kinetic advantage of the SC (see also Figure 4F).

Appendix 1—figure 1
Analysis of membrane properties.

Mismatch of the hydrophobic region of CI and CIII2 relative to the lipid membrane. (A) The hydrophobic region in CI (left) and CIII2 (right), with hydrophobic amino acids in gold, and hydrophilic amino acids in cyan. The maps were created using ChimeraX. (B) Ion distribution as a function of z position in a protein-free slab of simulation S1. (C) Dielectric constant (at z = 0.5 nm from the membrane–water interface) around the SC (left) and CI (right), with the SC leading to larger membrane surface with a perturbed ε. (D) Distribution of the membrane thickness in a POPC/POPE/CDL/QH2 membrane and a CDL membrane. (E) Dielectric profile (ε) computed perpendicular to the membrane plane for a POPC/POPE and a pure CDL membrane. The Stern layer, characterized by the local increase in the ε, extends ca. 1.5 nm from the membrane plane and is affected by the lipid composition. (F) Position of sampled dielectric profiles shown in (GL).

Appendix 1—figure 2
Concentration of lipids and quinones, and analysis of membrane thickness in atomistic molecular dynamics (aMD) simulations.

(A) Local membrane thickness in the SCI/CIII2 (top), CI (middle), and CIII2 (bottom). (B) Local cardiolipin concentration. (C) Local QH2 concentration. Difference in the membrane thickness around the SC relative to CI (left) or relative to CIII2 (right) from (D) aMD and (E) coarse-grained molecular dynamics (cgMD). (F–H) Error bars for the local CDL and Q/QH2 enrichment (data presented in Figure 2H, I) were estimated by averaging over the membrane in the vicinity of the protein (within 1.5 nm distance) divided into six domains (see panel F). The average concentrations for each domain and the calculated errors of (G) CDL and (H) Q/QH2.

Appendix 1—figure 3
Quinone dynamics in the membrane.

(A) Flip-flop motion of the Q headgroup in the membrane from atomistic molecular dynamics (aMD) simulations. Lifetime spent in upper/lower leaflets indicated by bars above. (B) Lifetime analysis of the Q and QH2 headgroup, with the half-life (t1/2) for Q/QH2 spent on a given membrane leaflet during consecutive steps during the simulations. (C) Representative aMD structure of QH2 headgroup interaction with water and lipid headgroups at the membrane interface. (D–F) Average diffusion of Q/QH2 around the SC/IIII2, CI, and CIII2 from coarse-grained molecular dynamics (cgMD) simulations. (G–I) Concentration gradient of Q/QH2 around the SCI/III2, CI, and CIII2 from cgMD simulations. CI – in blue, CIII2 – in beige. The Q species show a directed rotational diffusion around the OXPHOS complexes that could arise from conservation laws due to the Q depletion/increase in unique regions.

Appendix 1—figure 4
Shifts in the membrane leaflet from coarse-grained molecular dynamics (cgMD) simulations.

Shifts in the membrane leaflet viewed from the N-side of the (A) SCI/III2, (B) CI, and (C) CIII2. Shifts in the membrane leaflet viewed from the P-side of the membrane for (D) CIII2 and (E) SCI/III2. The location of the Qo site of the proximal CIII protomer is indicated by a red star. (F) The CIII2 angle was measured between the z-axis of the simulation box and the principal axis of inertia of the protein. (G) Time evolution of the membrane angle averaged over atomistic molecular dynamics (aMD) simulations of the SCI/III2 and CIII2. (H) Projection of the principal axis onto the xy-plane. (I) Visualization of the membrane distortion effect based on trajectory-averaged membrane position.

Appendix 1—figure 5
Analysis of membrane-induced distortion effects.

(A) Relative strain effect relative to a lipid membrane from atomistic MD simulations of the SCI/III2, CI, and CIII2, suggesting reduction of the membrane strain (blue patches) in the SC surroundings. The figure shows the non-bonded energies relative to the average non-bonded energies from membrane simulations (simulation M4, Appendix 1—table 2). (B) The lipid strain contribution for different lipids calculated from non-bonded interaction energies of the lipids relative to the average lipid interaction in an IMM membrane model (simulation M4). The figure shows the relative strain contribution for nearby lipids r < 2 Å, in color from panel (C), and lipids >5 Å from the OXPHOS proteins. (C) Selection of lipids (<2 Å) interacting with the OXPHOS proteins. (D) Potential of mean force (PMF) of membrane thickness derived from thickness distributions from coarse-grained molecular dynamics (cgMD) simulations of a membrane, the SCI/III2, CI, and CIII2. (E) Membrane thickness as a function of CDL concentration from cgMD simulations. (F) ΔGthick of the SC as a function of membrane thickness based on cgMD simulations. (G) Membrane curvature around the SCI/III2 (left), CI (middle), and CIII2 (right) from atomistic simulations. (H) Squared membrane curvature obtained from cgMD simulations, within a 20-nm radius around the center of the system. These maps correspond to the curvature field used in the calculation of the bending deformation energy term (Gcurv).

Appendix 1—figure 6
Analysis of lipid end-to-end distance from atomistic molecular dynamics (aMD) simulations of (A) SC, (B) CI, and (C) CIII2.
Appendix 1—figure 7
Local lipid concentrations from coarse-grained molecular dynamics (cgMD) simulations.

Concentration of (A) POPC, (B) POPE, (C) CDL, (D) QH2, and (E) Q from cgMD simulations of the SCI/III2 (left), CI (middle), and CIII (right).

Appendix 1—figure 8
Surface area per lipid and local distortion of membrane leaflet.

(A) The area per lipid on the P-side of the membrane. (B) The local leaflet shift in the P-side of the membrane. (C) The area per lipid on the N-side of the membrane. (D) The local leaflet shift in the N-side of the membrane.

Appendix 1—figure 9
Diffusion of Q molecules in coarse-grained molecular dynamics (cgMD) simulations.

Map of Q diffusion for (A) SCI/III2, (B) CI, and (C) CIII2.

Appendix 1—figure 10
Relative movement of lipids in coarse-grained molecular dynamics (cgMD) simulations.

Map of average movement of (A) CDL and (B) QH2 for SCI/III2 (left), CI (middle), and CIII2 (right). The movement is normalized relative to the local movement of POPC and POPE.

Appendix 1—figure 11
Correlated motion within the SC and allosteric network analysis.

(A, B) The Pearson correlation coefficients from Cα positions during atomistic molecular dynamics (aMD) simulations for (A) within the SCI/III2. (B) Correlation within the individual CI and CIII2. Core subunits are indicated by labels. Differences between A and B are indicated by yellow boxes. (C, D) Allosteric network analysis (see Materials and methods). (C) Path length difference of SCI/III2 network for the CI QH2 bound state and CI apo state. (D) Difference in average path length between the CI QH2 bound state and CI apo state.

Appendix 1—figure 12
Interaction area between CI and CIII2 from atomistic molecular dynamics (aMD) simulations.

(A) Structure of the SCI/III2 interface. Inset: Specific interactions at the CI and CIII2 interface. (B) Distribution of the pair distances within the SCI/III2. (C) Distribution of interface area between different subunits. (D) Distribution of the total interface area for simulations in the apo and QH2-bound states of CI. (E) Total interface area as a function of simulation time. (F) Multiple sequence alignment of UQCRC1, with the carboxylate motif forming the SC contact highlighted in red. (G) Standard deviation of the inter-subunit distance matrix for the CI apo state (upper triangle) and CI QH2 bound state (lower triangle). (H) Distance distribution between C-terminus of NDUFB7 and CIII2. No specific contacts form between these regions during the MD simulations, supporting that substitution of this region resulted in unaltered SCs.

Appendix 1—figure 13
Analysis of the TM3ND6 dihedral angle in CI.

(A) Position of the ND6 subunit within CI. Average conformation of ND6 from the last 100 ns – (B) the side view and (C) the top view of ND6. (D) Distribution of φ for SCI/III2 and CI in the apo and QH2-bound states of CI. (E) Position of the NDUFS2 subunit of CI. (F) Average conformation (over last 100 ns) of the β1–β2 loop of NDUFS2 from MD simulations for SCI/III2 (top) and CI (bottom). (G, H) Distance distribution between His59 and Asp160 of NDUFS2.

Appendix 1—figure 14
Interaction of cardiolipin and quinone/quinol with the SCI/III2.

(A) Overview of the CDL interaction sites shown in (C–H). (B) Overview of the Q/QH2 interaction sites shown in (I–L).

Appendix 1—figure 15
Enthalpy–entropy compensations from the lattice model.

(A) Free energy changes as a function of the introduced membrane strain term. (B) Enthalpic and entropic contributions to the free energy as a function of the introduced membrane strain term.

Videos

Video 1
Normal modes of the SC.
Video 2
Normal modes of the isolated CI.
Video 3
Normal modes of the isolated CIII2.
Video 4
Comparison of the normal modes between the SC and the isolated CI.
Video 5
Comparison of the normal modes between the SC and the isolated CIII2.

Tables

Appendix 1—table 1
Estimated membrane deformation energies associated with individual complexes and their supercomplex (SC).

Energies are reported as mean ± standard deviation from bootstrap analysis. The deformation energies are reported for a membrane patch with an area of ca. 1000 nm2.

ΔGcurv (kcal mol–1)ΔGthick (kcal mol–1)
CI123.6 ± 2.823.9 ± 0.6
CIII259.6 ± 4.428.0 ± 1.4
SC103.6 ± 1.649.1 ± 0.9
SC – (CI + CIII2)–79.2 ± 5.2–2.8 ± 2.0
Appendix 1—table 2
List of atomistic MD and coarse-grained molecular dynamics (cgMD) simulations.
cgMD
Simulation nameProtein modelCI ligand stateCIII proximal ligand stateCIII distal ligand stateSimulation time (µs)
A1SCI/III2apoapoapo2 × 0.5
A2apoQ in Qi, QH2 in Qoapo2 × 0.5
A3apoapoQ in Qi, QH2 in Qo2 × 0.5
A4QH2apoapo2 × 0.5
A5QH2Q in Qi, QH2 in Qoapo2 × 0.5
A6QH2apoQ in Qi, QH2 in Qo2 × 0.5
A7CIQH2n/an/a2 × 0.5
A8CIapon/an/a2 × 0.5
A9CIII2n/aapoapo2 × 0.5
A10n/aQ in Qi, QH2 in Qoapo2 × 0.5
A11n/aapoQ in Qi, QH2 in Qo2 × 0.5
Total11 μs
cgMD
Simulation nameProtein modelCI ligand stateCIII proximal ligand stateCIII distal ligand stateSimulation time (µs)
C1SCI/III2QH2Q in Qi, QH2 in QoQ in Qi, QH2 in Qo75 + 50
C2CIQH2n/an/a50 + 40
C3CIII2n/aQ in Qi, QH2 in QoQ in Qi, QH2 in Qo50 + 50
Total315 μs
Appendix 1—table 3
Simulation details of membrane models.
aMD
Simulation nameSystem size (Å3)Lipid compositionLipid/Q/QH2 moleculesSimulation time (µs)
M178 × 78 × 82POPC/POPE (1:1)2042 × 0.2
M278 × 78 × 82CDL1022 × 0.2
M3187 × 187 × 182POPC/POPE/CDL/Q (38:38:19:5)10002 × 0.35
M4187 × 187 × 182POPC/POPE/CDL/QH2 (38:38:19:5)10002 × 0.35
cgMD
Simulation nameSystemLipid bead distance (nm)Simulation time (µs)
cgM1/2/3SCI/III20.44/0.50/0.535/5/5
cgM4/5/6CI0.44/0.50/0.535/5/5
cgM7/8/9CIII20.44/0.50/0.535/5/5
cgM10–13Membrane0.44/0.47/0.50/0.535/23/5/5
Appendix 1—table 4
Non-standard protonation states used in atomistic molecular dynamics (aMD) simulations.

*The protonation states of H59NDUFS2 and Y108NDUFS2 were modeled in their neutral (His0)/deprotonated (TyrO) states with QH2 in CI.

SubunitResidues
ND1E192, E206, H247(ε), H287(ε)
ND2H25(ε), K46, H48(ε), H112(ε), H186(ε), H232(ε), K263
ND4H82(ε), H213(ε/δ), H220(ε), Lys283, H293(ε), H319(ε), H338(ε/δ), H419(ε), H422(ε)
ND5H27(ε), H56(ε), H109(ε), K119, H230(ε), H248(ε), H323(ε), H348(ε), K392, H484(ε), H509(ε), H605(ε)
NDUFS2H55(ε), H59(ε/δ)*, Y108*, D104, H150(ε), H157(ε), H190(ε), H200(ε), D292, E343, H348(ε)
ND3D66, E68, E105
NDUFS7D68, E154
NDUFS3H19(ε), H53(ε), H145(ε)
NDUFV2H9(ε), H42(ε), H99(ε)
NDUFV1H29(ε), D98, H113(ε), H116(ε), H261(ε/δ), H283(ε), H356(ε/δ), H437(ε)
NDUFS1H43(ε), D232, H255(ε), H293(ε/δ), D324, E347, H401(ε), H421(ε), H437(ε), H494(ε), H549(ε)
NDUFS8H65(ε/δ), K81, H144(ε/δ)
ND6E100
ND4LH52(ε/δ)
18 kDaH29(ε)
9 kDaH43(ε), H44(ε), H75(ε)
B8H21(ε)
B12E29
B17H67(ε), H74(ε), H83(ε), H89(ε), H127(ε/δ)
B18H3(ε), H60(ε), H81(ε/δ), H84(ε), Asp87, Glu90, H91(ε)
B22H11(ε), H25(ε), H32(ε), H50(ε), H72(ε), H75(ε), H107(ε)
AGGGH6(ε), H42 (ε/δ), H50(ε)
ASHIH66(ε), H78(ε), H104(ε), H155(ε)
ESSSH45(ε)
MNLLH10(ε), H13(ε)
Appendix 1—table 5
Estimation of volume changes from coarse-grained molecular dynamics (cgMD) simulations.

The volume change of the membrane with embedded OXPHOS proteins. The SCI/III2 formation leads to the membrane strain relative to the individual CI and CIII2.

SystemMembrane volume (nm3)ΔV (nm3)
Membrane4160±1.0
CI3946±0.8214
CIII24093±0.768
SCI/III23865±0.7296
ΔV14 nm3
Appendix 1—table 6
Estimation of protein copy numbers and effective surface area in the IMM.

The data is based on Schlame, 2021, see also Petrache et al., 2000; Morgenstern et al., 2021; Fedor and Hirst, 2018. The crowding model, with a square length of 163 nm, describes the flat membrane regions of the IMM, thus excluding the area of ATP synthase and 10 lipids bound to the c-ring, located at the cristae (cf. Petrache et al., 2000).

ComponentCopy numberSurface area/molecule (nm2)Total surface area (nm2)
CI131902470
CII1516240
CIII2181101980
CIV64644096
ATP synthase212004200
ATP/ADP carrier160101600
Lipids48,3000.716,905/leaflet
of which Q/QH2 1%483
291 proteins
48,300 membrane
14,586 nm2 protein
16,905 nm2 membrane
Total: 31,491 nm2
Appendix 1—table 7
Energy terms in the lattice model.

The protein–protein interaction is described by specific interactions term (Especific <0 kBT) and non-specific interactions (Enon-specific >0). The membrane–protein interaction determines the strain energy of the membrane (Estrain), based on the number of neighboring ‘lipid’ occupied grids that are in contact with proteins (Figure 4A). The interaction between the lipids was indirectly accounted for by the background energy of the model. The proteins can occupy four unique orientations on a grid ([North, East, South, West]). The table summarizes the unique energies linked to the respective microstates.

StateCI and CIII2 positiond(CI – CIII2)Energy
1Neighbors with specific interaction1Especific + 10 Estrain
2Neighbors with non-specific interaction1Enon-specific + 10 Estrain
3Diagonal neighbors212 Estrain
4Separated by one lattice position in either cardinal direction213 Estrain
5Separated by one lattice position in cardinal direction514 Estrain
6Diagonal neighbors separated by one lattice position2215 Estrain
7No interaction>2216 Estrain
Appendix 1—table 8
Monte Carlo (MC) simulations of the lattice model.

The conformational landscape was sampled by MC using 107 MC iterations with 100 replicas. Temperature effects were modeled by varying β, and the effect of different protein-to-lipid ratios by increasing the grid area. The following simulations were performed (energy units are given in kB × 310 K).

SimulationTemperature (K)Grid size, NEspecificEnon-specificEstrain
G1–G5300, 305, …, 3204−110.1
G6–G10300, 305, …, 3204−110.5
G11–G15300, 305, …, 3204−111
G16–G20300, 305, …, 3204−112

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Maximilian C Pöverlein
  2. Alexander Jussupow
  3. Hyunho Kim
  4. Ville RI Kaila
(2026)
Protein-induced membrane strain drives supercomplex formation
eLife 13:RP102104.
https://doi.org/10.7554/eLife.102104.4