A Mutation in Mouse MT-ATP6 Gene Induces Respiration Defects and Opposed Effects on the Cell Tumorigenic Phenotype

Abstract

As the last step of the OXPHOS system, mitochondrial ATP synthase (or complex V) is responsible for ATP production by using the generated proton gradient, but also has an impact on other important functions linked to this system. Mutations either in complex V structural subunits, especially in mtDNA-encoded ATP6 gene, or in its assembly factors, are the molecular cause of a wide variety of human diseases, most of them classified as neurodegenerative disorders. The role of ATP synthase alterations in cancer development or metastasis has also been postulated. In this work, we reported the generation and characterization of the first mt-Atp6 pathological mutation in mouse cells, an m.8414A>G transition that promotes an amino acid change from Asn to Ser at a highly conserved residue of the protein (p.N163S), located near the path followed by protons from the intermembrane space to the mitochondrial matrix. The phenotypic consequences of the p.N163S change reproduce the effects of MT-ATP6 mutations in human diseases, such as dependence on glycolysis, defective OXPHOS activity, ATP synthesis impairment, increased ROS generation or mitochondrial membrane potential alteration. These observations demonstrate that this mutant cell line could be of great interest for the generation of mouse models with the aim of studying human diseases caused by alterations in ATP synthase. On the other hand, mutant cells showed lower migration capacity, higher expression of MHC-I and slightly lower levels of HIF-1α, indicating a possible reduction of their tumorigenic potential. These results could suggest a protective role of ATP synthase inhibition against tumor transformation that could open the door to new therapeutic strategies in those cancer types relying on OXPHOS metabolism.

1. Introduction

Although ATP synthesis, using the electrochemical proton gradient generated in the mitochondrial electron transport chain, is the main function of F₁F_o ATP synthase (EC 7.1.2.2 or complex V) [1], this enzyme is also involved in other pathways that have an impact on relevant functions connected to the OXPHOS system [2]. For instance, ATP synthase has an important role in shaping the structure of mitochondrial cristae together with other inner mitochondrial membrane (IMM) proteins, such as OPA1 (optic atrophy 1) and some components of the MICOS (Mitochondrial contact site and cristae organizing system) complex [3,4,5]; it is involved in the permeabilization of the IMM under physiological and pathological conditions [6,7] and participates in the control of intracellular signaling pathways [8], regulating cell death and survival [9,10].

Mitochondrial ATP synthase is a rotary enzyme composed of 17 different protein subunits with two of them, subunits a and A6L, encoded within mtDNA [11]. This large complex is organized in two functional domains: the hydrophobic one (F_o), embedded in the IMM, and the catalytical domain (F₁), located in the matrix and responsible for ADP phosphorylation to produce ATP. Both domains are connected by two stalks that couple proton transfer through F_o from the intermembrane space (IMS) to the matrix, to ATP synthesis, in F₁ [1,12].

Due to the crucial role of ATP synthase for maintaining cellular functions and proliferation, mutations in either complex V structural subunits or in its assembly factors are the molecular cause of a wide variety of human diseases, mostly neurodegenerative disorders [13]. The majority of the disorders affecting mitochondrial complex V characterized to date are due to mutations in the mitochondrial gene encoding subunit a MT-ATP6. This is a transmembrane protein located in the F_o domain, where it participates in two connected functions: first, it takes part of the channel that conducts protons through the IMM to the matrix and second, it interacts with the ring of c-subunits promoting its rotation driven by proton pumping [1]. The resulting mechanical energy, conducted through the stalk, induces conformational changes in F₁ domain to allow ATP production [1,12,14,15]. Moreover, this subunit has also been suggested to have an important role in complex V dimerization and stabilization. This process takes place via the F_o domain, and subunit a constitutes the most important basis for dimerization as it presents a high number of transmembrane helixes [1,16,17]. The most frequent mutations in the MT-ATP6 gene affect the same amino acid residue of subunit a, p.L156R (m.8993T>G) and p.L156P (m.8993T>C) and have been identified in numerous patients who present the NARP (neuropathy, ataxia, and retinitis pigmentosa) or MILS (maternally inherited Leigh syndrome) syndromes, depending on the mutation load. The m.8993T>G mutation compromises severely ATP production (up to 90%), due probably to a block in proton translocation or to defects in coupling or in assembly of ATP synthase [15,18,19,20,21,22,23,24,25]. In contrast, the second mutation (m.8993T>C) promotes a less severe pathological phenotype with a decrease in ATP production of 70%, due mainly to a reduction of assembly or stability of subunit a [22,26,27,28,29]. Two mutations at amino acid 217, p.L217R (m.9176T>G) and p.L217P (m.9176T>C) are also frequent and promote similar diseases to the ones associated to position 156 changes. Whereas p.L217R is extremely severe as it blocks the assembly of subunit a, the second one affects ATP synthase function without compromising the assembly [30]. Finally, the m.9185T>C (p.L220P) mutation has been found in patients with MILS and other neurological disorders [31,32] and promotes a drop of 50–90% in ATP production, due to changes in the enzyme activity. Moreover, other less frequent mutations have been described in MT-ATP6 up to date, promoting a wide variety of clinical symptoms including Leber hereditary optic neuropathy (LHON) (for a review see [30]).

Recently, ATP synthase functions have been associated in the literature with different properties of cancer cells such as proliferation behavior, and cell death evasion (reviewed in [33]). Several mutations within the MT-ATP6 gene have been found in different human carcinomas, promoting tumor growth by reducing cell death [34,35]. In addition, changes in H⁺-ATP synthase subunit expression in cancer cells have been reported. For example, a downregulation of the catalytic subunit β-F₁-ATPase expression has been found in the most prevalent human carcinomas when compared to normal tissues [36,37]. Furthermore, several tumor types such as colon, lung, breast and ovarian carcinomas, among others [38], show an upregulation of the expression of the ATPase inhibitor IF1 in order to favor cell proliferation by preserving the ATP pool [38,39]. For these reasons, mitochondrial ATP synthase has been proposed as a possible therapeutic target against cancer progression and metastasis [40,41].

Here, we reported the generation and characterization of the first mt-Atp6 pathological mutation in mouse cells, an m.8414A>G transition that promotes an amino acid change from Asn to Ser. This change affects a highly conserved residue of the protein (p.N163S), located in the aH5 helix near the path of protons from the IMS to the mitochondrial matrix and in a position laying between two pathological mutations previously reported in humans [42,43]. The phenotypic consequences of the p.N163S substitution reproduce many of the effects of MT-ATP6 mutations in human diseases, such as increased dependence on glycolysis, defective OXPHOS activity, ATP synthesis impairment, increased reactive oxygen species (ROS) generation, mitochondrial membrane potential alteration and F₁-subcomplex accumulation. These observations demonstrate that this mutant cell line could be a useful tool to generate mouse models to study human diseases caused by alterations in ATP synthase. Modelling of the mutation effects on subunit a structure and interactions suggests an alteration in proton flow and in F_o rotation that could explain the ATP synthesis defect. On the other hand, the analysis of mutant cells revealed a lower migration capacity, higher expression of major histocompatibility complex I (MHC-I) and lower levels of hypoxia inducible factor 1α (HIF-1α), suggesting a reduced tumorigenic potential. These preliminary results might indicate a protective role of ATP synthase inhibition against tumor transformation that could be explored to find new therapeutic strategies in some types of cancer with an OXPHOS-dependent metabolism.

2. Results

2.1. Isolation of a Mouse Cell Line with an mtDNA Driven Mutation in ATP Synthase

A mutagenesis and selection strategy developed to generate mtDNA mutations [44] was applied to obtain new mouse cell lines that could represent models for human mtDNA-associated mutations. Briefly, a control cell line, TmBalb/cJ, obtained by transference of mitochondria from Balb/cJ mouse platelets to mt-DNA depleted cells (ρ⁰L929^neo) as previously described [45], was treated with trimethyl-psoralen (TMP) in order to randomly induce mutations in mtDNA [44,46]. After treatment with ethidium bromide (EthBr) to reduce the mtDNA copy number, subcloning and recovery to allow mtDNA repopulation, a selection in a galactose-containing medium was applied and a potential OXPHOS defective clone (BcTMPII.1) was selected for further analysis. Full mtDNA sequencing of the potentially mutant clone showed the presence of a single mutation, an A to G transition at position 8414 (Figure 1A) within the mt-Atp6 gene in homoplasmic form, as confirmed by RFLP analysis (Figure 1B). This mutation promotes an amino acid change from asparagine to serine at position 163, which is a highly conserved exposed and functional residue located in the aH5 helix of the F_o ATPase a subunit (MT-ATP6) (Figure 1C), forming part of a positively charged region around the highly conserved Arg159, essential for the transfer of protons from the IMS to the mitochondrial matrix [30] (Figure 1D). Multiple protein sequence alignment against 100 non-redundant UniProtKB/SwissProt sequences, including animals (vertebrates and invertebrates), plants, fungi and bacteria, showed 100% conservation in this position among the 99 species, presenting more than 165 amino acids in subunit a. As shown in Table S1, both mt-Atp6 gene and protein sequences present variable identity scores among different organisms. In addition, different amino acid substitutions in MT-ATP6 are responsible for human pathologies such as NARP/MILS, LHON or Leigh diseases (Figure 1E) among others, several of them being located adjacent or very close to the equivalent of mouse 163 position mutated in the clone BcTMPII.1 and affecting aH5 helix.

Although the mutagenesis and selection methods applied mainly lead to mtDNA mutations, they could also affect nuclear genes. To exclude the effects of potential nuclear gene mutations, we performed mitochondrial transfer from both control and BcTMPII.1 cell lines into a receptor cell line lacking mtDNA (ρ⁰L929^hygro), thus generating transmitochondrial clones harboring the same nuclear background, named TmControl and TmA6^MUT, respectively, and those were used for further analysis.

2.2. Cells Carrying the m.8414A>G mtDNA Variant Are Defective in OXPHOS Performance and Dependent on Glycolysis

Cells with mtDNA mutations usually present variable difficulties to grow in a culture medium where glucose is replaced by galactose as energy source, depending on the degree of the OXPHOS impairment. As shown in Figure 2A, TmA6^MUT cells had a parallel growth curve to control cells in glucose medium but were unable to grow in the presence of galactose, whereas the control cells showed similar doubling times in both glucose and galactose media. This result strongly supports the idea that mutant cells have a severe defect in OXPHOS and rely mainly on glycolysis for energy supply. The observation of a higher glucose consumption and lactate production in the glucose-containing culture medium by TmA6^MUT cells after an identical time in culture and equivalent cell density (Figure 2B; also confirmed by the faster acidification of the medium (Figure S1A)), reinforces this interpretation.

