Analysis of a mathematical model for the growth of tumors under the action of external inhibitors

 In this paper we make rigorous analysis to a mathematical model for the growth of nonnecrotic tumors under the action of external inhibitors. By external inhibitor we mean an inhibitor that is either developed from the immune system of the body or administered by medical treatment to distinguish with that secreted by tumor itself. The model modifies a similar model proposed by H. M. Byrne and M. A. J. Chaplain. After simply establishing the well-posedness of the model, we discuss the asymptotic behavior of its solutions by rigorous analysis. The result shows that an evolutionary tumor will finally disappear, or converge to a stationary state (dormant state), or expand unboundedly, depending on which of the four disjoint regions Δ ₁ , ..., Δ ₄ the parameter vector (A ₁ ,A ₂ ) belongs to, how large the scaled apoptosis number ˜σ is, and how large the initial radius R ₀ of the tumor is. Finally, we discuss some biological implications of the result, which reveals how a tumor varies when inhibitor supply is increased and nutrient supply is reduced. Received: 6 June 2000 / Revised version: 7 November 2001 / Published online: 8 May 2002     

Analysis of a mathematical model for the growth of tumors

 In this paper we study a recently proposed model for the growth of a nonnecrotic, vascularized tumor. The model is in the form of a free-boundary problem whereby the tumor grows (or shrinks) due to cell proliferation or death according to the level of a diffusing nutrient concentration. The tumor is assumed to be spherically symmetric, and its boundary is an unknown function r=s(t). We concentrate on the case where at the boundary of the tumor the birth rate of cells exceeds their death rate, a necessary condition for the existence of a unique stationary solution with radius r=R ₀ (which depends on the various parameters of the problem). Denoting by c the quotient of the diffusion time scale to the tumor doubling time scale, so that c is small, we rigorously prove that (i) lim inf _t→∞ s(t)>0, i.e. once engendered, tumors persist in time. Indeed, we further show that (ii) If c is sufficiently small then s(t)→R ₀ exponentially fast as t→∞, i.e. the steady state solution is globally asymptotically stable. Further, (iii) If c is not “sufficiently small” but is smaller than some constant γ determined explicitly by the parameters of the problem, then lim sup _t→∞ s(t)<∞; if however c is “somewhat” larger than γ then generally s(t) does not remain bounded and, in fact, s(t)→∞ exponentially fast as t→∞. Received: 25 February 1998 / Revised version: 30 April 1998     

A mathematical model of glioma growth: the effect of chemotherapy on spatio-temporal growth

During the past two decades computerized tomography (CT) and magnetic resonance imaging (MRI) have permitted the detection of tumours at much earlier stages in their development than was previously possible. In spite of this earlier diagnosis the effects of earlier and more extensive treatments have been difficult to document. This failure has led to an increasing awareness of the importance of infiltration of glioma cells into surrounding grossly normal brain tissue such that recurrence still occurs. In this paper a simple mathematical model for the proliferation and infiltration of such tumours is introduced, based in part on quantitative image analysis of histological sections of a human brain glioma and especially on cross-sectional area/volume measurements of serial CT images while the patient was undergoing chemotherapy. The model parameters were estimated using optimization techniques to give the best fit of the simulated tumour area to the CT scan data. Numerical solution of the model on a two-dimensional domain, which took into account the geometry of the brain and its natural barriers to diffusion, was used to determine the effect of chemotherapy on the spatiotemporal growth of the tumour.     

A mathematical study of a two-regional population growth model

The paper provides a mathematical study of a model of urban dynamics, adjusting to an ecological model proposed by Lotka and Volterra. The model is a system of two first-order non-linear ordinary differential equations. The study proposed here completes the original proof by using the main tools such as a Lyapunov function.     

A mathematical model of glioma growth: the effect of extent of surgical resection

Abstract. We have developed a mathematical model based on proliferation and infiltration of neoplastic cells that allows predictions to be made concerning the life expectancies following various extents of surgical resection of gliomas of all grades of malignancy. The key model parameters are the growth rate and the diffusion rate. These rates were initially derived from analysis of a case of recurrent anaplastic astrocytoma treated by chemotherapies. Numerical simulations allow us to estimate what would have happened to that patient if various extents of surgical resection, rather than chemotherapies, had been used. In each case, the shell of the infiltrating tumour that remains after 'gross total removal' or even a maximal excision continues to grow and regenerates the tumour mass remarkably rapidly. By developing a model that allows the growth and diffusion rates to define the distribution of cells at the time of diagnosis, and then varying these rates by about 50%, we created a hypothetical tumour patient population whose survival times show good agreement with the results recently reported by Kreth for treatments of glioblastomas. Tenfold decreases in the rates of growth and diffusion mimic the results reported by many other investigators with more slowly growing gliomas. Thus, the model quantitatively supports the ideas that (i) gliomas infiltrate so diffusely that they cannot be cured by resection alone, surgical or radiological, no matter how extensive that may be; (ii) the more extensive the resection, regardless of the degree of malignancy of the glioma, the greater the life expectancy; and (iii) measurements of the two rates, growth and diffusion, may be able to predict survival rates better than the current histological estimates of the type and grade of gliomas.     

