Systems modeling of anti-apoptotic pathways in prostate cancer: psychological stress triggers a synergism pattern switch in drug combination therapy.

Prostate cancer patients often have increased levels of psychological stress or anxiety, but the molecular mechanisms underlying the interaction between psychological stress and prostate cancer as well as therapy resistance have been rarely studied and remain poorly understood. Recent reports show that stress inhibits apoptosis in prostate cancer cells via epinephrine/beta2 adrenergic receptor/PKA/BAD pathway. In this study, we used experimental data on the signaling pathways that control BAD phosphorylation to build a dynamic network model of apoptosis regulation in prostate cancer cells. We then compared the predictive power of two different models with or without the role of Mcl-1, which justified the role of Mcl-1 stabilization in anti-apoptotic effects of emotional stress. Based on the selected model, we examined and quantitatively evaluated the induction of apoptosis by drug combination therapies. We predicted that the combination of PI3K inhibitor LY294002 and inhibition of BAD phosphorylation at S112 would produce the best synergistic effect among 8 interventions examined. Experimental validation confirmed the effectiveness of our predictive model. Moreover, we found that epinephrine signaling changes the synergism pattern and decreases efficacy of combination therapy. The molecular mechanisms responsible for therapeutic resistance and the switch in synergism were explored by analyzing a network model of signaling pathways affected by psychological stress. These results provide insights into the mechanisms of psychological stress signaling in therapy-resistant cancer, and indicate the potential benefit of reducing psychological stress in designing more effective therapies for prostate cancer patients.

A novel missense-mutation-related feature extraction scheme for 'driver' mutation identification.

MOTIVATION: It becomes widely accepted that human cancer is a disease involving dynamic changes in the genome and that the missense mutations constitute the bulk of human genetic variations. A multitude of computational algorithms, especially the machine learning-based ones, has consequently been proposed to distinguish missense changes that contribute to the cancer progression ('driver' mutation) from those that do not ('passenger' mutation). However, the existing methods have multifaceted shortcomings, in the sense that they either adopt incomplete feature space or depend on protein structural databases which are usually far from integrated. RESULTS: In this article, we investigated multiple aspects of a missense mutation and identified a novel feature space that well distinguishes cancer-associated driver mutations from passenger ones. An index (DX score) was proposed to evaluate the discriminating capability of each feature, and a subset of these features which ranks top was selected to build the SVM classifier. Cross-validation showed that the classifier trained on our selected features significantly outperforms the existing ones both in precision and robustness. We applied our method to several datasets of missense mutations culled from published database and literature and obtained more reasonable results than previous studies. AVAILABILITY: The software is available online at http://www.methodisthealth.com/software and https://sites.google.com/site/drivermutationidentification/. CONTACT: xzhou@tmhs.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A computational model for compressed sensing RNAi cellular screening.

BACKGROUND: RNA interference (RNAi) becomes an increasingly important and effective genetic tool to study the function of target genes by suppressing specific genes of interest. This system approach helps identify signaling pathways and cellular phase types by tracking intensity and/or morphological changes of cells. The traditional RNAi screening scheme, in which one siRNA is designed to knockdown one specific mRNA target, needs a large library of siRNAs and turns out to be time-consuming and expensive. RESULTS: In this paper, we propose a conceptual model, called compressed sensing RNAi (csRNAi), which employs a unique combination of group of small interfering RNAs (siRNAs) to knockdown a much larger size of genes. This strategy is based on the fact that one gene can be partially bound with several small interfering RNAs (siRNAs) and conversely, one siRNA can bind to a few genes with distinct binding affinity. This model constructs a multi-to-multi correspondence between siRNAs and their targets, with siRNAs much fewer than mRNA targets, compared with the conventional scheme. Mathematically this problem involves an underdetermined system of equations (linear or nonlinear), which is ill-posed in general. However, the recently developed compressed sensing (CS) theory can solve this problem. We present a mathematical model to describe the csRNAi system based on both CS theory and biological concerns. To build this model, we first search nucleotide motifs in a target gene set. Then we propose a machine learning based method to find the effective siRNAs with novel features, such as image features and speech features to describe an siRNA sequence. Numerical simulations show that we can reduce the siRNA library to one third of that in the conventional scheme. In addition, the features to describe siRNAs outperform the existing ones substantially. CONCLUSIONS: This csRNAi system is very promising in saving both time and cost for large-scale RNAi screening experiments which may benefit the biological research with respect to cellular processes and pathways.

Multi-scale agent-based brain cancer modeling and prediction of TKI treatment response: incorporating EGFR signaling pathway and angiogenesis.

