Bcl2-associated Athanogene 3 Interactome Analysis Reveals a New Role in Modulating Proteasome Activity*

Bcl2-associated athanogene 3 (BAG3), a member of the BAG family of co-chaperones, plays a critical role in regulating apoptosis, development, cell motility, autophagy, and tumor metastasis and in mediating cell adaptive responses to stressful stimuli. BAG3 carries a BAG domain, a WW domain, and a proline-rich repeat (PXXP), all of which mediate binding to different partners. To elucidate BAG3's interaction network at the molecular level, we employed quantitative immunoprecipitation combined with knockdown and human proteome microarrays to comprehensively profile the BAG3 interactome in humans. We identified a total of 382 BAG3-interacting proteins with diverse functions, including transferase activity, nucleic acid binding, transcription factors, proteases, and chaperones, suggesting that BAG3 is a critical regulator of diverse cellular functions. In addition, we characterized interactions between BAG3 and some of its newly identified partners in greater detail. In particular, bioinformatic analysis revealed that the BAG3 interactome is strongly enriched in proteins functioning within the proteasome-ubiquitination process and that compose the proteasome complex itself, suggesting that a critical biological function of BAG3 is associated with the proteasome. Functional studies demonstrated that BAG3 indeed interacts with the proteasome and modulates its activity, sustaining cell survival and underlying resistance to therapy through the down-modulation of apoptosis. Taken as a whole, this study expands our knowledge of the BAG3 interactome, provides a valuable resource for understanding how BAG3 affects different cellular functions, and demonstrates that biologically relevant data can be harvested using this kind of integrated approach.

Bcl2-associated athanogene 3 (BAG3) 1 is a member of the human BAG family of molecular co-chaperone proteins (1). BAG proteins share a common C-terminal BAG domain that mediates interactions with the ATPase domain of heat shock protein (Hsp) 70 (1), but they have divergent N termini that confer different properties and functions to each member (2). For example, in addition to the BAG domain, BAG3 contains a proline-rich (PXXP) region and an N-terminal WW domain (1) that mediate binding to partners other than Hsp70. Moreover, in addition to the BAG domain, BAG3 possesses one of the longest C-terminal sequences among the known BAG proteins. These multifaceted structural characteristics underlie the ability of BAG3 to modulate biological processes such as apoptosis, development, cytoskeleton organization, and autophagy and to mediate cell adaptive responses to stressful stimuli (3).
As a chaperone, BAG3 seems to influence cell survival by interacting with different molecular partners and thus activating multiple pathways. For example, BAG3 can potently suppress cell apoptosis via its interaction with Hsp70, a chaperone protein able to modulate apoptosis by interfering with cytochrome c release, apoptosome formation, and other events in the death process (4,5). In addition, BAG3 binds PLC-gamma (6) and Bcl-2, with which it synergizes to prevent cell death (7). Other cellular signaling molecules reported to be regulated by BAG3 include Raf-1, CDK-4, EGFR, MMP-2, Hsps B8 and B6, and focal adhesion kinase (3). Notably, recent reports have shown additional functions of this chaperone in the regulation of autophagy (8,9), virus replication (10,11), and NF-B activity (12). Moreover, in a previous study, we demonstrated that the up-regulation of BAG3 may be part of the cellular response that SARS-CoV relies on for replication (13).
Despite the multifunctional roles of this stress-and survivalrelated protein, little is known about the molecular mechanisms underlying these functions. Given that most proteins are part of one or more protein complexes, knowledge of protein-protein interactions is essential in order for one to understand the nature of protein-mediated biological processes. Some of the interacting partners of BAG3 are known, but the diverse functions identified for BAG3 suggest that this protein might have an even greater assortment of binding partners. Thus, the identification and characterization of the entire set of BAG3 interaction partners is critical for a comprehensive understanding of its biological functions. To this end, we have performed a global analysis of BAG3 interaction partners using the quantitative immunoprecipitation combined with knockdown (QUICK) and proteome microarray methods.
The QUICK method, recently developed by Mann and colleagues (14), is a relatively new method for identifying interactions between proteins at endogenous levels under physiological conditions via a combination of stable isotope labeling with amino acids in cell culture (SILAC), RNA interference (RNAi)-induced knockdown, co-immunoprecipitation, and quantitative MS. This highly sensitive and accurate approach for protein-protein interaction analysis has been applied to identify the interaction partners of ␤-catenin and Cbl (14), 14-3-3 (15), Lrrk2 (16), and Stat3 (17).
Proteome microarrays, also known as protein chips, are constructed by spotting hundreds to thousands of individually purified proteins at a high density on a solid surface (18,19). A major breakthrough in this technology came from a study by Zhu et al. in which a proteome-wide protein array consisting of 5,800 unique yeast proteins was constructed and used to identify calmodulin-and phospholipid-binding proteins (20). Since then, proteome microarrays have been employed to probe various types of protein biochemical properties, enabling high-throughput screening of thousands of proteins in a single experiment (21)(22)(23). Recently, a human proteome microarray comprising 17,000 purified, full-length proteins was constructed and used to identify highly specific monoclonal antibodies (24). Proteome microarrays provide information about direct biochemical and physical interactions among biomolecules and have been successfully used for the identification of protein interaction partners (20,25,26). They have become one of the most powerful multiplexed detection platforms for addressing diverse biological questions in a highthroughput fashion.
Although both QUICK and proteome microarrays have been successfully used for the identification of protein-protein interactions, each approach detects and reports these interactions in a distinct manner. For example, the proteome microarrays detect binary interactions, whereas the QUICK method purifies and identifies protein complexes. The QUICK method uses untagged endogenous bait proteins to pull down prey proteins followed by mass spectrometric detection of proteins, whereas the proteome microarray method relies on separate protein complementation schemes to ultimately report on whether a protein pair is interacting. In addition, the QUICK method identifies interactions at approximately physiological cellular protein concentrations, whereas the concentrations of the interacting partners in the proteome microarray screens are not necessarily comparable to those found in the native environment. However, the protein concentrations used are less diverse on the proteome microarray, and this may facilitate the identification of BAG3 interacting proteins expressed at low levels under physiological conditions. The QUICK method is presently not a high-throughput method, whereas the proteome microarray is an outstandingly flexible tool for high-throughput studies. Therefore, although each approach has distinct advantages, disadvantages, and limitations, the combination of these complementary tools provides a powerful methodology for proteomic investigations.
In this study, we employed both QUICK and proteome microarray methods to decipher the human BAG3 interactome. Our results confirm that BAG3 indeed interacts with many more proteins than previously thought and suggest a number of different functions for BAG3, which can now be further explored. This study also demonstrates that the complementary information obtained by combining the QUICK method and proteome microarrays can provide a comprehensive portrait of the protein-protein interactions with any protein of interest.

