Cell viability prediction and optimization in extrusion-based bioprinting via neural network-based Bayesian optimization models.

The fields of regenerative medicine and cancer modeling have witnessed tremendous growth in the application of 3D bioprinting. Maintaining high cell viability throughout the bioprinting process is crucial for the success of this technology, as it directly affects the accuracy of the 3D bioprinted models, the validity of experimental results, and the discovery of new therapeutic approaches. Therefore, optimizing bioprinting conditions, which include numerous variables influencing cell viability during and after the procedure, is of utmost importance to achieve desirable results. So far, these optimizations have been accomplished primarily through trial and error and repeating multiple time-consuming and costly experiments. To address this challenge, we initiated the process by creating a dataset of these parameters for gelatin and alginate-based bioinks and the corresponding cell viability by integrating data obtained in our laboratory and those derived from the literature. Then, we developed machine learning models to predict cell viability based on different bioprinting variables. The trained neural network yielded regressionR2value of 0.71 and classification accuracy of 0.86. Compared to models that have been developed so far, the performance of our models is superior and shows great prediction results. The study further introduces a novel optimization strategy that employs the Bayesian optimization model in combination with the developed regression neural network to determine the optimal combination of the selected bioprinting parameters to maximize cell viability and eliminate trial-and-error experiments. Finally, we experimentally validated the optimization model's performance.

Computational Multi-Scale Modeling of Drug Delivery into an Anti-Angiogenic Therapy-Treated Tumor.

The present study develops a numerical model, which is the most complex one, in comparison to previous research to investigate drug delivery accompanied by the anti-angiogenesis effect. This paper simulates intravascular blood flow and interstitial fluid flow using a dynamic model. The model accounts for the non-Newtonian behavior of blood and incorporates the adaptation of the diameter of a heterogeneous microvascular network derived from modeling the evolution of endothelial cells toward a circular tumor sprouting from two-parent vessels, with and without imposing the inhibitory effect of angiostatin on a modified discrete angiogenesis model. The average solute exposure and its uniformity in solid tumors of different sizes are studied by numerically solving the convection-diffusion equation. Three different methodologies are considered for simulating anti-angiogenesis: modifying the capillary network, updating the transport properties, and considering both microvasculature and transport properties modifications. It is shown that anti-angiogenic therapy decreases drug wash-out in the periphery of the tumor. Results show the decisive role of microvascular structure, particularly its distribution, and interstitial transport properties modifications induced via vascular normalization on the quality of drug delivery, such that it is improved by 39% in uniformity by the second approach in R = 0.2 cm.

Predicting intratumoral fluid pressure and liposome accumulation using physics informed deep learning.

Liposome-based anticancer agents take advantage of the increased vascular permeability and transvascular pressure gradients for selective accumulation in tumors, a phenomenon known as the enhanced permeability and retention(EPR) effect. The EPR effect has motivated the clinical use of nano-therapeutics, with mixed results on treatment outcome. High interstitial fluid pressure (IFP) has been shown to limit liposome drug delivery to central tumour regions. Furthermore, high IFP is an independent prognostic biomarker for treatment efficacy in radiation therapy and chemotherapy for some solid cancers. Therefore, accurately measuring spatial liposome accumulation and IFP distribution within a solid tumour is crucial for optimal treatment planning. In this paper, we develop a model capable of predicting voxel-by-voxel intratumoral liposome accumulation and IFP using pre and post administration imaging. Our approach is based on physics informed machine learning, a novel technique combining machine learning and partial differential equations. through application to a set of mouse data and a set of synthetically-generated tumours, we show that our approach accurately predicts the spatial liposome accumulation and IFP for an individual tumour while relying on minimal information. This is an important result with applications for forecasting tumour progression and designing treatment.

Controlled tumor heterogeneity in a co-culture system by 3D bio-printed tumor-on-chip model.

