Computational quantification of global effects induced by mutations and drugs in signaling networks of colorectal cancer cells.

Colorectal cancer (CRC) is one of the most deadly and commonly diagnosed tumors worldwide. Several genes are involved in its development and progression. The most frequent mutations concern APC, KRAS, SMAD4, and TP53 genes, suggesting that CRC relies on the concomitant alteration of the related pathways. However, with classic molecular approaches, it is not easy to simultaneously analyze the interconnections between these pathways. To overcome this limitation, recently these pathways have been included in a huge chemical reaction network (CRN) describing how information sensed from the environment by growth factors is processed by healthy colorectal cells. Starting from this CRN, we propose a computational model which simulates the effects induced by single or multiple concurrent mutations on the global signaling network. The model has been tested in three scenarios. First, we have quantified the changes induced on the concentration of the proteins of the network by a mutation in APC, KRAS, SMAD4, or TP53. Second, we have computed the changes in the concentration of p53 induced by up to two concurrent mutations affecting proteins upstreams in the network. Third, we have considered a mutated cell affected by a gain of function of KRAS, and we have simulated the action of Dabrafenib, showing that the proposed model can be used to determine the most effective amount of drug to be delivered to the cell. In general, the proposed approach displays several advantages, in that it allows to quantify the alteration in the concentration of the proteins resulting from a single or multiple given mutations. Moreover, simulations of the global signaling network of CRC may be used to identify new therapeutic targets, or to disclose unexpected interactions between the involved pathways.

Systems medicine in colorectal cancer: from a mathematical model toward a new type of clinical trial.

Current colorectal cancer (CRC) treatment guidelines are primarily based on clinical features, such as cancer stage and grade. However, outcomes may be improved using molecular treatment guidelines. Potentially useful biomarkers include driver mutations and somatically inherited alterations, signaling proteins (their expression levels and (post) translational modifications), mRNAs, micro-RNAs and long noncoding RNAs. Moving to an integrated system is potentially very relevant. To implement such an integrated system: we focus on an important region of the signaling network, immediately above the G1-S restriction point, and discuss the reconstruction of a Molecular Interaction Map and interrogating it with a dynamic mathematical model. Extensive model pretraining achieved satisfactory, validated, performance. The model helps to propose future target combination priorities, and restricts drastically the number of drugs to be finally tested at a cellular, in vivo, and clinical-trial level. Our model allows for the inclusion of the unique molecular profiles of each individual patient's tumor. While existing clinical guidelines are well established, dynamic modeling may be used for future targeted combination therapies, which may progressively become part of clinical practice within the near future. WIREs Syst Biol Med 2016, 8:314-336. doi: 10.1002/wsbm.1342 For further resources related to this article, please visit the WIREs website.

Systems Medicine in Oncology: Signaling Network Modeling and New-Generation Decision-Support Systems.

Two different perspectives are the main focus of this book chapter: (1) A perspective that looks to the future, with the goal of devising rational associations of targeted inhibitors against distinct altered signaling-network pathways. This goal implies a sufficiently in-depth molecular diagnosis of the personal cancer of a given patient. A sufficiently robust and extended dynamic modeling will suggest rational combinations of the abovementioned oncoprotein inhibitors. The work toward new selective drugs, in the field of medicinal chemistry, is very intensive. Rational associations of selective drug inhibitors will become progressively a more realistic goal within the next 3-5 years. Toward the possibility of an implementation in standard oncologic structures of technologically sufficiently advanced countries, new (legal) rules probably will have to be established through a consensus process, at the level of both diagnostic and therapeutic behaviors.(2) The cancer patient of today is not the patient of 5-10 years from now. How to support the choice of the most convenient (and already clinically allowed) treatment for an individual cancer patient, as of today? We will consider the present level of artificial intelligence (AI) sophistication and the continuous feeding, updating, and integration of cancer-related new data, in AI systems. We will also report briefly about one of the most important projects in this field: IBM Watson US Cancer Centers. Allowing for a temporal shift, in the long term the two perspectives should move in the same direction, with a necessary time lag between them.

Advances in dynamic modeling of colorectal cancer signaling-network regions, a path toward targeted therapies.

