- Research
- Open access
- Published:

# A mathematical model separates quantitatively the cytostatic and cytotoxic effects of a HER2 tyrosine kinase inhibitor

*Theoretical Biology and Medical Modelling*
**volume 4**, Article number: 14 (2007)

## Abstract

### Background

Oncogene signaling is known to deregulate cell proliferation resulting in uncontrolled growth and cellular transformation. Gene amplification and/or somatic mutations of the HER2/Neu (ErbB2) proto-oncogene occur in approximately 20% of breast cancers. A therapeutic strategy that has been used to block HER2 function is the small molecule tyrosine kinase inhibitor lapatinib. Using human mammary epithelial cells that overexpress HER2, we determined the anti-proliferative effect of lapatinib through measuring the total cell number and analyzing the cell cycle distribution. A mathematical model was used to interpret the experimental data.

### Results

The model suggests that lapatinib acts as expected by slowing the transition through G_{1} phase. However, the experimental data indicated a previously unreported late cytotoxic effect, which was incorporated into the model. Both effects depend on the dosage of the drug, which shows saturation kinetics.

### Conclusion

The model separates quantitatively the cytostatic and cytotoxic effects of lapatinib and may have implications for preclinical studies with other anti-oncogene therapies.

## Background

Molecule-targeted anti-cancer drugs have been developed as a result of our understanding of tumor cell and molecular biology. Compared to "traditional" cancer therapies, targeted drugs such as the tyrosine kinase inhibitors (TKIs) have higher specificity and relatively lower toxicity in selected patients with corresponding oncogene expression. Members of the type 1 receptor tyrosine kinase (RTK) family, which includes the epidermal growth factor receptor (EGFR), HER2 (ErbB2), HER3 and HER4 play a crucial role in growth and differentiation of both normal and malignant mammary epithelial cells [1, 2]. Binding of receptor-specific ligands to the ectodomain of EGFR, HER3 and HER4 results in the formation of receptor dimers and hetero-oligomers to which HER2 is recruited as the preferred heterodimerization partner [3]. HER2 gene amplification has been reported in approximately 20% of breast cancers, where it is associated with poor patient outcome [4]. Studies with HER2-overexpressing breast cancer cell lines and human tumors have shown constitutive phosphorylation of HER2 [5, 6]. Overexpression of HER2 is associated with transformation of mammary epithelial cells [7, 8] as well as shorter survival in patients with breast carcinoma [4, 9]. These facts make HER2 a rational therapeutic target in human breast cancer. One therapeutic approach against HER2-overexpressing breast cancers is the generation of trastuzumab, a humanized IgG1 that binds to residues 529–626 in domain IV of the HER2 ectodomain [2, 10]. However, many patients with HER2-overexpressing advanced disease do not respond clinically to trastuzumab and many that initially respond eventually relapse with antibody-resistant disease. Lapatinib (GW572016, GlaxoSmithKline) is a selective reversible inhibitor of both EGFR and HER2 tyrosine kinases. Lapatinib mimics ATP and binds to the ATP site in the tyrosine kinase domain of HER2, resulting in blockade of the receptor's catalytic activity [11].

Preclinical data have shown that tumor cells overexpressing EGFR or HER2 are growth inhibited by lapatinib both in vitro and in vivo [12–14]. Lapatinib inhibits the activation of cell proliferation effectors, Erk1/2 (also known as mitogen-activated protein kinase, or MAPK) and AKT, which are downstream of EGFR and HER2 [11, 15]. In another study in which over 30 breast cancer cell lines were tested for their responses to lapatinib, concentration-dependent antiproliferative effects of lapatinib were seen in all cells but varied significantly between individual cell lines [13]. Response to lapatinib is significantly correlated with HER2 expression and its ability to inhibit the phosphorylation of HER2 and downstream effectors. In phase II clinical trials, treatment with lapatinib resulted in objective tumor responses in 28% of patients with HER2-positive advanced breast cancer [12]. Modeling the antiproliferative effects of this oncogene inhibitor using mathematical tools will lead to novel insights into the functioning features and mechanisms of the inhibitor. The model may also provide constructive clinical implications, such as the predictive effects of the inhibitor in first-line therapy in combination with chemotherapy.