Cancer treatment resistance is a caused by presence of various types of cells and heterogeneity within the tumor. Tumor cell-cell and cell-microenvironment interactions play a significant role in the tumor progression and invasion, which have important implications for diagnosis, and resistance to chemotherapy. In this study, we develop 3D bioprinted in vitro models of the breast cancer tumor microenvironment made of co-cultured cells distributed in a hydrogel matrix with controlled architecture to model tumor heterogeneity. We hypothesize that the tumor could be represented by a cancer cell-laden co-culture hydrogel construct, whereas its microenvironment can be modeled in a microfluidic chip capable of producing a chemical gradient. Breast cancer cells (MCF7 and MDA-MB-231) and non-tumorigenic mammary epithelial cells (MCF10A) were embedded in the alginate-gelatine hydrogels and printed using a multi-cartridge extrusion bioprinter. Our approach allows for precise control over position and arrangements of cells in a co-culture system, enabling the design of various tumor architectures. We created samples with two different types of cells at specific initial locations, where the density of each cell type was carefully controlled. The cells were either randomly mixed or positioned in sequential layers to create cellular heterogeneity. To study cell migration toward chemoattractant, we developed a chemical microenvironment in a chamber with a gradual chemical gradient. As a proof of concept, we studied different migration patterns of MDA-MB-231 cells toward the epithelial growth factor gradient in presence of MCF10A cells in different ratios using this device. Our approach involves the integration of 3D bioprinting and microfluidic devices to create diverse tumor architectures that are representative of those found in various patients. This provides an excellent tool for studying the behavior of cancer cells with high spatial and temporal resolution.

Investigation of the evolution of tumor-induced microvascular network under the inhibitory effect of anti-angiogenic factor, angiostatin: A mathematical study.

Anti-angiogenesis as a treatment strategy for normalizing the microvascular network of tumors is of great interest among researchers, especially in combination with chemotherapy or radiotherapy. According to the vital role that angiogenesis plays in tumor growth and in exposing the tumor to therapeutic agents, this work develops a mathematical framework to study the influence of angiostatin, a plasminogen fragment that shows the anti-angiogenic function, in the evolutionary behavior of tumor-induced angiogenesis. Angiostatin-induced microvascular network reformation is investigated in a two-dimensional space by considering two parent vessels around a circular tumor by a modified discrete angiogenesis model in different tumor sizes. The effects of imposing modifications on the existing model, i.e., the matrix-degrading enzyme effect, proliferation and death of endothelial cells, matrix density function, and a more realistic chemotactic function, are investigated in this study. Results show a decrease in microvascular density in response to the angiostatin. A functional relationship exists between angiostatin's ability to normalize the capillary network and tumor size or progression stage, such that capillary density decreases by 55%, 41%, 24%, and 13% in tumors with a non-dimensional radius of 0.4, 0.3, 0.2, and 0.1, respectively, after angiostatin administration.

Predicting and elucidating the post-printing behavior of 3D printed cancer cells in hydrogel structures by integrating in-vitro and in-silico experiments.

A key feature distinguishing 3D bioprinting from other 3D cell culture techniques is its precise control over created structures. This property allows for the high-resolution fabrication of biomimetic structures with controlled structural and mechanical properties such as porosity, permeability, and stiffness. However, analyzing post-printing cellular dynamics and optimizing their functions within the 3D fabricated environment is only possible through trial and error and replicating several experiments. This issue motivated the development of a cellular automata model for the first time to simulate post-printing cell behaviour within the 3D bioprinted construct. To improve our model, we bioprinted a 3D construct using MDA-MB-231 cell-laden hydrogel and evaluated cellular functions, including viability and proliferation in 11 days. The results showed that our model successfully simulated the 3D bioprinted structure and captured in-vitro observations. We demonstrated that in-silico model could predict and elucidate post-printing biological functions for different initial cell numbers in bioink and different bioink formulations with gelatine and alginate, without replicating several costly and time-consuming in-vitro measurements. We believe such a computational framework will substantially impact 3D bioprinting's future application. We hope this study inspires researchers to further realize how an in-silico model might be utilized to advance in-vitro 3D bioprinting research.

Deep learning characterization of brain tumours with diffusion weighted imaging.

