Simulating the impact of a molecular 'decision-process' on cellular phenotype and multicellular patterns in brain tumors

Experimental evidence indicates that human brain cancer cells proliferate or migrate, yet do not display both phenotypes at the same time. Here, we present a novel computational model simulating this cellular decision-process leading up to either phenotype based on a molecular interaction network of genes and proteins. The model's regulatory network consists of the epidermal growth factor receptor (EGFR), its ligand transforming growth factor-alpha (TGF alpha), the downstream enzyme phospholipaseC-gamma (PLC gamma) and a mitosis-associated response pathway. This network is activated by autocrine TGF alpha secretion, and the EGFR-dependent downstream signaling this step triggers, as well as modulated by an extrinsic nutritive glucose gradient. Employing a framework of mass action kinetics within a multiscale agent-based environment, we analyse both the emergent multicellular behavior of tumor growth and the single-cell molecular profiles that change over time and space. Our results show that one can indeed simulate the dichotomy between cell migration and proliferation based solely on an EGFR decision network. It turns out that these behavioral decisions on the single cell level impact the spatial dynamics of the entire cancerous system. Furthermore, the simulation results yield intriguing experimentally testable hypotheses also on the sub-cellular level such as spatial cytosolic polarization of PLC gamma towards an extrinsic chemotactic gradient. Implications of these results for future works, both on the modeling and experimental side are discussed.     

Quantitative proteomics and phosphoproteomics reveal novel insights into complexity and dynamics of the EGFR signaling network

The epidermal growth factor receptor (EGFR/ErbB1/Her1) belongs to the ErbB family of receptor tyrosine kinases (RTKs) and is a key player in the regulation of cell proliferation, differentiation, survival, and migration. Overexpression and mutational changes of EGFR have been identified in a variety of human cancers and the regulation of EGFR signaling plays a critical role in tumor development and progression. Due to its biological significance the EGFR signaling network is a widely used model system for the development of analytical techniques. Novel quantitative proteomics and phosphoproteomics approaches play an important role in the characterization of signaling pathways in a time and stimulus dependent manner. Recent studies discussed in this review provide new insights into different aspects of EGFR signal transduction, such as regulation and dynamics of its phosphorylation sites, association with interaction partners and identification of regulated phosphoproteins. Correlation of data from functional proteomics studies with results from other fields of signal transduction research by systems biology will be necessary to integrate and translate these findings into successful clinical applications.     

Unligated epidermal growth factor receptor forms higher order oligomers within microclusters on A431 cells that are sensitive to tyrosine kinase inhibitor binding

Characterization of the association states of the unligated epidermal growth factor receptor (EGFR) is important in understanding the mechanism of EGFR tyrosine kinase activation in a tumor cell environment. We analyzed, in detail, the association states of unligated, immunotagged EGFR on the surface of intact epidermoid carcinoma A431 cells, using AlexaFluor488 and AlexaFluor546 anti-EGFR antibody, mAb528, as probes. Image correlation microscopy revealed the presence of unligated EGFR in submicron scale clusters containing an average of 10-30 receptors (mean cluster density = 32 +/- 9 clusters per square micron). Lifetime-based F?rster resonance energy transfer (FRET) techniques as a function of acceptor:donor labeling ratio disclosed a clustering of the unligated EGFR in clusters containing an average of four receptors on the nanometer (<10 nm) scale. The relationship between the nanoscale and submicron scale associations was determined using a new analysis that combines nanoscale information from lifetime-detected FRET imaging with submicron scale information obtained with image correlation microscopy. This analysis revealed the presence of monomers (or small oligomers) and larger clusters containing 15-30 receptors that were partially associated on the sub-10 nm scale. Pretreatment of the cells with the tyrosine kinase inhibitor AG1478 caused a partial dispersal of the submicron clusters (mean cluster density = 85 +/- 15 clusters per square micron; mean degree of association = 4-10 receptors per cluster) and reduced the level of FRET down to our limit of detection. These results are consistent with a higher order nanoscale receptor organization of the unligated receptor population that is partially controlled by the kinase domains. The ramifications of the results to mechanisms of EGFR activation in a tumor cell environment are discussed.     