In this study we used MCF10A human mammary epithelial cells engineered to overexpress HER2 in order to determine the anti-tumor effects of lapatinib. Compared to control MCF10A cells that do not overexpress HER2, MCF10A/HER2 cells exhibit a gain-of-function phenotype including increased proliferation and filling of the lumen when grown in three dimensions, as a result of oncogene overexpression [16]. Lapatinib inhibits the phosphorylation and function of HER2 in these cells and suppresses growth [16]. At the molecular level the functional mechanisms of HER2 inhibitors are evaluated by the activities of downstream signaling networks, which are often determined by immunoblots. However, signaling pathways such as the PI3K/Akt and the MEK/Erk pathways can converge at various levels of the signaling cascades, making it difficult to separate a combined effect on cell growth and survival. Quantitative models can separate the strengths of drug action on individual phases of the cell cycle. Previous molecular biological studies have shown that HER2 is associated with increases of both G_{1}-S-specific cyclins (cyclins D and E) and G_{2}-M-specific cyclin (cyclin A) [17, 18], which are crucial for G_{1}-S and G_{2}-M progression, respectively. Our objective in this study is to use quantitative models to determine if HER2 inhibitors abolish the function on both phase transitions and how this contributes to cell cycle blockage.

Mathematical modeling has been applied extensively to study the growth kinetics of tumors, with and without treatment; see [19–26] and the references therein. These authors have focused on phenomena such as decelerated growth, quiescence, homeostasis and chemotherapy scheduling.

It has been recognized that, apart from killing cells outright, anticancer drugs can also act by delaying the progression through the cell cycle. Moreover, this blocking effect can be phase specific [27]. Transition through one phase of the cell cycle may be delayed while transition through another phase is unaffected. Mathematical modeling here provides the tool to test possible alternative scenarios against each other and to gain new insight. In a series of papers, Ubezio and collaborators used a mathematical modeling approach to investigate phase-specific cytotoxic and cytostatic effects of drugs such as cisplatin, melphalan and topotecan *in vitro* [26–29]. A continuous model has been used by Agur and coworkers [22] to predict the effect of periodic treatments with cycle-specific cytotoxic drugs.

Our mathematical model consists of populations of proliferating and nonproliferating cells with individual cells distinguished by cell cycle position and is described in detail below. Numerical simulations of the model give good agreement with the experimental data. We find that the experimental data are consistent with a theory in which lapatinib preferentially affects cells growing in monolayer culture in G_{1}-phase in a dose-specific manner. As the dose of lapatinib is increased, however, our study indicates that other phases of the cell cycle are affected as well. Moreover, we see a gradual onset of the cytostatic effect as opposed to a sudden onset. We observe a simple functional relationship between the strength of the cytostatic effect and the drug concentration (see equation (13) for details). Finally, our study indicates that a cytotoxic effect is present after longer periods of exposure to the drug.

## Results

In the control scenario (without treatment) the cell counts showed an initial exponential increase of the population and then a leveling off (see Figure 1A). To explain this leveling off, the nonproliferating cell class was incorporated into the model. Nonlinear models with nonproliferating subpopulations have been used extensively to explain Gompertzian growth kinetics of tumors [20, 21]. Proliferating cells enter the nonproliferating class irreversibly at a rate dependent on their maturity and the total population count of both proliferating and nonproliferating cells. This nonlinearity in the model accounts for the confluence observed in the control study on day 6. Staining of cells with the marker for proliferation Ki-67 showed a dramatic decrease of the proliferating fraction from day 4 to day 6 (see Figure 2B), well borne out by the numerical simulation. In the model, nonproliferating cells arrested their maturity value at the moment of transition from proliferation (see equations (2)–(3)). No mortality of cells was assumed in the model for the control, since no decrease in cell numbers was observed. In addition, staining for the marker of apoptosis Caspase-3 was negative for the control (data not shown). The flow cytometric data from the control scenario are shown in Figure 1B. Owing to uncertainty in the experimental measurements, there were some discrepancies in the fit, particularly during days 1 – 3.

