Mathematical modeling of solid cancer growth with angiogenesis

<p>Abstract</p> <p>Background</p> <p>Cancer arises when within a single cell multiple malfunctions of control systems occur, which are, broadly, the system that promote cell growth and the system that protect against erratic growth. Additional systems within the cell must be corrupted so that a cancer cell, to form a mass of any real size, produces substances that promote the growth of new blood vessels. Multiple mutations are required before a normal cell can become a cancer cell by corruption of multiple growth-promoting systems.</p> <p>Methods</p> <p>We develop a simple mathematical model to describe the solid cancer growth dynamics inducing angiogenesis in the absence of cancer controlling mechanisms.</p> <p>Results</p> <p>The initial conditions supplied to the dynamical system consist of a perturbation in form of pulse: The origin of cancer cells from normal cells of an organ of human body. Thresholds of interacting parameters were obtained from the steady states analysis. The existence of two equilibrium points determine the strong dependency of dynamical trajectories on the initial conditions. The thresholds can be used to control cancer.</p> <p>Conclusions</p> <p>Cancer can be settled in an organ if the following combination matches: better fitness of cancer cells, decrease in the efficiency of the repairing systems, increase in the capacity of sprouting from existing vascularization, and higher capacity of mounting up new vascularization. However, we show that cancer is rarely induced in organs (or tissues) displaying an efficient (numerically and functionally) reparative or regenerative mechanism.</p

Tumor Growth Rate Determines the Timing of Optimal Chronomodulated Treatment Schedules

In host and cancer tissues, drug metabolism and susceptibility to drugs vary in a circadian (24 h) manner. In particular, the efficacy of a cell cycle specific (CCS) cytotoxic agent is affected by the daily modulation of cell cycle activity in the target tissues. Anti-cancer chronotherapy, in which treatments are administered at a particular time each day, aims at exploiting these biological rhythms to reduce toxicity and improve efficacy of the treatment. The circadian status, which is the timing of physiological and behavioral activity relative to daily environmental cues, largely determines the best timing of treatments. However, the influence of variations in tumor kinetics has not been considered in determining appropriate treatment schedules. We used a simple model for cell populations under chronomodulated treatment to identify which biological parameters are important for the successful design of a chronotherapy strategy. We show that the duration of the phase of the cell cycle targeted by the treatment and the cell proliferation rate are crucial in determining the best times to administer CCS drugs. Thus, optimal treatment times depend not only on the circadian status of the patient but also on the cell cycle kinetics of the tumor. Then, we developed a theoretical analysis of treatment outcome (TATO) to relate the circadian status and cell cycle kinetic parameters to the treatment outcomes. We show that the best and the worst CCS drug administration schedules are those with 24 h intervals, implying that 24 h chronomodulated treatments can be ineffective or even harmful if administered at wrong circadian times. We show that for certain tumors, administration times at intervals different from 24 h may reduce these risks without compromising overall efficacy

Dickkopf1 Regulates Fate Decision and Drives Breast Cancer Stem Cells to Differentiation: An Experimentally Supported Mathematical Model

BACKGROUND: Modulation of cellular signaling pathways can change the replication/differentiation balance in cancer stem cells (CSCs), thus affecting tumor growth and recurrence. Analysis of a simple, experimentally verified, mathematical model suggests that this balance is maintained by quorum sensing (QS). METHODOLOGY/PRINCIPAL FINDINGS: To explore the mechanism by which putative QS cellular signals in mammary stem cells (SCs) may regulate SC fate decisions, we developed a multi-scale mathematical model, integrating extra-cellular and intra-cellular signal transduction within the mammary tissue dynamics. Preliminary model analysis of the single cell dynamics indicated that Dickkopf1 (Dkk1), a protein known to negatively regulate the Wnt pathway, can serve as anti-proliferation and pro-maturation signal to the cell. Simulations of the multi-scale tissue model suggested that Dkk1 may be a QS factor, regulating SC density on the level of the whole tissue: relatively low levels of exogenously applied Dkk1 have little effect on SC numbers, whereas high levels drive SCs into differentiation. To verify these model predictions, we treated the MCF-7 cell line and primary breast cancer (BC) cells from 3 patient samples with different concentrations and dosing regimens of Dkk1, and evaluated subsequent formation of mammospheres (MS) and the mammary SC marker CD44(+)CD24(lo). As predicted by the model, low concentrations of Dkk1 had no effect on primary BC cells, or even increased MS formation among MCF-7 cells, whereas high Dkk1 concentrations decreased MS formation among both primary BC cells and MCF-7 cells. CONCLUSIONS/SIGNIFICANCE: Our study suggests that Dkk1 treatment may be more robust than other methods for eliminating CSCs, as it challenges a general cellular homeostasis mechanism, namely, fate decision by QS. The study also suggests that low dose Dkk1 administration may be counterproductive; we showed experimentally that in some cases it can stimulate CSC proliferation, although this needs validating in vivo

Exploiting Clinical Trial Data Drastically Narrows the Window of Possible Solutions to the Problem of Clinical Adaptation of a Multiscale Cancer Model

