Systems analysis of EGF receptor signaling dynamics with microwestern arrays.

We describe microwestern arrays, which enable quantitative, sensitive and high-throughput assessment of protein abundance and modifications after electrophoretic separation of microarrayed cell lysates. This method allowed us to measure 91 phosphosites on 67 proteins at six time points after stimulation with five epidermal growth factor (EGF) concentrations in A431 human carcinoma cells. We inferred the connectivities among 15 phosphorylation sites in 10 receptor tyrosine kinases (RTKs) and two sites from Src kinase using Bayesian network modeling and two mutual information-based methods; the three inference methods yielded substantial agreement on the network topology. These results imply multiple distinct RTK coactivation mechanisms and support the notion that small amounts of experimental data collected from phenotypically diverse network states may enable network inference.

Systems-level understanding of protein functions in biological processes remains a challenge. The western1 blot is a powerful protein analysis method because the electrophoretic separation step allows for reduction in sample complexity and the antibody detection step then results in signal amplitude proportional to the abundance of the immobilized antigen at a physical location on the detection membrane that can be related to molecular size standards. Because western blots require a relatively large amount of sample and a great deal of human labor, they have been of limited utility in large-scale protein studies. Reverse phase lysate arrays (RPAs), performed by arraying lysates directly on nitrocellulose-coated slides and probing them with antibodies, are useful for quantifying large numbers of proteins from limited amounts of material such as in biomarker discovery2,3. In contrast to western blots, however, RPAs lack confirmatory data for signal veracity; in a side-by-side comparison of measurements from RPAs and western blots, only 4 of 34 phospho-specific antibodies examined generated equivalent information with the two approaches4. Antibody cross-reactivity was concluded to contribute substantial noise to RPAs, confounding true protein measurements. Many antibodies have been validated for use with the Luminex xMAP bead-sorting system, although about 1,000-fold greater cell material is required per protein analysis than RPAs and the cost of detection reagents per protein about 30-fold higher. Flow cytometry permits a (relatively small) cohort of proteins to be examined simultaneously in individual cells; this multiplexing feature has been exploited with Bayesian network modeling to predict new signaling network causalities5.

In contrast to antibody-based methods, mass spectrometry (MS) can identify novel proteins. Using MS, thousands of peptides were assessed in lung cancers to identify commonly activated receptor tyrosine kinases and downstream signaling pathways6. Relative abundances can be examined quantitatively using isotopic labels across time points, cell types, or perturbations as in examination of phosphorylation dynamics of Hela7 and mammary epithelial cells8 following EGF or heregulin treatment. However, the large sample amount required by MS can limit the number of conditions that can be analyzed; ~108 cells are typically required for an MS experiment5 versus ~105 cells for an immunoblot or ~103 cells for RPAs9.

Here we describe micro-western arrays (MWA), which combine the scalability of RPAs while retaining vital attributes of western blots for highly multiplexed proteomic measurements: reduction of sample complexity and signals that can be related to protein size standards. In combination with suitable pan- and modification-specific antibodies, dynamics of protein abundance and modification may be simultaneously monitored across many samples. We demonstrate that MWA in combination with computational modeling techniques can yield useful systems-level biological insights for EGF receptor signaling dynamics.

RESULTS

Fabrication of micro-western arrays

Our strategy allows us to compare protein abundances and differences in post-translational modifications for cells stimulated under different conditions (Fig. 1). To interface the microscopic western blots with micro-titer-based liquid handling methods, cell lysates are printed via a non-contact microarrayer on gels in a series of 96 identical blocks with dimensions of a 96-well plate10 (Fig. 1). Using these dimensions 6 different lysates may be examined with 96 different antibodies or 24 different lysates may be examined with 24 different antibodies. An acetate running buffer is used to increase the migration rate of large proteins and slow the rate of smaller ones, obviating the need for a stacking gel. For each spot, 6 nl of sample was arrayed over the same gel position ten times, allowing for greater spotting density and signal than micro-depositing the entire 60 nl in a single dispense. One spot of size standard and six spots of experimental sample were arrayed at one millimeter pitch at the top edge of each block. After printing, the samples are subject to semidry electrophoresis for 12 min and then transferred to a nitrocellulose membrane. The membrane is placed in a 96-well gasket to isolate each set of 6 separated lysates; each block is then incubated with a different antibody. Following incubation with dye-labeled secondary antibody, the blot is scanned with an infrared fluorescence scanner. This format allows interrogation of 192 antibodies against 6 samples when two antibodies from different hosts (e.g. rabbit and mouse) are utilized. A total of 1,152 antibody-sample readouts is therefore possible per MWA device. Each measurement requires ~1,000 cells (equivalent to 250 ng of protein) and 16 ng of detection antibody and allows for the analysis of ~4,000 protein abundances.

