ABSTRACT
The incidence of neoplastic diseases shows that the search for more effective cancer treatments is still necessary. In addition to standard treatments such as chemo- and radiotherapy, new treatment modalities have been focused in recent advances in immunology. Since little has been discussed about the biological implications of chemotherapy with respect to its impact on the immune system and on the other normal, “healthy” cells, we devote the present mathematical modeling work to do so. First, we prove the invariance of the region where all state variables remain positive (i.e., the number of cancer cells, normal cells and immune cells, and also the amount of chemotherapy). Afterward, we analyze the model in terms of the linear stability of the system and establish a necessary and sufficient condition for the local stability of the cure equilibrium point. Moreover, we simulate some scenarios involving both the immune system and chemotherapy, showing that a reasonable treatment strategy occurs when these are combined suitably.
Keywords:
cancer; immune system; chemotherapy; mathematical modeling; ordinary differential equations
RESUMO
A incidência de doenças neoplásicas mostra que a busca por tratamentos mais eficazes contra o câncer ainda é necessária. Além dos tratamentos padrão, como quimioterapia e radioterapia, novas modalidades de tratamento têm sido focadas nos recentes avanços na imunologia. Como pouco tem sido discutido sobre as implicações biológicas da quimioterapia em relação ao seu impacto tanto no sistema imunológico como em outras células normais e “saudáveis”, dedicamos o presente trabalho de modelagem matemática a fazê-lo. Primeiro, provamos a invariância da região na qual todas as variáveis de estado permanecem positivas (ou seja, o número de células cancerosas, células normais e células imunes, e também a quantidade de quimioterapia). Em seguida, analisamos o modelo em termos de estabilidade linear do sistema e estabelecemos uma condição necessária e suficiente para a estabilidade local do ponto de equilíbrio de cura. Além disso, simulamos alguns cenários envolvendo tanto o sistema imunológico quanto a quimioterapia, mostrando que uma estratégia de tratamento razoável ocorre quando estes são combinados adequadamente.
Palavras-chave:
câncer; sistema imune; quimioterapia; modelagem matemática; equações diferenciais ordinárias
1 INTRODUCTION
In the near future, cancer has the potential to become the leading cause of premature death, overshadowing cardiovascular diseases 1515 World Health Organization. Cancer. URL URL http://www.who.int/topics/cancer/en/
. Http://www.who.int/topics/cancer/en/. Accessed: February 25, 2019.
http://www.who.int/topics/cancer/en/...
assical treatments for cancer include chemotherapy, radiotherapy, and surgery. More recently, immunotherapy has been also an alternative to fight against cancer growth. Although there is a consistent understanding of the immunological responses to various pathogens, many aspects of the tumor-immune interaction remain unknown not only from the biological point of view but also quantitatively.
Mathematical models of tumor growth under the influence of the immune system have been addressed in several papers by de Pillis and collaborators (e.g., 44 L.G. de Pillis & A. Radunskaya. A mathematical model of immune response to tumor invasion. In K.J. Bathe (editor), “Computational Fluid and Solid Mechanics 2003”. Elsevier Science Ltd, Oxford (2003), pp. 1661-1668. doi:10.1016/B978-008044046-0.50404-8. URL https://www.sciencedirect.com/science/article/pii/B9780080440460504048.
https://www.sciencedirect.com/science/ar...
),(55 L.G. de Pillis, A.E. Radunskaya & C.L. Wiseman. A validated mathematical model of cell-mediated immune response to tumor growth. Cancer Research, 65(17) (2005), 7950-7958. doi:10.1158/0008-5472.CAN-05-0564.
https://doi.org/10.1158/0008-5472.CAN-05...
e latter, the authors make use of estimated parameters through experimental data to propose different ways of tumor-immune interaction and tumor dynamics models, which are further validated. With respect to tumor-immune interaction with chemotherapy, the search for better strategies of treatment is also addressed by de Pillis et al. in 22 G.L. De Pillis & A. Radunskaya. A mathematical tumour model with immune resistance and drug therapy: An optimal control approach. Journal of Theoretical Medicine, 3(2) (2001), 79-100. doi:10.1080/10273660108833067. URL https://www.hindawi.com/journals/cmmm/2001/318436/abs/.
https://www.hindawi.com/journals/cmmm/20...
),(33 L.G. de Pillis, W. Gu, K.R. Fister, T. Head, K. Maples, A. Murugan, T. Neal & K. Yoshida. Chemotherapy for tumors: An analysis of the dynamics and a study of quadratic and linear optimal controls. Mathematical Biosciences, 209(1) (2007), 292-315. doi:10.1016/j.mbs.2006.05.003. URL http://www.sciencedirect.com/science/article/pii/S0025556406000836.
http://www.sciencedirect.com/science/art...
the authors report tests with different optimal control strategies for chemotherapy.
In the present paper, based on 22 G.L. De Pillis & A. Radunskaya. A mathematical tumour model with immune resistance and drug therapy: An optimal control approach. Journal of Theoretical Medicine, 3(2) (2001), 79-100. doi:10.1080/10273660108833067. URL https://www.hindawi.com/journals/cmmm/2001/318436/abs/.
https://www.hindawi.com/journals/cmmm/20...
), (66 A.C. Fassoni. “Modelos matemáticos em câncer abordando fase inicial e tratamento de tumor avascular (in portuguese)”. Ph.D. thesis, Universidade Estadual de Campinas, Campinas (2016). URL http://repositorio.unicamp.br/jspui/handle/REPOSIP/307199.
http://repositorio.unicamp.br/jspui/hand...
), (88 V.A. Kuznetsov, I.A. Makalkin, M.A. Taylor & A.S. Perelson. Nonlinear dynamics of immunogenic tumors: Parameter estimation and global bifurcation analysis. Bulletin of Mathematical Biology, 56(2) (1994), 295-321. doi:10.1016/S0092-8240(05)80260-5. URL http://www.sciencedirect.com/science/article/pii/S0092824005802605.
http://www.sciencedirect.com/science/art...
e propose a mathematical model that involves the influence of the immune system and chemotherapy on tumor cells, but also the normal cells and the unavoidable effect of chemotherapy on these. First, we show that the region of practical interest is invariant (once state variables are set positive they remain positive). Afterward, we analyze the model in terms of the linear stability of the system and establish a necessary and sufficient condition for the local stability of the cure equilibrium point. We also address some scenarios involving both the immune system and chemotherapy, showing that a reasonable treatment strategy occurs when these are combined suitably.
2 MODEL
Based on 66 A.C. Fassoni. “Modelos matemáticos em câncer abordando fase inicial e tratamento de tumor avascular (in portuguese)”. Ph.D. thesis, Universidade Estadual de Campinas, Campinas (2016). URL http://repositorio.unicamp.br/jspui/handle/REPOSIP/307199.
http://repositorio.unicamp.br/jspui/hand...
), (22 G.L. De Pillis & A. Radunskaya. A mathematical tumour model with immune resistance and drug therapy: An optimal control approach. Journal of Theoretical Medicine, 3(2) (2001), 79-100. doi:10.1080/10273660108833067. URL https://www.hindawi.com/journals/cmmm/2001/318436/abs/.
https://www.hindawi.com/journals/cmmm/20...
and 1111 D.S. Rodrigues, P.F. de Arruda Mancera, T. de Carvalho & L.F. Gonçalves. A mathematical model for chemoimmunotherapy of chronic lymphocytic leukemia. Applied Mathematics and Computation, 349 (2019), 118-133. doi:10.1016/j.amc.2018.12.008. URL https://www.sciencedirect.com/science/article/pii/S009630031831049X.
https://www.sciencedirect.com/science/ar...
, let us consider N
1 the number of tumor cells, N
2 the number of normal cells, I the number of cells of the immune system, and Q the amount of chemotherapeutic drug. Our proposed model is given by:
Although there are different mathematical models that are used to describe tumor growth, we consider the logistic model (see 1414 V.G. Vaidya & F.J. Alexandro Jr. Evaluation of some mathematical models for tumor growth. International Journal of Bio-Medical Computing, 13(1) (1982), 19-35. doi:10.1016/0020-7101(82)90048-4. URL http://www.sciencedirect.com/science/article/pii/0020710182900484.
http://www.sciencedirect.com/science/art...
, where r
1 is the intrinsic growth rate, with the tumor cells carrying capacity given by k
1. The coefficients of competition between the populations N
1 and N
2 are given by α
ij , which measure the effects of population j on the population . The parameter r
2 represents the total constant reproduction of normal cells (see 66 A.C. Fassoni. “Modelos matemáticos em câncer abordando fase inicial e tratamento de tumor avascular (in portuguese)”. Ph.D. thesis, Universidade Estadual de Campinas, Campinas (2016). URL http://repositorio.unicamp.br/jspui/handle/REPOSIP/307199.
http://repositorio.unicamp.br/jspui/hand...
).
The dynamics in the immune cell population is activated by the tumor population at a rate ρ, γ being the half-saturation constant of the Michaelis-Menten functional response given by (see 22 G.L. De Pillis & A. Radunskaya. A mathematical tumour model with immune resistance and drug therapy: An optimal control approach. Journal of Theoretical Medicine, 3(2) (2001), 79-100. doi:10.1080/10273660108833067. URL https://www.hindawi.com/journals/cmmm/2001/318436/abs/.
https://www.hindawi.com/journals/cmmm/20...
), and also there is a natural death rate of immune cells given by m. As we know that some immune cells can directly eliminate tumor cells (see 22 G.L. De Pillis & A. Radunskaya. A mathematical tumour model with immune resistance and drug therapy: An optimal control approach. Journal of Theoretical Medicine, 3(2) (2001), 79-100. doi:10.1080/10273660108833067. URL https://www.hindawi.com/journals/cmmm/2001/318436/abs/.
https://www.hindawi.com/journals/cmmm/20...
two more terms are in order, -c
1
IN
1 and- c
2
IN
1, the latter representing the inactivation of immune cells acting on tumor cells and the former is due to the death of tumor cells due to the action of the immune system. The term s describes a natural source of immune cells (see 88 V.A. Kuznetsov, I.A. Makalkin, M.A. Taylor & A.S. Perelson. Nonlinear dynamics of immunogenic tumors: Parameter estimation and global bifurcation analysis. Bulletin of Mathematical Biology, 56(2) (1994), 295-321. doi:10.1016/S0092-8240(05)80260-5. URL http://www.sciencedirect.com/science/article/pii/S0092824005802605.
http://www.sciencedirect.com/science/art...
).
In order to model the amount of chemotherapy injected in the system, the function models the infusion of the drug into the system and λ is the drug washout rate. As in 1212 D.S. Rodrigues & P.F.A. Mancera. Mathematical analysis and simulations involving chemotherapy and surgery on large human tumours under a suitable cell-kill functional response. Mathematical Biosciences & Engineering, 10 (2013), 221-234. doi:10.3934/mbe.2013.10.221. URL http://aimsciences.org/article/doi/10.3934/mbe.2013.10.221.
http://aimsciences.org/article/doi/10.39...
, the response of each cell population to the chemotherapy is considered to be of a Michaelis-Menten form, with half-saturation parameters a, b and c; μ is the treatment rate of the tumor cells; v is the mortality rate of normal cells due to treatment; and δ represents the mortality rate of immune cells due to the chemotherapeutic drug (see also 1313 D.S. Rodrigues, P.F.A. Mancera & S.T.R. Pinho. Understanding the antiangiogenic effect of metronomic chemotherapy through a simple mathematical model. Physica A: Statistical Mechanics and its Applications, 464 (2016), 251-266. doi:10.1016/j.physa.2016.07.076. URL http://www.sciencedirect.com/science/article/pii/S0378437116305167.
http://www.sciencedirect.com/science/art...
).
Throughout the paper, the drug administration q(t) takes two different forms:
-
1. Continuous administration
-
, and the drug infusion is at a constant rate.
-
2. Administration in cycles
-
Following 99 R. Martin & K.L. Teo. “Optimal Control of Drug Administration in Cancer Chemotherapy”. World Scientific Publishing Company, Singapore (1993). URL https://www.worldscientific.com/worldscibooks/10.1142/2048.
https://www.worldscientific.com/worldsci... , q(t) is a periodic function defined as
-
where T is the time between drug infusions, represents the instants of administration, and τ is the time taken for infusion.
3 RESULTS
We analyze the stability of the system firstly considering that there is no chemotherapy and then we address the chemotherapeutic case.
3.1 Stability analysis without chemotherapy
As the model represents the time evolution of non-negative state variables, it is desirable that these variables remain with this property. But even for a positive initial condition, there is no apparent reason for the solutions not to converge to an equilibrium with one or more negative coordinates. However, from this point on we will prove that this cannot occur given the dynamics on the planes N 1, N 2 and I.
-
On the plane the model is given by:
As the normal cell population is strictly increasing, that is, the trajectories will follow the direction from to . Then, if the initial condition comprises , then the values of , will always be positive (see Figure 1).
-
On the plane the system is written as:
As the immune cells increase, that is, any trajectory will follow the direction from to (see Figure 2).
-
When one gets the following system:
Note that this plane is an invariant subspace, since . As a consequence, given an initial condition in which , the cancer population will always be zero. Moreover, solving the equations on this plane with initial condition (0, N 20, I 0) one gets:
Considering the immune cells variation given by the third equation in (3.1), one sees that the straight line is an invariant set inside the plane . Moreover, since , where , we conclude that H(t) is an attractor, which means that for a large value of the population I converges to s/m and the population N 2 goes to infinity (see Figure 3). These subspaces are very illustrative because they are related to the elimination of cancer cells.
3.2 Stability analysis with chemotherapy
Now we add chemotherapy in our analysis, in order to study the full model (2.1). We have the same behavior in the subspaces and , as the situation without chemotherapy. Let and (constant) results in the following system:
The subspace is still invariant, since to the initial condition one gets
Since is the solution to the cancer cells population, then we write g(t) and f(t) as the solutions to the normal and immune cells populations, respectively. Let us consider the solution of the system of equations given by and W, the subspace where . For an initial condition , the solution is given by
Then, if the solution is in this subspace, it will remain there, resulting that W is invariant. Consequently, we can discard .
The equilibrium point of the chemotherapeutic equation is given by and it is stable, since , where . We have another invariant subspace , given by and . It is invariant because to , the solution given by . When and , we conclude that W 3 has a similar behavior as the invariant plane without chemotherapy, where the solutions are getting closer to the invariant line and . As q increases, , there is an equilibrium point on the invariant line with coordinates
on the subspace with . As N 1 does not change on the invariant subspace W in 3D, we analyze the stability on this subspace. The Jacobian matrix is:
Evaluating the Jacobian matrix on P inf it follows:
The eigenvalues of (3.2) are given by
and
implying that the equilibrium point is stable on this subspace. This point, P inf, occurs due to the killing of normal cells by the chemotherapeutic drug.
Now we study the cancer model in the full domain. The system (2.1) has five equilibrium points for ; four of them have the following structure
and the other is
The point P inf represents the extinction of cancer cells, and it was analyzed in the subspace W, with . The Jacobian matrix is
where
and
Consequently, J 2(P inf) takes the form:
where
and
The eigenvalues of this matrix are
and
The first three eigenvalues are negative since the parameters are positive. The fourth one will be negative under the following condition:
3.3 Numerical simulations: chemotherapy in cycles
To set the chemotherapeutic schedule, we define the function q based on a drug regimen already used in practice. In the specific case of chronic lymphocytic leukemia (CLL), one possible choice is a drug combination with pentostatin, cyclophosphamide and rituximab, usually known as PCR protocol 77 N.E. Kay, S.M. Geyer, T.G. Call, T.D. Shanafelt, C.S. Zent, D.F. Jelinek, R. Tschumper, N.D. Bone, G.W. Dewald, T.S. Lin, N.A. Heerema, L. Smith, M.R. Grever & J.C. Byrd. Combination chemoimmunotherapy with pentostatin, cyclophosphamide, and rituximab shows significant clinical activity with low accompanying toxicity in previously untreated B chronic lymphocytic leukemia. Blood, 109(2) (2007), 405-411. doi:10.1182/blood-2006-07-033274. URL http://www.bloodjournal.org/content/109/2/405.
http://www.bloodjournal.org/content/109/...
. In our numerical simulations, however, we only consider the use of cyclophosphamide 11 F.T. Awan. Cure for CLL? Blood, 127(3) (2016), 274 p. doi:10.1182/blood-2015-11-678532. URL http://www.bloodjournal.org/content/127/3/274.
http://www.bloodjournal.org/content/127/...
.
The chemotherapeutic regimen consists of a 21-day, 5- or 6-cycle schedule of 600 mg of cyclophosphamide per m2 of body surface per infusion. For a hypothetical human patient of height 1.70m and weight 65 kg, resulting in a body surface of 1.75m2 (1010 R.D. Mosteller. Simplified Calculation of Body-Surface Area. New England Journal of Medicine, 317(17) (1987), 1098. doi:10.1056/NEJM198710223171717. URL https://doi.org/10.1056/NEJM198710223171717.
https://doi.org/10.1056/NEJM198710223171...
, the dose per infusion is mg. We consider a drug infusion time of 3 hours or 1/8 day, resulting in mg/day. The other parameters are listed in Table 1.
Then
with for 5 infusions or for 6 infusions.
We have simulated two situations, shown in Figure 4, one with and another with . We observe how important is the value of s (how efficient the immune system is), resulting in the cure or disease relapse.
In Figure 5 we have taken and the number of infusions 5 or 6. We observe that one more infusion results in a change of the equilibrium point attained, from the failure of treatment (5 infusions) to cure (6 infusions).
Figure 6 illustrates the cytotoxic effects of the chemotherapeutic drug on the cells population. We have set new values for the parameters δ and v: and . We note that the cure is not anymore possible and the number of normal cells decreases. The action of the immune system is limited by chemotherapy.
We have increased the value of ρ, which is the rate of production of immune cells stimulated by the cancer cells. Keeping and with 5 or 6 infusions, the number of cancer cells does not go to zero as it can be seen in Figure 7. In this situation, the cancer growth is postponed, but not avoided (see Figure 8).
To show the importance of the immune system, we have fixed ρ as and chosen , with 5 infusions of the drug. The cancer cells go to zero, even with 5 infusions of the chemotherapeutic drug, as seen in Figure 9. This is a crucial situation since we are achieving the cure administrating a lesser amount of chemotherapy.
4 CONCLUSIONS
We developed an ODE mathematical model for tumor growth under the action of the immune system and chemotherapy. We established the positivity of the dependent variables of the model. As a result, we presented the stability analysis of the referred system, revealing the condition that must hold for the cure equilibrium to be attained. From the numerical simulations, we have shown both immune system and chemotherapy are important to the decreasing of tumor population, but a very strong immune system combined with chemotherapy is much more crucial to reducing the number of tumor cells to a low level. According to our results, it is also crucial to allow the immune system to recover from the chemotherapy before considering more eventual applications of oncological therapies. As a final remark, as it is shown in Figures 8 and 9, let us note that the success of chemotherapy critically depends on the intensity in which the tumor cells stimulate the immune cell production. In the near future, this could be explored in detail in some immunotherapeutic strategy.
REFERENCES
-
1F.T. Awan. Cure for CLL? Blood, 127(3) (2016), 274 p. doi:10.1182/blood-2015-11-678532. URL http://www.bloodjournal.org/content/127/3/274
» https://doi.org/10.1182/blood-2015-11-678532» http://www.bloodjournal.org/content/127/3/274 -
2G.L. De Pillis & A. Radunskaya. A mathematical tumour model with immune resistance and drug therapy: An optimal control approach. Journal of Theoretical Medicine, 3(2) (2001), 79-100. doi:10.1080/10273660108833067. URL https://www.hindawi.com/journals/cmmm/2001/318436/abs/
» https://doi.org/10.1080/10273660108833067» https://www.hindawi.com/journals/cmmm/2001/318436/abs/ -
3L.G. de Pillis, W. Gu, K.R. Fister, T. Head, K. Maples, A. Murugan, T. Neal & K. Yoshida. Chemotherapy for tumors: An analysis of the dynamics and a study of quadratic and linear optimal controls. Mathematical Biosciences, 209(1) (2007), 292-315. doi:10.1016/j.mbs.2006.05.003. URL http://www.sciencedirect.com/science/article/pii/S0025556406000836
» https://doi.org/10.1016/j.mbs.2006.05.003» http://www.sciencedirect.com/science/article/pii/S0025556406000836 -
4L.G. de Pillis & A. Radunskaya. A mathematical model of immune response to tumor invasion. In K.J. Bathe (editor), “Computational Fluid and Solid Mechanics 2003”. Elsevier Science Ltd, Oxford (2003), pp. 1661-1668. doi:10.1016/B978-008044046-0.50404-8. URL https://www.sciencedirect.com/science/article/pii/B9780080440460504048
» https://doi.org/10.1016/B978-008044046-0.50404-8» https://www.sciencedirect.com/science/article/pii/B9780080440460504048 -
5L.G. de Pillis, A.E. Radunskaya & C.L. Wiseman. A validated mathematical model of cell-mediated immune response to tumor growth. Cancer Research, 65(17) (2005), 7950-7958. doi:10.1158/0008-5472.CAN-05-0564.
» https://doi.org/10.1158/0008-5472.CAN-05-0564 -
6A.C. Fassoni. “Modelos matemáticos em câncer abordando fase inicial e tratamento de tumor avascular (in portuguese)”. Ph.D. thesis, Universidade Estadual de Campinas, Campinas (2016). URL http://repositorio.unicamp.br/jspui/handle/REPOSIP/307199
» http://repositorio.unicamp.br/jspui/handle/REPOSIP/307199 -
7N.E. Kay, S.M. Geyer, T.G. Call, T.D. Shanafelt, C.S. Zent, D.F. Jelinek, R. Tschumper, N.D. Bone, G.W. Dewald, T.S. Lin, N.A. Heerema, L. Smith, M.R. Grever & J.C. Byrd. Combination chemoimmunotherapy with pentostatin, cyclophosphamide, and rituximab shows significant clinical activity with low accompanying toxicity in previously untreated B chronic lymphocytic leukemia. Blood, 109(2) (2007), 405-411. doi:10.1182/blood-2006-07-033274. URL http://www.bloodjournal.org/content/109/2/405
» https://doi.org/10.1182/blood-2006-07-033274» http://www.bloodjournal.org/content/109/2/405 -
8V.A. Kuznetsov, I.A. Makalkin, M.A. Taylor & A.S. Perelson. Nonlinear dynamics of immunogenic tumors: Parameter estimation and global bifurcation analysis. Bulletin of Mathematical Biology, 56(2) (1994), 295-321. doi:10.1016/S0092-8240(05)80260-5. URL http://www.sciencedirect.com/science/article/pii/S0092824005802605
» https://doi.org/10.1016/S0092-8240(05)80260-5» http://www.sciencedirect.com/science/article/pii/S0092824005802605 -
9R. Martin & K.L. Teo. “Optimal Control of Drug Administration in Cancer Chemotherapy”. World Scientific Publishing Company, Singapore (1993). URL https://www.worldscientific.com/worldscibooks/10.1142/2048
» https://www.worldscientific.com/worldscibooks/10.1142/2048 -
10R.D. Mosteller. Simplified Calculation of Body-Surface Area. New England Journal of Medicine, 317(17) (1987), 1098. doi:10.1056/NEJM198710223171717. URL https://doi.org/10.1056/NEJM198710223171717
» https://doi.org/10.1056/NEJM198710223171717» https://doi.org/10.1056/NEJM198710223171717 -
11D.S. Rodrigues, P.F. de Arruda Mancera, T. de Carvalho & L.F. Gonçalves. A mathematical model for chemoimmunotherapy of chronic lymphocytic leukemia. Applied Mathematics and Computation, 349 (2019), 118-133. doi:10.1016/j.amc.2018.12.008. URL https://www.sciencedirect.com/science/article/pii/S009630031831049X
» https://doi.org/10.1016/j.amc.2018.12.008» https://www.sciencedirect.com/science/article/pii/S009630031831049X -
12D.S. Rodrigues & P.F.A. Mancera. Mathematical analysis and simulations involving chemotherapy and surgery on large human tumours under a suitable cell-kill functional response. Mathematical Biosciences & Engineering, 10 (2013), 221-234. doi:10.3934/mbe.2013.10.221. URL http://aimsciences.org/article/doi/10.3934/mbe.2013.10.221
» https://doi.org/10.3934/mbe.2013.10.221» http://aimsciences.org/article/doi/10.3934/mbe.2013.10.221 -
13D.S. Rodrigues, P.F.A. Mancera & S.T.R. Pinho. Understanding the antiangiogenic effect of metronomic chemotherapy through a simple mathematical model. Physica A: Statistical Mechanics and its Applications, 464 (2016), 251-266. doi:10.1016/j.physa.2016.07.076. URL http://www.sciencedirect.com/science/article/pii/S0378437116305167
» https://doi.org/10.1016/j.physa.2016.07.076» http://www.sciencedirect.com/science/article/pii/S0378437116305167 -
14V.G. Vaidya & F.J. Alexandro Jr. Evaluation of some mathematical models for tumor growth. International Journal of Bio-Medical Computing, 13(1) (1982), 19-35. doi:10.1016/0020-7101(82)90048-4. URL http://www.sciencedirect.com/science/article/pii/0020710182900484
» https://doi.org/10.1016/0020-7101(82)90048-4» http://www.sciencedirect.com/science/article/pii/0020710182900484 -
15World Health Organization. Cancer. URL URL http://www.who.int/topics/cancer/en/ Http://www.who.int/topics/cancer/en/. Accessed: February 25, 2019.
» http://www.who.int/topics/cancer/en/
Publication Dates
-
Publication in this collection
16 Sept 2019 -
Date of issue
May-Aug 2019
History
-
Received
08 Aug 2018 -
Accepted
28 Feb 2019