The development of computational models for simulating tumor growth and response to treatment has gained significant momentum during the last few decades. At the dawn of the era of personalized medicine, providing insight into complex mechanisms involved in cancer and contributing to patient-specific therapy optimization constitute particularly inspiring pursuits. The in silico oncology community is facing the great challenge of effectively translating simulation models into clinical practice, which presupposes a thorough sensitivity analysis, adaptation and validation process based on real clinical data. In this paper, the behavior of a clinically-oriented, multiscale model of solid tumor response to chemotherapy is investigated, using the paradigm of nephroblastoma response to preoperative chemotherapy in the context of the SIOP/GPOH clinical trial. A sorting of the model's parameters according to the magnitude of their effect on the output has unveiled the relative importance of the corresponding biological mechanisms; major impact on the result of therapy is credited to the oxygenation and nutrient availability status of the tumor and the balance between the symmetric and asymmetric modes of stem cell division. The effect of a number of parameter combinations on the extent of chemotherapy-induced tumor shrinkage and on the tumor's growth rate are discussed. A real clinical case of nephroblastoma has served as a proof of principle study case, demonstrating the basics of an ongoing clinical adaptation and validation process. By using clinical data in conjunction with plausible values of model parameters, an excellent fit of the model to the available medical data of the selected nephroblastoma case has been achieved, in terms of both volume reduction and histological constitution of the tumor. In this context, the exploitation of multiscale clinical data drastically narrows the window of possible solutions to the clinical adaptation problem

2012 Activity Report of the Regional Research Programme on Hadrontherapy for the ETOILE Center

2012 is the penultimate year of financial support by the CPER 2007-2013 for ETOILE's research program, sustained by the PRRH at the University Claude Bernard. As with each edition we make the annual review of the research in this group, so active for over 12 years now. Over the difficulties in the decision-making process for the implementation of the ETOILE Center, towards which all our efforts are focussed, some "themes" (work packages) were strengthened, others have progressed, or have been dropped. This is the case of the eighth theme (technological developments), centered around the technology for rotative beam distribution heads (gantries) and, after being synchronized with the developments of ULICE's WP6, remained so by ceasing its activities, coinciding also with the retirement of its historic leader at IPNL, Marcel Bajard. Topic number 5 ("In silico simulations") has suffered the departure of its leader, Benjamin Ribba, although the work has still been provided by Branka Bernard, a former postdoctoral fellow in Lyon Sud, and now back home in Croatia, still in contract with UCBL for the ULICE project. Aside from these two issues (and the fact that the theme "Medico-economical simulations" is now directly linked to the first one ("Medical Project"), the rest of the teams are growing, as evidenced by the publication statistics at the beginning of this report. This is obviously due to the financial support of our always faithful regional institutions, but also to the synergy that the previous years, the European projects, the arrival of the PRIMES LabEx, and the national France Hadron infrastructure have managed to impulse. The Rhone-Alpes hadron team, which naturally includes the researchers of LPC at Clermont, should also see its influence result in a strong presence in France Hadron's regional node, which is being organized. The future of this regional research is not yet fully guaranteed, especially in the still uncertain context of ETOILE, but the tracks are beginning to emerge to allow past and present efforts translate into a long future that we all want to see established. Each of the researchers in PRRH is aware that 2013 will be (and already is) the year of great challenge : for ETOILE, for the PRRH, for hadron therapy in France, for French hadrontherapy in Europe (after the opening and beginning of treatments in the German [HIT Heidelberg, Marburg], Italian [CNAO, Pavia] and Austrian [MedAustron, Wien Neuerstadt]) centers. Let us meet again in early 2014 for a comprehensive review of the past and a perspective for the future ..

A tumor cord model for Doxorubicin delivery and dose optimization in solid tumors

<p>Abstract</p> <p>Background</p> <p>Doxorubicin is a common anticancer agent used in the treatment of a number of neoplasms, with the lifetime dose limited due to the potential for cardiotoxocity. This has motivated efforts to develop optimal dosage regimes that maximize anti-tumor activity while minimizing cardiac toxicity, which is correlated with peak plasma concentration. Doxorubicin is characterized by poor penetration from tumoral vessels into the tumor mass, due to the highly irregular tumor vasculature. I model the delivery of a soluble drug from the vasculature to a solid tumor using a tumor cord model and examine the penetration of doxorubicin under different dosage regimes and tumor microenvironments.</p> <p>Methods</p> <p>A coupled ODE-PDE model is employed where drug is transported from the vasculature into a tumor cord domain according to the principle of solute transport. Within the tumor cord, extracellular drug diffuses and saturable pharmacokinetics govern uptake and efflux by cancer cells. Cancer cell death is also determined as a function of peak intracellular drug concentration.</p> <p>Results</p> <p>The model predicts that transport to the tumor cord from the vasculature is dominated by diffusive transport of free drug during the initial plasma drug distribution phase. I characterize the effect of all parameters describing the tumor microenvironment on drug delivery, and large intercapillary distance is predicted to be a major barrier to drug delivery. Comparing continuous drug infusion with bolus injection shows that the optimum infusion time depends upon the drug dose, with bolus injection best for low-dose therapy but short infusions better for high doses. Simulations of multiple treatments suggest that additional treatments have similar efficacy in terms of cell mortality, but drug penetration is limited. Moreover, fractionating a single large dose into several smaller doses slightly improves anti-tumor efficacy.</p> <p>Conclusion</p> <p>Drug infusion time has a significant effect on the spatial profile of cell mortality within tumor cord systems. Therefore, extending infusion times (up to 2 hours) and fractionating large doses are two strategies that may preserve or increase anti-tumor activity and reduce cardiotoxicity by decreasing peak plasma concentration. However, even under optimal conditions, doxorubicin may have limited delivery into advanced solid tumors.</p

A quantitative systems pharmacology consortium approach to managing immunogenicity of therapeutic proteins

Immunogenicity is a major challenge in drug development and patient care. Currently, most efforts are dedicated to the elimination of the unwanted immune responses through T‐cell epitope prediction and protein engineering. However, because it is unlikely that this approach will lead to complete eradication of immunogenicity, we propose that quantitative systems pharmacology models should be developed to predict and manage immunogenicity. The potential impact of such a mechanistic model‐based approach is precedented by applications of physiologically‐based pharmacokinetics