BACKGROUND: The epidermal growth factor receptor (EGFR) signaling pathway and angiogenesis in brain cancer act as an engine for tumor initiation, expansion and response to therapy. Since the existing literature does not have any models that investigate the impact of both angiogenesis and molecular signaling pathways on treatment, we propose a novel multi-scale, agent-based computational model that includes both angiogenesis and EGFR modules to study the response of brain cancer under tyrosine kinase inhibitors (TKIs) treatment. RESULTS: The novel angiogenesis module integrated into the agent-based tumor model is based on a set of reaction-diffusion equations that describe the spatio-temporal evolution of the distributions of micro-environmental factors such as glucose, oxygen, TGFα, VEGF and fibronectin. These molecular species regulate tumor growth during angiogenesis. Each tumor cell is equipped with an EGFR signaling pathway linked to a cell-cycle pathway to determine its phenotype. EGFR TKIs are delivered through the blood vessels of tumor microvasculature and the response to treatment is studied. CONCLUSIONS: Our simulations demonstrated that entire tumor growth profile is a collective behaviour of cells regulated by the EGFR signaling pathway and the cell cycle. We also found that angiogenesis has a dual effect under TKI treatment: on one hand, through neo-vasculature TKIs are delivered to decrease tumor invasion; on the other hand, the neo-vasculature can transport glucose and oxygen to tumor cells to maintain their metabolism, which results in an increase of cell survival rate in the late simulation stages.

A Spatio-Temporal Model of Macrophage-Mediated Drug Resistance in Glioma Immunotherapy.

The emergence of drug resistance is often an inevitable obstacle that limits the long-term effectiveness of clinical cancer chemotherapeutics. Although various forms of cancer cell-intrinsic mechanisms of drug resistance have been experimentally revealed, the role and the underlying mechanism of tumor microenvironment in driving the development of acquired drug resistance remain elusive, which significantly impedes effective clinical cancer treatment. Recent experimental studies have revealed a macrophage-mediated drug resistance mechanism in which the tumor microenvironment undergoes adaptation in response to macrophage-targeted colony-stimulating factor-1 receptor (CSF1R) inhibition therapy in gliomas. In this study, we developed a spatio-temporal model to quantitatively describe the interplay between glioma cells and CSF1R inhibitor-targeted macrophages through CSF1 and IGF1 pathways. Our model was used to investigate the evolutionary kinetics of the tumor regrowth and the associated dynamic adaptation of the tumor microenvironment in response to the CSF1R inhibitor treatment. The simulation result obtained using this model was in agreement with the experimental data. The sensitivity analysis revealed the key parameters involved in the model, and their potential impacts on the model behavior were examined. Moreover, we demonstrated that the drug resistance is dose-dependent. In addition, we quantitatively evaluated the effects of combined CSFR inhibition and IGF1 receptor (IGF1R) inhibition with the goal of designing more effective therapies for gliomas. Our study provides quantitative and mechanistic insights into the microenvironmental adaptation mechanisms that operate during macrophage-targeted immunotherapy and has implications for drug dose optimization and the design of more effective combination therapies. Mol Cancer Ther; 17(4); 814-24. ©2018 AACR.

Mathematical Modeling of Therapy-induced Cancer Drug Resistance: Connecting Cancer Mechanisms to Population Survival Rates.

Drug resistance significantly limits the long-term effectiveness of targeted therapeutics for cancer patients. Recent experimental studies have demonstrated that cancer cell heterogeneity and microenvironment adaptations to targeted therapy play important roles in promoting the rapid acquisition of drug resistance and in increasing cancer metastasis. The systematic development of effective therapeutics to overcome drug resistance mechanisms poses a major challenge. In this study, we used a modeling approach to connect cellular mechanisms underlying cancer drug resistance to population-level patient survival. To predict progression-free survival in cancer patients with metastatic melanoma, we developed a set of stochastic differential equations to describe the dynamics of heterogeneous cell populations while taking into account micro-environment adaptations. Clinical data on survival and circulating tumor cell DNA (ctDNA) concentrations were used to confirm the effectiveness of our model. Moreover, our model predicted distinct patterns of dose-dependent synergy when evaluating a combination of BRAF and MEK inhibitors versus a combination of BRAF and PI3K inhibitors. These predictions were consistent with the findings in previously reported studies. The impact of the drug metabolism rate on patient survival was also discussed. The proposed model might facilitate the quantitative evaluation and optimization of combination therapeutics and cancer clinical trial design.

Modeling of signaling crosstalk-mediated drug resistance and its implications on drug combination.

