Conserved cardiolipin-mitochondrial ADP/ATP carrier interactions assume distinct structural and functional roles that are clinically relevant

Abstract

The mitochondrial phospholipid cardiolipin (CL) promotes bioenergetics via oxidative phosphorylation (OXPHOS). Three tightly bound CLs are evolutionarily conserved in the ADP/ATP carrier (AAC in yeast; adenine nucleotide translocator, ANT in mammals) which resides in the inner mitochondrial membrane and exchanges ADP and ATP to enable OXPHOS. Here, we investigated the role of these buried CLs in the carrier using yeast Aac2 as a model. We introduced negatively charged mutations into each CL-binding site of Aac2 to disrupt the CL interactions via electrostatic repulsion. While all mutations disturbing the CL-protein interaction destabilized Aac2 monomeric structure, transport activity was impaired in a pocket-specific manner. Finally, we determined that a disease-associated missense mutation in one CL-binding site in ANT1 compromised its structure and transport activity, resulting in OXPHOS defects. Our findings highlight the conserved significance of CL in AAC/ANT structure and function, directly tied to specific lipid-protein interactions.

Fig. 1:. Generation and characterization of yeast Aac2 CL-binding mutants.

(A) Strategy to disrupt CL-binding sites in Aac2. Aac2 was modeled onto bovine ANT1 (PDB ID: 2C3E) using SWISS-MODEL. Negatively charged amino acids introduced (light orange) into Aac2 CL-binding motifs face toward the CL phosphate headgroup (orange). (B) Expression of WT and mutant Aac2 was detected in whole cell extracts by immunoblot; Kgd1 and Tom70 served as loading controls (n=5). Significant difference was obtained by one-way ANOVA with Dunnett’s multiple comparisons test (vs. WT) ***p<0.001. (C) Growth phenotype of Aac2 CL-binding mutants. Serial dilutions of indicated cells were spotted onto fermentable (YPD) and respiratory (YPEG) media and incubated at 30 or 37°C for 3–5 days (n=3). (D) Membrane topology of WT and mutant Aac2. Isolated mitochondria were osmotically ruptured and treated with or without proteinase K as indicated. Aac2 N-terminus was detected by a monoclonal antibody 6H8 recognizing the first 13 amino acids MSSNAQVKTPLPP (n=6). IMS, intermembrane space; IMM, inner mitochondrial membrane; OMM, outer mitochondrial membrane. In B-D, representative images from the indicated replicates are shown.

Fig. 2:. Native mass spectrometry analysis to detect lipid-protein interaction of Aac2.

(A) FlagAac2 affinity purified from WT mitochondria was associated with CATR and up to three CL molecules. (B) FlagAac2 affinity purified from crd1Δ mitochondria which lack CL did not co-purify other phospholipids. Insets in A and B show the fractional population of FlagAac2 relative to FlagAac2+CATR when immunoprecipitated from WT and crd1Δ mitochondria, respectively. Mean with SEM (n=3). Significant differences as determined by Student’s t-test indicated (****p<0.0001). (C) CL, PE, PG, or lyso-PC was added to FlagAac2 purified from crd1Δ mitochondria at the indicated concentrations. (D) Fractional population of FlagAac2, FlagAac2+CATR alone or associated with one to three CL molecules was determined for WT and Aac2 CL-binding mutants. Mean with SEM (n=5–6). Significant differences were obtained by one-way ANOVA with Dunnett’s multiple comparisons test (* vs. FlagAac2; † vs. FlagAac2+CATR+1CL); *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001.

Fig. 3:. CLs associated with Aac2 stabilize its monomeric structure.

Mitochondria (100 μg) from WT and Aac2 CL-binding mutants were mock-treated or instead incubated with either 40 μM CATR (A) or 10 μM BKA (B) and then solubilized with 1.5% (w/v) digitonin. The extracts were resolved by 6 to 16% blue native-PAGE and immunoblotted for Aac2 (n=6). Representative images from the indicated replicates are shown.

Fig. 4:. Distinct functional roles of the three tightly bound CLs.

(A-C) ADP/ATP exchange: The efflux of matrix ATP was detected with isolated mitochondria as NADPH formation (A340; absorbance at 340 nm) occurring coupled with an in vitro glycolysis reaction which contained glucose, hexokinase, and glucose-6-phosphate dehydrogenase. The reaction was initiated by adding ADP. The measurement was performed in the presence of 5 mM malate and 5 mM pyruvate (+Mal/Pyr). Where indicated, WT mitochondria were treated with 5 μM CATR prior to the efflux reaction. (A) The linear part of the initial velocity for the ATP efflux was plotted and curve fitting performed by nonlinear regression (mean with SEM, n=6). Plots of aac2Δ and WT are repeated in all panels. (B) The initial velocity following the addition of 33 μM ADP presented as scatter plots (mean with SEM, n=6). (C) Fitted Km and Vmax values were obtained using the Michaelis-Menten equation (mean). The range of determined values is shown in brackets (n.d., not detected). (D) Oxygen consumption rate (OCR) in isolated mitochondria with sequential injections of respiratory substrate (NADH or succinate), ADP, oligomycin, and CCCP was measured using Seahorse XF96e FluxAnalyzer. (E, F) ADP-stimulated respiration and respiratory control ratio (RCR) of WT and Aac2 mutant mitochondria under NADH (E) and succinate (F) were plotted. RCR was obtained by dividing OCR of ADP-stimulated respiration by that of basal respiration (see also fig. S5). WT mitochondria was treated with 50 μM CATR before the measurement. Mean with SEM, n=21–35. Significant differences obtained by two-way ANOVA followed by Tukey’s multiple comparisons test are shown as * for comparison with WT and † for comparison between pocket mutants; *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001.