Since ATP synthase activity is coupled to mitochondrial respiration, we measured oxygen consumption rate in intact control and mutant cells. As shown in Figure 2C, mutant cells present significantly impaired coupled mitochondrial respiration compared to control cells (left panel). Uncoupled respiration (Figure 2C, right panel), which is an index of maximal respiratory capacity independent of ATPase activity [47], confirmed this impairment. To determine which complexes were limiting respiration capacity in mutant cells, we performed polarographic measurements of respiratory chain complexes in digitonin-permeabilized cells. As illustrated in Figure 2D, the reduction in respiratory capacity in the mutant line is maintained at any entry point of electrons, suggesting that the whole electron transport chain is affected. However, the analysis of in-gel activity assay (IGA) for complex IV, that can reflect the isolated activity of this complex, did not reveal significant differences between control and mutant cells (Figure 2E).

The previous results strongly pointed to a higher dependence of the mutant cells on glycolysis due to damage in the mitochondrial function. To further test this possibility, we compared the sensitivity of control and mutant cells to drugs that inhibit glycolysis, testing different concentrations and determining their effect on cell viability. Thus, we estimated the inhibitory concentration 50 (IC₅₀) for sodium iodoacetate (IA), which irreversibly inhibits the glycolytic enzyme glyceraldehyde-3-posphate dehydrogenase (GAPDH) [48], and sodium dichloroacetate (DCA), which inhibits pyruvate dehydrogenase kinase, promoting pyruvate conversion into acetyl-CoA. As an additional control, we evaluated the cytotoxic effect of both inhibitors in ρ⁰ cells, that are completely dependent on glycolysis. As shown in Figure 2F, TmA6^MUT cells are significantly more sensitive to IA than control cells confirming their higher dependence on glucose metabolism and they are also more sensitive to DCA, confirming that their damaged OXPHOS system cannot fully compensate the reduced glycolytic route.

2.3. Mutant Cells Show Impaired ATP Synthesis and Differences in Mitochondrial Membrane Potential

In order to investigate if the mutation in mt-Atp6 gene promotes impairment in ATP synthase activity, we evaluated the mitochondrial ATP synthesis rate of mutant compared to control cells. ATP synthesis was assayed in digitonin-permeabilized TmControl and TmA6^MUT cells using pyruvate and malate as substrates. As it can be seen in Figure 3A and Figure S1B, ATP synthesis is severely hampered in TmA6^MUT cells, being produced at a rate of 11% relative to control cells (Figure 3A). The treatment with the specific inhibitor of the ATP synthase F_o fraction oligomycin, revealed that mutant cells were significantly more sensitive (lower IC₅₀) than control cells (Figure 3B). This result could be explained by the fact that their OXPHOS system is already partially damaged and hence its further inhibition causes that a threshold affecting viability (probably through essential routes connected with this activity such as nucleotide synthesis or redox balance) is reached at a lower oligomycin concentration than in cells with a stronger OXPHOS [49].

Mitochondrial membrane potential (ΔΨm) alterations can be caused by the accumulation of protons in the IMS or by an abnormal flux to the matrix coupled to ATP synthesis. Since TmA6^MUT cells maintained a substantial respiration capacity (although reduced when compared with control cells) and ATP synthesis was severely hampered, one could expect alterations in their ΔΨm. To test this possibility, we measured ΔΨm by cytometry, using DiOC₆ staining. As shown in Figure 3C, we observed a slight but not significant increase of ΔΨm in mutant cells. To further analyze this aspect, we cultivated both cell lines in the presence of different concentrations of DCA to force OXPHOS and we measured ΔΨm again. As shown in Figure 3D, although the mean fluorescence intensity (MFI) in mutants seems to decrease after treatment with 25 mM DCA, two different peaks can be distinguished: one of them matches with the size and complexity of normal-sized healthy cells, while the other one corresponds with shrunk cells, which are thus significantly affected by DCA treatment and probably dying. Considering exclusively the peak of healthy cells, the ΔΨm increases significantly when compared to that of mutant untreated cells or control cells grown in the presence of 25mM DCA (MFI: 513 vs. 345 and 484 respectively, Figure 3E). Therefore, this hyperpolarization observed when OXPHOS is forced and glycolysis inhibited, indicates an abnormal proton flux through the ATP F_o portion in the mutant cells. The measurement of ΔΨm by cytometry with a different dye, TMRE, and in the presence of an uncoupler, FCCP, confirmed the validity of the previous data and the existence of differences between both cell lines (Figure S1C).

2.4. m.8414A>G Mutation Induces Differences in Mitochondrial Mass and ROS Production

Due to defects in F₁F_o ATP synthase function, human cybrids harboring mt-ATP6 mutations have been reported to produce higher basal levels of mitochondrial ROS (mtROS) than their control counterparts [1,15,50,51]. For this reason, we sought to measure mitochondrial superoxide generation in our cell lines using MitoSOX^TM. As shown in Figure 4A, mutant cells generate around 30% more mitochondrial superoxide than control cells. To determine whether the increase in ROS could trigger a rise in the number of mitochondria or on mtDNA levels as an adaptive response, as has been reported in some cases [52,53], we evaluated mitochondrial mass. However, as shown in Figure 4B, both mtDNA copy number, measured by q-PCR, and mitochondrial mass, analyzed by MitoTracker Green staining followed by flow cytometry, indicated a significant reduction in mitochondrial mass in mutant cells. When normalizing ROS levels to mitochondrial mass, the differences between control and A6 mutant cells become much higher (Figure 4C). If the respiration values shown in Figure 2C,D are also normalized by mitochondrial content to better evaluate the respiratory chain efficiency, a significant reduction is still detected in mutant cells with respect to control cells in basal (coupled) respiration (Figure 4D). However, normalization of the respiration observed in permeabilized cells with complex-specific substrates, shows differences between mutant and control cells that become non-significant (Figure 4E), suggesting that the main defect induced by the mutation involves primary CV and not the respiratory chain efficiency, although a general reduction in respiration activity is developed in mutant cells.

2.5. m.8414A>G Mutation Does Not Affect ATP Synthase Assembly

Some described human ATP6 mutations affect only CV function [19,30,54], while others affect both assembly and function [12,23,47,54,55]. To understand the possible effect of mt-Atp6 mutation on complex V assembly, BN-PAGE followed by WB was performed, using digitonin-solubilized mitochondria from mutant and control cells. The F₁F_o band was assigned by both CV-IGA assays and Western blot (WB). As shown in Figure 5A,B, the complete monomeric F₁F_o fraction was found in similar amount in control and mutant samples. However, a lower molecular weight band with ATPase activity which could correspond to catalytically active F₁ subcomplexes (CV*, Figure 5A,B), previously reported in human mutants [12,23,47,54,55], was in proportion more intense in the mutant mitochondria. These observations could indicate a slightly more rapid turnover or slower assembly of CV in the mutant cell line. As seen in Figure 2E and Figure 5, the assembly and activity of other OXPHOS complexes was similar in both cell types.

2.6. m.8414A>G Mutation Decreases Migration Ability and Tumorigenic Phenotype in TmA6^MUT Cells

According to the literature, ATP synthase activity has been associated with proliferation, invasive capacity, and apoptosis resistance in cancer cells [33,38,40]. Indeed, most cancer cells rely on aerobic glycolysis and many of them overexpress the ATPase inhibitor factor 1 (IF1), which affects cancer development through different effects such as cristae shaping or inhibition of ATP hydrolysis, among others [56]. In addition, several mt-ATP6 mutations have been detected in different human cancer cell lines and solid tumors [57,58,59,60], next to the 8414 position and affecting aH5 helix (Figure 6A). All these aspects, together with the higher glycolysis rate, ROS production and sensitivity to DCA (a metabolic drug being tested for cancer treatment in humans) observed in TmA6^MUT cells, prompted us to evaluate aspects related to the potential tumorigenic phenotype of this cell line. First, we measured MHC-I surface expression in A6 mutant cells, as some tumor cells lose this factor to escape from immune system surveillance [61]. In this case, TmA6^MUT expression of MHC-I is increased two-fold compared to that of control cells (Figure 6B) and it is further increased by DCA treatment (Figure 6D). To test the effect of DCA on mutant cells, we measured cell viability and apoptosis after incubation with either 5 or 25 mM DCA for 72 h. As shown in Figure 6C, cell growth of mutant cells is inhibited 20% and 50% upon 5 or 25 mM DCA treatment, respectively, when compared to non-treated cells. In contrast, control cells’ growth rate is not affected by the treatment (Figure 6C, left panel). However, although DCA does not induce significant differences in cell death in either control or mutant cells at the lowest concentration tested, it promotes apoptosis in both cell lines when used at 25 mM, being significantly higher in complex V mutant cells (Figure 6C right panel). Re-evaluation of mt-ROS levels after incubation of 25 mM DCA for 72 h indicates a three-fold increase in superoxide levels of mutant cells vs the 1.4-fold increase observed in TmControl cells (Figure 6E). Finally, we analyzed the levels of the hypoxia inducible factor 1α (HIF-1α), a transcription factor that regulates the expression of genes implicated in glycolysis and plays a key role in cancer metabolism reprogramming [62], in both control and mutant cells. The analysis was performed by WB (not shown) and immunocytochemistry followed by flow cytometry in basal conditions and after DCA treatment (Figure 6F,G). As shown in Figure 6F, expression of HIF-1α under basal conditions in mutant cells, although statistically non-significant, is slightly lower than that observed in control cells, and DCA treatment did not differentially modify the expression levels of this protein among the two cell lines (Figure 6G).

To further analyse the tumor phenotype of TmA6^MUT cells, we evaluated the migration capacity of control and mutant cells using microfluidics technology. For this purpose, cells were 3D-cultured in collagen-based hydrogels within the microfluidic devices and time-lapse images were taken every 20 min for 24 h, ensuring that the tracked cells were embedded within the 3D network. Then, we calculated the average (Vmean: total distance travelled vs time) and effective (Veff: end to end distance travelled vs time) speed of both cell types. As shown in Figure 6H, although there was no difference in Veff between control and mutant cells, Vmean was significantly lower in TmA6^MUT cells, indicating lower migration ability.

3. Discussion

Mutations affecting mitochondrial complex V structural subunits or assembly factors are the molecular cause of a wide variety of mitochondrial disorders in humans. Especially relevant are those changes located in the mtDNA-encoded ATP6 gene which concentrates most of the mutations in mitochondrial ATP synthase characterized up to date [30]. One of the main challenges of mitochondrial diseases is the clinical heterogeneity observed in patients which may cause difficulties in diagnosis. The other one is the treatment of patients as, in most of the cases, there is no effective cure for mitochondrial disorders and treatments are focused on managing symptoms [63]. Although there are different cellular models that are very useful for unravelling the molecular mechanisms of mitochondrial diseases due to mutations in mtDNA, animal models are a more suitable option to study these disorders in the context of a complete organism. Since these mutations have not been found in nature in mice, to generate such models it is crucial to obtain mouse cell lines harboring mtDNA mutations.