The efficacy of pharmacological perturbation to the signaling transduction network depends on the network topology. However, whether and how signaling dynamics mediated by crosstalk contributes to the drug resistance are not fully understood and remain to be systematically explored. In this study, motivated by a realistic signaling network linked by crosstalk between EGF/EGFR/Ras/MEK/ERK pathway and HGF/HGFR/PI3K/AKT pathway, we develop kinetic models for several small networks with typical crosstalk modules to investigate the role of the architecture of crosstalk in inducing drug resistance. Our results demonstrate that crosstalk inhibition diminishes the response of signaling output to the external stimuli. Moreover, we show that signaling crosstalk affects the relative sensitivity of drugs, and some types of crosstalk modules that could yield resistance to the targeted drugs were identified. Furthermore, we quantitatively evaluate the relative efficacy and synergism of drug combinations. For the modules that are resistant to the targeted drug, we identify drug targets that can not only increase the relative drug efficacy but also act synergistically. In addition, we analyze the role of the strength of crosstalk in switching a module between drug-sensitive and drug-resistant. Our study provides mechanistic insights into the signaling crosstalk-mediated mechanisms of drug resistance and provides implications for the design of synergistic drug combinations to reduce drug resistance.

Genome-wide mutational spectra analysis reveals significant cancer-specific heterogeneity.

Cancer is widely recognized as a genetic disease in which somatic mutations are sequentially accumulated to drive tumor progression. Although genomic landscape studies are informative for individual cancer types, a comprehensive comparative study of tumorigenic mutations across cancer types based on integrative data sources is still a pressing need. We systematically analyzed ~10(6) non-synonymous mutations extracted from COSMIC, involving ~8000 genome-wide screened samples across 23 major human cancers at both the amino acid and gene levels. Our analysis identified cancer-specific heterogeneity that traditional nucleotide variation analysis alone usually overlooked. Particularly, the amino acid arginine (R) turns out to be the most favorable target of amino acid alteration in most cancer types studied (P < 10(-9), binomial test), reflecting its important role in cellular physiology. The tumor suppressor gene TP53 is mutated exclusively with the HYDIN, KRAS, and PTEN genes in large intestine, lung, and endometrial cancers respectively, indicating that TP53 takes part in different signaling pathways in different cancers. While some of our analyses corroborated previous observations, others indicated relevant candidates with high priority for further experimental validation. Our findings have many ramifications in understanding the etiology of cancer and the underlying molecular mechanisms in particular cancers.

Systems biology approach to studying proliferation-dependent prognostic subnetworks in breast cancer.

Tumor proliferative capacity is a major biological correlate of breast tumor metastatic potential. In this paper, we developed a systems approach to investigate associations among gene expression patterns, representative protein-protein interactions, and the potential for clinical metastases, to uncover novel survival-related subnetwork signatures as a function of tumor proliferative potential. Based on the statistical associations between gene expression patterns and patient outcomes, we identified three groups of survival prognostic subnetwork signatures (SPNs) corresponding to three proliferation levels. We discovered 8 SPNs in the high proliferation group, 8 SPNs in the intermediate proliferation group, and 6 SPNs in the low proliferation group. We observed little overlap of SPNs between the three proliferation groups. The enrichment analysis revealed that most SPNs were enriched in distinct signaling pathways and biological processes. The SPNs were validated on other cohorts of patients, and delivered high accuracy in the classification of metastatic vs non-metastatic breast tumors. Our findings indicate that certain biological networks underlying breast cancer metastasis differ in a proliferation-dependent manner. These networks, in combination, may form the basis of highly accurate prognostic classification models and may have clinical utility in guiding therapeutic options for patients.

In silico study on multidrug resistance conferred by I223R/H275Y double mutant neuraminidase.

It was recently reported that an I223R/H275Y double mutant of neuraminidase (NA) creates a multidrug-resistant form of the pandemic influenza A (H1N1) virus. However, a comprehensive understanding of the molecular mechanisms is still lacking. We conducted a systematic in silico study to explore the structural basis underlying this multidrug resistance. By molecular docking analyses and molecular dynamics (MD) simulations, we compared various biochemical and biophysical properties of the wild type, the I223R single mutant and the I223R/H275Y double mutant NA with two inhibitors, zanamivir (ZMR) and oseltamivir (G39). The binding free energy of oseltamivir with all types of NA was substantially lower than its zanamivir counterpart. On the other hand, the binding free energy of each inhibitor with wild type NA was generally higher than that with mutant NAs. MD simulation outcomes exemplify distinct patterns for oseltamivir and zanamivir with all types of NA. In particular, the stronger resistance of the double mutant NA relative to the wild and single mutant types can be ascribed to the overall looser but locally more compact structure of the former. Specifically, as a whole the double mutant NA adapts to the larger gyration radius and greater distance between charged atom groups, which is contrary to the pattern in the local binding site region. The enhanced resistance of all types of NA to oseltamivir rather than zanamivir might be accounted for similarly. We expect these findings to provide significant insights into improving inhibitors for the multidrug-resistant neuraminidase of H1N1 influenza viruses.