Glioblastoma multiforme (GBM) is one of the most deadly forms of cancer. Methods of characterizing these tumours are valuable for improving predictions of their progression and response to treatment. A mathematical model called the proliferation-invasion (PI) model has been used extensively in the literature to model the growth of these tumours, though it relies on known values of two key parameters: the tumour cell diffusivity and proliferation rate. Unfortunately, these parameters are difficult to estimate in a patient-specific manner, making personalized tumour forecasting challenging. In this paper, we develop and apply a deep learning model capable of making accurate estimates of these key GBM-characterizing parameters while simultaneously producing a full prediction of the tumour progression curve. Our method uses two sets of multi sequence MRI in order to produce estimations and relies on a preprocessing pipeline which includes brain tumour segmentation and conversion to tumour cellularity. We first apply our deep learning model to synthetic tumours to showcase the model's capabilities and identify situations where prediction errors are likely to occur. We then apply our model to a clinical dataset consisting of five patients diagnosed with GBM. For all patients, we derive evidence-based estimates for each of the PI model parameters and predictions for the future progression of the tumour, along with estimates of the parameter uncertainties. Our work provides a new, easily generalizable method for the estimation of patient-specific tumour parameters, which can be built upon to aid physicians in designing personalized treatments.

Temporal optimization of radiation therapy to heterogeneous tumour populations and cancer stem cells.

External beam radiation therapy is a key part of modern cancer treatments which uses high doses of radiation to destroy tumour cells. Despite its widespread usage and extensive study in theoretical, experimental, and clinical works, many questions still remain about how best to administer it. Many mathematical studies have examined optimal scheduling of radiotherapy, and most come to similar conclusions. Importantly though, these studies generally assume intratumoral homogeneity. But in recent years, it has become clear that tumours are not homogeneous masses of cancerous cells, but wildly heterogeneous masses with various subpopulations which grow and respond to treatment differently. One subpopulation of particular importance is cancer stem cells (CSCs) which are known to exhibit higher radioresistence compared with non-CSCs. Knowledge of these differences between cell types could theoretically lead to changes in optimal treatment scheduling. Only a few studies have examined this question, and interestingly, they arrive at apparent conflicting results. However, an understanding of their assumptions reveals a key difference which leads to their differing conclusions. In this paper, we generalize the problem of temporal optimization of dose distribution of radiation therapy to a two cell type model. We do so by creating a mathematical model and a numerical optimization algorithm to find the distribution of dose which leads to optimal cell kill. We then create a data set of optimization solutions and use data analysis tools to learn the relationships between model parameters and the qualitative behaviour of optimization results. Analysis of the model and discussion of biological importance are provided throughout. We find that the key factor in predicting the behaviour of the optimal distribution of radiation is the ratio between the radiosensitivities of the present cell types. These results can provide guidance for treatment in cases where clinicians have knowledge of tumour heterogeneity and of the abundance of CSCs.

An integrative systems biology approach to overcome venetoclax resistance in acute myeloid leukemia.

The over-expression of the Bcl-2 protein is a common feature of many solid cancers and hematological malignancies, and it is typically associated with poor prognosis and resistance to chemotherapy. Bcl-2-specific inhibitors, such as venetoclax, have recently been approved for the treatment of chronic lymphocytic leukemia and small lymphocytic lymphoma, and they are showing promise in clinical trials as a targeted therapy for patients with relapsed or refractory acute myeloid leukemia (AML). However, successful treatment of AML with Bcl-2-specific inhibitors is often followed by the rapid development of drug resistance. An emerging paradigm for overcoming drug resistance in cancer treatment is through the targeting of mitochondrial energetics and metabolism. In AML in particular, it was recently observed that inhibition of mitochondrial translation via administration of the antibiotic tedizolid significantly affects mitochondrial bioenergetics, activating the integrated stress response (ISR) and subsequently sensitizing drug-resistant AML cells to venetoclax. Here we develop an integrative systems biology approach to acquire a deeper understanding of the molecular mechanisms behind this process, and in particular, of the specific role of the ISR in the commitment of cells to apoptosis. Our multi-scale mathematical model couples the ISR to the intrinsic apoptosis pathway in venetoclax-resistant AML cells, includes the metabolic effects of treatment, and integrates RNA, protein level, and cellular viability data. Using the mathematical model, we identify the dominant mechanisms by which ISR activation helps to overcome venetoclax resistance, and we study the temporal sequencing of combination treatment to determine the most efficient and robust combination treatment protocol.