Imaging epidermal growth factor receptor phosphorylation in human colorectal cancer cells and human tissues

In tumor cells, high phosphorylation levels of receptor tyrosine kinases may occur in the absence of exogenous ligands due to autocrine signaling or enhanced tyrosine kinase activity. Here we show that the phosphorylation state of the endogenous epidermal growth factor receptor (EGFR) can be quantitatively imaged in tumor cells and tissues by detecting fluorescence resonance energy transfer between fluorophores conjugated to antibodies against the receptor and phosphotyrosine, respectively. Five different human colorectal cell lines were analyzed for activity and expression of EGFR. All cell lines exhibited basal EGFR phosphorylation under serum starvation conditions. Phosphorylation levels increased after stimulation with EGF or pervanadate, dependent on the level of basal EGFR phosphorylation in the respective cell lines. This basal activity correlated inversely with receptor expression. Using the acceptor photobleaching fluorescence resonance energy transfer imaging approach, a significantly higher phosphorylation state of EGFR was also found in resected human colorectal tumor samples as compared with adjacent healthy tissue. Imaging of EGFR phosphorylation may thus serve as a valuable tool to investigate the role of receptor tyrosine kinase activity in malignant cell growth.     

Epidermal growth factor receptor induced apoptosis: potentiation by inhibition of Ras signaling

Previous studies have shown that certain tumor cell lines which naturally express high levels of the epidermal growth factor receptor (EGFR) undergo apoptosis when exposed to epidermal growth factor. Whether this phenomenon is a direct result of receptor overexpression or some other genetic alteration renders these cells sensitive to apoptosis is yet to be established. We show that experimentally increasing the level of EGFR expression predictably leads to apoptosis in a variety of cell types which requires an active tyrosine kinase but not EGFR autophosphorylation sites. Expression of a dominant negative Ras mutant in EGFR overexpressing cells results in a significant potentiation of EGFR induced apoptosis suggesting that Ras activation is a key survival signal generated by the EGFR. We propose that potentiation of EGFR induced apoptosis by dominant negative Ras results, at least in part, by a block of Akt activation.     

Lyso-GM3, its dimer,and multimer: their synthesis,and their effect on epidermal growth factor-induced receptor tyrosine kinase

Glycosphingolipids, particularly gangliosides, are known to modulate growth factor receptor tyrosine kinase. A well-documented example is the inhibitory effect of GM3 on kinase associated with epidermal growth factor receptor (EGFR) in human epidermoid carcinoma A431 cells. Lyso-GM3 was detected as a minor component in A431 cells, and may function as an auxiliary factor in GM3-dependent inhibition of EGFR. We studied the inhibitory effect of chemically synthesized GM3, lyso-GM3, and its derivatives, on EGFR function, based on their interaction in membrane microdomain, with the following major findings: (1) GM3, EGFR, and caveolin coexist, but tetraspanins CD9 and CD82 are essentially absent, within the same low-density membrane fraction, separated by sucrose density gradient ultracentrifugation. (2) Strong interaction between EGFR and GM3 was indicated by increasing binding of EGFR to GM3-coated polystyrene beads, in a GM3 dose-dependent manner. Confocal microscopy results suggested that three components in the microdomain (GM3, EGFR, and caveolin) are closely associated. (3) Lyso-GM3 or lyso-GM3 dimer strongly inhibited EGFR kinase activity, in a dose-dependent manner, while lyso-GM3 trimer and tetramer did not. >50 μM lyso-GM3 was cytolytic, while >50 μM lyso-GM3 dimer was not cytolytic, yet inhibited EGFR kinase strongly. Thus, lyso-GM3 and its dimer exert an auxiliary effect on GM3-induced inhibition of EGFR kinase and cell growth, and lyso-GM3 dimer may be a good candidate for pharmacological inhibitor of epidermal tumor growth.     