The interconnected network of pathways downstream of the TGFß, WNT and EGF-families of receptor ligands play an important role in colorectal cancer pathogenesis.We studied and implemented dynamic simulations of multiple downstream pathways and described the section of the signaling network considered as a Molecular Interaction Map (MIM). Our simulations used Ordinary Differential Equations (ODEs), which involved 447 reactants and their interactions.Starting from an initial "physiologic condition", the model can be adapted to simulate individual pathologic cancer conditions implementing alterations/mutations in relevant onco-proteins. We verified some salient model predictions using the mutated colorectal cancer lines HCT116 and HT29. We measured the amount of MYC and CCND1 mRNAs and AKT and ERK phosphorylated proteins, in response to individual or combination onco-protein inhibitor treatments. Experimental and simulation results were well correlated. Recent independently published results were also predicted by our model.Even in the presence of an approximate and incomplete signaling network information, a predictive dynamic modeling seems already possible. An important long term road seems to be open and can be pursued further, by incremental steps, toward even larger and better parameterized MIMs. Personalized treatment strategies with rational associations of signaling-proteins inhibitors, could become a realistic goal.

Parameter space exploration within dynamic simulations of signaling networks.

We started offering an introduction to very basic aspects of molecular biology, for the reader coming from computer sciences, information technology, mathematics. Similarly we offered a minimum of information about pathways and networks in graph theory, for a reader coming from the bio-medical sector. At the crossover about the two different types of expertise, we offered some definition about Systems Biology. The core of the article deals with a Molecular Interaction Map (MIM), a network of biochemical interactions involved in a small signaling-network sub-region relevant in breast cancer. We explored robustness/sensitivity to random perturbations. It turns out that our MIM is a non-isomorphic directed graph. For non physiological directions of propagation of the signal the network is quite resistant to perturbations. The opposite happens for biologically significant directions of signal propagation. In these cases we can have no signal attenuation, and even signal amplification. Signal propagation along a given pathway is highly unidirectional, with the exception of signal-feedbacks, that again have a specific biological role and significance. In conclusion, even a relatively small network like our present MIM reveals the preponderance of specific biological functions over unspecific isomorphic behaviors. This is perhaps the consequence of hundreds of millions of years of biological evolution.

Potent synergistic interaction between the Nampt inhibitor APO866 and the apoptosis activator TRAIL in human leukemia cells.

OBJECTIVE: The nicotinamide phosphoribosyltransferase (Nampt) inhibitor APO866 depletes intracellular nicotinamide adenine dinucleotide (NAD(+)) and shows promising anticancer activity in preclinical studies. Tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) binds to plasma membrane receptors DR4 and DR5 and induces apoptosis via caspase-8 and -10. Here we have explored the interaction between APO866 and TRAIL in leukemia cell lines and in primary B-cell chronic lymphocytic leukemia cells. MATERIALS AND METHODS: Cells were treated with APO866, TRAIL, or their combination. Viability and mitochondrial transmembrane potential (ΔΨ(m)) were determined by cell staining with propidium iodide and tetramethylrhodamine ethyl ester, respectively, and flow cytometry. Nampt and γ-tubulin levels, as well as caspase-3 cleavage were detected by immunoblotting. DR4 and DR5 expression were assessed by immunostaining and flow cytometry. Caspases were inhibited with zVAD-FMK and zDEVD-FMK; autophagy with 3-methyladenine, LY294002, and wortmannin. Intracellular NAD(+) and adenosine triphosphate (ATP) were measured by cycling assays and high-performance liquid chromatography (HPLC), respectively. RESULTS: APO866 induced NAD(+) depletion, ΔΨ(m) dissipation, and ATP shortage in leukemia cells, thereby leading to autophagic cell death. TRAIL induced caspase-dependent apoptosis. TRAIL addition to APO866 synergistically increased its activity in leukemia cells by enhancing NAD(+) depletion, ΔΨ(m) dissipation, and ATP shortage. No DR5 upregulation at the cell surface in response to APO866 was observed. Remarkably, in healthy leukocytes APO866 and TRAIL were poorly active and failed to show any cooperation. CONCLUSIONS: Activation of the extrinsic apoptotic cascade with TRAIL selectively amplifies the sequelae of Nampt inhibition in leukemia cells, and appears as a promising strategy to enhance APO866 activity in hematological malignancies.

Ras-induced resistance to lapatinib is overcome by MEK inhibition.