The model for the control case was used as a reference for the treatment cases, with two separate effects of the drug added. The first was the cytostatic effect, which slowed maturation velocity. Our numerical simulations indicate that lapatinib preferentially blocks cells in G_{1}-phase. At higher dose (2 *μM*) the model also incorporates blocking effects in G_{2}/M phase. We find that the strength of the cytostatic effect saturates at higher doses (see Table 1 and Figure 3A). The second effect of the drug was a cytotoxic action. This was incorporated into the model to explain the decrease in cell counts from day 5 to day 6, which was not present in the control (Figure 4). In the model it was assumed that this cytotoxic action only set in after 5 days. This is supported by experimental observations, as staining for the marker of apoptosis Caspase-3 was negative before day 5 (Figure 5). The model simulations agreed substantially with the experimental data, both in the total population counts and the flow cytometric data (Figure 1C, D).

## Discussion and conclusion

The mathematical model provided a means to separate the cytostatic and cytotoxic action of the drug in the experiments. We summarize our findings as follows.

### Cell cycle specificity of cytostatic effects

The strength of the cytostatic effects depends on the drug dosage; the drug showed saturation kinetics (see the last section and Figure 3A for details). At low concentrations (0.1 – 1 *μM*) only cells in G_{1}-phase are delayed. At higher concentrations (2 *μM*) we hypothesize, on the basis of our model simulations, that cells in G_{2}-phase are delayed as well and may be prevented from entering mitosis. It was not necessary to introduce a cytostatic effect for cells in S-phase. We therefore suggest that cells in S-phase remain unaffected at all concentrations.

### Dynamic behavior

Our numerical simulations indicate a buildup phase for the drug action that stretches over several days. Initially we assumed a sudden onset of the cytostatic action after a certain time. The simulations showed pronounced oscillations in the fractions of cells in each phase that were not observed in the data (simulations not shown). We then assumed a gradual increase (with respect to time) of the cytostatic effects. As the cycle time of an average cell is less than one day, we conjecture that the initial effects of the drugs alters a cell's protein contents but still allows division. Only after several generations is the progression of the cells through the respective cell cycle phases fully retarded. The same holds for the loss of cells. This can be explained by the functional mechanism of the drug and the nature of cell mitosis. Oncogene inhibitors such as lapatinib affect the activities of oncogene downstream effectors, which usually include regulatory proteins crucial for proliferation and survival. As cytoplasmic division occurs during mitosis, the inhibitory effect of lapatinib on these crucial effectors can be "inherited" as protein concentration in the cytosol, where the drug effect can further accumulate. Once the concentrations of these crucial proteins exceed certain thresholds, physiological effects such as growth arrest and apoptosis will be induced in the descendant cells. The conclusions that the length of lapatinib treatment is crucial for the overall drug effects, and that it takes several cell generations for the drug to show a clear cellular effect, especially at a lower concentration, may have potential clinical implications.

### The cytotoxic effect

A model with cytostatic action alone cannot lead to a decrease in total cell number, which made it necessary to introduce a cytotoxic effect. To explain the additional mortality of cells after 120 hours, we propose that cells in which progression through the cell cycle has been retarded for too long become prone to apoptosis in the presence of lapatinib. We conjecture that the strength of the cytotoxic effect also saturates at higher doses, as is suggested by Figure 3B.

While it is difficult to extrapolate conclusions from our *in vitro* study to the *in vivo* situation, the following suggestions are plausible. Lapatinib acts chiefly through slowing the progression of proliferating cells in monolayer culture. Furthermore, it is advisable to combine lapatinib with cytotoxic therapeutic agents that kill not only proliferating cells but also quiescent cells, such as some alkylating agents. These drugs may complement the antitumor effect of lapatinib and therefore serve as good candidates to be tested in combined treatment in the future.

