Visualization and analysis of the complexome network of Saccharomyces cerevisiae

Most processes in the cell are delivered by protein complexes, rather than individual proteins. While the association of proteins has been studied extensively in protein-protein interaction networks (the interactome), an intuitive and effective representation of complex-complex connections (the complexome) is not yet available. Here, we describe a new representation of the complexome of Saccharomyces cerevisiae. Using the core-module-attachment data of Gavin et al. ( Nature 2006 , 440 , 631 - 6 ), protein complexes in the network are represented as nodes; these are connected by edges that represent shared core and/or module protein subunits. To validate this network, we examined the network topology and its distribution of biological processes. The complexome network showed scale-free characteristics, with a power law-like node degree distribution and clustering coefficient independent of node degree. Connected complexes in the network showed similarities in biological process that were nonrandom. Furthermore, clusters of interacting complexes reflected a higher-level organization of many cellular functions. The strong functional relationships seen in these clusters, along with literature evidence, allowed 44 uncharacterized complexes to be assigned putative functions using guilt-by-association. We demonstrate our network model using the GEOMI visualization platform, on which we have developed capabilities to integrate and visualize complexome data.     

Analysis of 953 Human Proteins from a Mitochondrial HEK293 Fraction by Complexome Profiling

Complexome profiling is a novel technique which uses shotgun proteomics to establish protein migration profiles from fractionated blue native electrophoresis gels. Here we present a dataset of blue native electrophoresis migration profiles for 953 proteins by complexome profiling. By analysis of mitochondrial ribosomal complexes we demonstrate its potential to verify putative protein-protein interactions identified by affinity purification – mass spectrometry studies. Protein complexes were extracted in their native state from a HEK293 mitochondrial fraction and separated by blue native gel electrophoresis. Gel lanes were cut into gel slices of even size and analyzed by shotgun proteomics. Subsequently, the acquired protein migration profiles were analyzed for co-migration via hierarchical cluster analysis. This dataset holds great promise as a comprehensive resource for de novo identification of protein-protein interactions or to underpin and prioritize candidate protein interactions from other studies. To demonstrate the potential use of our dataset we focussed on the mitochondrial translation machinery. Our results show that mitoribosomal complexes can be analyzed by blue native gel electrophoresis, as at least four distinct complexes. Analysis of these complexes confirmed that 24 proteins that had previously been reported to co-purify with mitoribosomes indeed co-migrated with subunits of the mitochondrial ribosome. Co-migration of several proteins involved in biogenesis of inner mitochondrial membrane complexes together with mitoribosomal complexes suggested the possibility of co-translational assembly in human cells. Our data also highlighted a putative ribonucleotide complex that potentially contains MRPL10, MRPL12 and MRPL53 together with LRPPRC and SLIRP.     

The mitochondrial complexome of Arabidopsis thaliana