Network Analysis Identifies Phase Transitions for Tumor With Interacting Cells.

Metastasis is the process by which cancer cells acquire the capability to leave the primary tumor and travel to distant sites. Recent experiments have suggested that the epithelial-mesenchymal transition can regulate invasion and metastasis. Another possible scenario is the collective motion of cells. Recent studies have also proposed a jamming-unjamming transition for epithelial cells based on physical forces. Here, we assume that there exists a short-range chemical attraction between cancer cells and employ the Brownian dynamics to simulate tumor growth. Applying the network analysis, we suggest three possible phases for a given tumor and study the transition between these phases by adjusting the attraction strength.

Circadian and Immunity Cycle Talk in Cancer Destination: From Biological Aspects to In Silico Analysis.

Cancer is the leading cause of death and a major problem to increasing life expectancy worldwide. In recent years, various approaches such as surgery, chemotherapy, radiation, targeted therapies, and the newest pillar, immunotherapy, have been developed to treat cancer. Among key factors impacting the effectiveness of treatment, the administration of drugs based on the circadian rhythm in a person and within individuals can significantly elevate drug efficacy, reduce adverse effects, and prevent drug resistance. Circadian clocks also affect various physiological processes such as the sleep cycle, body temperature cycle, digestive and cardiovascular processes, and endocrine and immune systems. In recent years, to achieve precision patterns for drug administration using computational methods, the interaction of the effects of drugs and their cellular pathways has been considered more seriously. Integrated data-derived pathological images and genomics, transcriptomics, and proteomics analyses have provided an understanding of the molecular basis of cancer and dramatically revealed interactions between circadian and immunity cycles. Here, we describe crosstalk between the circadian cycle signaling pathway and immunity cycle in cancer and discuss how tumor microenvironment affects the influence on treatment process based on individuals' genetic differences. Moreover, we highlight recent advances in computational modeling that pave the way for personalized immune chronotherapy.

The role of calcium oscillations in the phenotype selection in endothelial cells.

Angiogenesis is an important process in the formation and maintenance of tissues which is driven by a complex system of intracellular and intercellular signaling mechanisms. Endothelial cells taking part in early angiogenesis must select their phenotype as either a tip cells (leading, migratory) or a stalk cells (following). Recent experiments have demonstrated that rapid calcium oscillations within active cells characterize this phenotype selection process and that these oscillations play a necessary role in governing phenotype selection and eventual vessel architecture. In this work, we develop a mathematical model capable of describing these oscillations and their role in phenotype selection then use it to improve our understanding of the biological mechanisms at play. We developed a model based on two previously published and experimentally validated mathematical models of calcium and angiogenesis then use our resulting model to simulate various multi-cell scenarios. We are able to capture essential calcium oscillation dynamics and intercellular communication between neighboring cells. The results of our model show that although the late DLL4 (a transmembrane protein that activates Notch pathway) levels of a cell are connected with its initial IP3 (Inositol 1,4,5-trisphosphate) level, cell-to-cell communication determines its eventual phenotype.

The network structure affects the fixation probability when it couples to the birth-death dynamics in finite population.

The study of evolutionary dynamics on graphs is an interesting topic for researchers in various fields of science and mathematics. In systems with finite population, different model dynamics are distinguished by their effects on two important quantities: fixation probability and fixation time. The isothermal theorem declares that the fixation probability is the same for a wide range of graphs and it only depends on the population size. This has also been proved for more complex graphs that are called complex networks. In this work, we propose a model that couples the population dynamics to the network structure and show that in this case, the isothermal theorem is being violated. In our model the death rate of a mutant depends on its number of neighbors, and neutral drift holds only in the average. We investigate the fixation probability behavior in terms of the complexity parameter, such as the scale-free exponent for the scale-free network and the rewiring probability for the small-world network.

Reinforcement learning derived chemotherapeutic schedules for robust patient-specific therapy.