A study of the effect of radiation on megakaryocytopoiesis using a mathematical model

Mathematical modelling used in analysing the postirradiation changes in megakaryocytopoiesis permitted to determine the level of radiation-induced injury in each experiment conducted and to show that megakaryocytopoiesis regulation followed the same mechanism after irradiation as it does normally and after the effect of hydroxyurea and anti-thrombocyte serum. The analysis has demonstrated that after the stem cell death induced by ionizing radiation, the regeneration can be provided by the committed cells, and the level of regeneration is determined by the maturity of precursors.     

The effect of interfering substances on the disinfection process: a mathematical model

AIMS: To gain a greater understanding of the effect of interfering substances on the efficacy of disinfection. METHODS AND RESULTS: Current kinetic disinfection models were augmented by a term designed to quantify the deleterious effect of soils such as milk on the disinfection process of suspended organisms. The model was based on the assumption that inactivation by added soil occurred at a much faster rate than microbial inactivation. The new model, the fat-soil model, was also able to quantify the effect of changing the initial inoculum size (1 x 10(7)-5 x 10(7) ml(-1) of Staphylococcus aureus) on the outcome of the suspension tests. Addition of catalase to the disinfection of Escherichia coli by hydrogen peroxide, resulted in changes to the shape of the log survivor/time plots. These changes were modelled on the basis of changing biocide concentration commensurate with microbial inactivation. CONCLUSIONS: The reduction in efficacy of a disinfectant in the presence of an interfering substance can be quantified through the use of adaptations to current disinfection models. SIGNIFICANCE AND IMPACT OF THE STUDY: Understanding the effect of soil on disinfection efficacy allows us to understand the limitations of disinfectants and disinfection procedures. It also gives us a mechanism with which to investigate the soil tolerance of new biocides and formulations.     

A mathematical model of tumour-induced capillary growth

The corneal limbal vessels of an animal host respond to the presence of a source of Tumour Angiogenesis Factor (TAF) implanted in the cornea by the formation of new capillaries which grow towards the source. This neovasculature can be easily seen and studied and this paper describes a mathematical model of some of the important features of the growth. The model includes the diffusion of TAF, the formation of sprouts from pre-existing vessels and models the movement of these sprouts to form new capillaries as a chemotactic response to the presence of TAF. Numerical results are produced for various values of the parameters which characterize the model and it is suggested that the model might form the framework for further theoretical work on related phenomena such as wound healing or to develop strategies for the investigation of anti-angiogenesis.     

A simplified mathematical model of tumor growth

A one-dimensional model of tumor tissue growth is presented in which the source of mitotic inhibitor is nonuniformly distributed within the tissue (in contrast to many earlier models). As a result, stable and unstable regimes of growth become significantly modified from the uniform-source case, indicating that the model, schematic though it is, is very sensitive to the type of source term assumed, and this has implications for experimental and theoretical comparisons in more realistic geometries.     

Describing force-induced bone growth and adaptation by a mathematical model

This work proposes a mathematical model that qualitative describes the process of mechanically force-induced bone growth and adaptation. The mathematical model includes osteocytes as the key interfacing layer connecting tissue, cellular and molecular signaling levels. Specifically, in the presence of an increase in the mechanical stimuli, osteocytes respond by mechano-transduction releasing the local factors nitric oxide (NO) and prostaglandin E(2) (PGE(2)). These local factors act as the signaling recruitment signals for bone cells progenitors and influence the coupling activity among osteoblasts and osteoclasts during the process of bone remodeling. The model is in agreement with qualitative observations found in the literature concerning the process of bone adaptation and the cellular interactions during a local bone remodeling cycle induced by mechanical stimulation.     

Analysis of a mathematical model for transport of I-albumin

The hypotheses and results given are motivated by the study of the distribution of albumin in man which represents a class of delay-differential systems. The approach used is to study the behavior of the solutions of nonlinear delay-differential systems with variable coefficients under the assumptions of continuity and boundedness of coefficients. The criterions are conditions on the roots of a certain “quasi-polynomial”, i.e., a polynomial in a variable and exponential of that variable. These criterions bear a resemblance to the ones in the case of constant coefficients and retardations and are applicable to this case also. The method is based on Lyapunov type functional with appropriate properties.     

A mathematical model for the growth of bacterial microcolonies on marine sediment

Counts of bacterial microcolonies attached to deep-sea sediment particles showed 4-, 8-, 16-, and 32-celled microcolonies to be very rare. This was investigated with a mathematical model in which microcolonies grew from single cells at a constant growth rate (), detached from particles at constant rate (), and reattached as single cells. Terms for attachment of foreign bacteria (a) and death of single cells (d) were also included. The best method of fitting the model to the microcolony counts was a weighted least-squares approach by which(0.83 hour^–1) was estimated to be about 20 times greater than(0.038 hour^–1). This showed that the bacteria were very mobile between sediment particles and this mobility was explained in terms of attachment by reversible sorption. The implications of the results for the frequency of dividing cell method for estimating growth rates of sediment bacteria are discussed. The ratio of and was found to be very robust both in terms of the errors associated with the microcolony counts and the range of microcolony sizes used to obtain the solution.     

