Comparative Clustering (CompaCt) of eukaryote complexomes identifies novel interactions and sheds light on protein complex evolution.

Complexome profiling allows large-scale, untargeted, and comprehensive characterization of protein complexes in a biological sample using a combined approach of separating intact protein complexes e.g., by native gel electrophoresis, followed by mass spectrometric analysis of the proteins in the resulting fractions. Over the last decade, its application has resulted in a large collection of complexome profiling datasets. While computational methods have been developed for the analysis of individual datasets, methods for large-scale comparative analysis of complexomes from multiple species are lacking. Here, we present Comparative Clustering (CompaCt), that performs fully automated integrative analysis of complexome profiling data from multiple species, enabling systematic characterization and comparison of complexomes. CompaCt implements a novel method for leveraging orthology in comparative analysis to allow systematic identification of conserved as well as taxon-specific elements of the analyzed complexomes. We applied this method to a collection of 53 complexome profiles spanning the major branches of the eukaryotes. We demonstrate the ability of CompaCt to robustly identify the composition of protein complexes, and show that integrated analysis of multiple datasets improves characterization of complexes from specific complexome profiles when compared to separate analyses. We identified novel candidate interactors and complexes in a number of species from previously analyzed datasets, like the emp24, the V-ATPase and mitochondrial ATP synthase complexes. Lastly, we demonstrate the utility of CompaCt for the automated large-scale characterization of the complexome of the mosquito Anopheles stephensi shedding light on the evolution of metazoan protein complexes. CompaCt is available from https://github.com/cmbi/compact-bio.

The cytochrome b carboxyl terminal region is necessary for mitochondrial complex III assembly.

Mitochondrial bc 1 complex from yeast has 10 subunits, but only cytochrome b (Cytb) subunit is encoded in the mitochondrial genome. Cytb has eight transmembrane helices containing two hemes b for electron transfer. Cbp3 and Cbp6 assist Cytb synthesis, and together with Cbp4 induce Cytb hemylation. Subunits Qcr7/Qcr8 participate in the first steps of assembly, and lack of Qcr7 reduces Cytb synthesis through an assembly-feedback mechanism involving Cbp3/Cbp6. Because Qcr7 resides near the Cytb carboxyl region, we wondered whether this region is important for Cytb synthesis/assembly. Although deletion of the Cytb C-region did not abrogate Cytb synthesis, the assembly-feedback regulation was lost, so Cytb synthesis was normal even if Qcr7 was missing. Mutants lacking the Cytb C-terminus were non-respiratory because of the absence of fully assembled bc 1 complex. By performing complexome profiling, we showed the existence of aberrant early-stage subassemblies in the mutant. In this work, we demonstrate that the C-terminal region of Cytb is critical for regulation of Cytb synthesis and bc 1 complex assembly.

MRPS36 provides a structural link in the eukaryotic 2-oxoglutarate dehydrogenase complex.

The tricarboxylic acid cycle is the central pathway of energy production in eukaryotic cells and plays a key part in aerobic respiration throughout all kingdoms of life. One of the pivotal enzymes in this cycle is 2-oxoglutarate dehydrogenase complex (OGDHC), which generates NADH by oxidative decarboxylation of 2-oxoglutarate to succinyl-CoA. OGDHC is a megadalton protein complex originally thought to be assembled from three catalytically active subunits (E1o, E2o, E3). In fungi and animals, however, the protein MRPS36 has more recently been proposed as a putative additional component. Based on extensive cross-linking mass spectrometry data supported by phylogenetic analyses, we provide evidence that MRPS36 is an important member of the eukaryotic OGDHC, with no prokaryotic orthologues. Comparative sequence analysis and computational structure predictions reveal that, in contrast with bacteria and archaea, eukaryotic E2o does not contain the peripheral subunit-binding domain (PSBD), for which we propose that MRPS36 evolved as an E3 adaptor protein, functionally replacing the PSBD. We further provide a refined structural model of the complete eukaryotic OGDHC of approximately 3.45 MDa with novel mechanistic insights.

Regulation of mitochondrial proteostasis by the proton gradient.