Fig. 1

Micro-western array (MWA) procedure

Procedural schematic for the micro-western array method.

Validation of micro-western array method

We compared the resolution and linearity in signal of MWAs with macroscopic gels using a labeled molecular weight standard (LI-COR 928-40000) (Fig. 2a). For proteins of 150, 50, and 25 kDa, the intensity of each ladder spot was proportional to the fold dilution over a range of two orders of magnitude for both methods (Fig. 2b,c). The coefficient-of-variation from arraying, rehydration, and transfer of a single band of the LI-COR ladder across the area of the membrane was < 9%.

Fig. 2

MWA Validation of Linear Response

(a) (left) 5 μL of a series of 1/2 dilutions of LI-COR protein ladder were electrophoresed using traditional 10% SDS-PAGE. (right) 60 nl of the identical samples were run in the micro-western format. The difference in scale is noted below the figure. (b) The median net signal intensity was quantified for three bands (150 kDa, 50 kDa, and 25 kDa) of the LI-COR protein ladder in the traditional western blot as shown in part (a). The intensity vs. concentration of proteins displayed a linear relationship for all protein bands quantified. (c) The median net signal intensities quantified from micro-western dilution series. The intensity vs. concentration displayed a similar linear relationship as that from the traditional western. (d) 1/2 dilutions of lysates from A431 cells stimulated for 5 min with 200 ng ml-1 of EGF probed with seven rabbit primary antibodies directed at proteins spanning a range of molecular weights (175 kDa to 15 kDa) and detected with IR700-labeled secondary antibody. (e) The median net signal intensity for each band vs. relative concentration as shown in part (d). The graphs show a linear relationship between net fluorescence and concentration for all antibodies tested.

We then tested the linearity of signal response in quantifying proteins from A431 human carcinoma cell lysates using a two-stage fluorescent immuno-detection system (Fig. 2d,e). We used five phospho- and two pan-specific antibodies to analyze proteins from 15-175 kDa in EGF-stimulated A431 cells lysates. All MWAs showed a linear relationship between relative antigen concentration and signal intensity over their detectable range (from 100- to 1,000-fold). Assuming an expression level of 1.2 × 106 receptors per A431 cell11 EGFR was detectable down to one cell equivalent (~2 attomoles; ~340 femtograms). We assumed linearity for all further analysis.

Comparison of macro-western blots and micro-western arrays

To compare performance of MWAs with macro-western blots for monitoring phosphorylation dynamics, we selected a representative test-set of 11 antibodies. Four were previously shown to generate equivalent quantitative data by RPAs and western blots4; another four were shown to result in substantial compression of dynamic range by RPAs due to antibody cross-reactivity4.

Measurements obtained by MWAs were similar to macro-western blots for all antibodies (Fig. 3) and did not display the dynamic range compression observed for RPAs. For many protein phosphosites, including EGFR, IRS1, and AKT, we observed bands at the predicted size as well as at additional sizes that could obscure quantitative measurements by RPAs. The precision in estimating sizes of proteins > 100 kDa by MWAs was ~10 kDa, and for smaller proteins ~5 kDa. While protein sizes could be determined with precision approaching that of a standard western blot, proteins were not completely resolved unless they differed by more than the following: 75 kDa for proteins greater than 200 kDa; 50 kDa for proteins between 100 and 200 kDa; 25kDa for proteins between 50 and 100 kDa; and 10 kDa for proteins less than 50 kDa (Fig. 2d, Fig. 3), corresponding to a migration distance of about 1.5 mm, double the diameter of spotted protein (Fig. 2d, Fig. 3). Resolution equal to a macro-western blot could be obtained by electrophoresing the samples for twice the distance (Fig. 2a).

Fig. 3

Comparison of MWA to traditional western blot