We have shown that a mathematical model based on population dynamics can be applied to interpret the cytostatic and cytotoxic effects of lapatinib. Earlier mathematical models [27] used a discrete partition of the cell cycle into age compartments and also a discrete time scale. We find that continuum models are advantageous from the viewpoint of parametrization and computability. In particular, the number of free parameters in our model is significantly smaller than in previously proposed models and each has a straightforward biological interpretation. Both the cytostatic and cytotoxic effects of lapatinib are currently being investigated in various other HER2-overexpressing breast cancer cell lines. Our model can certainly be applied to other oncogene inhibitors that have cytostatic effects on cells during a specific phase of the cell cycle.

## Materials and methods

MCF10A/HER2 cells were generated and grown as described [30]. Lapatinib was kindly provided by Tona Gilmer (GlaxoSmithKline, Research Triangle Park, NC). Equal numbers of cells (3·10^{4}/well) were seeded on 6-well plates for growth assay and cell cycle analysis over a time course of 6 days. At 4 *h* post seeding when > 98% of cells were attached to the plates, the cells were either harvested by trypsinization as day 0 samples, or treated by different doses (0.1, 0.5, 1 and 2 *μM*) of lapatinib, which was added and left in the medium throughout. Media containing fresh drug were replenished every 2 days. Cells were harvested by trypsinization every 8 *h* after adding drug for the initial 3 days and every 24 *h* for days 4–6. One tenth of the cells collected (100 *μl* out of 1 *ml*) were subjected to total cell number counting using a Coulter counter. The rest of the cells (900 *μl* out of 1 *ml*) were subjected to cell cycle analysis. For flow cytometric analysis of cell cycle distribution, cells were fixed in 70% ethanol for 24 *h* at – 20°C, and rehydrated in cold phosphate buffered saline (PBS) for 30 *min* on ice. Redydrated cells were subsequently labeled with 50 *μg/ml* propidium iodide (Sigma, St. Louis, MO) containing 125 units/*ml* protease-free RNase (Calbiochem, San Diego, CA) in the dark for 30 *min* at room temperature and filtered through a 95 *μm* pore size nylon mesh (Small Parts, Miami Lakes, FL). A total of 10,000 stained nuclei was analyzed in a FACS/Calibur Flow Cytometer (BD Biosciences, Franklin Lakes, NJ). All treatments were done in triplicate.

Indirect immunofluorescence assays (IFA) were performed as described previously [16] to detect markers for proliferation and apoptosis. Fluorescent images were captured using a Zeiss inverted LSM510 confocal microscopy system. Primary antibodies include Ki-67 (Calbiochem) and cleaved caspase-3 (Cell Signaling, Danvers, MA). The fluorescent antibodies were Oregon Green-*α*-mouse IgG and Texas Red-*α*-rabbit IgG (Molecular Probes, Carlsbad, CA).

The mathematical model was designed to quantify the cytostatic and cytotoxic effects of the drug on the basis of the population dynamics observed in the experiments. The model consists of a system of differential equations describing these dynamics over the 6 day time course (see the last section for a detailed description). Cells are classified as proliferating or nonproliferating. In the model, proliferating cells are tracked according to position in the cell cycle by assigning to each cell a variable called maturity. Maturity in the control scenario (0 *μM* lapatinib) corresponds to cell age (*i.e*., time that has passed since the last mitosis). The maturity values in the control delimiting the phases of the cell cycle are set at 0 – 7 *h* (G_{1}), 7 – 11 *h* (S), and 11 – 30 *h* (G_{2}/M) (see Figure 6B and the last section for further explanation). The model takes into account the variability of intermitotic times with mean age of division approximately 19 *h* in the control (see Figure 6A). The mathematical model was programmed using MATLAB (version 7.1, The MathWorks, Inc., Natick, MA). The codes are available upon request from the corresponding author. A standard upwind scheme is used for the numerical solution of the partial differential equations [31].

### The mathematical model

Let *t* ≥ 0 denote the time since the beginning of the experiment and *a* ∈ [0, *a*_{
M
}] denote the maturity of a cell. This maturity variable can be thought of as the position of a cell in its cell cycle. We wish to emphasize that in the absence of cytostatic effects of drugs, maturity coincides with chronological age, the time since cell division. In experimental terms, maturity is measured by differential DNA-content. See the discussion in [26] for further information.

Let *p*(*a*, *t*) and *n*(*a*, *t*) denote the densities of proliferating and nonproliferating cells, respectively, of maturity *a* at time *t*. The total number of cells at time *t* is obtained by integrating both densities over the age space

We state our model equations, which balance the biological processes occurring in time

*p*(*a* 0) = *p*_{0}(*a*)*n*(*a*, 0) = *n*_{0}(*a*).

The left hand side of equation (2) describes the aging process for proliferating cells. On the right hand side of the same equation we find that cells of maturity *a* are lost due to three independent processes. Firstly, cells undergo mitosis, at a rate *β* depending on *a*. Such a cell shows up as two proliferating cells of maturity 0, hence the boundary condition (4). Secondly, proliferating cells are also lost due to a transition into the nonproliferating class. In contrast to proliferating cells, the nonproliferating cells do not mature and do not give rise to new cells. The rate at which the transition between the two classes occurs, \tilde{\mu}, depends on the maturity *a* of the respective cell as well as on the total number of cells *M*. We set

The function *μ*(*a*) is depicted in Figure 6B and its parameters are given in Table 1. The particular choice of a piecewise constant function is a result of the experimental observations for the control scenario (Figure 1**B**). Indeed we saw that the percentages of cells in specific phases change after day 4, as more and more cells enter the nonproliferating class. We find it plausible that a cell that has entered S-phase will finish it and therefore be less prone to entering nonproliferation. Thirdly, there is an additional time-dependent cytotoxic effect *ε* for both classes in the presence of drug. We assume the log-kill hypothesis, i.e. the cell kill is proportional to the instantaneous population [32]. Both equations are supplied with initial maturity distributions *p*_{0} and *n*_{0} at time 0.

The cytostatic action of the drug changes the maturation velocity 1 – *δ* of the proliferating cells. The ordinary differential equation

describes the characteristic curves of equation (2). Since maturation is irreversible, the function *δ* must satisfy 0 ≤ *δ* (*a*, *t*) ≤ 1. In the absence of cytostatic effects, we have *δ* = 0. Then *a* - *t* = *const*, that is cells age one-to-one with time, as stated earlier. On the other hand, if *δ* ≈ 1, cell maturation is (almost) completely blocked. Observe from equation (8) that *δ* is a dimensionless quantity.

The model predicts the total cell number (1) as well as the percentages of cells in the three stages of the cell cycle. (Cells in G_{2}- and M-phases are counted together.) These are defined as

where the boundaries between compartments {a}_{{G}_{1}} and *a*_{
S
}are chosen on the basis of experimental observations and *a*_{
M
}is the maximum age (see Table 1 and Figure 6B). More precisely, the control scenario data were used to fix the parameters {a}_{{G}_{1}} and *a*_{
S
}. We compared the fractions of cells in each phase over time, determined by analysis of DNA content, to the model simulation from equations (9). Following [33], we have chosen the distribution of intermitotic times *φ* to be a specific shifted Γ-distribution

Here *a*_{
m
}is the minimal maturity a cell has to reach before it can divide. The parameter *σ* determines the standard deviation of *ϕ* (see Table 1 and Figure 6A). The corresponding age-dependent proliferation rate is given by

where

is the fraction of cells that reach age *a* without division (see Figure 6A). We allow for a certain percentage to reach maximum age *a*_{
M
}without division. As argued by [34], the mean duration of the cell cycle in solid tumors is relatively constant and therefore not influenced by the total cell number *M*.

Another source of uncertainty lies in the initial age distributions *p*_{0} and *q*_{0} at time 0. The experimental data indicated that initially proliferating cells in all stages of the cell cycle were present while there were no nonproliferating cells. It should be remarked that in the absence of nonlinear crowding effects one would observe asynchronous exponential growth. That is, the total cell number would grow exponentially with a certain well-defined rate *λ* while the percentage of cells in each age bracket would approach a steady state. This phenomenon has been widely studied in the population dynamics literature; we refer here to [35] and the references therein.

How do we let the drug act on the cells? We want to test the hypothesis that only cells in G_{1}-phase are blocked while cells in other phases remain unaffected. Hence, we let

where *d* is the drug concentration, {\delta}_{{G}_{1}} (*d*) ≤ 1 is a dose-dependent rate corresponding to maximum blocking effect and *T =* 144 *h* is the duration of the experiment. We decided to let the action increase linearly with respect to time throughout the entire duration of the experiment. When we tried the simplest way to model time dependence, namely a sudden switch from 0 to a constant {\delta}_{{G}_{1}} > 0, the model predicted oscillations in the percentages that were not present in the experimental data. At the higher concentration of 2 *μM* it becomes necessary to introduce an action of the drug on cells in G_{2}/M-phase, similar to equation (11). In order to explain the decrease of cell counts from day 5 to day 6 we assume that the additional mortality is of the type*ε*(*t*) = max{0, *ε*_{0}(*d*)·(*t* - 120)}.

We list the numerical values of the parameters that are fixed throughout all scenarios in Table 1. The values of the parameters {\delta}_{{G}_{1}} (*d*) and *ε*_{0}(*d*) that differ for specific doses are stated in Table 2. We find that the dependence of {\delta}_{{G}_{1}} (*d*) on the drug dosage *d* is described very well by the equation

with *c*_{1} = 3.5 (Figure 3**A**). A similar behavior can be surmised for the dependence of the cytotoxic effect *ε*_{0}(*d*), which follows

with *c*_{2} = 4.7·10^{-3} and *c*_{3} = 2.1 (Figure 3**B**). The functional relationships (13) and (14) are usable for any dose within a certain range.

## References

Hynes NE, Stern DF: The biology of erbB-2/neu/HER-2 and its role in cancer. Biochim Biophys Acta. 1994, 1198: 165-184.

Roskoski R: The ErbB/HER receptor protein-tyrosine kinases and cancer. Biochem Biophys Res Commun. 2004, 319: 1-11. 10.1016/j.bbrc.2004.04.150.

Yarden Y, Sliwkowski MX: Untangling the ErbB signalling network. Nat Rev Mol Cell Biol. 2001, 2: 127-137. 10.1038/35052073.

Slamon DJ, Godolphin W, Jones LA, Holt JA, Wong SG, Keith DE, Levin WJ, Stuart S, Udove J, Ullrich A: Studies of the HER-2/neu proto-oncogene in human breast and ovarian cancer. Science. 1989, 244: 707-712. 10.1126/science.2470152.

Alimandi M, Romano A, Curia MC, Muraro R, Fedi P, Aaronson SA, Fiore PPD, Kraus MH: Cooperative signaling of ErbB3 and ErbB2 in neoplastic transformation and human mammary carcinomas. Oncogene. 1995, 10: 1813-1821.

Thor AD, Liu S, Edgerton S, Moore D, Kasowitz KM, Benz CC, Stern DF, DiGiovanna MP: Activation (tyrosine phosphorylation) of ErbB-2 (HER-2/neu): a study of incidence and correlation with outcome in breast cancer. J Clin Oncol. 2000, 18: 3230-3239.

Pierce JH, Arnstein P, DiMarco E, Artrip J, Kraus MH, Lonardo F, Fiore PPD, Aaronson SA: Oncogenic potential of erbB-2 in human mammary epithelial cells. Oncogene. 1991, 6: 1189-1194.

Muthuswamy SK, Li D, Lelievre S, Bissell MJ, Brugge JS: ErbB2, but not ErbB1, reinitiates proliferation and induces luminal repopulation in epithelial acini. Nat Cell Biol. 2001, 3: 785-792. 10.1038/ncb0901-785.

Ross JS, Fletcher JA: The HER-2/neu oncogene in breast cancer: prognostic factor, predictive factor, and target for therapy. Oncologist. 1998, 3: 237-252.

Carter P, Presta L, Gorman CM, Ridgway JB, Henner D, Wong WL, Rowland AM, Kotts C, Carver ME, Shepard HM: Humanization of an anti-p185HER2 antibody for human cancer therapy. Proc Natl Acad Sci USA. 1992, 89: 4285-4289. 10.1073/pnas.89.10.4285.

Spector NL, Xia W, Burris H, Hurwitz H, Dees EC, Dowlati A, O'Neil B, Overmoyer B, Marcom PK, Blackwell KL, Smith DA, Koch KM, Stead A, Mangum S, Ellis MJ, Liu L, Man AK, Bremer TM, Harris J, Bacus S: Study of the biologic effects of lapatinib, a reversible inhibitor of ErbB1 and ErbB2 tyrosine kinases, on tumor growth and survival pathways in patients with advanced malignancies. J Clin Oncol. 2005, 23: 2502-2512. 10.1200/JCO.2005.12.157.

Johnston SR, Leary A: Lapatinib: a novel EGFR/HER2 tyrosine kinase inhibitor for cancer. Drugs Today(Barc). 2006, 42: 441-453. 10.1358/dot.2006.42.7.985637.

Konecny GE, Pegram MD, Venkatesan N, Finn R, Yang G, Rahmeh M, Untch M, Rusnak DW, Spehar G, Mullin RJ, Keith BR, Gilmer TM, Berger M, Podratz KC, Slamon DJ: of the dual kinase inhibitor lapatinib (GW572016) against HER-2-overexpressing and trastuzumab-treated breast cancer cells. Cancer Res. 2006, 66: 1630-1639. 10.1158/0008-5472.CAN-05-1182.

Nelson MH, Dolder CR: Lapatinib: a novel dual tyrosine kinase inhibitor with activity in solid tumors. Ann Pharmacother. 2006, 40: 261-269. 10.1345/aph.1G387.

Xia W, Gerard CM, Liu L, Baudson NM, Ory TL, Spector NL: Combining lapatinib (GW572016), a small molecule inhibitor of ErbB1 and ErbB2 tyrosine kinases, with therapeutic anti ErbB2 antibodies enhances apoptosis of ErbB2-overexpressing breast cancer cells. Oncogene. 2005, 24: 6213-6221. 10.1038/sj.onc.1208774.

Wang S, Narasanna A, Perez-Torres M, Xiang B, Wu F, Yang S, Carpenter G, Gazdar A, Muthuswamy S, Arteaga C: HER2 kinase domain mutation results in constitutive phosphorylation and activation of HER2 and EGFR and resistance to EGFR tyrosine kinase inhibitors. Cancer Cell. 2006, 10: 25-38. 10.1016/j.ccr.2006.05.023.

Harari D, Yarden Y: Molecular mechanisms underlying ErbB2/HER2 action in breast cancer. Oncogene. 2000, 19: 6102-6114. 10.1038/sj.onc.1203973.

Le XF, Betrosian I, Mao W, Murray M, Lu Z, Keyomarsi K, Lee MH, Zhao J, Bast RG: Anti-HER2 antibody trastuzumab inhibits CDK2-mediated NPAT and histone H4 expression via the PI3K pathway. Cell Cycle. 2006, 5: 1654-1661.

Bertuzzi A, D'Onofrio A, Fasano A, Gandolfi A: Regression and regrowth of tumour cords following single-dose anticancer treatment. Bull Math Biol. 2003, 65: 903-931. 10.1016/S0092-8240(03)00054-5.

Gyllenberg M, Webb GF: Quiescence as an explanation of Gompertzian tumor growth. Growth Develop Aging. 1989, 53: 25-33.

Gyllenberg M, Webb GF: A nonlinear structured population model of tumor growth with quiescence. J Math Biol. 1990, 28: 671-694. 10.1007/BF00160231.

Kheifetz Y, Kogan Y, Agur Z: Long-range predictability in models of cell populations subjected to phase specific drugs: growth rate approximation using properties of positive compact operators. 2006, [Preprint],http://calvino.polito.it/~mcrtn/library.html

Kimmel M, Świerniak A: Using control theory to make cancer chemotherapy beneficial from phase dependence and resistant to drug resistance. Arch Control Sci. 2004, 14: 105-145.

Magni P, Simeoni M, Poggesi I, Rocchetti M, Nicolao GD: A mathematical model to study the effects of drugs administration on tumor growth dynamics. Math Biosci. 2006, 200: 127-151. 10.1016/j.mbs.2005.12.028.

Spinelli L, Torricelli A, Ubezio P, Basse B: Modelling the balance between quiescence and cell death in normal and tumour cell populations. Math Biosci. 2006, 202: 349-370. 10.1016/j.mbs.2006.03.016.

Ubezio P: Unraveling the complexity of cell cycle effects of anticancer drugs in cell populations. Discrete Contin Dyn Syst Ser B. 2004, 4: 323-335.

Montalenti F, Sena G, Cappella P, Ubezio P: Simulating cancer-cell kinetics after drug treatment: application to cisplatin on ovarian carcinoma. Phys Rev E. 1998, 57: 5877-5887. 10.1103/PhysRevE.57.5877.

Lupi M, Matera G, Branduardi D, D'Incalci M, Ubezio P: Cytostatic and cytotoxic effects of topotecan decoded by a novel mathematical simulation approach. Cancer Res. 2004, 64: 2825-2832. 10.1158/0008-5472.CAN-03-3810.

Lupi M, Cappella P, Matera G, Natoli C, Ubezio P: Interpreting cell cycle effects of drugs: the case of melphalan. Cancer Chemother Pharmacol. 2006, 57: 443-457. 10.1007/s00280-005-0044-1.

Ueda Y, Wang S, Dumont N, Yi JY, Koh Y, Arteaga CL: Overexpression of HER2 (erbB2) in human breast epithelial cells unmasks transforming growth factor

*β*-induced cell motility. J Biol Chem. 2004, 279: 24505-24513. 10.1074/jbc.M400081200.Quarteroni A, Sacco R, Saleri F: Numerical Mathematics. 2000, New York: Springer Verlag

Skipper HE, Schnabel FM, Wilcox WS: Experimental evaluation of potential anticancer agents. Cancer Chemother Rep. 1964, 35: 1-111.

Dibrov BF, Zhabotinsky AM, Neyfakh YA, Orlova MP, Churikova LI: Optimal scheduling for cell synchronization by cycle-phase-specific blockers. Math Biosci. 1983, 66: 167-185. 10.1016/0025-5564(83)90087-1.

Tubiana M: The kinetics of tumour cell proliferation and radiotherapy. Br J Radiol. 1971, 44: 325-347.

Dyson J, Villella-Bressan R, Webb GF: Asynchronous exponential growth in an age structured population of proliferating and quiescent cells. Math Biosci. 2002, 177/178: 73-83. 10.1016/S0025-5564(01)00097-9.

## Acknowledgements

This work was supported by the NCI Integrative Cancer Biology Program (U54-CA113007) and the Vanderbilt Integrative Cancer Biology Center (VICBC). Additional support was provided by R01-CA80195 and by Vanderbilt Breast Cancer SPORE grant P50-CA98131 (CLA). We would like to acknowledge fruitful discussions with Drs. Lourdes Estrada and Vito Quaranta.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Competing interests

The author(s) declare that they have no competing interests.

### Authors' contributions

PH and SEW jointly conceived the study. SEW carried out the experiments. PH created the mathematical model, performed the analysis of the experimental data and drafted the manuscript. CLA participated in the design and coordination of the study. GFW helped to create the mathematical model and participated in writing the manuscript. All authors read and approved the final manuscript.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

**Open Access**
This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License (
https://creativecommons.org/licenses/by/2.0
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Hinow, P., Wang, S.E., Arteaga, C.L. *et al.* A mathematical model separates quantitatively the cytostatic and cytotoxic effects of a HER2 tyrosine kinase inhibitor.
*Theor Biol Med Model* **4**, 14 (2007). https://doi.org/10.1186/1742-4682-4-14

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1742-4682-4-14