Ligand‐dependent responses of the ErbB signaling network: experimental and modeling analyses

Open Access
Publisher Website
Google Scholar
Abstract
Deregulation of ErbB signaling plays a key role in the progression of multiple human cancers. To help understand ErbB signaling quantitatively, in this work we combine traditional experiments with computational modeling, building a model that describes how stimulation of all four ErbB receptors with epidermal growth factor (EGF) and heregulin (HRG) leads to activation of two critical downstream proteins, extracellular‐signal‐regulated kinase (ERK) and Akt. Model analysis and experimental validation show that (i) ErbB2 overexpression, which occurs in approximately 25% of all breast cancers, transforms transient EGF‐induced signaling into sustained signaling, (ii) HRG‐induced ERK activity is much more robust to the ERK cascade inhibitor U0126 than EGF‐induced ERK activity, and (iii) phosphoinositol‐3 kinase is a major regulator of post‐peak but not pre‐peak EGF‐induced ERK activity. Sensitivity analysis leads to the hypothesis that ERK activation is robust to parameter perturbation at high ligand doses, while Akt activation is not. ### Synopsis The ErbB receptors, ErbB1 (EGFR), ErbB2 (HER2/NEU), ErbB3, and ErbB4, are widely expressed throughout the human body and control cell fate decisions, such as proliferation, apoptosis, and differentiation ([Yarden and Sliwkowski, 2001][1]). Given the wide physiological expression of ErbB receptors and their importance in controlling cell fate decisions, it is not surprising that deregulation of ErbB signaling is implicated in the progression of multiple human cancers ([Yarden and Sliwkowski, 2001][1]). Because the ErbB signaling network is a complex, dynamic system with multiple feedbacks, it is difficult to understand its overall signaling behavior through traditional, qualitative methods. In the current work, we combine traditional experimental methods with quantitative modeling to understand ErbB signaling in MCF‐7 breast cancer cells. We specifically focus on the short‐term (⩽30 min) dynamic behavior of all four ErbB receptors and the key downstream intermediates ERK and Akt in response to stimulation with the ligands epidermal growth factor (EGF) and heregulin (HRG) ([Figure 1][2]). Based on the schematic in [Figure 1][2], we build an ordinary differential equation model consisting of 117 species, 235 parameters, and 96 net reactions to represent the dynamics of ligand‐dependent ErbB network activation. We use this model to generate hypotheses about how the cell dynamically controls ERK and Akt activation in a ligand‐dependent fashion, and experimentally corroborate the predictions that EGF‐induced ERK activity is much more sensitive than HRG‐induced ERK activity to the MEK inhibitor U0126 ([Figure 7A and B][3]), and that the phosphoinositol‐3 kinase (PI‐3K) inhibitor wortmannin has a large affect on post‐ but not pre‐peak EGF‐induced ERK signaling. Gaining insight into the behavior of such a large model for the purposes of hypothesis generation is a difficult task, and in this work, we relied mainly on parametric sensitivity analysis to accomplish this. Naturally, the sensitivity analysis results revealed many obvious modes of regulation, such as control of ERK activation by MEK, and control of Akt activation by PI‐3K. However, these results also revealed non‐obvious, ligand‐dependent modes of regulation: (i) ERK activation is robust to parametric perturbations at high ligand doses, while Akt activation is not. This difference in robustness to parameter variation may be due to the involvement of ERK in multiple negative feedback loops, and the absence of such feedback from Akt. It is well known from control theory that negative feedback endows a system with robustness to disturbances ([Ogunnaike and Ray, 1994][4]; [Freeman, 2000][5]). (ii) PI‐3K abundance is a dominant factor for the post‐peak (declining portion of the ERK response curve) but not pre‐peak (rising portion of the ERK response curve) EGF‐induced ERK activation. This phenomenon is due to degradation of active ErbB1 homodimers. At short times, there are still many 1‐1 homodimers that can signal to ERK; as time progresses, however, these homodimers are degraded. Therefore, at longer times, signaling to ERK must rely on alternative mechanisms—in this case, the PI‐3K‐Grb2‐associated binder 1 (Gab1) pathway. The results of the sensitivity analysis allowed us to gain general understanding of the ligand‐dependent control of ERK and Akt activity. While this general understanding is valuable by itself, understanding how known, cancer‐correlated network abnormalities affect signaling is of primary and immediate interest. We used the model to investigate how such a network abnormality, ErbB2 overexpression, which occurs in approximately 25% of breast cancer, would affect the dynamics of ERK and Akt activation. Simulation results suggested that it should have a large effect on signaling at high EGF doses by transforming normally transient signals into more sustained signals. Significantly, this model prediction quantitatively agrees with the data of Wolf‐Yadlin et al (2006) and Kumar et al (2007), which shows that at 10 and 30 min after EGF stimulation, ERK activation is between 1.15‐ and 2‐fold higher for human mammary epithelial cells overexpressing ErbB2. This effect of ErbB2 overexpression is a direct consequence of ligand‐induced ErbB1 homodimer internalization and degradation. EGF transmits signals only through ErbB1 homodimers and ErbB1‐ErbB2 heterodimers, and ErbB1 homodimers undergo preferential ligand‐induced internalization and degradation relative to other ErbB dimers ([Baulida et al , 1996][6]). As ErbB2 overexpression shifts the ErbB dimer distribution toward more ErbB1‐ErbB2 heterodimers rather than ErbB1 homodimers, it...