A comparison of eleven traditional western blots is shown parallel to the identical samples run in triplicate in micro-western array format. Lysates from A431 cells stimulated with EGF (200 ng ml-1) and lysed at 0,1,5,15,30, and 60 minutes after stimulation along with the LI-COR protein ladder. β-actin monoclonal mouse primary antibody (detected with IR800 secondary antibody shown in green) was probed with each of the eleven rabbit primary antibodies (detected with IR700 secondary antibody shown in red) polyclonal antibodies to demonstrate equal loading of each sample. An arrow indicating the band quantified is indicated to the left of the blot along with the corresponding sizes of the LI-COR protein standard. Quantification of the fluorescence of the traditional western (red) is shown in comparison to the micro-western (blue). Error bars shown represent the standard error of the three technical replicates of the micro-westerns shown. Absolute sizes of the blots are indicated below to demonstrate the extent of the miniaturization of scale of the micro-western in comparison to the traditional western Blot.

Application of MWAs to analysis of EGFR signaling network

To examine EGF signaling dynamics using MWAs, we chose antibodies directed at a wide range of phosphosites to monitor many molecular biological processes (Supplementary Figure 1, Supplementary Table 1): 1) early positive growth factor response regulators; 2) negative signaling regulators; 3) downstream proliferation indicators; 4) nutrition sensing indicators; 5) adhesion and migration indicators; 6) phospho-lipid and calcium state indicators; 7) stress indicators; and 8) transcription and cell cycle indicators.

To observe signaling dynamics at doses approximating physiological levels12, we stimulated cells with 2, 50, 100, and 200 ng ml-1 EGF. A mock stimulation was performed to distinguish EGF-mediated signaling events from nutrition-related events. Figure 4 shows a MWA chip along with two magnified blocks demonstrating distinct temporal phosphorylation patterns of p-S6 Ribosomal Protein (Ser240,Ser244) and of p-EGFR (Tyr1068) (Supplementary Table 2). A mouse beta-actin antibody was used to probe all plate wells and demonstrates the reproducibility of the approach (Fig. 4, Supplementary Table 2). The coefficient-of-variation from arraying, rehydration, transfer, binding of primary antibody and secondary antibody was < 17%.

Fig. 4

An MWA containing 6 cell lysates probed with 192 antibodies

The red channel (700 nm laser) shows the stimulation of A431 cells with 200 ng/ml EGF probed with a panel of rabbit anti-human polyclonal antibodies detected with IR700-labeled secondary antibodies. The green channel (800 nm laser) reflects a scan of the samples probed with mouse monoclonal antihuman β-actin antibody detected with IR800-labeled secondary antibodies to demonstrate the consistency of printing across the area of the membrane. For antibody layout, see Supplementary Table 2.

We quantified 91 phosphosites from 67 proteins and 18 pan-specific protein abundances for a total of 75 proteins in technical triplicate resulting in ~9,800 signaling observations with sufficient residual lysates for further analyses using different antibody panels. Intensity and inferred sizes from spots detected with each antibody were recorded along with signal-to-background ratio (Supplementary Table 1, Supplementary Figures 2, 3). 17 of 91 phosphosites quantified here were previously quantified in one recent MS report using pan-phospho enrichment7 and 22 in another13 employing phosphotyrosine-specific enrichment (Supplementary Table 3). Many ubiquitous EGFR signaling proteins quantified by MWAs, including p-EGFR(Tyr845), p-SHP2(Tyr542), p-p70S6K(Ser371), p-Raf(c-)(Ser338), p90RSK(Ser380), and p-Stat3(Ser727), were not quantified in either MS study, suggesting that MS detects only a fraction of phosphorylation events elicited by EGF. 4 of 91 phosphosites quantified here were quantified by RPA and demonstrated to generate equivalent information as a western blot3.

Comparison of signaling network at different EGF input levels

We next asked whether biological insights could be revealed using the MWA method. Five clusters of signaling profiles were organized based on time of maximal phosphorylation (Fig. 5). Phosphosites within clusters were rank ordered by fold-change. At the 2 ng ml-1 input level, we observed several phosphosites from EGFR, ErbB2, PLC gamma, Gab1, Mek, p90RSK, p70S6K, and Crkl that were absent in mock treatment (Fig. 5, Supplementary Figure 4, 5).

Fig. 5

A clustered heatmap profile of fold changes for antibody bands representing specific phosphorylation sites of proteins in A431 cells over six time points for four stimulation conditions and one 0 ng ml-1 control

The net fold-change is color coded as indicated in the legend. Antibody bands are clustered into six clusters according to the time point (0, 1, 5, 15, 30, or 60 min) at which maximal fold change occurs. The antibodies are in descending order sorted in each cluster by the value of the fold change at the 200 ng ml-1 stimulation condition at the time point representative of that particular cluster. The antibody names are listed as given by the manufacturer followed in parentheses by an approximation of the size of each band. The EGF mock stimulation is shown on the right image block.