A Mathematical Model for the Effects of HER2 Overexpression on Cell Proliferation in Breast Cancer

We present a mathematical model to study the effects of HER2 over-expression on cell proliferation in breast cancer. The model illustrates the proliferative behavior of cells as a function of HER2 and EGFR receptors numbers, and the growth factor EGF. This mathematical model comprises kinetic equations describing the cell surface binding of EGF growth factor to EGFR and HER2 receptors, coupled to a model for the dependence of cell proliferation rate on growth factor receptors binding. The simulation results from this model predict: (1) a growth advantage associated with excess HER2 receptors; (2) that HER2-over-expression is an insufficient parameter to predict the proliferation response of cancer cells to epidermal growth factors; and (3) the EGFR receptor expression level in HER2-over-expressing cells plays a key role in mediating the proliferation response to receptor-ligand signaling. This mathematical model also elucidates the interaction and roles of other model parameters in determining cell proliferation rate of HER2-over-expressing cells.     

The linoleic acid metabolite, 13-HpODE augments the phosphorylation of EGF receptor and SHP-2 leading to their increased association.

In previous studies with Syrian hamster embryo fibroblasts, we found that a specific lipoxygenase metabolite of linoleic acid, 13(S)-hydroperoxyoctadecadienoic acid (HpODE), enhanced epidermal growth factor (EGF) signal transduction in a tumor suppressor gene plus phenotype (supB+); with a diminished response to 13(S)-HpODE in a tumor suppressor gene minus phenotype (supB-). This differential response was attributed to differences in the rate of EGF receptor (EGFR) dephosphorylation. To further define the molecular basis for these observations, in this report we examine the interaction of phosphorylated EGFR with the SH2 domain-containing protein tyrosine phosphatase, SHP-2, a positive modulator of EGF dependent cell growth. SHP-2 associated with phosphorylated EGFR to a greater extent in supB+ cells when compared to supB-. This differential association could not be accounted for by differences between suppressor gene phenotypes in SHP-2 protein level or mutations in the molecular sequence. The addition of 13(S)-HpODE stimulated a concentration-dependent increase in EGF-dependent phosphorylation of SHP-2 and its association with EGFR. A more dramatic response was observed in the supB+ cells. Differences in SHP-2 interaction with EGFR may account, in part, for phenotypic differences in the growth rates and responsiveness to EGF between the supB+ and supB- cells. EGFR-SHP-2 association appears to play an important role in the regulation of EGFR signal transduction.     

Plasma membrane domain organization regulates EGFR signaling in tumor cells

Macromolecular complexes exhibit reduced diffusion in biological membranes; however, the physiological consequences of this characteristic of plasma membrane domain organization remain elusive. We report that competition between the galectin lattice and oligomerized caveolin-1 microdomains for epidermal growth factor (EGF) receptor (EGFR) recruitment regulates EGFR signaling in tumor cells. In mammary tumor cells deficient for Golgi beta1,6N-acetylglucosaminyltransferase V (Mgat5), a reduction in EGFR binding to the galectin lattice allows an increased association with stable caveolin-1 cell surface microdomains that suppresses EGFR signaling. Depletion of caveolin-1 enhances EGFR diffusion, responsiveness to EGF, and relieves Mgat5 deficiency-imposed restrictions on tumor cell growth. In Mgat5(+/+) tumor cells, EGFR association with the galectin lattice reduces first-order EGFR diffusion rates and promotes receptor interaction with the actin cytoskeleton. Importantly, EGFR association with the lattice opposes sequestration by caveolin-1, overriding its negative regulation of EGFR diffusion and signaling. Therefore, caveolin-1 is a conditional tumor suppressor whose loss is advantageous when beta1,6GlcNAc-branched N-glycans are below a threshold for optimal galectin lattice formation.     