The in-silico development of a chemotherapeutic dosing schedule for treating cancer relies upon a parameterization of a particular tumour growth model to describe the dynamics of the cancer in response to the dose of the drug. In practice, it is often prohibitively difficult to ensure the validity of patient-specific parameterizations of these models for any particular patient. As a result, sensitivities to these particular parameters can result in therapeutic dosing schedules that are optimal in principle not performing well on particular patients. In this study, we demonstrate that chemotherapeutic dosing strategies learned via reinforcement learning methods are more robust to perturbations in patient-specific parameter values than those learned via classical optimal control methods. By training a reinforcement learning agent on mean-value parameters and allowing the agent periodic access to a more easily measurable metric, relative bone marrow density, for the purpose of optimizing dose schedule while reducing drug toxicity, we are able to develop drug dosing schedules that outperform schedules learned via classical optimal control methods, even when such methods are allowed to leverage the same bone marrow measurements.

Systems biology informed neural networks (SBINN) predict response and novel combinations for PD-1 checkpoint blockade.

Anti-PD-1 immunotherapy has recently shown tremendous success for the treatment of several aggressive cancers. However, variability and unpredictability in treatment outcome have been observed, and are thought to be driven by patient-specific biology and interactions of the patient's immune system with the tumor. Here we develop an integrative systems biology and machine learning approach, built around clinical data, to predict patient response to anti-PD-1 immunotherapy and to improve the response rate. Using this approach, we determine biomarkers of patient response and identify potential mechanisms of drug resistance. We develop systems biology informed neural networks (SBINN) to calculate patient-specific kinetic parameter values and to predict clinical outcome. We show how transfer learning can be leveraged with simulated clinical data to significantly improve the response prediction accuracy of the SBINN. Further, we identify novel drug combinations and optimize the treatment protocol for triple combination therapy consisting of IL-6 inhibition, recombinant IL-12, and anti-PD-1 immunotherapy in order to maximize patient response. We also find unexpected differences in protein expression levels between response phenotypes which complement recent clinical findings. Our approach has the potential to aid in the development of targeted experiments for patient drug screening as well as identify novel therapeutic targets.

A mean-field approach for modeling the propagation of perturbations in biochemical reaction networks.

Often, the time evolution of a biochemical reaction network is crucial for determining the effects of combining multiple pharmaceuticals. Here we illustrate a mathematical framework for modeling the dominant temporal behaviour of a complicated molecular pathway or biochemical reaction network in response to an arbitrary perturbation, such as resulting from the administration of a therapeutic agent. The method enables the determination of the temporal evolution of a target protein as the perturbation propagates through its regulatory network. The mathematical approach is particularly useful when the experimental data that is available for characterizing or parameterizing the regulatory network is limited or incomplete. To illustrate the method, we consider the examples of the regulatory networks for the target proteins c-Myc and Chop, which play an important role in venetoclax resistance in acute myeloid leukemia. First we show how the networks that regulate each target protein can be reduced to a mean-field model by identifying the distinct effects that groups of proteins in the regulatory network have on the target protein. Then we show how limited protein-level data can be used to further simplify the mean-field model to pinpoint the dominant effects of the network perturbation on the target protein. This enables a further reduction in the number of parameters in the model. The result is an ordinary differential equation model that captures the temporal evolution of the expression of a target protein when one or more proteins in its regulatory network have been perturbed. Finally, we show how the dominant effects predicted by the mathematical model agree with RNA sequencing data for the regulatory proteins comprising the molecular network, despite the model not having a priori knowledge of this data. Thus, while the approach gives a simplified model for the expression of the target protein, it allows for the interpretation of the effects of the perturbation on the regulatory network itself. This method can be easily extended to sets of target proteins to model components of a larger systems biology model, and provides an approach for partially integrating RNA sequencing data and protein expression data. Moreover, it is a general approach that can be used to study drug effects on specific protein(s) in any disease or condition.

The impact of random microenvironmental fluctuations on tumor control probability.