A mathematical model of the dynamics of odontogenic cyst growth

OBJECTIVE: To formulate a mathematical model of odontogenic cyst growth and establish the dynamics of cyst enlargement and role of osmotic pressure forces throughout its growth. STUDY DESIGN: The model assumed a spherical cyst with a semipermeable lining of living cells and a core consisting of degraded cellular material, including generic osmotic material, fed by the continuous death of epithelial cells in the lining. The lining cells were assumed to have both elastic and viscous properties, reflecting the action of physical stresses by the surrounding cyst capsule, composed of fibroblasts and collagen fibers. The model couples the cyst radius and osmotic pressure differences resulting in a system of 2 nonlinear ordinary differential equations. RESULTS: The model predicts that in all parameter regimens the long-time behavior of the cyst is the same and that linear radial expansion results. CONCLUSION: In the early and intermediate stages of cystic growth, osmotic pressure differences play an important role; however, in very large cysts, this role becomes negligible, and cell birth in the lining dominates growth.     

A mathematical model for the growth of mycelial pellet populations

In liquid culture, filamentous organisms often grow in the form of pellets. Growth result in an increase in radius, whereas shear forces result in release of hyphal fragments which act as centers for further pellet growth and development. A previously published model for pellet growth of filamentous microorganisms has been examined and is found to be unstable for certain parameter values. This instability has been identified as being due to inaccuracies in estimating the numbers of fragments which seed the pellet population. A revised model has been formulated, based on similar premises, but adopting a finite element approach. This considers the population of pellets to be distributed in a range of size classes. Growth results in movement to classes of increasing pellet size, while fragments enter the smallest size class, from which they grow to form further pellets. The revised model is stable and predicts changes in the distribution of pellet sizes within a population growing in liquid batch culture. It considers pellet growth and death, with fragmentation providing new centers of growth within the pellet population, and predicts the effects of shear forces on pellet growth and size distribution. Predictions of pellet size distributions are tested using previously published data on the growth of fungal pellets and further predictions are generated which are suitable for experimental testing using cultures of filamentous fungi or actinomycetes. (c) 1995 John Wiley & Sons, Inc.     

Analysis of a time-delayed mathematical model for tumour growth with an almost periodic supply of external nutrients

In this paper, the existence, uniqueness and exponential stability of almost periodic solutions for a mathematical model of tumour growth are studied. The establishment of the model is based on the reaction–diffusion dynamics and mass conservation law and is considered with a delay in the cell proliferation process. Using a fixed-point theorem in cones, the existence and uniqueness of almost periodic solutions for different parameter values of the model is proved. Moreover, by the Gronwall inequality, sufficient conditions are established for the exponential stability of the unique almost periodic solution. Results are illustrated by computer simulations.     

Temperature dependence of vegetative growth and dark respiration: a mathematical model

A mathematical model of the processes involved in carbon metabolism is described that predicts the influence of temperature on the growth of plants. The model assumes that the rate of production of dry matter depends both on the temperature and the level of nonstructural carbohydrate. The level of nonstructural carbohydrate is determined by the rates of photosynthesis, growth, and maintenance respiration. The model describes the rate of growth and dark respiration, and the levels of carbohydrate seen in vegetative growth of carnation and tomato. The model suggests that the growth of plants at low temperatures is limited by a shortage of respiratory energy, whereas at high temperatures growth is limited by the shortage of carbohydrate. Thermoperiodism, wherein a warm day and cool night results in faster growth than does constant temperature, is explained by the model as an increase in the level of nonstructural carbohydrate which promotes the rate of growth relative to the rate of maintenance respiration.     

A mathematical model for the growth of the abdominal aortic aneurysm

We present the first mathematical model to account for the evolution of the abdominal aortic aneurysm. The artery is modelled as a two-layered, cylindrical membrane using nonlinear elasticity and a physiologically realistic constitutive model. It is subject to a constant systolic pressure and a physiological axial prestretch. The development of the aneurysm is assumed to be a consequence of the remodelling of its material constituents. Microstructural recruitment and fibre density variables for the collagen are introduced into the strain energy density functions. This enables the remodelling of collagen to be addressed as the aneurysm enlarges. An axisymmetric aneurysm, with axisymmetric degradation of elastin and linear differential equations for the remodelling of the fibre variables, is simulated numerically. Using physiologically determined parameters to model the abdominal aorta and realistic remodelling rates for its constituents, the predicted dilations of the aneurysm are consistent with those observed in vivo. An asymmetric aneurysm with spinal contact is also modelled, and the stress distributions are consistent with previous studies.