1. Introduction
Mitochondria possess a double-membrane structure, consisting of the inner mitochondrial membrane (IMM) that encloses the matrix, and the outer mitochondrial membrane (OMM) that separates the intermembrane space from the cell cytoplasm. These double-membrane organelles form a dynamic network that constantly changes shape through fusion and fission processes. Maintaining a delicate balance between these two processes, collectively known as mitochondrial dynamics, is crucial for mitochondrial function in cellular energy generation and overall cell health [
1]. In mammalian cells, the elongation of mitochondria through fusion is orchestrated by the spatio-temporal coordination of two key transmembrane proteins from the dynamin superfamily: Mitofusins (Mfn1 and Mfn2) and OPA1, which mediate the fusion of OMM and IMM, respectively [
2]. Mitofusins are composed of an N-terminal GTPase domain, followed by a first heptad repeat domain (HR1), a transmembrane domain (TMD), and a second C-terminal heptad repeat domain (HR2). These different domains are essential for Mitofusin function but their exact mode of action in mitochondrial fusion remains not fully understood [
3,
4,
5,
6].
Based on recent structural studies [
7,
8,
9,
10], the current working hypothesis for the molecular mechanisms of OMM fusion involves a
cis- (
i.e., within the same membrane) and/or
trans- (
i.e., across two different membranes) oligomerization of Mitofusin proteins. This initial event would lead to a long-distance (20–30 nm) membrane docking step [
8,
11]. Through GTP hydrolysis, membrane-bridging trans-Mitofusin complexes would then transit from an open to a closed conformation bringing OMM in close apposition. The HR2 domain could stabilize this short-distance (5–10 nm) membrane docking step by forming a homodimeric antiparallel coiled-coil complex [
5,
8,
11].
Bringing membranes in close apposition is the first step in membrane fusion, but it is not enough on its own. For fusion to occur, the membrane structure must also be destabilized to facilitate the merging of lipid bilayers. Two recent studies have emphasized the crucial role of amphipathic helices within the Mitofusin sequence in triggering OMM fusion [
12,
13]. One study found that the HR1 domain mediated liposome fusion
in vitro and was essential for mitochondrial fusion
in situ [
12]. The fusion activity was attributed to a conserved amphipathic helix located at the C-terminus of HR1, suggesting that HR1 induces fusion by interacting with and perturbing the lipid bilayer structure. A similar fusion mechanism has been described for the C-terminal amphipathic tail of Atlastin, another dynamin-like transmembrane protein involved in endoplasmic reticulum (ER) membrane fusion [
14,
15,
16]. Interestingly, another study demonstrated that when the TMD of Mitofusin was replaced with that of Atlastin, the resulting chimeric protein localized to ER membranes and was capable of mediating ER fusion [
13]. Furthermore, an amphipathic helix identified between the TMD and the HR2 domain of Mitofusin could effectively replace the C-terminal amphipathic tail of Atlastin in both
in vitro liposome fusion and
in situ ER fusion.
The lipid composition of mitochondrial membranes also plays a crucial role in facilitating mitochondrial fusion. The successive stages leading to membrane fusion involve the formation of energy-demanding intermediate membrane structures with high curvature [
17]. Lipids that can relieve this energy stress by inducing favorable membrane bending therefore facilitate fusion. Mitochondrial membranes contain specific lipids such as phosphatidylethanolamine (PE) and phosphatidic acid (PA), both of which possess a cone-shaped structure with a small headgroup area compared to the cross-sectional area of their hydrophobic chains. As a result, they can induce negative membrane curvatures when present in the outer leaflet of lipid bilayers. Additionally, there is a unique lipid found in mitochondria known as cardiolipin (CL). When bound to divalent cations like Ca
2+ or Mg
2+, CL can also adopt a conical shape. These three lipids thus have the potential to facilitate mitochondrial fusion events.
Indeed, a previous study found that high concentrations of CL were required for
in vitro liposome fusion mediated by OPA1 [
18]. CL also allows the generation of PA at the OMM through its cleavage by the mitochondria-localized phospholipase D (MitoPLD) enzyme. Depletion of MitoPLD in mammalian cells through RNA interference (RNAi) led to reduced mitochondrial fusion, suggesting that PA production at the OMM is important for fusion. Similarly, reduction of PE levels in mammalian cells by RNAi silencing of the enzyme phosphatidylserine decarboxylase (Pisd), which converts phosphatidylserine (PS) into PE, resulted in mitochondrial fragmentation, indicating inhibition of mitochondrial fusion [
19]. Yeast cells lacking Psd1, the homolog of mammalian Pisd, also exhibited impaired mitochondrial fusion [
20]. Additionally, yeast cells which lacked both Psd1 and the CL synthase Crd1 displayed an even more fragmented mitochondrial network [
21]. These studies demonstrate the critical importance of PE, CL and PA lipids in promoting efficient mitochondrial fusion. However, the specific molecular mechanisms by which these lipids facilitate mitochondrial fusion events remain to be established.
In this study, we investigate the role of PE, CL and PA lipids in membrane fusion mediated by the HR1 domain of Mitofusin. We also examine the interplay between these lipids and the divalent cations Ca2+ and Mg2+ in facilitating membrane perturbation and fusion via the amphipathic helix of HR1.
2. Materials and Methods
2.1. Chemicals
N-2-hydroxyethylpiperazine-N'-2-ethanesulfonic acid (HEPES, OmniPur grade), potassium hydroxide solution 47% (KOH 47%, EMSURE grade for analysis), potassium chloride (KCl, OmniPur grade), calcium chloride dihydrate (CaCl2, OmniPur grade), magnesium chloride hexahydrate (MgCl2, OmniPur grade), glycerol (Molecular Biology grade), tris(2-carboxyethyl)phosphine hydrochloride (TCEP, ≥98% GC), n-octyl-β-D-glucopyranoside (β-OG, ≥98% GC), n-dodecyl β-D-maltoside (DDM, ULTROL grade) and sodium dithionite (Analytical grade) were purchased from Merck. 5-(N-2,3-dihydroxypropylacetamido)-2,4,6-tri-iodo-N,N'-bis-(2,3-dihydroxypropyl)isophthalamide (Nycodenz, ≥98%) was purchased from Proteogenix.
All aqueous solutions were prepared using 18.2 MΩ ultra-pure water and filtered with sterile 0.22 µm polyethersulfone (PES) membranes.
1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (PC), 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (PE), L-α-phosphatidylinositol (Liver, Bovine) (sodium salt) (PI), 1',3'-bis[1,2-dioleoyl-sn-glycero-3-phospho]-glycerol (sodium salt) (CL), 1,2-dioleoyl-sn-glycero-3-phosphate (sodium salt) (PA), 1,2-dioleoyl-sn-glycero-3-phospho-L-serine-N-(7-nitro-2-1,3-benzoxadiazol-4-yl) (ammonium salt) (NBD PS), 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-(lissamine rhodamine B sulfonyl) (ammonium salt) (Rho PE), 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-[4-(p-maleimidophenyl)butyramide] (sodium salt) (MAL), and 1,2-dioleoyl-sn-glycero-3-[(N-(5-amino-1-carboxypentyl)iminodiacetic acid)succinyl] (nickel salt) (NTA-Ni) were purchased from Avanti Polar Lipids as chloroform solutions.
2.2. Peptides
The heptad repeat domain 1 (HR1) of human Mitofusin 1 (Mfn1) used in this study was produced by Fmoc solid-phase peptide synthesis and one-step purification by reverse-phase HPLC (Proteogenix, purity > 95%). The produced sequence was Mfn1-HR1 (T350-L420; C411S, C418S; with a C-terminal Cys or His6 tag). Lyophilized samples (1 mg aliquots) were put on ice and solubilized by slowing adding 1 mL of ice-cold buffer H (25 mM HEPES/KOH, pH 7.4; 150 mM KCl; 10% (v/v) Glycerol) containing 0.25 mM TCEP in the case of Cys-tagged HR1. Peptide solutions were vortexed for 2 min at room temperature followed by 20 sec of sonication on ice to remove any potential aggregates (2 cycles of 10 sec on, at 10 W, and 10 sec off, using the ultrasonic homogenizer UP200St from Hielscher equipped with a 2 mm sonotrode). Samples were snap-frozen in liquid nitrogen and stored at -80 °C as aliquots of 50 µL.
2.3. Liposomes
Liposomes were prepared by the detergent-assisted method [
22]. 1.2 µmol of the appropriate lipid mixtures in chloroform were dried in glass tubes for 10 min under a gentle stream of argon, followed by 2 hours under vacuum. The dried lipid films were resuspended in 400 µL of buffer H containing 1% (w/v) β-OG by vigorously vortexing for 30 min at room temperature. The detergent concentration was next reduced below the critical micellar concentration, 0.33% (w/v), by slowly adding 800 µL of buffer H, and then removed by overnight flow dialysis against 4 L of buffer H. Liposomes (final lipid concentration of 1 mM) were stored on ice and protected from light for up to 2–3 weeks.
2.4. Multi-angle dynamic light scattering (MADLS)
5 µL of liposomes at 1 mM and 95 µL of buffer H were mixed in a low volume quartz batch cuvette (ZEN2112, Malvern Panalytical) and their size distribution was determined at 37 °C with the Zetasizer Ultra Red instrument (Malvern Panalytical) using the MADLS mode, which measures the correlation function in three scattering directions: back scatter (173 degrees), side scatter (90 degrees) and forward scatter (13 degrees).
2.5. Cryogenic transmission electron microscopy (cryo-TEM)
A 5 μL drop of the initial sample solution was deposited on “quantifoil” carbon membrane grids (Quantifoil Micro Tools GmbH, Germany). The excess of liquid on the grid was absorbed with a filter paper and the grid was quench-frozen quickly in liquid ethane to form a thin vitreous ice film using an homemade mechanical cryo plonger. Once placed in a Gatan 626 cryo-holder cooled with liquid nitrogen, the samples were transferred in the microscope and observed at low temperature (-180 °C). Cryo-TEM images were recorded on ultrascan 1000, 2k x 2k pixels CCD camera (Gatan, USA), using a LaB6 JEOL JEM2100 (JEOL, Japan) cryo-microscope operating at 200 kV with a JEOL low dose system (Minimum Dose System, MDS, JEOL, Japan) to protect the thin ice film from any irradiation before imaging and reduce the irradiation during the image capture.
2.6. Liposome co-flotation assay
To assess membrane binding of HR1, 100 µL of HR1 at 50 µM and 100 µL of liposomes at 1 mM were incubated together for 1 hour at 37 °C under intermittent gentle mixing (1 min at 300 rpm every 9 min). 50 µL was taken out to serve as input control (unfloated fraction) and the remaining 150 µL was mixed with 150 µL of Nycodenz 80% (w/v) in buffer H. The mixture was transferred to a centrifuge tube (0.8 mL, Open-Top Thinwall Ultra-Clear Tube, 5 x 41 mm, Beckman Coulter) and overlaid with 250 μL of Nycodenz 30% (w/v) in buffer H followed by 50 μL of buffer H alone. The resulting gradient was centrifuged at 192,000 g for 4 hours at 4 °C in a SW 55 Ti Swinging-Bucket rotor (Beckman Coulter) and 37.5 μL of liposomes was collected from the top layer (floated fraction). 12 µL of unfloated or floated fraction was mixed with 4 µL of sample buffer (NuPAGE LDS Sample Buffer 4X, Invitrogen) and separated by electrophoresis on a polyacrylamide gel (NuPAGE 4–12%, Bis-Tris, 1 mm, Mini Protein Gel, Invitrogen) stained with Coomassie G-250 (PageBlue Protein Staining Solution, Thermo Scientific). Images were acquired using a ChemiDoc Touch Imaging System (Bio-Rad). The band intensities in the unfloated and floated fractions were measured with the software ImageJ.
2.7. FRET-based lipid mixing assay
For each condition to be tested, two sub-populations of liposomes with the same lipid composition were generated, except that one sub-population, referred to as the fluorescent donor liposomes, was labeled with the Fluorescence Resonance Energy Transfer (FRET) pair of fluorescent lipids NBD PS and Rho PE. These fluorescent lipids were added at a concentration of 1.5 mol% each at the expense of PC lipids. 27 µL of non-fluorescent acceptor liposomes at 1 mM and 21 µL of buffer H (or 19 µL of buffer H and 2 µL of CaCl
2 or MgCl
2 at 30 mM in buffer H for experiments performed in the presence of divalent cations) were added to a flat bottom 96-well white polystyrene plate (Thermo Scientific) and pre-warmed at 37 °C for 7 min. 6 µL of fluorescent donor liposomes at 500 µM were carefully added to one side of the well. 6 µL of HR1 at 250 µM were added to another side of the well. The fusion reaction was initiated by shaking the plate in order to mix the three different solutions. Lipid mixing was measured by following fluorescence dequenching of the NBD probes from the donor liposomes resulting from their dilution into the acceptor liposomes. The NBD fluorescence was monitored at 1-min intervals for 90 min (excitation at 460 nm, emission at 535 nm, cutoff at 530 nm) by the SpectraMax M5 microplate reader (Molecular Devices) equilibrated to 37 °C. After 90 min, 10 µL of 2.5% (w/v) DDM was added to completely dissolve the liposomes and thus measure the NBD fluorescence at infinite dilution, Fmax(NBD). The data were then normalized using the following equation that gives the percentage of NBD fluorescence increase at time t, %F(NBD, t):
where Fmin(NBD) is the lowest NBD fluorescence value from all time points.
2.8. Sodium dithionite assay
To determine the proportion of hemifused versus fused liposomes, we employed a fluorescence quenching method using sodium dithionite. 5 µL of freshly prepared sodium dithionite solution at 100 mM in buffer H was incubated with 33.3 µL of fluorescent donor liposomes at 1 mM for 15 min at 37 °C. 28.3 µL of buffer H was then added to the mixture to dilute the liposomes to a final concentration of 500 µM. The FRET-based lipid mixing assay was then performed as described above. As sodium dithionite was shown to completely lose its activity after 10 min at 37 °C, this ensures that the fluorescence signal from NBD lipids of the inner leaflets remains unquenched, even when DDM is added to solubilize the liposomes.
The percentage of liposomes that underwent hemifusion at time t, H(t), is given by the equation:
where F
T and F
I are respectively the normalized fluorescence dequenching signals without and with prior sodium dithionite treatment (total lipid mixing and inner monolayer lipid mixing, respectively), and α is the proportion of lipids residing in the inner monolayer of liposomes.
4. Discussion
In a previous study, we identified an amphipathic helix at the C-terminal end of the HR1 domain. We proposed that this amphipathic helix induces HR1-mediated fusion by perturbing the lipid bilayer structure, especially in membrane regions presenting lipid packing defects, caused by either high membrane curvature or the presence of lipids with a cone-shaped structure [
12]. In this study, through cryo-EM observations, we noticed that pure PC liposomes displayed locally disappearing membrane structure after incubation with HR1 (
Figure 1d). This finding is consistent with a perturbation of the lipid bilayer structure when HR1 interacts with the liposome membrane.
The primary objective of this study was to examine the influence of membrane lipid composition, both in the presence and absence of divalent cations (Ca
2+ or Mg
2+), on HR1-mediated liposome fusion. In the absence of cations, the inclusion of 30 mol% PE lipids in the liposome membrane resulted in a more than 2-fold increase in the extent of lipid mixing mediated by HR1 compared to liposomes composed of pure PC lipids (
Figure 1b,c). We attribute this effect to the cone-shaped structure of PE, which induces packing defects in the membrane structure and, consequently, enhances membrane binding and perturbation by the amphipathic helix of HR1. The presence of PE in the liposome membrane also promoted hemifusion events over full fusion mediated by HR1. The percentage of liposomes undergoing hemifusion increased from 50 to 80% when the liposomes contained 30 mol% PE lipids (
Figure 5b). This effect can also be associated with the cone-shaped structure of PE, which is known to induce the transition from bilayer to non-bilayer inverted hexagonal phase structures of high local curvature [
42], a process believed to take place during the formation of the stalk/hemifused fusion intermediate [
43].
Unlike PE, the inclusion of the cone-shaped PA lipid at either 10 or 30 mol% in the liposome membrane did not activate fusion by HR1 in comparison to pure PC liposomes; instead, it led to a slight concentration-dependent inhibition (
Figure 4a,b). This can be explained by the negative charge of the PA headgroup, which hinders the binding of HR1 – carrying a net negative charge at pH 7.4 – with the membrane. Similarly, the reduction of HR1-mediated fusion when liposomes contain CL lipids can be attributed to the electrostatic repulsion between HR1 and the negatively charged CL lipids (
Figure 3a,b). Furthermore, the presence of PA in the liposome membrane did not increase the number of hemifusion events by HR1 compared to pure PC liposomes (
Figure 5b). This effect is likely due to electrostatic repulsions between PA headgroups, preventing clustering of PA lipids into domains of high local curvature.
In the presence of Ca
2+, HR1-mediated lipid mixing of liposomes containing 30 mol% PA was strongly activated. Specifically, we observed a 3-fold increase compared to liposomes with 30 mol% PA in the absence of Ca
2+ and a 2-fold increase compared to pure PC liposomes in the presence of Ca
2+ (
Figure 4a,b). This activation can be attributed to the formation of PA-rich membrane domains facilitated by the presence of Ca
2+, which allows the attraction of PA headgroups [
44]. At the boundary of these membrane domains, the hydrophobic chains of lipids are exposed, enabling the interaction of the amphipathic helix of HR1 with the membrane structure. This interaction can be further reinforced by electrostatic attractions between the negatively charged PA headgroup and clusters of positive residues within the HR1 sequence. The amphipathic helix of HR1 notably includes a cluster of 3 consecutive positively charged Lysine residues (395-KKK-397). Similar recognition motifs have been identified in PA-binding proteins, including those with an amphipathic helix [
45,
46,
47]. As observed in the absence of Ca
2+, hemifusion events were not promoted for liposomes containing 30 mol% PA compared to PC liposomes, suggesting that zones of high curvatures enriched in PA do not form, even when PA-PA headgroup attraction is facilitated by the presence of Ca
2+. Surprisingly, in the presence of Mg
2+ ions, HR1-mediated fusion of liposomes with 30 mol% PA was not activated. This was unexpected because Ca
2+ and Mg
2+ ions were found to have a similar capacity to induce protein-free fusion of pure PA liposomes [
48,
49]. However, Mg
2+ ions are known to be more hydrated than Ca
2+ ions [
50], which may limit their accessibility to PA headgroups that are buried deep below the PC headgroups in mixed PC:PA bilayers [
45,
51]. Such an accessibility issue would not occur in protein-free fusion of pure PA bilayers. It would also explain why Ca
2+ and Mg
2+ have comparable effects on HR1-mediated fusion of PC and PC:PE bilayers in our study (
Figure S4), as PC and PE headgroups are fully accessible. In these cases, the activation of fusion can be explained by the presence of a small fraction of negative charges in PC bilayers [
52], which may be sufficient to induce their interaction with divalent cations and promote the formation of packing defects facilitating HR1-membrane interaction.
In conclusion, we found that the interplay of divalent cations and specific cone-shaped lipids allows the formation of membrane regions with molecular packing defects, which, in turn, facilitates membrane perturbation and fusion mediated by the amphipathic helix of HR1. Future studies will need to address the role of lipids and divalent cations in the mode of action of the other functional domains of Mitofusin and their impact within the context of the full-length protein. It will also be important to consider the contribution of ER-mitochondria contact sites in regulating mitochondrial fusion events [
53,
54]. These sites, known for facilitating lipid and Ca
2+ exchange between the two organelles [
55], may act as hotspots, leading to local increases in the concentration of fusogenic lipids and Ca
2+ ions, thus promoting efficient Mitofusin-mediated fusion.
Author Contributions
Conceptualization, D.T. and M.C.; methodology, D.T.; validation, A.V., K.N. and D.T.; formal analysis, A.V. and D.T.; investigation, A.V., K.N., H.F., J.-M.G. and D.T.; resources, D.T. and M.C.; writing—original draft preparation, A.V. and D.T.; writing—review and editing, A.V., M.C. and D.T.; visualization, A.V., K.N. and D.T.; supervision, D.T.; funding acquisition, D.T. and M.C. All authors have read and agreed to the published version of the manuscript.