In this work, we reported the generation and characterization of the first mt-Atp6 gene pathological mutation in mouse cells, a G to A transition at position 8414 that induces an amino acid change from Asn to Ser in a highly conserved position. The substitution affects an exposed and functional residue of the protein (p.N163S) which is located in aH5 helix and appears to have an impact on the interaction between subunit a, and the c-ring, as well as on proton translocation from the mitochondrial IMS to the matrix. The phenotypic consequences of this mutation reproduce the effects of MT-ATP6 mutations in human diseases: (i) increased dependence on glycolysis demonstrated by a higher glucose consumption and lactate production, and by the increased sensibility to glycolysis inhibitors such as sodium iodoacetate; (ii) defective OXPHOS activity, shown by their incapacity to survive in non-fermentative culture media, by the reduction in their oxygen consumption rate and the higher sensibility to DCA which forces OXPHOS; (iii) ATP synthesis impairment and higher sensitivity to oligomycin; (iv) increased ROS generation and decreased mitochondrial mass; (v) mitochondrial membrane potential increase when forced to use OXPHOS; and (vi) F₁-subcomplexes accumulation. Thus, this cell line could be of great interest to develop mouse models of human diseases due to mutations in MT-ATP6 gene.

The relevance of Helix 5 of subunit a in the function of ATP-synthase is evidenced by the high number of pathological mutations that have been found affecting this domain. As aH6, aH5 is intimately associated with the c-ring that constitutes the membrane region of the enzyme’s rotor (Figure 1D right panel) [64]. aH5 is kinked due to the presence of a proline at position 153 that enables the helix to follow the curvature of the c-ring and to seal the two hydrophilic pockets that connect the interface between subunit a and the ring to the IMS and the matrix [30]. In our mutant cell line, a substitution of an Asn by Ser at position 163 is the molecular cause for the altered phenotype observed. Although such substitution does not seem to constitute a very drastic change, as both are non-charged polar amino acids and their replacement only leads to small hydrophobicity changes, it is located in a highly conserved position, flanked by two amino acids that are mutated in a very similar way (p.A162V [43] and p.I164V [42]) in LHON patients. The molecular mechanism by which these mutations induce pathology in humans is not yet known [30] but it could be similar to that of p.N163S. Thus, TmA6^MUT cells could be an appropriate model to study disorders due to alterations in this region of mitochondrial ATP synthase subunit a.

As shown in Figure 7A, the p.N163S substitution disrupts an intramolecular H-bond between the side chains of N163 and Q210 of the aH6 that contributes to stabilize and pack both helices in an optimal conformation for defining the proton pathway. Moreover, the lateral chain of Asn 163 faces E59 of the c-ring, contributing to define the subunit a/c-ring interface. In addition, the substitution of an Asn for Ser, an amino acid with a smaller side chain, locally diminishes the solvent exposition of subunit a inducing an increase in the distance with the c-ring interface that can be detrimental for the suitable fitting of both subunits and, in consequence for c-ring rotation (Figure 7A), reducing in this way the synthesis of ATP. On the other hand, this mutation locally diminishes the positive electrostatic surface potential in a highly conserved region that also includes R159, and is essential for preventing leaks of proton by electrostatic repulsion during their transfer from the IMS entry channel to the matrix exit channel (Figure 7B) [65]. Thus, the reduction of positive surface potential may promote a retention of protons in the channel and a defective release to the matrix what could hinder c-ring rotation and therefore, ATP synthesis.

In recent years, complex V functions have been corelated with hallmarks of cancer cells such as proliferation, invasive behavior, and evasion of cell death [33]. Indeed, some carcinomas present MT-ATP6 mutations [34,35] while others downregulate the expression of the catalytic subunit β-F₁-ATPase [36,37] reducing in that way OXPHOS activity to favor the metabolic switch to a Warburg phenotype that promotes growth and cell survival [68]. In addition, it has been reported that several tumor types such as colon, lung, breast, and ovarian carcinomas upregulate the expression of ATPase inhibitory factor 1 (IF1) [38] which favors cell proliferation by blocking ATP hydrolysis and, consequently, preserving the amount of ATP [38,39]. More recently, different effects of IF1 on cancer development have been proposed [56]. In this sense, IF1 has been reported to have effects on metastasis promotion and cristae shaping [68,69], as well as on inhibition of ATP hydrolysis and escape from apoptosis. In addition, IF1 seems to modulate both ROS levels and ATP synthesis [70,71] during oxidative phosphorylation. In this work, although, our mutant cells show features which are typical of many cancer cells such as a higher sensitivity to DCA treatment and increased ROS levels, we observed that the p.N163S substitution in MT-APT6 reduces several cell features related to the tumorigenic and invasive phenotype of the cells. Thus, (i) the MHC-I expression, contrary to the downregulation observed in many cancer cells, is higher than in controls both in basal conditions and after DCA treatment; (ii) the HIF-1α levels are slightly reduced in basal conditions compared to control cells; and (iii) the migration capacity is diminished, as shown by a decrease in the average speed of mutant cells measured in vivo using microfluidic devices. This apparent reduction of tumorigenic potential and, especially, of migration ability observed in mutant cells could be explained by the decrease of ATP availability due to the severe reduction in its synthesis by complex V, although other causes cannot be ruled out. Our model suggests a protective role of ATP synthesis inhibition against tumor transformation and migration which can support the new therapeutic strategies proposing mitochondrial ATP synthase as a target in the treatment of some types of cancer with an OXPHOS-dependent metabolism [40,41].

4. Materials and Methods

4.1. Cell Lines and Media

All cell lines were grown in DMEM (GibcoBRL) supplemented with 5% FBS (fetal bovine serum, Gibco BRL). MtDNA-less mouse cells ρ⁰L929^hygro were generated by long-term exposure of L929 mouse cell line with high concentrations of Ethidium Bromide (EthBr) and transfection with the hygrocassette-containing plasmid pcDNA3.1 (Invitrogen, Waltham, MA, USA), as previously described [45,72]. Transmitochondrial cells were generated by transference of mitochondria from platelets or enucleated cell lines to mtDNA depleted cells, as described elsewhere [45,72]. Afterwards, they were isolated by growing them in DMEM supplemented with 5% dialyzed FBS and 10 µg/mL of puromycin (SIGMA) or 500 µg/mL of hygromycin (Gibco BRL). BcTMPII.1 cells were derived by random mutagenesis of TmBalb/cJ cells as previously described [44,46] and harbor an A to G transition at position 8414 within the mt-Atp6 gene.

4.2. DNA Analysis

Total DNA from cell lines was extracted using standard procedures. The complete mtDNA was amplified and sequenced as described before [73]. Primers were designed using the reference sequence (NC_005089) [74].

4.3. RFLP Analysis

To confirm the presence of the mutation RFLP analysis was achieved. Thus, a 571 bp fragment was amplified by PCR with the following primers: (i) Fw, CCAAAAGGACGAACATGAACCC (positions 8113–8134) and (ii) Rev, AGGAGGGCTGAAAAGGCTCC (positions 8664–8683, reverse primer).

The G8414 mutant version creates a recognition site for NheI and produces two bands of 298 and 273 bp upon digestion with this enzyme. Fragments were analyzed by electrophoresis in 5% polyacrylamide-TBE gel.

4.4. Growth Measurements

Growth capacity in galactose containing medium was determined by plating 5 × 10⁴ cells on 12 wells test plates in 2 mL of the appropriate medium (DMEM, which contains 4.5 mg of glucose/mL supplemented with 5% FBS, or DMEM lacking glucose and containing 0.9 mg of galactose/mL, supplemented with 5% dialyzed FBS), incubating them at 37 °C for 5 days. Cell counts were performed at daily intervals [52].

4.5. Oxygen Consumption Measurements

O₂ consumption determinations in intact or in digitonin-permeabilized cells were carried out in an Oxytherm Clark-type electrode (Hansatech), as previously described [75], with small modifications [45].

4.6. Glucose and Lactate Measurements

Glucose consumption and lactate production were measured spectrophotometrically in cell culture media collected from cultures at around 80% of confluence 24 h after medium renewal.

Glucose concentration was determined using the glucose oxidase/peroxidase kit from BioSystems. Briefly, 10 µL of diluted cell culture medium was mixed with 1 mL of reagent and absorbance at 500 nm was measured after 15 min of RT incubation. Concentration was calculated using standards of known concentration. For lactate determination, the colorimetric “L-Lactate Assay Kit” from Abcam was used. Samples were prepared following the manufacturer’s instructions and absorbance at 450 nm was measured in a 96-well plate reader.

4.7. Cell Growth Assays

The effect of m.8414A>G mtDNA mutation on sodium iodoacetate, oligomycin or DCA cytostatic ability was evaluated using the MTT reduction assay according to Mosmann et al. [76]. Briefly, cells were plated into flat-bottomed 96-well microtiter plates at a density of 5 × 10³ cells per ml. Then, cells were cultured for three days to be allowed to reach exponential growth rate before drug additions. Afterwards, cells were cultured in the presence of different inhibitors concentrations for 48 h. After drug exposure, cells were fed with fresh medium and allowed to grow for 2 population doubling times. At the end of the recovery period, plates were incubated with fresh medium and MTT for 4 h in a humidified atmosphere at 37 °C, formazan crystals were solubilized with DMSO and OD at 550 nm was read. The results obtained were given as relative values to the untreated control in percentage, and inhibitory concentration 50 (IC₅₀) was determined as the drug concentration required to reduce the absorbance to half that of the control. All experiments were performed at least in triplicate.

4.8. ATP Synthesis Measurements

ATP synthesis driven by the different substrates was measured in digitonin-permeabilized cells (2 × 10⁶ cells) using a kinetic luminescence assay, as previously described [77].

4.9. Flow Cytometry Analysis

Mitochondrial superoxide production, mitochondrial mass, mitochondrial membrane potential (MMP), MHC-I surface expression and cell death were measured in basal conditions and after DCA treatment by flow cytometry after cell staining with specific fluorescent probes, using a FACScalibur flow cytometer (BD Biosciences).

For mitochondrial superoxide production, cells were incubated with MitoSOX™ (5 µM, ThermoFisher, Waltham, MA, USA) for 30 min at 37 °C.

For mitochondrial mass, cells were stained with MitoTracker Green FM (200 nM, ThermoFisher), a mitochondrial dye that stains mitochondria regardless of Ψm, following the manufacturer’s recommendations.

