Mathematical modeling of the cancer cell’s control circuitry : paving the way to individualized therapeutic strategies

Cancer is a disease of signal transduction in which the dysregulation of the network of intracellular and extracellular signaling cascades is sufficient to thwart the cells finely-tuned biochemical control mechanisms. A keen interest in the mathematical modeling of cell signaling networks and the regulation of signal transduction has emerged in recent years, and has produced a glimmer of insight into the sophisticated feedback control and network regulation operating within cells. In this review, we present an overview of published theoretical studies on the control aspects of signal transduction, emphasizing the role and importance of mechanisms such as ‘ultrasensitivity’ and feedback loops. We emphasize that these exquisite and often subtle control strategies represent the key to orchestrating ‘simple’ signaling behaviors within the complex intracellular network, while regulating the trade-off between sensitivity and robustness to internal and external perturbations. Through a consideration of these apparent paradoxes, we explore how the basic homeostasis of the intracellular signaling network, in the face of carcinogenesis, can lead to neoplastic progression rather than cell death. A simple mathematical model is presented, furnishing a vivid illustration of how ‘control-oriented’ models of the deranged signaling networks in cancer cells may enucleate improved treatment strategies, including patient-tailored combination therapies, with the potential for reduced toxicity and more robust and potent antitumor activity.

On the mathematical modeling of wound healing angiogenesis in skin as a reaction-transport process

Over the last 30 years, numerous research groups have attempted to provide mathematical descriptions of the skin wound healing process. The development of theoretical models of the interlinked processes that underlie the healing mechanism has yielded considerable insight into aspects of this critical phenomenon that remain difficult to investigate empirically. In particular, the mathematical modeling of angiogenesis, i.e., capillary sprout growth, has offered new paradigms for the understanding of this highly complex and crucial step in the healing pathway. With the recent advances in imaging and cell tracking, the time is now ripe for an appraisal of the utility and importance of mathematical modeling in wound healing angiogenesis research. The purpose of this review is to pedagogically elucidate the conceptual principles that have underpinned the development of mathematical descriptions of wound healing angiogenesis, specifically those that have utilized a continuum reaction-transport framework, and highlight the contribution that such models have made toward the advancement of research in this field. We aim to draw attention to the common assumptions made when developing models of this nature, thereby bringing into focus the advantages and limitations of this approach. A deeper integration of mathematical modeling techniques into the practice of wound healing angiogenesis research promises new perspectives for advancing our knowledge in this area. To this end we detail several open problems related to the understanding of wound healing angiogenesis, and outline how these issues could be addressed through closer cross-disciplinary collaboration.

The Mathematical modeling of inhomogeneites in ventricular tissue

Cardiac arrhythmias such as ventricular tachycardia (VT) or ventricular fibrillation (VF) are the leading cause of death in the industrialised world. There is a growing consensus that these arrhythmias arise because of the formation of spiral waves of electrical activation in cardiac tissue; unbroken spiral waves are associated with VT and broken ones with VF. Several experimental studies have been carried out to determine the effects of inhomogeneities in cardiac tissue on such arrhythmias. We give a brief overview of such experiments, and then an introduction to partial-differential-equation models for ventricular tissue. We show how different types of inhomogeneities can be included in such models, and then discuss various numerical studies, including our own, of the effects of these inhomogeneities on spiral-wave dynamics. The most remarkable qualitative conclusion of our studies is that the spiral-wave dynamics in such systems depends very sensitively on the positions of these inhomogeneities.

Mathematical Modeling of Near-Wall Flows of Two-Phase Mixture with Evaporating Droplets

