Spatio-temporal modeling of signaling protein recruitment to EGFR
© Hsieh et al; licensee BioMed Central Ltd. 2010
Received: 14 July 2009
Accepted: 6 May 2010
Published: 6 May 2010
A stochastic simulator was implemented to study EGFR signal initiation in 3D with single molecule detail. The model considers previously unexplored contributions to receptor-adaptor coupling, such as receptor clustering and diffusive properties of both receptors and binding partners. The agent-based and rule-based approach permits consideration of combinatorial complexity, a problem associated with multiple phosphorylation sites and the potential for simultaneous binding of adaptors.
The model was used to simulate recruitment of four different signaling molecules (Grb2, PLCγ1, Stat5, Shc) to the phosphorylated EGFR tail, with rules based on coarse-grained prediction of spatial constraints. Parameters were derived in part from quantitative immunoblotting, immunoprecipitation and electron microscopy data. Results demonstrate that receptor clustering increases the efficiency of individual adaptor retainment on activated EGFR, an effect that is overridden if crowding is imposed by receptor overexpression. Simultaneous docking of multiple proteins is highly dependent on receptor-adaptor stability and independent of clustering.
Overall, we propose that receptor density, reaction kinetics and membrane spatial organization all contribute to signaling efficiency and influence the carcinogenesis process.
The ErbB or Epidermal Growth Factor Receptor (EGFR) family of receptor tyrosine kinases consists of four members: EGFR (ErbB1), ErbB2, ErbB3, and ErbB4. Under normal physiological conditions, they propagate signals regulating cell proliferation, differentiation, motility and apoptosis. Changes in expression and aberrant activation, especially of EGFR and ErbB2, are associated with a variety of cancers . Upon ligand binding, EGFR undergoes a conformational change that leads to the formation of homodimers (EGFR-EGFR) and heterodimers (i.e., EGFR-ErbB2) . Dimerization induces kinase activation and transphosphorylation of multiple tyrosine residues in receptor cytoplasmic tails [3–5]. The phosphotyrosine residues serve as docking sites for a large number of cytoplasmic adaptor proteins and enzymes . For a given cell type, the specificity and potency of EGFR-mediated intracellular signaling is mediated by the cell's repertoire of phosphotyrosine-binding proteins recruited to the EGFR cytoplasmic tail.
In this work, we use an agent-based model to evaluate the effects of reaction kinetics, steric constraints and receptor clustering on the docking of four EGFR binding partners: Grb2, Shc, Stat5 and PLCγ1. The adaptor Grb2 lacks enzymatic activity and consists of one Src homology (SH) 2 domain and two SH3 domains . Its SH2 domain docks to specific EGFR phosphotyrosine residues and its SH3 domains bind to a Ras guanine nucleotide exchange factor, Sos [8, 9]. The adaptor Shc also binds directly to activated EGFR by two distinct phosphotyrosine interaction domains, an NH2-terminal phosphotyrosine binding (PTB) domain and a COOH-terminal SH2 domain [10, 11]. Recruitment of Grb2 and Shc lead to activation of ERK (extracellular signal regulated kinase) , which translocates into the nucleus and induces gene expression . The transcription factor Stat5 is activated by phosphorylation after docking to EGFR or indirectly through Src-mediated EGFR signaling [14, 15]. Activated Stat5 translocates into the nucleus where it regulates the transcription of selected genes involved in oncogenesis [16, 17]. PLCγ1 has two SH2 domains, one SH3 domain and two pleckstrin homology (PH) domains . It is recruited to phosphorylated EGFR through its SH2 domains, where it serves as a substrate for EGFR kinase activity. Tyrosine phosphorylation of PLCγ1 then leads to an increase in its enzyme activity . PLCγ1 pathway plays a significant role in EGFR-mediated cell signaling, including calcium signaling , receptor endocytosis  and cell motility . Overexpression and hyperactivation of PLCγ1 has been implicated in breast and prostate cancers, and has especially been linked to cancer cell invasion [23, 24].
The process of signaling through ErbB receptors involves highly connected networks of interacting components. Improved understanding of receptor signaling through systems biology approaches has a number of potential practical applications, such as the rational design of drugs to treat cancer . The accuracy of mathematical models relies heavily on quantitative characterization of signaling components and their interactions, such as measurement of expression levels and reaction rate constants. However, the acquisition of quantitative information is no small task, in part because signaling proteins contain multiple phosphorylation sites and may interact with multiple binding partners. Many groups have studied the affinity between EGFR phosphopeptides and the binding domains of Grb2, Shc, STATs, and PLCγ1 using protein microarrays , Surface Plasmon Resonance (SPR) [27–29] and Isothermal Titration Calorimetry (ITC) [30, 31]. These studies provided estimates of dissociation equilibrium constants (Kd) but association and dissociation rate constants of the reactions were typically either not measured or derived indirectly. Moreover, none of these measurements were based upon whole EGFR within lipid bilayers. To understand distinct recruitment behaviors for the different signaling proteins, it is important to arrive at better estimates of their association and dissociation kinetics. This will require new experimental and computational approaches. In an recent experimental development, Morimatsu and colleagues applied single molecule analysis to measure the reaction rate constants of Grb2 with membranes bearing intact, phosphorylated EGFR . In this study, we combined several quantitative experimental approaches, including western blotting analysis and semi-quantitative electron microscopy, to evaluate the kinetics of EGFR phosphorylation and adaptor recruitment to the plasma membrane of EGF-stimulated cells. Rate constants for EGFR phosphorylation/dephosphorylation and adaptor docking are estimated by fitting this data to simulations in our agent-based stochastic model, Signaling Pathways Simulator (SPS) .
Our model specifically considers the phenomenon referred to as combinational complexity, which has been a challenging problem for deterministic mathematical models that employ differential equations to describe cell signaling pathways [34, 35]. For example, because the EGFR becomes phosphorylated on at least nine tyrosine residues during signaling, there are more than 260,000 distinct combinations of these phosphoforms for a dimer of EGFR. Additional molecular diversity can arise when accounting for potential simultaneous interactions of receptor tails with multiple cytoplasmic adaptors. Previous models of ErbB signaling reduce the problem of combinatorial complexity by making several assumptions, including simultaneous phosphorylation and dephosphorylation of receptor tyrosine residues, representation of all tyrosine residues as a single 'virtual phosphorylation site', and exclusion of multiple cytoplasmic adaptors on the same receptor tail based upon competitive binding [34, 36–39]. In an important advance, Blinov and colleagues developed a rules-based model of early EGFR signaling events, capable of evaluating more than 300 molecular species connected through ~4000 unidirectional reactions . Our spatial stochastic model is also "rules-based" and specifically designed to consider largely unexplored contributions of 1) EGFR clustering [33, 40, 41] and anomalous diffusion , 2) distinct temporal patterns of EGFR tyrosine phosphorylation and 3) the potential for multiple adaptors to bind to the same phosphorylated EGFR tail. We refer to the latter concept as "sharing" and base our simulation rules upon the results of coarse-grain molecular docking modeling. In SPS, receptors diffuse in the two dimensional plasma membrane. Rules established for diffusion in and out of defined subdomains of the membrane (protein islands or rafts) provide a mechanism for receptor clustering . We demonstrate that the agent-based spatial model can effectively address problems associated with combinational complexity and make testable predictions about adaptor binding and signaling output that are consistent with novel, quantitative experimental data sets. The simulation results suggest that adaptor sharing is highly dependent on reaction kinetics. The spatial model also predicts receptor clustering results in more efficient adaptor retainment, particularly at normal receptor expression levels.
Coarse-grained molecular docking simulations establish rules for competitive or simultaneous adaptor recruitment to the EGFR cytoplasmic tail
For convenience's sake, conventional models typically assume that the docking of adaptor proteins is a competitive process [34, 36, 37, 43]. Nevertheless, it is possible that neighboring phosphotyrosine residues on the EGFR tail can recruit distinct proteins at the same time , a phenomenon we refer to hereafter as "sharing". In theory, the ability of multiple proteins to dock on the same tail could influence signal transduction efficiency. To address this in our stochastic model, we first sought to establish docking "rules" based upon coarse-grained molecular docking methods.
Docking rules for adaptors on EGFR cytoplasmic tails, as established by coarse-grained molecular docking modeling simulations
Simulation of EGFR phosphorylation/dephosphorylation kinetics
Estimated rate constants for EGFR phosphorylation/dephosphorylation and adaptor recruitment
Receptor tyrosine phosphorylation rates
0.055 (nM × s)-1
0.055 (nM × s)-1
0.0063 (nM × s)-1
0.055 (nM × s)-1
Multiplier when a receptor is phosphorylated.
Multiplier when a receptor is ligand-bound.
Receptor tyrosine dephosphorylation rates
Multiplier when ligand is removed.
Adaptor docking rates
0.0072 (nM × s)-1
0.0055 (nM × s)-1
0.00216 (nM × s)-1
For Shc docking to pY1148
0.00936 (nM × s)-1
For Shc docking to pY1173
0.0056 (nM × s)-1
Adaptor dissociation rates
For Shc dissociation from pY1148
For Shc dissociation from pY1173
Kd of adaptor recruitment to EGFR
For Shc dissociation from pY1148
For Shc dissociation from pY1173
We used the PottersWheel parameter fitting toolbox  to estimate the "fast" and "slow" kinetics of tyrosine phosphorylation. Details of this ODE-based approach are found in Computational Methods. The estimated phosphorylation and dephosphorylation rate constants of the "fast" kinetics category are 0.055 (nM×s)-1 and 0.013/s, respectively. The rate constant estimated for "slow" kinetics phosphorylation is 0.0063 (nM×s)-1, with a dephosphorylation rate 0.0014/s (Table 2).
Stochastic simulations by SPS  were used to validate these estimated parameter values. Based on a value of 4 million receptors per A431 cell, the simulated cell membrane area of 0.49 μm2 contained 1592 EGFR particles. To set up the initial condition with 14% EGFR predimerized in resting A431 cells, we used the receptor conformational flux model from our previous work . In this model, collision between two transiently "open", dimerization-competent receptors leads to ligand-independent dimerization. We used a simulation time step of 25 μs, a random distribution of EGFR, a diffusion rate of 0.09 μm2/s for receptors, and 20 nM EGF. Ligand binding, dissociation, and receptor dimerization rate constants came from Shankaran's model  and are found in Computational Methods. As shown in Figures 2F and 2H, simulations run with the estimated parameter values for the "fast" and "slow" kinetics show close agreement with the experimental data.
Simulation of kinetics of EGFR with its adaptor proteins
Binding constants reported in literature
EGFR (32 PY peptides)
Protein microarray & SPRa
EGFR PY1148 peptide
TRK PY490 peptide
EGFR PY1148 peptide
EGFR PY1173 peptide
EGFR PY peptides
PLCγ1 N+C SH2
EGFR PY peptides
EGFR PY1068 peptide
EGFR PY1086 peptide
EGFR PY1068 peptide
EGFR PY1068 peptide
Activated EGFR in A431 membranes (fractional)
Single molecule analysis
PDGFR PY315 peptide
p85 2 SH2
(YMXM) peptides polyoma middle tumor antigen
Lck Y505 peptide
Effect of receptor clustering on efficiency of signal transduction
Using this approach, we compared the efficiency of adaptors that are retained on EGFR in the randomly distributed and clustered topography. Under both conditions, simulations included 20nM EGF, 100 EGFR and 280 Grb2 (equivalent to approximately 50,000 EGFR and 141,000 Grb2 per cell). In the clustered condition, we compared cluster sizes of 6.57 and 100 receptors per cluster. Figure 4C shows histograms plotting the number of events where Grb2 docked to another receptor within a 50 second interval after dissociating from a previous binding event, with comparisons in the three spatial environments. Receptor clusters of 100 increase the efficiency of Grb2 rebinding to a second EGFR by 6 fold, compared to randomly distributed receptors at this normal expression level. Overall efficiency of receptor coupling during the first 60 seconds is markedly higher in the clustered state (Figure 5B). This is consistent with the concept that an adaptor dissociating from a receptor has two possible outcomes: it can diffuse back into the cytosol or collide with the membrane. According to this scenario, receptor clustering creates a local density, increasing the likelihood that the adaptor will collide productively with a nearby membrane receptor.
Results in Figure 4D show that the increased efficiency of coupling contributed by receptor clustering would be obscured if overexpression creates a high overall density in the membrane. When using conditions applicable to the highly aggressive A431 cancer cell line (4 million receptors), plots for docking efficiency are essentially identical in the random and clustered state. This supports the hypothesis developed from our previous work that both receptor density and membrane spatial organization may be important factors in the carcinogenesis process .
Sharing docking model may suggest more efficient and diverse signaling output
In Figure 5A-B, we conclude by examining the frequency of simultaneous adaptor recruitment to activated EGFR under conditions of normal (50,000) and high (4,000,000) levels of receptor expression. Rules for permitted combinations of adaptors and signaling molecules on a single phosphorylated receptor tail were based upon our coarse-grained docking approach (Figure 1, Table 1). This work builds on that described in previous figures, including rate constants for docking and dissociation, as well as phosphorylation/dephosphorylation. The cytosolic simulation space was populated with identical numbers of the four signaling proteins (56 Grb2, 55 Shc, 59 Stat5, and 154 PLCγ1), while receptors varied from 20 (representing 50,000 receptors in the whole cell) to 1592 (4 million receptors/cell, as in A431 cells) for the case of random topography. To simulate clustered topography for the normal case of 50,000 receptors/cell, the simulation space was expanded 5-fold, with a corresponding increase in receptors and adaptors to match the new cellular volume representation. For each level of receptor, we compared results based upon the "sharing" docking model with results generated using a "competing" docking model. In the latter case, occupancy of a receptor tail at any given time step excluded another signaling protein or adaptor from binding to the same tail.
Simulation results in Figure 5A-B are intriguing in that they predict that the capability for "shared" docking does not significantly affect overall recruitment of adaptors, even in simulations using high density of receptors (either through clustering or overexpression). This result is intuitive in the case of EGFR overexpression, where the number of receptors is five times that the total number of the four adaptors (Figure 5A). This result was initially unexpected in the case where receptors fall well below the level of adaptors (50,000 receptors/cell, Figure 5B). However, this is explained by the prediction that the shared docking model should be profoundly dependent upon the rate constants applied in the simulation for dissociation of proteins from EGFR. This is illustrated in Figure 5C, where we substituted our docking and dissociation rate constants for Grb2, Shc and PLCγ1 with those of Kholodenko's model . For Stat5, we substituted the Kd estimated by Shao et al for docking of Stat3 to the phosphododecapeptide Y1068 . The simulations had 20 receptors in clustered topography in the simulated space (50,000 receptors per cell), 20nM EGF and the same numbers of the four adaptors as before. Figure 5C shows that use of these slower dissociation rate constants results in very large differences in total numbers of adaptors recruited using the two docking models. There are up to 7.5-fold increases in adaptors docked to EGFR at steady state using the sharing model, compared to the competitive model. If the dissociation rate of proteins bound to EGFR is slow, the formation of large complexes on the EGFR tail is likely to be a frequent event.
In this work, we apply agent-based, stochastic model to investigate mechanisms of adaptor proteins recruitment to EGFR as functions of time, receptor conformation, density and spatial distribution. Unique features include the inclusion of docking rules to consider the problem of combinatorial complexity and the consideration of cell membrane heterogeneity. Parameter values used in the model were estimated by fitting to our own western blotting and immunoelectron microscopy data from A431 cells, which provide evidence of distinct phosphorylation kinetics for different EGFR tyrosine residues and distinct behavior for adaptors recruited to phosphorylated receptors. Previous models for studying ErbB receptor signaling represented all tyrosine residues as a single 'virtual phosphorylation site' and assumed competition among cytoplasmic adaptors for receptor binding [34, 36–39]. The rule-based model of Blinov et al  can account more fully for potential molecular diversity; however, like others it is a differential equation model that cannot well describe cell surface heterogeneities, such as microdomains  or anomalous diffusion of surface receptors .
Our experimental data for resting cells were derived from serum-starved and batimastat-treated A431 cells, conditions that control for serum factors or autocrine stimulation (by shedding of EGFR ligands). It is notable that, at t = 0, there are already detectable levels of EGFR tyrosine phosphorylation and adaptors docked to EGFR (Figure 2). We hypothesize that these are contributed from the ligand-independent EGFR dimers as observed by many groups [33, 40, 49, 54–57]. Therefore, we made some assumptions when quantifying the western blotting data of receptor tyrosine residues phosphorylation: First, the observed phosphorylation at resting is entirely contributed from transient dimers formed by encounters between conformationally fluxing receptors. This process is density dependent and estimates of these "constitutive" and unstable dimers in A431 cells range from 5%  to 14% ; we use the latter value for our simulations. Second, we assume that the rise in tyrosine phosphorylation is equivalent in both receptors within the dimers that form (60% participation at peak values after treatment with 20 nM EGF). This may be an overestimate, since it is unknown whether both tails in the asymmetric dimer  are equally capable of phosphorylation. Indeed, there is evidence that receptor phosphorylation achieves values of only 10-35% in mammary epithelial cells , which would be consistent with unequal transphosphorylation by the two kinases in an EGFR homodimer. We also ignored internalization of receptors in the present work, although we acknowledge that this may be a component of the "fast" and "slow" kinetics for the four tyrosine residues whose phosphorylation kinetics we studied. The deterministic models of Wiley [38, 60] and Kholodenko  have considered the importance of EGFR endocytosis, particularly in the contexts of dimer composition and EGFR mutations, and we anticipate adding this feature to the SPS simulation model as we learn additional molecular details about the kinetics and docking characteristics of AP2 and clathrin recruitment.
Based upon parameter fitting, we estimated that the association (kon) and dissociation (koff) rate constants between Grb2 and EGFR are 0.0072 (nM × s)-1 and 9.34/s, respectively. Remarkably, these values are very close to the reported rate constants measured by single-molecule analysis (kon = 0.0022-0.016 (nM × s)-1; koff = 7.5-8.1/s) in the same cell line . We assumed that Shc PDB domain is the predominant means for recruitment to the EGFR, with pY1148 as a preferred site and pY1173 as the secondary binding site [61, 62]. We arrived at estimated parameter values for Shc to these two sites, with a higher kon for pY1148 (= 0.00936 (nM × s)-1) than pY1173 (= 0.0056 (nM × s)-1) and the same koff (13.39/s) for both. The estimated Kd values of adaptor interactions with EGFR in our studies are of the order of 1 μM (Table 2), and these values are close to those in recent reports [26, 32]. This relatively low affinity is consistent with the estimates that koff is large, and kon is small, such that adaptors can form complexes with EGFR and still can rapidly dissociate to limit signaling duration (Grb2/Shc) or propagate signals (Stat5, PLCγ). In this context, receptor clustering would promote rebinding to another active EGFR and provide a mechanism to enhance signaling efficiency despite high dissociation rates (Figure 4 and 5) A high overall density of receptors, typical of cancer cells that overexpress EGFR, also creates conditions of enhanced signaling efficiency.
It is important to note that, although simulations with these values agree well with the experimental observations, they are up to 100 fold higher (ie., lower affinity) than Kd estimates based upon the binding properties of recombinant SH2 domain and PDB domain to target phosphopeptides (see Table 3 and references therein). This is potentially due to "multi-state" interactions between intact receptors and adaptors, as suggested for EGFR-Grb2 interactions measured by single-molecule analysis . Thus, we caution that the estimated reaction rate constants in this work served solely as references for stochastic modeling. More precise rate constants need to be determined based upon novel experimental techniques for measuring binding between intact molecules in a biological context. As shown in Figure 5, the reaction rates have a large impact on the potential for assembly of large complexes on the EGFR tail. If dissociation rates are fast, combinatorial complexity will be low. If dissociation rates are slow, there is much greater chance for a single EGFR tail to recruit multiple signaling components where permitted by the constraints of steric hindrance.
To our knowledge, this is the first attempt to try to establish rules for simultaneous or competitive docking of adaptors to phosphorylated receptor tails, followed by application of these rules in simulations of single molecule behavior. To accomplish this, we relied upon coarse-grained docking predictions. These methods incorporate information about known or predicted structural domains, using homology and protein structure prediction methods. Domains were linked together by an approach previously applied to modeling of large RNA . Each coarse-grained protein structure was then docked, alone or in combination with the other adaptors, to a flexible C-terminus tail of EGFR based upon their shape complementary using PatchDock. The predicted rules are summarized in Table 1. This inexact approach provided several useful exclusionary rules for testing by simulation. For example, PLCγ1 appears to be sufficiently large that two cannot dock to a single receptor tail, despite the availability of 2 distinct binding sites. PLCγ1 docked at pY1173 likely also prevents Shc docking on the same receptor. We emphasize that the limitations of the coarse grain approach raise many uncertainties about these rules. In addition, there are many unknowns related to the asymmetric model of EGFR dimer, particularly the assumption that both kinases in the dimer become activated by conformational switching . All of these assumptions point to the need for additional experimental validation in order to develop more accurate models.
Mathematical modeling is most useful when it suggests new priorities for coupling of experimentation and simulation. From this work, we identify several areas for future development. First, we seek new insight into the mechanisms that drive receptor clustering and microdomain formation. Our simulation results suggest that, when clustering is introduced, adaptors are retained more readily at the plasma membrane. We predict that this efficiency would increase as the receptor cluster size increases (Figure 4C and 5B). At "normal" receptor expression levels, receptor clustering should create high local density and enhance the probability that dissociated adaptor proteins collide quickly with another receptor. However, simulations also suggest that overexpression of receptors have the same effect as receptor clustering, presumably because an overall high density of receptors improves the chances for dissociating adaptors to rebind. This supports the prediction from our previous work that both receptor density and membrane spatial organization contribute to the carcinogenesis process . Second, the coarse-grained docking results suggest that simultaneous adaptor docking to a single receptor should be possible but that the significance of this possibility is highly dependent upon the stability of receptor-adaptor complexes. Related to this is the need for accurate rate constants for receptor-adaptor complexes. We propose that this effort should take precedence over attempts to confirm the coarse grain docking results, since the problem of combinatorial complexity appears to be minimized if the lifetime of receptor-adaptor binding is short. The SPS platform is applicable to modeling these early signaling events and offers capabilities for extending cascades through the cytoplasm and nucleus. It should also be readily adaptable to other protein assembly problems in cells that rely on diffusion and conformational switches. It can explicitly consider aspects of membrane heterogeneity, combinatorial complexity and hierarchy of adaptor binding.
Coarse-grained Protein Modeling
Our primary goal was to 1) approximate the size and general shape of the EGFR tail and four of its docking partners and 2) use this information to arrive at reasonable predictions of combinations that might be accommodated on the EGFR tail simultaneously. The first step was to account for known domain structures within each protein, using SWISS-MODEL for automated homology-based modeling. This approach relies on sequence relationships of protein domains with one or more of known structure. Core models for individual domains (SH2, PTB, PH, etc) were based upon structural templates . Next, protein prediction methods were used to build regions of the structure not available from the templates. These methods combine machine learning methods, evolutionary information in the form of profiles, fragment libraries extracted from the Protein Data Bank (PDB) , and energy functions to predict protein structural features. The protein structure predictors used in this work were 3Dpro from Scratch , MaxSprout  and I-TASSER . 3Dpro predicts protein tertiary structure and outputs PDB files as a Carbon Alpha trace. Maxsprout was used to add backbone and side chain coordinates to obtain an all-atom model. I-TASSER combines both template-based and template-free modeling approaches. All the servers used in this work have been evaluated by the CASP (Critical Assessment of Structure Prediction) experiment . The overall protein structure was then generated by linking these models of sub-domains together by the software, Insight II http://accelrys.com/products/insight/, using an approach similar to that of  for modeling large RNA 3D structure. PatchDock  was then used to simulate docking of two or more structures. No claim is made regarding the accuracy of these predicted structures.
The PDB files of all the predicted protein structures are available upon request. Brief strategies for each structure follow: 1) For the EGFR tail, the target sequence (L1001 to A1210) was imported into SWISS-MODEL and the First Approach Mode was performed without preselected template files. The resulting 3D model structure contains one domain (A1118 to F1176) built with 1F4H as template structure (37% identities) and has missing sequences from L1001 to P1117 and from F1177 to A1210. These two sequences, with additional amino acids for linking (L1001 to V1133 and L1167 to A1210), were then modeled by 3Dpro and MaxSprout and were linked to the domain obtained from SWISS-MODEL by Insight II. The full-length c-terminal of EGFR was then linked to the asymmetric model of the kinase domains of EGFR dimer by Insight II. 2) Shc is a 583-amino acid protein, and the crystal structure of two domains of Shc (M111 to R317 and Q482 to L583) has been solved (1MIL and 1N3H). The protein structure prediction of the two missing sequences was performed with 3Dpro and Maxsprout. Insight II was then used to link these two models with 1MIL and 1N3H together. 3) Stat5 is 794 amino acid residues long and majority of its crystal structure has been solved (1Y1U: S138 to A690). The missing sequences with additional amino acids for linking, from M1 to E150 and from K681 to S794, were modeled by the protein structure predictor, 3Dpro and MaxSprout. The resulting model structures were then linked to 1Y1U by Insight II. 4) PLCγ1 is composed of 1290 amino acids and the crystal structure of its SH3 domain has been solved (1HSQ: 790 to 851). I-TASSER was used to model the two missing sequences (M1 to F800 and F841 to L1290), and the resulting model structures were linked to 1HSQ by Insight II. Grb2 has most of its crystal structure solved with one small missing sequence from L28 to D33. The full Grb2 was modeled by the first approach mode of SWISS-MODEL and the template was the known structure of Grb2 itself (1GRI).
We used a Matlab parameter-fitting toolbox, PottersWheel , which provides interactive modeling including multi-experiment fitting with highly optimized model integration. An ODE (ordinary differential equations) model was first derived to represent features of the EGFR network. Experimental data was then imported and fit to the data by automatically or manually adjusting model parameters to optimize matching between model trajectory and experimental data points . This strategy was used to determine the phosphorylation rate constant and the dephosphorylation rate constant for tyrosine residues 992, 1068, 1148 and 1173), using the following ODE series:
d(L)/dt = -1*k1*R*L + k2*RL;
d(R)/dt = -1*k1*R*L + k2*RL -2*k3*R*R + 2*5*k4*RR - k3*RL*R + k4*RLR;
d(RL)/dt = -1*k3*RL*R + k4*RLR - 2*k3*RL*RL + 2*k4*RLRL + k1*R*L - k2*RL;
d(RR)/dt = k3*R*R - 5*k4*RR - k1*RR*L + k2*RLR - k5*RR + 4*k6*pRR;
d(RLR)/dt = k3*RL*R - k4*RLR - k1*RLR*L + k2*RLRL - 2*k5*RLR + k6*pRLR;
d(RLRL)/dt = k3*RL*RL - k4*RLRL + k1*RLR*L - k2*RLRL -2*k5*RLRL + k6*pRLRL;
d(pRR)/dt = k5*RR - 4*k6*pRR - 3*k5*pRR + 4*k6*pRpR;
d(pRLR)/dt = k5*RLR - k6*pRLR - 2*3*k5*pRLR + k6*pRLpR;
d(pRLRL)/dt = 2*k5*RLRL - k6*pRLRL - 3*2*k5*pRLRL + k6*pRLpRL;
d(pRpR)/dt = 3*k5*pRR - 4*k6*pRpR;
d(pRLpR)/dt = 2*3*k5*pRLR - k6*pRLpR;
d(pRLpRL)/dt = 3*2*k5*pRLRL - k6*pRLpRL;
where k1 is the ligand binding rate constant (= 0.00000331 (# × sec/simspace)-1, fixed), k2 is the ligand dissociation rate constant (= 0.004/s, fixed), k3 is the dimerization rate constant (= 0.014 (# × sec/simspace)-1, fixed), k4 is the dimer dissociation rate constant (= 0.01/s, fixed), k5 is the phosphorylation rate constant (= 0.01 (# × sec/simspace)-1 initially), and k6 is the dephosphorylation rate constant (= 0.005/s initially). L is the number of ligands in the simulated system (= 780, fixed), and R is the number of receptors in the simulated system (= 1592, fixed). LR is the number of ligand-bound monomer, RR is the number of ligand-free dimer, pRR is the number of dimers that have only one phosphorylated receptor, pRpR is the number of dimers of which both receptors are phosphorylated, and so on. The initial number of these complexes is set to be zero in the model. The rate constants of ligand binding, dissociation and receptor dimerization as well as the reaction rate multiplier came from Shankaran's model (Table 2). The simulated space, simspace, represents ~1/2512 of a typical epithelial cell. In the model, 1 (# × sec/simspace)-1 equals to 0.72 (nM × s) -1 based on the assumption that the cytoplasmic water volume of a cell is 3 × 10-12 l . The ODE model was loaded into PottersWheel, and k5 and k6 of each tyrosine residue are estimated by fitting to the corresponding western blotting data (Figure 2F and 2H) both automatically and manually to avoid being trapped in the local optimum.
We also used this approach to estimate the parameters for binding of the four adaptor proteins (Shc, Grb2, Stat5, and PLCγ1) to target tyrosine residues. The ODE series for parameter-fitting was:
d(L)/dt = -1*k1*R*L + k2*RL;
d(R)/dt = -1*k1*R*L + k2*RL -2*k3*R*R + 2*5*k4*RR - k3*RL*R + k4*RLR;
d(RL)/dt = -1*k3*RL*R + k4*RLR - 2*k3*RL*RL + 2*k4*RLRL + k1*R*L - k2*RL;
d(RR)/dt = k3*R*R - 5*k4*RR - k1*RR*L + k2*RLR - k5*RR + 4*k6*pRR;
d(RLR)/dt = k3*RL*R - k4*RLR - k1*RLR*L + k2*RLRL - 2*k5*RLR + k6*pRLR;
d(RLRL)/dt = k3*RL*RL - k4*RLRL + k1*RLR*L - k2*RLRL -2*k5*RLRL + k6*pRLRL;
d(pRR)/dt = k5*RR - 4*k6*pRR - 3*k5*pRR + 4*k6*pRpR - 1*k7*pRR*AP + k8*pRRAP;
d(pRLR)/dt = k5*RLR - k6*pRLR - 2*3*k5*pRLR + k6*pRLpR - 1*k7*pRLR*AP + k8*pRLRAP;
d(pRLRL)/dt = 2*k5*RLRL - k6*pRLRL - 3*2*k5*pRLRL + k6*pRLpRL - 1*k7*pRLRL*AP + k8*pRLRLAP;
d(pRpR)/dt = 3*k5*pRR - 4*k6*pRpR - 1*k7*pRpR*AP + k8*pRpRAP;
d(pRLpR)/dt = 2*3*k5*pRLR - k6*pRLpR - 1*k7*pRLpR*AP + k8*pRLpRAP;
d(pRLpRL)/dt = 3*2*k5*pRLRL - k6*pRLpRL - 1*k7*pRLpRL*AP + k8*pRLpRLAP;
d(AP)/dt = -1*k7*pRR*AP + k8*pRRAP - 1*k7*pRLR*AP + k8*pRLRAP - 1*k7*pRLRL*AP + k8*pRLRLAP - 1*k7*pRpR*AP + k8*pRpRAP - 1*k7*pRLpR*AP + k8*pRLpRAP - 1*k7*pRLpRL*AP + k8*pRLpRLAP;
d(pRRAP)/dt = -1*k8*pRRAP + k7*pRR*AP;
d(pRLRAP)/dt = -1*k8*pRLRAP + k7*pRLR*AP;
d(pRLRLAP)/dt = -1*k8*pRLRLAP + k7*pRLRL*AP;
d(pRpRAP)/dt = -1*k8*pRpRAP + k7*pRpR*AP;
d(pRLpRAP)/dt = -1*k8*pRLpRAP + k7*pRLpR*AP;
d(pRLpRLAP)/dt = -1*k8*pRLpRLAP + k7*pRLpRL*AP;
where k7 is the docking rate constant (= 0.01 (# × sec/simspace)-1 initially), k8 is the rate constant of adaptor proteins dissociated from a receptor (= 0.005/s initially), and the values of k5 and k6 of a tyrosine residue are determined from the previous fitting and are fixed here. AP is the number of adaptor proteins in the simulated space which is determined from experimental data and was estimated to be 56 for Grb2, 55 for Shc, 59 for Stat5 and 154 for PLCγ1. pRLRAP is the number of adaptor protein-bound dimer of which one receptor is ligand-bound and the other receptor is phosphorylated, and so on. The experimental data points used for parameter fitting were derived by converting the values obtained by immunoelectron microscopy for a 3 μm2 area, as shown in Figure 3C-F, to corresponding values in our 0.49 μm2 simulated space. The rate constants of adaptor proteins docking and dissociation are estimated by adjusting k7 and k8 to minimize the distance between the model trajectory (sum of the adaptor protein-bound dimer; i.e. 'pRRAP + pRLRAP + pRLRLAP + pRpRAP +pRLpRAP + pRLpRLAP') and the experimental data points. Rate constants derived by this method are reported in Table 2.
The Signaling Pathways Simulator
The agent-based, stochastic simulator was previously described . In brief, the 3D simulation space is composed of an extracellular domain, plasma membrane and cytosol. In this work, ligands were simulated as a single concentration; however SPS can also treat ligands as individual species. Proteins in the 2D membrane and 3D cytosolic space are represented by sphere-like particles with a radius determined from experimental data and their coarse-grained molecular models. At each time step (typically 25 μs), these particles diffuse and have the potential to react with neighbors. SPS is designed for flexible model development and deployment by a modularized and rule-based approach. It tracks the individual reactions of multistate molecules and accommodates complex situations.
Modeling bio-reactions between EGFR and its adaptor proteins
where k on is the phosphorylation rate constant, k off is the dephosphorylation rate constant (different tyrosine residues have different phosphorylation and dephosphorylation rates) and Δt is the time step.
where k a in this case is a diffusion-limited reaction rate constant, V is the volume of the spatially homogeneous chemical system of interest, and f is a scaling factor used to match the well-mixed chemical system to our agent-based spatially heterogeneous model. Its values are approximated by matching the simulation results with the western blotting and EM experimental data (Figure 3J to 3M) and are approximated to be in the range of 10,000 to 30,000 for our modeling of EGFR-adaptor protein bio-reactions.
where k a is the association rate constant and k d is the dissociation rate constant, neither of which plays a role in DistFac(r) after normalization in equation (5).
The diameters and the diffusion coefficients of the four adaptor proteins of EGFR (parameter a and D in equations (4) to (9)) are determined as follows: the diameters of the four adaptor proteins, Grb2, Shc, PLCγ1, and Stat5, are determined from their protein modeling and are determined to be 6 nm, 11 nm, 11.5 nm, and 15.5 nm, respectively. It has been suggested that the diffusion coefficient of molecules, D, is related to the molecular weight, M, by the relation, . The molecular weight of Grb2, Shc, PLCγ1, and Stat5 are 25 kDa, 50 kDa, 145 kDa, and 80 kDa, respectively http://www.ncbi.nlm.nih.gov/. Assuming the diffusion coefficient of Grb2 to be 100 μm2/s , the diffusion coefficients of Shc, PLCγ1, and Stat5 can then be derived and are approximated to be 80, 56, and 68 μm2/s, respectively.
Diffusion of cytoplasmic proteins is simulated by Brownian motion while diffusion of receptor agents in the membrane is based upon the Constrained Brownian Motion Algorithm . In CBM, overlap between two molecules is not permitted; this is the basis for modeling ligand-receptor, receptor-receptor, and receptor-adaptor protein bio-reactions. The algorithm is a modification of the Gilliespie's approach and is summarized as follows:
The probabilities of reactions are then scaled by , denoted as S. S is defined such that the sum of all probabilities of reactions scaled by S equals to (1- P NR ), which is the probability of one reaction occurring. A random number, X, between 0 and 1 is generated to select an event for the receptor agent during the time step, permitting a specific bio-reaction to occur or else the agent undergoes constrained Brownian motion.
Cell treatments and western blotting analysis
A431 cells were obtained from ATCC and cultured according to ATCC recommendations. Epidermal growth factor (EGF) was from Biomedical Technologies (Stoughton, MA). Prior to experimentation, cells were incubated for 3 hrs in serum-free medium with or without batimastat (a gift of P. McGuire, UNM), followed by lysis in ice cold 1% NP-40 buffer (150 mM NaCl, 50 mM Tris/HCl pH 7.2 with protease inhibitors). Protein concentrations in clarified lysates were measured using the BCA protein assay (Pierce, Rockford, IL). Supernatants were mixed with 6x sample buffer. Proteins were separated by SDS-PAGE and transferred to nitrocellulose. Membranes were blocked and sequentially probed with primary and HRP-conjugated secondary antibodies. EGFR antibodies were from Santa Cruz (La Jolla, CA) and phosphotyrosine-specific EGFR antibodies were from Cell Signaling (Beverly, MA). Shc antibodies were from BD Bioscience (Los Angeles, CA), PLCγ1 antibodies were from Millopore (Bedford, MA), and Stat5 and Grb2 antibodies were from Santa Cruz (Santa Cruz, CA). Immuno-reactive bands were detected by enhanced chemiluminescence (Pierce) and their intensity digitally analyzed following densitometry using MultiGauge software.
Preparation of plasma membrane sheets and gold labeling for TEM
Methods for labeling proteins membrane sheets are described in Yang et al . Cells were grown on glass coverslips, treated as described in legends, then fixed in 0.5% PFA for 10 min. Coverslips were inverted onto EM grids (glow discharged, formvar and poly-L-lysine coated) and ripped to leave plasma membranes on the grid, cytoplasmic face up. After fixation (2% PFA, 20 min), grids were incubated sequentially with primary antibodies and gold-conjugated secondary antibodies. Samples were post-fixed with 2% glutaraldehyde, stained with tannic acid and examined using a Hitachi H-7500 transmission electron microscope (TEM) equipped with a 6.8 megapixel digital camera. A customized plugin for ImageJ was used to acquire positions of gold particles and clustering was analyzed using the Hopkins spatial statistic .
Cells were serum starved for 3 hours in the presence of batimastat, followed by EGF stimulation as indicated in legends. Reactions were halted by transfer to 4°C and rinsing with cold PBS. Cells were scraped off and briefly sonicated; intact cells and debris were sedimented by microcentrifugation (10 minutes). Supernatants were subjected to ultracentrifugation (100,000 g, 1 hr, 4°C) to yield membrane and cytosol fractions. Membrane pellets were dissolved in cold NP-40 lysis buffer. Protein concentrations in fractions were determined by BCA assay (Pierce) to normalize samples prepared for SDS-PAGE.
This work was supported by NIH grants R01 CA119232 (to BW) and P50GM085273 (the New Mexico Spatiotemporal Modeling Center). The SPS platform is accessible by request to the authors. We thank Alexandre Chigaev (UNM) and Chang-Shung Tung (LANL) for insightful discussions.
- Blume-Jensen P, Hunter T: Oncogenic kinase signalling. Nature. 2001, 411: 355-365. 10.1038/35077225View ArticlePubMedGoogle Scholar
- Schlessinger J: Ligand-induced, receptor-mediated dimerization and activation of EGF receptor. Cell. 2002, 110: 669-672. 10.1016/S0092-8674(02)00966-2View ArticlePubMedGoogle Scholar
- Yarden Y, Schlessinger J: Epidermal growth factor induces rapid, reversible aggregation of the purified epidermal growth factor receptor. Biochemistry. 1987, 26: 1443-1451. 10.1021/bi00379a035View ArticlePubMedGoogle Scholar
- Carpenter G, King L, Cohen S: Epidermal growth factor stimulates phosphorylation in membrane preparations in vitro. Nature. 1978, 276: 409-410. 10.1038/276409a0View ArticlePubMedGoogle Scholar
- Ushiro H, Cohen S: Identification of phosphotyrosine as a product of epidermal growth factor-activated protein kinase in A-431 cell membranes. J Biol Chem. 1980, 255: 8363-8365.PubMedGoogle Scholar
- Schulze WX, Deng L, Mann M: Phosphotyrosine interactome of the ErbB-receptor kinase family. Molecular Systems Biology. 2005, 1: 2005.0008- 10.1038/msb4100012PubMed CentralView ArticlePubMedGoogle Scholar
- Lowenstein EJ, Daly RJ, Batzer AG, Li W, Margolis B, Lammers R, Ullrich A, Skolnik EY, Bar-Sagi D, Schlessinger J: The SH2 and SH3 domain-containing protein GRB2 links receptor tyrosine kinases to ras signaling. Cell. 1992, 70: 431-442. 10.1016/0092-8674(92)90167-BView ArticlePubMedGoogle Scholar
- Chardin P, Camonis JH, Gale NW, van Aelst L, Schlessinger J, Wigler MH, Bar-Sagi D: Human Sos1: a guanine nucleotide exchange factor for Ras that binds to GRB2. Science. 1993, 260: 1338-1343. 10.1126/science.8493579View ArticlePubMedGoogle Scholar
- Rozakis-Adcock M, Fernley R, Wade J, Pawson T, Bowtell D: The SH2 and SH3 domains of mammalian Grb2 couple the EGF receptor to the Ras activator mSos1. Nature. 1993, 363: 83-85. 10.1038/363083a0View ArticlePubMedGoogle Scholar
- Rozakis-Adcock M, McGlade J, Mbamalu G, Pelicci G, Daly R, Li W, Batzer A, Thomas S, Brugge J, Pelicci PG, et al.: Association of the Shc and Grb2/Sem5 SH2-containing proteins is implicated in activation of the Ras pathway by tyrosine kinases. Nature. 1992, 360: 689-692. 10.1038/360689a0View ArticlePubMedGoogle Scholar
- Bonfini L, Migliaccio E, Pelicci G, Lanfrancone L, Pelicci PG: Not all Shc's roads lead to Ras. Trends Biochem Sci. 1996, 21: 257-261.View ArticlePubMedGoogle Scholar
- Burack WR, Shaw AS: Live Cell Imaging of ERK and MEK: simple binding equilibrium explains the regulated nucleocytoplasmic distribution of ERK. J Biol Chem. 2005, 280: 3832-3837. 10.1074/jbc.M410031200View ArticlePubMedGoogle Scholar
- Yang SH, Whitmarsh AJ, Davis RJ, Sharrocks AD: Differential targeting of MAP kinases to the ETS-domain transcription factor Elk-1. Embo J. 1998, 17: 1740-1749. 10.1093/emboj/17.6.1740PubMed CentralView ArticlePubMedGoogle Scholar
- Kloth MT, Catling AD, Silva CM: Novel activation of STAT5b in response to epidermal growth factor. J Biol Chem. 2002, 277: 8693-8701. 10.1074/jbc.M111884200View ArticlePubMedGoogle Scholar
- Olayioye MA, Neve RM, Lane HA, Hynes NE: The ErbB signaling network: receptor heterodimerization in development and cancer. Embo J. 2000, 19: 3159-3167. 10.1093/emboj/19.13.3159PubMed CentralView ArticlePubMedGoogle Scholar
- Darnell JE: STATs and gene regulation. Science. 1997, 277: 1630-1635. 10.1126/science.277.5332.1630View ArticlePubMedGoogle Scholar
- Bowman T, Garcia R, Turkson J, Jove R: STATs in oncogenesis. Oncogene. 2000, 19: 2474-2488. 10.1038/sj.onc.1203527View ArticlePubMedGoogle Scholar
- Wahl M, Carpenter G: Selective phospholipase C activation. Bioessays. 1991, 13: 107-113. 10.1002/bies.950130303View ArticlePubMedGoogle Scholar
- Spivak-Kroizman T, Rotin D, Pinchasi D, Ullrich A, Schlessinger J, Lax I: Heterodimerization of c-erbB2 with different epidermal growth factor receptor mutants elicits stimulatory or inhibitory responses. J Biol Chem. 1992, 267: 8056-8063.PubMedGoogle Scholar
- Wahl MI, Nishibe S, Kim JW, Kim H, Rhee SG, Carpenter G: Identification of two epidermal growth factor-sensitive tyrosine phosphorylation sites of phospholipase C-gamma in intact HSC-1 cells. J Biol Chem. 1990, 265: 3944-3948.PubMedGoogle Scholar
- Choi JH, Park JB, Bae SS, Yun S, Kim HS, Hong WP, Kim IS, Kim JH, Han MY, Ryu SH, et al.: Phospholipase C-gamma1 is a guanine nucleotide exchange factor for dynamin-1 and enhances dynamin-1-dependent epidermal growth factor receptor endocytosis. J Cell Sci. 2004, 117: 3785-3795. 10.1242/jcs.01220View ArticlePubMedGoogle Scholar
- Chen P, Xie H, Sekar MC, Gupta K, Wells A: Epidermal growth factor receptor-mediated cell motility: phospholipase C activity is required, but mitogen-activated protein kinase activity is not sufficient for induced cell movement. J Cell Biol. 1994, 127: 847-857. 10.1083/jcb.127.3.847View ArticlePubMedGoogle Scholar
- Thomas SM, Coppelli FM, Wells A, Gooding WE, Song J, Kassis J, Drenning SD, Grandis JR: Epidermal growth factor receptor-stimulated activation of phospholipase Cgamma-1 promotes invasion of head and neck squamous cell carcinoma. Cancer Res. 2003, 63: 5629-5635.PubMedGoogle Scholar
- Wells A, Grandis JR: Phospholipase C-gamma1 in tumor progression. Clin Exp Metastasis. 2003, 20: 285-290. 10.1023/A:1024088922957View ArticlePubMedGoogle Scholar
- Yarden Y, Sliwkowski MX: Untangling the ErbB signalling network. Nat Rev Mol Cell Biol. 2001, 2: 127-137. 10.1038/35052073View ArticlePubMedGoogle Scholar
- Jones RB, Gordus A, Krall JA, MacBeath G: A quantitative protein interaction network for the ErbB receptors using protein microarrays. Nature. 2006, 439: 168-174. 10.1038/nature04177View ArticlePubMedGoogle Scholar
- Batzer AG, Rotin D, Urena JM, Skolnik EY, Schlessinger J: Hierarchy of binding sites for Grb2 and Shc on the epidermal growth factor receptor. Mol Cell Biol. 1994, 14: 5192-5201.PubMed CentralView ArticlePubMedGoogle Scholar
- Chattopadhyay A, Vecchi M, Ji Q, Mernaugh R, Carpenter G: The role of individual SH2 domains in mediating association of phospholipase C-gamma1 with the activated EGF receptor. J Biol Chem. 1999, 274: 26091-26097. 10.1074/jbc.274.37.26091View ArticlePubMedGoogle Scholar
- Chook YM, Gish GD, Kay CM, Pai EF, Pawson T: The Grb2-mSos1 complex binds phosphopeptides with higher affinity than Grb2. J Biol Chem. 1996, 271: 30472-30478. 10.1074/jbc.271.48.30472View ArticlePubMedGoogle Scholar
- Mandiyan S, Schumacher C, Cioffi C, Sharif H, Yuryev A, Lappe R, Monia B, Hanson S, Goff S, Wennogle L: Molecular and cellular characterization of baboon C-Raf as a target for antiproliferative effects of antisense oligonucleotides. Antisense Nucleic Acid Drug Dev. 1997, 7: 539-548.View ArticlePubMedGoogle Scholar
- Lemmon MA, Ladbury JE, Mandiyan V, Zhou M, Schlessinger J: Independent binding of peptide ligands to the SH2 and SH3 domains of Grb2. J Biol Chem. 1994, 269: 31653-31658.PubMedGoogle Scholar
- Morimatsu M, Takagi H, Ota KG, Iwamoto R, Yanagida T, Sako Y: Multiple-state reactions between the epidermal growth factor receptor and Grb2 as observed by using single-molecule analysis. Proc Natl Acad Sci USA. 2007, 104: 18013-18018. 10.1073/pnas.0701330104PubMed CentralView ArticlePubMedGoogle Scholar
- Hsieh M, Yang S, Raymond-Stinz M, Steinberg S, Vlachos D, Shu W, Wilson B, Edwards JS: Stochastic Simulations of ErbB Homo and Hetero-Dimerization: Potential Impacts of Receptor Conformational State and Spatial Segregation. IET Systems Biology. 2008, 2: 256-272. 10.1049/iet-syb:20070073View ArticlePubMedGoogle Scholar
- Birtwistle MR, Hatakeyama M, Yumoto N, Ogunnaike BA, Hoek JB, Kholodenko BN: Ligand-dependent responses of the ErbB signaling network: experimental and modeling analyses. Mol Syst Biol. 2007, 3: 144- 10.1038/msb4100188PubMed CentralView ArticlePubMedGoogle Scholar
- Blinov ML, Faeder JR, Goldstein B, Hlavacek WS: A network model of early events in epidermal growth factor receptor signaling that accounts for combinatorial complexity. Biosystems. 2006, 83: 136-151. 10.1016/j.biosystems.2005.06.014View ArticlePubMedGoogle Scholar
- Kholodenko BN, Demin OV, Moehren G, Hoek JB: Quantification of Short Term Signaling by the Epidermal Growth Factor Receptor. J Biol Chem. 1999, 274: 30169-30181. 10.1074/jbc.274.42.30169View ArticlePubMedGoogle Scholar
- Schoeberl B, Eichler-Jonsson C, Gilles ED, Muller G: Computational modeling of the dynamics of the MAP kinase cascade activated by surface and internalized EGF receptors. Nat Biotechnol. 2002, 20: 370-375. 10.1038/nbt0402-370View ArticlePubMedGoogle Scholar
- Hendriks BS, Opresko LK, Wiley HS, Lauffenburger D: Quantitative analysis of HER2-mediated effects on HER2 and epidermal growth factor receptor endocytosis: distribution of homo- and heterodimers depends on relative HER2 levels. J Biol Chem. 2003, 278: 23343-23351. 10.1074/jbc.M300477200View ArticlePubMedGoogle Scholar
- Shankaran H, Wiley HS, Resat H: Modeling the effects of HER/ErbB1-3 coexpression on receptor dimerization and biological response. Biophys J. 2006, 90: 3993-4009. 10.1529/biophysj.105.080580PubMed CentralView ArticlePubMedGoogle Scholar
- Sako Y, Minoghchi S, Yanagida T: Single-molecule imaging of EGFR signalling on the surface of living cells. Nat Cell Biol. 2000, 2: 168-172. 10.1038/35004044View ArticlePubMedGoogle Scholar
- Yang S, Raymond-Stintz MA, Ying W, Zhang J, Lidke DS, Steinberg SL, Williams L, Oliver JM, Wilson BS: Mapping ErbB receptors on breast cancer cell membranes during signal transduction. J Cell Sci. 2007, 120: 2763-2773. 10.1242/jcs.007658View ArticlePubMedGoogle Scholar
- Smith PR, Morrison IE, Wilson KM, Fernandez N, Cherry RJ: Anomalous diffusion of major histocompatibility complex class I molecules on HeLa cells determined by single particle tracking. Biophys J. 1999, 76: 3331-3344. 10.1016/S0006-3495(99)77486-2PubMed CentralView ArticlePubMedGoogle Scholar
- Resat H, Ewald JA, Dixon DA, Wiley HS: An Integrated Model of Epidermal Growth Factor Receptor Trafficking and Signal Transduction. Biophys J. 2003, 85: 730-743. 10.1016/S0006-3495(03)74516-0PubMed CentralView ArticlePubMedGoogle Scholar
- Schneidman-Duhovny D, Inbar Y, Nussinov R, Wolfson HJ: PatchDock and SymmDock: servers for rigid and symmetric docking. Nucleic Acids Res. 2005, 33: W363-367. 10.1093/nar/gki481PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang X, Gureasko J, Shen K, Cole PA, Kuriyan J: An allosteric mechanism for activation of the kinase domain of epidermal growth factor receptor. Cell. 2006, 125: 1137-1149. 10.1016/j.cell.2006.05.013View ArticlePubMedGoogle Scholar
- Lee NY, Koland JG: Conformational changes accompany phosphorylation of the epidermal growth factor receptor C-terminal domain. Protein Sci. 2005, 14: 2793-2803. 10.1110/ps.051630305PubMed CentralView ArticlePubMedGoogle Scholar
- Wilson BS, Steinberg SL, Liederman K, Pfeiffer JR, Surviladze Z, Zhang J, Samelson LE, Yang LH, Kotula PG, Oliver JM: Markers for detergent-resistant lipid rafts occupy distinct and dynamic domains in native membranes. Mol Biol Cell. 2004, 15: 2580-2592. 10.1091/mbc.E03-08-0574PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang J, Leiderman K, Pfeiffer JR, Wilson BS, Oliver JM, Steinberg SL: Characterizing the topography of membrane receptors and signaling molecules from spatial patterns obtained using nanometer-scale electron-dense probes and electron microscopy. Micron. 2006, 37: 14-34. 10.1016/j.micron.2005.03.014View ArticlePubMedGoogle Scholar
- Martin-Fernandez M, Clarke DT, Tobin MJ, Jones SV, Jones GR: Preformed oligomeric epidermal growth factor receptors undergo an ectodomain structure change during signaling. Biophys J. 2002, 82: 2415-2427. 10.1016/S0006-3495(02)75585-9PubMed CentralView ArticlePubMedGoogle Scholar
- Maiwald T, Timmer J: Dynamical modeling and multi-experiment fitting with PottersWheel. Bioinformatics. 2008, 24: 2037-2043. 10.1093/bioinformatics/btn350PubMed CentralView ArticlePubMedGoogle Scholar
- Shao H, Cheng HY, Cook RG, Tweardy DJ: Identification and characterization of signal transducer and activator of transcription 3 recruitment sites within the epidermal growth factor receptor. Cancer Res. 2003, 63: 3923-3930.PubMedGoogle Scholar
- Lillemeier BF, Pfeiffer JR, Surviladze Z, Wilson BS, Davis MM: Plasma membrane-associated proteins are clustered into islands attached to the cytoskeleton. Proc Natl Acad Sci USA. 2006, 103: 18992-18997. 10.1073/pnas.0609009103PubMed CentralView ArticlePubMedGoogle Scholar
- Murase K, Fujiwara T, Umemura Y, Suzuki K, Iino R, Yamashita H, Saito M, Murakoshi H, Ritchie K, Kusumi A: Ultrafine membrane compartments for molecular diffusion as revealed by single molecule techniques. Biophys J. 2004, 86: 4075-4093. 10.1529/biophysj.103.035717PubMed CentralView ArticlePubMedGoogle Scholar
- Yu X, Sharma KD, Takahashi T, Iwamoto R, Mekada E: Ligand-independent dimer formation of epidermal growth factor receptor (EGFR) is a step separable from ligand-induced EGFR signaling. Mol Biol Cell. 2002, 13: 2547-2557. 10.1091/mbc.01-08-0411PubMed CentralView ArticlePubMedGoogle Scholar
- Lidke DS, Nagy P, Barisas BG, Heintzmann R, Post JN, Lidke KA, Clayton AH, Arndt-Jovin DJ, Jovin TM: Imaging molecular interactions in cells by dynamic and static fluorescence anisotropy (rFLIM and emFRET). Biochem Soc Trans. 2003, 31: 1020-1027. 10.1042/BST0311020View ArticlePubMedGoogle Scholar
- Liu P, Sudhaharan T, Koh RM, Hwang LC, Ahmed S, Maruyama IN, Wohland T: Investigation of the dimerization of proteins from the epidermal growth factor receptor family by single wavelength fluorescence cross-correlation spectroscopy. Biophys J. 2007, 93: 684-698. 10.1529/biophysj.106.102087PubMed CentralView ArticlePubMedGoogle Scholar
- Saffarian S, Li Y, Elson EL, Pike LJ: Oligomerization of the EGF receptor investigated by live cell fluorescence intensity distribution analysis. Biophys J. 2007, 93: 1021-1031. 10.1529/biophysj.107.105494PubMed CentralView ArticlePubMedGoogle Scholar
- Wiley HS: Anomalous binding of epidermal growth factor to A431 cells is due to the effect of high receptor densities and a saturable endocytic system. J Cell Biol. 1988, 107: 801-810. 10.1083/jcb.107.2.801View ArticlePubMedGoogle Scholar
- Shankaran H, Zhang Y, Opresko L, Resat H: Quantifying the effects of EGFR-HER2 co-expression on HER activation and trafficking. Biochem Biophys Res Commun. 2008, 371: 220-224. 10.1016/j.bbrc.2008.04.043PubMed CentralView ArticlePubMedGoogle Scholar
- Wiley HS, Shvartsman SY, Lauffenburger DA: Computational modeling of the EGF-receptor system: a paradigm for systems biology. Trends Cell Biol. 2003, 13: 43-50. 10.1016/S0962-8924(02)00009-0View ArticlePubMedGoogle Scholar
- Sakaguchi K, Okabayashi Y, Kido Y, Kimura S, Matsumura Y, Inushima K, Kasuga M: Shc phosphotyrosine-binding domain dominantly interacts with epidermal growth factor receptors and mediates Ras activation in intact cells. Mol Endocrinol. 1998, 12: 536-543. 10.1210/me.12.4.536View ArticlePubMedGoogle Scholar
- Okabayashi Y, Sugimoto Y, Totty NF, Hsuan J, Kido Y, Sakaguchi K, Gout I, Waterfield MD, Kasuga M: Interaction of Shc with adaptor protein adaptins. J Biol Chem. 1996, 271: 5265-5269. 10.1074/jbc.271.9.5265View ArticlePubMedGoogle Scholar
- Tung CS, Joseph S, Sanbonmatsu KY: All-atom homology model of the Escherichia coli 30S ribosomal subunit. Nat Struct Biol. 2002, 9: 750-755. 10.1038/nsb841View ArticlePubMedGoogle Scholar
- Moult J: A decade of CASP: progress, bottlenecks and prognosis in protein structure prediction. Curr Opin Struct Biol. 2005, 15: 285-289. 10.1016/j.sbi.2005.05.011View ArticlePubMedGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28: 235-242. 10.1093/nar/28.1.235PubMed CentralView ArticlePubMedGoogle Scholar
- Cheng J, Randall AZ, Sweredoski MJ, Baldi P: SCRATCH: a protein structure and structural feature prediction server. Nucleic Acids Res. 2005, 33: W72-76. 10.1093/nar/gki396PubMed CentralView ArticlePubMedGoogle Scholar
- Holm L, Sander C: Database algorithm for generating protein backbone and side-chain co-ordinates from a C alpha trace application to model building and detection of co-ordinate errors. J Mol Biol. 1991, 218: 183-194. 10.1016/0022-2836(91)90883-8View ArticlePubMedGoogle Scholar
- Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008, 9: 40- 10.1186/1471-2105-9-40PubMed CentralView ArticlePubMedGoogle Scholar
- Moult J, Fidelis K, Zemla A, Hubbard T: Critical assessment of methods of protein structure prediction (CASP): round IV. Proteins. 2001, 2-7. Suppl 5,
- Stoeber B, Liepmann D: Microfabricated microdialysis microneedles for continuous medicalmonitoring. 1st Ann Int IEEE-EMBS Special Topic Conf on Microtechnologies in Medicine and Biology; October; Lyon, France. 2000, 224-228. full_text.Google Scholar
- Mandiyan V, O'Brien R, Zhou M, Margolis B, Lemmon MA, Sturtevant JM, Schlessinger J: Thermodynamic studies of SHC phosphotyrosine interaction domain recognition of the NPXpY motif. J Biol Chem. 1996, 271: 4770-4775. 10.1074/jbc.271.9.4770View ArticlePubMedGoogle Scholar
- Ladbury JE, Lemmon MA, Zhou M, Green J, Botfield MC, Schlessinger J: Measurement of the binding of tyrosyl phosphopeptides to SH2 domains: a reappraisal. Proc Natl Acad Sci USA. 1995, 92: 3199-3203. 10.1073/pnas.92.8.3199PubMed CentralView ArticlePubMedGoogle Scholar
- Felder S, Zhou M, Hu P, Urena J, Ullrich A, Chaudhuri M, White M, Shoelson SE, Schlessinger J: SH2 domains exhibit high-affinity binding to tyrosine-phosphorylated peptides yet also exhibit rapid dissociation and exchange. Mol Cell Biol. 1993, 13: 1449-1455.PubMed CentralView ArticlePubMedGoogle Scholar
- Lemmon MA, Ladbury JE: Thermodynamic studies of tyrosyl-phosphopeptide binding to the SH2 domain of p56lck. Biochemistry. 1994, 33: 5070-5076. 10.1021/bi00183a010View ArticlePubMedGoogle Scholar