For differences in Ψm due to the presence of the m.8414A>G mutation or to DCA treatment, cells were incubated with DiOC₆ at 20 nM for 30 min or TMRE at 60 nM for 15 min at 37 °C prior to flow cytometry analysis. When indicated, the uncoupler FCCP (10 µM) was added and incubated together with the fluorescent probe.

For MHC-I surface expression, cells were incubated with an anti-mouse MHC-I antibody (FITC Mouse Anti-Mouse H-2Kk, BD Biosciences, East Rutherford, NJ, USA), in PBS with 5% FBS, for 30 min at 4 °C. Next, cells were washed and, its expression was determined by flow cytometry.

Finally, cell death induced by DCA treatment was analyzed after incubation with annexin-V-FITC, which binds to the phosphatidylserine exposed in the cell surface, and 7-AAD (7-aminoactinomycin) which binds DNA, indicating cell damage, in annexin binding buffer (140 mM NaCl, 2.5 mM CaCl₂, 10 mM HEPES/NaOH pH 7.4) for 15 min at room temperature.

All data from the cytometer were analyzed using the FlowJo software.

4.10. mtDNA Copy Number Measurement

mtDNA copy number was measured by qPCR with a LightCycler^® 2.0 Instrument (Roche, Basel, Switzerland) and LightCycler Fast-Start DNA MasterPLUS SYBR Green I (Roche) in a 20 µL final volume, as recommended by the manufacturer. Total cellular DNA was used as a template and was amplified with specific oligodeoxynucleotide for mt-Co2 (NC_005089) and Sdha (AK049441) which was used as a reference for nuclear DNA content, as previously described [52].

4.11. Blue Native Polyacrylamide Electrophoresis

Mitochondria were isolated from cultured cell lines according to Schägger (1995) [78], with some modifications [79]. Digitonin-solubilized mitochondrial proteins (100 µg) were separated on blue native gradient gels (3–13% acrylamide, Invitrogen^TM, Waltham, MA, USA).

4.12. Respiratory Complex and Supercomplex Analysis

After electrophoresis, gels were further processed for either western blotting or In Gel Activity IGA assays.

In order to analyze the assembly status of respiratory complexes and supercomplexes, proteins separated by BN were transferred onto PVDF membranes and probed with specific antibodies that recognize the different respiratory complexes (anti-NDUFA9, anti-SDHA, anti-Uqcrc1 and anti-COI for complexes I, II, III, IV respectively, all of them from Invitrogen, and anti-β-F₁-ATPase, for complex V, kindly gifted by Dr. J.M. Cuezva). Then, they were incubated with their corresponding secondary antibody, tagged with peroxidase (ThermoFisher Scientific, Waltham, MA, USA). Finally, protein levels were evaluated by immunoblot, using an AmershamTM Imager 600 to document the signal.

Complex I IGA assays were performed after BN-PAGE by incubating the gels at room temperature in the presence of 0.1 M Tris-HCl pH 7.4, 0.14 mM NADH and 1 mg/mL of nitro-blue tetrazolium (NBT) as previously described [80]. Complex V IGA staining was performed by incubating the gels in 34 mM Tris, 270 mM glycine, 14 mM MgSO₄, 0.2% Pb(NO₃)₂ and 8 mM ATP at pH 7.8, as described before [81]. After incubation, the gels were photographed with a digital camera using as background a white translucid screen for Complex I and Complex IV and a black screen for Complex V [81].

4.13. HIF-1α Expression Evaluation

