Mathematical modeling of the cancer cell’s control circuitry : paving the way to individualized therapeutic strategies

Cancer is a disease of signal transduction in which the dysregulation of the network of intracellular and extracellular signaling cascades is sufficient to thwart the cells finely-tuned biochemical control mechanisms. A keen interest in the mathematical modeling of cell signaling networks and the regulation of signal transduction has emerged in recent years, and has produced a glimmer of insight into the sophisticated feedback control and network regulation operating within cells. In this review, we present an overview of published theoretical studies on the control aspects of signal transduction, emphasizing the role and importance of mechanisms such as ‘ultrasensitivity’ and feedback loops. We emphasize that these exquisite and often subtle control strategies represent the key to orchestrating ‘simple’ signaling behaviors within the complex intracellular network, while regulating the trade-off between sensitivity and robustness to internal and external perturbations. Through a consideration of these apparent paradoxes, we explore how the basic homeostasis of the intracellular signaling network, in the face of carcinogenesis, can lead to neoplastic progression rather than cell death. A simple mathematical model is presented, furnishing a vivid illustration of how ‘control-oriented’ models of the deranged signaling networks in cancer cells may enucleate improved treatment strategies, including patient-tailored combination therapies, with the potential for reduced toxicity and more robust and potent antitumor activity.

Development of reverse phase protein microarrays for clinical applications and patient-tailored therapy

While genomics provide important information about the somatic genetic changes, and RNA transcript profiling can reveal important expression changes that correlate with outcome and response to therapy, it is the proteins that do the work in the cell. At a functional level, derangements within the proteome, driven by post-translational and epigenetic modifications, such as phosphorylation, is the cause of a vast majority of human diseases. Cancer, for instance, is a manifestation of deranged cellular protein molecular networks and cell signaling pathways that are based on genetic changes at the DNA level. Importantly, the protein pathways contain the drug targets in signaling networks that govern overall cellular survival, proliferation, invasion and cell death. Consequently, the promise of proteomics resides in the ability to extend analysis beyond correlation to causality. A critical gap in the information knowledge base of molecular profiling is an understanding of the ongoing activity of protein signaling in human tissue: what is activated and “in use” within the human body at any given point in time. To address this gap, we have invented a new technology, called reverse phase protein microarrays, that can generate a functional read-out of cell signaling networks or pathways for an individual patient obtained directly from a biopsy specimen. This “wiring diagram” can serve as the basis for both, selection of a therapy and patient stratification.

Phosphoprotein pathway mapping: Akt/Mammalian target of rapamycin activation is negatively associated with childhood rhabdomyosarcoma survival

Mapping of protein signaling networks within tumors can identify new targets for therapy and provide a means to stratify patients for individualized therapy. Despite advances in combination chemotherapy, the overall survival for childhood rhabdomyosarcoma remains ∼60%. A critical goal is to identify functionally important protein signaling defects associated with treatment failure for the 40% nonresponder cohort. Here, we show, by phosphoproteomic network analysis of microdissected tumor cells, that interlinked components of the Akt/mammalian target of rapamycin (mTOR) pathway exhibited increased levels of phosphorylation for tumors of patients with short-term survival. Specimens (n = 59) were obtained from the Children's Oncology Group Intergroup Rhabdomyosarcoma Study (IRS) IV, D9502 and D9803, with 12-year follow-up. High phosphorylation levels were associated with poor overall and poor disease-free survival: Akt Ser473 (overall survival P < 0.001, recurrence-free survival P < 0.0009), 4EBP1 Thr37/46 (overall survival P < 0.0110, recurrence-free survival P < 0.0106), eIF4G Ser1108 (overall survival P < 0.0017, recurrence-free survival P < 0.0072), and p70S6 Thr389 (overall survival P < 0.0085, recurrence-free survival P < 0.0296). Moreover, the findings support an altered interrelationship between the insulin receptor substrate (IRS-1) and Akt/mTOR pathway proteins (P < 0.0027) for tumors from patients with poor survival. The functional significance of this pathway was tested using CCI-779 in a mouse xenograft model. CCI-779 suppressed phosphorylation of mTOR downstream proteins and greatly reduced the growth of two different rhabdomyosarcoma (RD embryonal P = 0.00008; Rh30 alveolar P = 0.0002) cell lines compared with controls. These results suggest that phosphoprotein mapping of the Akt/mTOR pathway should be studied further as a means to select patients to receive mTOR/IRS pathway inhibitors before administration of chemotherapy.

