Resistance to targeted cancer therapies such as trastuzumab is a frequent clinical problem not solely because of insufficient expression of HER2 receptor but also because of the overriding activation states of cell signaling pathways. Systems biology approaches lend themselves to rapid in silico testing of factors, which may confer resistance to targeted therapies. Inthis study, we aimed to develop a new kinetic model that could be interrogated to predict resistance to receptor tyrosine kinase (RTK) inhibitor therapies and directly test predictions in vitro and in clinical samples. The new mathematical model included RTK inhibitor antibody binding, HER2/HER3 dimerization and inhibition, AKT/mitogen-activated protein kinase cross-talk, and the regulatory properties of PTEN. The model was parameterized using quantitative phosphoprotein expression data from cancer cell lines using reverse-phase protein microarrays. Quantitative PTEN protein expression was found to be the key determinant of resistance to anti-HER2 therapy in silico, which was predictive of unseen experiments in vitro using the PTEN inhibitor bp(V). When measured in cancer cell lines, PTEN expression predicts sensitivity to anti-HER2 therapy; furthermore, this quantitative measurement is more predictive of response (relative risk, 3.0; 95% confidence interval, 1.6–5.5; P < 0.0001) than other pathway components taken in isolation and when tested by multivariate analysis in a cohort of 122 breast cancers treated with trastuzumab. For the first time, a systems biology approach has successfully been used to stratify patients for personalized therapy in cancer and is further compelling evidence that PTEN, appropriately measured in the clinical setting, refines clinical decision making in patients treated with anti-HER2 therapies. [Cancer Res 2009;69(16):6713–20]

Photodynamic therapy (PDT) is a technique developed to treat the ever-increasing global incidence of cancer. This technique utilises singlet oxygen (1O2) generation via a laser excited photosensitiser (PS) to kill cancer cells. However, prolonged sensitivity to intensive light (6–8 weeks for lung cancer), relatively low tissue penetration by activating light (630 nm up to 4 mm), and the cost of PS administration can limit progressive PDT applications. The development of quantum-dot laser diodes emitting in the highest absorption region (1268 nm) of triplet oxygen (3O2) presents the possibility of inducing apoptosis in tumour cells through direct 3O2 → 1O2 transition. Here we demonstrate that a single laser pulse triggers dose-dependent 1O2 generation in both normal keratinocytes and tumour cells and show that tumour cells yield the highest 1O2 far beyond the initial laser pulse exposure. Our modelling and experimental results support the development of direct infrared (IR) laser-induced tumour treatment as a promising approach in tumour PDT.

Cells are constantly exposed to Reactive Oxygen Species (ROS) produced both endogenously to meet physiological requirements and from exogenous sources. While endogenous ROS are considered as important signalling molecules, high uncontrollable ROS are detrimental. It is unclear how cells can achieve a balance between maintaining physiological redox homeostasis and robustly activate the antioxidant system to remove exogenous ROS. We have utilised a Systems Biology approach to understand how this robust adaptive system fulfils homeostatic requirements of maintaining steady-state ROS and growth rate, while undergoing rapid readjustment under challenged conditions. Using a panel of human ovarian and normal cell lines, we experimentally quantified and established interrelationships between key elements of ROS homeostasis. The basal levels of NRF2 and KEAP1 were cell line specific and maintained in tight correlation with their growth rates and ROS. Furthermore, perturbation of this balance triggered cell specific kinetics of NRF2 nuclear-cytoplasmic relocalisation and sequestration of exogenous ROS. Our experimental data were employed to parameterise a mathematical model of the NRF2 pathway that elucidated key response mechanisms of redox regulation and showed that the dynamics of NRF2-H2O2 regulation defines a relationship between half-life, total and nuclear NRF2 level and endogenous H2O2 that is cell line specific.

High levels of variability in cancer-related cellular signalling networks and a lack of parameter identifiability in large-scale network models hamper translation of the results of modelling studies into the process of anti-cancer drug development. Recently global sensitivity analysis (GSA) has been recognised as a useful technique, capable of addressing the uncertainty of the model parameters and generating valid predictions on parametric sensitivities. Here we propose a novel implementation of model-based GSA specially designed to explore how multi-parametric network perturbations affect signal propagation through cancer-related networks. We use area-under-the-curve for time course of changes in phosphorylation of proteins as a characteristic for sensitivity analysis and rank network parameters with regard to their impact on the level of key cancer-related outputs, separating strong inhibitory from stimulatory effects. This allows interpretation of the results in terms which can incorporate the effects of potential anti-cancer drugs on targets and the associated biological markers of cancer. To illustrate the method we applied it to an ErbB signalling network model and explored the sensitivity profile of its key model readout, phosphorylated Akt, in the absence and presence of the ErbB2 inhibitor pertuzumab. The method successfully identified the parameters associated with elevation or suppression of Akt phosphorylation in the ErbB2/3 network. From analysis and comparison of the sensitivity profiles of pAkt in the absence and presence of targeted drugs we derived predictions of drug targets, cancer-related biomarkers and generated hypotheses for combinatorial therapy. Several key predictions have been confirmed in experiments using human ovarian carcinoma cell lines. We also compared GSA-derived predictions with the results of local sensitivity analysis and discuss the applicability of both methods. We propose that the developed GSA procedure can serve as a refining tool in combinatorial anti-cancer drug discovery.

Abstract: Computer simulation can be used to inform in vivo and in vitro experimentation, enabling rapid, low-cost hypothesis generation and directing experimental design in order to test those hypotheses. In this way, in silico models become a scientific instrument for investigation, and so should be developed to high standards, be carefully calibrated and their findings presented in such that they may be reproduced. Here, we outline a framework that supports developing simulations as scientific instruments, and we select cancer systems biology as an exemplar domain, with a particular focus on cellular signalling models. We consider the challenges of lack of data, incomplete knowledge and modelling in the context of a rapidly changing knowledge base. Our framework comprises a process to clearly separate scientific and engineering concerns in model and simulation development, and an argumentation approach to documenting models for rigorous way of recording assumptions and knowledge gaps. We propose interactive, dynamic visualisation tools to enable the biological community to interact with cellular signalling models directly for experimental design. There is a mismatch in scale between these cellular models and tissue structures that are affected by tumours, and bridging this gap requires substantial computational resource. We present concurrent programming as a technology to link scales without losing important details through model simplification. We discuss the value of combining this technology, interactive visualisation, argumentation and model separation to support development of multi-scale models that represent biologically plausible cells arranged in biologically plausible structures that model cell behaviour, interactions and response to therapeutic interventions.