In the framework of the two-continuum approach, using the matched asymptotic expansion method, the equations of a laminar boundary layer in mist flows with evaporating droplets were derived and solved. The similarity criteria controlling the mist flows were determined. For the flow along a curvilinear surface, the forms of the boundary layer equations differ from the regimes of presence and absence of the droplet inertia deposition. The numerical results were presented for the vapor-droplet boundary layer in the neighborhood of a stagnation point of a hot blunt body. It is demonstrated that, due to evaporation, a droplet-free region develops near the wall inside the boundary layer. On the upper edge of this region, the droplet radius tends to zero and the droplet number density becomes much higher than that in the free stream. The combined effect of the droplet evaporation and accumulation results in a significant enhancement of the heat transfer on the surface even for small mass concentration of the droplets in the free stream. 在双连续介质理论框架下,采用匹配渐进展开方法导出并求解了具有蒸发液滴的汽雾流中层流边界层方程,给出了控制汽雾流的相似判据。对于沿曲面的流动,边界层方程的形式取决于是否存在液滴的惯性沉积。给出了热钝体验点附近蒸汽。液滴边界层的数值计算结果。它们表明：由于蒸发,在边界层内近壁处形成了一个无液滴区域；在该区上边界处,液滴半径趋于零而液滴数密度急剧增高。液滴蒸发及聚集的联合效应造成了表面热流的显著增加,甚至在自由来流中液滴质量浓度很低时此效应依然存在。

Mathematical modeling of colony formation in algal blooms: phenotypic plasticity in cyanobacteria

In this paper, we analyzed a mathematical model of algal-grazer dynamics, including the effect of colony formation, which is an example of phenotypic plasticity. The model consists of three variables, which correspond to the biomasses of unicellular algae, colonial algae, and herbivorous zooplankton. Among these organisms, colonial algae are the main components of algal blooms. This aquatic system has two stable attractors, which can be identified as a zooplankton-dominated (ZD) state and an algal-dominated (AD) state, respectively. Assuming that the handling time of zooplankton on colonial algae increases with the colonial algae biomass, we discovered that bistability can occur within the model system. The applicability of alternative stable states in algae-grazer dynamics as a framework for explaining the algal blooms in real lake ecosystems, thus, seems to depend on whether the assumption mentioned above is met in natural circumstances.

Mathematical modeling of the ester oil-refrigerant R134A mixture two-phase flow with foam formation through a small diameter tube

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Basic ingredients for mathematical modeling of tumor growth in vitro: Cooperative effects and search for space

Based on the literature data from HT-29 cell monolayers, we develop a model for its growth, analogous to an epidemic model, mixing local and global interactions. First, we propose and solve a deterministic equation for the progress of these colonies. Thus, we add a stochastic (local) interaction and simulate the evolution of an Eden-like aggregate by using dynamical Monte Carlo methods. The growth curves of both deterministic and stochastic models are in excellent agreement with the experimental observations. The waiting times distributions, generated via our stochastic model, allowed us to analyze the role of mesoscopic events. We obtain log-normal distributions in the initial stages of the growth and Gaussians at long times. We interpret these outcomes in the light of cellular division events: in the early stages, the phenomena are dependent each other in a multiplicative geometric-based process, and they are independent at long times. We conclude that the main ingredients for a good minimalist model of tumor growth, at mesoscopic level, are intrinsic cooperative mechanisms and competitive search for space. © 2013 Elsevier Ltd.

Mathematical modeling of a continuous alcoholic fermentation process in a two-stage tower reactor cascade with flocculating yeast recycle

Experiments of continuous alcoholic fermentation of sugarcane juice with flocculating yeast recycle were conducted in a system of two 0.22-L tower bioreactors in series, operated at a range of dilution rates (D (1) = D (2) = 0.27-0.95 h(-1)), constant recycle ratio (alpha = F (R) /F = 4.0) and a sugar concentration in the feed stream (S (0)) around 150 g/L. The data obtained in these experimental conditions were used to adjust the parameters of a mathematical model previously developed for the single-stage process. This model considers each of the tower bioreactors as a perfectly mixed continuous reactor and the kinetics of cell growth and product formation takes into account the limitation by substrate and the inhibition by ethanol and biomass, as well as the substrate consumption for cellular maintenance. The model predictions agreed satisfactorily with the measurements taken in both stages of the cascade. The major differences with respect to the kinetic parameters previously estimated for a single-stage system were observed for the maximum specific growth rate, for the inhibition constants of cell growth and for the specific rate of substrate consumption for cell maintenance. Mathematical models were validated and used to simulate alternative operating conditions as well as to analyze the performance of the two-stage process against that of the single-stage process.