Multimodal therapy for complete regression of malignant melanoma using constrained nonlinear optimal dynamic inversion

Computational systems biology approach for permanent tumor elimination and normal tissue protection using negative biasing: Experimental validation in malignant melanoma as case study

Complete spontaneous tumor regression (without treatment) is well documented to occur in animals and humans as epidemiological analysis show, whereby the malignancy is permanently eliminated. We have developed a novel computational systems biology model for this unique phenomenon to furnish insight into the possibility of therapeutically replicating such regression processes on tumors clinically, without toxic side effects. We have formulated oncological informatics approach using cell-kinetics coupled differential equations while protecting normal tissue. We investigated three main tumor-lysis components: (ⅰ) DNA blockade factors, (ⅱ) Interleukin-2 (IL-2), and (ⅲ) Cytotoxic T-cells (CD8⁺ T). We studied the temporal variations of these factors, utilizing preclinical experimental investigations on malignant tumors, using mammalian melanoma microarray and histiocytoma immunochemical assessment. We found that permanent tumor regression can occur by: 1) Negative-Bias shift in population trajectory of tumor cells, eradicating them under first-order asymptotic kinetics, and 2) Temporal alteration in the three antitumor components (DNA replication-blockade, Antitumor T-lymphocyte, IL-2), which are respectively characterized by the following patterns: (a) Unimodal Inverted-U function, (b) Bimodal M-function, (c) Stationary-step function. These provide a time-wise orchestrated tri-phasic cytotoxic profile. We have also elucidated gene-expression levels corresponding to the above three components: (ⅰ) DNA-damage G2/M checkpoint regulation [genes: CDC2-CHEK], (ⅱ) Chemokine signaling: IL-2/15 [genes: IL2RG-IKT3], (ⅲ) T-lymphocyte signaling (genes: TRGV5-CD28). All three components quantitatively followed the same activation profiles predicted by our computational model (Smirnov-Kolmogorov statistical test satisfied, α = 5%). We have shown that the genes CASP7-GZMB are signatures of Negative-bias dynamics, enabling eradication of the residual tumor. Using the negative-biasing principle, we have furnished the dose-time profile of equivalent therapeutic agents (DNA-alkylator, IL-2, T-cell input) so that melanoma tumor may therapeutically undergo permanent extinction by replicating the spontaneous tumor regression dynamics.

Tumor-microenvironmental-response Bi-functional molecules for efficient imaging and anti-tumor activity therapy

To improve the visualization and potency of anticancer agents, the diagnosis and treatment integration bi-functional molecules were constructed based on active candidate BD7, approved drug Linifanib, and monoclonal antibody Bevacizumab. Commercial available Rhodamine B was inducted to realize imaging-aided diagnosis and target efficiency monitoring for cancer cells. In order to maintain the anticancer activity of drugs, disulfide bond was incorporated as releasable group based on tumor microenviroment. After design, synthesis and structure characterization of title compounds, various biological evaluation and cancer cell imaging analysis were carried out. The results indicated that these title diagnosis and treatment integration bi-functional molecules exhibited comparable potency with that of corresponding parent drug. Meanwhile, these agents afforded good performance in cell imaging and could be used to differentiate cancer cells from normal ovarian cells in real time. Further optimization of these bi-functional molecules is ongoing to improve the potency and precision and will be reported in due course. Our findings are expected to achieve efficient screening and real-time prognostic monitoring under the premise of high anti-tumor activity for clinical application.

Differences in predictions of ODE models of tumor growth: a cautionary example

BACKGROUND

While mathematical models are often used to predict progression of cancer and treatment outcomes, there is still uncertainty over how to best model tumor growth. Seven ordinary differential equation (ODE) models of tumor growth (exponential, Mendelsohn, logistic, linear, surface, Gompertz, and Bertalanffy) have been proposed, but there is no clear guidance on how to choose the most appropriate model for a particular cancer.

METHODS

We examined all seven of the previously proposed ODE models in the presence and absence of chemotherapy. We derived equations for the maximum tumor size, doubling time, and the minimum amount of chemotherapy needed to suppress the tumor and used a sample data set to compare how these quantities differ based on choice of growth model.

RESULTS

We find that there is a 12-fold difference in predicting doubling times and a 6-fold difference in the predicted amount of chemotherapy needed for suppression depending on which growth model was used.

CONCLUSION

Our results highlight the need for careful consideration of model assumptions when developing mathematical models for use in cancer treatment planning.