Mitochondria adapt to different energetic demands reshaping their proteome. Mitochondrial proteases are emerging as key regulators of these adaptive processes. Here, we use a multiproteomic approach to demonstrate the regulation of the m-AAA protease AFG3L2 by the mitochondrial proton gradient, coupling mitochondrial protein turnover to the energetic status of mitochondria. We identify TMBIM5 (previously also known as GHITM or MICS1) as a Ca2+ /H+ exchanger in the mitochondrial inner membrane, which binds to and inhibits the m-AAA protease. TMBIM5 ensures cell survival and respiration, allowing Ca2+ efflux from mitochondria and limiting mitochondrial hyperpolarization. Persistent hyperpolarization, however, triggers degradation of TMBIM5 and activation of the m-AAA protease. The m-AAA protease broadly remodels the mitochondrial proteome and mediates the proteolytic breakdown of respiratory complex I to confine ROS production and oxidative damage in hyperpolarized mitochondria. TMBIM5 thus integrates mitochondrial Ca2+ signaling and the energetic status of mitochondria with protein turnover rates to reshape the mitochondrial proteome and adjust the cellular metabolism.

AIFM1 is a component of the mitochondrial disulfide relay that drives complex I assembly through efficient import of NDUFS5.

The mitochondrial intermembrane space protein AIFM1 has been reported to mediate the import of MIA40/CHCHD4, which forms the import receptor in the mitochondrial disulfide relay. Here, we demonstrate that AIFM1 and MIA40/CHCHD4 cooperate beyond this MIA40/CHCHD4 import. We show that AIFM1 and MIA40/CHCHD4 form a stable long-lived complex in vitro, in different cell lines, and in tissues. In HEK293 cells lacking AIFM1, levels of MIA40 are unchanged, but the protein is present in the monomeric form. Monomeric MIA40 neither efficiently interacts with nor mediates the import of specific substrates. The import defect is especially severe for NDUFS5, a subunit of complex I of the respiratory chain. As a consequence, NDUFS5 accumulates in the cytosol and undergoes rapid proteasomal degradation. Lack of mitochondrial NDUFS5 in turn results in stalling of complex I assembly. Collectively, we demonstrate that AIFM1 serves two overlapping functions: importing MIA40/CHCHD4 and constituting an integral part of the disulfide relay that ensures efficient interaction of MIA40/CHCHD4 with specific substrates.

Ablation of mitochondrial DNA results in widespread remodeling of the mitochondrial complexome.

So-called ρ0 cells lack mitochondrial DNA and are therefore incapable of aerobic ATP synthesis. How cells adapt to survive ablation of oxidative phosphorylation remains poorly understood. Complexome profiling analysis of ρ0 cells covered 1,002 mitochondrial proteins and revealed changes in abundance and organization of numerous multiprotein complexes including previously not described assemblies. Beyond multiple subassemblies of complexes that would normally contain components encoded by mitochondrial DNA, we observed widespread reorganization of the complexome. This included distinct changes in the expression pattern of adenine nucleotide carrier isoforms, other mitochondrial transporters, and components of the protein import machinery. Remarkably, ablation of mitochondrial DNA hardly affected the complexes organizing cristae junctions indicating that the altered cristae morphology in ρ0 mitochondria predominantly resulted from the loss of complex V dimers required to impose narrow curvatures to the inner membrane. Our data provide a comprehensive resource for in-depth analysis of remodeling of the mitochondrial complexome in response to respiratory deficiency.

Molecular characterization of a complex of apoptosis-inducing factor 1 with cytochrome c oxidase of the mitochondrial respiratory chain.

Combining mass spectrometry-based chemical cross-linking and complexome profiling, we analyzed the interactome of heart mitochondria. We focused on complexes of oxidative phosphorylation and found that dimeric apoptosis-inducing factor 1 (AIFM1) forms a defined complex with ∼10% of monomeric cytochrome c oxidase (COX) but hardly interacts with respiratory chain supercomplexes. Multiple AIFM1 intercross-links engaging six different COX subunits provided structural restraints to build a detailed atomic model of the COX-AIFM12 complex (PDBDEV_00000092). An application of two complementary proteomic approaches thus provided unexpected insight into the macromolecular organization of the mitochondrial complexome. Our structural model excludes direct electron transfer between AIFM1 and COX. Notably, however, the binding site of cytochrome c remains accessible, allowing formation of a ternary complex. The discovery of the previously overlooked COX-AIFM12 complex and clues provided by the structural model hint at potential roles of AIFM1 in oxidative phosphorylation biogenesis and in programmed cell death.

Composition and stage dynamics of mitochondrial complexes in Plasmodium falciparum.