Lapatinib, a dual HER2 and EGFR tyrosine kinase inhibitor is highly active in HER2+ breast cancer. However, its efficacy is limited by either primary or acquired resistance. Although mutations in ras genes are rarely found in breast cancer, H-ras overexpression is frequently observed. Moreover, genetic alterations that do not directly involve ras such as Brk amplification, ultimately result in increased ras signaling. Using SKBR3 cells, a HER2+ breast cancer cell line that is naturally devoid of mutations in PI3KCA, PTEN, BRAF, and ras we show that both H-ras overexpression and expression of an oncogenic ras allele (ras V12) reduce susceptibility to lapatinib in analogy to what observed with activating PI3KCA mutations and with a constitutively active form of Akt. Importantly, we found that resistance to lapatinib due to ras overexpression or to ras V12 is overcome by MEK inhibition with U0126, suggesting a key role for the MEK-Erk pathway in ras-induced resistance. Similar results were obtained in BT474 cells, another HER+ breast cancer cell line. Therefore, our data indicate that overexpressed/mutated ras may act as a biological modifier of the response to lapatinib. Combining MEK inhibitors with lapatinib may help overcome this form of resistance and increase the efficacy of lapatinib in these tumors.

Apoptosis of B-cell chronic lymphocytic leukemia cells induced by a novel BH3 peptidomimetic.

B-cell chronic lymphocytic leukemia (B-CLL) is the most common leukemia in human adults of the Western world and no definitive cure is yet available. The disease is characterized by accumulation of clonal malignant B lymphocytes resistant to apoptosis. Strategies to hit the anti-apoptotic drift of the Bcl-2 family in B-CLL cells are being explored. A novel peptidomimetic based on the BH3 domain of the pro-apoptotic protein Bim and recently shown to exert significant apoptotic activity on acute myeloid leukemia cells, both in vitro and in vivo, was assayed on ex-vivo derived leukemic cells from untreated B-CLL patients (n = 7). We found that this peptide, named 072RB, induced apoptosis of B-CLL samples at a concentration that does not affect viability of peripheral and bone marrow derived lymphocytes from healthy donors. Apoptosis was demonstrated by activation of Bak and Bax, externalization of plasma membranes phosphadydilserines, appearance of hypodiploid events in DNA flow cytometry histograms and was accompanied by dissipation of the mitochondrial transmembrane potential. Before the onset of marked apoptotic signs a progressive decline of the relevant anti-apoptotic proteins Bcl-X(L) and Mcl-1 could be observed. The negative control peptide 072RBL94A was ineffective for B-CLL cells, supporting the sequence specificity of 072RB activity. No relationship was found between responsiveness to 072RB and Mcl-1/Bcl-X(L) basal levels or decrease magnitude, possibly because of the limited sample size of the study. Altogether, we demonstrate that 072RB induces significant apoptosis of B-CLL cells subsequent to Bcl-X(L) and Mcl-1 downregulation.

A novel Bim-BH3-derived Bcl-XL inhibitor: biochemical characterization, in vitro, in vivo and ex-vivo anti-leukemic activity.

BH3-only members of the Bcl-2 family exert a fundamental role in apoptosis induction. This work focuses on the development of a novel peptidic molecule based on the BH3 domain of Bim. The antiapoptotic molecule Bcl-X(L), involved in cancer development/progression and tumour resistance to cytotoxic drugs, is a target for Bim. According to a rational study of the structural interactions between wt Bim-BH3 and Bcl-X(L), we replaced specific residues of Bim-BH3 with natural and non-natural aminoacids and added an internalizing sequence, thus increasing dramatically the inhibitory activity of our modified Bim-BH3 peptide, called 072RB. Confocal microscopy and flow cytometry demonstrated cellular uptake and internalization of 072RB, followed by co-localization with mitochondria. Multiparameter flow cytometry demonstrated that the 072RB dose-dependent growth inhibition of leukaemia cell lines was due to apoptotic cell death. No effect was observed when cells were treated with the internalizing vector alone or a mutated control peptide (single aminoacid substitution L94A). Ex-vivo derived leukemic cells from acute myeloid leukaemia (AML) patients underwent cell death when cultured in vitro in the presence of 072RB. Conversely, no significant cytotoxic effect was observed when 072RB was administered to cultures of peripheral blood mononuclear cells, either resting or PHA-stimulated, and bone marrow cells of normal donors. Xenografts of human AML cells in NOD/SCID mice displayed a significant delay of leukemic cell growth upon treatment with 072RB administered intravenously (15 mg/Kg three times, 48 hours after tumour cell injection). Altogether, these observations support the therapeutic potentials of this novel BH3 mimetic.