Control of cell motility by interaction of gangliosides, tetraspanins, and epidermal growth factor receptor in A431 versus KB epidermoid tumor cells

Growth of epidermoid carcinoma cell lines, A431 and KB, has been known to be controlled by the interaction of epidermal growth factor (EGF) and its receptor (EGFR) with tyrosine kinase. Ganglioside GM3 was previously found to interact with EGFR and to inhibit EGFR tyrosine kinase. However, motility of these cells, controlled by EGFR and ganglioside, was not studied. The present study is focused on the control mechanism of the motility of these cells through interaction of ganglioside, tetraspanin (TSP), and EGFR. Key results are as follows: (i) The level of EGFR expressed in A431 cells is 6 times higher than that expressed in KB cells, and motility of A431 cells is also much higher than that of KB cells, yet growth of A431 cells is either not affected or is inhibited by EGF. In contrast, growth of KB cells is enhanced by EGF. (ii) Levels of TSPs (CD9, CD82, and CD81) expressed in A431 cells are much higher than those expressed in KB cells, and TSPs expressed in A431 cells are reduced by treatment of cells with EtDO-P4, which inhibits the synthesis of glycosphingolipids (GSLs) and gangliosides. (iii) These TSPs are co-immunoprecipitated with EGFR in both A431 and KB cells, indicating that TSPs are closely associated with EGFR. (iv) High motility of A431 cells is greatly reduced, while low motility of KB cells is not affected, by treatment of cells with EtDO-P4. These results, taken together, suggest that there is a close correlation between high motility of A431 cells and high expression of EGFR and TSPs, and between ganglioside GM3/GM2 and TSP. A similar correlation was suggested between the low motility of KB cells and low levels of EGFR and TSP. The correlation between high motility and high level of EGFR with the ganglioside–TSP complex in A431 cells is unique. This is in contrast to our previous studies that indicate that motility of many types of tumor cells is inhibited by a high level of CD9 or CD82, together with growth factor receptors and integrins.     

Epidermal growth factor receptor (EGFR) signaling in cancer

The epidermal growth factor receptor (EGFR) belongs to the ErbB family of receptor tyrosine kinases (RTK). These trans-membrane proteins are activated following binding with peptide growth factors of the EGF-family of proteins. Evidence suggests that the EGFR is involved in the pathogenesis and progression of different carcinoma types. The EGFR and EGF-like peptides are often over-expressed in human carcinomas, and in vivo and in vitro studies have shown that these proteins are able to induce cell transformation. Amplification of the EGFR gene and mutations of the EGFR tyrosine kinase domain have been recently demonstrated to occur in carcinoma patients. Interestingly, both these genetic alterations of the EGFR are correlated with high probability to respond to anti-EGFR agents. However, ErbB proteins and their ligands form a complex system in which the interactions occurring between receptors and ligands affect the type and the duration of the intracellular signals that derive from receptor activation. In fact, proteins of the ErbB family form either homo- or hetero-dimers following ligand binding, each dimer showing different affinity for ligands and different signaling properties. In this regard, evidence suggests that cooperation of multiple ErbB receptors and cognate ligands is necessary to induce cell transformation. In particular, the growth and the survival of carcinoma cells appear to be sustained by a network of receptors/ligands of the ErbB family. This phenomenon is also important for therapeutic approaches, since the response to anti-EGFR agents might depend on the total level of expression of ErbB receptors and ligands in tumor cells.     

Effect of epidermal growth factor receptor internalization on regulation of the phospholipase C-gamma1 signaling pathway

