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.

Use of reverse phase protein microarrays and reference standard development for molecular network analysis of metastatic ovarian carcinoma

Cancer can be defined as a deregulation or hyperactivity in the ongoing network of intracellular and extracellular signaling events. Reverse phase protein microarray technology may offer a new opportunity to measure and profile these signaling pathways, providing data on post-translational phosphorylation events not obtainable by gene microarray analysis. Treatment of ovarian epithelial carcinoma almost always takes place in a metastatic setting since unfortunately the disease is often not detected until later stages. Thus, in addition to elucidation of the molecular network within a tumor specimen, critical questions are to what extent do signaling changes occur upon metastasis and are there common pathway elements that arise in the metastatic microenvironment. For individualized combinatorial therapy, ideal therapeutic selection based on proteomic mapping of phosphorylation end points may require evaluation of the patient's metastatic tissue. Extending these findings to the bedside will require the development of optimized protocols and reference standards. We have developed a reference standard based on a mixture of phosphorylated peptides to begin to address this challenge.

JK1 (FAM134B) represses cell migration in colon cancer : a functional study of a novel gene

Background JK1 is a novel cancer-related gene with unknown functional role in carcinogenesis. The aim of this study is to investigate the role of JK1 gene in carcinogenesis in an in vitro cell proliferation and migration analysis model. Methods Small hairpin RNAs (shRNA) were designed to knock-down JK1 expression in colon cancer cell line (SW480) using transduction ready lentiviral particles. Cell proliferation and cell migration assays were performed on multiple extracellular matrices to investigate the cellular effects of JK1 in colon cancer cells. A non-cancer colonic epithelial cell line (FHC) was used to compare the expression of JK1 in cancer cell line. Results JK1 knock-down did not affect cellular proliferation or survival in colon cancer. However, the manipulation increased cancer cell migration rates on collagen and fibronectin substrates. Conclusions JK1 was shown for the first time to have a functional role in the pathogenesis of colon cancer. The results imply that JK1 represses the capacity of cancer cells to migrate within their tissue. They also concurred with the previous findings of JK1 activity correlations with clinical and pathological features in colon cancer. The capacity may have utility as a means to prevent cancer cells forming metastases.

Correlation between BRAF mutation and the clinicopathological parameters in papillary thyroid carcinoma with particular reference to follicular variant

Mutation of the BRAF gene is common in thyroid cancer. Follicular variant of papillary thyroid carcinoma is a variant of papillary thyroid carcinoma that has created continuous diagnostic controversies among pathologists. The aims of this study are to (1) investigate whether follicular variant of papillary thyroid carcinoma has a different pattern of BRAF mutation than conventional papillary thyroid carcinoma in a large cohort of patients with typical features of follicular variant of papillary thyroid carcinoma and (2) to study the relationship of clinicopathological features of papillary thyroid carcinomas with BRAF mutation. Tissue blocks from 76 patients with diagnostic features of papillary thyroid carcinomas (40 with conventional type and 36 with follicular variant) were included in the study. From these, DNA was extracted and BRAF V600E mutations were detected by polymerase chain reaction followed by restriction enzyme digestion and sequencing of exon 15. Analysis of the data indicated that BRAF V600E mutation is significantly more common in conventional papillary thyroid carcinoma (58% versus 31%, P = .022). Furthermore, the mutation was often noted in female patients (P = .017), in high-stage cancers (P = .034), and in tumors with mild lymphocytic thyroiditis (P = .006). We concluded that follicular variant of papillary thyroid carcinoma differs from conventional papillary thyroid carcinoma in the rate of BRAF mutation. The results of this study add further information indicating that mutations in BRAF play a role in thyroid cancer development and progression.

MAS NMR measurements of intact articular bovine cartilage

Articular cartilage (AC), an avascular connective tissue lining articulating surfaces of the long bones, comprises extracellular biopolymers. In functionally compromised states such as osteoarthritis, thinned or lost AC causes reduced mobility and increased health-care costs. Understanding of the characteristics responsible for the load bearing efficiency of AC and the factors leading to its degradation are incomplete. DTI shows the structural alignment of collagen in AC [1] and T2 relaxation measurements suggest that the average director of reorientational motion of water molecules depends on the degree of alignment of collagen in AC [2]. Information on the nature of the chemical interactions involved in functional AC is lacking. The need for AC structural integrity makes solid state NMR an ideal tool to study this tissue. We examined the contribution of water in different functional ‘compartments’ using 1H-MAS, 13C-MAS and 13C-CPMAS NMR of bovine patellar cartilage incubated in D2O. 1H-MAS spectra signal intensity was reduced due to H/D exchange without a measureable redistribution of relative signal intensity. Chemical shift anisotropy was estimated by lineshape analysis of multiple peaks in the 1H-MAS spinning sidebands. These asymmetrical sidebands suggested the presence of multiple water species in AC. Therefore, water was added in small aliquots to D2O saturated AC and the influence of H2O and D2O on organic components was studied with 13C-MAS-NMR and 13C-CPMAS-NMR. Signal intensity in 13C-MAS spectra showed no change in relative signal intensity throughout the spectrum. In 13C-CPMAS spectra, displacement of water by D2O resulted in a loss of signal in the aliphatic region due to a reduction in proton availability for cross-polarization. These results complement dehydration studies of cartilage using osmotic manipulation [3] and demonstrate components of cartilage that are in contact with mobile water.