EXPERIMENTAL PROCEDURES
Cell Culture-Cells from human cervical cancer cell line HeLa, lung cancer cell line A549, and multiple myeloma cell line U266 were purchased from the American Type Culture Collection (Rockville, MD). HeLa and A549 cells were routinely maintained in Dulbecco modified Eagle's medium (DMEM) (Sigma-Aldrich, Saint Louis, MO) supplemented with 1% penicillin/streptomycin and 10% fetal bovine serum (FBS) (Sigma-Aldrich). All cells were grown at 37°C in a humidified incubator with an atmosphere of 5% CO 2.
QUICK Identification of BAG3 Interacting Proteins SILAC Labeling and RNA Interference-The BAG3 interaction proteins in HeLa cells were identified using the QUICK approach. SILAC labeling was performed essentially as described elsewhere (15). A SILAC Protein Quantitation Kit (Pierce Biotechnology, Rockford, IL) was used according to the manufacturer's instructions. In brief, HeLa cells were grown in SILAC DMEM (Pierce) containing 10% v/v dialyzed FBS and 0.1 mg/ml heavy [ 13 C 6 ] or light [ 12 C 6 ] L-lysine. To ensure full incorporation of heavy-and light-labeled amino acids, cells were grown for at least six cell doublings prior to treatment. For BAG3 RNA interference, small interfering RNAs (siRNAs) specific for BAG3 (siBAG3: 5Ј-GCAUGCCAGAAACCACUCA-3Ј) or negative control siRNA (siNC: 5Ј-UUCUCCGAACGUGUCACGUTT-3Ј) were purchased from Shanghai GenePharma (Shanghai, China). Heavy [ 13 C 6 ]-labeled HeLa cells were transfected with siBAG3 using Lipofectamine 2000 (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. A control experiment was performed by transfecting light [ 12 C 6 ]-labeled HeLa cells with negative control siNC. A549 cells were transfected using the same protocol. Changes in BAG3 protein levels were determined via Western blotting.
In Vivo Formaldehyde Cross-linking-72 h post-transfection, cells were harvested, washed with PBS, and incubated in 1.0% w/v paraformaldehyde in PBS for 10 min at 37°C. To stop the cross-linking reaction, 1.25 M glycine was added to a final concentration of 125 mM for 5 min at room temperature. Cells were harvested, washed three times with ice-cold washing buffer (10 M Tris-HCl, 250 M sucrose, pH 7.0), and transferred to a clean 1.5-ml Eppendorf tube. Cells were lysed with RIPA lysis buffer, and cellular debris was removed via centrifugation for 30 min at 13,200 ϫ g at 4°C. Supernatants were collected, and protein concentrations were measured in duplicate using a BCA protein assay kit (Pierce).
Antibody Coupling and Immunoprecipitation-In order to prevent interference from an excess of antibodies with mass spectrometric detection, the BAG3 antibody was coupled to Protein A/G Sepharose (GE Healthcare). Protein A/G Sepharose was mixed with antibody in 0.1 M borate buffer, pH 8.2, for 30 min at room temperature with gentle shaking, after which the Sepharose beads were washed with excess borate buffer. The Sepharose was washed with 0.2 M triethanolamine, pH 8.2, and then resuspended in 20 volumes of dimethyl pimelimidate dihydrochloride (Pierce) freshly made up in 0.2 M triethanolamine, readjusting the pH to 8.2. The mixture was agitated gently at room temperature for 1 h, and the reaction was stopped by centrifuging the beads (500 ϫ g for 1 min) and resuspending in an equal volume of ethanolamine (pH 8.2) of the same molarity as the dimethyl pimelimidate dihydrochloride. After 5 min, the cross-linked beads were centrifuged and washed three times with lysis buffer. Prior to immunoprecipitation, 2 mg of protein from light siNC transfected HeLa cells was mixed with 2 mg of protein from heavy siBAG3 transfected HeLa cells and precleared with 0.2 ml Protein A/G Sepharose for 2 h. Protein A-Sepharose beads coupled with antibodies were mixed with precleared mixtures and rotated for 18 h at 4°C. The beads were washed twice with 2 ml of lysis buffer containing 0.1% Triton X-100 and twice with 2 ml of lysis buffer without detergent. After complete removal of the supernatant with a microliter syringe, immunoprecipitated proteins were eluted with 3 ϫ 200 l of 5% acetic acid pH 3 and dried in a speed-vacuum. Laemlli sample buffer, pH 7.4 (10 mM Tris, 10% v/v glycerol, 2% w/v SDS, 5 mM EDTA, 0.02% bromphenol blue, and 6% ␤-mercaptoethanol), was added and boiled for 20 min at 95°C. Proteins were then separated on a 10% SDS-PAGE gel and visualized with Coomassie Blue staining.
Protein Separation and In-gel Digestion-Protein bands were excised from the SDS-PAGE gel and cut into 15 sections for in-gel tryptic digestion. The excised sections were chopped into small particles, washed in water, and completely destained using 100 mM ammonium bicarbonate in 50% acetonitrile. A reduction step was performed by adding 100 l of 10 mM DTT at 37°C for 1 h. The proteins were alkylated by the addition of 100 l of 50 mM iodoacetamide and allowed to react in the dark at room temperature for 1 h. Gel sections were washed first in water and then in acetonitrile, and finally they were dried using a SpeedVac (Thermo Fisher Scientific, Waltham, MA) for 30 min. Digestion was carried out using 20 g/ml sequencing-grade modified trypsin (Promega) in 50 mM ammonium bicarbonate. Sufficient trypsin solution was added to swell the gel pieces, which were kept at 4°C for 45 min and then incubated at 37°C overnight. The supernatants were transferred into a 200-l microcentrifuge tube, and the gels were extracted once with extraction buffer (67% acetonitrile containing 2.5% trifluoroacetic acid). The peptide extract and the supernatant of the gel slices were combined and then completely dried in the SpeedVac.
Mass Spectrometry and Protein Identification and Quantification-Dried peptides were reconstituted in 5% acetonitrile/0.1% formic acid. LC-MS/MS analyses were performed on an LTQ-Orbitrap mass spectrometer (Thermo Fisher) equipped with an Agilent 1100 highperformance liquid chromatography (HPLC) pump (Agilent Technologies) and a Famos autosampler (LC Packings). The peptide mixtures were loaded onto a C18 column (100 m inner diameter, 10 cm long, 5 m resin, from Michrom Bioresources, Auburn, CA) using an autosampler. Peptides were eluted with a 0%-35% gradient (Buffer A, 0.1% formic acid and 5% acetonitrile; Buffer B, 0.1% formic acid and 95% acetonitrile) over 90 min and detected online in an LTQ-Orbitrap using a data-dependent method (27). All samples were analyzed in triplicate. The general mass spectrometric conditions were as follows: spray voltage, 1.80 kV; no sheath and auxiliary gas flow; ion transfer tube temperature, 200°C. Ion selection thresholds were 1,000 counts for MS 2 and 500 counts for MS 3 . An activation q ϭ 0.25 and an activation time of 30 ms were applied in MS 2 acquisitions. The mass spectrometers were operated in positive ion mode, employing a data-dependent automatic switch between MS and MS 2 acquisition modes. For each cycle, one full MS scan of the five most intense ions in the Orbitrap at 1 ϫ 10 6 AGC target was followed by five MS 2 scans in the LTQ at 5,000 AGC target. Selected ions were excluded from further selection for 90 s. The maximum ion accumulation time was 500 ms for full MS scans and 100 ms for MS 2 scans. All MS/MS spectra were collected using normalized collision energy (a setting of 35%), an isolation window of 3 m/z, and one micro-scan. The application of mass spectrometer scan functions and HPLC solvent gradients was controlled by an XCalibur data system (Thermo Fisher Scientific, Waltham, MA). The acquired raw files were processed using BioWorks 3.3.1 (Thermo Finnigan, San Jose, CA), and the MS/MS spectra were searched using the SEQUEST search engine (version 28.13, Thermo Fisher Scientific, Waltham, MA) against a catenated database combining the International Protein Index human database (version 3.63, 84,118 sequences), the reversed sequences of all sequences in the International Protein Index human database, and 262 commonly observed contaminants. The following search criteria were employed: full tryptic specificity was required; two missed cleavages were allowed; parent ion mass tolerance ϭ 20 ppm; mass tolerance for fragment ion ϭ 0.5 Da; minimum peptide length of six residues; Cys (ϩ57.0215 Da, carbamidomethylation) set as a fixed modification; Met (ϩ15.9949 Da, oxidation) and Lys (ϩ6.0201 Da, SILAC heavy amino acid) considered as variable modifications. Peptides were accepted if their Xcorr was greater than or equal to 1.9, 2.5, 3.5, or 4.5 for ϩ1, ϩ2, ϩ3, and ϩ4 charged ions, respectively, and if their ⌬Cn was at least 0.1. The global false discovery rate for both peptides and proteins was set to 0.01. If the peptide charge and peptide sequence were the same, we selected only the peptide with the highest Xcorr value. The selected peptides were converted to dtaselect-filter.txt by the dtaselect-filter file format and quantified by Census (28). Only proteins with a minimum of two quantifiable peptides were included in our final dataset. The quantification was based on three independent QUICK experiments, and only proteins identified and quantified in at least two of three experiments were considered for further analysis. Census normalized protein ratios were utilized for the downstream analyses, and a p value of 0.05 was selected as the threshold for significant enrichment of BAG3 interaction proteins. Only those proteins with p Յ 0.05 were considered as BAG3-interacting proteins.
Proteome Microarray-based BAG3 Interactome Screen-Human Proteome Microarray Construction-N-terminal GSTtagged human proteins (ϳ17,000) were purified from 6 ml yeast culture as described elsewhere (24). Briefly, yeast transformants were grown at 30°C and inoculated into SC-ura/raffinose media. The cultures were grown overnight and induced at an A 600 of 0.6 to 0.8 to a final concentration of 2% galactose. Cells were then harvested via centrifugation after 6 h. The cells were lysed via shaking with glass beads, and the supernatant was incubated with glutathione-Sepharose (GE Healthcare). The beads were washed and the proteins were eluted with 40 mM glutathione. The purified proteins were aliquoted into 384-well plates and printed onto Fullmoon slides (Fullmoon Biosystems, Sunnyvale, CA) using a 48-pin contact microarrayer (ChipWriter Pro, Bio-Rad, Hercules, CA). Each protein was spotted in duplicate.
Probe BAG3 on Human Proteome Microarray-Proteome microarrays were blocked with blocking buffer (1% BSA in 0.1% Tween 20, TBST) for 1 h at room temperature with gentle agitation. BAG3 protein (Abnova, Taipei, Taiwan) was diluted to 0.01 g/l in blocking buffer and incubated on the blocked proteome microarray at room temperature for 1 h. The microarrays were washed three times for 5 min each time with TBST, incubated with 1:1,000 diluted BAG3 antibody (Abcam Inc., Cambridge, MA) in blocking buffer for 1 h at room temperature, and then washed three more times (5-min washes) in TBST. The microarrays were further incubated with a Cy5 conjugated second antibody at 1:1,000 dilution (Jackson ImmunoResearch Laboratories, Inc., West Grove, PA) for 1 h at room temperature and underwent three more 5-min washes. The microarrays were spun dry at 1,500 rpm for 3 min and subjected to scanning with a Genepix 4200A (Axon Instruments, Sunnyvale, CA) in order for results to be visualized and recorded. A GenePix Pro 6.0 was used for data analysis.
Proteome Microarray Data Processing-The data were extracted by GenePix Pro 6.0 from the microarray images. To generate the candidate list of BAG3 binding proteins, the signal-to-noise ratios (SNR ϭ F635 median/B635 median) were first calculated for all the spots. The SNR of a protein was averaged from two duplicated spots. The SNRs of all the proteins were calculated for microarrays with or without BAG3 as SNR(ϩ) or SNR(Ϫ), respectively. Another index, the fold change, was defined as (Signal(ϩ) Ϫ Background(ϩ))/(Signal(Ϫ) Ϫ Background(Ϫ)) and was generated for all the proteins. To determine the final candidates, the cutoff was set as SNR(ϩ) Ն 2 and fold change Ն 3.
Western Blotting Analysis-Protein extracts (50 g) prepared with RIPA lysis buffer were resolved on a 10% SDS-PAGE gel and transferred to an Immobilon-P PVDF transfer membrane (Millipore, Bedford, MA) via electroblotting. After blocking with 5% nonfat milk, the membranes were probed with rabbit monoclonal anti-Bag1, goat polyclonal anti-Bag2, rabbit monoclonal anti-BAG3, rabbit polyclonal anti-Bag4, rabbit polyclonal anti-Bag5 (Abcam Inc., Cambridge, MA), mouse anti-Actin monoclonal (Santa Cruz Biotechnology, Santa Cruz, CA), and mouse anti-ubiquitin monoclonal antibodies (P4D1, Cell Signaling). Blots were then incubated with a streptavidin alkaline phosphatase-conjugate secondary antibody (Amersham Biosciences, Piscataway, NJ) for 1 h at room temperature. The signal was visualized with ECL detection reagent (GE Healthcare).
Bioinformatics Analysis of the BAG3 Interactome-The BAG3 interactome was classified based on the PANTHER (Protein Analysis through Evolutionary Relationships) system, which is a unique resource that classifies genes and proteins according to their functions (29). Functional diversity among proteins present in two datasets was analyzed as previously described (30). The relative distribution of all proteins from Homo sapiens, as annotated according to the PANTHER functional classes, is shown by the line labeled "Expected" in Fig. 3. It denotes the relative fraction (coverage) of all human proteins that belong to a given category. The deviations of the QUICK and proteome microarray datasets from this distribution are shown by different-colored bars and values in Fig. 3.
To determine whether any types of proteins were overrepresented, enrichment analysis of Gene Ontology (GO) terms (31), Kyoto Encyclopedia of Genes and Genomes pathways (32), and Pfam domain families (33) was performed using the Web-accessible program DAVID 6.7 (34,35). The default human proteome was used as the background list. The significance of the enrichments was statistically evaluated with a modified Fisher's exact test (EASE score), and a p value for each term was calculated by applying a Benjamini-Hochberg false discovery rate correction (34,35). For GO term enrichments, the GO FAT annotation available in DAVID was used. GO FAT is a subset of the GO term set created by filtering out the broadest ontology terms so as to not overshadow more specific ones. The enrichment of GO biological process terms was also analyzed using Cytoscape and its Plugin BiNGO 2.3 (36), using the complete GO term set and a hypergeometric statistical test with Benjamini-Hochberg false discovery rate correction. The GO categories GO SLim Generic assignment, the distribution of cellular components, molecular functions, and biological processes of the BAG3 interactome were analyzed.
The BAG3 interactome was searched against the STRING database, version 8 (37), for protein-protein interactions. Only interactions between proteins belonging to the BAG3 interactome dataset were selected. STRING defines a metric called a "confidence score" to define interaction confidence; we fetched all interactions for our BAG3 interactome dataset that had a confidence score of Ն0.4 (medium confidence). A network of protein-protein interactions was generated in this manner and then visualized by Cytoscape v2.8.1. This BAG3 interaction network was further analyzed for densely connected regions using a graph theoretic clustering algorithm, "Molecular Complex Detection" (MCODE) (38). MCODE is part of the plug-in toolkit in Cytoscape (39). The top four high-ranking modules were selected for further analysis and rendering.
Surface Plasmon Resonance Analysis-Surface plasmon resonance (SPR) experiments were carried out at 25°C using a BIAcore 3000 instrument (BIAcore AB, Uppsala, Sweden). SPR buffers, regeneration solutions, and sensor chips were purchased from GE Healthcare. Chemicals were from Sigma (St. Louis, MO) unless otherwise stated. Purchased proteins were BAG3 (Abnova, Taipei, Taiwan), 14-3-3, Hsp70, MAPK7, ATG3, PSMB1, HSF1, and ENO1 (OriGene Technologies, Rockville, MD). The surface of the CM5 chip was activated following a standard 1-ethyl-3(3-dimethylaminopropyl)-carbodiimide hydrochloride/N-hydroxyl-succinimide amine coupling Biacore protocol. For the immobilization of proteins, the BAG3 protein was injected in 100 l of 10 mM sodium acetate (flow rate: 10 l/min). After immobilization, each surface was blocked by 1 M ethanolamine at pH 8.5 for 6 min. For interaction measurements, recombinant human proteins were dissolved according to the manufacturer's instructions (OriGene Technologies, Rockville, MD). Various concentrations of proteins were injected over this interaction surface for 5 min in a buffer containing 150 mM NaCl, 5 mM CaCl 2 , 0.005% (v/v) Tween 20, and 20 mM Hepes (pH 7.4) with a flow rate of 5 l/min. The sensor surface was regenerated between sample injections by washing with 6 M guanidine-HCl for 0.5 min at a flow rate of 50 l/min. An equivalent volume of each protein sample (analyte) was injected over the chip surface with no protein immobilized to serve as a blank phase for background subtraction. ENO1, which was identified as a background protein in this study, was injected as a negative control using the same conditions. All experiments were repeated at least three times. The kinetic parameters (k a , association rate constant; k d , dissociation rate constant; K D ϭ k d /k a ϭ equilibrium dissociation constant) for each interaction were determined by globally fitting the experimental data using BIAevaluation software 4.1 (GE Healthcare).
Plasmid Constructs and Transient Transfection-A cDNA encoding human BAG3 was generated via PCR from a human brain cDNA library (Invitrogen, Carlsbad, CA) and subcloned into the EcoRI/XhoI sites of the eukaryotic expression plasmid pcDNA3 tagged with a Flag epitope. HeLa or A549 cells were transfected with Lipofectamine 2000 (Invitrogen) according to the manufacturer's instructions.
Proteasome Activity Assay-Cytosolic extracts (without protease inhibitors) were used to measure proteasome activity using a 20S proteasome assay kit (Chemicon, Temecula, CA) following the manufacturer's instructions. The proteasome has three distinct peptidase activities: chymotryptic, tryptic, and caspase-like activities (40). The specific proteasome activities were measured as described previously, with a few modifications (41,42). Treated or untreated HeLa and A549 cells were washed with PBS and pelleted via centrifugation. Cell pellets were lysed via repeated freeze-thaw cycles in 0.25 M sucrose, 25 mM Hepes (pH 7.8), 10 mM MgCl 2 , 1 mM EDTA, and 1 mM dithiothreitol. Lysates were centrifuged at 12,000 ϫ g for 10 min, and supernatants were collected. Protein concentrations were measured in duplicate using a BCA protein assay kit (Pierce). Cell lysates were diluted with proteolysis assay buffer (50 mM Hepes (pH 7.8), 20 mM KCl, 5 mM MgCl 2 , and 1 mM dithiothreitol) to a final concentration of 0.5 mg/ml. The specific proteasome activity was measured via the addition of 170 l proteolysis assay buffer and 30 l Suc-LLVY-AMC (chymotryptic substrate; Biomol International, Plymouth Meeting, PA), Z-ARR-AMC (tryptic substrate; Calbiochem), or Z-LLE-AMC (caspase-like substrate; Biomol International) (1 mM stock solution in DMSO) to 100 l of the diluted cell lysate. Assays of 26S proteasome complexes were carried out in the presence of 5 mM ATP, which stimulates 26S proteasome activity. 20S proteasome activity was measured in the presence of 0.02% SDS without ATP, which activates 20S proteasomes and partially inhibits 26S proteasomes. The mixture was incubated at 37°C for 2 h. Following proteasome cleavage, the fluorescent group, 7-amido-4-methylcoumarin (AMC), was released and measured on a Perkin-Elmer fluorometer using an excitation wavelength of 380 nm and an emission wavelength of 460 nm. The data are plotted as mean Ϯ standard error and are expressed as a percentage of the control, which was considered to be 100%. All samples were assayed in three independent experiments. The statistical significance of the data was tested using Student's t test for unpaired samples. Differences were considered to be significant if p Յ 0.05 (*).
Detection of Apoptosis-Apoptosis was detected by using Annexin V/PI staining. In brief, cells were washed once in PBS and then stained with Annexin V-FITC and PI (2 mg/ml) according to the manufacturer's instructions. Samples were acquired on a FACScan flow cytometer (Becton Dickinson, San Jose, CA) and analyzed with WinMDI 2.8. All samples were assayed in three independent experiments. The statistical significance of the data was tested using Student's t test for unpaired samples. Differences were considered to be significant if p Յ 0.05 (*).