The tumor control probability (TCP) is a metric used to calculate the probability of controlling or eradicating tumors through radiotherapy. Cancer cells vary in their response to radiation, and although many factors are involved, the tumor microenvironment is a crucial one that determines radiation efficacy. The tumor microenvironment plays a significant role in cancer initiation and propagation, as well as in treatment outcome. We have developed stochastic formulations to study the impact of arbitrary microenvironmental fluctuations on TCP and extinction probability (EP), which is defined as the probability of cancer cells removal in the absence of treatment. Since the derivation of analytical solutions are not possible for complicated cases, we employ a modified Gillespie algorithm to analyze TCP and EP, considering the random variations in cellular proliferation and death rates. Our results show that increasing the standard deviation in kinetic rates initially enhances the probability of tumor eradication. However, if the EP does not reach a probability of 1, the increase in the standard deviation subsequently has a negative impact on probability of cancer cells removal, decreasing the EP over time. The greatest effect on EP has been observed when both birth and death rates are being randomly modified and are anticorrelated. In addition, similar results are observed for TCP, where radiotherapy is included, indicating that increasing the standard deviation in kinetic rates at first enhances the probability of tumor eradication. But, it has a negative impact on treatment effectiveness if the TCP does not reach a probability of 1.

A mathematical study of the impact of cell plasticity on tumour control probability.

The tumour control probability (TCP) is a treatment planning tool that evaluates the probability of tumour eradication and helps in the assessment of the relative efficacy of different radiotherapy regimens. The response of tumours to radiation differs greatly even between patients with same types of cancers. Tumour heterogeneity or cellular diversity among cancer cells has a pronounced impact on the success of administered radiotherapy protocols. Tumour heterogeneity can be explained using the cancer stem cells (CSCs) hypothesis, which posits that CSCs are responsible for tumour initiation and propagation as well as therapeutic resistance. Moreover, the existence of plasticity or bidirectional transition between CSCs and non-CSCs indicates that, sometimes, non-CSCs appear to mimic CSC phenotypes, resulting in an increase in resistance. Here, we have developed a stochastic model to investigate the impact of plasticity on the efficacy of radiotherapy. The effect of plasticity on TCP is explored by applying the model to standard and hyper-fractionated schedules for a three week period of treatment as well as standard, hyper-fractionated, and accelerated hyper-fractionated schedules with an equal total dose of 30 Gy. Our results confirm that tumour control becomes increasingly difficult in the presence of plasticity as well as for the most resistant tumours. For the case with equal total dose, it is observed that increasing fractionation, at first enhances the probability of CSCs and tumour removal, but ultimately results in lower TCPS+P and TCPS. In addition, the combination of radiotherapy and targeted therapy (with increasing CSC differentiation) improves both the probability of CSC and tumour removal, in the absence of plasticity. However, in the presence of plasticity, the impact of combination therapy is not significant.

Nanoengineered Disruption of Heat Shock Protein 90 Targets Drug-Induced Resistance and Relieves Natural Killer Cell Suppression in Breast Cancer.

Drug-induced resistance, or tolerance, is an emerging yet poorly understood failure of anticancer therapy. The interplay between drug-tolerant cancer cells and innate immunity within the tumor, the consequence on tumor growth, and therapeutic strategies to address these challenges remain undescribed. Here, we elucidate the role of taxane-induced resistance on natural killer (NK) cell tumor immunity in triple-negative breast cancer (TNBC) and the design of spatiotemporally controlled nanomedicines, which boost therapeutic efficacy and invigorate "disabled" NK cells. Drug tolerance limited NK cell immune surveillance via drug-induced depletion of the NK-activating ligand receptor axis, NK group 2 member D, and MHC class I polypeptide-related sequence A, B. Systems biology supported by empirical evidence revealed the heat shock protein 90 (Hsp90) simultaneously controls immune surveillance and persistence of drug-treated tumor cells. On the basis of this evidence, we engineered a "chimeric" nanotherapeutic tool comprising taxanes and a cholesterol-tethered Hsp90 inhibitor, radicicol, which targets the tumor, reduces tolerance, and optimally reprimes NK cells via prolonged induction of NK-activating ligand receptors via temporal control of drug release in vitro and in vivo. A human ex vivo TNBC model confirmed the importance of NK cells in drug-induced death under pressure of clinically approved agents. These findings highlight a convergence between drug-induced resistance, the tumor immune contexture, and engineered approaches that consider the tumor and microenvironment to improve the success of combinatorial therapy. SIGNIFICANCE: This study uncovers a molecular mechanism linking drug-induced resistance and tumor immunity and provides novel engineered solutions that target these mechanisms in the tumor and improve immunity, thus mitigating off-target effects.