Conversely, many phosphosites related to phosphoinositide signaling displayed substantial fold-change in mock stimulation but not EGF treatment, including sites from PDK1 and its downstream targets AKT, PKC gamma, and PKC delta; downstream targets of AKT including mTOR, FOXO1; and mTor substrate p70S6K and its downstream target S6 Ribosomal Protein. We speculate that activation of PLC gamma following EGF stimulation leads to the hydrolysis of phosphatidylinositol 4,5-bisphosphate, causing down-regulation of PDK1 and AKT. Reduced AKT activity could produce the observed A431 cell cycle inhibition14 through decreased phosphorylation of cyclin dependent kinase (CDK) inhibitors, including p21CIP1(Thr145) and p27Kip1(Thr157). Consistent with this notion, Insulin-like Growth Factor (IGF), a potent PI3K stimulator, is also a potent mitogen for A431 cells12.

We then asked how the dynamic range and timing of the EGF signaling network were influenced by EGF input level. The first wave of phosphorylation peaking at the 1 minute point included 33 tyrosines from EGFR and other RTKs and membrane-localized proteins (Fig 5, Supplementary Figures 3-5, Supplementary Table 1). At 5 minutes, we observed serine and threonine sites from downstream kinases and transcription factors including Raf, MEK, p70S6 kinase, mTor, and ATF2. At 15 minutes, we observed phosphosites from Erk, P38 MAPK, and cell cycle-related kinases and substrates. Sites peaking at 30 minutes included those of the Crkl adapter protein and MAPKAPK2, a substrate of P38 MAPK. Proteins with sites peaking at 60 minutes included the PDK1 substrates AKT and PKCD, and the AKT substrate 4EBP1, among others. The timing of most phosphorylation events was not affected by EGF concentrations.

Bayesian network modeling of receptor layer connectivity

In order to try to elucidate the directional influences among phosphosites, we applied Bayesian network modeling approaches to the receptor layer of the EGF signaling network. This permitted us to verify known influences while identifying new directional relationships underlying receptor-level crosstalk. Bayesian networks are graphical representations of conditional independencies in a probability distribution over a set of variables15, and can potentially be inferred from experimental data such as that generated by MWAs. The network analyzed here comprised seventeen phosphosites: two from the Src kinase and fifteen from the ten RTKs for which fold-change measurements were observed specifically in all four EGF treatment input levels and for which antibodies were predicted to result in no cross-reactivity with each other (Fig. 6, Supplementary Table 4). Each time-point was considered as an independent sample of the EGF-stimulated network state, giving 20 samples for each phosphosite (4 conditions across 5 non-zero time points of one biological replicate). All data were normalized to the zero time point.

Fig. 6

Consensus model of EGFR receptor level influences modeled by Bayesian network inference with comparison to ARACNe and CLR. (a) A consensus model of the EGF signaling network obtained by exact Bayesian model averaging following Bayesian network inference (Supplementary Note 1). Significant (p < 0.001) positive edges (green), significant (p < 0.05) negative edges (red blunt edges), and interactions with a non-significant correlation coefficient (black) are shown. Edges for which the directionality could not be determined using equivalence class analysis are shown as undirected. (b) Heatmaps show the undirected adjacency matrices comparing the Bayesian network to the ARACNe and CLR networks. An edge between node i and node j is represented by matrix value (i, j). Because the undirected networks are compared, the adjacency matrix is symmetric across the diagonal, and thus only the lower triangular matrix of the adjacency matrix is shown. Edge weight thresholds were set to > 0.3 for the Bayesian network and ARACNe (using ARACNe Data Processing Inequality parameter τ= 0.03) and to Z > 1.13 for CLR. Eight of 11 edges present only in the Bayesian network and not in the ARACNe network would induce three-node triplets in the ARACNe network, which is precisely what ARACNe is designed to prune out (see Supplementary Figure 8). (c) Venn diagram comparing edges across the three networks. The ARACNe network forms a complete subnetwork of the CLR network and a near complete subnetwork of the Bayesian network, which forms a near complete subnetwork of the CLR network.