Our current understanding of mitochondrial functioning is largely restricted to traditional model organisms, which only represent a fraction of eukaryotic diversity. The unusual mitochondrion of malaria parasites is a validated drug target but remains poorly understood. Here, we apply complexome profiling to map the inventory of protein complexes across the pathogenic asexual blood stages and the transmissible gametocyte stages of Plasmodium falciparum. We identify remarkably divergent composition and clade-specific additions of all respiratory chain complexes. Furthermore, we show that respiratory chain complex components and linked metabolic pathways are up to 40-fold more prevalent in gametocytes, while glycolytic enzymes are substantially reduced. Underlining this functional switch, we find that cristae are exclusively present in gametocytes. Leveraging these divergent properties and stage dynamics for drug development presents an attractive opportunity to discover novel classes of antimalarials and increase our repertoire of gametocytocidal drugs.

CEDAR, an online resource for the reporting and exploration of complexome profiling data.

Complexome profiling is an emerging 'omics' approach that systematically interrogates the composition of protein complexes (the complexome) of a sample, by combining biochemical separation of native protein complexes with mass-spectrometry based quantitation proteomics. The resulting fractionation profiles hold comprehensive information on the abundance and composition of the complexome, and have a high potential for reuse by experimental and computational researchers. However, the lack of a central resource that provides access to these data, reported with adequate descriptions and an analysis tool, has limited their reuse. Therefore, we established the ComplexomE profiling DAta Resource (CEDAR, www3.cmbi.umcn.nl/cedar/), an openly accessible database for depositing and exploring mass spectrometry data from complexome profiling studies. Compatibility and reusability of the data is ensured by a standardized data and reporting format containing the "minimum information required for a complexome profiling experiment" (MIACE). The data can be accessed through a user-friendly web interface, as well as programmatically using the REST API portal. Additionally, all complexome profiles available on CEDAR can be inspected directly on the website with the profile viewer tool that allows the detection of correlated profiles and inference of potential complexes. In conclusion, CEDAR is a unique, growing and invaluable resource for the study of protein complex composition and dynamics across biological systems.

Investigation of central energy metabolism-related protein complexes of ANME-2d methanotrophic archaea by complexome profiling.

The anaerobic oxidation of methane is important for mitigating emissions of this potent greenhouse gas to the atmosphere and is mediated by anaerobic methanotrophic archaea. In a 'Candidatus Methanoperedens BLZ2' enrichment culture used in this study, methane is oxidized to CO2 with nitrate being the terminal electron acceptor of an anaerobic respiratory chain. Energy conservation mechanisms of anaerobic methanotrophs have mostly been studied at metagenomic level and hardly any protein data is available at this point. To close this gap, we used complexome profiling to investigate the presence and subunit composition of protein complexes involved in energy conservation processes. All enzyme complexes and their subunit composition involved in reverse methanogenesis were identified. The membrane-bound enzymes of the respiratory chain, such as F420H2:quinone oxidoreductase, membrane-bound heterodisulfide reductase, nitrate reductases and Rieske cytochrome bc1 complex were all detected. Additional or putative subunits such as an octaheme subunit as part of the Rieske cytochrome bc1 complex were discovered that will be interesting targets for future studies. Furthermore, several soluble proteins were identified, which are potentially involved in oxidation of reduced ferredoxin produced during reverse methanogenesis leading to formation of small organic molecules. Taken together these findings provide an updated, refined picture of the energy metabolism of the environmentally important group of anaerobic methanotrophic archaea.

Novel defect in phosphatidylinositol 4-kinase type 2-alpha (PI4K2A) at the membrane-enzyme interface is associated with metabolic cutis laxa.

Inherited cutis laxa, or inelastic, sagging skin is a genetic condition of premature and generalised connective tissue ageing, affecting various elastic components of the extracellular matrix. Several cutis laxa syndromes are inborn errors of metabolism and lead to severe neurological symptoms. In a patient with cutis laxa, a choreoathetoid movement disorder, dysmorphic features and intellectual disability we performed exome sequencing to elucidate the underlying genetic defect. We identified the amino acid substitution R275W in phosphatidylinositol 4-kinase type IIα, caused by a homozygous missense mutation in the PI4K2A gene. We used lipidomics, complexome profiling and functional studies to measure phosphatidylinositol 4-phosphate synthesis in the patient and evaluated PI4K2A deficient mice to define a novel metabolic disorder. The R275W residue, located on the surface of the protein, is involved in forming electrostatic interactions with the membrane. The catalytic activity of PI4K2A in patient fibroblasts was severely reduced and lipid mass spectrometry showed that particular acyl-chain pools of PI4P and PI(4,5)P2 were decreased. Phosphoinositide lipids play a major role in intracellular signalling and trafficking and regulate the balance between proliferation and apoptosis. Phosphatidylinositol 4-kinases such as PI4K2A mediate the first step in the main metabolic pathway that generates PI4P, PI(4,5)P2 and PI(3,4,5)P3 . Although neurologic involvement is common, cutis laxa has not been reported previously in metabolic defects affecting signalling. Here we describe a patient with a complex neurological phenotype, premature ageing and a mutation in PI4K2A, illustrating the importance of this enzyme in the generation of inositol lipids with particular acylation characteristics.