Modeling vascularized bone regeneration within a porous biodegradable CaP scaffold loaded with growth factors.

Osteogenetic microenvironment is a complex constitution in which extracellular matrix (ECM) molecules, stem cells and growth factors each interact to direct the coordinate regulation of bone tissue development. Importantly, angiogenesis improvement and revascularization are critical for osteogenesis during bone tissue regeneration processes. In this study, we developed a three-dimensional (3D) multi-scale system model to study cell response to growth factors released from a 3D biodegradable porous calcium phosphate (CaP) scaffold. Our model reconstructed the 3D bone regeneration system and examined the effects of pore size and porosity on bone formation and angiogenesis. The results suggested that scaffold porosity played a more dominant role in affecting bone formation and angiogenesis compared with pore size, while the pore size could be controlled to tailor the growth factor release rate and release fraction. Furthermore, a combination of gradient VEGF with BMP2 and Wnt released from the multi-layer scaffold promoted angiogenesis and bone formation more readily than single growth factors. These results demonstrated that the developed model can be potentially applied to predict vascularized bone regeneration with specific scaffold and growth factors.

A 3D multiscale model of cancer stem cell in tumor development.

BACKGROUND: Recent reports indicate that a subgroup of tumor cells named cancer stem cells (CSCs) or tumor initiating cells (TICs) are responsible for tumor initiation, growth and drug resistance. This subgroup of tumor cells has self-renewal capacity and could differentiate into heterogeneous tumor cell populations through asymmetric proliferation. The idea of CSC provides informative insights into tumor initiation, metastasis and treatment. However, the underlying mechanisms of CSCs regulating tumor behaviors are unclear due to the complex cancer system. To study the functions of CSCs in the complex tumor system, a few mathematical modeling studies have been proposed. Whereas, the effect of microenvironment (mE) factors, the behaviors of CSCs, progenitor tumor cells (PCs) and differentiated tumor cells (TCs), and the impact of CSC fraction and signaling heterogeneity, are not adequately explored yet. METHODS: In this study, a novel 3D multi-scale mathematical modeling is proposed to investigate the behaviors of CSCsin tumor progressions. The model integrates CSCs, PCs, and TCs together with a few essential mE factors. With this model, we simulated and investigated the tumor development and drug response under different CSC content and heterogeneity. RESULTS: The simulation results shown that the fraction of CSCs plays a critical role in driving the tumor progression and drug resistance. It is also showed that the pure chemo-drug treatment was not a successful treatment, as it resulted in a significant increase of the CSC fraction. It further shown that the self-renew heterogeneity of the initial CSC population is a cause of the heterogeneity of the derived tumors in terms of the CSC fraction and response to drug treatments. CONCLUSIONS: The proposed 3D multi-scale model provides a new tool for investigating the behaviors of CSC in CSC-initiated tumors, which enables scientists to investigate and generate testable hypotheses about CSCs in tumor development and drug response under different microenvironments and drug perturbations.

Cytokine combination therapy prediction for bone remodeling in tissue engineering based on the intracellular signaling pathway.

The long-term performance of tissue-engineered bone grafts is determined by a dynamic balance between bone regeneration and resorption. We proposed using embedded cytokine slow-releasing hydrogels to tune this balance toward a desirable final bone density. In this study we established a systems biology model, and quantitatively explored the combinatorial effects of delivered cytokines from hydrogels on final bone density. We hypothesized that: 1) bone regeneration was driven by transcription factors Runx2 and Osterix, which responded to released cytokines, such as Wnt, BMP2, and TGFß, drove the development of osteoblast lineage, and contributed to bone mass generation; and 2) the osteoclast lineage, on the other hand, governed the bone resorption, and communications between these two lineages determined the dynamics of bone remodeling. In our model, Intracellular signaling pathways were represented by ordinary differential equations, while the intercellular communications and cellular population dynamics were modeled by stochastic differential equations. Effects of synergistic cytokine combinations were evaluated by Loewe index and Bliss index. Simulation results revealed that the Wnt/BMP2 combinations released from hydrogels showed best control of bone regeneration and synergistic effects, and suggested optimal dose ratios of given cytokine combinations released from hydrogels to most efficiently control the long-term bone remodeling. We revealed the characteristics of cytokine combinations of Wnt/BMP2 which could be used to guide the design of in vivo bone scaffolds and the clinical treatment of some diseases such as osteoporosis.