{"id":2729,"date":"2015-04-28T08:20:50","date_gmt":"2015-04-28T08:20:50","guid":{"rendered":"http:\/\/biomedpharmajournal.org\/?p=2729"},"modified":"2020-04-26T07:29:35","modified_gmt":"2020-04-26T07:29:35","slug":"structurally-steady-control-of-tumor-cell-populations","status":"publish","type":"post","link":"https:\/\/biomedpharmajournal.org\/staging\/vol6no2\/structurally-steady-control-of-tumor-cell-populations\/","title":{"rendered":"Structurally Steady Control of Tumor Cell Populations"},"content":{"rendered":"<p><strong>Introduction<\/strong><\/p>\n<p>Modern medicine is providing reliable approaches for screening, diagnosis, treatment, prevention, and early intervention of cancer. Fortunately, current research levels are high enough to identify many cancer types as curable and studies are underway to find those cures.<\/p>\n<p>In this work, we study mathematical models for tumor growth. Some previous studies have developed mathematical models for biological immune systems. For example, Adimya and Crauste considered a nonlinear mathematical model of hematopoietic stem cell dynamics in which proliferation and apoptosis are controlled by growth factor concentrations [1].<\/p>\n<p>Structured models were found by Roeder et al. [2] to consistently explain a broad variety of phenomena in patients with chronic myeloid leukemia. Furthermore, Kim et al. [3, 4] developed a model for describing the dynamics during the treatment of chronic myelogenous leukemia. This model is intended to replace the recent agent-based model.<\/p>\n<p>By changing the initial susceptible population after each epidemic generation, models developed by Roeder et al. [5] incorporate the effects of vaccination.<\/p>\n<p>In this work, we model tumor growth by considering interactions between tumor cells and cells that resist cancer growth, that is, cells that regulate antitumor activity. The latter can be viewed as a control system applied to tumors, and the creation of models that combine control laws with a class of structurally steady mappings constitutes a new approach to cancer study. The work may be important for both its practical applications and its theoretical novelty.<\/p>\n<p>To develop and test new models in the theory of control in immunology, we apply classical and modern theories of control, theoretical and applied methods in the study of dynamic systems, modern catastrophe theory and its application to population models, and Lyapunov stability theory.<\/p>\n<p>We view cancer as uncontrolled proliferation of tumor cells in an organism. Our task is to suppress outbreaks and stabilize the organism (through the influence of a vaccine) when outbreaks do occur.<\/p>\n<p><strong>Mathematical Models<\/strong><\/p>\n<p>To track the growth of tumor and tumor-resistant cell populations, we start with a relatively simple model (\u00a7 2.1). Then we add terms to account for the control of tumor growth (\u00a7 2.2) and determine the stability of that extended model (\u00a7 2.3). Then we extend the model further to account for the suppression of the boomerang effect (\u00a7 2.4).<\/p>\n<p><strong>Initial Model<\/strong><\/p>\n<p>To model cell populations, including their interactions and influence, we start with two dynamical variables: the concentration of tumor cells and the concentration of cells that have specific resistance to tumor growth. Cells with natural resistance have some level of activity, which is included in the model as a parameter.<\/p>\n<p>In particular, the temporal evolution of the two concentrations is modeled by these two differential equations [6]:<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9017\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for1.jpg\" alt=\"Vol-6No1_STRU_JANA_for1\" width=\"512\" height=\"138\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for1-300x81.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for1.jpg 512w\" sizes=\"(max-width: 512px) 100vw, 512px\" \/><\/p>\n<p>Here <em>x<\/em><sub>1<\/sub> is the concentration of free cells (effector cells), <em>x<\/em><sub>2<\/sub> is the concentration of tumor cells, <em>j<\/em> is the growth rate of effector cells, <em>b<\/em><em>x<\/em><sub>1<\/sub><em>x<\/em><sub>2<\/sub> is the death rate of effector cells due to their interactions with tumor cells, <em>x<\/em><sub>1<\/sub><em>x<\/em><sub>2<\/sub> is the death rate of tumor cells due to their interactions with effector cells, <em>a<\/em> controls the growth rate of effector cells due to interactions with tumor cells, <em>d<\/em> controls the growth rate of tumor cells, <em>g<\/em> controls the natural death rate of effector cells, and <em>m<\/em> controls the natural death rate of tumor cells. Thus, our model is a system of two differential equations involving six parameters.<\/p>\n<p>The phase portrait [7, 8] for system (1) is shown in Fig. 1. The stationary points of the model are determined from the conditions <em>x<\/em><sub>1<\/sub> = 0 and <em>x<\/em><sub>2<\/sub> = 0. Therefore, from (1), model 1 has three equilibrium states:<img decoding=\"async\" class=\"alignnone size-full wp-image-9018\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f1.jpg\" alt=\"Vol-6No1_STRU_JANA_f1\" width=\"392\" height=\"75\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f1-300x57.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f1.jpg 392w\" sizes=\"(max-width: 392px) 100vw, 392px\" \/><\/p>\n<p>These are marked by the filled circles in Fig. 1.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td>\u00a0<img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-9061\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig1-150x150.jpg\" alt=\"Figure 1: Phase Portrait of System (1).\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig1-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig1.jpg 328w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 1:\u00a0Phase Portrait of System (1).<\/strong><\/p>\n<p>&nbsp;<\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig1.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>Note that the Jacobi matrix of system (1) is<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9019\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f2.jpg\" alt=\"Vol-6No1_STRU_JANA_f2\" width=\"371\" height=\"126\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f2-300x102.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f2.jpg 371w\" sizes=\"(max-width: 371px) 100vw, 371px\" \/><\/p>\n<p><strong>Adding a Control Law to the Initial Model<\/strong><\/p>\n<p>To promote the stabilization of an organism in a predictable manner, we present a new method for controlling cell populations through interactions between effector and tumor cells. We call this approach a \u201ccompensatory effect.\u201d<\/p>\n<p>To model system (1), we add a synchronic controlling D-factor, which is a function of structurally steady mappings known as a \u201cfold catastrophe.\u201d One of the most important points in catastrophe theory is the behavior of smooth functions that depend on a parameter. The D-factor used here represents a fold form (one of the 7 forms of catastrophe theory) [9, 10]. Then model (1) becomes<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9020\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f3.jpg\" alt=\"Vol-6No1_STRU_JANA_f3\" width=\"333\" height=\"127\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f3-300x114.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f3.jpg 333w\" sizes=\"(max-width: 333px) 100vw, 333px\" \/><\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9021\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for2.jpg\" alt=\"Vol-6No1_STRU_JANA_for2\" width=\"512\" height=\"135\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for2-300x79.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for2.jpg 512w\" sizes=\"(max-width: 512px) 100vw, 512px\" \/><\/p>\n<p>Negative values indicate a specific influence on the organism. The controlling parameter can be thought of as the influence of the therapy (vaccination) being received.<\/p>\n<p>Stationary points are determined from the conditions , and therefore, system (2) has these equilibrium states:<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9022\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f4.jpg\" alt=\"Vol-6No1_STRU_JANA_f4\" width=\"337\" height=\"192\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f4-300x171.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f4.jpg 337w\" sizes=\"(max-width: 337px) 100vw, 337px\" \/><\/p>\n<p>The phase portrait of system (2) is shown in Fig 2.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td>\u00a0<img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-9062\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig2-150x150.jpg\" alt=\"Figure 2: Phase Portrait of System (2).\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig2-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig2.jpg 317w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 2:\u00a0Phase Portrait of System (2).<\/strong><\/p>\n<p>&nbsp;<\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig2.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>The Jacobi matrix of system (2) is<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9023\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f5.jpg\" alt=\"Vol-6No1_STRU_JANA_f5\" width=\"445\" height=\"114\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f5-300x77.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f5.jpg 445w\" sizes=\"(max-width: 445px) 100vw, 445px\" \/><\/p>\n<p><strong>Region of Stability for the Extended Model<\/strong><\/p>\n<p>We now determine the conditions at which the system of nonlinear equations (2) exhibits stability. To this end, we construct the Lyapunov function. We also use the central Morse theorem with fold catastrophe from the classification of Thom\u2019s theorem.<\/p>\n<p>Stability is a fundamental notion in the qualitative theory of differential equations and is essential for many applications. In turn, Lyapunov functions are a basic instrument for studying stability; however, there is no universal method for constructing Lyapunov functions. Nevertheless, in some special cases, a function can be constructed by applying special techniques [11, 12].<\/p>\n<p>We construct the Lyapunov function for system (2) and then use geometric interpretation to find the region of stability. The geometric identification of stable states reduces to creating a family of closed surfaces that surround the zero equilibrium of coordinates. The system state moves across contour curves; that is, each integrated curve can cross each of these surfaces, as illustrated in Fig. 3.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td>\u00a0<img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-9063\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig3-150x150.jpg\" alt=\"Figure 3: Geometrical Interpretation of Stability\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig3-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig3-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig3.jpg 354w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 3:\u00a0Geometrical Interpretation of Stability<\/strong><\/p>\n<p>&nbsp;<\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig3.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>To apply the geometric interpretation, we must find the Lyapunov function, which is a positive definite state function <em>V(x)<\/em> &gt; 0 whose time derivative is locally negative definite \u00a0&lt; 0 or<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9024\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f6.jpg\" alt=\"Vol-6No1_STRU_JANA_f6\" width=\"194\" height=\"68\" \/><\/p>\n<p>It is known that this system can be described by these dynamical equations:<img decoding=\"async\" class=\"alignnone size-full wp-image-9025\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for3.jpg\" alt=\"Vol-6No1_STRU_JANA_for3\" width=\"579\" height=\"69\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for3-300x36.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for3.jpg 579w\" sizes=\"(max-width: 579px) 100vw, 579px\" \/><\/p>\n<p>It is also known that the time derivative of the corresponding Lyapunov function has the form<img decoding=\"async\" class=\"alignnone size-full wp-image-9026\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for4.jpg\" alt=\"Vol-6No1_STRU_JANA_for4\" width=\"557\" height=\"77\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for4-300x41.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for4.jpg 557w\" sizes=\"(max-width: 557px) 100vw, 557px\" \/><\/p>\n<p>or<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9027\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for5.jpg\" alt=\"Vol-6No1_STRU_JANA_for5\" width=\"575\" height=\"67\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for5-300x35.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for5.jpg 575w\" sizes=\"(max-width: 575px) 100vw, 575px\" \/><\/p>\n<p>The gradient of a function \u00a0is a vector defined by<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9028\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f7.jpg\" alt=\"Vol-6No1_STRU_JANA_f7\" width=\"224\" height=\"76\" \/><\/p>\n<p>We can introduce a vector \u00a0with projections corresponding to (3), namely<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9029\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f8.jpg\" alt=\"Vol-6No1_STRU_JANA_f8\" width=\"291\" height=\"67\" \/><\/p>\n<p>The vector \u00a0is a velocity vector of a point M, as depicted in Fig 4.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td>\u00a0<img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-9064\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig4-150x150.jpg\" alt=\"Figure 4: Gradient of a Function.\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig4-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig4-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig4.jpg 315w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 4:\u00a0Gradient of a Function.<\/strong><\/p>\n<p>&nbsp;<\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig4.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>According to (4) and (5), we have<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9030\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f9.jpg\" alt=\"Vol-6No1_STRU_JANA_f9\" width=\"281\" height=\"70\" \/><\/p>\n<p>where.The time derivative of the Lyapunov function is a scalar product of the gradient and the phase velocity vector. Note that in Fig. 4, the vector is perpendicular to the surface \u00a0and is in the direction of the maximum increase in .<\/p>\n<p>Now, we construct the vector function<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9031\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f10.jpg\" alt=\"Vol-6No1_STRU_JANA_f10\" width=\"320\" height=\"60\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f10-300x56.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f10.jpg 320w\" sizes=\"(max-width: 320px) 100vw, 320px\" \/><\/p>\n<p>whose gradient has a magnitude equal to the absolute value of the velocity vector for system (2) and whose direction is opposite to that of the velocity vector. For system (2), we form components of the gradient-vector Lyapunov function and find<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9032\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f11.jpg\" alt=\"Vol-6No1_STRU_JANA_f11\" width=\"294\" height=\"145\" \/><\/p>\n<p>and<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9033\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f12.jpg\" alt=\"Vol-6No1_STRU_JANA_f12\" width=\"340\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f12-300x132.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f12.jpg 340w\" sizes=\"(max-width: 340px) 100vw, 340px\" \/><\/p>\n<p>Then the full derivative of the Lyapunov function is<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9034\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f13.jpg\" alt=\"Vol-6No1_STRU_JANA_f13\" width=\"330\" height=\"212\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f13-300x193.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f13.jpg 330w\" sizes=\"(max-width: 330px) 100vw, 330px\" \/><\/p>\n<p>This expression shows that <em>dV\/dt<\/em> &lt; 0.<\/p>\n<p>To find the integral components, we integrate<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9035\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f14.jpg\" alt=\"Vol-6No1_STRU_JANA_f14\" width=\"82\" height=\"71\" \/><\/p>\n<p>only over \u00a0and integrate<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9036\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f15.jpg\" alt=\"Vol-6No1_STRU_JANA_f15\" width=\"91\" height=\"67\" \/><\/p>\n<p>only over .\u00a0We also integrate<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9037\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f16.jpg\" alt=\"Vol-6No1_STRU_JANA_f16\" width=\"83\" height=\"73\" \/><\/p>\n<p>and<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9038\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f17.jpg\" alt=\"Vol-6No1_STRU_JANA_f17\" width=\"90\" height=\"73\" \/><\/p>\n<p>first over \u00a0and then over . Therefore, we have [13]\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9039\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f18.jpg\" alt=\"Vol-6No1_STRU_JANA_f18\" width=\"204\" height=\"94\" \/><\/p>\n<p>and<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9040\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f19.jpg\" alt=\"Vol-6No1_STRU_JANA_f19\" width=\"199\" height=\"80\" \/><\/p>\n<p>We find<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9041\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f20.jpg\" alt=\"Vol-6No1_STRU_JANA_f20\" width=\"384\" height=\"190\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f20-300x148.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f20.jpg 384w\" sizes=\"(max-width: 384px) 100vw, 384px\" \/><\/p>\n<p>and<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9043\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f21.jpg\" alt=\"Vol-6No1_STRU_JANA_f21\" width=\"484\" height=\"130\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f21-300x81.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f21.jpg 484w\" sizes=\"(max-width: 484px) 100vw, 484px\" \/><\/p>\n<p>Thus, the Lyapunov function takes this scalar form:<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9044\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f22.jpg\" alt=\"Vol-6No1_STRU_JANA_f22\" width=\"222\" height=\"41\" \/><\/p>\n<p>However, the Lyapunov function <em>V(x)<\/em> is not positive definite because it contains negative components (terms of odd degree). Hence, we need to find the conditions under which the Lyapunov function is positive definite.<\/p>\n<p>To do so, we next build the Hessian matrix of the Lyapunov function:<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9045\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f23.jpg\" alt=\"Vol-6No1_STRU_JANA_f23\" width=\"419\" height=\"138\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f23-300x99.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f23.jpg 419w\" sizes=\"(max-width: 419px) 100vw, 419px\" \/><\/p>\n<p>and<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9046\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f241.jpg\" alt=\"Vol-6No1_STRU_JANA_f24\" width=\"491\" height=\"140\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f241-300x86.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f241.jpg 491w\" sizes=\"(max-width: 491px) 100vw, 491px\" \/><\/p>\n<p>Thus, the Hessian matrix is<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9047\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f25.jpg\" alt=\"Vol-6No1_STRU_JANA_f25\" width=\"356\" height=\"92\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f25-300x78.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f25.jpg 356w\" sizes=\"(max-width: 356px) 100vw, 356px\" \/><\/p>\n<p>At point (0;0), the Lyapunov function in square form is<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9048\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f26.jpg\" alt=\"Vol-6No1_STRU_JANA_f26\" width=\"329\" height=\"52\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f26-300x47.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f26.jpg 329w\" sizes=\"(max-width: 329px) 100vw, 329px\" \/><\/p>\n<p>and the conditions for being positive definite are given by<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9049\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f27.jpg\" alt=\"Vol-6No1_STRU_JANA_f27\" width=\"157\" height=\"79\" \/><\/p>\n<p>At another stationary point (<em>M;0<\/em>), the matrix Hessian is<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9050\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f28.jpg\" alt=\"Vol-6No1_STRU_JANA_f28\" width=\"409\" height=\"104\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f28-300x76.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f28.jpg 409w\" sizes=\"(max-width: 409px) 100vw, 409px\" \/><\/p>\n<p>and at this point, the Lyapunov function is<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9051\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f29.jpg\" alt=\"Vol-6No1_STRU_JANA_f29\" width=\"348\" height=\"82\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f29-300x71.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f29.jpg 348w\" sizes=\"(max-width: 348px) 100vw, 348px\" \/><\/p>\n<p>Then the conditions for positive definiteness of the square form are<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9052\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f30.jpg\" alt=\"Vol-6No1_STRU_JANA_f30\" width=\"230\" height=\"108\" \/><\/p>\n<p>From the condition \u00a0and the definition of <em>M<\/em> below (2), the parameters satisfy \u00a0and . They are also real.<\/p>\n<p>Thus, we have found conditions for the parameters under which our Lyapunov function is positive definite and its time derivative is negative definite. These conditions identify an area of robust stability for dynamical system (2) [14]. This means that system (2) captures system states in which tumors do not continue to grow.<\/p>\n<p><strong>Boomerang<\/strong> <strong>Effect<\/strong><\/p>\n<p>Although the model in (2) includes situations in which tumors do not grow, in clinical practice, we know that many tumors initially respond well to treatment, but after a short time, they start to grow again. That is, once a chemical treatment has started, some tumors may adapt; hence, the treatment is no longer effective. We call this the \u201cboomerang effect\u201d [7]. Clinically, protective measures must be repeated regularly to avoid the boomerang effect. If protective measures are not regular, there will be time for tumor cells to adapt to chemical treatments.<\/p>\n<p>To exclude the \u201cboomerang effect\u201d in our mathematical model, we consider the following model:<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9053\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for6.jpg\" alt=\"Vol-6No1_STRU_JANA_for6\" width=\"454\" height=\"109\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for6-300x72.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for6.jpg 454w\" sizes=\"(max-width: 454px) 100vw, 454px\" \/><\/p>\n<p>where the fractional terms describe loss changes. Now we add D-fa\u0441tor terms to system (6):<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9054\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for7.jpg\" alt=\"Vol-6No1_STRU_JANA_for7\" width=\"646\" height=\"119\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for7-300x55.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_for7.jpg 646w\" sizes=\"(max-width: 646px) 100vw, 646px\" \/><\/p>\n<p>Thus, with the D-factors, we exclude the boomerang effect and can control the growth of cell numbers resulting from uncontrolled cell proliferation. Phase portraits for system (7) appear in Fig. 7.<\/p>\n<p><strong>Results and Discussion<\/strong><strong>\u00a0<\/strong><\/p>\n<p>In this section, we illustrate the behavior of each of the above models for selected values of their parameters. Figure 5 shows results from system (1) when its six parameters have these values:<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9055\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f31.jpg\" alt=\"Vol-6No1_STRU_JANA_f31\" width=\"382\" height=\"43\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f31-300x34.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f31.jpg 382w\" sizes=\"(max-width: 382px) 100vw, 382px\" \/><\/p>\n<p>The stationary point at <em>x<\/em><sub>2<\/sub> = 0, <em>x<\/em><sub>1<\/sub> \u2260 0 represents zero concentration of cancer cells, that is, complete recovery of the organism.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td>\u00a0<img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-9065\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig5-150x150.jpg\" alt=\"Figure 5: Phase Portrait of System (1).\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig5-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig5.jpg 235w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 5:\u00a0Phase Portrait of System (1).<\/strong><\/p>\n<p>&nbsp;<\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig5.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>Figure 6 shows calculated results for system (2) using parameter values<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9056\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f32.jpg\" alt=\"Vol-6No1_STRU_JANA_f32\" width=\"441\" height=\"36\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f32-300x24.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f32.jpg 441w\" sizes=\"(max-width: 441px) 100vw, 441px\" \/><\/p>\n<p>The values for the first six parameters are the same as those used for system (1) above; the difference is only in the seventh parameter. Figure 6 also shows a stationary state that represents a complete cure (<em>x<\/em><sub>2<\/sub> = 0).<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td>\u00a0<img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-9066\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig6-150x150.jpg\" alt=\"Figure 6: Phase Portrait of System (2).\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig6-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig6.jpg 229w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 6:\u00a0Phase Portrait of System (2).<\/strong><\/p>\n<p>&nbsp;<\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig6.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>Results for system (6) are also shown in Fig. 5 for the following parameter values:<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9058\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f33.jpg\" alt=\"Vol-6No1_STRU_JANA_f33\" width=\"435\" height=\"59\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f33-300x41.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f33.jpg 435w\" sizes=\"(max-width: 435px) 100vw, 435px\" \/><\/p>\n<p>Figure 7 contains results for system (7). Figure 7a pertains to these parameter values:<\/p>\n<p><img decoding=\"async\" class=\"alignnone size-full wp-image-9059\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f34.jpg\" alt=\"Vol-6No1_STRU_JANA_f34\" width=\"431\" height=\"63\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f34-300x44.jpg 300w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_f34.jpg 431w\" sizes=\"(max-width: 431px) 100vw, 431px\" \/><\/p>\n<p>In Fig. 7b, the results are for the same parameter values, except that the value of parameter <em>J<\/em> was increased from 10 to 20.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td>\u00a0<img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-9067\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig7-150x150.jpg\" alt=\"Figure 7: Phase Portrait of System (7).\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig7-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig7.jpg 501w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 7:\u00a0Phase Portrait of System (7).<\/strong><\/p>\n<p>&nbsp;<\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2015\/04\/Vol-6No1_STRU_JANA_fig7.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>In both Figs. 7a and 7b, a complete cure is possible, but in Fig. 7b, the cure requires a higher concentration of effector cells (response of the immune system).<\/p>\n<p><strong>Conclusions<\/strong><\/p>\n<p>This paper has reported on mathematical methods for modeling and controlling growth of tumor cells in living cancerous tissues. The models are pairs of nonlinear differential equations that include terms to account for growth of tumor and effector cells, control of cancerous growth by effector cells, natural mortality of both kinds of cells, and suppression of the boomerang effect in cancer treatments. Control of tumor growth and exclusion of the boomerang effect were modeled using one-parameter mappings (D-factors).<\/p>\n<p>In general, phase portraits from the models exhibit one or two stationary states (see Fig. 2). The first represents cancer treatment that reduces the concentration of tumor cells to zero: complete recovery of the organism. The second stationary state is a \u201cfrozen\u201d state that represents treatment in which the tumor is not eliminated, but its growth can be controlled at the expense of effector cells (the immune system).<\/p>\n<p>In the case of model system (2), we could construct the Lyapunov function and identify the region of stability in terms of model parameters. In <em>physical terms<\/em>, this means that the parameter controlling tumor growth somewhat decreases the population of effector cells. In <em>biological <\/em><em>terms<\/em>, it means that for the system to be stable, the growth of effector cells () must be positive, i.e., there should be enough protein in the cells. Furthermore, interactions between parameters that represent natural mortality rates of the tumor and effective cells (and) are essential. Therefore, effector cells have to be chemo-stable.<\/p>\n<p>We also demonstrated that by introducing a D-factor into models, we can exclude the boomerang effect. The boomerang effect occurs in some cancer treatments in which tumor growth is initially reduced, but the tumor adapts to the treatment, and the tumor resumes growth.<\/p>\n<p>More generally, this paper illustrates that the \u201ccompensatory effect\u201d can be applied to control the growth of tumor cells in an organism [15, 16], fixing the immune system at a predictable level [17, 18]. Furthermore, with the help of regulating factors, it is possible to evaluate the level to which populations of tumor cells can be controlled [19, 20].<\/p>\n<p>These results should motivate new theoretical studies on the medical-biological treatment of tumor cells and on the regulation of the immune system. These studies are especially important for the prevention and early intervention of cancer.<\/p>\n<p><strong>Acknowledgments<\/strong><\/p>\n<p>Special thanks Professor Khlebopros R.G., Professor Beisenbi M.A, As. Prof. Nikulin V.V., As. Prof, Ph.D. Omarov A.N. for fruitful discussions on the biological and mathematical implications and possible extensions of the work.<\/p>\n<p><strong>References<\/strong><\/p>\n<ol>\n<li>Adimya, M. and Crauste, F. Mathematical model of hematopoiesis dynamics with growth factor-dependent apoptosis and proliferation regulations. Comput. Model. <u>49<\/u> (2009) 2128-2137.<\/li>\n<li>Roeder, I., Herberg, M. and Horn, M. An \u201cage\u201d-structured model of hematopoietic stem cell organization with application to chronic myeloid leukemia. B. Math. Biol. <u>71<\/u>(3) (2009) 602-626.<\/li>\n<li>Kim, P.S., Lee, P.P. and Levy, D. Modeling imatinib-treated chronic myelogenous leukemia: reducing the complexity of agent-based models. B. Math. Biol. <u>70<\/u>(3) (2008) 728-744<\/li>\n<li>Kim, P., Lee, P. and Levy, D. A PDE model for imatinib-treated chronic myelogenous leukemia. B. Math. Biol. <u>70<\/u> (7) (2008) 1994-2016.<\/li>\n<li>Roeder, I., Herberg, M. and Horn, M. Modelling repeated epidemics with general infection kernels. Res. Lett. Inf. Math. Sci. <u>6<\/u> (2004) 159-170.<\/li>\n<li>Bazykyn, A.D. Nonlinear Dynamics of Interaction of Population, <em>Moscow<\/em> 2003.<\/li>\n<li>Khlebopros, R.G., Isaev, A.S. and Nedorezov, L.V. Dynamics of wood insects number. <em>Novosibirsk. Sciene,<\/em> 1984.<\/li>\n<li>Khlebopros, R.G., Okhonin, V.A. and Fet, A.I. Catastrophes in Nature and Society, Volume 1, 2010.<\/li>\n<li>Nikolis, G. and Prigogine, I. Exploring complexity: an introduction, <em>Moscow<\/em>, 1990.<\/li>\n<li>Poston, T. and Stewart, I. Catastrophe Theory and its Applications, <em>Dover Publications<\/em>, 1996.<\/li>\n<li>Boris, T.P. Introduction to optimization. Optimization Software, <em>Inc. Publications Division, New York<\/em>, 2010.<\/li>\n<li>Liapunov Functions and Stability in Control Theory. Bacciotti, Andrea, Rosier, Lionel. Originally published as volume 267 in the series \u201cLecture Notes Control\u201d by <em>Springer, London,<\/em> 2001, 2nd ed., 2005, XIV, 238 p. 20 illus.<\/li>\n<li>Schwartz, L. Analyse Mathematique. <em>Hermann,<\/em> Vol. 1. Part IV. \u21168, 1967.<\/li>\n<li>Beisenbi, M. A., Kulniyazova, K. S. Determining the region of robust stability of control systems via Lyapunov functions, in: The Proceedings of the VIth International Asian School-Seminar. Siberian Branch Russian Academic Science Novosibirsk, 2010.<\/li>\n<li>Yermekbayeva, \u00a0\u00a0 J.J. and Omarov, A.N. Compensatory effect for controlling the population of phytophages and enthomophages, Eurasian Entomological Journal <u>9<\/u>(3) (2010). <em>http:\/\/www.euroentjourn.narod.ru\/indexrus.html<\/em><\/li>\n<li>Yermekbayeva, J.J. The doctoral thesis. The analysis and research of complex behaviour population models. ISBN 9965-31-334-4, L.N.Gumilyov Eurasian National University, <em>Reviewer:<\/em> <em>R.G.Khlebopros (2010).<\/em><\/li>\n<li>Yermekbayeva, J.J. Use of one-parametrical structural-steady steady mapping for control of system predator-prey. <em>Bulletin of L.N.Gumilyov Eurasian National University, 2010, Vol.4.<\/em> <em>http:\/\/repository.enu.kz\/xmlui\/bitstream\/handle\/123456789\/83\/ispolzovanie_odnoparametricheskih.pdf?sequence=1<\/em>.<\/li>\n<li>Yermekbayeva, J.J. The Influence of structural and steady displays in competition models. <em>The proceedings of the Perspektywiczne opracowania s\u0105 nauk\u0105 i technikami \u2013 2011: materia\u0142y vii i\u0119dzynarodowej naukowi-praktycznej konferencji<\/em>, Przemy\u015bl, 2011.-V.11. S.15-20.\u00a0<em>http:\/\/www.rusnauka.com\/28_PRNT_2011\/Economics\/8_94027.doc.htm<\/em><\/li>\n<li>Yermekbayeva, J.J. Control of tumoral cages growth with the help of \u00abcompensatory effect\u00bb. <em>The proceedings of the XVI All-Russia symposium with international participation \u00abComplex systems in extreme conditions\u00bb,<\/em> Krasnoyarsk scientific center of the Siberian Branch of the Russian Academy of Science, 2012.<\/li>\n<li>Abitova, G., Nikulin, V., Beisenbi, M. and Ainagulova, A. Control System with High Robust Stability Characteristics Based On Catastrophe Function. The proceedings of the ICECCS\u2019 12. 2012.<\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"<p>Introduction Modern medicine is providing reliable approaches for screening, diagnosis,  [&#8230;]<\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[17],"tags":[],"class_list":["post-2729","post","type-post","status-publish","format-standard","hentry","category-vol6no2"],"_links":{"self":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/2729","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/comments?post=2729"}],"version-history":[{"count":5,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/2729\/revisions"}],"predecessor-version":[{"id":33227,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/2729\/revisions\/33227"}],"wp:attachment":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/media?parent=2729"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/categories?post=2729"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/tags?post=2729"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}