Differential expression of nuclear lamins in normal and cancerous prostate tissues.

The process of carcinogenesis is characterized by definite changes in the protein composition of the nuclear matrix. We have recently found that lamins form, in addition to the nuclear lamina, an intranuclear web of thin fibrils. This finding prompted us to address the question of whether changes in the expression of lamins occur in the course of tumor development. In prostate cancer, lamin B undergoes a significant increase; interestingly, its nuclear content strongly correlates with tumor differentiation. Moreover, all the lamins show reproducible alterations in the distribution of the isoelectric variants, suggesting that dephosphorylation events could trigger changes in the pattern of gene expression by inducing structural rearrangements of the nuclear scaffold.

SCN- binding to the charged lysines of histones end domains mimics acetylation and shows the major histone-DNA interactions involved in eu and heterochromatin stabilization.

SCN- binds to the charged amino group of lysines, inducing local changes in the electrostatic free energy of histones. We exploited this property to selectively perturb the histone-DNA interactions involved in the stabilization of eu and heterochromatin. Differential scanning calorimetry (DSC) was used as leading technique in combination with trypsin digestion that selectively cleaves the histone end domains. Euchromatin undergoes progressive destabilization with increasing KSCN concentration from 0 to 0.3 M. Trypsin digestion in the presence of 0.2 M KSCN show that the stability of the linker decreases as a consequence of the competitive binding of SCN- to the amino groups located in the C and N-terminal domain of H1 and H3, respectively; likewise, the release of the N-terminal domain of H4 induces an appreciable depression in both the temperature and enthalpy of melting of core particle DNA. Unfolding of heterochromatin requires, in addition to further cleavage of H4, extensive digestion of H2A and H2B, strongly suggesting that these histones stabilize the higher order structure by forming a protein network which extends throughout the heterochromatin domain.

A new assay for the discovery of Bcl-XL inhibitors.

The Bcl-2 family of antiapoptotic proteins is commonly over expressed in many types of human cancer and remains one of the few validated targets. Antiapoptotic family proteins such as Bcl-2 and Bcl-XL function, at least in part, by binding proapoptotic members such as Bax and Bak and thereby prevent release of the apoptotic cascade of events. "BH3-only" members of the family disrupt this interaction by binding, via their BH3 domain, to a hydrophobic pocket on the surface of the antiapoptotic members. Disruption of heterodimerization could be used to modulate cell death reinstating apoptosis in cancer cells. An affinity displacement assay based on Bcl-XL/BH3 interaction has been developed. This assay makes use of soluble His-tagged Bcl-XL and fluorescein tagged BH3. Binding is measured as fluorescence associated with magnetic beads. The assay was miniaturized to 96-well microtiter plates and can be employed in high throughput screening (HTS), in addition it is robust enough to be applied to microbial fermentation extracts.

Meta-analysis of the role of p53 status in isogenic systems tested for sensitivity to cytotoxic antineoplastic drugs.

The role of p53 in modifying sensitivity to cytotoxic drugs has been commonly studied by creating transfection pairs of wt p53 parental cells and altered p53 daughter cells, or vice versa. Authors inevitably tended to extrapolate and generalize their experimental observations, and conflicting reports have been more the rule than the exception. We have performed a meta-analysis of 356 independent studies. Average changes of drug sensitivity after a change of p53 status were observed. E6 transfection predominantly induces sensitization to cytotoxic drugs, whereas p53-/- knockout cells are more drug-resistant than their normal p53+/+ counterpart. Unexpectedly, transfection with a mutated p53 does not change much the drug sensitivity of most wt p53 cancer lines, with the notable exception of A2780, a predominant cell line in the studies analyzed (A2780 cells show increased resistance after transfection with a mutated p53). Rather interestingly, mitotic spindle poisons acted differently from other classes of cytotoxic drugs. A crucial indication of our findings is that the role of p53 alone in determining sensitivity/resistance to cytotoxic drugs is limited: the individual molecular pathology and differentiation of a given cancer line prevail over any average trend, and are causal to a broad spreading of the data. We also identify major "confounding factors", alias independent categorical variables, capable of affecting the average outcome.