Given typically limited amounts of data, a variety of graph structures can be generated by Bayesian inference modeling that describe the data reasonably well, so a consensus model is often sought rather than aiming to find a unique best-scoring graph15. Accordingly, we illustrate here (Fig. 6) a consensus model containing only edges with a score > 0.3, derived from exact Bayesian network model averaging over all directed acyclic graph (DAG) structures having at most three parents per node16,17. By considering only those DAG structures in the equivalence class of the consensus model with a directed edge from p-SRC(Tyr416) to p-EGFR(Tyr845), directionality of the remaining compelled edges was determined18 (Supplementary Note 1). Signs of directional influences (i.e. positive vs. negative) could also be discerned. EGFR(Tyr845) is a known Src kinase substrate which is not phosphorylated by the EGFR kinase19. This prior knowledge was only used to distinguish edge directionality in the equivalence class; no prior structural knowledge was used to derive the consensus model.

The three linked root nodes from which most downstream influences in the graph structure were derived included p-SRC(Tyr416), p-EGFR(Tyr845), and p-PDGFRB(Tyr1009). The model suggests that the EGFR and PDGFRA,B influence one another, with p-EGFR(Tyr1068), p-EGFR(Tyr1173), p-ERBB2(Tyr1221,Tyr1222), and p-KIT(C)(Tyr719) depicted directly downstream of both p-PDGFRB(Tyr1009) and p-EGFR(Tyr845). Notably, PDGFRB has previously been described to heterodimerize and transactivate the EGFR20 in response to PDGF, even in the presence of a PDGFR inhibitor. While others have previously suggested A431 cells lack PDGFR expression, we witnessed bands at the predicted molecular weights using several phospho- and pan-specific antibodies directed at the intracellular region of the receptor (Supplementary Figures 1, 6) as well as with RT-PCR using primers directed at the juxtamembrane region of PDGFR (data not shown).

Interestingly, p-PDGFRA(Tyr849),PDGFRB(Tyr857), the activation loop phosphosite of PDGFRA,B was depicted by the model to lie downstream of p-MET(Tyr1349), a root node, and p-EGFR(Tyr1173), which was downstream of the root nodes p-EGFR(Tyr845) and p-PDGFRB(Tyr1009). p-EGFR(Tyr1173) began to display robust phosphorylation at the 100 ng ml-1 EGF input level, the same level that the activation loop of PDGFRA,B began to display phosphorylation; at low EGF levels Src kinase may mediate the phosphorylation of a subset of PDGFR sites other than Tyr849,Tyr857, while at higher levels the PDGFR kinase itself becomes activated through a mechanism involving or concurrent with the phosphorylation of p-EGFR(Tyr1173).

p-EGFR(Tyr1068), modeled to be upstream of p-EGFR(Tyr1086), p-ERBB4(Tyr1284) and both p-FGFR1(Tyr653,Tyr654) activation loop isoforms, was distinct among EGFR sites in displaying maximal phosphorylation at 5 minutes and sustained phosphorylation amplitude for the duration of the time course. The edge directed from p-EGFR(Tyr1068) to p-FGFR1(Tyr653,Tyr654)(145 kDa) displayed a relatively high edge score (.80; see Supplementary Figure 7 for all consensus Bayesian network edge weights), similar to that displayed between p-EGFR(Tyr1068) and p-EGFR(Tyr1086) (which displayed an edge score of .89), suggesting that EGFR is able to mediate the activation of FGFR1. We speculate that the 145 kDa and 100 kDa forms of FGFR1 represent hyper- and hypoglycosylated forms of the receptor, respectively. Hyper-glycosylation of FGFR1 has been shown to inhibit its interaction with both FGF2 and heparin-derived oligosaccharides21, which was predicted to decrease its activity. Only the short form phosphosite was modeled to have subsequent downstream targets. The only site negatively regulated in the model was p-PDGFRA(Tyr754), which is known to recruit the SHP2 phosphatase22 resulting in dephosphorylation of RASGAP recruitment sites on PDGFRA and B and increased MAPK signaling. Therefore down-regulation of p-PDGFRA(Tyr754) would be predicted to decrease MAPK signaling. Consistent with previous reports, the model suggests that p-SRC(Tyr527), a known inhibitory site of Src kinase, is disconnected from the EGF network23.