Solid-state MAS NMR measurements of intact articular bovine cartilage

Articular cartilage (AC), an avascular connective tissue lining articulating surfaces of the long bones, comprises extracellular biopolymers. In functionally compromised states such as osteoarthritis, thinned or lost AC causes reduced mobility and increased health-care costs. Understanding of the characteristics responsible for the load bearing efficiency of AC and the factors leading to its degradation are incomplete. DTI shows the structural alignment of collagen in AC [1] and T2 relaxation measurements suggest that the average director of reorientational motion of water molecules depends on the degree of alignment of collagen in AC [2]. Information on the nature of the chemical interactions involved in functional AC is lacking. The need for AC structural integrity makes solid state NMR an ideal tool to study this tissue. We examined the contribution of water in different functional ‘compartments’ using 1H-MAS, 13C-MAS and 13C-CPMAS NMR of bovine patellar cartilage incubated in D2O. 1H-MAS spectra signal intensity was reduced due to H/D exchange without a measureable redistribution of relative signal intensity. Chemical shift anisotropy was estimated by lineshape analysis of multiple peaks in the 1H-MAS spinning sidebands. These asymmetrical sidebands suggested the presence of multiple water species in AC. Therefore, water was added in small aliquots to D2O saturated AC and the influence of H2O and D2O on organic components was studied with 13C-MAS-NMR and 13C-CPMAS-NMR. Signal intensity in 13C-MAS spectra showed no change in relative signal intensity throughout the spectrum. In 13C-CPMAS spectra, displacement of water by D2O resulted in a loss of signal in the aliphatic region due to a reduction in proton availability for cross-polarization. These results complement dehydration studies of cartilage using osmotic manipulation [3] and demonstrate components of cartilage that are in contact with mobile water.

Role of the N-terminal region in G protein-coupled receptor functions : negative modulation revealed by 5-HT2B receptor polymorphisms

The putative role of the N-terminal region of rhodopsin-like 7 transmembrane biogenic amine receptors in agonist-induced signaling has not yet been clarified despite recent advances in 7 transmembrane receptor structural biology. Given the existence of N-terminal nonsynonymous polymorphisms (R6G;E42G) within the HTR2B gene in a drug-abusing population, we assessed whether these polymorphisms affect 5-hydroxytryptamine 2B (5-HT2B) receptor in vitro pharmacologic and coupling properties in transfected COS-7 cells. Modification of the 5-HT2B receptor N terminus by the R6G;E42G polymorphisms increases such agonist signaling pathways as inositol phosphate accumulation as assessed by either classic or operational models. The N-terminal R6G;E42G mutations of the 5-HT2B receptor also increase cell proliferation and slow its desensitization kinetics compared with the wild-type receptor, further supporting a role for the N terminus in transduction efficacy. Furthermore, by coexpressing a tethered wild-type 5-HT2B receptor N terminus with a 5-HT2B receptor bearing a N-terminal deletion, we were able to restore original coupling. This reversion to normal activity of a truncated 5-HT2B receptor by coexpression of the membrane-tethered wild-type 5-HT2B receptor N terminus was not observed using a membrane-tethered 5-HT2B receptor R6G;E42G N terminus. These data suggest that the N terminus exerts a negative control over basal as well as agonist-stimulated receptor activity that is lost in the R6G;E42G mutant. Our findings reveal a new and unanticipated role of the 5-HT2B receptor N terminus as a negative modulator, affecting both constitutive and agonist-stimulated activity. Moreover, our data caution against excluding the N terminus and extracellular loops in structural studies of this 7 transmembrane receptor family

5-HT(2B) receptors are required for serotonin-selective antidepressant actions

The therapeutic effects induced by serotonin-selective reuptake inhibitor (SSRI) antidepressants are initially triggered by blocking the serotonin transporter and rely on long-term adaptations of pre- and post-synaptic receptors. We show here that long-term behavioral and neurogenic SSRI effects are abolished after either genetic or pharmacological inactivation of 5-HT(2B) receptors. Conversely, direct agonist stimulation of 5-HT(2B) receptors induces an SSRI-like response in behavioral and neurogenic assays. Moreover, the observation that (i) this receptor is expressed by raphe serotonergic neurons, (ii) the SSRI-induced increase in hippocampal extracellular serotonin concentration is strongly reduced in the absence of functional 5-HT(2B) receptors and (iii) a selective 5-HT(2B) agonist mimics SSRI responses, supports a positive regulation of serotonergic neurons by 5-HT(2B) receptors. The 5-HT(2B) receptor appears, therefore, to positively modulate serotonergic activity and to be required for the therapeutic actions of SSRIs. Consequently, the 5-HT(2B) receptor should be considered as a new tractable target in the combat against depression.