Molecular interaction maps--a diagrammatic graphical language for bioregulatory networks.

Molecular interaction maps (MIMs) use a clear, accurate, and versatile graphical language to depict complex biological processes. Here, we discuss the main features of the MIM language and its potential uses. MIMs can be used as database resources and simulation guides, and can serve to generate new hypotheses regarding the roles of specific molecules in the bioregulatory networks that control progression through the cell cycle, differentiation, and cell death.

Structural investigation of proapoptotic peptide by CD and NMR spectroscopy.

We have performed a systematic investigation of the structural features of the peptides Int (a sequence able to cross cell membranes) and Int-H1(S6A,F8A) (which shows interesting antitumoral properties). After screening in aqueous solution at different ionic strength and pH values, we analyzed the structures of the peptides in different water/trifluoroethanol mixtures by Circular Dichroism and Nuclear Magnetic Resonance techniques.

DNA supercoiling in apoptotic chromatin.

In a previous paper, we have reported that in rat thymocyte apoptosis chromatin undergoes a specific structural change as well as an appreciable increase in the unacetylated forms of histones H3 and H4. Here, we show that H3 and H4 deacetylation bears no relation to chromatin condensation, and present new ultrastructural and topological observations that largely clarify the organization of the condensed state. The texture of the latter corresponds to a closely woven network of negatively supercoiled 11 nm fibers, as shown by both ultrastructural observations and relaxation experiments using ethidium bromide. Circularly closed chromatin loops undergoing apoptotic condensation, clearly showing nucleosome compact dimers or higher oligomers, as well as long stretches of supercoiled DNA, have also been detected. All of these modifications are strongly reminiscent of the alterations induced in nucleosome bearing plasmids by the chromatin remodeling factors SWI/SNF and RSC.

Nuclear matrix proteins changes in cancerous prostate tissues and their prognostic value in clinically localized prostate cancer.

BACKGROUND: After the discovery that nuclear matrix (NM) directs the spatial organization of DNA transcription and replication, there has been an increasing interest in studying NM changes associated with malignant transformation and their potential usefulness in the clinical setting. METHODS: High-resolution two-dimensional gel electrophoresis was used to analyze the NM proteins (NMP) of specimens of prostate cancer tissue obtained from the prostates of 75 patients undergoing retropubic prostatectomy. RESULTS: Nine NMP with different molecular weights and isoelectric points have been identified. They were expressed differently by prostate cancer tissues. An increasing trend toward the expression of such proteins like NMP 6-8 was evident with increasing tumor stage and dedifferentiation. NMP 6-8 were also significantly correlated with the risk of biochemical progression. However, Gleason score was the only significant discriminant in this regard in multiparametric models. CONCLUSIONS: This study confirms that prostate cancer progression is related to profound changes in NMP expression patterns. However, due to the complexity of the methods required to define these latter, the clinical relevance of NMP appears to be still limited.

An intranuclear frame for chromatin compartmentalization and higher-order folding.

Recent ultrastructural, immunoelectron, and confocal microscopy observations done in our laboratory [Barboro et al. [2002] Exp. Cell. Res. 279:202-218] have confirmed that lamins and the nuclear mitotic apparatus protein (NuMA) are localized inside the interphase nucleus in a polymerized form. This provided evidence of the existence of a RNA stabilized lamin/NuMA frame, consisting of a web of thin ( approximately 3 and approximately 5 nm) lamin filaments to which NuMA is anchored mainly in the form of discrete islands, which might correspond to the minilattices described by Harborth et al. [1999] (EMBO. J. 18:1689-1700). In this article we propose that this scaffold is involved in the compartmentalization of both chromatin and functional domains and further determines the higher-order nuclear organization. This hypothesis is strongly supported by the scrutiny of different structural transitions which occur inside the nucleus, such as chromatin displacement and rearrangements, the collapse of the internal nuclear matrix after RNA digestion and the disruption of chromosome territories induced by RNase A and high salt treatment. All of these destructive events directly depend on the loss of the stabilizing effect exerted on the different levels of structural organization by the interaction of RNA with lamins and/or NuMA. Therefore, the integrity of nuclear RNA must be safeguarded as far as possible to isolate the matrix in the native form. This material will allow for the first time the unambiguous ultrastructural localization inside the INM of the components of the functional domains, so opening new avenues of investigation on the mechanisms of gene expression in eukaryotes.