The epidermal growth factor receptor (EGFR) ligands, epidermal growth factor (EGF), and transforming growth factor-alpha (TGFalpha) elicit differential postendocytic processing of ligand and receptor molecules, which impacts long-term cell signaling outcomes. These differences arise from the higher affinity of the EGF-EGFR interaction versus that of TGFalpha-EGFR in the acidic conditions of sorting endosomes. To determine whether EGFR occupancy in endosomes might also affect short-term signaling events, we examined activation of the phospholipase C-gamma1 (PLC-gamma1) pathway, an event shown to be essential for growth factor-induced cell motility. We found that EGF continues to stimulate maximal tyrosine phosphorylation of EGFR following internalization, while, as expected, TGFalpha stimulates markedly less. The resulting higher level of receptor activation by EGF, however, did not yield higher levels of phosphatidylinositol (4,5)-bisphosphate (PIP2) hydrolysis over those stimulated by TGFalpha. By altering the ratio of activated receptors between the cell surface and the internalized compartment, we found that only cell surface receptors effectively participate in PLC function. In contrast to PIP2 hydrolysis, PLC-gamma1 tyrosine phosphorylation correlated linearly with the total level of Tyr(P)-EGFR stimulated by either ligand, indicating that the functional deficiency of internal EGFR cannot be attributed to an inability to interact with and phosphorylate signaling proteins. We conclude that EGFR signaling through the PLC pathway is spatially restricted at a point between PLC-gamma1 phosphorylation and PIP2 hydrolysis, perhaps because of limited access of EGFR-bound PLC-gamma1 to its substrate in endocytic trafficking organelles.     

Escape of autocrine ligands into extracellular medium: experimental test of theoretical model predictions

We have developed an experimental system for testing mathematical model predictions concerning escape of autocrine ligands into the extracellular bulk medium. This system employs anti-receptor blocking antibodies against the epidermal growth factor receptor (EGFR)/transforming growth factor alpha (TGFalpha) receptor/ligand pair. TGFalpha was expressed under the control of a tetracycline-repressed promoter, together with a constitutively expressed human EGFR in B82 mouse fibroblast cells. This expression system allowed us to vary TGFalpha synthesis rates over a roughly 300-fold range by adjusting tetracycline concentration. TGFalpha accumulation in the extracellular bulk medium was then measured as a function of cell density, TGFalpha synthesis rate, and anti-EGFR blocking antibody concentration. Consistent with model predictions, amounts of ligand in the medium on a per cell basis were found to diminish as cell density was increased but with reduced dependence on cell density at higher ligand synthesis rates. Similarly consistent with model predictions, higher ligand synthesis rates also decreased the effect of anti-receptor blocking antibodies. Our investigation has established that we can successfully analyze and understand autocrine ligand secretion behavior from the basis of our theoretical model.     

p21-activated kinase 1 (PAK1) interacts with the Grb2 adapter protein to couple to growth factor signaling

A variety of intracellular signaling pathways are linked to cell surface receptor signaling through their recruitment by Src homology 2 (SH2)/SH3-containing adapter molecules. p21-activated kinase 1 (PAK1) is an effector of Rac/Cdc42 GTPases that has been implicated in the regulation of cytoskeletal dynamics, proliferation, and cell survival signaling. In this study, we describe the specific interaction of PAK1 with the Grb2 adapter protein both in vitro and in vivo. We identify the site of this interaction as the second proline-rich SH3 binding domain of PAK1. Stimulation of the epidermal growth factor receptor (EGFR) in HaCaT cells enhances the level of EGFR-associated PAK1 and Grb2, although the PAK1-Grb2 association is itself independent of this stimulation. A cell-permeant TAT-tagged peptide encompassing the second proline-rich SH3 binding domain of PAK1 simultaneously blocked Grb2 and activated EGFR association with PAK1, in vitro and in vivo, indicating that Grb2 mediates the interaction of PAK1 with the activated EGFR. Blockade of this interaction decreased the epidermal growth factor-induced extension of membrane lamellae. Thus Grb2 may serve as an important mechanism for linking downstream PAK signaling to various upstream pathways.