To corroborate the Bayesian network results, we also inferred network connectivities using ‘Algorithm for the Reconstruction of Accurate Cellular Networks’ (ARACNe) and ‘Context Likelihood of Relatedness’ (CLR) algorithms24,25. 22 of 24 edges in the Bayesian network were also identified by ARACNe and/or CLR, though as undirected edges since these latter methods are based on mutual information notions (Fig. 6b,c, Supplementary Figures 8, 9). Consistency across network inference methods has been shown experimentally to characterize accurate edge inference26,27, and in our case here data permutation studies show that the topologies inferred by the Bayesian network, ARACNe, and CLR are significant (Supplementary Note 1 and Supplementary Figure 10). Because in the context of proteomic signaling networks it is problematic to make broad assumptions about edge directionality absent extensive prior knowledge (e.g., concerning particular kinase-substrate relationships), we believe that predicting edge directionality using methods such as Bayesian network modeling offers appealing advantage.

DISCUSSION

In contrast to RPAs, MWAs have the ability to reduce the complexity of lysates after arraying, minimizing effort in experimental scale-up. Most of the information of a traditional western blot can be obtained, using 200-fold less protein and antibody. MWAs should be useful for analysis of proteins from cell lines and tissues where lysates are sufficient to print hundreds of MWAs that could be distributed en masse in an analogous manner to spotted DNA microarrays for interrogation with the user's choice of antibodies. The only devices required following printing are commercially available 96-well gaskets and an imager. The ability to obtain information regarding hundreds of proteins with the MWA method should allow advances in our understanding of cell-context-specific networks underlying human disease when combined with appropriate computational modeling methods.

MWAs could also be very useful for large-scale, systematic validation of antibodies. Antibody collections could be systematically verified for selectivity by examining lysates from cells transfected with a cDNA or depleted for the cognate protein by RNAi. The amount of antibody obtained from a single rabbit immunization (~5 mg) would be sufficient for over 100,000 MWAs, thus minimizing lot-tolot variability of polyclonal antibodies. MWAs could be useful for current efforts to build a human protein atlas; samples from tissues used for in situ analyses could be examined with MWAs to verify that signals observed with each antibody resulted from proteins of the predicted molecular weight(s).

The ability to gather dynamic information regarding hundreds of proteins under many conditions poses new challenges for computational modeling. The Bayesian network described here represents direct and/or indirect effects of a given node on other nodes as indicated by high-probability connecting arcs, which are hypothesized to represent relationships of influence among the phosphoproteins in the network. Using prior knowledge to restrict edge directionality across a Bayesian network equivalence class, one can bolster the case for assigning directionality to these edges. To further support a case for interpreting network connections as causal, one could explicitly model the temporal data28 and/or use interventional data15,17 which will be the source of future inquiry.

The timing and amplitude of phosphorylation dynamics observed here coupled with the connectivities modeled in the Bayesian network suggest several candidate sources of RTK-coactivation, each of which may be important in specific cancer contexts: 1) direct dimerization and/or phosphorylation by EGFR or other downstream tyrosine kinases as suggested by the rapid phosphorylation kinetics of Src, ErbB2 and 4, coupled with their close proximity at the top of the network; 2) activation of proteases that activate precursor growth factors or latent RTKs as might be predicted from the delayed phosphorylation amplitudes of FGFR1(100kDa) and MET activation loop sites coupled with their distance from EGFR in the network and; 3) inactivation of tyrosine phosphatases through oxidation by reactive oxygen species29. Phosphorylation of Tyr542 of Shp2 phosphatase displayed the highest fold-change of any site in our analysis; this site has been suggested to relieve inhibition of phosphatase30 activity. The sustained phosphorylation of this and other tyrosine sites at EGF levels ≥ 50 ng/ml suggests that it (and other cysteine-based tyrosine phosphatases) may be inactivated thus unmasking many tyrosine kinase activities. Each of these mechanisms may play distinct roles in the context of cancers that have become resistant to single kinase inhibitors; systems level analysis of other tyrosine kinase-driven cancers may be helpful in revealing appropriate therapeutic targets.

EdSumm AOP

Micro-western arrays combine the advantages of scalability of reverse phase protein arrays and the information content of western blotting, for analyzing protein abundance and modification state with high sensitivity and throughput. The method is demonstrated for analyzing phosphorylation state changes in the EGF receptor signaling network using Bayesian network modeling.

EdSumm issue

Micro-western arrays combine the advantages of scalability of reverse phase protein arrays and the information content of western blotting, for analyzing protein abundance and modification state with high sensitivity and throughput. The method is demonstrated for analyzing phosphorylation state changes in the EGF receptor signaling network using Bayesian network modeling.