TMEM70 functions in the assembly of complexes I and V.

Protein complexes from the oxidative phosphorylation (OXPHOS) system are assembled with the help of proteins called assembly factors. We here delineate the function of the inner mitochondrial membrane protein TMEM70, in which mutations have been linked to OXPHOS deficiencies, using a combination of BioID, complexome profiling and coevolution analyses. TMEM70 interacts with complex I and V and for both complexes the loss of TMEM70 results in the accumulation of an assembly intermediate followed by a reduction of the next assembly intermediate in the pathway. This indicates that TMEM70 has a role in the stability of membrane-bound subassemblies or in the membrane recruitment of subunits into the forming complex. Independent evidence for a role of TMEM70 in OXPHOS assembly comes from evolutionary analyses. The TMEM70/TMEM186/TMEM223 protein family, of which we show that TMEM186 and TMEM223 are mitochondrial in human as well, only occurs in species with OXPHOS complexes. Our results validate the use of combining complexome profiling with BioID and evolutionary analyses in elucidating congenital defects in protein complex assembly.

A salvage pathway maintains highly functional respiratory complex I.

Regulation of the turnover of complex I (CI), the largest mitochondrial respiratory chain complex, remains enigmatic despite huge advancement in understanding its structure and the assembly. Here, we report that the NADH-oxidizing N-module of CI is turned over at a higher rate and largely independently of the rest of the complex by mitochondrial matrix protease ClpXP, which selectively removes and degrades damaged subunits. The observed mechanism seems to be a safeguard against the accumulation of dysfunctional CI arising from the inactivation of the N-module subunits due to attrition caused by its constant activity under physiological conditions. This CI salvage pathway maintains highly functional CI through a favorable mechanism that demands much lower energetic cost than de novo synthesis and reassembly of the entire CI. Our results also identify ClpXP activity as an unforeseen target for therapeutic interventions in the large group of mitochondrial diseases characterized by the CI instability.

Complexome analysis of the nitrite-dependent methanotroph Methylomirabilis lanthanidiphila.

The atmospheric concentration of the potent greenhouse gases methane and nitrous oxide (N2O) has increased drastically during the last century. Methylomirabilis bacteria can play an important role in controlling the emission of these two gases from natural ecosystems, by oxidizing methane to CO2 and reducing nitrite to N2 without producing N2O. These bacteria have an anaerobic metabolism, but are proposed to possess an oxygen-dependent pathway for methane activation. Methylomirabilis bacteria reduce nitrite to NO, and are proposed to dismutate NO into O2 and N2 by a putative NO dismutase (NO-D). The O2 produced in the cell can then be used to activate methane by a particulate methane monooxygenase. So far, the metabolic model of Methylomirabilis bacteria was based mainly on (meta)genomics and physiological experiments. Here we applied a complexome profiling approach to determine which of the proposed enzymes are actually expressed in Methylomirabilis lanthanidiphila. To validate the proposed metabolic model, we focused on enzymes involved in respiration, as well as nitrogen and carbon transformation. All complexes suggested to be involved in nitrite-dependent methane oxidation, were identified in M. lanthanidiphila, including the putative NO-D. Furthermore, several complexes involved in nitrate reduction/nitrite oxidation and NO reduction were detected, which likely play a role in detoxification and redox homeostasis. In conclusion, complexome profiling validated the expression and composition of enzymes hypothesized to be involved in the energy, methane and nitrogen metabolism of M. lanthanidiphila, thereby further corroborating their unique metabolism involved in the environmentally relevant process of nitrite-dependent methane oxidation.

COmplexome Profiling ALignment (COPAL) reveals remodeling of mitochondrial protein complexes in Barth syndrome.