Comparative analysis of the Metarhizium anisopliae secretome in response to exposure to the greyback cane grub and grub cuticles

Metarhizium anisopliae is a well-characterized biocontrol agent of a wide range of insects including cane grubs. In this study, a two-dimensional (2D) electrophoresis was used to display secreted proteins of M. anisopliae strain FI-1045 growing on the whole greyback cane grubs and their isolated cuticles. Hydrolytic enzymes secreted by M. anisopliae play a key role in insect cuticle-degradation and initiation of the infection process. We have identified all the 101 protein spots displayed by cross-species identification (CSI) from the fungal kingdom. Among the identified proteins were 64-kDa serine carboxypeptidase, 1,3 beta-exoglucanase, Dynamin GTPase, THZ kinase, calcineurin like phosphoesterase, and phosphatidylinositol kinase secreted by M. ansiopliae (FI-1045) in response to exposure to the greyback cane grubs and their isolated cuticles. These proteins have not been previously identified from the culture supernatant of M. anisopliae during infection. To our knowledge, this the first proteomic map established to study the extracellular proteins secreted by M. ansiopliae (FI-1045) during infection of greyback cane grubs and its cuticles.

Encapsulation of Hydrocortisone and Mesalazine in Zein Microparticles

Zein was investigated for use as an oral-drug delivery system by loading prednisolone into zein microparticles using coacervation. To investigate the adaptability of this method to other drugs, zein microparticles were loaded with hydrocortisone, which is structurally related to prednisolone; or mesalazine, which is structurally different having a smaller LogP and ionizable functional groups. Investigations into the in vitro digestibility, and the electrophoretic profile of zein, and zein microparticles were conducted to shed further insight on using this protein as a drug delivery system. Hydrocortisone loading into zein microparticles was comparable with that reported for prednisolone, but mesalazine loading was highly variable. Depending on the starting quantities of hydrocortisone and zein, the average amount of microparticles equivalent to 4 mg hydrocortisone, (a clinically used dose), ranged from 60-115 mg, which is realistic and practical for oral dosing. Comparatively, an average of 2.5 g of microparticles was required to deliver 250 mg of mesalazine (a clinically used dose), so alternate encapsulation methods that can produce higher and more precise mesalazine loading are required. In vitro protein digestibility revealed that zein microparticles were more resistant to digestion compared to the zein raw material, and that individual zein peptides are not preferentially coacervated into the microparticles. In combination, these results suggest that there is potential to formulate a delivery system based on zein microparticles made using specific subunits of zein that is more resistant to digestion as starting material, to deliver drugs to the lower gastrointestinal tract.

The infectivity and host range of Orgyia anartoides nucleopolyhedrovirus

The painted apple moth (PAM), Teia anartoides (Walker) (Lepidoptera: Lymantriidae) made a recent incursion into New Zealand. A nucleopolyhedrovirus (NPV), Orgyia anartoides NPV (OranNPV), originally isolated from PAM in Australia, was tested for its pathogenicity to PAM and a range of non-target insect species found in New Zealand, to evaluate its suitability as a microbial control for this insect invader. Dosage-mortality tests showed that OranNPV was highly pathogenic to PAM larvae; mean LT50 values for third instars ranged from 17.9 to 8.1 days for doses from 102 to 105 polyhedral inclusion bodies/larva, respectively. The cause of death in infected insects was confirmed as OranNPV. Molecular analysis established that OranNPV can be identified by PCR and restriction digestion, and this process complemented microscopic examination of infected larvae. No lymantriid species occur in New Zealand; however, the virus had no significant effects on species from five other lepidopteran families (Noctuidae, Tortricidae, Geometridae, Nymphalidae and Plutellidae) or on adult honeybees. Thus, all indications from this initial investigation are that OranNPV would be an important tool in the control of PAM in a future incursion of this species into New Zealand.

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.

Activation of forkhead box O transcription factors by oncogenic BRAF promotes p21cip1-dependent senescence

Oncogene-induced senescence (OIS) is a potent tumor-suppressive mechanism that is thought to come at the cost of aging. The Forkhead box O (FOXO) transcription factors are regulators of life span and tumor suppression. However, whether and how FOXOs function in OIS have been unclear. Here, we show a role for FOXO4 in mediating senescence by the human BRAFV600E oncogene, which arises commonly in melanoma. BRAFV600E signaling through mitogen-activated protein kinase/extracellular signal-regulated kinase kinase resulted in increased reactive oxygen species levels and c-Jun NH 2 terminal kinase-mediated activation of FOXO4 via its phosphorylation on Thr223, Ser226, Thr447, and Thr451. BRAFV600E-induced FOXO4 phosphorylation resulted in p21cip1-mediated cell senescence independent of p16 ink4a or p27kip1. Importantly, melanocyte-specific activation of BRAFV600E in vivo resulted in the formation of skin nevi expressing Thr223/Ser226-phosphorylated FOXO4 and elevated p21cip1. Together, these findings support a model in which FOXOs mediate a trade-off between cancer and aging.

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.