Fig. 5. Disturbed Aac2-CL interaction is a pathological mechanism.

(A) Mitochondrial myopathy patient mutation L155F was introduced into yeast Aac2 as Fig. 1A. (B) Native MS analysis obtained the fractional population of WT and indicated FlagAac2 mutants associated with CATR and one to three CL molecules. (C) Expression of WT and mutant Aac2 was detected in whole cell extract by immunoblot. Kgd1 and Tom70 served as loading controls (n=4). (D) Growth phenotype of Aac2 CL-binding mutants. Serial dilutions of indicated cells were spotted onto fermentable (YPD) and respiratory (YPEG) media and incubated at 30 or 37°C for 3–5 days (n=4). (E) Aac2 monomeric structure: 100 μg of mitochondria from WT and Aac2-binding mutants were mock-treated or instead incubated with either 40 μM CATR (C) or 10 μM BKA (B) and then solubilized with 1.5% (w/v) digitonin. The extracts were resolved by 6 to 16% blue native-PAGE and immunoblotted for FlagAac2 (n=6). (F-H) ADP/ATP exchange: The efflux of matrix ATP was detected with isolated mitochondria as NADPH formation (A340; absorbance at 340 nm) occurring coupled with in vitro glycolysis reaction as Fig. 4. The measurement was performed in the presence of 5 mM malate and 5 mM pyruvate (n=6). (F) The linear part of the initial velocity following the addition of ADP at indicated concentrations was plotted (mean with SEM). Curve fitting was performed by nonlinear regression. (G) The linear part of velocity when 33 μM ADP was added shown as scatter plots (mean with SEM). (H) Fitted Km and Vmax values were obtained by the Michaelis-Menten equation from the replicated experiments (mean). Significant differences were obtained by one-way ANOVA with Dunnett’s multiple comparisons test (vs. WT); *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001. (I) Predicted roles of the buried CLs within Aac2. If CLs are dissociated, Aac2 monomeric structure is destabilized and ADP/ATP transport is compromised, which disrupts energy production via OXPHOS. In C-E, representative images from the indicated replicates are shown.

Fig. 6. Human ANT1 L141F are structurally and functionally compromised.

(A) Human ANT1 was modeled onto bovine ANT1 (PDB ID: 2C3E) using SWISS-MODEL. (B-F) FlagANT1 and the indicated ANT1 mutants were induced by 0.25 μg/ml doxycycline in ant^null T-REx-293 cells with three ANT isoforms (ANT1, 2, and 3) knocked out. (B) Expression of ANT1 was detected in whole cell extracts by immunoblot. β-actin and GRP75 were loading controls (n=3). (C) ANT1 monomeric structure: 80 μg of mitochondria from WT and ANT1 mutants were mock-treated or instead incubated with either 40 μM CATR or 10 μM BKA and then solubilized with 1.5% (w/v) digitonin. The extracts were resolved by 6 to 16% blue native-PAGE and immunoblotted for FlagANT1 (n=3). (D-E) ADP/ATP exchange: The efflux of matrix ATP was detected with isolated mitochondria as NADPH formation (A340; absorbance at 340 nm) as Fig. 4. The measurement was performed in the presence of 5 mM malate and 5 mM pyruvate. 5 μM CATR was added to WT mitochondria prior to stimulating the efflux (n=3). (D) The initial velocity following the addition of ADP at indicated concentrations was plotted (mean with SEM). Curve fitting was performed by nonlinear regression. Significant differences obtained by one-way ANOVA with Dunnett’s multiple comparisons test are shown as *, L141E; †, L141F; #, WT+CATR (vs. WT). *p<0.05, **p<0.01, ***p<0.001. (E) Fitted Km and Vmax values from the Michaelis-Menten equation (mean). (F) Cellular oxygen consumption rate (OCR) was measured using a Seahorse XF96e Flux Analyzer with the Mito Stress Test kit under indicated conditions. Basal and maximal OCR were obtained under glucose stimulation after FCCP treatment to uncouple mitochondria. ATP production-coupled respiration is defined as basal OCR subtracted by post-oligomycin OCR. Significant differences were determined by one-way ANOVA with Dunnett’s multiple comparisons test (vs. WT), ***p<0.001 ****p < 0.0001. Means with SEM (n=16). In B and C, representative images from the indicated replicates are shown.

Fig. 7:. MD simulations predict reduced CL-ANT1 affinity at pocket 2 for L141 mutants.

(A-B). Time-averaged 2D density maps of CL lipids for the equilibrium prebound (A) and unbound (B) simulations. The CL binding pockets are defined by the blue (pocket 1), cyan (pocket 2) and green (pocket 3) dots, which represent the amino acids that comprise each binding site. The density surfaces were calculated using vmd volmap with 1 Å resolution. Occupancy scale bar on right applies to all images. (C-D) Radial distribution function profile for the identification of CL lipids around residue 141 for prebound (C) and unbound (D) simulation. (E) The estimated relative free energy required to decouple CL from binding pocket 2 of ANT1 for WT and mutants L141F and L141E. Using Wilcoxson rank sum test between WT and L141E indicates statistical significance (p< 0.05), while the difference between WT and L141F is not significant. Note: Binding site residues considered in the present study are based on the protein-CL lipid interactions and are the selected amino acid residues within 8 Å of the bound CL molecules: Pocket 1: 36, 53, 54, 55, 271, 272, 273, 274, 275, and 276; Pocket 2: 71, 72, 73, 74, 75, 141, 152, 155, 156, 157, and 158; Pocket 3: 251, 252, 253, 254, 255, 174, 175, 176, 177, and 178, respectively.