A control theoretic paradigm for cell signaling networks : a simple complexity for a sensitive robustness

The fields of molecular biology and cell biology are being flooded with complex genomic and proteomic datasets of large dimensions. We now recognize that each molecule in the cell and tissue can no longer be viewed as an isolated entity. Instead, each molecule must be considered as one member of an interacting network. Consequently, there is an urgent need for mathematical models to understand the behavior of cell signaling networks in health and in disease.

Innate transcriptional networks activated in bladder in response to uropathogenic Escherichia coli drive diverse biological pathways and rapid synthesis of IL-10 for defense against bacterial urinary tract infection

Early transcriptional activation events that occur in bladder immediately following bacterial urinary tract infection (UTI) are not well defined. In this study, we describe the whole bladder transcriptome of uropathogenic Escherichia coli (UPEC) cystitis in mice using genome-wide expression profiling to define the transcriptome of innate immune activation stemming from UPEC colonization of the bladder. Bladder RNA from female C57BL/6 mice, analyzed using 1.0 ST-Affymetrix microarrays, revealed extensive activation of diverse sets of innate immune response genes, including those that encode multiple IL-family members, receptors, metabolic regulators, MAPK activators, and lymphocyte signaling molecules. These were among 1564 genes differentially regulated at 2 h postinfection, highlighting a rapid and broad innate immune response to bladder colonization. Integrative systems-level analyses using InnateDB (http://www.innatedb.com) bioinformatics and ingenuity pathway analysis identified multiple distinct biological pathways in the bladder transcriptome with extensive involvement of lymphocyte signaling, cell cycle alterations, cytoskeletal, and metabolic changes. A key regulator of IL activity identified in the transcriptome was IL-10, which was analyzed functionally to reveal marked exacerbation of cystitis in IL-10–deficient mice. Studies of clinical UTI revealed significantly elevated urinary IL-10 in patients with UPEC cystitis, indicating a role for IL-10 in the innate response to human UTI. The whole bladder transcriptome presented in this work provides new insight into the diversity of innate factors that determine UTI on a genome-wide scale and will be valuable for further data mining. Identification of protective roles for other elements in the transcriptome will provide critical new insight into the complex cascade of events that underpin UTI.

Clinical Escherichia coli isolates utilize alpha-hemolysin to inhibit in vitro epithelial cytokine production

Uropathogenic Escherichia coli is the primary cause of urinary tract infections, which affects over 60% of women during their lifetime. UPEC exhibits a number of virulence traits that facilitate colonization of the bladder, including inhibition of cytokine production by bladder epithelial cells. The goal of this study was to identify the mechanism of this inhibition. We observed that cytokine suppression was associated with rapid cytotoxicity toward epithelial cells. We found that cytotoxicity, cytokine suppression and alpha-hemolysin production were all tightly linked in clinical isolates. We screened a UPEC fosmid library and identified clones that gained the cytotoxicity and cytokine-suppression phenotypes. Both clones contained fosmids encoding a PAI II(J96)-like domain and expressed the alpha-hemolysin (hlyA) encoded therein. Mutation of the fosmid-encoded hly operon abolished cytotoxicity and cytokine suppression. Similarly, mutation of the chromosomal hlyCABD operon of UPEC isolate F11 also abolished these phenotypes, and they could be restored by introducing the PAI II(J96)-like domain-encoding fosmid. We also examined the role of alpha-hemolysin in cytokine production both in the murine UTI model as well as patient specimens. We conclude that E. coli utilizes alpha-hemolysin to inhibit epithelial cytokine production in vitro. Its contribution to inflammation during infection requires further study.

Migraine and glutamate: Modulators of glutamatergic signalling as potential treatments of neuropsychiatric disorders

Migraine is a complex neurological disorder with a well-documented genetic basis. Migraine is a product of allelic variation in genes of neurological, vascular and hormonal origin interacting with environmental triggers. Presentation can include attacks of head pain with symptoms of nausea, emesis, photophobia, phonophobia, and occasionally, visual sensory disturbances, known as aura. Migraine pain is difficult to ignore, associated with a deep sense of malaise and manifests as a throbbing, pulsatile headache, localized to one side of the head that intensifies with physical activity and that can last from 4-72 hours. Migraine is diagnosed according to criteria developed by the International Headache Society (IHS) and is subdivided into two main types based on the occurrence of aura symptoms that may be present in the early stages of the headache: migraine with aura (MA) and migraine without aura (MO). The majority (about 70%) of migraineurs are diagnosed with the MO subtype whilst the remaining 30% experience MA accompanied by neurological symptoms that manifest as fully reversible, visual, sensory and/or dysphasic speech disturbances in conjunction with their headache. Glutamate is the primary excitatory neurotransmitter in the central nervous system (CNS) and over-excitation of glutamate receptors is regarded as a contributing factor, through various mechanisms, to the pathology of migraine. In this chapter we present an overview of the pathophysiology and co-morbidity of migraine with other psychiatric disorders and discuss the role of the glutamatergic system in migraine, its molecular components as potential drug targets, in addition to the current treatments and progress of modulators of glutamatergic signaling.

A fiber Bragg grating sensor system for train axle counting

Railway signaling facilitates two main functions, namely, train detection and train control, in order to maintain safe separations among the trains. Track circuits are the most commonly used train detection means with the simple open/close circuit principles; and subsequent adoption of axle counters further allows the detection of trains under adverse track conditions. However, with electrification and power electronics traction drive systems, aggravated by the electromagnetic interference in the vicinity of the signaling system, railway engineers often find unstable or even faulty operations of track circuits and axle counting systems, which inevitably jeopardizes the safe operation of trains. A new means of train detection, which is completely free from electromagnetic interference, is therefore required for the modern railway signaling system. This paper presents a novel optical fiber sensor signaling system. The sensor operation, field setup, axle detection solution set, and test results of an installation in a trial system on a busy suburban railway line are given.

Generalized binomial Tau-leap method for biochemical kinetics incorporating both delay and intrinsic noise

The delay stochastic simulation algorithm (DSSA) by Barrio et al. [Plos Comput. Biol.2, 117–E (2006)] was developed to simulate delayed processes in cell biology in the presence of intrinsic noise, that is, when there are small-to-moderate numbers of certain key molecules present in a chemical reaction system. These delayed processes can faithfully represent complex interactions and mechanisms that imply a number of spatiotemporal processes often not explicitly modeled such as transcription and translation, basic in the modeling of cell signaling pathways. However, for systems with widely varying reaction rate constants or large numbers of molecules, the simulation time steps of both the stochastic simulation algorithm (SSA) and the DSSA can become very small causing considerable computational overheads. In order to overcome the limit of small step sizes, various τ-leap strategies have been suggested for improving computational performance of the SSA. In this paper, we present a binomial τ- DSSA method that extends the τ-leap idea to the delay setting and avoids drawing insufficient numbers of reactions, a common shortcoming of existing binomial τ-leap methods that becomes evident when dealing with complex chemical interactions. The resulting inaccuracies are most evident in the delayed case, even when considering reaction products as potential reactants within the same time step in which they are produced. Moreover, we extend the framework to account for multicellular systems with different degrees of intercellular communication. We apply these ideas to two important genetic regulatory models, namely, the hes1 gene, implicated as a molecular clock, and a Her1/Her 7 model for coupled oscillating cells.

Exome sequencing identifies recurrent somatic RAC1 mutations in melanoma

We characterized the mutational landscape of melanoma, the form of skin cancer with the highest mortality rate, by sequencing the exomes of 147 melanomas. Sun-exposed melanomas had markedly more ultraviolet (UV)-like C>T somatic mutations compared to sun-shielded acral, mucosal and uveal melanomas. Among the newly identified cancer genes was PPP6C, encoding a serine/threonine phosphatase, which harbored mutations that clustered in the active site in 12% of sun-exposed melanomas, exclusively in tumors with mutations in BRAF or NRAS. Notably, we identified a recurrent UV-signature, an activating mutation in RAC1 in 9.2% of sun-exposed melanomas. This activating mutation, the third most frequent in our cohort of sun-exposed melanoma after those of BRAF and NRAS, changes Pro29 to serine (RAC1P29S) in the highly conserved switch I domain. Crystal structures, and biochemical and functional studies of RAC1P29S showed that the alteration releases the conformational restraint conferred by the conserved proline, causes an increased binding of the protein to downstream effectors, and promotes melanocyte proliferation and migration. These findings raise the possibility that pharmacological inhibition of downstream effectors of RAC1 signaling could be of therapeutic benefit.

COX-derived prostanoid pathways in gastrointestinal cancer development and progression : novel targets for prevention and intervention

Arachidonic acid metabolism through cyclooxygenase (COX) pathways leads to the generation of biologically active eicosanoids. Eicosanoid expression levels vary during development and progression of gastrointestinal (GI) malignancies. COX-2 is the major COX-isoform responsible for G.I. cancer development/progression. COX-2 expression increases during progression from a normal to cancerous state. Evidence from observational studies has demonstrated that chronic NSAID use reduces the risk of cancer development, while both incidence and risk of death due to G.I. cancers were significantly reduced by daily aspirin intake. A number of randomized controlled trials (APC trial, Prevention of Sporadic Adenomatous Polyps trial, APPROVe trial) have also shown a significant protective effect in patients receiving selective COX-2 inhibitors. However, chronic use of selective COX-2 inhibitors at high doses was associated with increased cardiovascular risk, while NSAIDs have also been associated with increased risk. More recently, downstream effectors of COX-signaling have been investigated in cancer development/progression. PGE 2, which binds to both EP and PPAR receptors, is the major prostanoid implicated in the carcinogenesis of G.I. cancers. The role of TXA 2 in G.I. cancers has also been examined, although further studies are required to uncover its role in carcinogenesis. Other prostanoids investigated include PGD 2 and its metabolite 15d-PGJ2, PGF 1α and PGI 2. Targeting these prostanoids in G.I. cancers has the promise of avoiding cardiovascular toxicity associated with chronic selective COX-2 inhibition, while maintaining anti-tumor reactivity.A progressive sequence from normal to pre-malignant to a malignant state has been identified in G.I. cancers. In this review, we will discuss the role of the COX-derived prostanoids in G.I. cancer development and progression. Targeting these downstream prostanoids for chemoprevention and/or treatment of G.I. cancers will also be discussed. Finally, we will highlight the latest pre-clinical technologies as well as avenues for future investigation in this highly topical research field. © 2011 Elsevier B.V.

A soluble activin Type IIA receptor induces bone formation and improves skeletal integrity

Diseases that affect the regulation of bone turnover can lead to skeletal fragility and increased fracture risk. Members of the TGF-superfamily have been shown to be involved in the regulation of bone mass. Activin A, a TGF-� signaling ligand, is present at high levels in bone and may play a role in the regulation of bone metabolism. Here we demonstrate that pharmacological blockade of ligand signaling through the high affinity receptor for activin, type II activin receptor (ActRIIA), by administration of the soluble extracellular domain of ActRIIA fused to a murine IgG2a-Fc, increases bone formation, bone mass, and bone strength in normal mice and in ovariectomized mice with established bone loss. These observations support the development of this pharmacological strategy for the treatment of diseases with skeletal fragility.

Heregulin regulates the actin cytoskeleton and promotes invasive properties in breast cancer cell lines

The metastatic process requires changes in tumor cell adhesion properties, cell motility and remodeling of the extracellular matrix. The erbB2 proto-oncogene is overexpressed in approximately 30% of breast cancers and is a major prognostic parameter when present in invasive disease. A ligand for the erbB2 receptor has not yet been identified but it can be activated by heterodimerization with heregulin (HRG)-stimulated erbB3 and erbB4 receptors. The HRGs are a family of polypeptide growth factors that have been shown to play a role in embryogenesis, tumor formation, growth and differentiation of breast cancer cells. The erbB3 and erbB4 receptors are involved in transregulation of erbB2 signaling. The work presented here suggests biological roles for HRG including regulation of the actin cytoskeleton and induction of motility and invasion in breast cancer cells. HRG-expressing breast cancer cell lines are characterized by low erbB receptor levels and a high invasive and metastatic index, while those which overexpress erbB2 demonstrate minimal invasive potential in vitro and are non-tumorigenic in vivo. Treatment of the highly tumorigenic and metastatic HRG-expressing breast cancer cell line MDA-MB-231 with an HRG-neutralizing antibody significantly inhibited proliferation in culture and motility in the Boyden chamber assay. Addition of exogenous HRG to non-invasive erbB2 overexpressing cells (SKBr-3) at low concentrations induced formation of pseudopodia, enhanced phagocytic activity and increased chemomigration and invasion in the Boyden chamber assay. The specificity of the chemomigration response to HRG is demonstrated by inhibition with the anti-HRG neutralizing antibody. These results suggest that either HRG can act as an autocrine or paracrine ligand to promote the invasive behavior of breast cancer cells in vitro or thus may enhance the metastatic process in vivo.

A computational model for BMP movement in sea urchin embryos

Bone morphogen proteins (BMPs) are distributed along a dorsal-ventral (DV) gradient in many developing embryos. The spatial distribution of this signaling ligand is critical for correct DV axis specification. In various species, BMP expression is spatially localized, and BMP gradient formation relies on BMP transport, which in turn requires interactions with the extracellular proteins Short gastrulation/Chordin (Chd) and Twisted gastrulation (Tsg). These binding interactions promote BMP movement and concomitantly inhibit BMP signaling. The protease Tolloid (Tld) cleaves Chd, which releases BMP from the complex and permits it to bind the BMP receptor and signal. In sea urchin embryos, BMP is produced in the ventral ectoderm, but signals in the dorsal ectoderm. The transport of BMP from the ventral ectoderm to the dorsal ectoderm in sea urchin embryos is not understood. Therefore, using information from a series of experiments, we adapt the mathematical model of Mizutani et al. (2005) and embed it as the reaction part of a one-dimensional reaction–diffusion model. We use it to study aspects of this transport process in sea urchin embryos. We demonstrate that the receptor-bound BMP concentration exhibits dorsally centered peaks of the same type as those observed experimentally when the ternary transport complex (Chd-Tsg-BMP) forms relatively quickly and BMP receptor binding is relatively slow. Similarly, dorsally centered peaks are created when the diffusivities of BMP, Chd, and Chd-Tsg are relatively low and that of Chd-Tsg-BMP is relatively high, and the model dynamics also suggest that Tld is a principal regulator of the system. At the end of this paper, we briefly compare the observed dynamics in the sea urchin model to a version that applies to the fly embryo, and we find that the same conditions can account for BMP transport in the two types of embryos only if Tld levels are reduced in sea urchin compared to fly.

Transcriptional pathway signatures predict MEK addiction and response to selumetinib (AZD6244)

Selumetinib (AZD6244, ARRY-142886) is a selective, non-ATP-competitive inhibitor of mitogen-activated protein/extracellular signal-regulated kinase kinase (MEK)-1/2. The range of antitumor activity seen preclinically and in patients highlights the importance of identifying determinants of response to this drug. In large tumor cell panels of diverse lineage, we show that MEK inhibitor response does not have an absolute correlation with mutational or phospho-protein markers of BRAF/MEK, RAS, or phosphoinositide 3-kinase (PI3K) activity. We aimed to enhance predictivity by measuring pathway output through coregulated gene networks displaying differential mRNA expression exclusive to resistant cell subsets and correlated to mutational or dynamic pathway activity. We discovered an 18-gene signature enabling measurement of MEK functional output independent of tumor genotype. Where the MEK pathway is activated but the cells remain resistant to selumetinib, we identified a 13-gene signature that implicates the existence of compensatory signaling from RAS effectors other than PI3K. The ability of these signatures to stratify samples according to functional activation of MEK and/or selumetinib sensitivity was shown in multiple independent melanoma, colon, breast, and lung tumor cell lines and in xenograft models. Furthermore, we were able to measure these signatures in fixed archival melanoma tumor samples using a single RT-qPCR-based test and found intergene correlations and associations with genetic markers of pathway activity to be preserved. These signatures offer useful tools for the study of MEK biology and clinical application of MEK inhibitors, and the novel approaches taken may benefit other targeted therapies.