Qiongyu
Li
a,
Yixuan
Xie
a,
Rachel
Rice
a,
Emanual
Maverakis
b and
Carlito B.
Lebrilla
*ac
aDepartment of Chemistry, University of California Davis, Davis, California, USA. E-mail: cblebrilla@ucdavis.edu
bDepartment of Dermatology, School of Medicine, University of California Davis, Davis, California, USA
cDepartment of Biochemistry, University of California Davis, Davis, California, USA
First published on 30th April 2022
Antibodies targeting specific antigens are widely utilized in biological research to investigate protein interactions or to quantify target antigens. Here, we introduce antigen–antibody proximity labeling (AAPL), a novel method to map the antigen interaction sites as well as interactors of antibody-targeted proteins. As a proof of concept, AAPL was demonstrated using sodium/potassium transporting ATPase (ATP1A1) and epidermal growth factor receptor 2 (ERBB2)-specific antibodies that were modified with an Fe(III) catalytic probe. Once bound to their target proteins, Fe(III)-induced catalytic oxidation occurred in proximity of the antigen's epitope. Oxidative proteomic analysis was then used to determine the degree of oxidation, the site of oxidation within the targeted antigen, and the interacting proteins that were in close proximity to the targeted antigen. An AAPL score was generated for each protein yielding the specificity of the oxidation and proximity of the interacting protein to the target antigen. As a final demonstration of its utility, the AAPL approach was applied to map the interactors of liver–intestine-cadherin (CDH17) in colon cancer cells.
Mass spectrometry-based proximity labeling coupled with proteomic techniques is another widely used approach for probing protein–protein interactions. By fusing the targeted proteins with specific enzymes or catalytic probes, proteins localized near the interaction site are labeled by catalytic reactions that can then be characterized by quantitative proteomics.15–17 With proximity labeling proteomics, interactors with both strong or weak interactions of the target are identified. Several studies have employed this approach to characterize protein interactions. Lobingier et al.18 have developed a proximity labeling method combining with the engineered ascorbic acid peroxidase (APEX) to investigate the interaction network of the protein β adrenergic receptor (B2AR). In another study, a biotinylation reaction was used for the proximity labeling of cell membrane proteins.19 However, biotinylation-based approaches require more complicated sample preparation and streptavidin enrichment does not provide information about the site of the interaction. Hydroxyl radical protein footprinting (HRPF) is commonly used for the determination of protein structure and protein–ligand interactions, where the interaction sites are characterized by the site-specific oxidation of proteins and mass spectrometry.20,21 Combining proximity labeling and hydroxyl radical protein footprinting, Li and coworkers developed a labeling approach by immobilizing a catalytic Fe(III) probe on cell membrane sialic acids to map the proteins interacting with sialylated glycans through site-specific oxidation.22 The enrichment of interacting proteins was not required with this approach, which assess the extent of the protein–protein interactions by quantifying the degree of oxidation at interaction sites.
In the current study, we developed antigen–antibody proximity labeling (AAPL) as a technique to map the epitopes of a targeted antigen as well as the interactors of the antigen by proximity-dependent oxidation. Two separate approaches were employed to modify the antigen-specific targeting antibody to create a molecule with catalytic function. In one method, the N-glycans of the antibodies were modified by changing terminal galactoses (Gal) to N-azidoacetylgalactosamines (GalNAz), followed by the conjugation of a synthesized clickable probe, dibenzocyclooctyne-functionalized bromoacetamidobenzyl-EDTA iron(III) chelate (DBCO-FeBABE) to the GalNAz groups. In another approach, the antibody was directly modified with the iron(III) chelate of bromoacetamidobenzyl-EDTA (FeBABE) molecule at cysteine residues via their free thiol group. Modified antibodies directed against sodium/potassium transporting ATPase (Na+/K+-ATPase) and human epidermal growth factor receptor 2 (ERBB2) were used as prototypical examples, with the goal of labeling the targeted cell membrane proteins and their interactors in the immortalized prostate epithelium cell line, PNT2, and the human breast cancer cell line, SKBR3, respectively. The oxidations of targeted antigens and their associated proteins were monitored through the quantitative oxidative proteomic analysis to optimize the reaction conditions for specific labeling. The AAPL method with optimized condition was then applied to study the interactors of liver–intestine cadherin (LI-cadherin/CDH17), which is a cadherin-like protein expressed mainly in liver and intestine, as well as the involvement of CDH17 glycosylation playing in its interactome.
Both versions of probes (FeDBAb and FeAb), were tested in parallel using cell lines expressing their respective antigens. After antibody binding, hydrogen peroxide (H2O2) was added to the culture medium to produce hydroxyl radicals in the proximity of the attached antibody (Fig. 1b). Optimal conditions to minimize overoxidation have been previously reported.22 Here, controls to assess background oxidation included cell line only (C1), cell line treated with unmodified antibody (C2), and cell line treated with H2O2 (C3). For each of the engineered antigen-specific oxidizing probes, treatment conditions were optimized to obtain the optimal oxidation conditions. For FeAb, various ratios of antibody to FeBABE were tested to optimize Fe(III) conjugation (M1–M3, ESI Table 1†). The resulting oxidizing probes were also assessed under various oxidation conditions to optimize target antigen oxidation. Similarly, FeDBAb optimization also involved varying concentrations of FeDBAb to optimize the oxidation conditions (M4–M6, ESI Table 1†).
Membrane proteins were enriched after H2O2 treatment and subjected to quantitative oxidative proteomics analysis (Fig. 1c). The degree of oxidation for each site on a protein was determined by calculating the ratio of the oxidized peptide intensity to the summed unmodified and modified peptide intensities for each tryptic peptide. Additionally, the extent of protein oxidation (EPO) for each protein was obtained by calculating the ratio of the summed oxidized peptide intensities to the summed total peptide intensities of each protein. The optimal oxidation conditions were determined based on the EPO. The oxidation sites on ATP1A1 and ERBB2 and their positions in the protein domains showed the specificity of each antibody (Fig. 2a and b). It was observed that among conditions M1–M6 and controls, conditions under M5 yielded the highest EPO for both ATP1A1 and ERBB2 (Fig. 2c and d) with their respective antibody.
Antibody cross reactivity was assessed by examining the oxidation of ATP1A1 after treatment with ERBB2-specific oxidizing antibody probe. We noted that three oxidation sites were observed on ATP1A1 after incubation with the ERBB2-specific oxidizing antibody (Fig. 2a); however, none were highly oxidized under M5 conditions (Fig. 2e), demonstrating the high specificity of the AAPL method. To illustrate the EPO of all oxidized proteins, a heatmap was generated for proteins under the M1–M6 and three control conditions (Fig. S3†) using both ATP1A1 and ERBB2-specific oxidizing antibodies. The heatmap demonstrated that similar conditions clustered together, and condition M1 was optimal for FeAb and M5 for FeDBAb.
A method for calculating topological scores (TopS) of protein interactions in affinity purification mass spectrometry (AP-MS) was adapted here to generate an AAPL score.23 In TopS, the spectral counts of each bait-enriched prey protein are used to calculate a TopS score to evaluate the extent of interaction of the prey protein and the bait protein. Similarly, in AAPL the extent of interaction was evaluated based on an AAPL score calculated using the spectral count of each oxidized protein together with its EPO. Detailed equations for the calculations are included in the Method section of the ESI.† The proteins with AAPL scores larger than 0 were considered as more specifically oxidized relative to controls. This method was applied to identify ATP1A1 and ERBB2 interactors, which generated over 40 proteins with AAPL scores pass the 0 threshold under M5 condition.
To generate a more stringent selection criteria for the antigen interactors, we first calculated the average number of oxidation sites for all proteins with AAPL values above 0 and below 1 after ATP1A1- and ERBB2-targeted oxidation. Results differed depending on the targeted antigen. For ATP1A1, a score of 0 corresponded to 2.8 sites, while a score of 1 corresponded to 4.8 oxidation sites. In ERBB2-targeted proximity labeling proteins with AAPL score above 0 had on average 2 oxidation sites, while those with AAPL score above 1 had on average 2.4 oxidation sites. The number of average oxidation sites suggest that proteins with AAPL scores above 1 are more oxidized, indicating a stronger interaction with the target protein.
To validate the AAPL scores, the classification of interactors was compared to STRING, which provides scores based on the extent of interaction between proteins using derived literature data. The oxidized proteins were examined with STRING. Oxidized proteins that were known to interact directly with the target protein were classified as Type 1 (T1) and Type 2 (T2) with T1 defined as high confidence interactors, having STRING scores above 0.7, and T2 as medium level interactors, having STRING scores in the range of 0.4–0.7. Secondary interactors, those interacting with T1 and T2 proteins and not directly with ATP1A1 or ERBB2, were classified as T3 interactors. The remaining oxidized proteins were further classified as T4 interactors. The numbers of each type of interactors for ATP1A1 and ERBB2 based on the STRING analysis are summarized in Fig. 3a and b, respectively.
The distributions of AAPL scores for T1 and T2 interactors were then examined. For ATP1A1, more than 50% of the interactors (T1 and T2 combined) had AAPL scores larger than 1. Nearly 15% had AAPL scores in the range of 0–1, and the remainder were below 0 (Fig. 3c). For ERBB2 around 50% of T1 and T2 interactors had AAPL scores >1, with 8% in the range of 0–1, and the remainder below 0 (Fig. 3d). Thus, for the T1 and T2 interactors with AAPL scores larger than 0, a large fraction was with AAPL scores above 1. In addition to AAPL scores, the extent of protein oxidation fold changes (EPO FC) for T1 and T2 interactors was also determined. It was found that for ATP1A1, all the T1 and T2 interactors (corresponding to 9 proteins) with AAPL scores larger than 0 had total EPO FC from 1.5 to 5. For ERBB2, all the T1 and T2 interactors (6 proteins) with AAPL scores larger than 0 had total EPO FC from 1.5 to 3, with one exception. Thus, the STRING and AAPL methods identified similar sets of proteins as strong interactors with the target antigen.
We therefore classified the oxidized proteins based on the AAPL scores with the following criteria. Proteins with AAPL scores above 1 were classified as Level 1 (L1) interactors, while those with AAPL scores in the range of 0–1 and EPO FC above 1.5 were classified as Level 2 (L2) interactors. With this classification criteria, several of STRING classified T3 and T4 oxidized proteins were also found as L1 interactors of ATP1A1 and ERBB2 in PNT2 and SKBR3 cell lines, respectively; while not all the T1 and T2 proteins were primary interactors. We should note interactors identified by the AAPL scores may be dependent on the specific cell lines (or tissue) and culture conditions used in the analysis. In total, four T3 proteins and thirteen T4 proteins were found as AAPL L1 interactors of ATP1A1. Twelve T3 proteins and three T4 proteins were found as AAPL L1 interactors of ERBB2. ATP1A1 and ERBB2 L1 and L2 interactors and corresponding STRING-type associations are summarized in Fig. 3e and f.
To evaluate how well the AAPL identified interactors of ATP1A1 and ERBB2 were related to the corresponding antigen, gene ontology enrichment analysis was conducted. For ATP1A1 interactors, the analysis of molecular function demonstrated that several genes belong to the groups of proteins that were involved in active transmembrane transporter activity and wide pore channel activity, which were related to the function of ATP1A1.24,25 All of enriched molecular function terms had adjusted p-value below 0.01, indicating a significant level of enrichment (Fig. 4c).
Interactors of ERBB2 were involved in four molecular functions, including lyase activity, oligosaccharyl transferase activity, apolipoprotein binding, and macromolecule transmembrane transporter activity (Fig. 4d). Among the STRING T1 interactors, G3P (glyceraldehyde-3-phosphate dehydrogenase) and FAS (fatty acid synthase) were identified as AAPL L1 interactors. Four STRING T2 interactors were also identified as AAPL L1 interactors (S100A14, TFRC, PTPN1, and CPNE3). Among these, CPNE3 (copine-3) was demonstrated to interact with ERBB2 through a receptor dependent manner and it may be a novel factor regulating the ERBB2-dependent cancer cell motility.26
The specificity of the AAPL method was also investigated by comparing the assigned interactor types of proteins that were found commonly in both ATP1A1 and ERBB2 AAPL networks. A merged interaction network was also generated and is illustrated in Fig. S4.† In total, 5 proteins including GAPDH, TFRC, VDAC1, ATP2A2, and ATP1A1 were found in both interaction networks. It was found that AAPL interactors of ERBB2, including GAPDH and TFRC were not AAPL interactors of ATP1A1, corresponding to their types, they were not increasingly oxidized under ATP1A1 AAPL condition. Similarly, ATP2A2 was a member of the ATP1A1 AAPL interaction network but not the ERBB2 AAPL interaction network. Since ATP2A2 was increasingly oxidized in response to ATP1A1 AAPL while not ERBB2 AAPL, we hypothesize that the degree of protein oxidation is proportional to the degree of its interaction with the target protein. Lastly, VDAC1 was identified as AAPL interactor of both ATP1A1 and ERBB2, corresponding to the observation that the protein was increasingly oxidized under the AAPL condition of both ATP1A1 and ERBB2. These observations demonstrate the specificity of the AAPL method to the targeted antigen and its interactors.
The employment of AAPL to differentiated Caco-2 cells generated 30 AAPL interactors of LI-cadherin, with 25 L1 interactors and 5 L2 interactors. STRING was used again to organize the oxidized proteins based on known interactions. T1 interactors predicted by STRING, such as catenin subunits and desmoplakin, were found in the cell but were not labeled as L1 or L2 interactors, which is consistent with the independence of the LI-cadherin mediated cell aggregation to the commonly known cadherin-family interacting proteins.30 Among the STRING T2 proteins, PTPN1 and VTNC were identified as AAPL interactors, and they have been shown to interact with E-cadherin interaction network.31,32 Besides, 12 STRING T3 proteins and 15 STRING T4 proteins were determined as AAPL interactors of LI-cadherin in the differentiated Caco2 cell line. The complete interaction network is shown in Fig. 5a. The gene ontology enrichment result is shown in Fig. 5b. It demonstrates that AAPL interactors of LI-cadherin are involved in several molecular functions including integrin binding. These results demonstrate that AAPL successfully labeled proteins that are known to be related to cadherin. The results for the oxidized peptides and proteins of LI-cadherin in differentiated Caco-2 cells are summarized in ESI_Table_III.†
Cell surface glycosylation was modified with various inhibitors, including kifunensine (Kif), an α-mannosidase inhibitor,33 and 3Fax-feracetylNeu5Ac, a sialyltransferase inhibitor (SI).34 These inhibitors are known to dramatically alter cell glycosylation. For example, Kif treatment of Caco-2 cells drastically increases high-mannose N-glycan surface expression. In contrast, SI extensively decreased the relative abundances of sialylated glycans. To investigate the effects of glycans on the protein interaction networks, differentiated Caco-2 cells were treated with the different glycosylation inhibitors, followed by incubation with the LI-cadherin-specific oxidizing probe and proximity labeling with H2O2. The oxidized sequences were then quantified by nanoLC-MS/MS.
In total, there were 30 oxidized proteins that were also glycosylated based on the glycoproteomics analysis. Among these proteins, six of them were identified as interactors of LI-cadherin, with three classified as L1 interactors including CEAM1, ATLA3, and STT3A, and three as L2 interactors including 1A02, VTNC, and NOMO2. To examine whether the extents of the interactions were altered by the inhibitors, we measured the extent of protein oxidation of the native and the glycan modified cells and used the variations to determine extent of interaction. It was found that the oxidation of CEAM1, ATLA3, 1AO2, and VTNC were all decreased when cells were pretreated with kifunensine, suggesting that these proteins interacted less with LI-cadherin. However, there were little to no decrease in the extent of oxidation for STT3A and NOMO2. An identical behavior was observed when the sialic acid expression was decreased by SI (Fig. 6d). Importantly, based on glycoproteomic analysis the former four proteins are all sialylated, while the latter two proteins have primarily high mannose type glycans (Fig. 6e). These results strongly suggested that the interactions of LI-cadherin to the group of four proteins occurred mainly through sialic acid interactions while the two were mainly through high mannose interactions.
In addition to the specific oxidation of the target proteins, proteins interacting with the target protein are also oxidized, albeit to a lesser extent. These proteins were primarily those that associated with the target antigen as verified by comparing them to published interactors. The STRING algorithm collates all published interactions and classifies them based on the extent of the interactions. The AAPL approach obtained similar results with proteins identified that matched those obtained by STRING. There were differences between those that were oxidized and those that were previously identified as interactors. The oxidized secondary proteins overlapped with the STRING interactors approximately 50% of the time, when they were present. The major source of the discrepancy could be that protein–protein interactions on the cell membrane may depend on several factors including the specific cell line and culture conditions used. Other proteins were also more highly oxidized than the known STRING interactors, suggesting stronger interactions in the systems studied here than those previously catalogued by STRING. Other factors may also be involved in the oxidation of the secondary proteins including non-specific interactions between the antibody and these proteins. It is worth noting that the interactors characterized through AAPL were in close proximity to the target protein. However, additional methods would be required to further validate the possible direct physical interaction. Nonetheless, AAPL would be well suited for discovering new interactions that can be further validated.
The variations in the oxidation of the AAPL interactors in the cells observed with LI-cadherin presents a new opportunity to elucidate the role of glycosylation in protein–protein interactions. Altering the glycan expression of some AAPL interactors from complex type glycans to high mannose or from sialylated type glycans to non-sialylated glycans affect the extent of oxidation of AAPL interactors, demonstrating the suppressed interactions of proteins to LI-cadherin. Interestingly, LI-cadherin and its AAPL interactors with their extent of interaction affected by altered glycosylation are post-translationally modified by complex type highly sialylated glycans. In contrast, AAPL assessment of cell surface proteins typically modified by high mannose glycans, which contain no sialic acid, yielded similar EPO results before and after Kif or SI treatment. Hence, the AAPL method can be applied to the study of glycosylation mediated protein–protein interactions.
Footnote |
† Electronic supplementary information (ESI) available. See https://doi.org/10.1039/d1sc06898a |
This journal is © The Royal Society of Chemistry 2022 |