QUICK Identification of BAG3 Interacting Proteins-We
identified BAG3 interacting proteins first using the QUICK approach, a SILAC-based quantitative strategy to capture endogenous protein-protein interactions (Fig. 1A). A prerequisite for the QUICK assay is the availability of RNAi cells exhibiting reduced expression of the protein of interest. BAG3 gene silencing was achieved by transfecting HeLa cells with BAG3 siRNA (siBAG3). BAG3 knockdown efficiency was monitored by means of Western blotting using a BAG3-specific antibody. Levels of BAG3 proteins were markedly reduced in siBAG3-treated HeLa cells realative to HeLa cells transfected with negative control siRNA (siNC) (Fig. 1B). The depletion of intracellular levels of BAG3 showed no appreciable effect on the expression of the other BAG isoforms or on the internal control, actin (Fig. 1B). Proteins were immunoprecipitated, and peptides from pooled eluates were subjected to LC-MS/MS. As expected, the level of BAG3 was 4.89-fold lower in the knock-down treatment than in the control (supplemental Table S1). From among the quantified proteins, we include here only the results for those proteins that could be quantified in two or more experiments. Using the criteria described in "Experimental Procedures," we identified a total of 174 proteins as potential BAG3 interaction partners (supplemental Table S1). Twenty-five proteins had p values greater than 0.05, indicating that under our experimental conditions, they are likely to interact nonspecifically with BAG3 and were therefore considered as background proteins (supplemental Table S2). Detailed information about the quantified proteins and the ratios for each peptide are summarized in supplemental Table S2.
Identification of BAG3 Interacting Proteins Using a Human Proteome Microarray-We next used a human proteome microarray to identify putative BAG3 interacting proteins ( Fig. 2A). As shown in Fig. 2B, a protein microarray composed of ϳ17,000 individually purified human proteins as N-terminal GST fusions was incubated with affinity-purified human BAG3 protein. To detect binding events, we applied a BAG3-specific monoclonal antibody to the microarray and followed with incubation with a Cy5-labeled secondary antibody (Fig. 2B,  left panel). Because anti-BAG3 and Cy5-labeled secondary antibodies can potentially recognize proteins other than BAG3 on the human proteome microarray, we also performed a control binding experiment in the absence of BAG3 on the microarray (Fig. 2B, right panel). Using stringent criteria as described in "Experimental Procedures," 289 proteins were identified as potential BAG3 interacting proteins (supplemental Table S3), 11 of which are shown in Fig. 2C. To make an estimation of the false positive rate, we adopted a statistical model developed for estimating DNA microarray false positive rates by Wang et al. (43). Based on our experience, for proteome microarray experiments, the discordant rate of two replicated microarrays could be easily controlled under 10%.
If we assume that the discordant rate of this study was 10%, then the false positive rate is (10% Functional Diversity in the Two Datasets and Construction of the BAG3 Interactome-A comparison of the QUICK and proteome microarray results revealed that a substantial fraction of the proteins were shared between the two datasets. Of all 382 proteins identified as BAG3 interacting proteins, 81 (21.2%) were identified by both methods (Fig. 3A). Some of the known BAG3 interacting proteins, such as Hsp70 (5), sequestosome 1 (9), CCT (44), and HspB6 (45), were recovered in our datasets (supplemental Table S4), supporting the reliability of our data.
Although 20% of the proteins in the two datasets overlapped, each dataset contained a distinct group of proteins, and functional classification of these proteins permitted characterization of the different properties of the two datasets. Fig.   FIG. 2. Identification of BAG3 interacting proteins with a human proteome microarray. A, workflow of the proteome microarray strategy for the identification of protein-protein interactions. B, human proteome microarrays were probed with BAG3 (left), followed by a BAG3 antibody and a Cy5-labeled secondary antibody. A control experiment was carried out without BAG3 (right). Many stronger signals can be observed from the microarray with BAG3 than from the microarray without BAG3. Enlarged areas from the same locations of the two micorarrays with and without BAG3 are illustrated in the middle. When the signals from the two microarrays were compared, 289 candidate BAG3 interacting proteins were identified. C, representative BAG3 interacting proteins.

FIG. 3. Comparison of BAG3 interacting proteins identified via the QUICK and proteome microarray methods. A, Venn diagrams
showing the degree of overlap between the QUICK and proteome microarray data sets. B, functional diversity among proteins present in two datasets. The relative distribution of all proteins from Homo sapiens, annotated according to the PANTHER functional classes, is shown by the line labeled "Expected." This denotes the relative fraction (coverage) of all human proteins that belong to a given category. Deviations of the QUICK and proteome microarray datasets from this distribution are shown by different colored bars; for example, proteins labeled "transcription factor" are relatively underrepresented in the two datasets.
3B shows the relative distribution of functional classes for all Homo sapiens proteins in 22 PANTHER functional classes (29) (labeled "Expected") (30) and the interacting proteins from the QUICK and proteome microarray datasets. Across all 22 PANTHER functional classes, we found the greatest deviations in the QUICK dataset to be associated with an overrepresentation in the chaperone, transfer/carrier protein, and protease categories and an underrepresentation of transcription factors and receptors. In the proteome microarray dataset, adaptor and cytoskeletal proteins were the two most overrepresented categories, and transcription factors and receptors were the most underrepresented. It is noteworthy that the transcription factor and receptor categories in the proteome microarray dataset were much less underrepresented than in QUICK (Fig. 3B). Therefore, there are large, clear differences between the two datasets, particularly for proteins with chaperone, transcription factor, or receptor functions. A plausible explanation for this difference is that transcription factors and receptors are generally expressed at low levels (46), making them hard to identify with QUICK, because the affinity purification step will most likely result in a loss of these proteins. However, this is not a problem with the protein microarray because each protein is overexpressed in yeast and individually purified. Similarly, proteins that function in organized assemblies, such as chaperones, would be preferentially included in the QUICK dataset, whereas the binary binding assay performed on the proteome microarrays is not designed to identify interactions involving more than two proteins. In order to take full advantage of the complementary nature of data obtained via the QUICK and proteome microarray methods, we combined the proteins identified by each method, giving a nonredundant list of 382 human proteins. This list, referred to as the BAG3 interactome (supplemental Table S4), was used in all subsequent analyses.
Categorization and GO Analysis of the BAG3 Interactome-To evaluate the biological relevance of the BAG3 interactome, the PANTHER classification system was used to classify these interacting proteins according to their functions and biological processes. The PANTHER classification analysis revealed that the interactome can be classified into 17 groups according to biological processes and 23 groups according to molecular functions (supplemental Fig. S1). The largest group within the BAG3 interactome is the transferase family (12%). Other significant functional categories include nucleic acid binding (10%), transcription factors (8%), enzyme modulators (8%), proteases (7%), and chaperones (6%). These results confirm that BAG3 is a critical regulator of diverse cellular functions across widespread biological processes (3). It should be noted that many of these proteins are multifunctional and are assigned to more than one functional class.
We performed a GO analysis to gain insights into the functional roles of BAG3, identifying enriched terms in the BAG3 interactome. First, the assignment of GO categories gave us an overview of the GO distribution (supplemental Fig. S2). Next, we performed GO biological process, molecular function, and cellular component analysis (Table I and  supplemental Tables S5-S7). Consistent with the fact that BAG3 is a molecular chaperone regulator, 47 proteins were annotated in the "catabolic processes" term (GO term: 9057, p ϭ 8.33 ϫ 10 Ϫ12 ). The biological process terms are presented as a graph in which direct links describe the hier- archy and relationships between terms (supplemental Fig.  S3A). Colored nodes are those that are significantly overrepresented. Examination of the closest branch points of the overrepresented GO biological processes indicated that four major functional groups were strongly enriched in the BAG3 interactome, including two clusters of biological processes involved in the control of the proteasome-ubiquitination system and metabolic processes and two other clusters related to the regulation of cellular processes such as cell cycle and cell death (supplemental Fig. S3A). In the GO molecular function category, we found that the most overrepresented functions were involved in regulating enzyme activity, such as threonine-type endopeptidases and peptidases, glutathione transferases, and small protein-activating enzymes (supplemental Fig. S3B). For example, threonine-type endopeptidase and peptidase activities were significantly overrepresented with a p value of 2.17 ϫ 10 Ϫ9 . In addition, functions including cytoskeletal protein binding, calcium-dependent phospholipid binding, proteasome activator activity, and coenzyme binding were also significant (p Ͻ 0.01). In the GO cellular component category, we found that most of the BAG3 interactome proteins identified were localized in non-membrane-bounded organelles (p ϭ 8.83 ϫ 10 Ϫ4 ), the cytosol (p ϭ 6.23 ϫ 10 Ϫ15 ), and the cytoskeleton (p ϭ 2.99 ϫ 10 Ϫ4 ) (supplemental Fig. S3C). Notably, 21 BAG3 interacting proteins were attributed to the proteasome complex (p ϭ 7.25 ϫ 10 Ϫ20 ), suggesting that that BAG3 might interact with the proteasome and regulate its activity.
The subcellular localization and function of proteins are often dictated by their domains. We performed a similar analysis at the structural level by investigating the enrichment of protein domains in the BAG3 interactome. 93.4% of the BAG3 interacting proteins contain at least one Pfam-annotated domain; those domains that are significantly enriched in the BAG3 interactome are listed in Table II. Pfam domains associated with the proteasome and 14-3-3 proteins were most overrepresented in the BAG3 interactome (Table II). It is well established that Bag3 has a modular structure that contains many domains to mediate potential interactions with chaperones and/or other proteins (47). For example, it has been reported that the PXXP region of BAG3 can bind to Src homology 3 (SH3)-containing proteins (6). Consistent with these reports, nine BAG3 interactors are SH3-containing proteins, and SH3 domain families are enriched in the BAG3 interactome (p Ͻ 10 Ϫ6 ) (Table II). Therefore, it is likely that these proteins may use the SH3 domain to interact with the PXXP region of BAG3.
Signaling Pathway and Protein Interaction Network Analysis-To determine whether BAG3 plays a role in signaling pathways, we superimposed the human Kyoto Encyclopedia of Genes and Genomes pathway database on the BAG3 interactome. As shown in supplemental Table S8, many iden- tified BAG3 interacting proteins are involved in pathways responsible for the control of key physiological and pathological processes. In particular, 27 BAG3 partners were found to be involved in the regulation of protein degradation. Of these, 16 proteins are mapped to the proteasome, 8 to ubiquitin-mediated proteolysis, and 3 to the regulation of autophagy. Furthermore, we identified 12 BAG3 interacting proteins, including CHUK, HspA8, JUND, MAPK7, and MAPK11, which take part in MAPK signaling pathways involved in various cellular functions, including cell proliferation, differentiation, and migration (supplemental Table S8). In addition, five proteins, including Hsp90AA1, NFKBIB, and MAPK11, were assigned to the NOD-like receptor signaling pathway, a crucial pathway in the recognition of intracellular ligands. Furthermore, some of the proteins in our BAG3 interactome, such as Hsp90, MAPK11, CHUK, and 14-3-3, serve as key molecular switches in these signaling pathways. In summary, pathway analysis of the BAG3 interactome implied direct involvement of the BAG3 interactome in complex physiological and pathological systems such as the proteasome-ubiquitination system, immune responses, cancer, and neurodegenerative disorders.
To generate a BAG3 protein interaction network, we used the STRING database of physical and functional interactions (48). 283 (74%) of the identified BAG3 interacting proteins mapped to the BAG3 interaction network either directly or indirectly (Fig. 4A). Recent large-scale, high-throughput studies have established that most, if not all, cellular functions and processes are carried out by protein complexes rather than by individual proteins. In attempts to better understand biological processes, mapping interactions between protein complexes is as important as mapping interactions between individual proteins. To characterize protein complexes that associate with BAG3, therefore, we analyzed the BAG3 interactome for highly connected regions using MCODE (38) and determined the interconnectivity between these complexes and BAG3 by means of protein network analysis. As shown in Figs. 4B-4E, several complexes and cellular functions formed prominent, highly connected clusters. For example, many proteasome subunits are connected in a dense protein-protein interaction network that forms the densest cluster (Cluster 1), whereas BAG3 forms a complex with Hsp90AA1, HspA1A, HSF1, and ST13 (Cluster 4). It is therefore tempting to speculate that the physiological interaction between these two complexes might contribute to the cooperation and/or coordination of their functions in the control of numerous intracellular signaling and regulatory pathways.
Validation of BAG3 Interacting Proteins by Reciprocal Immunoprecipitation-To validate some of the proteins identified via QUICK and proteome microarrays, co-immunoprecipitation experiments and Western blot analysis were performed. Based on how well their biological functions and importance were known, we selected two proteins identified only by QUICK, 14-3-3 and ATG3, and two proteins identified only by proteome microarrays, MAPK7 and HSF1, for validation. We also included two proteins, PSMB1 and HSP70, that were identified by both methods for validation. 14-3-3 belongs to the 14-3-3 protein family, which mediates signal transduction by binding to phosphoserine-containing proteins (49). ATG3, the yeast Apg3 homolog, is a protein-conjugating enzyme essential for autophagy (50). PSMB1 is the beta subunit of the highly ordered ring-shaped 20S core of the proteasome (51). MAPK7 plays a role in diverse cellular processes such as proliferation, differentiation, and cell survival (52). HSF1 is a DNA-binding protein that specifically binds heat shock promoter elements and activates transcription (53). Hsp70 was included as a positive control for the co-immunoprecipitation experiments because it was previously reported to physically interact with BAG3 (5). As shown in Fig. 5A, endogenously expressed 14-3-3, ATG3, PSMB1, MAPK7, HSF1, and Hsp70 proteins could be readily detected in the BAG3 complex immunoprecipitated with anti-BAG3 antibodies. These interactions are specific, as no detectable signals were observed in IgG negative controls. Furthermore, reciprocal immunoprecipitation-Western blot analysis confirmed their specific interactions in cells (Figs. 5B-5G). Thus, compelling evidence shows that 14-3-3, ATG3, PSMB1, MAPK7, HSF1, and Hsp70 proteins are bona fide BAG3 interacting proteins.
Validation of BAG3 Interacting Proteins by SPR-As an independent approach, we employed SPR to determine the binding affinity of BAG3 interacting proteins. We included one protein, ENO1, which was identified as a background protein as a negative control for these SPR experiments. This in vitro assay system clearly showed that 14-3-3, ATG3, PSMB1, MAPK7, HSF1, and Hsp70 interact with BAG3 immobilized on the SPR surface. As expected, we did not observe any appreciable interaction between ENO1 and BAG3 under the same conditions (Table III). Using the SPR data, we further determined the affinity constants (K D ) for the interaction between immobilized BAG3 and 14-3-3, ATG3, PSMB1, MAPK7, HSF1, and Hsp70 as 51, 170, 1,200, 1,100, 34, and 17 nM, respectively (Table III).
BAG3 Regulates Proteasome Activity-GO analysis revealed that the most overrepresented biological process in the BAG3 interactome is involved in the control of the proteasome function, suggesting that BAG3 is a critical regulator of proteasome activity. To evaluate the role of BAG3 in regulating proteasome activity, we knocked down BAG3 using RNAi and examined its potential impact in HeLa cells. Western  blotting using anti-BAG3 antibodies in RNAi-treated cells confirmed a significant decrease in the protein level of BAG3 relative to the negative control (Fig. 6A). Because previous studies have shown that BAG3 is induced by the proteasome inhibitor MG-132 at the transcriptional level, and because the induction of BAG3 by proteasome inhibition is anti-apoptotic (54,55), we performed MG-132 treatment together with BAG3 knockdown. As shown in Fig. 6B, knockdown of BAG3 alone did not affect apoptosis in HeLa cells. However, in the BAG3knockdown background, MG-132-induced apoptosis was significantly higher (49.5 Ϯ 5.1%) than that in control cells treated with a control siRNA (36.1 Ϯ 2.4%; p Ͻ 0.05). We further measured the effects of BAG3 knockdown on 20S proteasome activity in HeLa cells. As shown in Fig. 6C, knock-down of BAG3 resulted in a significant decrease in 20S proteasome activity (p Ͻ 0.05). In cells transfected with BAG3 siRNA, proteasome activity of MG-132-treated cells (42 Ϯ 3.1%) was significantly lower than that in cells transfected with control siRNA (58 Ϯ 6.4%; p Ͻ 0.05). Inhibition of proteasome activity should prevent the rapid degradation of ubiquitin-tagged proteins. Our results (Fig. 6D) show that BAG3 knockdown caused an accumulation of ubiquitintagged proteins in HeLa cells, as determined via Western blotting using anti-ubiquitin antibody. The 20 and 26S proteasomes have at least three distinct kinds of proteolytic activities, including chymotrypsin-, trypsin-, and caspase-like activities (40). To monitor specific protease activities of the 20 and 26S proteasomes, three fluorogenic substrates (Suc- LLVY-AMC for chymotrypsin-like, Z-ARR-AMC for trypsinlike, and Z-LLE-AMC for caspase-like activity) were incubated with cell lysates from treated or untreated cultures. As shown in Fig. 6E, knockdown of BAG3 resulted in a significant decrease in chymotrypsin-like activity of the 26S proteasome (p Ͻ 0.05). The caspase-like and trypsin-like activities of the 26S proteasome were slightly, but not significantly, less in BAG3 silenced HeLa cells. In cells transfected with BAG3 siRNA, chymotrypsin-like activity of MG-132 treated cells (43 Ϯ 5.4%) was significantly lower than that in cells transfected with control siRNA (56 Ϯ 4.5%) (p Ͻ 0.05), whereas the caspase-like and trypsin-like activities of the 26S proteasome were not significantly different. Similar to its effects on the 26S proteasome, BAG3 knockdown resulted in a significant (p Ͻ 0.05) decrease in the chymotrypsin-like activity of the 20S proteasome in HeLa cells with or without MG-132 treatment (Fig. 6F). The effects of BAG3 on the specific protease activities of the 20 and 26S proteasomes were further determined using additional A549 (human lung adenocarcinoma epithelial cell) cancer cells. Similar to its effects on apoptosis and the proteasome in HeLa cells, BAG3 knockdown resulted a significant (p Ͻ 0.05) increase in MG-132 sensitivity and decrease in the chymotrypsin-like activity of the 26S and 20S proteasome in A549 cells (supplemental Fig. S4).
These results prompted us to investigate the effect of BAG3 overexpression on proteasome activity and apoptosis. BAG3 was transiently overexpressed in HeLa cells, as confirmed by Western blotting (Fig. 7A). As shown in Fig. 7B, overexpression of BAG3 reduced MG-132 sensitivity in HeLa cells relative to cells transfected with a control vector. Quantitative analysis revealed that the percentage of apoptotic cells when BAG3 was overexpressed together with MG-132 treatment (24.9% Ϯ 3.1%) was significantly lower than that in the vector control (34.2% Ϯ 3.4%; p Ͻ 0.05). When BAG3 was overexpressed, the 20S proteasome activity of MG-132 treated cells (77% Ϯ 7.1%) was significantly higher (p Ͻ 0.01) than that of MG-132-treated cells transfected with a control vector (62% Ϯ 5.4%) (Fig. 7C), and the accumulation of ubiquitintagged proteins in MG-132-treated HeLa cells was reduced (Fig. 7D). We further examined the specific effect of BAG3 on proteasome activity. As shown in Figs. 7E and 7F, the chymotrypsin-like activity of MG-132-treated cells (68% Ϯ 4.2% for 26S and 71% Ϯ 6.1% for 20S) was significantly higher (p Ͻ 0.05) than that of MG-132-treated cells transfected with a control vector (54% Ϯ 5.4% for 26S and 58% Ϯ 4.5% for 20S). In addition, similar effects of BAG3 overexpression on apoptosis and proteasome activities were also observed in A549 cells (supplemental Fig. S5).
Together, these results suggest that BAG3 plays an important role in regulating proteasome activity in cancer cells. Our results also support the notion that the inhibition of BAG3 could in principle be a useful therapeutic strategy for the treatment of cancer.

DISCUSSION
BAG3 has a multifaceted ability to form complexes with many proteins; it is involved in a number of major biological processes including apoptosis, development, cytoskeleton organization, virus replication, and autophagy. Recently, we found that BAG3 plays an important role in SARS-CoV repli- cation (13). However, the mechanisms underlying BAG3's functions remain somewhat enigmatic and controversial. A key step in understanding the physiological function of BAG3 is assigning BAG3 to appropriate biochemical pathways through the identification of its interacting partners. In this study, we performed a global analysis of the interactome of BAG3 using the QUICK and proteome microarray methods. Our results clearly demonstrate that datasets generated via QUICK and proteome microarrays form distinct, partially overlapping sets with different protein-protein interaction characteristics (Fig. 3). The differences between the QUICK and proteome microarray results were not unexpected, as many studies have shown that different proteomic methods possess distinct strengths and weaknesses and provide complementary types of information (56 -58). Therefore, we consolidated these two datasets into a single set that was designated as the BAG3 interactome. This dataset should provide a starting point for more in-depth studies on the functions of BAG3.
Because our bioinformatics analysis indicated that the proteasome-ubiquitination process and proteasome complex were strongly enriched in the BAG3 interactome, we investigated their relationship in greater depth. We report for the first time that BAG3 interacts with the proteasome and regulates its activity. This agrees well with reports that proteasomeinteracting partners dynamically associate with proteasomes and alter proteasome function (59 -64). In this study, downregulation of BAG3 was found to decrease proteasome activity and increase cell death in response to proteasome inhibition. This effect can be attributed to the fact that BAG3 inhibition suppresses a wide constellation of anti-apoptotic protective cellular pathways, thereby rendering cells more sensitive to various pro-apoptotic stimuli. This indicates that BAG3 might reduce the apoptotic response to proteasome inhibition, probably by maintaining proteasome activity. Proteasomes degrade a multitude of protein substrates in the cytosol and nucleus and are thus essential for many aspects of cellular function (65). Proteasome inhibitors represent a novel class of anticancer drugs with very promising clinical activity in multiple myeloma and possibly in various solid tumors (66). However, proteasome inhibitors also activate survival signals, and thus it seems plausible that the activation of survival signaling cascades might compromise their antitumoral efficacy (67). The synergistic effect between BAG3 inhibition and MG-132 observed here confirms the functional significance of BAG3 and suggests the possibility that combination treatments that include BAG3 inhibitors may augment clinical efficacy and overcome clinical refractoriness to proteasome inhibitors.
Heat shock proteins, known as molecular chaperones that assist in protein folding, are the products of several distinct gene families that are required for cell survival during stress and are named according to their molecular weights (68). Different classes of heat shock proteins play diverse roles in governing proper protein assembly, folding, and translocation (69). Among them, the Hsp70s control protein homeostasis in the cell during and following applied stress and, together with other chaperones, regulate physiological processes, including vesicular trafficking and signaling for growth, differentiation, and apoptosis (70). It is well established that BAG3 contains a conserved BAG domain able to interact with Hsp70 proteins and that the influence of BAG3 on cell survival is mediated at least in part by its binding to Hsp70 proteins (3,71). Our data therefore support the idea that BAG3 interacts with both the proteasome and Hsp70 and raise the possibility that BAG3 might form a complex in vivo with Hsp70 and the proteasome. It is likely that the interaction of BAG3 with the proteasome and Hsp70 may regulate the stability and dynamics of targeted proteins. It is also possible that such interactions may be dynamically regulated by stimuli in coordination with the diverse functions of these proteins during cellular processes such as cell survival (stress response), proliferation, migration, and apoptosis. The interaction of BAG3 with the proteasome and Hsp70 might therefore represent a novel mechanism for maintaining the stability and dynamics of targeted proteins in cancer cells. Further work that combines individual knockouts of BAG3, Hsp70, and proteasome subunits with quantitative proteomics is ongoing to clarify the role of these proteins in cancer cells.
Seventeen of the BAG3 interacting proteins belong to the transfer/carrier protein class, including five ubiquitin-conjugating enzymes and two autophagy-related proteins (Atg3 and Atg7). The interaction of Atg3 with BAG3 was confirmed via reciprocal co-immunoprecipitation and SPR analysis ( Fig.  5 and Table III). ATG genes are highly conserved among eukaryotes and are required for macroautophagy and its related processes. The discovery of ATG genes, originally in yeast and subsequently in multicellular organisms, has provided an important breakthrough in the understanding of autophagy and of its functions in physiology and disease (72). To date, more than 30 ATG genes have been identified. The corresponding gene products constitute the core machinery that coordinates specific steps in the autophagic pathway (73). It is well established that protein degradation through both the proteasome and autophagy pathways is regulated by BAG3 and its interacting proteins (3,9). BAG3 participates, along with HspB8, a member of the HspB family of molecular chaperones, in the degradation of misfolded and aggregated proteins via autophagy (74,75). BAG3 may act as a molecular switch to determine whether proteins are degraded in a proteasomal or autophagic manner (3,9). Additionally, BAG3 was found to act in concert with the ubiquitin-binding protein sequestosome 1 to increase autophagic activity (3,9). In this study, we have shown clearly that Atg3 and Atg7 are BAG3associated proteins. This suggests that, through binding to Atg proteins, BAG3 may regulate the formation and maturation of autophagosomes and their fusion with lysosomes. Further experiments are necessary to establish the functional consequences and detailed mechanisms of these interactions, as well as to define further the role of BAG3 in autophagy signaling. These data, combined with data from previous studies, will provide a new window into the autophagy machinery and a possible new target for therapeutic intervention.
To our knowledge, this is the first study that has combined QUICK and proteome microarrays to obtain a comprehensive profile of protein-protein interactions. Both of these technologies have their own unique advantages and limitations, and thus combining the two provides a powerful approach for proteomic investigations. With BAG3 as an example, we have shown that QUICK and proteome microarrays are complementary approaches (Fig. 3A). In addition, well-validated BAG3 interacting proteins such as Hsp70 (5), sequestosome 1 (9), and HspB6 (45) were successfully identified via both of these approaches. The lists of BAG3 interacting proteins generated by the two approaches overlapped to a high degree; 47% (81/174) of the QUICK candidates were also identified by the proteome microarray. Furthermore, we constructed a very comprehensive interactome with 382 BAG3 interacting proteins by combining the independent results obtained using these two methodologies. This complementary strategy of using both QUICK and proteome microarrays could easily be adopted for interactome construction for other proteins of interest, such as transcriptional factors, the protein-protein interactions of which cannot theoretically be studied using traditional yeast two-hybrid assays.
The complete BAG3 interactome obtained using these complementary approaches provides an essential step toward our understanding of the system-wide role of BAG3 in the regulation of various cellular processes such as signal transduction, transcription, proteasome-ubiquitination, and autophagy. It is now important to further characterize the interactions of BAG3 with different binding partners in different physiological and pathophysiological contexts.