HIF-1α expression was evaluated by flow cytometry using a specific antibody anti- HIF-1α. 10⁵ cells were seeded in a 48-well-plate and let them attach for 24 h at 37 °C and 5% CO₂. Next day, cell culture media was removed, cells were trypsinized, transferred into a flow cytometry tube and washed with PBS three times by centrifugation. Then, cells were fixed using 4% PFA for 10 min at 4 °C. After performing three washes with PBS, cells were permeabilized with 0.1% saponin and incubated for 20 min at RT. Subsequently, saponin was removed by washing three more times and primary antibody anti- HIF-1α (1:700, Novus Bio #NB100-449) was incubated for 1h at RT. Secondary anti-rabbit antibody Alexa Fluor^® 488 (1:1000, Invitrogen #A11034) was incubated for 10 min at RT after removing primary antibody and washing three times. Finally, fluorescence intensity was assessed by flow cytometry. To avoid unspecific interactions, control cells were only incubated with the secondary labeled-antibody and compared fluorescence intensities between cells stained with both primary and secondary antibodies. Data obtained was analyzed using FlowJo software.

4.14. Cell Migration Analysis

The ability of cells to migrate was evaluated using microfluidic devices. Microdevices were fabricated in PDMS following the methodology described by Shin et al. [82]. TmControl and TmA6^MUT cells were cultured within the microfluidic chip embedded in 3D-collagen hydrogels and migration was evaluated by time-lapse imaging carried out with a Nikon D-Eclipse Ti Microscope with a 10X objective, acquiring phase contrast images every 20 min for 24 h as described before [83]. Briefly, cells were mixed with collagen type I hydrogel (2.5 mg/mL BD Biosciences) to a final concentration of 2 × 10⁵ cells/mL. Then, the hydrogel dilution was pipetted into the device’s central gel chamber and confined by surface tension. Once in place, collagen gel solution was polymerized in a humidity chamber at 37 °C and 5% CO₂ for 20 min. After that, the gel was hydrated with culture medium and incubated overnight for matrix stabilization and cell adhesion. Time-lapse imaging was obtained as described above, at 37 °C, 5% CO₂ and 95% of humidity. Cell tracking was performed using a hand coded semi-automatic MATLAB script described in previous studies [83,84]. Acquired data were then analyzed in terms of mean and effective speeds. By comparing pixel intensities and using matrix convolution techniques, the software was able to find and track cell centroids, request visual correction from the user, and post-process the migration results. These trajectories were used to extract the cell mean (Vmean) and effective (Veff) velocities. Vmean and Veff definitions are shown in Figure 5H.

4.15. Bioinformatic Tools

The structural model of MT-ATP6 from mouse was obtained from AlphaFold protein structure database [66,67] and used to create N163S mutant using mutagenesis tool embedded in PyMOL Molecular Graphics System [85] Alpha fold models conservation scores were calculated by the ConSurf server [86] using solely the protein 513 sequence of interest as input. ESP for WT and N163 variant was calculated using the APBS-PDB2PQR software suite (https://www.poissonboltzmann.org/ (accessed on 17 November 2022)) at pH 7.4 and 150 mM salt concentration. The PyMOL Molecular Graphics System was used to produce the structural figures.

4.16. Statistical Analysis

Results, displayed as mean ± SEM, were statistically analyzed by Fisher’s PLSD post hoc test from ANOVA. Data were analyzed with the statistical package StatView 5.0 (SAA Institute, Mesa, AZ, USA) or with GraphPad Prism 8.1.2 software. In all cases, differences were considered statistically significant at p ≤ 0.05. The significance of the observed differences in migration assays, was obtained via Welch’s test of Vmean and Mann–Whitney test of Veff.

Author Contributions

Conceptualization: P.F.-S., R.M.-L. and J.A.E.; Data curation: P.F.-S., R.M.-L., J.M.-B. and P.F.; Funding acquisition: P.F.-S. and R.M.-L.; Investigation: R.M.-L. (cell culture, respiration and viability analysis), N.M. (cell culture, ATP synthesis, migration analysis), J.M.-B. (Cytometer measurements), R.S.-A. (WB and immunocytochemistry), P.F. (protein modelling) and P.F.-S. (Blue-Native analysis) Project administration: P.F.-S. and R.M.-L.; Supervision: P.F.-S. and R.M.-L.; Validation: P.F.-S. and R.M.-L.; Writing—original draft: R.M.-L.; review & editing: all authors. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in supplementary material.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures

Figure 1. Analysis of m.8414A>G mutation in mt-Atp6 gene and protein: (A) chromatogram showing the m.8414A>G mutation within the mt-Atp6 gene in mutant BcTMP11.1 cells; (B) RFLP analysis of the mutation in transmitochondrial mutant (TmA6MUT) and control (TmControl) cells. The presence of the mutation creates a recognition site for NheI; (C) ConSurf analysis of FoF1 ATP synthase subunit a (MT-ATP6) from Mus musculus shown as cartoon colored according to the conservation score calculated by ConSurf server (left). The protein is folded in six transmembrane α-helixes connected by loops. The N163 (represented with sticks, black arrow) is an exposed residue located in a highly conserved region of helix aH5. As shown in the aH5 sequence (right), N163 is a highly conserved amino acid; (D) structure of mitochondrial ATP synthase (left) and interaction of subunit a with the c-ring (right). Left: Cartoon representation of subunit composition of the bovine ATP synthase monomer (PDB 6zpo). The catalytic subunits α and β (upper part) are lemon yellow and pale green, respectively. The asymmetrical central stalk subunits ƴ, δ and ε are red, light blue and violet, respectively. The peripheral stalk subunits OSCP, F6, b and d are blue, violet, pink and orange, respectively. The membrane domain composed of the c-ring in contact with subunit a, and the supernumerary subunit A6L are gray, wheat and green forest, respectively. To simplify the representation, the rest of supernumerary subunits e, f, g, j and k have been omitted. Right: View of the entire c-ring and subunit a from the intermembrane space and matrix (upper and bottom panels, respectively) and the path of protons, showing subunit a and c-ring (PDB 6zpo). The side chains of the essential residues for proton transfer, E59 and R159 (c-ring and subunit a, respectively), and the mutated N163 residue are shown as CPK colored sticks with C atoms in grey for c-ring, and in wheat for the subunit a; and (E) alignment of human and murine MT-ATP6 protein sequences. Reported human pathological mutations and the mouse mutation (pink) are highlighted in the sequence with the following color code: green for LOHN causing mutations, blue for NARP, MILS or Leigh syndromes (dark blue: frequent mutations; light blue: rare mutations), purple for autism or schizophrenia-related mutations and yellow for “other mutations”. Triangles indicate mutations creating a stop codon (red: insertion; grey: deletion). Red dots indicate the presence of mutations in the human protein that introduce an amino acid change identical to that of the murine protein.

Figure 2. Analysis of mutant cells glucose dependence and OXPHOS performance: (A) representative growing curve of TmControl and TmA6MUT in two different culture media: Glucose containing medium and galactose containing medium in which only OXPHOS competent cells are able to grow; (B) glucose consumption and lactate production by TmControl and TmA6MUT cell semi-confluent cultures in 24 h. The graph represents the concentrations relationship between lactate produced and glucose consumed, measured in cell culture media collected from a 24 h culture (80% of confluence); (C) endogenous (left) and maximal (right) oxygen consumption rate measurement in intact cells (n = 9 and 13 for TmControl and TmA6MUT, respectively; p < 0.0001 for endogenous respiration and p < 0.05 for maximal respiration); (D) oxygen consumption of digitonin-permeabilized cells in the presence of specific electron donors for complex I (glutamate + malate), complex II + III (succinate + G3P) and complex IV (TMPD) (n= 2 and 3 for TmControl and TmA6MUT, respectively; p = 0.0071, p = 0.0072 and p = 0.048 for complex I + III + IV, complex II + III + IV and complex IV, respectively). (E) Pattern of mitochondrial complexes and supercomplexes analysed by BN-PAGE followed by enzymatic In Gel Activity for complex IV (left) or by coomassie staining (right); and (F) IC50 expressed as % of that of TmControl cells for the indicated inhibitors in TmA6MUT and ρ0 cells (n = 3 in all cases for IA and n = 3, 3 and 2 for TmControl, TmA6MUT and ρ0 cells, respectively, for DCA; p < 0.05 for all the pairs and inhibitors). All values are given as mean ± SD of the mean. Asterisks indicate significant differences respect to control cells, tested by ANOVA post-hoc Fisher PLSD (p < 0.05 or as indicated in each panel description).

Figure 3. Analysis of mutant cells mitochondrial ATP synthesis activity and membrane potential: (A) ATP synthesis rate measurement in control and mutant cells (n = 9 and 6 for control and mutant cells, respectively; p < 0.0001); (B) IC50 for olygomycin in TmControl and TmA6MUT cells (n = 4 and p = 0.0008); (C) analysis of mitochondrial membrane potential (ΔΨm) in control and mutant cells. Changes in ΔΨm were determined by staining with DiOC6 and analyzed by flow cytometry. As indicated in the legend, red histogram corresponds to TmControl cells, while blue histogram corresponds to mutant cells. The left panel shows a representative histogram of ΔΨm measurement in both cell lines under basal conditions (the further to the right the peak is, the higher the membrane potential) whereas the right panel represents the average ΔΨm of mutant cells relative to controls from two independent experiments; (D) ΔΨm evaluation in both cell lines after treatment with different concentrations of DCA; and (E) Mutant cells were further separated into different populations, depending on their size (normal size and shrunken, left), and their corresponding ΔΨm was compared (right). All values are given as mean ± SD of the mean. Asterisks indicate significant differences respect to control cells, tested by ANOVA post-hoc Fisher PLSD (p < 0.05 or as indicated).

Figure 4. Analysis of mitochondrial ROS (mtROS) production and mitochondrial mass: (A) mtROS production in control and mutant cells was analyzed by flow cytometry after staining with MitoSOXTM (n = 2 for both cell lines, p = 0.0208); (B) mitochondrial mass was estimated by two different ways: mtDNA copy number quantification by qPCR (left) and flow cytometry (right) after mitochondria staining with MitoTracker GreenTM, a fluorescent probe that specifically stains mitochondria regardless of ΔΨm (n = 5 for mtDNA copy number quantification in both cell lines, p < 0.001 and n = 2 for mitochondrial mass estimation using MitoTracker Green; p = 0.0002); (C) mtROS production per mitochondria. Estimation of mitochondrial superoxide production normalized to mitochondrial mass (n = 3 in both cases, p = 0.025); (D) oxygen consumption (coupled respiration) normalized per mitochondrial mass (n = 9 and 12 for control and mutant respectively, p = 0.046); and (E) oxygen consumption of digitonin-permeabilized cells with specific substrates normalized per mitochondrial mass (n = 2 and 3 for control and mutant respectively and p > 0.05 in all cases). In all cases, the normalization by mitochondrial mass was done using the data from MitoTracker GreenTM staining. All values are expressed as mean ± SD of the mean. Asterisks indicate significant differences respect to control cells, tested by ANOVA post-hoc Fisher PLSD (p < 0.05 or as indicated).

Figure 5. OXPHOS complex and supercomplex assembly analysis: (A) pattern of mitochondrial complexes and supercomplexes containing CI (left panel) and CV (right panel) analyzed by BN-PAGE followed by enzymatic IGA of the indicated complexes; and (B) immunodetection of the different OXPHOS complexes and supercomplexes in TmControl and TmA6MUT mitochondria. Antibodies: anti-NDUFA9 for complex I, anti-Uqcrc1 for complex III, anti-Cox1 for complex IV, anti-β-F1 ATPase for complex V and anti-SDHA for complex II. Notice the higher presence of a band corresponding to CV subcomplexes (CV*) and of low molecular weight antibody-reacting material in mutant mitochondria.

Figure 6. Tumor phenotype analysis: (A) MT-ATP6 mutations found in human cancers are highlighted in green and the corresponding amino acid change is indicated under the human sequence [33]. Red dots indicate that the mutation promotes a change that is identical to the murine sequence. Blue dot indicates two different changes in the same position; (B) MHC-I expression analysis in untreated TmControl and TmA6MUT cells determined by flow cytometry; (C) Cell growth (left) and cell death (right) after treatment with different DCA concentrations. Cell growth effect was measured using the MTT reduction method (n = 2 for each cell line and treatment; p = 0.0109 and 0.0108 for 5 mM and 25 mM DCA, respectively). To evaluate cell death, cells were simultaneously stained with annexin-V-FITC and 7 AAD and fluorescence was measured by flow cytometry (n= 3 for each cell line and treatment); (D,E) effect of different DCA concentrations on MHC-I expression (D) and on mitochondrial superoxide production (E) in TmControl (upper panels) and TmA6MUT (lower panels) cell lines by flow cytometry; (F,G) analysis of HIF-1α expression levels in basal conditions; (F) and after treatment with different concentrations of DCA; (G) by immunocytochemistry followed by flow cytometry. CoCl2 (500 µM) was included as a control of hypoxia induction effect; and (H) Evaluation of cell migration using microfluidic devices. Definition (left panels) and measurement (right panels) of effective (Veff) and mean (Vmean) speed for TmControl and TmA6MUT cells (n = 3 independent experiments with 3 technical replicas each; p < 0.001). Asterisks indicate statistically significant differences between control and mutant cells: * indicates p < 0.05 and *** indicates p < 0.001.

Figure 6. Tumor phenotype analysis: (A) MT-ATP6 mutations found in human cancers are highlighted in green and the corresponding amino acid change is indicated under the human sequence [33]. Red dots indicate that the mutation promotes a change that is identical to the murine sequence. Blue dot indicates two different changes in the same position; (B) MHC-I expression analysis in untreated TmControl and TmA6MUT cells determined by flow cytometry; (C) Cell growth (left) and cell death (right) after treatment with different DCA concentrations. Cell growth effect was measured using the MTT reduction method (n = 2 for each cell line and treatment; p = 0.0109 and 0.0108 for 5 mM and 25 mM DCA, respectively). To evaluate cell death, cells were simultaneously stained with annexin-V-FITC and 7 AAD and fluorescence was measured by flow cytometry (n= 3 for each cell line and treatment); (D,E) effect of different DCA concentrations on MHC-I expression (D) and on mitochondrial superoxide production (E) in TmControl (upper panels) and TmA6MUT (lower panels) cell lines by flow cytometry; (F,G) analysis of HIF-1α expression levels in basal conditions; (F) and after treatment with different concentrations of DCA; (G) by immunocytochemistry followed by flow cytometry. CoCl2 (500 µM) was included as a control of hypoxia induction effect; and (H) Evaluation of cell migration using microfluidic devices. Definition (left panels) and measurement (right panels) of effective (Veff) and mean (Vmean) speed for TmControl and TmA6MUT cells (n = 3 independent experiments with 3 technical replicas each; p < 0.001). Asterisks indicate statistically significant differences between control and mutant cells: * indicates p < 0.05 and *** indicates p < 0.001.

Figure 7. The impact of N163S mutation in the structural model of the mouse ATP synthase subunit a: (A) detail of the interface view between the c-ring and a subunit for WT (left panel) and p.N163S (right panel). Chain-N and chain-M of c-ring are displayed with cartoon representation (light and dark gray, respectively). Subunit a is displayed as surface with cartoon representation of the secondary structure elements, helices aH5 and aH6, implicated in the interaction. Relevant residues are represented as sticks and CPK colored; and (B) impact of p.N163S mutation on the electrostatic surface potential (ESP) of mouse subunit a. ESP for wild-type (left panel)) and its N163S (right panel) variant was calculated at pH 7.4 and 150 mM of salt using the APBS-PDB2PQR software suite (https://www.poissonboltzmann.org/ (accessed on 17 November 2022)) and then plotted using PyMOL. Position of N163S mutation is indicated by a black arrow. The range of change potential is −5 (red) to +5 (blue). The structural model of mouse subunit a (AF-P00848) was obtained from AlphaFold protein structure database [66,67].

Figure 7. The impact of N163S mutation in the structural model of the mouse ATP synthase subunit a: (A) detail of the interface view between the c-ring and a subunit for WT (left panel) and p.N163S (right panel). Chain-N and chain-M of c-ring are displayed with cartoon representation (light and dark gray, respectively). Subunit a is displayed as surface with cartoon representation of the secondary structure elements, helices aH5 and aH6, implicated in the interaction. Relevant residues are represented as sticks and CPK colored; and (B) impact of p.N163S mutation on the electrostatic surface potential (ESP) of mouse subunit a. ESP for wild-type (left panel)) and its N163S (right panel) variant was calculated at pH 7.4 and 150 mM of salt using the APBS-PDB2PQR software suite (https://www.poissonboltzmann.org/ (accessed on 17 November 2022)) and then plotted using PyMOL. Position of N163S mutation is indicated by a black arrow. The range of change potential is −5 (red) to +5 (blue). The structural model of mouse subunit a (AF-P00848) was obtained from AlphaFold protein structure database [66,67].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24021300/s1.

References

1. Jonckheere, A.I.; Smeitink, J.A.; Rodenburg, R.J. Mitochondrial ATP synthase: Architecture, function and pathology. J. Inherit. Metab. Dis.; 2012; 35, pp. 211-225. [DOI: https://dx.doi.org/10.1007/s10545-011-9382-9]

2. Dominguez-Zorita, S.; Romero-Carraminana, I.; Cuezva, J.M.; Esparza-Molto, P.B. The ATPase Inhibitory Factor 1 is a Tissue-Specific Physiological Regulator of the Structure and Function of Mitochondrial ATP Synthase: A Closer Look Into Neuronal Function. Front. Physiol.; 2022; 13, 868820. [DOI: https://dx.doi.org/10.3389/fphys.2022.868820]

3. Eydt, K.; Davies, K.M.; Behrendt, C.; Wittig, I.; Reichert, A.S. Cristae architecture is determined by an interplay of the MICOS complex and the F1FO ATP synthase via Mic27 and Mic10. Microb. Cell; 2017; 4, pp. 259-272. [DOI: https://dx.doi.org/10.15698/mic2017.08.585] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28845423]

4. Kuhlbrandt, W. Structure and Mechanisms of F-Type ATP Synthases. Annu. Rev. Biochem.; 2019; 88, pp. 515-549. [DOI: https://dx.doi.org/10.1146/annurev-biochem-013118-110903] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30901262]

5. Quintana-Cabrera, R.; Quirin, C.; Glytsou, C.; Corrado, M.; Urbani, A.; Pellattiero, A.; Calvo, E.; Vazquez, J.; Enriquez, J.A.; Gerle, C. et al. The cristae modulator Optic atrophy 1 requires mitochondrial ATP synthase oligomers to safeguard mitochondrial function. Nat. Commun.; 2018; 9, 3399. [DOI: https://dx.doi.org/10.1038/s41467-018-05655-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30143614]

6. Carraro, M.; Jones, K.; Sartori, G.; Schiavone, M.; Antonucci, S.; Kucharczyk, R.; di Rago, J.P.; Franchin, C.; Arrigoni, G.; Forte, M. et al. The Unique Cysteine of F-ATP Synthase OSCP Subunit Participates in Modulation of the Permeability Transition Pore. Cell Rep.; 2020; 32, 108095. [DOI: https://dx.doi.org/10.1016/j.celrep.2020.108095]

7. Mnatsakanyan, N.; Jonas, E.A. ATP synthase c-subunit ring as the channel of mitochondrial permeability transition: Regulator of metabolism in development and degeneration. J. Mol. Cell. Cardiol.; 2020; 144, pp. 109-118. [DOI: https://dx.doi.org/10.1016/j.yjmcc.2020.05.013]

8. Esparza-Molto, P.B.; Nuevo-Tapioles, C.; Cuezva, J.M. Regulation of the H⁺-ATP synthase by IF1: A role in mitohormesis. Cell. Mol. Life Sci.; 2017; 74, pp. 2151-2166. [DOI: https://dx.doi.org/10.1007/s00018-017-2462-8]

9. Martinez-Reyes, I.; Cuezva, J.M. The H⁺-ATP synthase: A gate to ROS-mediated cell death or cell survival. Biochim. Biophys. Acta; 2014; 1837, pp. 1099-1112. [DOI: https://dx.doi.org/10.1016/j.bbabio.2014.03.010]

10. Mnatsakanyan, N.; Jonas, E.A. The new role of F1Fo ATP synthase in mitochondria-mediated neurodegeneration and neuroprotection. Exp. Neurol.; 2020; 332, 113400. [DOI: https://dx.doi.org/10.1016/j.expneurol.2020.113400]

11. Anderson, S.; Bankier, A.T.; Barrell, B.G.; de Bruijn, M.H.; Coulson, A.R.; Drouin, J.; Eperon, I.C.; Nierlich, D.P.; Roe, B.A.; Sanger, F. et al. Sequence and organization of the human mitochondrial genome. Nature; 1981; 290, pp. 457-465. [DOI: https://dx.doi.org/10.1038/290457a0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/7219534]

12. Capiau, S.; Smet, J.; De Paepe, B.; Yildiz, Y.; Arslan, M.; Stevens, O.; Verschoore, M.; Stepman, H.; Seneca, S.; Vanlander, A. Clinical Heterogeneity in MT-ATP6 Pathogenic Variants: Same Genotype-Different Onset. Cells; 2022; 11, 489. [DOI: https://dx.doi.org/10.3390/cells11030489]

13. Fernandez-Vizarra, E.; Zeviani, M. Mitochondrial disorders of the OXPHOS system. FEBS Lett.; 2021; 595, pp. 1062-1106. [DOI: https://dx.doi.org/10.1002/1873-3468.13995] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33159691]

14. Kucharczyk, R.; Ezkurdia, N.; Couplan, E.; Procaccio, V.; Ackerman, S.H.; Blondel, M.; di Rago, J.P. Consequences of the pathogenic T9176C mutation of human mitochondrial DNA on yeast mitochondrial ATP synthase. Biochim. Biophys. Acta; 2010; 1797, pp. 1105-1112. [DOI: https://dx.doi.org/10.1016/j.bbabio.2009.12.022] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20056103]

15. Sgarbi, G.; Baracca, A.; Lenaz, G.; Valentino, L.M.; Carelli, V.; Solaini, G. Inefficient coupling between proton transport and ATP synthesis may be the pathogenic mechanism for NARP and Leigh syndrome resulting from the T8993G mutation in mtDNA. Biochem. J.; 2006; 395, pp. 493-500. [DOI: https://dx.doi.org/10.1042/BJ20051748]

16. Wittig, I.; Schagger, H. Structural organization of mitochondrial ATP synthase. Biochim. Biophys. Acta; 2008; 1777, pp. 592-598. [DOI: https://dx.doi.org/10.1016/j.bbabio.2008.04.027]

17. Wittig, I.; Velours, J.; Stuart, R.; Schagger, H. Characterization of domain interfaces in monomeric and dimeric ATP synthase. Mol. Cell. Proteom.; 2008; 7, pp. 995-1004. [DOI: https://dx.doi.org/10.1074/mcp.M700465-MCP200]

18. Baracca, A.; Sgarbi, G.; Mattiazzi, M.; Casalena, G.; Pagnotta, E.; Valentino, M.L.; Moggio, M.; Lenaz, G.; Carelli, V.; Solaini, G. Biochemical phenotypes associated with the mitochondrial ATP6 gene mutations at nt8993. Biochim. Biophys. Acta; 2007; 1767, pp. 913-919. [DOI: https://dx.doi.org/10.1016/j.bbabio.2007.05.005]

19. Cortes-Hernandez, P.; Vazquez-Memije, M.E.; Garcia, J.J. ATP6 homoplasmic mutations inhibit and destabilize the human F1F0-ATP synthase without preventing enzyme assembly and oligomerization. J. Biol. Chem.; 2007; 282, pp. 1051-1058. [DOI: https://dx.doi.org/10.1074/jbc.M606828200]

20. Houstĕk, J.; Klement, P.; Hermanská, J.; Houstková, H.; Hansíková, H.; Van den Bogert, C.; Zeman, J. Altered properties of mitochondrial ATP-synthase in patients with a T→G mutation in the ATPase 6 (subunit a) gene at position 8993 of mtDNA. Biochim. Biophys. Acta; 1995; 1271, pp. 349-357. [DOI: https://dx.doi.org/10.1016/0925-4439(95)00063-A]

21. Mattiazzi, M.; Vijayvergiya, C.; Gajewski, C.D.; DeVivo, D.C.; Lenaz, G.; Wiedmann, M.; Manfredi, G. The mtDNA T8993G (NARP) mutation results in an impairment of oxidative phosphorylation that can be improved by antioxidants. Hum. Mol. Genet.; 2004; 13, pp. 869-879. [DOI: https://dx.doi.org/10.1093/hmg/ddh103] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14998933]

22. Morava, E.; Rodenburg, R.J.; Hol, F.; de Vries, M.; Janssen, A.; van den Heuvel, L.; Nijtmans, L.; Smeitink, J. Clinical and biochemical characteristics in patients with a high mutant load of the mitochondrial T8993G/C mutations. Am. J. Med. Genet. A; 2006; 140, pp. 863-868. [DOI: https://dx.doi.org/10.1002/ajmg.a.31194] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16532470]

23. Nijtmans, L.G.; Henderson, N.S.; Attardi, G.; Holt, I.J. Impaired ATP synthase assembly associated with a mutation in the human ATP synthase subunit 6 gene. J. Biol. Chem.; 2001; 276, pp. 6755-6762. [DOI: https://dx.doi.org/10.1074/jbc.M008114200] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11076946]

24. Trounce, I.; Neill, S.; Wallace, D.C. Cytoplasmic transfer of the mtDNA nt 8993 T→G (ATP6) point mutation associated with Leigh syndrome into mtDNA-less cells demonstrates cosegregation with a decrease in state III respiration and ADP/O ratio. Proc. Natl. Acad. Sci. USA; 1994; 91, pp. 8334-8338. [DOI: https://dx.doi.org/10.1073/pnas.91.18.8334] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8078883]

25. Vazquez-Memije, M.E.; Shanske, S.; Santorelli, F.M.; Kranz-Eble, P.; DeVivo, D.C.; DiMauro, S. Comparative biochemical studies of ATPases in cells from patients with the T8993G or T8993C mitochondrial DNA mutations. J. Inherit. Metab. Dis.; 1998; 21, pp. 829-836. [DOI: https://dx.doi.org/10.1023/A:1005418718299]

26. Aure, K.; Dubourg, O.; Jardel, C.; Clarysse, L.; Sternberg, D.; Fournier, E.; Laforet, P.; Streichenberger, N.; Petiot, P.; Gervais-Bernard, H. et al. Episodic weakness due to mitochondrial DNA MT-ATP6/8 mutations. Neurology; 2013; 81, pp. 1810-1818. [DOI: https://dx.doi.org/10.1212/01.wnl.0000436067.43384.0b] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24153443]

27. Craig, K.; Elliott, H.R.; Keers, S.M.; Lambert, C.; Pyle, A.; Graves, T.D.; Woodward, C.; Sweeney, M.G.; Davis, M.B.; Hanna, M.G. et al. Episodic ataxia and hemiplegia caused by the 8993T->C mitochondrial DNA mutation. J. Med. Genet.; 2007; 44, pp. 797-799. [DOI: https://dx.doi.org/10.1136/jmg.2007.052902]

28. Kucharczyk, R.; Rak, M.; di Rago, J.P. Biochemical consequences in yeast of the human mitochondrial DNA 8993T>C mutation in the ATPase6 gene found in NARP/MILS patients. Biochim. Biophys. Acta; 2009; 1793, pp. 817-824. [DOI: https://dx.doi.org/10.1016/j.bbamcr.2009.02.011]

29. Vilarinho, L.; Barbot, C.; Carrozzo, R.; Calado, E.; Tessa, A.; Dionisi-Vici, C.; Guimaraes, A.; Santorelli, F.M. Clinical and molecular findings in four new patients harbouring the mtDNA 8993T>C mutation. J. Inherit. Metab. Dis.; 2001; 24, pp. 883-884. [DOI: https://dx.doi.org/10.1023/A:1013908728445]

30. Dautant, A.; Meier, T.; Hahn, A.; Tribouillard-Tanvier, D.; di Rago, J.P.; Kucharczyk, R. ATP Synthase Diseases of Mitochondrial Genetic Origin. Front. Physiol.; 2018; 9, 329. [DOI: https://dx.doi.org/10.3389/fphys.2018.00329]

31. Castagna, A.E.; Addis, J.; McInnes, R.R.; Clarke, J.T.; Ashby, P.; Blaser, S.; Robinson, B.H. Late onset Leigh syndrome and ataxia due to a T to C mutation at bp 9,185 of mitochondrial DNA. Am. J. Med. Genet. A; 2007; 143A, pp. 808-816. [DOI: https://dx.doi.org/10.1002/ajmg.a.31637] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17352390]

32. Childs, A.M.; Hutchin, T.; Pysden, K.; Highet, L.; Bamford, J.; Livingston, J.; Crow, Y.J. Variable phenotype including Leigh syndrome with a 9185T>C mutation in the MTATP6 gene. Neuropediatrics; 2007; 38, pp. 313-316. [DOI: https://dx.doi.org/10.1055/s-2008-1065355] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18461509]

33. Galber, C.; Acosta, M.J.; Minervini, G.; Giorgio, V. The role of mitochondrial ATP synthase in cancer. Biol. Chem.; 2020; 401, pp. 1199-1214. [DOI: https://dx.doi.org/10.1515/hsz-2020-0157] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32769215]

34. Petros, J.A.; Baumann, A.K.; Ruiz-Pesini, E.; Amin, M.B.; Sun, C.Q.; Hall, J.; Lim, S.; Issa, M.M.; Flanders, W.D.; Hosseini, S.H. et al. mtDNA mutations increase tumorigenicity in prostate cancer. Proc. Natl. Acad. Sci. USA; 2005; 102, pp. 719-724. [DOI: https://dx.doi.org/10.1073/pnas.0408894102]

35. Shidara, Y.; Yamagata, K.; Kanamori, T.; Nakano, K.; Kwong, J.Q.; Manfredi, G.; Oda, H.; Ohta, S. Positive contribution of pathogenic mutations in the mitochondrial genome to the promotion of cancer by prevention from apoptosis. Cancer Res.; 2005; 65, pp. 1655-1663. [DOI: https://dx.doi.org/10.1158/0008-5472.CAN-04-2012]

36. Cuezva, J.M.; Krajewska, M.; de Heredia, M.L.; Krajewski, S.; Santamaria, G.; Kim, H.; Zapata, J.M.; Marusawa, H.; Chamorro, M.; Reed, J.C. The bioenergetic signature of cancer: A marker of tumor progression. Cancer Res.; 2002; 62, pp. 6674-6681.

37. Isidoro, A.; Martinez, M.; Fernandez, P.L.; Ortega, A.D.; Santamaria, G.; Chamorro, M.; Reed, J.C.; Cuezva, J.M. Alteration of the bioenergetic phenotype of mitochondria is a hallmark of breast, gastric, lung and oesophageal cancer. Biochem. J.; 2004; 378, pp. 17-20. [DOI: https://dx.doi.org/10.1042/bj20031541]

38. Esparza-Molto, P.B.; Cuezva, J.M. The Role of Mitochondrial H⁺-ATP Synthase in Cancer. Front. Oncol.; 2018; 8, 53. [DOI: https://dx.doi.org/10.3389/fonc.2018.00053]

39. Garcia-Bermudez, J.; Cuezva, J.M. The ATPase Inhibitory Factor 1 (IF1): A master regulator of energy metabolism and of cell survival. Biochim. Biophys. Acta; 2016; 1857, pp. 1167-1182. [DOI: https://dx.doi.org/10.1016/j.bbabio.2016.02.004]

40. Fiorillo, M.; Ozsvari, B.; Sotgia, F.; Lisanti, M.P. High ATP Production Fuels Cancer Drug Resistance and Metastasis: Implications for Mitochondrial ATP Depletion Therapy. Front. Oncol.; 2021; 11, 740720. [DOI: https://dx.doi.org/10.3389/fonc.2021.740720]

41. Fiorillo, M.; Scatena, C.; Naccarato, A.G.; Sotgia, F.; Lisanti, M.P. Bedaquiline, an FDA-approved drug, inhibits mitochondrial ATP production and metastasis in vivo, by targeting the gamma subunit (ATP5F1C) of the ATP synthase. Cell Death Differ.; 2021; 28, pp. 2797-2817. [DOI: https://dx.doi.org/10.1038/s41418-021-00788-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33986463]

42. Povalko, N.; Zakharova, E.; Rudenskaia, G.; Akita, Y.; Hirata, K.; Toyojiro, M.; Koga, Y. A new sequence variant in mitochondrial DNA associated with high penetrance of Russian Leber hereditary optic neuropathy. Mitochondrion; 2005; 5, pp. 194-199. [DOI: https://dx.doi.org/10.1016/j.mito.2005.03.003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16050984]

43. Shidara, K.; Wakakura, M. Leber’s hereditary optic neuropathy with the 3434, 9011 mitochondrial DNA point mutation. Jpn. J. Ophthalmol.; 2012; 56, pp. 175-180. [DOI: https://dx.doi.org/10.1007/s10384-011-0106-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22183138]

44. Bayona-Bafaluy, M.P.; Movilla, N.; Perez-Martos, A.; Fernandez-Silva, P.; Enriquez, J.A. Functional genetic analysis of the mammalian mitochondrial DNA encoded peptides: A mutagenesis approach. Methods Mol. Biol.; 2008; 457, pp. 379-390. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19066042]

45. Acin-Perez, R.; Bayona-Bafaluy, M.P.; Bueno, M.; Machicado, C.; Fernandez-Silva, P.; Perez-Martos, A.; Montoya, J.; Lopez-Perez, M.J.; Sancho, J.; Enriquez, J.A. An intragenic suppressor in the cytochrome c oxidase I gene of mouse mitochondrial DNA. Hum. Mol. Genet.; 2003; 12, pp. 329-339. [DOI: https://dx.doi.org/10.1093/hmg/ddg021]

46. Acin-Perez, R.; Bayona-Bafaluy, M.P.; Fernandez-Silva, P.; Moreno-Loshuertos, R.; Perez-Martos, A.; Bruno, C.; Moraes, C.T.; Enriquez, J.A. Respiratory complex III is required to maintain complex I in mammalian mitochondria. Mol. Cell; 2004; 13, pp. 805-815. [DOI: https://dx.doi.org/10.1016/S1097-2765(04)00124-8]

47. D’Aurelio, M.; Vives-Bauza, C.; Davidson, M.M.; Manfredi, G. Mitochondrial DNA background modifies the bioenergetics of NARP/MILS ATP6 mutant cells. Hum. Mol. Genet.; 2010; 19, pp. 374-386. [DOI: https://dx.doi.org/10.1093/hmg/ddp503]

48. Sabri, M.I.; Ochs, S. Inhibition of glyceraldehyde-3-phosphate dehydrogenase in mammalian nerve by iodoacetic acid. J. Neurochem.; 1971; 18, pp. 1509-1514. [DOI: https://dx.doi.org/10.1111/j.1471-4159.1971.tb00013.x]

49. Manfredi, G.; Gupta, N.; Vazquez-Memije, M.E.; Sadlock, J.E.; Spinazzola, A.; De Vivo, D.C.; Schon, E.A. Oligomycin induces a decrease in the cellular content of a pathogenic mutation in the human mitochondrial ATPase 6 gene. J. Biol. Chem.; 1999; 274, pp. 9386-9391. [DOI: https://dx.doi.org/10.1074/jbc.274.14.9386]

50. Lopez-Gallardo, E.; Emperador, S.; Solano, A.; Llobet, L.; Martin-Navarro, A.; Lopez-Perez, M.J.; Briones, P.; Pineda, M.; Artuch, R.; Barraquer, E. et al. Expanding the clinical phenotypes of MT-ATP6 mutations. Hum. Mol. Genet.; 2014; 23, pp. 6191-6200. [DOI: https://dx.doi.org/10.1093/hmg/ddu339]

51. Sikorska, M.; Sandhu, J.K.; Simon, D.K.; Pathiraja, V.; Sodja, C.; Li, Y.; Ribecco-Lutkiewicz, M.; Lanthier, P.; Borowy-Borowski, H.; Upton, A. et al. Identification of ataxia-associated mtDNA mutations (m.4052T>C and m.9035T>C) and evaluation of their pathogenicity in transmitochondrial cybrids. Muscle Nerve; 2009; 40, pp. 381-394. [DOI: https://dx.doi.org/10.1002/mus.21355] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19626676]

52. Moreno-Loshuertos, R.; Acin-Perez, R.; Fernandez-Silva, P.; Movilla, N.; Perez-Martos, A.; de Cordoba, S.R.; Gallardo, M.E.; Enriquez, J.A. Differences in reactive oxygen species production explain the phenotypes associated with common mouse mitochondrial DNA variants. Nat. Genet.; 2006; 38, pp. 1261-1268. [DOI: https://dx.doi.org/10.1038/ng1897] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17013393]

53. Moreno-Loshuertos, R.; Ferrín, G.; Acín-Pérez, R.; Gallardo, M.E.; Viscomi, C.; Pérez-Martos, A.; Zeviani, M.; Fernández-Silva, P.; Enríquez, J.A. Evolution Meets Disease: Penetrance and Functional Epistasis of Mitochondrial tRNA Mutations. PLoS Genet.; 2011; 7, e1001379. [DOI: https://dx.doi.org/10.1371/journal.pgen.1001379] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21533077]

54. Bugiardini, E.; Bottani, E.; Marchet, S.; Poole, O.V.; Beninca, C.; Horga, A.; Woodward, C.; Lam, A.; Hargreaves, I.; Chalasani, A. et al. Expanding the molecular and phenotypic spectrum of truncating MT-ATP6 mutations. Neurol. Genet.; 2020; 6, e381. [DOI: https://dx.doi.org/10.1212/NXG.0000000000000381] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32042910]

55. Kenvin, S.; Torregrosa-Muñumer, R.; Reidelbach, M.; Pennonen, J.; Turkia, J.J.; Rannila, E.; Kvist, J.; Sainio, M.T.; Huber, N.; Herukka, S.K. et al. Threshold of heteroplasmic truncating MT-ATP6 mutation in reprogramming, Notch hyperactivation and motor neuron metabolism. Hum. Mol. Genet.; 2022; 31, pp. 958-974. [DOI: https://dx.doi.org/10.1093/hmg/ddab299]

56. Gatto, C.; Grandi, M.; Solaini, G.; Baracca, A.; Giorgio, V. The F1Fo-ATPase inhibitor protein IF1 in pathophysiology. Front. Physiol.; 2022; 13, 917203. [DOI: https://dx.doi.org/10.3389/fphys.2022.917203]

57. Chatterjee, A.; Mambo, E.; Sidransky, D. Mitochondrial DNA mutations in human cancer. Oncogene; 2006; 25, pp. 4663-4674. [DOI: https://dx.doi.org/10.1038/sj.onc.1209604]

58. Grzybowska-Szatkowska, L.; Slaska, B.; Rzymowska, J.; Brzozowska, A.; Florianczyk, B. Novel mitochondrial mutations in the ATP6 and ATP8 genes in patients with breast cancer. Mol. Med. Rep.; 2014; 10, pp. 1772-1778. [DOI: https://dx.doi.org/10.3892/mmr.2014.2471]

59. Jimenez-Morales, S.; Perez-Amado, C.J.; Langley, E.; Hidalgo-Miranda, A. Overview of mitochondrial germline variants and mutations in human disease: Focus on breast cancer (Review). Int. J. Oncol.; 2018; 53, pp. 923-936. [DOI: https://dx.doi.org/10.3892/ijo.2018.4468]

60. Wu, S.; Akhtari, M.; Alachkar, H. Characterization of Mutations in the Mitochondrial Encoded Electron Transport Chain Complexes in Acute Myeloid Leukemia. Sci. Rep.; 2018; 8, 13301. [DOI: https://dx.doi.org/10.1038/s41598-018-31489-0]

61. Marco-Brualla, J.; Al-Wasaby, S.; Soler, R.; Romanos, E.; Conde, B.; Justo-Mendez, R.; Enriquez, J.A.; Fernandez-Silva, P.; Martinez-Lostao, L.; Villalba, M. et al. Mutations in the ND2 Subunit of Mitochondrial Complex I Are Sufficient to Confer Increased Tumorigenic and Metastatic Potential to Cancer Cells. Cancers; 2019; 11, 1027. [DOI: https://dx.doi.org/10.3390/cancers11071027] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31330915]

62. Semenza, G.L. HIF-1: Upstream and downstream of cancer metabolism. Curr. Opin. Genet. Dev.; 2010; 20, pp. 51-56. [DOI: https://dx.doi.org/10.1016/j.gde.2009.10.009] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19942427]

63. Russell, O.M.; Gorman, G.S.; Lightowlers, R.N.; Turnbull, D.M. Mitochondrial Diseases: Hope for the Future. Cell; 2020; 181, pp. 168-188. [DOI: https://dx.doi.org/10.1016/j.cell.2020.02.051] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32220313]

64. Vinothkumar, K.R.; Montgomery, M.G.; Liu, S.; Walker, J.E. Structure of the mitochondrial ATP synthase from Pichia angusta determined by electron cryo-microscopy. Proc. Natl. Acad. Sci. USA; 2016; 113, pp. 12709-12714. [DOI: https://dx.doi.org/10.1073/pnas.1615902113] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27791192]

65. Spikes, T.E.; Montgomery, M.G.; Walker, J.E. Structure of the dimeric ATP synthase from bovine mitochondria. Proc. Natl. Acad. Sci. USA; 2020; 117, pp. 23519-23526. [DOI: https://dx.doi.org/10.1073/pnas.2013998117] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32900941]

66. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Zidek, A.; Potapenko, A. et al. Highly accurate protein structure prediction with AlphaFold. Nature; 2021; 596, pp. 583-589. [DOI: https://dx.doi.org/10.1038/s41586-021-03819-2]

67. Varadi, M.; Anyango, S.; Deshpande, M.; Nair, S.; Natassia, C.; Yordanova, G.; Yuan, D.; Stroe, O.; Wood, G.; Laydon, A. et al. AlphaFold Protein Structure Database: Massively expanding the structural coverage of protein-sequence space with high-accuracy models. Nucleic Acids Res.; 2022; 50, pp. D439-D444. [DOI: https://dx.doi.org/10.1093/nar/gkab1061]

68. Campanella, M.; Casswell, E.; Chong, S.; Farah, Z.; Wieckowski, M.R.; Abramov, A.Y.; Tinker, A.; Duchen, M.R. Regulation of mitochondrial structure and function by the F1Fo-ATPase inhibitor protein, IF1. Cell Metab.; 2008; 8, pp. 13-25. [DOI: https://dx.doi.org/10.1016/j.cmet.2008.06.001]

69. Faccenda, D.; Tan, C.H.; Seraphim, A.; Duchen, M.R.; Campanella, M. IF1 limits the apoptotic-signalling cascade by preventing mitochondrial remodelling. Cell Death Differ.; 2013; 20, pp. 686-697. [DOI: https://dx.doi.org/10.1038/cdd.2012.163]

70. Barbato, S.; Sgarbi, G.; Gorini, G.; Baracca, A.; Solaini, G. The inhibitor protein (IF1) of the F1F0-ATPase modulates human osteosarcoma cell bioenergetics. J. Biol. Chem.; 2015; 290, pp. 6338-6348. [DOI: https://dx.doi.org/10.1074/jbc.M114.631788]

71. Garcia-Bermudez, J.; Sanchez-Arago, M.; Soldevilla, B.; Del Arco, A.; Nuevo-Tapioles, C.; Cuezva, J.M. PKA Phosphorylates the ATPase Inhibitory Factor 1 and Inactivates Its Capacity to Bind and Inhibit the Mitochondrial H⁺-ATP Synthase. Cell Rep.; 2015; 12, pp. 2143-2155. [DOI: https://dx.doi.org/10.1016/j.celrep.2015.08.052] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26387949]

72. Chomyn, A.; Lai, S.T.; Shakeley, R.; Bresolin, N.; Scarlato, G.; Attardi, G. Platelet-mediated transformation of mtDNA-less human cells: Analysis of phenotypic variability among clones from normal individuals—And complementation behavior of the tRNALys mutation causing myoclonic epilepsy and ragged red fibers. Am. J. Hum. Genet.; 1994; 54, pp. 966-974. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8198140]

73. Gallardo, M.E.; Moreno-Loshuertos, R.; Lopez, C.; Casqueiro, M.; Silva, J.; Bonilla, F.; Rodriguez de Cordoba, S.; Enriquez, J.A. m.6267G>A: A recurrent mutation in the human mitochondrial DNA that reduces cytochrome c oxidase activity and is associated with tumors. Hum. Mutat.; 2006; 27, pp. 575-582. [DOI: https://dx.doi.org/10.1002/humu.20338] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16671096]

74. Bayona-Bafaluy, M.P.; Acin-Perez, R.; Mullikin, J.C.; Park, J.S.; Moreno-Loshuertos, R.; Hu, P.; Perez-Martos, A.; Fernandez-Silva, P.; Bai, Y.; Enriquez, J.A. Revisiting the mouse mitochondrial DNA sequence. Nucleic Acids Res.; 2003; 31, pp. 5349-5355. [DOI: https://dx.doi.org/10.1093/nar/gkg739]

75. Hofhaus, G.; Shakeley, R.M.; Attardi, G. Use of polarography to detect respiration defects in cell cultures. Methods Enzymol.; 1996; 264, pp. 476-483. [DOI: https://dx.doi.org/10.1016/S0076-6879(96)64043-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8965720]

76. Mosmann, T. Rapid colorimetric assay for cellular growth and survival: Application to proliferation and cytotoxicity assays. J. Immunol. Methods; 1983; 65, pp. 55-63. [DOI: https://dx.doi.org/10.1016/0022-1759(83)90303-4]

77. Vives-Bauza, C.; Yang, L.; Manfredi, G. Assay of mitochondrial ATP synthesis in animal cells and tissues. Methods Cell Biol.; 2007; 80, pp. 155-171. [DOI: https://dx.doi.org/10.1016/S0091-679X(06)80007-5]

78. Schagger, H. Native electrophoresis for isolation of mitochondrial oxidative phosphorylation protein complexes. Methods Enzymol.; 1995; 260, pp. 190-202. [DOI: https://dx.doi.org/10.1016/0076-6879(95)60137-6]

79. Acin-Perez, R.; Fernandez-Silva, P.; Peleato, M.L.; Perez-Martos, A.; Enriquez, J.A. Respiratory active mitochondrial supercomplexes. Mol. Cell; 2008; 32, pp. 529-539. [DOI: https://dx.doi.org/10.1016/j.molcel.2008.10.021]

80. Perales-Clemente, E.; Fernandez-Vizarra, E.; Acin-Perez, R.; Movilla, N.; Bayona-Bafaluy, M.P.; Moreno-Loshuertos, R.; Perez-Martos, A.; Fernandez-Silva, P.; Enriquez, J.A. Five entry points of the mitochondrially encoded subunits in mammalian complex I assembly. Mol. Cell. Biol.; 2010; 30, pp. 3038-3047. [DOI: https://dx.doi.org/10.1128/MCB.00025-10]

81. Van Coster, R.; Smet, J.; George, E.; De Meirleir, L.; Seneca, S.; Van Hove, J.; Sebire, G.; Verhelst, H.; De Bleecker, J.; Van Vlem, B. et al. Blue native polyacrylamide gel electrophoresis: A powerful tool in diagnosis of oxidative phosphorylation defects. Pediatr. Res.; 2001; 50, pp. 658-665. [DOI: https://dx.doi.org/10.1203/00006450-200111000-00020] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11641463]

82. Shin, Y.; Han, S.; Jeon, J.S.; Yamamoto, K.; Zervantonakis, I.K.; Sudo, R.; Kamm, R.D.; Chung, S. Microfluidic assay for simultaneous culture of multiple cell types on surfaces or within hydrogels. Nat. Protoc.; 2012; 7, pp. 1247-1259. [DOI: https://dx.doi.org/10.1038/nprot.2012.051] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22678430]

83. Movilla, N.; Borau, C.; Valero, C.; Garcia-Aznar, J.M. Degradation of extracellular matrix regulates osteoblast migration: A microfluidic-based study. Bone; 2018; 107, pp. 10-17. [DOI: https://dx.doi.org/10.1016/j.bone.2017.10.025]

84. Moreno-Arotzena, O.; Borau, C.; Movilla, N.; Vicente-Manzanares, M.; Garcia-Aznar, J.M. Fibroblast Migration in 3D is Controlled by Haptotaxis in a Non-muscle Myosin II-Dependent Manner. Ann. Biomed. Eng.; 2015; 43, pp. 3025-3039. [DOI: https://dx.doi.org/10.1007/s10439-015-1343-2]

85. De Lano, W.L. PyMOL: An open-source molecular graphics tool. CCP4 725 Newsl. Protein Crystallogr.; 2002; 40, pp. 82-92.

86. Ashkenazy, H.; Abadi, S.; Martz, E.; Chay, O.; Mayrose, I.; Pupko, T.; Ben-Tal, N. ConSurf 2016: An improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res.; 2016; 44, pp. W344-W350. [DOI: https://dx.doi.org/10.1093/nar/gkw408] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27166375]