MOTIVATION: Complexome profiling combines native gel electrophoresis with mass spectrometry to obtain the inventory, composition and abundance of multiprotein assemblies in an organelle. Applying complexome profiling to determine the effect of a mutation on protein complexes requires separating technical and biological variations from the variations caused by that mutation. RESULTS: We have developed the COmplexome Profiling ALignment (COPAL) tool that aligns multiple complexome profiles with each other. It includes the abundance profiles of all proteins on two gels, using a multi-dimensional implementation of the dynamic time warping algorithm to align the gels. Subsequent progressive alignment allows us to align multiple profiles with each other. We tested COPAL on complexome profiles from control mitochondria and from Barth syndrome (BTHS) mitochondria, which have a mutation in tafazzin gene that is involved in remodeling the inner mitochondrial membrane phospholipid cardiolipin. By comparing the variation between BTHS mitochondria and controls with the variation among either, we assessed the effects of BTHS on the abundance profiles of individual proteins. Combining those profiles with gene set enrichment analysis allows detecting significantly affected protein complexes. Most of the significantly affected protein complexes are located in the inner mitochondrial membrane (mitochondrial contact site and cristae organizing system, prohibitins), or are attached to it (the large ribosomal subunit). AVAILABILITY AND IMPLEMENTATION: COPAL is written in python and is available from http://github.com/cmbi/copal. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cryo-EM structure of respiratory complex I at work.

Mitochondrial complex I has a key role in cellular energy metabolism, generating a major portion of the proton motive force that drives aerobic ATP synthesis. The hydrophilic arm of the L-shaped ~1 MDa membrane protein complex transfers electrons from NADH to ubiquinone, providing the energy to drive proton pumping at distant sites in the membrane arm. The critical steps of energy conversion are associated with the redox chemistry of ubiquinone. We report the cryo-EM structure of complete mitochondrial complex I from the aerobic yeast Yarrowia lipolytica both in the deactive form and after capturing the enzyme during steady-state activity. The site of ubiquinone binding observed during turnover supports a two-state stabilization change mechanism for complex I.

Locking loop movement in the ubiquinone pocket of complex I disengages the proton pumps.

Complex I (proton-pumping NADH:ubiquinone oxidoreductase) is the largest enzyme of the mitochondrial respiratory chain and a significant source of reactive oxygen species (ROS). We hypothesized that during energy conversion by complex I, electron transfer onto ubiquinone triggers the concerted rearrangement of three protein loops of subunits ND1, ND3, and 49-kDa thereby generating the power-stoke driving proton pumping. Here we show that fixing loop TMH1-2ND3 to the nearby subunit PSST via a disulfide bridge introduced by site-directed mutagenesis reversibly disengages proton pumping without impairing ubiquinone reduction, inhibitor binding or the Active/Deactive transition. The X-ray structure of mutant complex I indicates that the disulfide bridge immobilizes but does not displace the tip of loop TMH1-2ND3. We conclude that movement of loop TMH1-2ND3 located at the ubiquinone-binding pocket is required to drive proton pumping corroborating one of the central predictions of our model for the mechanism of energy conversion by complex I proposed earlier.

Barth syndrome cells display widespread remodeling of mitochondrial complexes without affecting metabolic flux distribution.

Barth syndrome (BTHS) is a rare X-linked disorder that is characterized by cardiac and skeletal myopathy, neutropenia and growth abnormalities. The disease is caused by mutations in the tafazzin (TAZ) gene encoding an enzyme involved in the acyl chain remodeling of the mitochondrial phospholipid cardiolipin (CL). Biochemically, this leads to decreased levels of mature CL and accumulation of the intermediate monolysocardiolipin (MLCL). At a cellular level, this causes mitochondrial fragmentation and reduced stability of the respiratory chain supercomplexes. However, the exact mechanism through which tafazzin deficiency leads to disease development remains unclear. We therefore aimed to elucidate the pathways affected in BTHS cells by employing proteomic and metabolic profiling assays. Complexome profiling of patient skin fibroblasts revealed significant effects for about 200 different mitochondrial proteins. Prominently, we found a specific destabilization of higher order oxidative phosphorylation (OXPHOS) supercomplexes, as well as changes in complexes involved in cristae organization and CL trafficking. Moreover, the key metabolic complexes 2-oxoglutarate dehydrogenase (OGDH) and branched-chain ketoacid dehydrogenase (BCKD) were profoundly destabilized in BTHS patient samples. Surprisingly, metabolic flux distribution assays using stable isotope tracer-based metabolomics did not show reduced flux through the TCA cycle. Overall, insights from analyzing the impact of TAZ mutations on the mitochondrial complexome provided a better understanding of the resulting functional and structural consequences and thus the pathological mechanisms leading to Barth syndrome.