- 1Center for Molecular and Biomolecular Informatics, Radboud Institute for Molecular Life Sciences, Radboud University Medical Center, Nijmegen, Netherlands
- 2Department of Pediatrics, Radboud Center for Mitochondrial Medicine, Radboud University Medical Center, Nijmegen, Netherlands
- 3Department of Medical Microbiology, Radboud Institute for Molecular Life Sciences, Radboud University Medical Center, Nijmegen, Netherlands
- 4Biomolecular Mass Spectrometry and Proteomics, University of Utrecht, Utrecht, Netherlands
- 5Bijvoet Center for Biomolecular Research, University of Utrecht, Utrecht, Netherlands
- 6Utrecht Institute for Pharmaceutical Sciences, University of Utrecht, Utrecht, Netherlands
- 7Netherlands Proteomics Center, Utrecht, Netherlands
- 8University Children’s Research@Kinder-UKE, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
Complexome profiling (CP) is a state-of-the-art approach that combines separation of native proteins by electrophoresis, size exclusion chromatography or density gradient centrifugation with tandem mass spectrometry identification and quantification. Resulting data are computationally clustered to visualize the inventory, abundance and arrangement of multiprotein complexes in a biological sample. Since its formal introduction a decade ago, this method has been mostly applied to explore not only the composition and abundance of mitochondrial oxidative phosphorylation (OXPHOS) complexes in several species but also to identify novel protein interactors involved in their assembly, maintenance and functions. Besides, complexome profiling has been utilized to study the dynamics of OXPHOS complexes, as well as the impact of an increasing number of mutations leading to mitochondrial disorders or rearrangements of the whole mitochondrial complexome. Here, we summarize the major findings obtained by this approach; emphasize its advantages and current limitations; discuss multiple examples on how this tool could be applied to further investigate pathophysiological mechanisms and comment on the latest advances and opportunity areas to keep developing this methodology.
1 Introduction
Proteins, the so-called “workhorses” of life, are the tools that make living machines work (Adams, 2008). The vast majority of biological processes are thus structured, mediated and executed by these macromolecules. Even though many single proteins do perform specific functions by their own, most proteins habitually interact with other proteins, DNA, RNA and lipid molecules rather than acting individually for achieving their biological tasks. The resultant protein complexes can form transient or steady interactions, which also correlate with the type of biological processes they are involved in (De Las Rivas and Fontanillo, 2010). For example, housekeeping cell processes are likely performed by a large fraction of steady protein complexes, whereas in signaling, cell migration, membrane trafficking, metabolic response and other highly dynamic processes, involvement of transient complexes is definitely required for rapid responses and adaptation. The entire set of multi-protein complexes in a cell or compartment is referred to as the complexome (Deshaies et al., 2002; Ceulemans et al., 2006; Lasserre et al., 2006).
In situations where impaired interaction of the elements of a protein complex affect its proper formation; e.g., due to genetic mutations, the related cell process(es) could be compromised and result in biological dysfunction. Higher complexity in multicellular organisms is of course accompanied by larger complexomes than those from unicellular species. Therefore, alterations in protein complexes may affect not only multiple cell processes, but also lead to severe pathophysiological issues at the tissue/organ level. An integral elucidation of both complexomes and protein interaction networks; i.e., interactome (Vidal et al., 2011) under different cellular scenarios, becomes crucial to fully comprehend the molecular mechanisms behind cell physiology and disease.
Numerous biochemical, biophysical, structural, genetic, microscopy and mass spectrometry (MS) approaches have been applied to characterize protein complexes. Although these methods have different principles, in most cases, they require experimental interventions, cell/tissue fractionation, protein extraction and/or time-consuming protocols prior to data collection and analyses. Besides, the amount of information obtained is often limited to one or a small set of protein complexes. The recent breakthrough in quantitative high-throughput technologies and bioinformatic tools has substantially increased the efficiency and quality of the large-scale identification of protein interactors (Iacobucci et al., 2021; Low et al., 2021). Accordingly, an outstanding volume of evidence on the composition, 3D structure, interactions and molecular roles of hundreds of protein complexes is now easily accessible through multiple repositories, such as RCSB PDB (Burley et al., 2021), STRING (Szklarczyk et al., 2021), CORUM (Giurgiu et al., 2019), BioGRID (Oughtred et al., 2021), Pfam (Mistry et al., 2021), Complex Portal (Meldal et al., 2021), NCBI (Coordinators, 2016) and UniProt (UniProt, 2021). Yet, a substantial fraction of the protein interactors reported in those databases still lack full validation of their occurrence in vivo by using novel and more reliable methods.
A systematic approach known as complexome profiling (CP) has recently been introduced as a straightforward and unbiased way to explore protein complexes (Heide et al., 2012; Wessels et al., 2013). CP combines fractionation of protein complexes under native conditions with digestion and quantitative tandem MS-based identification of the derived peptides. Resultant data are clustered to visualize the inventory, abundance and arrangement of multi-protein complexes in a biological sample. Since its formal introduction a decade ago, CP has been primarily used in mitochondrial research (Kiirika et al., 2013; Senkler et al., 2017; Giese et al., 2021; Guerrero-Castillo et al., 2021b; Wittig and Malacarne, 2021); albeit CP can be applied to any kind of cellular fraction.
2 Complexome Profiling
2.1 General Overview
The first requisite to carry out a CP experiment is the collection of biological material; e.g., tissue pieces/biopsies, cultured cells, enriched cell fractions, purified organelles, etc. These materials might also come from previous experimental interventions. Homogenization and cell fractionation methods should aim to keep the native state of protein complexes. For solubilization of membrane proteins, it is recommended to use mild, non-ionic detergents; e.g., Triton X-100, NP-40, dodecyl-maltoside, octyl-glucoside or digitonin (Eubel et al., 2005; Wittig et al., 2006). Next, native protein extracts are separated by a so-called “untargeted” method (Iacobucci et al., 2021), such as native electrophoresis, size-exclusion chromatography or density gradient ultracentrifugation (Figure 1). Regardless of the type of protein separation, CP follows a well-defined protocol that slightly differs among variants. In all cases, protein complexes are separated by size, hence, by shape and molecular mass. Each fraction is further digested with high-purity and specific proteases (e.g., trypsin, chymotrypsin and endoproteinase Lys-C) followed by MS/MS identification of the resulting peptides (Figure 1); i.e., bottom-up proteomics strategy (Zhang et al., 2013).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796128/fcell-09-796128-HTML/image_m/fcell-09-796128-g001.jpg)
FIGURE 1. Overall workflow of complexome profiling (CP). After collection, homogenization and fractionation of biological materials, proteins are separated by either native polyacrylamide gel electrophoresis (PAGE), size exclusion chromatography (SEC) or density gradient ultracentrifugation (DGU) for CP studies. The obtained protein-containing fractions are individually digested with specific proteases. Resultant peptides are extracted, cleaned and usually separated by reversed-phase high-performance liquid chromatography (RP-HPLC) followed by tandem mass spectrometry (MS/MS) analysis. MS data can be acquired in data-dependent or data-independent modes, DDA or DIA, respectively. Next, a proteomic search is performed to match obtained MS spectra against proteome databases using a variety of available software. Icons of the most popular tools used for CP studies are shown (see Section 2.1 for details). Protein abundance profiles are further obtained by plotting LFQ/iBAQ values against the number of fractions. The list of identified protein groups is computationally sorted based on similarities of the abundance patterns across fractions; e.g., hierarchical clustering. Prior to analysis, complexome data are pre-processed to account for protein loading/MS sensitivity differences. Data correction and normalization between samples are regularly applied. In SEC- and native-PAGE-based CP, a mass calibration can be implemented for a more meaningful biological dimension. Data can be analyzed by using available bioinformatic tools specifically created for CP. Some of these programs are shown in the figure (see Section 2.4 for details). For reusing and further analysis, large CP datasets can be deposited and publicly shared through PRIDE (Proteome Xchange consortium) or the new dedicated website for CP, CEDAR.
Previous to MS analysis, peptides are usually separated based on their hydrophobicity by reversed-phase high-performance liquid chromatography (RP-HPLC) followed by electrospray ionization (ESI) to produce gas phase ions (Zhang et al., 2013). The majority of MS data for CP studies has been acquired in data-dependent mode (DDA); i.e., collection of a predefined number of precursor ions for fragmentation during the MS2 cycle is done according to their charge states and abundance (Hu et al., 2016). In conventional DDA modes, the possibility to identify low abundant peptides thus becomes significantly limited. To avoid this issue, acquisition strategies in data-independent mode (DIA) (Krasny and Huang, 2021), such as sequential window acquisition of all theoretical mass spectra (SWATH-MS) (Gillet et al., 2012), in which fragmentation of all precursor ions identified in the MS1 cycle can be analyzed in the MS2 cycle, have recently been introduced to CP-like workflows (Heusel et al., 2019; Bludau et al., 2020; Calvo et al., 2020).
MS spectra are routinely matched against a proteome database of the species of interest by using a variety of software/search engines (Figure 1), such as Mascot (Perkins et al., 1999), openMS (Rost et al., 2016), PEAKS DB (Zhang et al., 2012), Proteome discoverer (Orsburn, 2021), Protein Prospector (Chalkley et al., 2005) and MaxQuant (Tyanova et al., 2016a). In case that the proteome of the studied organism is not yet available or fully annotated, metagenome or transcriptome data are useful to generate a list of proteins and perform the search. Alternatively, peptide sequences might be deduced directly from tandem MS spectra by using database-independent computational approaches; i.e., de novo peptide sequencing (Tran et al., 2017).
Protein abundance profiles are further obtained by plotting, for example, label-free quantification (LFQ) (Cox et al., 2014) or intensity-based absolute quantification (iBAQ) (Schwanhausser et al., 2011; Tyanova et al., 2016a) values against the number of fraction. While LFQ intensities are typically utilized for comparing relative amounts of proteins in multiple samples, iBAQ values offer a more stoichiometric impression of the identified protein groups as those are proportional to their molar quantities. iBAQ values are calculated as the sum of all individual peptide intensities of a given protein group divided by its number of theoretical identifiable peptides. It is thus not surprising that the majority of CP studies use iBAQ since it enables a fair comparison between different protein groups in the same and multiple samples. Additional information regarding LFQ and other alternative methods have been reviewed in Fabre et al. (2014) and Wittig and Malacarne (2021). Ultimately, the list of hundreds or even thousands of identified proteins is mainly sorted based on similarities of the abundance patterns across fractions by hierarchical clustering analysis (Figure 1). Proteins that are part of the same complex consistently cluster together since these show co-migration in the same fraction(s) and similar abundance profiles.
2.2 Different Complexome Profiling Setups, One Common Goal
As a ground-breaking “-omics” method, CP has been rapidly developing and spreading among the scientific community. It is thus not striking that in less than a decade, multiple strategies have complemented or even improved its scope. At the same time, several limitations of CP have been circumvented gradually by including novel MS-based strategies, specific adaptions in sample processing and new tools for data analysis. All variants of CP do share a common goal nonetheless: to unravel the composition of protein complexes, as well as their abundance, stabilities, apparent molecular masses and stoichiometries.
The most used CP workflow, currently referred to as “classic CP” (Wittig and Malacarne, 2021), uses Blue Native-polyacrylamide gel electrophoresis (BN-PAGE) followed by MS identification. The main advantages of BN-PAGE are the relatively low amounts of biological material required and high-resolution separation of native proteins. BN-PAGE conditions keep proteins in similar native states to those occurring in vivo (Wittig et al., 2006). However, if the presence of Coomassie blue dye affects the stability of one or more protein complexes of interest, milder dye-free versions of native-PAGE can easily be used instead; e.g., high-resolution clear native-PAGE (hrCN-PAGE) (Wittig et al., 2007; Ladig et al., 2011). Native-PAGE is suitable for separation of proteins between ∼0.02–10 MDa. If the study requires interrogation of bigger protein complexes (>10 MDa), large pore BN-PAGE (LP-BN-PAGE) can be applied instead (Strecker et al., 2010; Heide et al., 2012).
After electrophoresis, entire gel strips are fixed and often cut into 32–70 slices followed by LC-MS/MS identification (Heide et al., 2012; Senkler et al., 2017; Vidoni et al., 2017; Giese et al., 2021). For CP analysis, slice numbers can be transformed into apparent molecular mass values by carrying out a calibration using standard proteins as exemplified in Wittig et al. (2010) and Heide et al. (2012). Molecular mass accuracy is optimized by using different sets of standard protein complexes, one for water-soluble proteins and the other for membrane proteins. On the whole, the more fractions collected the better resolution of the resultant profiles will be. In 2016, Müller and co-workers developed cryo-slicing Blue Native-Mass Spectrometry (csBN-MS), which helped increase the resolution, accuracy and quantification of the complexome profiles by slicing the BN-lanes in 230 pieces as well as optimizing MS analysis and protein identification (Muller et al., 2016). The resolution of complexome profiles obtained from ≤30 slices would not be suitable to unambiguously assign protein interactions.
Besides native-PAGE, size exclusion chromatography (SEC) and density gradient ultracentrifugation (DGU) have been implemented in different CP setups. Using SEC, native proteins are separated by filtration through a gel matrix (resin), which consists of spherical porous beads of different sizes depending on the desired range of molecular masses (Burgess, 2018). Elution of heavier proteins is faster than the lighter ones by this method. Apparent molecular masses of eluted fractions can be determined with proper standard mixes; although the mass calibration is more accurate for globular proteins (Hong et al., 2012; Korepanova and Matayoshi, 2012). SEC resins are suitable for separating proteins in a wide range of molecular masses; e.g., 1–700 kDa (Sephadex™, Superdex™); 5–5,000 kDa (Superose™) and 0.01–40 MDa (Sepharose™). The number of fractions that can be collected after SEC is comparable to those aforementioned. For example, in the recently introduced SEC-SWATH-MS approach (Heusel et al., 2019), 81 fractions were collected, subjected to MS identification and examined by complex-centric proteome analysis. This CP setup led to higher sensitivity and accuracy in protein quantification, less noise, and validation of protein interactions. SEC has been particularly useful in CP studies characterizing large protein complexes, such as nuclear components (Connelly et al., 2018). However, key limitations of SEC-involving CP setups would be the larger amounts of biological material required, loss of protein interactions by dilution, formation of self-oligomers and inadequate identification of membrane and low abundant protein complexes (Burgess, 2018; Heusel et al., 2019; Iacobucci et al., 2021).
On the other hand, DGU uses solutions of different densities made of glycerol, sucrose, cesium chloride, iodixanol or Ficoll® through which protein complexes are separated based on their sedimentation rates, where heavier complexes sediment faster. DGU is an excellent technique for separating large protein complexes, ribosomes, membrane vesicles and subcellular organelles. DGU is also suitable to separate cleared cell lysates as well as immunoprecipitation (IP)-captured protein complexes (Lee and Skalnik, 2013; Caudron-Herger et al., 2019). After separation, most of resolved protein complexes remain in near-native states. In a recent study, a reliable DGU-based CP variant, referred to as quantitative density gradient analysis by mass spectrometry (qDGMS) has been developed to study the human mitoribosome (Palenikova et al., 2021a). qDGMS combines stable isotope labeling by amino acids in cell culture (SILAC) and DGU followed by fractionation and LC-MS/MS analysis. SILAC is a MS-based technique that quantifies the differences in protein abundance/expression among biological samples (Ong et al., 2002). In conventional SILAC, two different cell lines/strains are cultured using media supplemented with either “heavy” or “light” essential aminoacids that are labeled with non-radioactive isotopes (e.g., 13C, 2H, 18O, 15N) or unlabeled, respectively (Geiger et al., 2011). The specific labeled aminoacids are hence incorporated into all cell proteins during cell growth. Proteins of cell lysates/fractions obtained from the two samples are mixed (1:1), digested and analyzed together by LC-MS/MS. The isotope-labeled peptides appear in MS spectra as pairs with identic chemical composition but different masses. Ratios of intensities for the many identified peptide pairs thus denote the respective changes in protein abundance between samples. Incorporation of SILAC not only in qDGMS but also in classic CP allows duplexing, MS time-saving and, most important, higher accuracy of quantification in proteomic analysis of two experimental conditions (Palenikova et al., 2021a; Palenikova et al., 2021b).
A major challenge of DGU relies on retrieving of fractions without manual disturbance of the resolved layers. To account for this issue, several strategies and devices have been developed, including commercially available automatic fractionators and freezing the gradient after centrifugation followed by cryo-slicing to fractionate samples consistently (Yu et al., 2016). Location of resolved proteins by DGU is considerably more spread when compared to native-PAGE or SEC, which means that increasing the number of fractions does not necessarily lead to higher resolution. Furthermore, sedimentation rates do not only depend on molecular masses of protein particles but also on shape and densities from both the particles and fluid used for making the gradient (Cole et al., 2008). For these reasons, co-migration of identified proteins by this CP setup does not immediately represent actual associations rather than merely similar sedimentation rates.
2.3 Other Mass Spectrometry-Based Approaches Related to Complexome Profiling
Additional untargeted MS-based assays have been described in the literature, which share the same principles of a CP experiment. The group of Matthias Mann developed protein correlation profiling (PCP) by combining DGU and quantitative MS to identify and validate several components of the human centrosome (Andersen et al., 2003). A few years later, Harner and co-workers applied PCP to discover the mitochondrial contact site (MICOS) complex (Harner et al., 2011). The further inclusion of hierarchical clustering and bioinformatic analysis of PCP data improved substantially this method leading to a broader characterization of both organelle and cell complexomes (Foster et al., 2006; Havugimana et al., 2012). PCP has also been related to co-fractionation coupled to MS (CoFrac-MS) experiments (Bludau, 2021). In CoFrac-MS, ion-exchange chromatography (IEX) or SEC have however been used for a higher resolution separation of protein complexes by their charges or hydrodynamic radii, respectively (Bludau and Aebersold, 2020). The combination of PCP/CoFrac-MS with other techniques, such as SILAC or stable isotope labelling in mammals (SILAM) have been proven valuable in proteome-scale interactome studies across cells (Kristensen et al., 2012) and tissues (Skinnider et al., 2021).
CP-like approaches using BN-PAGE followed by MS analysis have also been applied in different studies. A systematic proteomic analysis of series of 15 protein bands excised from individual BN-gel lanes allowed the identification of subunits of the mitochondrial respiratory complexes and was helpful to correctly re-assign NDUFA4, previously misidentified as a complex I subunit, as a bona fide component of the cytochrome c oxidase (Balsa et al., 2012). Likewise, Cogliati and co-workers separated mitochondrial proteins from mice by BN-PAGE and cut the strips in 26 slices before MS identification (Cogliati et al., 2016a) using a data-independent scanning (DiS) method (Guaras et al., 2016). The latter workflow has recently been standardized and referred to as Blue-DiS (Calvo et al., 2020).
2.4 Complexome Profiling Data Analysis and Visualization
Huge output files are obtained in CP studies after protein searches, which usually contain all data necessary for further analysis, including identified protein groups, unique peptides, sequence coverage, MS/MS counts, scores, LFQ/iBAQ values from each fraction, etc. The most common visualizations of a complexome profile are heatmaps accompanied by line charts plotting the LFQ/iBAQ values throughout the fractions (Figure 1). Complexome profiles of a short set of proteins can be manually generated and analyzed. Yet, the large volume of protein identifications in CP datasets makes full manual inspection impractical. In the last years, several tools have been specifically designed for automated processing and exploration of CP datasets.
ComplexomeMap has been developed as an online public platform for mining the mitochondrial complexome of plants Arabidopsis thaliana and Viscum album (Senkler et al., 2017). NOVA is a user-friendly software to perform cluster analysis, mass calibration, normalization, visual inspection, links to protein databases and comparison of experimental conditions (Giese et al., 2015). The software COPAL has proven helpful for analyzing multiple CP datasets, aligning experimental replicates and detecting significantly affected protein complexes (Van Strien et al., 2019). It also generates files that can be directly used for gene set enrichment analysis. ComPrAn, a Shiny R app, has been developed for analyzing qDGMS data; this tool enables analysis of peptide-level data, normalization, clustering, visualization options and a graphical user interface (Palenikova et al., 2021a). In addition, ComPrAn is particularly useful to analyze proteomic data obtained from SILAC-treated samples. ComplexFinder, a Python-based software suit, has recently been released to analyze fractionation of native protein complexes, particularly from BN-PAGE- or SEC-based CP experiments (Nolte and Langer, 2021). This tool allows machine learning-based prediction of potential protein-protein interactions (PPIs) and high flexibility in CP data analysis. ComplexFinder also provides improved peak-centric quantification and kinetic modelling, protein connectivity networks and compatibility with different quantification strategies. CCprofiler is a robust R package for analyzing co-fractionation MS datasets, which has been originally designed for SEC-SWATH-MS (Heusel et al., 2019). CCprofiler and its web interface, SECexplorer, offer multiple functions, such as quality control and filtering for less erroneous assignment of protein interactors, usage of curated reference datasets, protein quantification, protein- and complex-centric data analysis and visualization.
Other software available for proteomic and interactome analysis can also be used for CP data processing. For instance, free tools such as Perseus (Tyanova et al., 2016b), PrInCE (Stacey et al., 2017) and EPIC (Hu et al., 2019) may provide suitable options for visualization, protein quantification, statistical analysis, prediction of PPIs from co-fractionation data, obtention of supporting info from public repositories and/or cross-omics comparisons.
2.5 Storage, Sharing and Examination of Complexome Profiling Data
Although CP datasets contain a long list of identified proteins, the majority of studies have focused only on a small subset of the identified protein complexes. Naturally, the amount of information contained in those datasets could be useful not only to unveil uncharacterized PPIs, but also to validate previously reported protein complexes. All these under-explored results thus represent a valuable source of information, which could be reused and further mined by other researchers. To achieve this purpose, the ComplexomE profiling DAta Resource (CEDAR) has recently been established as the first website for depositing and sharing CP data (van Strien et al., 2021). A standardized format containing the minimum information required for a CP experiment (MIACE) ensures compatibility and correct retrieval of data. CEDAR is able to link access to raw MS data and other files deposited parallelly in the known repository PRIDE (ProteomeXchange consortium) (Perez-Riverol et al., 2019). Additionally, CEDAR has a profile viewer tool that allows users to examine CP data directly on the website. Since its release in June 2020, 24 complete CP experiments encompassing a total of 146 samples have been submitted to CEDAR (web V1.1, accessed October 6, 2021). CP data from 11 species covering different life domains are already accessible (Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796128/fcell-09-796128-HTML/image_m/fcell-09-796128-g002.jpg)
FIGURE 2. Different species used in complexome profiling studies. Phylogenetic tree including all the species reported in complexome profiling (CP) studies until November 2021. As it can be seen, CP has been used to study multiple species covering the three different life domains. Where a tree-like icon is shown next to the species name, CP data are currently available on CEDAR (https://www3.cmbi.umcn.nl/cedar/browse/). The template of the tree was generated online using Interactive Tree of Life (iTOL) v6.3.3 (https://itol.embl.de/) and further edited manually.
It should be mentioned that the data available through CEDAR are just a fraction of the CP studies reported in the literature. At present, ∼145 publications mentioning the term “complexome” (or complexomic/complexomics) have been found in PubMed (NCBI) records; ∼62 publications correspond to CP-related studies, reviews, methods and software tools (Figure 3). In ∼40 and ∼18 publications, CP has been used for mitochondrial research or to study protein complexes involved in disease, respectively (Figure 3). As the number of CP publications continues accelerating, it is desirable that future studies using this approach would also share resultant CP datasets publicly through CEDAR to simplify their reuse by the scientific community and keep expanding its valuable offer.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796128/fcell-09-796128-HTML/image_m/fcell-09-796128-g003.jpg)
FIGURE 3. Reported studies in the literature including the complexome profiling approach. Since the introduction of the term “complexome” back in 2002 (*), ∼145 publications referring to this term have been reported in PubMed (NCBI) (blue line). Formal introduction of complexome profiling (CP) in 2012 (**). In less than 10 years, ∼62 studies have applied CP (orange line), where ∼40 and ∼18 publications are related to mitochondrial research (green line) or to protein complexes involved in disease (purple line), respectively. The positive tendency of multiple groups around the world implementing CP to study protein-protein interactions in both fundamental- and disease-related studies as a standard method is incontrovertible.
3 Complexome Profiling as a Tool to Investigate Mitochondrial Protein Complexes
Mitochondria are traditionally known as the eukaryotes’ “powerhouses” since they contain all the enzymes required to generate ATP via the oxidative phosphorylation (OXPHOS) pathway. This process is catalyzed canonically by four respiratory chain complexes (I, II, III and IV) and a F1FO-ATP synthase (complex V). Complexes I, III and IV couple the energy released during electron transfer to oxygen to the generation of an electrochemical gradient of protons across the inner mitochondrial membrane (IMM). The proton gradient not only drives the synthesis of ATP, but also other processes such as metabolite transport, protein import, redox and Ca2+ homeostasis, fusion/fission, signaling and cell death. These organelles do also contain their own genome, also known as mitochondrial DNA (mtDNA), which encodes for a few but essential subunits of OXPHOS complexes. The roles of mitochondria are thus not limited to their energy duties. These organelles are multi-functional “hubs” critical to almost every cell process.
Despite of the great progress in molecular characterization and structure elucidation of numerous mitochondrial protein complexes, mostly OXPHOS- and mitoribosome-related proteins, a large number of mitochondrial PPIs remain elusive. In the last years, however, CP has proven valuable for identifying novel protein interactors, validating previous results and shedding light on intricate assembly pathways. In this section, we summarize these findings and describe how CP has boosted the analysis of protein complexes in the mitochondrial research field. Expectedly, most of these findings are again related to components or mediators of the assembly of OXPHOS complexes, energy metabolism-linked proteins and mitoribosomes.
The first groups that established CP have also been interested in the many features of OXPHOS and in particular of complex I (CI). CI is the largest redox enzyme of the mitochondrial respiratory chain constituted by ∼45 subunits depending on the species. CI generates proton-motive force driven by the transfer of electrons from NADH to ubiquinone (Hirst, 2013). Although the redox features and composition of mitochondrial CI in several species were already known by the start of the 2010s, major queries on this enzyme had yet unresolved: its entire 3D structure, its energy-conversion mechanism and how it assembles. To help tackle the last one, the formal introduction of CP by Heide and co-workers was useful to identify TMEM126B interacting with the known CI assembly factors NDUFAF1, ECSIT and ACAD9 (Heide et al., 2012). TMEM126B knockdown in 143B osteosarcoma cells led to ∼95% specific decrease in CI-containing supercomplexes, i.e., supramolecular associations of complexes I, III and IV. These results thus proposed TMEM126B as a CI assembly factor. In an earlier report using HEK293 cells, Wessels and co-workers identified with a similar approach two other assembly factors of CI: C6ORF66 and C3ORF60 (Wessels et al., 2009), currently known as NDUFAF4 and NDUFAF3, respectively. In this study, mitochondrial proteins were separated by BN-PAGE followed by LC-MS/MS identification, but comparison of migration patterns was limited to the use of PCP.
A few years later, a dynamic CP strategy was successfully implemented to describe the step-by-step integration of the subunits, assembly factors and the different assembly intermediates of human CI (Guerrero-Castillo et al., 2017a). The authors used CP data from time-based mitochondrial translation recovery to describe a number of assembly intermediates of CI accumulating at various timepoints after removal of chloramphenicol, a reversible inhibitor of mitoribosomes (Ugalde et al., 2004). These data corroborated the long time proposed modular assembly pathway of CI, which basically involves the coordinated formation and pre-assembly of its functional modules, N, Q, PP and PD, before forming the entire enzyme. This study also provided insight into the specific involvement of earlier reported assembly factors and novel interactors. Complex IV (CIV) assembly-related chaperone COA1 and TMEM186 were clearly found associated with membrane arm intermediates that also interact with TMEM126B, NDUFAF1 and ECSIT. At present, all these proteins have been recognized as true components of the so-called mitochondrial CI assembly (MCIA) complex (Formosa et al., 2020). ATP5SL clustered with FOXRED1, another CI assembly factor. However, ATP5SL disruption did not lead to CI deficiency (Andrews et al., 2013) and its involvement in the assembly pathway is thereby unclear. Complex V assembly-involved protein TMEM70 was also identified in assembly intermediates of CI. To further explore the putative role of TMEM70 on CI assembly, Sánchez-Caballero and co-workers performed a study including proximity-dependent biotin identification (BioID), co-evolution analyses and CP (Sanchez-Caballero et al., 2020). In this study, TMEM70 was found in close contact with subunits of complexes I and V, as well as the small mitoribosome subunit. The absence of TMEM70 resulted in slight depletion of fully assembled CI and accumulation of assembly intermediates of the distal part of the membrane arm of CI, whereas mature complex V (CV) content was ∼70% lower than in the control cell line. These authors showed TMEM70 as a non-essential assembly factor of complexes I and V, and proposed a possible role in tethering mitoribosomes to the IMM during translation.
Apart from studies in human cells, CP has also been used to study the role of several accessory subunits of CI in different models. For example, Kmita and co-workers analyzed the assembly of CI in a deletion strain of the yeast Yarrowia lipolytica lacking the Zn2+-containing subunit NUMM/NDUFS6 (Kmita et al., 2015). Absence of this subunit did not prevent formation of the entire CI, since a lower enzyme content and higher fraction of peripheral arm subcomplexes were observed. Unexpectedly, the almost fully assembled CI contained the assembly factor N7BML/NDUFAF2 bound instead of accessory subunit N7BM/NDUFAF12. Moreover, Angerer et al. implemented CP in a study on the mitochondrial acyl carrier proteins (ACPM) 1 and 2, which are accessory subunits of CI and contain LYR motifs (Angerer et al., 2014). Although ACPM1/2 are predominantly associated to CI, ACPM1 was also found as a free protein in the matrix and forming a complex with LYRM4(ISD11)/NFS1 involved in iron-sulfur cluster biosynthesis. Furthermore, two different groups in parallel unveiled the existence of two isoforms of the peripheral arm subunit NDUFV3 of CI in mammals (Bridges et al., 2017; Guerrero-Castillo et al., 2017b). An extra exon present in gene NDUFV3 can be alternative spliced, hence generating short and long isoforms of ∼10 and ∼50 kDa, respectively. CP helped find a different expression of these isoforms in different tissues of bovine, mouse and rat as well as in cultured human cells. The canonical short isoform was predominantly identified in heart and skeletal muscle, whereas the large isoform was the foremost isoform in liver, brain and lung tissues. Both NDUFV3 isoforms can also be expressed at the same time and correctly assembled onto CI in some tissues; yet, one isoform predominated in each case.
Other OXPHOS complexes and their chaperones have also been explored by CP. Singhal et al. identified the product of ORF YDR381C-A as an assembly factor of yeast complexes III and IV (Singhal et al., 2017). This protein was renamed as cytochrome c oxidase interacting protein 1 (Coi1), which occurs only in fungi. Deletion of Coi1 resulted in severe alteration of mitochondrial function, diminished amounts of CIV-associated heme, defective assembly of complexes III, IV and their supercomplexes, as well as accumulation of assembly intermediates. Vidoni and co-workers reported that the short isoform of myofibrillogenesis regulator 1 (MR-1S) associates with chaperones PET100 and PET117 to mediate human CIV assembly (Vidoni et al., 2017). Authors implemented not only CP to study assembly intermediates accumulated in control and MT-CO3 mutant cybrids, but also SILAC and quantitative MS.
CP has been helpful to better understand the formation of supercomplexes or respirasomes. The specific factors for mediating this process remain unclear. COX7A2L, also known as SCAFI, has been shown critical for the association of complexes III and IV (Perez-Perez et al., 2016). Its putative involvement in larger respirasome formation has also been reported (Lapuente-Brun et al., 2013). To better understand the role of SCAFI, Férnandez-Vizarra et al. explored the mitochondrial complexome of a SCAFI knockout human cell line by SILAC-based CP (Fernández-Vizarra et al., 2021). Absence of SCAFI resulted in marked loss of supercomplex III2-IV, whereas CI-containing respirasomes were not affected. In contrast, authors showed that ∼70% of respirasomes contained COX7A2 in either human cell line analyzed. This discrepancy has apparently been related to a tissue-specific expression of the two isoforms (Lapuente-Brun et al., 2013). On the other hand, Protasoni and co-workers demonstrated that absence of a fully assembled complex III (CIII) in a MT-CYB-deficient human cell line stalled CI biogenesis by preventing the integration of its N-module (Protasoni et al., 2020). Although substantial accumulation of partially assembled Q/P intermediate was observed by SILAC-based CP, a slight fraction of CI could still be detected. Reasonably, respirasome assembly was totally lost in the mutant. Besides, assembly of CIV was affected since several COX subunits were found associated with CIII sub-assemblies; hence interfering their correct assembly. It has thus been proposed that human complexes I, III and IV assemble in a cooperative fashion, where CIII seems to be central.
Mitochondria are characterized by a highly folded IMM. These folds, called cristae, have emerged as dynamic compartments whose shape and dimensions influence structure and functioning of the OXPHOS system (Cogliati et al., 2016b). A key player in shaping cristae appears to be the MICOS complex, which also interacts with the sorting and assembly machinery (SAM) complex from the outer mitochondrial membrane (OMM) (Huynen et al., 2016). Interaction of these two complexes constitutes the mitochondrial intermembrane space bridging (MIB) complex (Ott et al., 2015). Since its relatively recent discovery (Harner et al., 2011), much of what we know about the MICOS complex has been derived from CP-based studies.
For instance, Weber and co-workers identified apolipoprotein O (APOO/MIC26) and apolipoprotein O-like protein (APOOL/MIC27) as potential components of the human MIB complex (Weber et al., 2013). Huynen and co-workers described the composition and apparent masses of the fully assembled MIB complex (2.2–2.8 MDa) and other putative assembly intermediates, e.g., the free MICOS complex (∼700 kDa) (Huynen et al., 2016). Anand et al. identified MIC13 (QIL1) as a novel MICOS complex subunit (Anand et al., 2016). Deletion of MIC13 resulted in a smaller albeit still assembled MICOS complex while having no effect on integrity of OXPHOS complexes. They also demonstrated a disruptive effect of MIC13 deletion on cristae morphology accompanied by reduced respiratory capacity. The same group used a similar approach to characterize MIC26 and MIC27 (Anand et al., 2020). They showed negative effects on CV stability as well as cristae morphology defects both of which were much more pronounced in the double knockout than either single knockout, suggesting overlapping roles. Assembly of other MICOS components was however unimpeded by MIC26/27 knockouts, proposing these proteins are not essential for its assembly or stability. The link of MICOS components to CV was further explored by Eydt et al. (2017), showing that MIC10 comigrates with CV dimers. Additionally, they presented evidence that MIC10 acts antagonistically to MIC27 in negatively controlling CV oligomerization, while MIC27 appears to be a positive regulator. A recent paper by Bock and co-workers reported the interaction of PGC-1- and ERR-induced regulator in muscle 1 (PERM1) with multiple components of the MICOS/MIB complex as well as vimentin and ankyrin B in skeletal muscle from mouse (Bock et al., 2021). These findings suggested a novel mechanism to help understand not only the interconnection of mitochondria and sarcolemma but also mitochondria-cytoskeleton associations as well as the organization of a functional mitochondrial network in this tissue. Taken together CP has greatly aided our understanding of MICOS/MIB complex composition, the role of individual subunits and its interaction with OXPHOS complexes.
The IMM contains other key protein complexes involved in proteolytic events, fusion/fission of mitochondrial membranes, regulation of cristae morphology and cell signaling. Some of these proteins belong to the SPFH (stomatin, prohibitin, flotillin and HflC/K) family and usually arrange as large scaffolds. To have a clearer picture, Wai and co-workers implemented CP to unveil the interaction partners of stomatin-like protein 2 (SLP2) and at least two proteases, PARL and YME1L, in immortalized embryonic fibroblasts mitochondria (Wai et al., 2016). SLP2 works as a regulator of PARL, which also modulates other proteins such as PGAM5, PINK1 and OMA1. These interactions were suggested to occur in defined sites of the IMM and related to mitochondrial proteostasis, dynamics and cell survival. Similarly, Konig and co-workers identified MAIP1 (C2Orf47) as a novel m-AAA protease-binding protein in complexome profiles (Konig et al., 2016). MAIP1 was shown essential for regulating the assembly of subunit EMRE into mitochondrial Ca2+ uniporter (MCU) complexes.
CP has also played an important role in the field of plant and parasite mitochondrial biology. Using CP, Senkler and co-workers performed a systematic characterization of mitochondrial complexes in the plant A. thaliana (Senkler et al., 2017). CP data were also useful to update the subunit composition of OXPHOS complexes and identifying respective assembly intermediates. This group then used the same approach on mitochondria from the mistletoe V. album, revealing a highly unusual OXPHOS system (Senkler et al., 2018). This species completely lacks CI, the only multicellular eukaryote to date with this distinction; instead, V. album contains alternative NAD(P)H oxidoreductases. This species also expresses an alternative oxidase and its complexes III and IV are firmly associated as supercomplexes. In comparison to A. thaliana, the abundance of complexes II and V was particularly low, suggesting a shift in stoichiometry of the OXPHOS complexes in the IMM of this plant. Rugen and co-workers took a more focused CP strategy to investigate the composition of the mitoribosome in A. thaliana (Rugen et al., 2019). Utilizing LP-BN-PAGE- and DGU-based CP setups, several non-conventional proteins were found attached to the mitoribosome, mostly from the class of pentatricopeptide repeat proteins that seem to be involved in RNA processing and protein maturation. Presence of these additional interactors results in a larger mitoribosome with unusual large and small subunits of ∼3 and ∼5.7 MDa, respectively. Additional functions are thus likely incorporated into the plant mitoribosome.
Due to their extreme divergence, apicomplexan parasites of the genera Plasmodium and Toxoplasma that cause the infectious diseases malaria and toxoplasmosis, respectively, remain poorly understood. In fact, more than one-third of their genes still lack any functional annotation (Aurrecoechea et al., 2009; Harb et al., 2020). To help narrow this gap, Evers et al. and Maclean et al. obtained evidence for highly divergent composition of the mitochondrial OXPHOS complexes in Plasmodium falciparum and Toxoplasma gondii, respectively (Evers et al., 2021; Maclean et al., 2021). More than 30 novel subunits across complexes II, III, IV and V with no recognizable orthologs outside of the myzozoan phylum could be identified by CP. These novel subunits did not only replace subunits typically observed in standard models, but also increased the size of P. falciparum OXPHOS complexes by around 50, 50, 130 and 70% as compared to complexes II, III, IV and V from mammalian mitochondria, respectively. This was consistent with the observations in T. gondii, thus making apicomplexan OXPHOS complexes the largest described to date. Furthermore, abundance of OXPHOS complexes was ∼32 fold increased in the transmissible gametocyte stages compared to the pathogenic asexual stages of P. falciparum (Evers et al., 2021). This finding offered protein-level support for the long-standing hypothesis that malaria parasites undergo a metabolic switch towards mitochondrial catabolism to facilitate their transmission back to the insect vector (MacRae et al., 2013). Despite its relatively recent introduction to the field of parasite research, CP has already provided valuable data for mapping a number of previously unknown interactors to protein complexes and biological processes. As these novel additions are unlike what is known from standard models, it opens the door for new discoveries, even in pathways as fundamental as respiration.
4 The Mitochondrial Complexome and Human Disease
In the last years, several metabolic diseases have been characterized along with accompanying complexome studies of human subjects carrying mutations in subunits of the OXPHOS complexes, assembly factors or proteins required for efficient mitochondrial protein biosynthesis. Most of these studies involved patients with metabolic disorders that affect mitochondrial processes leading to heterogeneous clinical manifestations ranging from mild to severe phenotypes. These studies analyzed subcellular preparations from primary skin fibroblasts from patients and controls, or from cell lines with genetic mutations/deletions used as disease models. CP analysis revealed either changes in abundance of specific protein complexes, indicating alterations of their relative steady-state levels in comparison to the controls, or molecular mass shifts resulting from the loss of protein components or from destabilization of the complexes. Beyond the described main alteration of the complex directly involved with the respective genetic mutation, an enormous amount of information can be retrieved from the deposited complexome data of these studies that could potentially be re-analyzed and investigated, for instance, in relation to the clinical features of these patients. Such in-depth characterization of the mitochondrial proteome at the level of macromolecular complexes in disease states may conceivably help uncover molecular basis and consequences of metabolic disorders and, importantly, could serve to identify molecular targets for developing novel therapeutic strategies. In this section, we summarize studies of genetic defects causing metabolic diseases in human subjects and mutations in model organisms analyzed to date by CP (Table 1).
4.1 Complex I Deficiency
CI deficiencies may potentially originate from genetic defects in either the seven genes encoded in the mtDNA or in the 37 subunits and the increasing number of assembly factors encoded in the nuclear DNA (Fassone and Rahman, 2012; Sánchez-Caballero et al., 2016a). Pathogenic mutations in patients have been found in the majority of CI subunits and in many assembly factors (Rodenburg, 2016), and the number of mutations identified continuously grows. Since TMEM126B was the first assembly factor of CI identified by CP, it was not unexpected that the first complexomes from patient fibroblasts with genetic alterations described in the literature belonged to subjects harboring mutations in this gene. Independently, two research groups described six (Alston et al., 2016) and three (Sánchez-Caballero et al., 2016b) subjects, respectively, presenting isolated CI deficiency with muscle weakness in most cases and with recessive mutations, evidenced by whole exome sequencing, in the gene encoding TMEM126B. Complexome profiles of three subjects revealed aberrant migration patterns of CI subunits with decreased abundance of mature CI-containing supercomplexes at ∼1.7 MDa and accumulation of assembly intermediates. In both studies, increased steady-state levels of subunits belonging to the Q-module plus assembly factors NDUFAF3 and NDUFAF4, and to the P-distal ND4 module indicated stalled assembly intermediates of these unaffected modules, suggesting that TMEM126B takes part, together with the other MCIA complex components, ACAD9, NDUFAF1 and ECSIT, in the assembly of the P-proximal ND2 module. Complementation of patient fibroblasts using lentiviral expression of wild-type TMEM126B partially restored the abundance of CI-containing supercomplexes and decreased the abundance of Q- and P-module assembly intermediates, validating its role as assembly factor. Notably, specific ubiquitination and further degradation of this assembly factor has been implicated in a regulatory mechanism to decrease CI content and consequently oxygen consumption under hypoxic conditions (Fuhrmann et al., 2018).
Complexome profiles of patient fibroblasts with mutations either in a subunit or in an assembly factor of the Q-module, NDUFA6 and NDUFAF8, respectively, proved useful to corroborate the pathogenicity of the respective mutation and for pointing out the stage in which CI assembly was stalled. Bi-allelic mutations in NDUFA6, a LYR-motif subunit important for proper functioning of ubiquinone redox chemistry (Galemou Yoga et al., 2020), were found responsible for defects in CI assembly (Alston et al., 2018). In one of the three subjects studied, it has been found a decrease not only in the subunits of the Q-module, NDUFA6, NDUFA7 and NDUFA12, but also in subunits belonging to the N-module, correlating with the more severe phenotype observed in this subject. Of note, in these three patients, assembly factor NDUFAF2, which is known to interchange for subunit NDUFA12 in a late assembly step (Pereira et al., 2013), remained to some extent bound to CI intermediates, indicating that also in the other two subjects N-module incorporation was inefficient. On the other hand, NDUFAF8 has been identified to physically interact with NDUFAF5 (Floyd et al., 2016), an assembly factor of CI required to hydroxylate an arginine residue of the Q-module subunit NDUFS7 (Rhein et al., 2016). Alston and co-workers identified three subjects presenting Leigh syndrome with bi-allelic mutations in NDUFAF8 causing isolated CI deficiency (Alston et al., 2020). Complexome profile of one of the patients revealed decrease abundance of CI-containing supercomplexes. Consistent with the role of this protein in early CI assembly, NDUFAF8-mutated patient fibroblasts showed incomplete formation of the Q-module since neither its signal nor the fraction of this module holding assembly factors NDUFAF3/4 were observed. These fractions are visible even in control cells where CI assembly process occurs normally; in particular the more prominent one including the assembly factors. In contrast, accumulation of the earliest Q-module sub-assembly formed only by subunits NDUFA5, NDUFS2 and NDUFS3 has been described. In addition, accumulation of the P-distal ND4 module has also been reported.
Extending the genotype of CI deficiency, pathogenic mutations in subunit NDUFC2 have been identified in two unrelated children presenting Leigh syndrome (Alahmad et al., 2020). CP analysis from both subjects indicated severe impairment of CI assembly since abundance of CI-containing supercomplexes decreased prominently. In addition, intermediates composed by subunits of the Q-module, the Q/P-proximal module and the P-distal ND4 module accumulated, suggesting that stalling of CI assembly in these patients occurred at the intermediate stage (Guerrero-Castillo et al., 2017a). Comparably, mutations in the gene encoding TMEM70 have been identified in patients with deficiencies in complexes I and V (Jonckheere et al., 2011). The importance of TMEM70 in CV assembly has become clear from many mutations identified in human subjects (Cízková et al., 2008; Honzík et al., 2010), but recently in-depth analyses of complexome profiles delineated the role of TMEM70 in the formation of the F1-c intermediate (Sanchez-Caballero et al., 2020). In regard to CI however, accumulations of assembly intermediates corresponding to the Q/P-proximal module as well as to the P-distal ND5 module, evidenced by complexome profiling of TMEM70 knockout HAP1 cells, indicated stalling of CI assembly at the binding of P-distal ND4 module to these intermediates (Sanchez-Caballero et al., 2020).
CI deficiency may not exclusively originate from pathogenic mutations in structural components and assembly factors, but also from mutations in genes involved in the maintenance and repair of damaged subunits. In an exhaustive study including 29 families (Stenton et al., 2021), Stenton and co-workers identified mutations in DNAJC30 as a cause of LHON syndrome. Co-migration of DNAJC30 with OXPHOS supercomplexes was identified in the control but not in the patient complexome profile. The abundance of CI assembled into supercomplexes, slightly higher in the patient samples, excluded a role of DNAJC30 in the assembly of CI. Instead, turnover of subunits of the N-module of CI that are more frequently interchanged during optimal maintenance was slowed down, implicating DNAJC30 in a repair mechanism of CI still associated into the supercomplex. Interestingly, the ATP-dependent serine protease CLPP, whose mutations in humans are known to cause Perrault syndrome (Jenkinson et al., 2013), has also been implicated in CI repair by selectively replacing subunits of the N-module that are exposed to oxidative damage (Szczepanowska et al., 2020).
4.2 Complex III Deficiency
Proper assembly of CIII is not only important for generation of proton-motive force driven by the enzymatic oxidation of ubiquinol molecules produced by multiple dehydrogenases localized in the IMM, but also, at the core of respiratory supercomplexes, for structural stabilization of complexes I and IV. This phenomenon has been observed in several CIII-deficient patients where the activity and content of CI and CIV were concomitantly decreased (recently reviewed by Vercellino and Sazanov, 2021). In addition, a SILAC-duplexed CP study of cells lacking subunit cytochrome b due to a 4-bp deletion in MT-CYB and therefore unable to form holo-complex III, was consistent with a role of CIII-containing supercomplexes in CI maturation, since a late assembly intermediate containing assembly factor NDUFAF2 accumulated in cells with mutated cytochrome b (Protasoni et al., 2020). In contrast, mutations that partially affect the enzymatic activity of CIII, but not its assembly, do not compromise the stability and activity of CI and CIV. For instance, using a mouse model of CIII deficiency by deletion of TTC19, a protein involved in the removal of inhibitory proteolytic fragments of the Rieske iron-sulfur protein from the holo-complex, Bottani and coworkers showed that while the activity of CIII decreased in the tissues tested, the content and activity of CI and CIV were hardly altered (Bottani et al., 2017).
4.3 Complex IV Deficiency
Respiratory complex IV is a terminal oxidase that transfers electrons from soluble cytochrome c to oxygen, which also generates proton-motive force. Although no complexomes from patients with mutations in subunits or assembly factors of CIV have been described so far, multiple cell lines representing disease models of its deficiency have lately been characterized by CP. A complexome study of COX4I1/2 double knockout using HEK293 cells illustrated the interdependency of OXPHOS complexes in terms of their assembly and stability (Cunatova et al., 2021). The lack of subunits COX4I1/2 completely abolished formation of CIV. This is consistent with the role of this subunit forming an early assembly intermediate with COX5. In the absence of CIV, a pronounced decrease in the abundance of the CIV-containing supercomplexes have certainly been evidenced, whereas supercomplex I-III2 (S0) and complexes II, III and V remained unaffected. Notably, even if complex V forms large macromolecular assemblies that are independent of respirasomes, its interconnection to other OXPHOS complexes became apparent by the accumulation in COX4I1/2 double knockout cells not only of the F1 segment, but also of the early FO intermediate containing subunits ATP5MF, ATP5PB, ATP5ME and ATP5MG.
Similarly, Lobo-Jarne and co-workers analyzed the effects of the lack of CIV holoenzyme on the stability of respiratory complexes and respirasomes using trans-mitochondrial cybrids of human osteosarcoma 143B cells harboring nonsense mutations in mitochondrial genes MT-CO1 and MT-CO2 that resulted in truncated versions of subunits COX1 or COX2, respectively (Lobo-Jarne et al., 2020). In contrast to control 143B cells that, besides monomeric CIV, contained larger assemblies of this complex corresponding to IV2, III2-IV and respirasome I-III2-IV (S1), in cybrids with truncated COX1 and COX2 no CIV was formed but S0 remained stable. In the absence of COX1, the only detected assembly intermediate contained subunits COX4 and COX5, whereas the lack of signals of all other unassembled CIV subunits indicated their degradation. Conversely, cells lacking COX2 displayed, in addition to the heterodimer intermediate formed by COX4 and COX5, a larger intermediate containing also subunit COX1. Importantly, these authors described the formation of a heavier S0 containing several subunits and assembly factors of CIV, which has also been corroborated by bidirectional co-IP assays. The association of this partially assembled unconventional intermediate of CIV bound to S0 exemplified how alternative biogenesis of respirasomes could occur in disease-related conditions.
4.4 Combined Oxidative Phosphorylation Complexes Deficiency
The intricate assembly process of OXPHOS complexes requires not only expression and import of functional subunits as well as ancillary proteins produced by cytosolic ribosomes, but also proper expression and insertion of the mtDNA-encoded components. Therefore, diseases featured by energy metabolism alterations may also result from mitochondrial translation deficits. In mammals, the 13 mtDNA-encoded subunits of OXPHOS complexes are synthesized by the mitochondrial translation machinery (Hällberg and Larsson, 2014). This includes the small and large subunits of the mitoribosome [mt-SSU (28S) and mt-LSU (39S), respectively], a set of tRNAs and many regulatory proteins required for the proper functioning and fine tuning of this process (Kummer and Ban, 2021). Taking into account all the factors involved in mitochondrial translation, it is thus not surprising that an increasing number of genetic mutations causing multiple OXPHOS complex deficiencies due to impaired mitochondrial protein biosynthesis have been identified (Pearce et al., 2013).
Perrault syndrome is a rare genetic mitochondrial disorder with a relatively mild phenotype characterized by sensorineural deafness and ovarian dysgenesis. Genetic mutations in several genes encoding mitochondrial and peroxisomal proteins have been identified as direct causes for this disease, including CLPP, HSD17B4, LARS2, HARS2 and TWNK. Chatzispyrou et al. identified two unrelated female patients presenting Perrault syndrome symptoms where whole exome sequencing revealed no mutations in any of the enlisted genes causing this syndrome (Chatzispyrou et al., 2017). Instead, a predicted pathogenic mutation in ERAL1, a chaperone of the 12S mitochondrial rRNA was identified. Complexome profiles of the two patients evidenced a pronounced decrease in the abundance of the mt-SSU, while the mt-LSU remained unaffected; however, very slight changes in OXPHOS complexes intensities have been detected. This correlated with the milder phenotypes of Perrault syndrome patients, contrary to the clinical features observed when the mutations occurred in structural components of the mitoribosome, where multiple OXPHOS complex deficiencies manifested (Gardeitchik et al., 2018).
A study of two subjects with mutations in MRPS2, a structural component of the mt-SSU, has linked these mutations to combined OXPHOS complex defects with sensorineural hearing loss (Gardeitchik et al., 2018). The complexome profile of one of the subjects has been analyzed in parallel with fibroblasts from a previously studied patient with mutations in MRPS22 exhibiting similar symptoms (Smits et al., 2011). Complexomes of both patients revealed a marked decrease in the mt-SSU, while the large subunit remained unaffected. Decrease of mt-SSU abundance was also reflected in the patterns of OXPHOS complexes, which, in both patients, showed largely decreased abundance, particularly for CI and CIV. The abundance of monomeric CV remained unaltered in fibroblasts with mutations in MRPS2 or MRPS22. However, accumulations of assembly intermediates of CV, corresponding to the F1 segment and to FO early subassemblies that are not detected in control cell lines also indicated mitochondrial translation impairment in these patients.
A comparable decrease in the abundance of the mt-SSU has been found in a knockout cell line of METTL15, a N4-methylcytidine methyltransferase responsible for modification of the 12S mitochondrial rRNA (Van Haute et al., 2019). The absence of METTL15 impaired the translation of OXPHOS complex subunits encoded in the mtDNA and decreased the steady-state levels subunits of complexes I, III and IV. The DGU-based complexome profiles showed an interaction of METTL15 with mitoribosome components. In contrast to the aforementioned genetic defects where only the mt-SSU was affected, deletion of METTL15 resulted in a parallel decrease of the mt-LSU as well.
Shifting the energy metabolism from oxidative to glycolytic involves remodeling of biochemical pathways and structural rearranges of protein complexes inside and outside mitochondria. In this regard, Rho 0 (ρ0) cells, devoid of mtDNA, are incapable to respire due to the lack of essential subunits of OXPHOS complexes thus depend entirely on anaerobic ATP synthesis (King and Attardi, 1989). These cells represent an optimal system to study mtDNA-related defects and OXPHOS deficiencies. Beyond the impaired assembly of OXPHOS complexes I, III, IV and V, which did retain partially assembled intermediates formed by nuclear-encoded subunits, complexomes of ρ0 cells displayed unexpected modifications in migration patterns of a wide range of membrane-embedded and water-soluble protein complexes (Guerrero-Castillo et al., 2021b). Changes in expression of, for example, mitochondrial inner membrane carrier proteins, TCA cycle enzymes, mtDNA maintenance proteins and protein translocases localized to both mitochondrial membranes were described. Moreover, the absence of the mtDNA-encoded 12S and 16S rRNAs hindered the formation of mitochondrial ribosomal subunits, from which only smaller parts, composed exclusively of proteins, were preserved. This extensive re-organization of the mitochondrial complexome greatly illustrated the adaptability of mitochondria to self-adjust and compensate the absence of aerobic energy metabolism.
4.5 Pathogenic Mutations of Mitochondrial Proteins Affecting Indirectly the Assembly or Stability of Oxidative Phosphorylation Complexes
Reorganization of the mitochondrial energy metabolism can be triggered not only by genetic defects in OXPHOS complexes, assembly factors, mitochondrial translation and/or import machineries, but also by changes in the lipid environment of the inner mitochondrial membrane. This is particularly evident in patients with Barth syndrome, an inherited disorder of cardiolipin remodeling characterized by cardiomyopathy, neutropenia, 3-methylglutaconic aciduria and growth retardation (Neustein et al., 1979; Barth et al., 1983). This syndrome is caused by mutations in TAZ (Bione et al., 1996), encoding a mitochondrial acyltransferase called tafazzin, responsible for a CoA-independent exchange of (poly)unsaturated acyl chains from phosphatidyl-choline or phosphatidyl-ethanolamine to cardiolipin (Xu et al., 2006). The polyunsaturated acyl chains shape cardiolipin into a conical form, required for stabilization of membrane curvatures that optimize bioenergetics functions and are especially abundant at the inner mitochondrial membrane of multiple tissues, including heart. Complexome profiles of four Barth syndrome patients with altered saturation patterns of cardiolipin acyl chains revealed marked destabilization of several protein complexes involved in different mitochondrial processes. Significant disruption of supercomplexes containing 2–4 copies of CIV (I-III2-IV2–4), but not supercomplexes S0 and S1, has been found as well as partial destabilization of the complexes α-ketoglutarate dehydrogenase and branched-chain α-keto acid dehydrogenase (Chatzispyrou et al., 2018). Importantly, protein complexes involved in maintenance of mitochondrial morphology (e.g., MICOS) and in protein-driven mitochondrial apoptosis were substantially increased in patient fibroblasts. Overall, more than 200 proteins displayed altered migration patterns as a consequence of the membrane curvature destabilization, indicating a widespread reorganization of the mitochondrial complexome. These changes have not been limited to components of the inner mitochondrial membrane; but unexpectedly also spread among complexes located in the other three mitochondrial compartments.
In-depth understanding of remodeling of the mitochondrial complexome might be helpful to explain the links between the biochemical and clinical features observed in patients. For example, to delineate the molecular mechanisms underlying 3-methylglutaconic aciduria or neutropenia. These two clinical signs are also present in subjects with MEGCANN (3-methylglutaconic aciduria with cataracts, neurologic involvement and neutropenia), which is caused by mutations in CLPB, an ATP-dependent refoldase (dissagregase) of the mitochondrial intermembrane space. In a recent study comparing two subjects with bi-allelic mutations in CLPB and three subjects with mono-allelic mutations, different patterns of oligomeric CLPB and HAX1, whose genetics variants have been associated to neutropenia, were interpreted as a retention of HAX1 in inefficient CLPB oligomers in the patients (Wortmann et al., 2021). Besides this effect, aberrant migration patterns of multiple mitochondrial proteins were noticed. Comprehensively re-analysing these complexome datasets together with those from Barth syndrome patients and TMEM70 knockout cells, in which accumulation of 3-methylglutaconic acid is commonly observed, could potentially shed light into the causes and consequences of this type of aciduria.
5 Complexome Profiling: Perspectives
The implementation of CP analysis to help tackle the many open questions regarding PPIs is not limited to mitochondrial biology. CP is swiftly spreading to other scientific fields and already utilized to investigate protein complexes from other cell compartments and different species. However, there are still key limitations in the current setups, which open the door for further development and expanding its reach. In this final section we discussed not only the most recent advances in CP, but also its application to study protein complexes beyond the mitochondrion.
5.1 Improving Separation of Protein Complexes, Mass Spectrometry Quantification and Data Acquisition Time
The biochemical native separation techniques used to date for complexome analysis cover wide molecular mass ranges. Although each of these methods can easily be adapted to enhance the separation resolution of a certain mass window of interest, alternative strategies or modifications to the currently applied methods are envisioned, especially to improve separation resolution of large complexes.
Replacement of detergents by amphiphilic polymers (amphipols) could also help improve the analysis of membrane protein complexes. Amphipol-enclosed transmembrane complexes are more stable during separation than the ones solubilized with detergents (Althoff et al., 2011). Moreover, field flow fractionation (FFF) is a method to separate large (>1 nm) protein complexes. Proteins in solution are passed through a long narrow channel without a stationary phase and separated by their mobility upon a hydraulic field applied perpendicularly (Reschiglian and Moon, 2008), where small particles migrate faster than the bigger ones due to their smaller hydrodynamic sizes and higher diffusion coefficients. This technique has proven effective for analyzing ribosome profiles of the plant Nicothiana benthamiana with a similar resolution to the one of density gradients to distinguish ribosomal subunits from monosomes and free proteins (Pitkanen et al., 2014). FFF could thus represent an alternative separation technique in future CP setups.
Resolution of complexome profiles might be significantly improved by increasing the number of fractions in native-PAGE- and SEC-based CP setups. The more fractions the more MS data acquisition time is required however. Hence, a major limitation of CP is LC-MS time; especially if analysis of multiple experimental conditions is intended. Duplexing and triplexing with conventional SILAC/SILAM as well as multiplexing using isobaric labeling compounds are promising strategies for increasing the throughput of CP while reducing MS measuring times (Guerrero-Castillo et al., 2021a; Palenikova et al., 2021b). Incorporation of these strategies in CP studies also offers a higher level of proteomic quantification.
On the one hand, SILAC/SILAM-based quantifications enable more precise comparisons between samples or conditions than label-free quantification methods, and software tools have already been developed to facilitate data analysis (Palenikova et al., 2021a). SILAC is a suitable, affordable and compatible method with practically any cell culture condition. Accuracy and reproducibility are not hampered in this method since there are no precursor interferences as it occurs in isobaric tagging strategies (Chen et al., 2015). Other related variants, such as pulsed (p)SILAC and neutron-encoding (NeuCode) SILAC/SILAM, may also expand the possibilities this strategy can be used for in combination with CP. In pSILAC, the labeled essential aminoacids are added only for a brief period of time during cell growth to monitor differences in newly synthesized proteins and their half-life (Schwanhausser et al., 2009). As aforesaid, pSILAC and CP have proven useful to study the roles of mitochondrial protease ClpXP in turnover of N-module subunits of CI (Szczepanowska et al., 2020) and the chaperone DNAJC30 in the repair mechanism of impaired CI (Stenton et al., 2021). Furthermore, NeuCode SILAC/SILAM allows higher sensitivity and multiplexing with metabolic labeling using heavy stable isotopologues with additional neutrons (Hebert et al., 2013; Merrill et al., 2014; Overmyer et al., 2018). The subtle mass differences caused by neutron-binding energy variations are revealed under ultrahigh-resolution analysis that is only offered by Fourier-transform (FT)-MS or orbitrap analyzer systems (Overmyer et al., 2018). Thereby, implementation of this strategy in CP setups seems plausible.
Multiplexed CP (MCP), on the other hand, employs a set of isobaric labeling compounds, such as tandem mass tags that allows pooling a higher number of samples together after tryptic digestion. With commercially available tags, it is now possible to combine peptide solutions from up to 16 samples together before LC-MS/MS analyses. The quantitative advantage of using multiplexing isobaric compounds is that these molecules contain a reporter ion that is cleaved off during high collision dissociation together with precursor ion fragmentation. Therefore, quantifications of the reporter ions are performed at the MS2 level in the same spectrum used for peptide identification. Nevertheless, state-of-the-art equipment capable of MS3 or quantifications based on complement reporter ions are required to avoid quantification distortion due to interfering co-isolated ions (McAlister et al., 2014; Sonnett et al., 2018).
5.2 Combined Approaches for Selective Structural Analysis of Protein Complexes
Cross-linking mass spectrometry (XL-MS) has emerged as a powerful technology for interactome and structural biology studies in mitochondria (Schweppe et al., 2017; Liu et al., 2018; Ryl et al., 2020; Hevler et al., 2021b). Chemical cross-linkers react covalently with residues of proteins that are in close proximity to stabilize native PPIs for further characterization (Gomes and Gozzo, 2010; Mintseris and Gygi, 2020). After scrutiny of the identified intra- or inter-cross-linked peptides, a list of protein interactors and distance restraints is generated. This information can be used to guide computational modeling of single proteins and complex interfaces, most often in combination with data derived by classical structural methods, such as X-ray crystallography, nuclear magnetic resonance (NMR) or cryogenic electron microscopy (cryo-EM) (Herzog et al., 2012; Kornberg et al., 2015; O'Reilly and Rappsilber, 2018; Steigenberger et al., 2020).
Before structure analysis, determination of the optimal concentrations of cross-linkers and protein is strictly required to abate artifacts. To help circumvent this, a hybrid workflow termed in-gel cross-linking (IGX-MS) has recently been developed (Hevler et al., 2021a). In this method, proteins are separated by BN-PAGE and the gel spots of interest are cut and incubated with a cross-linker, digested and further analyzed by LC-MS/MS.
IGX-MS not only makes time-consuming experimental optimization steps (e.g., determining concentration of cross-linkers, optimal buffer system, etc.) nearly obsolete but also decreases the required protein amounts as well as undesired over-length cross-links. In contrast to classical in-solution XL-MS workflows, IGX-MS allows the differentiation of conformation- and interaction-specific distance restraints as it has been shown for various protein complexes, among them complexes I and V from bovine heart mitochondria. IGX-MS seems specially promising for future modeling studies that aim at characterizing co-occurring protein complexes with different stoichiometries or assembly states.
It has lately been demonstrated that CP and XL-MS are highly complementary and could provide valuable insights into the macromolecular organization of, for instance, the mitochondrial complexome. By combining CP and XL-MS data, Hevler and co-workers build a detailed atomic model of bovine CIV associated with dimeric apoptosis-inducing factor 1 (AIFM1) (Hevler et al., 2021b). CP analysis was indeed helpful to determine the stoichiometry of such a complex and validating this interaction. Combination of both approaches offers the possibility to cross-validate PPIs and to identify protein complexes that could be overlooked otherwise. Although actual scalability and throughput have not yet been addressed, combination of CP and XL-MS holds a great potential not only to better define the composition and states of protein complexes but also to help stabilize transient PPIs.
XL-MS has also been used to validate another CP-like setup that includes protein separation by SEC, quantitative MS analysis, cryo-electron microscopy (EM) and computational modeling. This structural proteomics approach allowed simultaneous description of the abundance, PPIs and structure profiling of 1/3 of the proteome of Chaetomium thermophilum (Kastritis et al., 2017). These results were integrated in a network map comprising 48 protein complexes and communities. This approach was also suitable to resolve the structure of the fatty acid synthase complex and its arrangements. Recent inclusion of image-processing workflows based on machine-learning methods opens the door to a much more robust data analysis, improved identification of PPIs, higher resolution in structure models and multi-scale molecular description of protein communities in situ (Kyrilis et al., 2021).
5.3 Investigating DNA-/RNA-Protein Complexes by Complexome Profiling
CP has not been popular for the analysis of DNA-interacting protein complexes. These have been addressed with more targeted approaches like chromatin immunoprecipitation (ChIP) (Das et al., 2004), electrophoretic mobility shift assay (EMSA) (Hellman and Fried, 2007) and a variety of Co-IP techniques (Sahr and Buchrieser, 2013), which have been extensively reviewed by Ferraz et al. (2021). Although both genomic and mitochondrial DNA molecules that interact with proteins are too large to enter standard native gels and do not exist in populations separable by size, like RNPs, CP could be useful for decomposing DNA-associated protein complexes. For example, Munawar and co-workers studied the complexome of a chromatin-enriched fraction separated by BN-PAGE after enzymatic digestion of DNA (Munawar et al., 2015). The resulting profile contained >50% of known chromatin complexes represented by at least half their subunits and contained information about assemblies for some of the complexes, such as PRC2 and NuRD. This illustrated the potential of CP for massive characterization of DNA-associated complexes. This fact is also corroborated by appearance of a variety of DNA-related complexes in existing CP datasets, such as histones, DNA repairing, replication and transcription factors.
Two research groups aiming to characterize the protein composition of chromatin complexes have made use of separation of protein complexes by size and LC-MS/MS analysis followed by correlation profiling. One utilized sucrose density gradients for separation of complexes containing chromatin-associated protein Wdr82 (Lee and Skalnik, 2013). The other used BN-PAGE to separate different populations of a chromatin remodeling complex NuRD derived from mouse embryonic stem cells. This approach led, among other things, to discovery of a new subunit of NuRD-associated protein, Wdr5 (Bode et al., 2016). Although these works were performed using affinity-purified complexes, they illustrated the potential of CP for characterization of chromatin-associated protein complexes as both allowed identification of different populations and assemblies in native state in contrast to pull-down techniques.
Conversely, CP variants to study RNA-protein complexes have been developed. Yet, these approaches have been particularly used for studying ribonucleoproteins (RNPs). Applications of RNP complexomics are generously reviewed in Gerovac et al. (2021). Complexomics studies based on protein co-migration profiles obtained from density gradients are commonly used to investigate ribosomal protein components and associated factors (Yu et al., 2005; Aviner et al., 2017; Van Haute et al., 2019; Palenikova et al., 2021a). Besides ribosomes, other RNPs are likely to have specific sedimentation patterns on density gradients due to their complex composition (Caudron-Herger et al., 2019). The simplicity and gentleness of fractionation by density gradient allows parallel identification of both proteins and RNA molecules by LC-MS/MS and RNA-seq, respectively. This setup has been established in a robust RNP complexomics method called Grad-seq that enabled global classification of stable RNA-protein complexes as well as identification of novel RNPs (Smirnov et al., 2016; Hor et al., 2020).
SEC is another method that has successfully been used for separation of RNPs from complex mixtures (Yoshikawa et al., 2018; Mallam et al., 2019). It provides fast and reproducible automated separation of particles with an impressive upper limit size of ∼30–40 MDa. SEC has been applied to specifically separate human 80S monosomes (∼4.3 MDa) and n-polysomes from HeLa cells (Yoshikawa et al., 2018). Moreover, hierarchical clustering revealed complexes involved in translation associated with ribosomal subunits, which makes this method a good alternative to the time-consuming and less specific DGU.
Classic CP has also been useful for dissecting protein composition of mitoribosomes as they enter regular BN-gels (Wessels et al., 2013; Chatzispyrou et al., 2017; Gardeitchik et al., 2018; Van Strien et al., 2019; Cunatova et al., 2021). For larger ribosomal species however, LP-native gels can be used (Strecker et al., 2010). As aforementioned, LP-BN-PAGE-based CP helped identify a RNP complex containing the mt-SSU associated with several translation and transcription factors at ∼5.7 MDa (Rugen et al., 2019). LP-gels could thus be very useful for further dissection of larger RNPs or those that are associated with many other factors as it has been commonly reported for translation and RNA-processing machineries (Acestor et al., 2009; Hocine et al., 2010; Blombach et al., 2011; Doetsch et al., 2011; Sprink et al., 2016; Duss et al., 2019; Gopalakrishna et al., 2019; Hilander et al., 2021).
A potential issue for characterization of RNA-protein complexes may result from the poor stability of RNA molecules, which are also prone to degradation by RNases. In this case, RNA-protein interactions would be lost and leading to identification of not only sub-assemblies but also artifacts. To avoid this, prior to CP, separation of these complexes should be thus performed using RNase-free materials and RNase inhibitors to keep stable both the protein and RNA interactors (López De Heredia and Jansen, 2004; Mili and Steitz, 2004; Wang et al., 2016; Aarum et al., 2020). Alternatively, to stabilize weak interactions of fragile RNA-protein complexes during separation, cross-linking strategies can be implemented. Treatment with chemicals (Rugen et al., 2019; Herrmannova et al., 2020; Patton et al., 2020), such as formaldehyde or UV irradiation (Trendel et al., 2019; Urdaneta et al., 2019; Urdaneta and Beckmann, 2020; van Esveld and Spelbrink, 2021) create short cross-links between a nucleic acid and close-contact protein interactors. In contrast, different types of enzymatic digestion of RNA can be used purposely to affect structure and stability of RNA-protein complexes to build additional levels of information for CP analysis, which can enable identification of previously unknown RNA-binding proteins (Caudron-Herger et al., 2019; Mallam et al., 2019; Gerovac et al., 2021).
5.4 Filling the Void: Complexome Profiling as a Workflow to Identify Protein-Protein Interactions and Improve Functional Annotation of Protein Complexes Across Species
Many species still lack functional annotation for a large portion of their genomes despite their significance as pathogens or drivers of major ecological processes. This is often owed to poor conservation, missing orthologues and the additional pitfall that the most unique biology that is not found in standard models, remains enigmatic. Systematic knowledge of PPIs can allow us to place these proteins into pathways, compartments and interaction networks. CP represents a clear choice to help narrow this gap and has the added benefit that no genetic intervention is required, which allows immediate assessment of species without a well-developed toolkit.
For example, Hillier et al. applied CP to study schizont stages of three different Plasmodium species. Integration with machine-learning led to identification of ∼20,000 putative PPIs arranged into ∼600 protein clusters, which were used to map the interactome (Hillier et al., 2019). Next to creating a valuable resource for Plasmodium research, the authors shed light on the interaction network of a group of parasite-specific transcription factors as well as identifying novel parasite-specific interactors. The latter could serve as targets for drug development due to their divergence from host biology. These data are sorely needed in the face of emerging drug resistance that is threatening our attempts to eradicate malaria (Haldar et al., 2018).
Prokaryotes also deviate considerably from the standard protein complexes found in eukaryotes. These organisms cover virtually every ecological niche on earth and have adapted their protein complexes accordingly. This is indeed exemplified in energy metabolic pathways, where prokaryotes utilize a wide variety of electron donors and acceptors, each requiring their own adapted set of classical and alternative respiratory complexes (Poole and Cook, 2000; Kaila and Wikstrom, 2021), most of which remain undescribed. While genomic analyses led to the hypothetical annotation of many homologous proteins involved in prokaryotic respiration, this approach is ill suited to the identification of novel components and interactors. CP complements this shortcoming and stands up as a good approach for characterizing PPIs in these organisms. This was illustrated beautifully by de Almeida and co-workers who investigated the electron transport system of the anaerobic ammonium-oxidizing bacterium Kuenenia stuttgartiensis (de Almeida et al., 2016). These bacteria oxidize ammonium with nitrite to produce N2 gas. It has been estimated that this process account for ∼50% of all N2 emitted into the atmosphere (Devol, 2015). Although this enormous ecological significance, the underlying respiratory complexes were largely theoretical (Kartal et al., 2013). By applying CP, these authors did not only confirm the presence of nearly all predicted respiratory complexes but also found evidence for novel protein complexes. One striking finding was that this bacterium appears to assemble a set of respiratory complexes that utilize Na+ besides the expected proton-pumping complexes, which suggested a role for a sodium-motive force to drive, for instance, the unfavorable reduction of ferredoxin. A set of three Rieske/cytb complexes has also been found that shed light on how these bacteria can capture energy from hydrazine oxidation, challenging a previous hypothesis suggesting a quinone-reducing enzyme (Kartal et al., 2013). Similarly, CP has been used to analyze the nitrite-dependent methanotroph Methylomirabilis lanthanidiphila (Versantvoort et al., 2019). This species may potentially contribute to reduction of methane and nitrous oxide; i.e., potent greenhouse gases, by oxidizing methane to CO2 and reducing nitrate without producing nitrous oxide. To explain this process, a metabolic model based on genomic information had previously been proposed (Ettwig et al., 2010). Versantvoort and co-workers could confirm this model by identifying all proposed underlying protein complexes. Interestingly, three protein complexes thought to be involved in nitrate reduction were also identified, suggesting metabolic capabilities which had not yet been described in Methylomirabilis species. On the other hand, application of CP to study the marine sulfate reducer Desulfobacula toluolica, also led to detection of not only the expected membrane complexes but also new assemblies, unexpected associations, multimeric states, redox complexes involved in Na+-based bioenergetics and indications of respiratory supercomplexes (Wöhlbrand et al., 2016).
To shed more light on the biogenesis of thylakoid membrane complexes, CP was used to study auxiliary factors involved in the assembly of photosystem II (PSII) from the single-cell green alga Chlamydomonas reinhardtii (Spaniol et al., 2021). Authors demonstrated that the homolog for low PSII accumulation 2 (LPA2) protein is critical for PSII assembly and further supercomplex formation. Besides, novel interactors were found and possibly involved in regulation of PSII. This study exemplifies the great potential CP has for further studying photosynthesis-related protein complexes in algae, cyanobacteria and plants.
Finally, even in microorganisms that are not currently obtainable in axenic cultures, such as the archaeon Candidatus Methanoperedens, CP has proven invaluable (Berger et al., 2021). In this study, Berger et al. confirmed the presence of several predicted membrane-bound respiratory chain complexes as well as novel associations. Taken together these results demonstrated the suitability of CP to explore species previously inaccessible to most protein research. Overall, it is clear that CP complements (meta)genomic studies perfectly, by identifying in one experiment how these gene products assemble to form protein complexes, what subset of genome is really expressed and consequently update our understanding of metabolic capabilities and how these organisms achieve them.
5.5 Exploring Effects of Pathogenic Mutations in Non-Mitochondrial Proteins by Complexome Profiling
Although to date the vast majority of CP analyses in patients have been focused on disorders affecting mitochondria, its application is broadening nonetheless. The effects of pathogenic missense mutations in structural subunits of the V-type ATPase on the stability of the holoenzyme have been characterized by this method (Van Damme et al., 2017). By whole exome sequencing, predicted pathogenic mutations have been identified in genes ATP6V1E1 and ATP6V1A, encoding subunits E and A from the V1 segment of the V-type ATPase, respectively. Primary skin fibroblasts from three patients, one with a mutation in E1 and two with mutations in subunit A were compared to controls and to an earlier described subject presenting cutis laxa with a mutation in subunit ATP6V0A2 of the VO segment. To further investigate this, enriched Golgi preparations have been analyzed by BN-PAGE-based CP. Migration patterns of healthy controls revealed three assembly intermediates of the V1 segment, the VO segment and a highly abundant, fully assembled V1VO-ATPase complex. In contrast, in patient fibroblasts with mutations in subunits of the V1 segment, only the VO intermediate and traces of the holo-complex have been detected. Comparably, the mutation in ATP6V0A2 decreased the abundance of intermediates of both domains as well as the holo-complex. The destabilization of the fully assembled V1VO-ATPase evidenced by CP contributed to detailing the molecular consequences and pathogenicity of these mutations.
In a recent CP study of peroxisome-deficient cells using peroxins deletion strains of S. cerevisiae as a model system and fibroblasts from a Zellweger syndrome patient, Nübel and co-workers linked accumulation of miss-targeted peroxisomal membrane proteins into mitochondria to metabolic and morphological abnormalities of this organelle (Nuebel et al., 2021). Complexome profiles of yeast strains lacking peroxisomes showed formation of the docking subassembly of the peroxisomal importomer complex in mitochondria, explaining also the incorporation of peroxisomal matrix proteins to this subcellular compartment. It would be interesting to also explore the complexome profiles of fibroblasts from patients with Zellweger spectrum disorder to investigate in detail the rearrangements of mitochondrial complexes caused by peroxisomes deficit in humans.
These two studies exemplified how CP offers a straightforward option to shed light on the pathophysiological roles of critical proteins located in cell compartments other than mitochondria and the specific consequences of impairments in protein complexes and PPIs resulting from genetic and/or biochemical defects. CP thus holds a great potential to help not only unveil the molecular mechanisms of diseases, but also develop novel treatments.
6 Concluding Remarks
In a very short time, CP has revolutionized the way we study mitochondrial protein complexes and significantly increased the amount of data collected from single unbiased experiments. Research findings obtained by this methodology include, for instance, identification of novel OXPHOS-related subunits and assembly factors; stepwise description of intricate assembly pathways; molecular evidence on the role of accessory subunits and existence of tissue-specific subunit isoforms of CI; better insight on the formation of respirasomes and the supramolecular organization of the MICOS/MIB complexes as well as the unique composition of OXPHOS complexes and mitoribosomes in different species. Notably, CP has also been invaluable to better understand the effects of mutations found in mitochondrial genes from patients as well as in disease models. Recent developments in duplexing and multiplexing CP experiments not only enable accurate proteomic quantifications, but also greatly increase the feasibility of including this approach in future clinical diagnosis. We have also discussed innovative strategies exemplifying how CP easily complements other state-of-the-art methods such as XL-MS, RNAseq and cryo-EM. However, there is still a lot of room for improving this method and keep expanding its reach. It is projected that CP will further evolve with the advent of the next generation of LC-MS instruments and proteomic strategies as well as upcoming complexome data analysis software. Finally, we also encourage the scientific community to apply CP for studying protein complexes from cell compartments other than the mitochondrion and share publicly the resultant complexome datasets through CEDAR.
Author Contributions
AC-O and SG-C contributed to the conception and design of this review. AC-O and SG-C wrote the first draft of the manuscript with contributions of AP, FE and JH in specific subsections. SG-C generated Table 1. AC-O and AP generated the set of figures. All authors contributed to discussion, literature revision, edition of figures as well as manuscript revision. All authors approved the submitted version.
Funding
AC-O is supported by the Netherlands Organization for Health Research and Development (ZonMW TOP-Grant 91217009); AP is supported by the European Union’s Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement 721757; FE is supported by the Netherlands Organization for Scientific Research (NWO-VIDI 864.13.009); JH is supported by the NWO funding the Netherlands Proteomics Centre through the X-omics Road Map program (Project 184.034.019), the TOP project 714.017.004 and the EU Horizon 2020 program Epic-XS (Project 823839). This work was also supported by the Else Kröner-Fresenius-Stiftung grant 2019_A135 to SG-C.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
Authors thank Prof. Ulrich Brandt for carefully revising this manuscript and helpful suggestions.
References
Aarum, J., Cabrera, C. P., Jones, T. A., Rajendran, S., Adiutori, R., Giovannoni, G., et al. (2020). Enzymatic Degradation of RNA Causes Widespread Protein Aggregation in Cell and Tissue Lysates. EMBO Rep. 21. doi:10.15252/embr.201949585
Acestor, N., Panigrahi, A. K., Carnes, J., Zíková, A., and Stuart, K. D. (2009). The MRB1 Complex Functions in Kinetoplastid RNA Processing. RNA 15, 277–286. doi:10.1261/RNA.1353209
Adams, J. (2008). The Proteome: Discovering the Structure and Function of Proteins. Nat. Education 1, 6.
Alahmad, A., Nasca, A., Heidler, J., Thompson, K., Oláhová, M., Legati, A., et al. (2020). Bi‐allelic Pathogenic Variants in NDUFC2 Cause Early‐onset Leigh Syndrome and Stalled Biogenesis of Complex I. EMBO Mol. Med. 12, e12619. doi:10.15252/emmm.202012619
Alston, C. L., Compton, A. G., Formosa, L. E., Strecker, V., Oláhová, M., Haack, T. B., et al. (2016). Biallelic Mutations in TMEM126B Cause Severe Complex I Deficiency with a Variable Clinical Phenotype. Am. J. Hum. Genet. 99, 217–227. doi:10.1016/j.ajhg.2016.05.021
Alston, C. L., Heidler, J., Dibley, M. G., Kremer, L. S., Taylor, L. S., Fratter, C., et al. (2018). Bi-allelic Mutations in NDUFA6 Establish its Role in Early-Onset Isolated Mitochondrial Complex I Deficiency. Am. J. Hum. Genet. 103, 592–601. doi:10.1016/j.ajhg.2018.08.013
Alston, C. L., Veling, M. T., Heidler, J., Taylor, L. S., Alaimo, J. T., Sung, A. Y., et al. (2020). Pathogenic Bi-allelic Mutations in NDUFAF8 Cause Leigh Syndrome with an Isolated Complex I Deficiency. Am. J. Hum. Genet. 106, 92–101. doi:10.1016/j.ajhg.2019.12.001
Althoff, T., Mills, D. J., Popot, J.-L., and Kühlbrandt, W. (2011). Arrangement of Electron Transport Chain Components in Bovine Mitochondrial Supercomplex I1III2IV1. EMBO J. 30, 4652–4664. doi:10.1038/emboj.2011.324
Anand, R., Kondadi, A. K., Meisterknecht, J., Golombek, M., Nortmann, O., Riedel, J., et al. (2020). MIC26 and MIC27 Cooperate to Regulate Cardiolipin Levels and the Landscape of OXPHOS Complexes. Life Sci. Alliance 3, e202000711. doi:10.26508/lsa.202000711
Anand, R., Strecker, V., Urbach, J., Wittig, I., and Reichert, A. S. (2016). Mic13 Is Essential for Formation of Crista Junctions in Mammalian Cells. PLOS ONE 11, e0160258. doi:10.1371/journal.pone.0160258
Andersen, J. S., Wilkinson, C. J., Mayor, T., Mortensen, P., Nigg, E. A., and Mann, M. (2003). Proteomic Characterization of the Human Centrosome by Protein Correlation Profiling. Nature 426, 570–574. doi:10.1038/nature02166
Andrews, B., Carroll, J., Ding, S., Fearnley, I. M., and Walker, J. E. (2013). Assembly Factors for the Membrane Arm of Human Complex I. Proc. Natl. Acad. Sci. 110, 18934–18939. doi:10.1073/pnas.1319247110
Angerer, H., Radermacher, M., Ma kowska, M., Steger, M., Zwicker, K., Heide, H., et al. (2014). The LYR Protein Subunit NB4M/NDUFA6 of Mitochondrial Complex I Anchors an Acyl Carrier Protein and Is Essential for Catalytic Activity. Proc. Natl. Acad. Sci. 111, 5207–5212. doi:10.1073/pnas.1322438111
Aurrecoechea, C., Brestelli, J., Brunk, B. P., Dommer, J., Fischer, S., Gajria, B., et al. (2009). PlasmoDB: a Functional Genomic Database for Malaria Parasites. Nucleic Acids Res. 37, D539–D543. doi:10.1093/nar/gkn814
Aviner, R., Hofmann, S., Elman, T., Shenoy, A., Geiger, T., Elkon, R., et al. (2017). Proteomic Analysis of Polyribosomes Identifies Splicing Factors as Potential Regulators of Translation during Mitosis. Nucleic Acids Res. 45, 5945–5957. doi:10.1093/nar/gkx326
Balsa, E., Marco, R., Perales-Clemente, E., Szklarczyk, R., Calvo, E., Landázuri, M. O., et al. (2012). NDUFA4 Is a Subunit of Complex IV of the Mammalian Electron Transport Chain. Cel Metab. 16, 378–386. doi:10.1016/j.cmet.2012.07.015
Barth, P. G., Scholte, H. R., Berden, J. A., Van Der Klei-Van Moorsel, J. M., Luyt-Houwen, I. E. M., Van'T Veer-Korthof, E. T., et al. (1983). An X-Linked Mitochondrial Disease Affecting Cardiac Muscle, Skeletal Muscle and Neutrophil Leucocytes. J. Neurol. Sci. 62, 327–355. doi:10.1016/0022-510x(83)90209-5
Berger, S., Cabrera-Orefice, A., Jetten, M. S. M., Brandt, U., and Welte, C. U. (2021). Investigation of central Energy Metabolism-Related Protein Complexes of ANME-2d Methanotrophic Archaea by Complexome Profiling. Biochim. Biophys. Acta (Bba) - Bioenerg. 1862, 148308. doi:10.1016/j.bbabio.2020.148308
Bione, S., D'adamo, P., Maestrini, E., Gedeon, A. K., Bolhuis, P. A., and Toniolo, D. (1996). A Novel X-Linked Gene, G4.5. Is Responsible for Barth Syndrome. Nat. Genet. 12, 385–389. doi:10.1038/ng0496-385
Blombach, F., Brouns, S. J. J., and van der Oost, J. (2011). Assembling the Archaeal Ribosome: Roles for Translation-Factor-Related GTPases. Biochem. Soc. Trans. 39, 45–50. doi:10.1042/BST0390045
Bludau, I., and Aebersold, R. (2020). Proteomic and Interactomic Insights into the Molecular Basis of Cell Functional Diversity. Nat. Rev. Mol. Cel Biol 21, 327–340. doi:10.1038/s41580-020-0231-2
Bludau, I. (2021). Discovery-Versus Hypothesis-Driven Detection of Protein-Protein Interactions and Complexes. Ijms 22, 4450. doi:10.3390/ijms22094450
Bludau, I., Heusel, M., Frank, M., Rosenberger, G., Hafen, R., Banaei-Esfahani, A., et al. (2020). Complex-centric Proteome Profiling by SEC-SWATH-MS for the Parallel Detection of Hundreds of Protein Complexes. Nat. Protoc. 15, 2341–2386. doi:10.1038/s41596-020-0332-6
Bock, T., Türk, C., Aravamudhan, S., Keufgens, L., Bloch, W., Rozsivalova, D. H., et al. (2021). PERM1 Interacts with the MICOS-MIB Complex to Connect the Mitochondria and Sarcolemma via Ankyrin B. Nat. Commun. 12, 4900. doi:10.1038/s41467-021-25185-3
Bode, D., Yu, L., Tate, P., Pardo, M., and Choudhary, J. (2016). Characterization of Two Distinct Nucleosome Remodeling and Deacetylase (NuRD) Complex Assemblies in Embryonic Stem Cells. Mol. Cell Proteomics 15, 878–891. doi:10.1074/mcp.M115.053207
Bottani, E., Cerutti, R., Harbour, M. E., Ravaglia, S., Dogan, S. A., Giordano, C., et al. (2017). TTC19 Plays a Husbandry Role on UQCRFS1 Turnover in the Biogenesis of Mitochondrial Respiratory Complex III. Mol. Cel 67, 96–105. doi:10.1016/j.molcel.2017.06.001
Bridges, H. R., Mohammed, K., Harbour, M. E., and Hirst, J. (2017). Subunit NDUFV3 Is Present in Two Distinct Isoforms in Mammalian Complex I. Biochim. Biophys. Acta (Bba) - Bioenerg. 1858, 197–207. doi:10.1016/j.bbabio.2016.12.001
Burgess, R. R. (2018). A Brief Practical Review of Size Exclusion Chromatography: Rules of Thumb, Limitations, and Troubleshooting. Protein Expr. Purif. 150, 81–85. doi:10.1016/j.pep.2018.05.007
Burley, S. K., Bhikadiya, C., Bi, C., Bittrich, S., Chen, L., Crichlow, G. V., et al. (2021). RCSB Protein Data Bank: Powerful New Tools for Exploring 3D Structures of Biological Macromolecules for Basic and Applied Research and Education in Fundamental Biology, Biomedicine, Biotechnology, Bioengineering and Energy Sciences. Nucleic Acids Res. 49, D437–D451. doi:10.1093/nar/gkaa1038
Calvo, E., Cogliati, S., Hernansanz-Agustín, P., Loureiro-López, M., Guarás, A., Casuso, R. A., et al. (2020). Functional Role of Respiratory Supercomplexes in Mice: SCAF1 Relevance and Segmentation of the Q Pool. Sci. Adv. 6, eaba7509. doi:10.1126/sciadv.aba7509
Caudron-Herger, M., Rusin, S. F., Adamo, M. E., Seiler, J., Schmid, V. K., Barreau, E., et al. (2019). R-DeeP: Proteome-wide and Quantitative Identification of RNA-dependent Proteins by Density Gradient Ultracentrifugation. Mol. Cel 75, 184–199. doi:10.1016/j.molcel.2019.04.018
Ceulemans, H., Beke, L., and Bollen, M. (2006). Approaches to Defining the Ancestral Eukaryotic Protein Complexome. Bioessays 28, 316–324. doi:10.1002/bies.20373
Chalkley, R. J., Baker, P. R., Huang, L., Hansen, K. C., Allen, N. P., Rexach, M., et al. (2005). Comprehensive Analysis of a Multidimensional Liquid Chromatography Mass Spectrometry Dataset Acquired on a Quadrupole Selecting, Quadrupole Collision Cell, Time-Of-Flight Mass Spectrometer. Mol. Cell Proteomics 4, 1194–1204. doi:10.1074/mcp.D500002-MCP200
Chatzispyrou, I. A., Alders, M., Guerrero-Castillo, S., Zapata Perez, R., Haagmans, M. A., Mouchiroud, L., et al. (2017). A Homozygous Missense Mutation in ERAL1, Encoding a Mitochondrial rRNA Chaperone, Causes Perrault Syndrome. Hum. Mol. Genet. 26, 2541–2550. doi:10.1093/hmg/ddx152
Chatzispyrou, I. A., Guerrero-Castillo, S., Held, N. M., Ruiter, J. P. N., Denis, S. W., Ijlst, L., et al. (2018). Barth Syndrome Cells Display Widespread Remodeling of Mitochondrial Complexes without Affecting Metabolic Flux Distribution. Biochim. Biophys. Acta (Bba) - Mol. Basis Dis. 1864, 3650–3658. doi:10.1016/j.bbadis.2018.08.041
Chen, X., Wei, S., Ji, Y., Guo, X., and Yang, F. (2015). Quantitative Proteomics Using SILAC: Principles, Applications, and Developments. Proteomics 15, 3175–3192. doi:10.1002/pmic.201500108
Čížková, A., Stránecký, V., Mayr, J. A., Tesařová, M., Havlíčková, V., Paul, J., et al. (2008). TMEM70 Mutations Cause Isolated ATP Synthase Deficiency and Neonatal Mitochondrial Encephalocardiomyopathy. Nat. Genet. 40, 1288–1290. doi:10.1038/ng.246
Cogliati, S., Calvo, E., Loureiro, M., Guaras, A. M., Nieto-Arellano, R., Garcia-Poyatos, C., et al. (2016a). Mechanism of Super-assembly of Respiratory Complexes III and IV. Nature 539, 579–582. doi:10.1038/nature20157
Cogliati, S., Enriquez, J. A., and Scorrano, L. (2016b). Mitochondrial Cristae: Where Beauty Meets Functionality. Trends Biochem. Sci. 41, 261–273. doi:10.1016/j.tibs.2016.01.001
Cole, J. L., Lary, J. W., P. Moody, T., and Laue, T. M. (2008). Analytical Ultracentrifugation: Sedimentation Velocity and Sedimentation Equilibrium. Methods Cel Biol 84, 143–179. doi:10.1016/S0091-679X(07)84006-4
Connelly, K. E., Hedrick, V., Paschoal Sobreira, T. J., Dykhuizen, E. C., and Aryal, U. K. (2018). Analysis of Human Nuclear Protein Complexes by Quantitative Mass Spectrometry Profiling. Proteomics 18, 1700427. doi:10.1002/pmic.201700427
Coordinators, N. R. (2016). Database Resources of the National Center for Biotechnology Information. Nucleic Acids Res. 44, D7–D19. doi:10.1093/nar/gkv1290
Cox, J., Hein, M. Y., Luber, C. A., Paron, I., Nagaraj, N., and Mann, M. (2014). Accurate Proteome-wide Label-free Quantification by Delayed Normalization and Maximal Peptide Ratio Extraction, Termed MaxLFQ. Mol. Cell Proteomics 13, 2513–2526. doi:10.1074/mcp.M113.031591
Čunátová, K., Reguera, D. P., Vrbacký, M., Fernández-Vizarra, E., Ding, S., Fearnley, I. M., et al. (2021). Loss of COX4I1 Leads to Combined Respiratory Chain Deficiency and Impaired Mitochondrial Protein Synthesis. Cells 10, 369. doi:10.3390/cells10020369
Das, P. M., Ramachandran, K., Vanwert, J., and Singal, R. (2004). “Chromatin Immunoprecipitation Assay,” in BioTechniques. Biotechniques). doi:10.2144/04376rv01
De Almeida, N. M., Wessels, H. J. C. T., De Graaf, R. M., Ferousi, C., Jetten, M. S. M., Keltjens, J. T., et al. (2016). Membrane-bound Electron Transport Systems of an Anammox Bacterium: A Complexome Analysis. Biochim. Biophys. Acta (Bba) - Bioenerg. 1857, 1694–1704. doi:10.1016/j.bbabio.2016.07.006
De Las Rivas, J., and Fontanillo, C. (2010). Protein-protein Interactions Essentials: Key Concepts to Building and Analyzing Interactome Networks. Plos Comput. Biol. 6, e1000807. doi:10.1371/journal.pcbi.1000807
Deshaies, R. J., Seol, J. H., Mcdonald, W. H., Cope, G., Lyapina, S., Shevchenko, A., et al. (2002). Charting the Protein Complexome in Yeast by Mass Spectrometry. Mol. Cell Proteomics 1, 3–10. doi:10.1074/mcp.r100001-mcp200
Devol, A. H. (2015). Denitrification, Anammox, and N2 Production in Marine Sediments. Annu. Rev. Mar. Sci. 7, 403–423. doi:10.1146/annurev-marine-010213-135040
Doetsch, M., Schroeder, R., and Fürtig, B. (2011). Transient RNA-Protein Interactions in RNA Folding. Febs J. 278, 1634–1642. doi:10.1111/J.1742-4658.2011.08094.X
Duss, O., Stepanyuk, G. A., Puglisi, J. D., and Williamson, J. R. (2019). Transient Protein-RNA Interactions Guide Nascent Ribosomal RNA Folding. Cell 179, 1357–1369. e1316. doi:10.1016/J.CELL.2019.10.035
Ettwig, K. F., Butler, M. K., Le Paslier, D., Pelletier, E., Mangenot, S., Kuypers, M. M. M., et al. (2010). Nitrite-driven Anaerobic Methane Oxidation by Oxygenic Bacteria. Nature 464, 543–548. doi:10.1038/nature08883
Eubel, H., Braun, H.-P., and Millar, A. (2005). Blue-native PAGE in Plants: a Tool in Analysis of Protein-Protein Interactions. Plant Methods 1, 11. doi:10.1186/1746-4811-1-11
Evers, F., Cabrera-Orefice, A., Elurbe, D. M., Kea-Te Lindert, M., Boltryk, S. D., Voss, T. S., et al. (2021). Composition and Stage Dynamics of Mitochondrial Complexes in Plasmodium Falciparum. Nat. Commun. 12, 3820. doi:10.1038/s41467-021-23919-x
Eydt, K., Davis, K. M., Behrendt, C., Wittig, I., and Reichert, A. S. (2017). Cristae Architecture Is Determined by an Interplay of the MICOS Complex and the F1Fo ATP Synthase via Mic27 and Mic10. Microb. Cel 4, 259–272. doi:10.15698/mic2017.08.585
Fabre, B., Lambour, T., Bouyssié, D., Menneteau, T., Monsarrat, B., Burlet-Schiltz, O., et al. (2014). Comparison of Label-free Quantification Methods for the Determination of Protein Complexes Subunits Stoichiometry. EuPA Open Proteomics 4, 82–86. doi:10.1016/j.euprot.2014.06.001
Fassone, E., and Rahman, S. (2012). Complex I Deficiency: Clinical Features, Biochemistry and Molecular Genetics. J. Med. Genet. 49, 578–590. doi:10.1136/jmedgenet-2012-101159
Fernández-Vizarra, E., López-Calcerrada, S., Formosa, L. E., Pérez-Pérez, R., Ding, S., Fearnley, I. M., et al. (2021). SILAC-based Complexome Profiling Dissects the Structural Organization of the Human Respiratory Supercomplexes in SCAFIKO Cells. Biochim. Biophys. Acta (Bba) - Bioenerg. 1862, 148414. doi:10.1016/j.bbabio.2021.148414
Ferraz, R. a. C., Lopes, A. L. G., Da Silva, J. a. F., Moreira, D. F. V., Ferreira, M. J. N., and De Almeida Coimbra, S. V. (2021). “DNA–protein Interaction Studies: a Historical and Comparative Analysis,” in Plant Methods. Plant Methods).
Floyd, B. J., Wilkerson, E. M., Veling, M. T., Minogue, C. E., Xia, C., Beebe, E. T., et al. (2016). Mitochondrial Protein Interaction Mapping Identifies Regulators of Respiratory Chain Function. Mol. Cel 63, 621–632. doi:10.1016/j.molcel.2016.06.033
Formosa, L. E., Muellner-Wong, L., Reljic, B., Sharpe, A. J., Jackson, T. D., Beilharz, T. H., et al. (2020). Dissecting the Roles of Mitochondrial Complex I Intermediate Assembly Complex Factors in the Biogenesis of Complex I. Cel Rep. 31, 107541. doi:10.1016/j.celrep.2020.107541
Foster, L. J., De Hoog, C. L., Zhang, Y., Zhang, Y., Xie, X., Mootha, V. K., et al. (2006). A Mammalian Organelle Map by Protein Correlation Profiling. Cell 125, 187–199. doi:10.1016/j.cell.2006.03.022
Fuhrmann, D. C., Wittig, I., Dröse, S., Schmid, T., Dehne, N., and Brüne, B. (2018). Degradation of the Mitochondrial Complex I Assembly Factor TMEM126B under Chronic Hypoxia. Cell. Mol. Life Sci. 75, 3051–3067. doi:10.1007/s00018-018-2779-y
Galemou Yoga, E., Parey, K., Djurabekova, A., Haapanen, O., Siegmund, K., Zwicker, K., et al. (2020). Essential Role of Accessory Subunit LYRM6 in the Mechanism of Mitochondrial Complex I. Nat. Commun. 11, 6008. doi:10.1038/s41467-020-19778-7
Gardeitchik, T., Mohamed, M., Ruzzenente, B., Karall, D., Guerrero-Castillo, S., Dalloyaux, D., et al. (2018). Bi-allelic Mutations in the Mitochondrial Ribosomal Protein MRPS2 Cause Sensorineural Hearing Loss, Hypoglycemia, and Multiple OXPHOS Complex Deficiencies. Am. J. Hum. Genet. 102, 685–695. doi:10.1016/j.ajhg.2018.02.012
Geiger, T., Wisniewski, J. R., Cox, J., Zanivan, S., Kruger, M., Ishihama, Y., et al. (2011). Use of Stable Isotope Labeling by Amino Acids in Cell Culture as a Spike-In Standard in Quantitative Proteomics. Nat. Protoc. 6, 147–157. doi:10.1038/nprot.2010.192
Gerovac, M., Vogel, J., and Smirnov, A. (2021). The World of Stable Ribonucleoproteins and its Mapping with Grad-Seq and Related Approaches. Front. Mol. Biosci. 8, 661448. doi:10.3389/fmolb.2021.661448
Giese, H., Ackermann, J., Heide, H., Bleier, L., Drose, S., Wittig, I., et al. (2015). NOVA: a Software to Analyze Complexome Profiling Data. Bioinformatics 31, 440–441. doi:10.1093/bioinformatics/btu623
Giese, H., Meisterknecht, J., Heidler, J., and Wittig, I. (2021). Mitochondrial Complexome Profiling. Methods Mol. Biol. 2192, 269–285. doi:10.1007/978-1-0716-0834-0_19
Gillet, L. C., Navarro, P., Tate, S., Röst, H., Selevsek, N., Reiter, L., et al. (2012). Targeted Data Extraction of the MS/MS Spectra Generated by Data-independent Acquisition: a New Concept for Consistent and Accurate Proteome Analysis. Mol. Cell Proteomics 11, O111. doi:10.1074/mcp.O111.016717
Giurgiu, M., Reinhard, J., Brauner, B., Dunger-Kaltenbach, I., Fobo, G., Frishman, G., et al. (2019). CORUM: the Comprehensive Resource of Mammalian Protein Complexes-2019. Nucleic Acids Res. 47, D559–D563. doi:10.1093/nar/gky973
Gomes, A. F., and Gozzo, F. C. (2010). Chemical Cross-Linking with a Diazirine Photoactivatable Cross-Linker Investigated by MALDI- and ESI-MS/MS. J. Mass. Spectrom. 45, 892–899. doi:10.1002/jms.1776
Gopalakrishna, S., Pearce, S. F., Dinan, A. M., Schober, F. A., Cipullo, M., Spåhr, H., et al. (2019). C6orf203 Is an RNA-Binding Protein Involved in Mitochondrial Protein Synthesis. Nucleic Acids Res. 47, 9386–9399. doi:10.1093/NAR/GKZ684
Guarás, A., Perales-Clemente, E., Calvo, E., Acín-Pérez, R., Loureiro-Lopez, M., Pujol, C., et al. (2016). The CoQH2/CoQ Ratio Serves as a Sensor of Respiratory Chain Efficiency. Cel Rep. 15, 197–209. doi:10.1016/j.celrep.2016.03.009
Guerrero‐Castillo, S., Van Strien, J., Brandt, U., and Arnold, S. (2021b). Ablation of Mitochondrial DNA Results in Widespread Remodeling of the Mitochondrial Complexome. EMBO J. 40, e108648. doi:10.15252/embj.2021108648
Guerrero-Castillo, S., Baertling, F., Kownatzki, D., Wessels, H. J., Arnold, S., Brandt, U., et al. (2017a). The Assembly Pathway of Mitochondrial Respiratory Chain Complex I. Cel Metab. 25, 128–139. doi:10.1016/j.cmet.2016.09.002
Guerrero-Castillo, S., Cabrera-Orefice, A., Huynen, M. A., and Arnold, S. (2017b). Identification and Evolutionary Analysis of Tissue-specific Isoforms of Mitochondrial Complex I Subunit NDUFV3. Biochim. Biophys. Acta (Bba) - Bioenerg. 1858, 208–217. doi:10.1016/j.bbabio.2016.12.004
Guerrero-Castillo, S., Krisp, C., Küchler, K., Arnold, S., Schlüter, H., and Gersting, S. W. (2021a). Multiplexed Complexome Profiling Using Tandem Mass Tags. Biochim. Biophys. Acta (Bba) - Bioenerg. 1862, 148448. doi:10.1016/j.bbabio.2021.148448
Haldar, K., Bhattacharjee, S., and Safeukui, I. (2018). Drug Resistance in Plasmodium. Nat. Rev. Microbiol. 16, 156–170. doi:10.1038/nrmicro.2017.161
Hällberg, B. M., and Larsson, N.-G. (2014). Making Proteins in the Powerhouse. Cel Metab. 20, 226–240. doi:10.1016/j.cmet.2014.07.001
Harb, O. S., Kissinger, J. C., and Roos, D. S. (2020). “ToxoDB: the Functional Genomic Resource for Toxoplasma and Related Organisms,” in Toxoplasma Gondii (Elsevier), 1021–1041. doi:10.1016/b978-0-12-815041-2.00023-2
Harner, M., Körner, C., Walther, D., Mokranjac, D., Kaesmacher, J., Welsch, U., et al. (2011). The Mitochondrial Contact Site Complex, a Determinant of Mitochondrial Architecture. EMBO J. 30, 4356–4370. doi:10.1038/emboj.2011.379
Haute, L. V., Hendrick, A. G., D’Souza, A. R., Powell, C. A., Rebelo-Guiomar, P., Harbour, M. E., et al. (2019). METTL15 Introduces N4-Methylcytidine into Human Mitochondrial 12S rRNA and Is Required for Mitoribosome Biogenesis. Nucleic Acids Res. 47, 10267–10281. doi:10.1093/nar/gkz735
Havugimana, P. C., Hart, G. T., Nepusz, T., Yang, H., Turinsky, A. L., Li, Z., et al. (2012). A Census of Human Soluble Protein Complexes. Cell 150, 1068–1081. doi:10.1016/j.cell.2012.08.011
Hebert, A. S., Merrill, A. E., Bailey, D. J., Still, A. J., Westphall, M. S., Strieter, E. R., et al. (2013). Neutron-encoded Mass Signatures for Multiplexed Proteome Quantification. Nat. Methods 10, 332–334. doi:10.1038/nmeth.2378
Heide, H., Bleier, L., Steger, M., Ackermann, J., Dröse, S., Schwamb, B., et al. (2012). Complexome Profiling Identifies TMEM126B as a Component of the Mitochondrial Complex I Assembly Complex. Cel Metab. 16, 538–549. doi:10.1016/j.cmet.2012.08.009
Hellman, L. M., and Fried, M. G. (2007). Electrophoretic Mobility Shift Assay (EMSA) for Detecting Protein-Nucleic Acid Interactions. Nat. Protoc. 2, 1849–1861. doi:10.1038/nprot.2007.249
Herrmannová, A., Prilepskaja, T., Wagner, S., Šikrová, D., Zeman, J., Poncová, K., et al. (2020). Adapted Formaldehyde Gradient Cross-Linking Protocol Implicates Human eIF3d and eIF3c, K and L Subunits in the 43S and 48S Pre-initiation Complex Assembly, Respectively. Nucleic Acids Res. 48, 1969–1984. doi:10.1093/nar/gkz1185
Herzog, F., Kahraman, A., Boehringer, D., Mak, R., Bracher, A., Walzthoeni, T., et al. (2012). Structural Probing of a Protein Phosphatase 2A Network by Chemical Cross-Linking and Mass Spectrometry. Science 337, 1348–1352. doi:10.1126/science.1221483
Heusel, M., Bludau, I., Rosenberger, G., Hafen, R., Frank, M., Banaei‐Esfahani, A., et al. (2019). Complex‐centric Proteome Profiling by SEC ‐ SWATH ‐ MS. Mol. Syst. Biol. 15, e8438. doi:10.15252/msb.20188438
Hevler, J. F., Lukassen, M. V., Cabrera‐Orefice, A., Arnold, S., Pronker, M. F., Franc, V., et al. (2021a). Selective Cross‐linking of Coinciding Protein Assemblies by In‐gel Cross‐linking Mass Spectrometry. EMBO J. 40, e106174. doi:10.15252/embj.2020106174
Hevler, J. F., Zenezeni Chiozzi, R., Cabrera-Orefice, A., Brandt, U., Arnold, S., and Heck, A. J. R. (2021b). Molecular Characterization of a Complex of Apoptosis-Inducing Factor 1 with Cytochrome C Oxidase of the Mitochondrial Respiratory Chain. Proc. Natl. Acad. Sci. USA 118, e2106950118. doi:10.1073/pnas.2106950118
Hilander, T., Jackson, C. B., Robciuc, M., Bashir, T., and Zhao, H. (2021). The Roles of Assembly Factors in Mammalian Mitoribosome Biogenesis. Mitochondrion 60, 70–84. doi:10.1016/J.MITO.2021.07.008
Hillier, C., Pardo, M., Yu, L., Bushell, E., Sanderson, T., Metcalf, T., et al. (2019). Landscape of the Plasmodium Interactome Reveals Both Conserved and Species-specific Functionality. Cel Rep. 28, 1635–1647. e1635. doi:10.1016/j.celrep.2019.07.019
Hirst, J. (2013). Mitochondrial Complex I. Annu. Rev. Biochem. 82, 551–575. doi:10.1146/annurev-biochem-070511-103700
Hocine, S., Singer, R. H., and Grünwald, D. (2010). RNA Processing and Export. Cold Spring Harbor Perspect. Biol. 2, a000752. doi:10.1101/CSHPERSPECT.A000752
Hong, P., Koza, S., and Bouvier, E. S. P. (2012). A Review Size-Exclusion Chromatography for the Analysis of Protein Biotherapeutics and Their Aggregates. J. Liquid Chromatogr. Relat. Tech. 35, 2923–2950. doi:10.1080/10826076.2012.743724
Honzík, T., Tesarová, M., Mayr, J. A., Hansíková, H., Jesina, P., Bodamer, O., et al. (2010). Mitochondrial Encephalocardio-Myopathy with Early Neonatal Onset Due to TMEM70 Mutation. Arch. Dis. Child. 95, 296–301. doi:10.1136/adc.2009.168096
Hör, J., Di Giorgio, S., Gerovac, M., Venturini, E., Förstner, K. U., and Vogel, J. (2020). Grad-seq Shines Light on Unrecognized RNA and Protein Complexes in the Model Bacterium Escherichia coli. Nucleic Acids Res. 48, 9301–9319. doi:10.1093/nar/gkaa676
Hu, A., Noble, W. S., and Wolf-Yadlin, A. (2016). Technical Advances in Proteomics: New Developments in Data-independent Acquisition. F1000Res 5, 419. doi:10.12688/f1000research.7042.1
Hu, L. Z., Goebels, F., Tan, J. H., Wolf, E., Kuzmanov, U., Wan, C., et al. (2019). EPIC: Software Toolkit for Elution Profile-Based Inference of Protein Complexes. Nat. Methods 16, 737–742. doi:10.1038/s41592-019-0461-4
Huynen, M. A., Mühlmeister, M., Gotthardt, K., Guerrero-Castillo, S., and Brandt, U. (2016). Evolution and Structural Organization of the Mitochondrial Contact Site (MICOS) Complex and the Mitochondrial Intermembrane Space Bridging (MIB) Complex. Biochim. Biophys. Acta (Bba) - Mol. Cel Res. 1863, 91–101. doi:10.1016/j.bbamcr.2015.10.009
Iacobucci, I., Monaco, V., Cozzolino, F., and Monti, M. (2021). From Classical to New Generation Approaches: An Excursus of -omics Methods for Investigation of Protein-Protein Interaction Networks. J. Proteomics 230, 103990. doi:10.1016/j.jprot.2020.103990
Jenkinson, E. M., Rehman, A. U., Walsh, T., Clayton-Smith, J., Lee, K., Morell, R. J., et al. (2013). Perrault Syndrome Is Caused by Recessive Mutations in CLPP, Encoding a Mitochondrial ATP-dependent Chambered Protease. Am. J. Hum. Genet. 92, 605–613. doi:10.1016/j.ajhg.2013.02.013
Jonckheere, A. I., Huigsloot, M., Lammens, M., Jansen, J., Van Den Heuvel, L. P., Spiekerkoetter, U., et al. (2011). Restoration of Complex V Deficiency Caused by a Novel Deletion in the Human TMEM70 Gene Normalizes Mitochondrial Morphology. Mitochondrion 11, 954–963. doi:10.1016/j.mito.2011.08.012
Kaila, V. R. I., and Wikström, M. (2021). Architecture of Bacterial Respiratory Chains. Nat. Rev. Microbiol. 19, 319–330. doi:10.1038/s41579-020-00486-4
Kartal, B., De Almeida, N. M., Maalcke, W. J., Op Den Camp, H. J. M., Jetten, M. S. M., and Keltjens, J. T. (2013). How to Make a Living from Anaerobic Ammonium Oxidation. FEMS Microbiol. Rev. 37, 428–461. doi:10.1111/1574-6976.12014
Kastritis, P. L., O'reilly, F. J., Bock, T., Li, Y., Rogon, M. Z., Buczak, K., et al. (2017). Capturing Protein Communities by Structural Proteomics in a Thermophilic Eukaryote. Mol. Syst. Biol. 13, 936. doi:10.15252/msb.20167412
Kiirika, L. M., Behrens, C., Braun, H.-P., and Colditz, F. (2013). The Mitochondrial Complexome of Medicago Truncatula. Front. Plant Sci. 4, 84. doi:10.3389/fpls.2013.00084
King, M. P., and Attardi, G. (1989). Human Cells Lacking mtDNA: Repopulation with Exogenous Mitochondria by Complementation. Science 246, 500–503. doi:10.1126/science.2814477
Kmita, K., Wirth, C., Warnau, J., Guerrero-Castillo, S., Hunte, C., Hummer, G., et al. (2015). Accessory NUMM (NDUFS6) Subunit Harbors a Zn-Binding Site and Is Essential for Biogenesis of Mitochondrial Complex I. Proc. Natl. Acad. Sci. USA 112, 5685–5690. doi:10.1073/pnas.1424353112
König, T., Tröder, S. E., Bakka, K., Korwitz, A., Richter-Dennerlein, R., Lampe, P. A., et al. (2016). The M -AAA Protease Associated with Neurodegeneration Limits MCU Activity in Mitochondria. Mol. Cel 64, 148–162. doi:10.1016/j.molcel.2016.08.020
Korepanova, A., and Matayoshi, E. D. (2012). HPLC‐SEC Characterization of Membrane Protein‐Detergent Complexes. Curr. Protoc. Protein Sci. 68–12. doi:10.1002/0471140864.ps2905s68
Krasny, L., and Huang, P. H. (2021). Data-independent Acquisition Mass Spectrometry (DIA-MS) for Proteomic Applications in Oncology. Mol. Omics 17, 29–42. doi:10.1039/d0mo00072h
Kristensen, A. R., Gsponer, J., and Foster, L. J. (2012). A High-Throughput Approach for Measuring Temporal Changes in the Interactome. Nat. Methods 9, 907–909. doi:10.1038/nmeth.2131
Kummer, E., and Ban, N. (2021). Mechanisms and Regulation of Protein Synthesis in Mitochondria. Nat. Rev. Mol. Cel Biol 22, 307–325. doi:10.1038/s41580-021-00332-2
Kyrilis, F. L., Belapure, J., and Kastritis, P. L. (2021). Detecting Protein Communities in Native Cell Extracts by Machine Learning: A Structural Biologist's Perspective. Front. Mol. Biosci. 8, 660542. doi:10.3389/fmolb.2021.660542
Ladig, R., Sommer, M. S., Hahn, A., Leisegang, M. S., Papasotiriou, D. G., Ibrahim, M., et al. (2011). A High-Definition Native Polyacrylamide Gel Electrophoresis System for the Analysis of Membrane Complexes. Plant J. 67, 181–194. doi:10.1111/j.1365-313X.2011.04577.x
Lapuente-Brun, E., Moreno-Loshuertos, R., Acín-Pérez, R., Latorre-Pellicer, A., Colás, C., Balsa, E., et al. (2013). Supercomplex Assembly Determines Electron Flux in the Mitochondrial Electron Transport Chain. Science 340, 1567–1570. doi:10.1126/science.1230381
Lasserre, J.-P., Beyne, E., Pyndiah, S., Lapaillerie, D., Claverol, S., and Bonneu, M. (2006). A Complexomic Study ofEscherichia Coli Using Two-Dimensional Blue Native/SDS Polyacrylamide Gel Electrophoresis. Electrophoresis 27, 3306–3321. doi:10.1002/elps.200500912
Lee, J.-H., and Skalnik, D. (2013). Simple and Efficient Identification of Chromatin Modifying Complexes and Characterization of Complex Composition. Methods Mol. Biol. 977, 289–298. doi:10.1007/978-1-62703-284-1_23
Liu, F., Lössl, P., Rabbitts, B. M., Balaban, R. S., and Heck, A. J. R. (2018). The Interactome of Intact Mitochondria by Cross-Linking Mass Spectrometry Provides Evidence for Coexisting Respiratory Supercomplexes. Mol. Cell Proteomics 17, 216–232. doi:10.1074/mcp.RA117.000470
Lobo‐Jarne, T., Pérez‐Pérez, R., Fontanesi, F., Timón‐Gómez, A., Wittig, I., Peñas, A., et al. (2020). Multiple Pathways Coordinate Assembly of Human Mitochondrial Complex IV and Stabilization of Respiratory Supercomplexes. EMBO J. 39, e103912. doi:10.15252/embj.2019103912
López De Heredia, M., and Jansen, R. P. (2004). RNA Integrity as a Quality Indicator during the First Steps of RNP Purifications: A Comparison of Yeast Lysis Methods. BMC Biochem. 5, 14. doi:10.1186/1471-2091-5-14
Low, T. Y., Syafruddin, S. E., Mohtar, M. A., Vellaichamy, A., A Rahman, N. S., Pung, Y.-F., et al. (2021). Recent Progress in Mass Spectrometry-Based Strategies for Elucidating Protein-Protein Interactions. Cell. Mol. Life Sci. 78, 5325–5339. doi:10.1007/s00018-021-03856-0
Maclean, A. E., Bridges, H. R., Silva, M. F., Ding, S., Ovciarikova, J., Hirst, J., et al. (2021). Complexome Profile of Toxoplasma Gondii Mitochondria Identifies Divergent Subunits of Respiratory Chain Complexes Including New Subunits of Cytochrome Bc1 Complex. Plos Pathog. 17, e1009301. doi:10.1371/journal.ppat.1009301
Macrae, J. I., Dixon, M. W., Dearnley, M. K., Chua, H. H., Chambers, J. M., Kenny, S., et al. (2013). Mitochondrial Metabolism of Sexual and Asexual Blood Stages of the Malaria Parasite Plasmodium Falciparum. BMC Biol. 11, 67. doi:10.1186/1741-7007-11-67
Mallam, A. L., Sae-Lee, W., Schaub, J. M., Tu, F., Battenhouse, A., Jang, Y. J., et al. (2019). Systematic Discovery of Endogenous Human Ribonucleoprotein Complexes. Cel Rep. 29, 1351–1368. e1355. doi:10.1016/j.celrep.2019.09.060
Mcalister, G. C., Nusinow, D. P., Jedrychowski, M. P., Wühr, M., Huttlin, E. L., Erickson, B. K., et al. (2014). MultiNotch MS3 Enables Accurate, Sensitive, and Multiplexed Detection of Differential Expression across Cancer Cell Line Proteomes. Anal. Chem. 86, 7150–7158. doi:10.1021/ac502040v
Meldal, B. H. M., Perfetto, L., Combe, C., Lubiana, T., Ferreira Cavalcante, J. V., Bye-A-Jee, H., et al. (2021). Complex Portal 2022: New Curation Frontiers. Nucleic Acids Res., gkab991. doi:10.1093/nar/gkab991
Merrill, A. E., Hebert, A. S., Macgilvray, M. E., Rose, C. M., Bailey, D. J., Bradley, J. C., et al. (2014). NeuCode Labels for Relative Protein Quantification. Mol. Cell Proteomics 13, 2503–2512. doi:10.1074/mcp.M114.040287
Mili, S., and Steitz, J. A. (2004). Evidence for Reassociation of RNA-Binding Proteins after Cell Lysis: Implications for the Interpretation of Immunoprecipitation Analyses. RNA 10, 1692–1694. doi:10.1261/rna.7151404
Mintseris, J., and Gygi, S. P. (2020). High-density Chemical Cross-Linking for Modeling Protein Interactions. Proc. Natl. Acad. Sci. USA 117, 93–102. doi:10.1073/pnas.1902931116
Mistry, J., Chuguransky, S., Williams, L., Qureshi, M., Salazar, G. A., Sonnhammer, E. L. L., et al. (2021). Pfam: The Protein Families Database in 2021. Nucleic Acids Res. 49, D412–D419. doi:10.1093/nar/gkaa913
Mohamed, M., Gardeitchik, T., Balasubramaniam, S., Guerrero‐Castillo, S., Dalloyaux, D., Kraaij, S., et al. (2020). Novel Defect in Phosphatidylinositol 4‐kinase Type 2‐alpha ( PI4K2A ) at the Membrane‐enzyme Interface Is Associated with Metabolic Cutis Laxa. Jrnl Inher Metab. Disea 43, 1382–1391. doi:10.1002/jimd.12255
Müller, C. S., Bildl, W., Haupt, A., Ellenrieder, L., Becker, T., Hunte, C., et al. (2016). Cryo-slicing Blue Native-Mass Spectrometry (csBN-MS), a Novel Technology for High Resolution Complexome Profiling. Mol. Cell Proteomics 15, 669–681. doi:10.1074/mcp.M115.054080
Munawar, N., Olivero, G., Jerman, E., Doyle, B., Streubel, G., Wynne, K., et al. (2015). Native Gel Analysis of Macromolecular Protein Complexes in Cultured Mammalian Cells. Proteomics 15, 3603–3612. doi:10.1002/pmic.201500045
Neustein, H. B., Lurie, P. R., Dahms, B., and Takahashi, M. (1979). An X-Linked Recessive Cardiomyopathy with Abnormal Mitochondria. Pediatrics 64, 24–29.
Nolte, H., and Langer, T. (2021). ComplexFinder: A Software Package for the Analysis of Native Protein Complex Fractionation Experiments. Biochim. Biophys. Acta (Bba) - Bioenerg. 1862, 148444. doi:10.1016/j.bbabio.2021.148444
Nuebel, E., Morgan, J. T., Fogarty, S., Winter, J. M., Lettlova, S., Berg, J. A., et al. (2021). The Biochemical Basis of Mitochondrial Dysfunction in Zellweger Spectrum Disorder. EMBO Rep. 22, e51991. doi:10.15252/embr.202051991
Ong, S.-E., Blagoev, B., Kratchmarova, I., Kristensen, D. B., Steen, H., Pandey, A., et al. (2002). Stable Isotope Labeling by Amino Acids in Cell Culture, SILAC, as a Simple and Accurate Approach to Expression Proteomics. Mol. Cell Proteomics 1, 376–386. doi:10.1074/mcp.m200025-mcp200
O’Reilly, F. J., and Rappsilber, J. (2018). Cross-linking Mass Spectrometry: Methods and Applications in Structural, Molecular and Systems Biology. Nat. Struct. Mol. Biol. 25, 1000–1008. doi:10.1038/s41594-018-0147-0
Orsburn, B. C. (2021). Proteome Discoverer-A Community Enhanced Data Processing Suite for Protein Informatics. Proteomes 9, 15. doi:10.3390/proteomes9010015
Ott, C., Dorsch, E., Fraunholz, M., Straub, S., and Kozjak-Pavlovic, V. (2015). Detailed Analysis of the Human Mitochondrial Contact Site Complex Indicate a Hierarchy of Subunits. PLoS One 10, e0120213. doi:10.1371/journal.pone.0120213
Oughtred, R., Rust, J., Chang, C., Breitkreutz, B. J., Stark, C., Willems, A., et al. (2021). TheBioGRIDdatabase: A Comprehensive Biomedical Resource of Curated Protein, Genetic, and Chemical Interactions. Protein Sci. 30, 187–200. doi:10.1002/pro.3978
Overmyer, K. A., Tyanova, S., Hebert, A. S., Westphall, M. S., Cox, J., and Coon, J. J. (2018). Multiplexed Proteome Analysis with Neutron-Encoded Stable Isotope Labeling in Cells and Mice. Nat. Protoc. 13, 293–306. doi:10.1038/nprot.2017.121
Páleníková, P., Harbour, M. E., Ding, S., Fearnley, I. M., Van Haute, L., Rorbach, J., et al. (2021a). Quantitative Density Gradient Analysis by Mass Spectrometry (qDGMS) and Complexome Profiling Analysis (ComPrAn) R Package for the Study of Macromolecular Complexes. Biochim. Biophys. Acta (Bba) - Bioenerg. 1862, 148399. doi:10.1016/j.bbabio.2021.148399
Páleníková, P., Harbour, M. E., Prodi, F., Minczuk, M., Zeviani, M., Ghelli, A., et al. (2021b). Duplexing Complexome Profiling with SILAC to Study Human Respiratory Chain Assembly Defects. Biochim. Biophys. Acta (Bba) - Bioenerg. 1862, 148395. doi:10.1016/j.bbabio.2021.148395
Patton, R. D., Sanjeev, M., Woodward, L. A., Mabin, J. W., Bundschuh, R., and Singh, G. (2020). Chemical Crosslinking Enhances RNA Immunoprecipitation for Efficient Identification of Binding Sites of Proteins that Photo-Crosslink Poorly with RNA. RNA 26, 1216–1233. doi:10.1261/rna.074856.120
Pearce, S., Nezich, C. L., and Spinazzola, A. (2013). Mitochondrial Diseases: Translation Matters. Mol. Cell Neurosci. 55, 1–12. doi:10.1016/j.mcn.2012.08.013
Pereira, B., Videira, A., and Duarte, M. (2013). Novel Insights into the Role of Neurospora Crassa NDUFAF2, an Evolutionarily Conserved Mitochondrial Complex I Assembly Factor. Mol. Cel Biol 33, 2623–2634. doi:10.1128/MCB.01476-12
Pérez-Pérez, R., Lobo-Jarne, T., Milenkovic, D., Mourier, A., Bratic, A., García-Bartolomé, A., et al. (2016). COX7A2L Is a Mitochondrial Complex III Binding Protein that Stabilizes the III2+IV Supercomplex without Affecting Respirasome Formation. Cel Rep. 16, 2387–2398. doi:10.1016/j.celrep.2016.07.081
Perez-Riverol, Y., Csordas, A., Bai, J., Bernal-Llinares, M., Hewapathirana, S., Kundu, D. J., et al. (2019). The PRIDE Database and Related Tools and Resources in 2019: Improving Support for Quantification Data. Nucleic Acids Res. 47, D442–D450. doi:10.1093/nar/gky1106
Perkins, D. N., Pappin, D. J. C., Creasy, D. M., and Cottrell, J. S. (1999). Probability-based Protein Identification by Searching Sequence Databases Using Mass Spectrometry Data. Electrophoresis 20, 3551–3567. doi:10.1002/(sici)1522-2683(19991201)20:18<3551::aid-elps3551>3.0.co;2-2
Pitkänen, L., Tuomainen, P., and Eskelin, K. (2014). Analysis of Plant Ribosomes with Asymmetric Flow Field-Flow Fractionation. Anal. Bioanal. Chem. 406, 1629–1637. doi:10.1007/s00216-013-7454-4
Poole, R. K., and Cook, G. M. (2000). Redundancy of Aerobic Respiratory Chains in Bacteria? Routes, Reasons and Regulation. Adv. Microb. Physiol. 43, 165–224. doi:10.1016/s0065-2911(00)43005-5
Protasoni, M., Pérez‐Pérez, R., Lobo‐Jarne, T., Harbour, M. E., Ding, S., Peñas, A., et al. (2020). Respiratory Supercomplexes Act as a Platform for Complex III ‐mediated Maturation of Human Mitochondrial Complexes I and IV. EMBO J. 39, e102817. doi:10.15252/embj.2019102817
Reschiglian, P., and Moon, M. H. (2008). Flow Field-Flow Fractionation: a Pre-analytical Method for Proteomics. J. Proteomics 71, 265–276. doi:10.1016/j.jprot.2008.06.002
Rhein, V. F., Carroll, J., Ding, S., Fearnley, I. M., and Walker, J. E. (2016). NDUFAF5 Hydroxylates NDUFS7 at an Early Stage in the Assembly of Human Complex I. J. Biol. Chem. 291, 14851–14860. doi:10.1074/jbc.M116.734970
Robinson, P. J., Trnka, M. J., Pellarin, R., Greenberg, C. H., Bushnell, D. A., Davis, R., et al. (2015). Molecular Architecture of the Yeast Mediator Complex. Elife 4, e08719. doi:10.7554/eLife.08719
Rodenburg, R. J. (2016). Mitochondrial Complex I-Linked Disease. Biochim. Biophys. Acta (Bba) - Bioenerg. 1857, 938–945. doi:10.1016/j.bbabio.2016.02.012
Röst, H. L., Sachsenberg, T., Aiche, S., Bielow, C., Weisser, H., Aicheler, F., et al. (2016). OpenMS: a Flexible Open-Source Software Platform for Mass Spectrometry Data Analysis. Nat. Methods 13, 741–748. doi:10.1038/nmeth.3959
Rugen, N., Straube, H., Franken, L. E., Braun, H.-P., and Eubel, H. (2019). Complexome Profiling Reveals Association of PPR Proteins with Ribosomes in the Mitochondria of Plants. Mol. Cell Proteomics 18, 1345–1362. doi:10.1074/mcp.RA119.001396
Ryl, P. S. J., Bohlke-Schneider, M., Lenz, S., Fischer, L., Budzinski, L., Stuiver, M., et al. (2020). In Situ Structural Restraints from Cross-Linking Mass Spectrometry in Human Mitochondria. J. Proteome Res. 19, 327–336. doi:10.1021/acs.jproteome.9b00541
Sahr, T., and Buchrieser, C. (2013). Co-immunoprecipitation: Protein-RNA and Protein-DNA Interaction. Methods Mol. Biol. 954, 583–593. doi:10.1007/978-1-62703-161-5_36
Sánchez-Caballero, L., Elurbe, D. M., Baertling, F., Guerrero-Castillo, S., Van Den Brand, M., Van Strien, J., et al. (2020). TMEM70 Functions in the Assembly of Complexes I and V. Biochim. Biophys. Acta (Bba) - Bioenerg. 1861, 148202. doi:10.1016/j.bbabio.2020.148202
Sánchez-Caballero, L., Guerrero-Castillo, S., and Nijtmans, L. (2016a). Unraveling the Complexity of Mitochondrial Complex I Assembly: A Dynamic Process. Biochim. Biophys. Acta (Bba) - Bioenerg. 1857, 980–990. doi:10.1016/j.bbabio.2016.03.031
Sánchez-Caballero, L., Ruzzenente, B., Bianchi, L., Assouline, Z., Barcia, G., Metodiev, M. D., et al. (2016b). Mutations in Complex I Assembly Factor TMEM126B Result in Muscle Weakness and Isolated Complex I Deficiency. Am. J. Hum. Genet. 99, 208–216. doi:10.1016/j.ajhg.2016.05.022
Schwanhäusser, B., Busse, D., Li, N., Dittmar, G., Schuchhardt, J., Wolf, J., et al. (2011). Global Quantification of Mammalian Gene Expression Control. Nature 473, 337–342. doi:10.1038/nature10098
Schwanhäusser, B., Gossen, M., Dittmar, G., and Selbach, M. (2009). Global Analysis of Cellular Protein Translation by Pulsed SILAC. Proteomics 9, 205–209. doi:10.1002/pmic.200800275
Schweppe, D. K., Chavez, J. D., Lee, C. F., Caudal, A., Kruse, S. E., Stuppard, R., et al. (2017). Mitochondrial Protein Interactome Elucidated by Chemical Cross-Linking Mass Spectrometry. Proc. Natl. Acad. Sci. USA 114, 1732–1737. doi:10.1073/pnas.1617220114
Senkler, J., Rugen, N., Eubel, H., Hegermann, J., and Braun, H.-P. (2018). Absence of Complex I Implicates Rearrangement of the Respiratory Chain in European Mistletoe. Curr. Biol. 28, 1606–1613. e1604. doi:10.1016/j.cub.2018.03.050
Senkler, J., Senkler, M., Eubel, H., Hildebrandt, T., Lengwenus, C., Schertl, P., et al. (2017). The Mitochondrial Complexome of Arabidopsis thaliana. Plant J. 89, 1079–1092. doi:10.1111/tpj.13448
Singhal, R. K., Kruse, C., Heidler, J., Strecker, V., Zwicker, K., Düsterwald, L., et al. (2017). Coi1 Is a Novel Assembly Factor of the Yeast Complex III-Complex IV Supercomplex. MBoC 28, 2609–2622. doi:10.1091/mbc.E17-02-0093
Skinnider, M. A., Scott, N. E., Prudova, A., Kerr, C. H., Stoynov, N., Stacey, R. G., et al. (2021). An Atlas of Protein-Protein Interactions across Mouse Tissues. Cell 184, 4073–4089. e4017. doi:10.1016/j.cell.2021.06.003
Smirnov, A., Förstner, K. U., Holmqvist, E., Otto, A., Günster, R., Becher, D., et al. (2016). Grad-seq Guides the Discovery of ProQ as a Major Small RNA-Binding Protein. Proc. Natl. Acad. Sci. USA 113, 11591–11596. doi:10.1073/pnas.1609981113
Smits, P., Saada, A., Wortmann, S. B., Heister, A. J., Brink, M., Pfundt, R., et al. (2011). Mutation in mitochondrial ribosomal protein MRPS22 leads to Cornelia de Lange-like phenotype, brain abnormalities and hypertrophic cardiomyopathy. Eur. J. Hum. Genet. 19, 394–399. doi:10.1038/ejhg.2010.214
Sonnett, M., Yeung, E., and Wühr, M. (2018). Accurate, Sensitive, and Precise Multiplexed Proteomics Using the Complement Reporter Ion Cluster. Anal. Chem. 90, 5032–5039. doi:10.1021/acs.analchem.7b04713
Spaniol, B., Lang, J., Venn, B., Schake, L., Sommer, F., Mustas, M., et al. (2021). Complexome Profiling on the Chlamydomonas Lpa2 Mutant Reveals Insights into PSII Biogenesis and New PSII Associated Proteins. J. Exp. Bot., erab390. doi:10.1093/jxb/erab390
Sprink, T., Ramrath, D. J. F., Yamamoto, H., Yamamoto, K., Loerke, J., Ismer, J., et al. (2016). Structures of Ribosome-Bound Initiation Factor 2 Reveal the Mechanism of Subunit Association. Sci. Adv. 2, e1501502. doi:10.1126/SCIADV.1501502
Stacey, R. G., Skinnider, M. A., Scott, N. E., and Foster, L. J. (2017). A Rapid and Accurate Approach for Prediction of Interactomes from Co-elution Data (PrInCE). BMC Bioinformatics 18, 457. doi:10.1186/s12859-017-1865-8
Steigenberger, B., Albanese, P., Heck, A. J. R., and Scheltema, R. A. (2020). To Cleave or Not to Cleave in XL-MS? J. Am. Soc. Mass. Spectrom. 31, 196–206. doi:10.1021/jasms.9b00085
Stenton, S. L., Sheremet, N. L., Catarino, C. B., Andreeva, N. A., Assouline, Z., Barboni, P., et al. (2021). Impaired Complex I Repair Causes Recessive Leber's Hereditary Optic Neuropathy. J. Clin. Invest. 131. doi:10.1172/JCI138267
Strecker, V., Wumaier, Z., Wittig, I., and Schägger, H. (2010). Large Pore Gels to Separate Mega Protein Complexes Larger Than 10 MDa by Blue Native Electrophoresis: Isolation of Putative Respiratory Strings or Patches. Proteomics 10, 3379–3387. doi:10.1002/pmic.201000343
Szczepanowska, K., Senft, K., Heidler, J., Herholz, M., Kukat, A., Höhne, M. N., et al. (2020). A Salvage Pathway Maintains Highly Functional Respiratory Complex I. Nat. Commun. 11, 1643. doi:10.1038/s41467-020-15467-7
Szklarczyk, D., Gable, A. L., Nastou, K. C., Lyon, D., Kirsch, R., Pyysalo, S., et al. (2021). The STRING Database in 2021: Customizable Protein-Protein Networks, and Functional Characterization of User-Uploaded Gene/measurement Sets. Nucleic Acids Res. 49, D605–D612. doi:10.1093/nar/gkaa1074
Tran, N. H., Zhang, X., Xin, L., Shan, B., and Li, M. (2017). De Novo peptide Sequencing by Deep Learning. Proc. Natl. Acad. Sci. USA 114, 8247–8252. doi:10.1073/pnas.1705691114
Trendel, J., Schwarzl, T., Horos, R., Prakash, A., Bateman, A., Hentze, M. W., et al. (2019). The Human RNA-Binding Proteome and its Dynamics during Translational Arrest. Cell 176, 391–403. e319. doi:10.1016/j.cell.2018.11.004
Tyanova, S., Temu, T., and Cox, J. (2016a). The MaxQuant Computational Platform for Mass Spectrometry-Based Shotgun Proteomics. Nat. Protoc. 11, 2301–2319. doi:10.1038/nprot.2016.136
Tyanova, S., Temu, T., Sinitcyn, P., Carlson, A., Hein, M. Y., Geiger, T., et al. (2016b). The Perseus Computational Platform for Comprehensive Analysis of (Prote)omics Data. Nat. Methods 13, 731–740. doi:10.1038/nmeth.3901
Ugalde, C., Vogel, R., Huijbens, R., Van Den Heuvel, B., Smeitink, J., and Nijtmans, L. (2004). Human Mitochondrial Complex I Assembles through the Combination of Evolutionary Conserved Modules: a Framework to Interpret Complex I Deficiencies. Hum. Mol. Genet. 13, 2461–2472. doi:10.1093/hmg/ddh262
Uniprot, C. (2021). UniProt: the Universal Protein Knowledgebase in 2021. Nucleic Acids Res. 49, D480–D489. doi:10.1093/nar/gkaa1100
Urdaneta, E. C., and Beckmann, B. M. (2020). Fast and Unbiased Purification of RNA-Protein Complexes after UV Cross-Linking. Methods 178, 72–82. doi:10.1016/J.YMETH.2019.09.013
Urdaneta, E. C., Vieira-Vieira, C. H., Hick, T., Wessels, H.-H., Figini, D., Moschall, R., et al. (2019). Purification of Cross-Linked RNA-Protein Complexes by Phenol-Toluol Extraction. Nat. Commun. 10, 990. doi:10.1038/s41467-019-08942-3
Van Damme, T., Gardeitchik, T., Mohamed, M., Guerrero-Castillo, S., Freisinger, P., Guillemyn, B., et al. (2017). Mutations in ATP6V1E1 or ATP6V1A Cause Autosomal-Recessive Cutis Laxa. Am. J. Hum. Genet. 100, 216–227. doi:10.1016/j.ajhg.2016.12.010
Van Esveld, S. L., and Spelbrink, J. N. (2021). RNA Crosslinking to Analyze the Mitochondrial RNA-Binding Proteome. Methods Mol. Biol. 2192, 147–158. doi:10.1007/978-1-0716-0834-0_12
Van Strien, J., Guerrero-Castillo, S., Chatzispyrou, I. A., Houtkooper, R. H., Brandt, U., and Huynen, M. A. (2019). COmplexome Profiling ALignment (COPAL) Reveals Remodeling of Mitochondrial Protein Complexes in Barth Syndrome. Bioinformatics 35, 3083–3091. doi:10.1093/bioinformatics/btz025
Van Strien, J., Haupt, A., Schulte, U., Braun, H.-P., Cabrera-Orefice, A., Choudhary, J. S., et al. (2021). CEDAR, an Online Resource for the Reporting and Exploration of Complexome Profiling Data. Biochim. Biophys. Acta (Bba) - Bioenerg. 1862, 148411. doi:10.1016/j.bbabio.2021.148411
Vercellino, I., and Sazanov, L. A. (2021). The Assembly, Regulation and Function of the Mitochondrial Respiratory Chain. Nat. Rev. Mol. Cel Biol. doi:10.1038/s41580-021-00415-0
Versantvoort, W., Guerrero-Castillo, S., Wessels, H. J. C. T., Van Niftrik, L., Jetten, M. S. M., Brandt, U., et al. (2019). Complexome Analysis of the Nitrite-dependent Methanotroph Methylomirabilis Lanthanidiphila. Biochim. Biophys. Acta (Bba) - Bioenerg. 1860, 734–744. doi:10.1016/j.bbabio.2019.07.011
Vidal, M., Cusick, M. E., and Barabási, A.-L. (2011). Interactome Networks and Human Disease. Cell 144, 986–998. doi:10.1016/j.cell.2011.02.016
Vidoni, S., Harbour, M. E., Guerrero-Castillo, S., Signes, A., Ding, S., Fearnley, I. M., et al. (2017). MR-1S Interacts with PET100 and PET117 in Module-Based Assembly of Human Cytochrome C Oxidase. Cel Rep. 18, 1727–1738. doi:10.1016/j.celrep.2017.01.044
Wai, T., Saita, S., Nolte, H., Müller, S., König, T., Richter‐Dennerlein, R., et al. (2016). The Membrane Scaffold SLP2 Anchors a Proteolytic Hub in Mitochondria Containing PARL and the I ‐AAA Protease YME1L. EMBO Rep. 17, 1844–1856. doi:10.15252/embr.201642698
Wang, X., Teferedegne, B., Shatzkes, K., Tu, W., and Murata, H. (2016). Endogenous RNase Inhibitor Contributes to Stability of RNA in Crude Cell Lysates: Applicability to RT-qPCR. Anal. Biochem. 513, 21–27. doi:10.1016/J.AB.2016.08.011
Weber, T. A., Koob, S., Heide, H., Wittig, I., Head, B., Van Der Bliek, A., et al. (2013). APOOL Is a Cardiolipin-Binding Constituent of the Mitofilin/MINOS Protein Complex Determining Cristae Morphology in Mammalian Mitochondria. PLoS One 8, e63683. doi:10.1371/journal.pone.0063683
Wessels, H. J. C. T., Vogel, R. O., Lightowlers, R. N., Spelbrink, J. N., Rodenburg, R. J., Van Den Heuvel, L. P., et al. (2013). Analysis of 953 Human Proteins from a Mitochondrial HEK293 Fraction by Complexome Profiling. PLoS One 8, e68340. doi:10.1371/journal.pone.0068340
Wessels, H. J. C. T., Vogel, R. O., Van Den Heuvel, L., Smeitink, J. A., Rodenburg, R. J., Nijtmans, L. G., et al. (2009). LC-MS/MS as an Alternative for SDS-PAGE in Blue Native Analysis of Protein Complexes. Proteomics 9, 4221–4228. doi:10.1002/pmic.200900157
Wittig, I., Beckhaus, T., Wumaier, Z., Karas, M., and Schägger, H. (2010). Mass Estimation of Native Proteins by Blue Native Electrophoresis. Mol. Cell Proteomics 9, 2149–2161. doi:10.1074/mcp.M900526-MCP200
Wittig, I., Braun, H.-P., and Schägger, H. (2006). Blue Native PAGE. Nat. Protoc. 1, 418–428. doi:10.1038/nprot.2006.62
Wittig, I., Karas, M., and Schägger, H. (2007). High Resolution clear Native Electrophoresis for In-Gel Functional Assays and Fluorescence Studies of Membrane Protein Complexes. Mol. Cell Proteomics 6, 1215–1225. doi:10.1074/mcp.M700076-MCP200
Wittig, I., and Malacarne, P. F. (2021). Complexome Profiling: Assembly and Remodeling of Protein Complexes. Ijms 22, 7809. doi:10.3390/ijms22157809
Wöhlbrand, L., Ruppersberg, H. S., Feenders, C., Blasius, B., Braun, H.-P., and Rabus, R. (2016). Analysis of Membrane-Protein Complexes of the marine Sulfate Reducer Desulfobacula Toluolica Tol2 by 1D Blue Native-PAGE Complexome Profiling and 2D Blue Native-/SDS-PAGE. PROTEOMICS 16, 973–988. doi:10.1002/pmic.201500360
Wortmann, S. B., Ziętkiewicz, S., Guerrero-Castillo, S., Feichtinger, R. G., Wagner, M., Russell, J., et al. (2021). Neutropenia and Intellectual Disability Are Hallmarks of Biallelic and De Novo CLPB Deficiency. Genet. Med. 23, 1705–1714. doi:10.1038/s41436-021-01194-x
Xu, Y., Malhotra, A., Ren, M., and Schlame, M. (2006). The Enzymatic Function of Tafazzin. J. Biol. Chem. 281, 39217–39224. doi:10.1074/jbc.M606100200
Yoshikawa, H., Larance, M., Harney, D. J., Sundaramoorthy, R., Ly, T., Owen-Hughes, T., et al. (2018). Efficient Analysis of Mammalian Polysomes in Cells and Tissues Using Ribo Mega-SEC. Elife 7. doi:10.7554/eLife.36530
Yu, H., Lu, J. J., Rao, W., and Liu, S. (20162016). Capitalizing Resolving Power of Density Gradient Ultracentrifugation by Freezing and Precisely Slicing Centrifuged Solution: Enabling Identification of Complex Proteins from Mitochondria by Matrix Assisted Laser Desorption/Ionization Time-Of-Flight Mass Spectrometry. J. Anal. Methods Chem. 2016, 1–7. doi:10.1155/2016/8183656
Yu, Y., Ji, H., Doudna, J. A., and Leary, J. A. (2005). Mass Spectrometric Analysis of the Human 40S Ribosomal Subunit: Native and HCV IRES-Bound Complexes. Protein Sci. 14, 1438–1446. doi:10.1110/ps.041293005
Zhang, J., Xin, L., Shan, B., Chen, W., Xie, M., Yuen, D., et al. (2012). PEAKS DB: De Novo Sequencing Assisted Database Search for Sensitive and Accurate Peptide Identification. Mol. Cell Proteomics 11, M111. doi:10.1074/mcp.M111.010587
Keywords: complexome profiling, mitochondria, protein complex, protein-protein interaction (PPI), oxidative phosphorylation, disease, proteomics, mass spectrometry
Citation: Cabrera-Orefice A, Potter A, Evers F, Hevler JF and Guerrero-Castillo S (2022) Complexome Profiling—Exploring Mitochondrial Protein Complexes in Health and Disease. Front. Cell Dev. Biol. 9:796128. doi: 10.3389/fcell.2021.796128
Received: 15 October 2021; Accepted: 08 December 2021;
Published: 12 January 2022.
Edited by:
Erika Fernandez-Vizarra, Veneto Institute of Molecular Medicine (VIMM), ItalyReviewed by:
Ilka Wittig, Goethe University Frankfurt, GermanyMassimo Zeviani, University of Padua, Italy
Copyright © 2022 Cabrera-Orefice, Potter, Evers, Hevler and Guerrero-Castillo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alfredo Cabrera-Orefice, alfredbiomed@gmail.com; Sergio Guerrero-Castillo, s.guerrerocastillo@uke.de
†These authors have contributed equally to this work