Mitochondria are central to cellular metabolism and energy conversion. In plants they also enable photosynthesis through additional components and functional flexibility. A majority of those processes relies on the assembly of individual proteins to larger protein complexes, some of which operate as large molecular machines. There has been a strong interest in the makeup and function of mitochondrial protein complexes and protein–protein interactions in plants, but the experimental approaches used typically suffer from selectivity or bias. Here, we present a complexome profiling analysis for leaf mitochondria of the model plant Arabidopsis thaliana for the systematic characterization of protein assemblies. Purified organelle extracts were separated by 1D Blue native (BN) PAGE, a resulting gel lane was dissected into 70 slices (complexome fractions) and proteins in each slice were identified by label free quantitative shot‐gun proteomics. Overall, 1359 unique proteins were identified, which were, on average, present in 17 complexome fractions each. Quantitative profiles of proteins along the BN gel lane were aligned by similarity, allowing us to visualize protein assemblies. The data allow re‐annotating the subunit compositions of OXPHOS complexes, identifying assembly intermediates of OXPHOS complexes and assemblies of alternative respiratory oxidoreductases. Several protein complexes were discovered that have not yet been reported in plants, such as a 530 kDa Tat complex, 460 and 1000 kDa SAM complexes, a calcium ion uniporter complex (150 kDa) and several PPR protein complexes. We have set up a tailored online resource ( https://complexomemap.de/at_mito_leaves ) to deposit the data and to allow straightforward access and custom data analyses.     

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

Sulfate‐reducing bacteria (SRB) obtain energy from cytoplasmic reduction of sulfate to sulfide involving APS‐reductase (AprAB) and dissimilatory sulfite reductase (DsrAB). These enzymes are predicted to obtain electrons from membrane redox complexes, i.e. the quinone‐interacting membrane‐bound oxidoreductase (QmoABC) and DsrMKJOP complexes. In addition to these conserved complexes, the genomes of SRB encode a large number of other (predicted) membrane redox complexes, the function and actual formation of which is unknown. This study reports the establishment of 1D Blue Native‐PAGE complexome profiling and 2D BN‐/SDS‐PAGE for analysis of the membrane protein complexome of the marine sulfate reducer Desulfobacula toluolica Tol2. Analysis of normalized score profiles of >800 proteins in combination with hierarchical clustering and identification of 2D BN‐/SDS‐PAGE separated spots demonstrated separation of membrane complexes in their native form, e.g. ATP synthase. In addition to the QmoABC and DsrMKJOP complexes, other complexes were detected that constitute the basic membrane complexome of D. toluolica Tol2, e.g. transport proteins (e.g. sodium/sulfate symporters) or redox complexes involved in Na⁺‐based bioenergetics (RnfABCDEG). Notably, size estimation indicates dimer and quadruple formation of the DsrMKJOP complex in vivo. Furthermore, cluster analysis suggests interaction of this complex with a rhodanese‐like protein (Tol2_C05230) possibly representing a periplasmic electron transfer partner for DsrMKJOP.     

Characterization of multiprotein complexes of the Borrelia burgdorferi outer membrane vesicles

Among bacterial cell envelopes, the Borrelia burgdorferi outer membrane (OM) is structurally unique in that the identities of many protein complexes remain unknown; however, their characterization is the first step toward our understanding of membrane protein interactions and potential functions. Here, we used two-dimensional blue native/SDS-PAGE/mass spectrometric analysis for a global characterization of protein-protein interactions as well as to identify protein complexes in OM vesicles isolated from multiple infectious sensu stricto isolates of B. burgdorferi. Although we uncovered the existence of at least 10 distinct OM complexes harboring several unique subunits, the complexome is dominated by the frequent occurrence of a limited diversity of membrane proteins, most notably P13, outer surface protein (Osp) A, -B, -C, and -D and Lp6.6. The occurrence of these complexes and specificity of subunit interaction were further supported by independent two-dimensional immunoblotting and coimmunoprecipitation assays as well as by mutagenesis studies, where targeted depletion of a subunit member (P66) selectively abolished a specific complex. Although a comparable profile of the OM complexome was detected in two major infectious isolates, such as B31 and 297, certain complexes are likely to occur in an isolate-specific manner. Further assessment of protein complexes in multiple Osp-deficient isolates showed loss of several protein complexes but revealed the existence of additional complex/subunits that are undetectable in wild-type cells. Together, these observations uncovered borrelial antigens involved in membrane protein interactions. The study also suggests that the assembly process of OM complexes is specific and that the core or stabilizing subunits vary between complexes. Further characterization of these protein complexes including elucidation of their biological significance may shed new light on the mechanism of pathogen persistence and the development of preventative measures against the infection.     

Francisella tularensis membrane complexome by blue native/SDS-PAGE

The study of membrane proteins and membrane protein complexes (MPC) provides crucial information in the field of bacterial physiology and pathogenesis research. The method of blue native polyacrylamide gel electrophoresis and its combination with SDS-PAGE (BN/SDS-PAGE) were here employed to study the membrane complexome of an intracellular bacterium Francisella tularensis, the causative agent of a severe disease tularemia. In the presented study we describe the subunit composition and stoichiometry of several MPC involved in various cell functions (oxidative phosphorylation, membrane transport, cell division, membrane or periplasmic proteins folding, iron storage, phospholipid and cell envelope biosynthesis). Moreover, some undocumented or hypothetical MPC with possible connection to virulence factors were also proposed and some newly detected subunits were assigned to known complexes. The BN/SDS-PAGE combined with mass spectrometry appeared to be a strong tool in the investigation of membrane proteins and complexes and thus successfully complements the traditional electrophoresis approaches.     

Functional coexpression of serine protein kinase SRPK1 and its substrate ASF/SF2 in Escherichia coli

Mammalian proteins expressed in Escherichia coli are used in a variety of applications. A major drawback in producing eukaryotic proteins in E.coli is that the bacteria lack most eukaryotic post-translational modification systems, including serine/threonine protein kinase(s). Here we show that a eukaryotic protein can be phosphorylated in E.coli by simultaneous expression of a mammalian protein kinase and its substrate. We show that in bacteria expressing SRPK1, ASF/SF2 becomes phosphorylated to a degree resembling native ASF/SF2 present in interphase HeLa cell nuclei. The E.coli phosphorylated ASF/SF2 is functional in splicing and, contrary to the unphosphorylated protein, soluble under native conditions.     

Immobilization of the first dimension in 2D blue native/SDS-PAGE allows the relative quantification of membrane proteomes

In biological membranes many proteins are organized in complexes. The method of choice for the global analysis of the subunits of these complexes is two-dimensional blue native (2D BN)/SDS-PAGE. In the 1st dimension complexes are separated by BN-PAGE, and in the 2nd dimension their subunits are resolved by SDS-PAGE. In the currently available protocols the 1st dimension BN gel lanes get distorted during their transfer to the 2nd dimension separation gels. This leads to low reproducibility and high variation of 2D BN/SDS-gels, rendering them unsuitable for comparative analysis. We have developed a 2D BN/SDS-PAGE protocol where the 1st dimension BN gel is cast on a GelBond PAG film. Immobilization prevents distortion of BN gel lanes, which lowers variation and greatly improves reproducibility of 2D BN/SDS-gels. 2D BN/SDS-PAGE with an immobilized 1st dimension was used for the comparative analysis of the cytoplasmic membrane proteomes of Escherichia coli cells overexpressing a membrane protein and to create a 2D BN/SDS-PAGE reference map of the E. coli cytoplasmic membrane proteome with 143 identified proteins from 165 different protein spots.     

The E. coli monothiol glutaredoxin GrxD forms homodimeric and heterodimeric FeS cluster containing complexes

Monothiol glutaredoxins (mono-Grx) represent a highly evolutionarily conserved class of proteins present in organisms ranging from prokaryotes to humans. Mono-Grxs have been implicated in iron sulfur (FeS) cluster biosynthesis as potential scaffold proteins and in iron homeostasis via an FeS-containing complex with Fra2p (homologue of E. coli BolA) in yeast and are linked to signal transduction in mammalian systems. However, the function of the mono-Grx in prokaryotes and the nature of an interaction with BolA-like proteins have not been established. Recent genome-wide screens for E. coli genetic interactions reported the synthetic lethality (combination of mutations leading to cell death; mutation of only one of these genes does not) of a grxD mutation when combined with strains defective in FeS cluster biosynthesis (isc operon) functions [Butland, G., et al. (2008) Nature Methods 5, 789-795]. These data connected the only E. coli mono-Grx, GrxD to a potential role in FeS cluster biosynthesis. We investigated GrxD to uncover the molecular basis of this synthetic lethality and observed that GrxD can form FeS-bound homodimeric and BolA containing heterodimeric complexes. These complexes display substantially different spectroscopic and functional properties, including the ability to act as scaffold proteins for intact FeS cluster transfer to the model [2Fe-2S] acceptor protein E. coli apo-ferredoxin (Fdx), with the homodimer being significantly more efficient. In this work, we functionally dissect the potential cellular roles of GrxD as a component of both homodimeric and heterodimeric complexes to ultimately uncover if either of these complexes performs functions linked to FeS cluster biosynthesis.     

Multiplexed complexome profiling using tandem mass tags

Complexome profiling is a rapidly spreading, powerful technique to gain insight into the nature of protein complexes. It identifies and quantifies protein complexes separated into multiple fractions of increasing molecular mass using mass spectrometry-based, label-free bottom-up proteomics. Complexome profiling enables a sophisticated and thorough characterization of the composition, molecular mass, assembly, and interactions of protein complexes. However, in practice, its application is limited by the large number of samples it generates and the related time of mass spectrometry analyses. Here, we report an improved process workflow that implements tandem mass tags for multiplexing complexome profiling. This workflow substantially reduces the number of samples and measuring time without compromising protein identification or quantification reliability. In profiles from mitochondrial fractions of cells recovering from chloramphenicol treatment, tandem mass tags-multiplexed complexome profiling exhibited migration patterns of mature ATP synthase (complex V) and assembly intermediates that were consistent in composition and abundance with profiles obtained by the label-free approach. Reporter ion quantifications of proteins and complexes unaffected by the chloramphenicol treatment presented less variation in comparison to the label-free method. Incorporation of tandem mass tags enabled an efficient and robust complexome profiling analysis and may foster broader application for protein complex profiling in biomedical research and diagnostics.     

YdgT, the Hha paralogue in Escherichia coli, forms heteromeric complexes with H-NS and StpA

In enteric bacteria, proteins of the Hha/YmoA family play a role in the regulation of gene expression in response to environmental factors. Interaction of both Hha and YmoA with H-NS has been reported, and an Hha/H-NS complex has been shown to modulate expression in Escherichia coli of the haemolysin operon of plasmid pHly152. In addition to the hns gene, the chromosome of E. coli and other enteric bacteria also includes the stpA gene that encodes the StpA protein, an H-NS paralogue. We report here the identification of the Hha paralogue in E. coli, the YdgT protein. As Hha paralogue, YdgT appears to fulfil some of the functions reported for StpA as H-NS paralogue: YdgT is overexpressed in hha mutants and can compensate, at least partially, some of the hha-induced phenotypes. We also demonstrate that YdgT interacts both with H-NS and with StpA. Protein cross-linking studies showed that YdgT/H-NS heteromeric complexes are generated within the bacterial cell. The StpA protein, which is subjected to Lon-mediated turnover, was less stable in the absence of Hha or YdgT. Our findings suggest that Hha, YdgT and StpA may form complexes in vivo.     

Systematic identification of the subproteome of Escherichia coli cell envelope reveals the interaction network of membrane proteins and membrane-associated peripheral proteins

Membrane proteins of Gram-negative bacteria are key molecules that interface the cells with the environment. Despite recent proteomic identification of numerous oligomer proteins in the Escherichia coli cell envelope, the protein complex of E. coli membrane proteins and their peripherally associated proteins remain ill-defined. In the current study, we systematically analyze the subproteome of E. coli cell envelope enriched in sarcosine-insoluble fraction (SIF) and sarcosine-soluble fraction (SSF) by using proteomic methodologies. One hundred and four proteins out of 184 spots on 2D electrophoresis gels are identified, which includes 31 outer membrane proteins (OMPs). Importantly, our further proteomic studies reveal a number of previously unrecognized membrane-interacting protein complexes, such as the complex consisting of OmpW and fumarate reductase. This established complete proteomic profile of E. coli envelope also sheds new insight into the function(s) of E. coli outer envelope.     

Structural models of the supramolecular organization of AQP0 and connexons in junctional microdomains

Membrane proteins perform many essential cellular functions. Over the last years, substantial advances have been made in our understanding of the structure and function of isolated membrane proteins. However, like soluble proteins, many membrane proteins assemble into supramolecular complexes that perform specific functions in specialized membrane domains. Since supramolecular complexes of membrane proteins are difficult to study by conventional approaches, little is known about their composition, organization and assembly. The high signal-to-noise ratio of the images that can be obtained with an atomic force microscope (AFM) makes this instrument a powerful tool to image membrane protein complexes within native membranes. Recently, we have reported high-resolution topographs of junctional microdomains in native eye lens membranes containing two-dimensional (2D) arrays of aquaporin-0 (AQP0) surrounded by connexons. While both proteins are involved in cell adhesion, AQP0 is a specific water channel whereas connexons form cell–cell communication channels with broad substrate specificity. Here, we have performed a detailed analysis of the supramolecular organization of AQP0 tetramers and connexon hexamers in junctional microdomains in the native lens membrane. We present first structural models of these junctional microdomains, which we generated by docking atomic models of AQP0 and connexons into the AFM topographs. The AQP0 2D arrays in the native membrane show the same molecular packing of tetramers seen in highly ordered double-layered 2D crystals obtained through reconstitution of purified AQP0. In contrast, the connexons that surround the AQP0 arrays are only loosely packed. Based on our AFM observations, we propose a mechanism that may explain the supramolecular organization of AQP0 and connexons in junctional domains in native lens membranes.     

TMEM70 functions in the assembly of complexes I and V

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

Systematic analysis of native membrane protein complexes in Escherichia coli

The cell envelope of Escherichia coli is an essential structure that modulates exchanges between the cell and the extra-cellular milieu. Previous proteomic analyses have suggested that it contains a significant number of proteins with no annotated function. To gain insight into these proteins and the general organization of the cell envelope proteome, we have carried out a systematic analysis of native membrane protein complexes. We have identified 30 membrane protein complexes (6 of which are novel) and present reference maps that can be used for cell envelope profiling. In one instance, we identified a protein with no annotated function (YfgM) in a complex with a well-characterized periplasmic chaperone (PpiD). Using the guilt by association principle, we suggest that YfgM is also part of the periplasmic chaperone network. The approach we present circumvents the need for engineering of tags and protein overexpression. It is applicable for the analysis of membrane protein complexes in any organism and will be particularly useful for less-characterized organisms where conventional strategies that require protein engineering (i.e., 2-hybrid based approaches and TAP-tagging) are not feasible.     

Quantitative Brightness Analysis of Fluorescence Intensity Fluctuations in E. Coli

The brightness measured by fluorescence fluctuation spectroscopy specifies the average stoichiometry of a labeled protein in a sample. Here we extended brightness analysis, which has been mainly applied in eukaryotic cells, to prokaryotic cells with E. coli serving as a model system. The small size of the E. coli cell introduces unique challenges for applying brightness analysis that are addressed in this work. Photobleaching leads to a depletion of fluorophores and a reduction of the brightness of protein complexes. In addition, the E. coli cell and the point spread function of the instrument only partially overlap, which influences intensity fluctuations. To address these challenges we developed MSQ analysis, which is based on the mean Q-value of segmented photon count data, and combined it with the analysis of axial scans through the E. coli cell. The MSQ method recovers brightness, concentration, and diffusion time of soluble proteins in E. coli. We applied MSQ to measure the brightness of EGFP in E. coli and compared it to solution measurements. We further used MSQ analysis to determine the oligomeric state of nuclear transport factor 2 labeled with EGFP expressed in E. coli cells. The results obtained demonstrate the feasibility of quantifying the stoichiometry of proteins by brightness analysis in a prokaryotic cell.     

The steroid-binding properties of recombinant glucocorticoid receptor: a putative role for heat shock protein hsp90

The steroid-binding domain of the human glucocorticoid receptor was expressed in Escherichia coli either as a fusion protein with protein A or under control of the T7 RNA polymerase promoter. The recombinant proteins were found to bind steroids with the normal specificity for a glucocorticoid receptor but with reduced affinity (Kd for triamcinolone acetonide approximately 70 nM). Glycerol gradient analysis of the E. coli lystate containing the recombinant protein indicated no interaction between the glucocorticoid receptor fragment and heat shock proteins. However, synthesis of the corresponding fragments of glucocorticoid receptor in vitro using rabbit reticulocyte lystate resulted in the formation of proteins that bound triamcinolone acetonide with high affinity (Kd 2nM). Glycerol gradient analysis of these proteins, with and without molybdate, indicated that the in vitro synthesised receptor fragments formed complexes with hsp90 as previously shown for the full-length rat glucocorticoid receptor. Radiosequence analysis of the recombinant steroid-binding domain expressed in E. coli and affinity labelled with dexamethasone mesylate identified binding of the steroid to Cys-638 predominantly. However, all cysteine residues within the steroid-binding domain were affinity labelled to a certain degree indicating that the recombinant protein has a structure similar to the native receptor but more open and accessible.     

Proteome-wide analysis of chaperonin-dependent protein folding in Escherichia coli

The E. coli chaperonin GroEL and its cofactor GroES promote protein folding by sequestering nonnative polypeptides in a cage-like structure. Here we define the contribution of this system to protein folding across the entire E. coli proteome. Approximately 250 different proteins interact with GroEL, but most of these can utilize either GroEL or the upstream chaperones trigger factor (TF) and DnaK for folding. Obligate GroEL-dependence is limited to only approximately 85 substrates, including 13 essential proteins, and occupying more than 75% of GroEL capacity. These proteins appear to populate kinetically trapped intermediates during folding; they are stabilized by TF/DnaK against aggregation but reach native state only upon transfer to GroEL/GroES. Interestingly, substantially enriched among the GroEL substrates are proteins with (betaalpha)8 TIM-barrel domains. We suggest that the chaperonin system may have facilitated the evolution of this fold into a versatile platform for the implementation of numerous enzymatic functions.