Modeling the effect of APC truncation on destruction complex function in colorectal cancer cells

Cancer‑associated fibroblast‑derived LRRC15 promotes the migration and invasion of triple‑negative breast cancer cells via Wnt/β‑catenin signalling pathway regulation

Triple-negative breast cancer (TNBC) is a highly aggressive tumour subtype associated with poor prognosis. The function of leucine-rich repeat-containing protein 15 (LRRC15), a member of the leucine-rich repeat superfamily, in TNBC has not yet been elucidated. The aim of this study was to identify the combined role of LRRC15 and Wnt/β-catenin signalling pathway in the development of TNBC. The expression of LRRC15 in TNBC tissues was analysed using data from The Cancer Genome Atlas. Cell migration and invasion assays were conducted to study the function of LRRC15 in TNBC. The expression of Wnt/β-catenin signalling proteins was analysed via western blotting. The effect of LRRC15 on β-catenin nuclear localisation was measured by performing western blotting and luciferase assays. It was found that high LRRC15 expression was associated with poor prognosis in patients with TNBC. High expression of LRRC15 in cancer-associated fibroblasts (CAFs) promoted cell migration and invasion in TNBC cells. In addition, TNBC cells with LRRC15 overexpression in CAFs showed an aberrant increase in β-catenin activity concomitant with nuclear localisation of β-catenin, which inhibited its degradation. These results showed that LRRC15 promoted tumour migration and invasion in TNBC cells by regulating the Wnt/β-catenin signalling pathway.

Reconstitution of the destruction complex defines roles of AXIN polymers and APC in β-catenin capture, phosphorylation, and ubiquitylation

The Wnt/β-catenin pathway is a highly conserved, frequently mutated developmental and cancer pathway. Its output is defined mainly by β-catenin's phosphorylation- and ubiquitylation-dependent proteasomal degradation, initiated by the multi-protein β-catenin destruction complex. The precise mechanisms underlying destruction complex function have remained unknown, largely because of the lack of suitable in vitro systems. Here we describe the in vitro reconstitution of an active human β-catenin destruction complex from purified components, recapitulating complex assembly, β-catenin modification, and degradation. We reveal that AXIN1 polymerization and APC promote β-catenin capture, phosphorylation, and ubiquitylation. APC facilitates β-catenin's flux through the complex by limiting ubiquitylation processivity and directly interacts with the SCFβ-TrCP E3 ligase complex in a β-TrCP-dependent manner. Oncogenic APC truncation variants, although part of the complex, are functionally impaired. Nonetheless, even the most severely truncated APC variant promotes β-catenin recruitment. These findings exemplify the power of biochemical reconstitution to interrogate the molecular mechanisms of Wnt/β-catenin signaling.

Analysis of Fn14-NF-κB signaling response dynamics using a mechanistic model

Fn14 is a transmembrane receptor protein belonging to the tumor necrosis factor receptor (TNFR) superfamily. Many experimental reports have shown that crosslinking of the receptor by its extracellular ligand TWEAK induces prolonged activation of transcription factor NF-κB. This behavior is distinct from TNF-α receptor, which is a more well-characterized member of the TNFR family. TNF-α receptor, despite sharing many similar molecular interactions with Fn14, only transiently activates NF-κB in response to TNF-α stimulation. Here, we investigate molecular mechanisms that enable Fn14 to display such distinctive behavior. In particular, we focus on two specific features of the Fn14 pathway that potentially give rise to a positive feedback regulation and differentiate it from the TNF-α receptor signaling. By developing a mechanistic model, we analyze how these features may determine the dynamics of an Fn14-NF-κB response. Our analysis reveals that stimulation of Fn14 by TWEAK may generate highly non-linear dynamics, including stable limit cycles and bistable responses. The type of response depends both on the strength and duration of a TWEAK signal. Our predictions and analyses also show that the molecular interactions underlying the positive feedback explain the prolonged activation of NF-κB under certain parameter regimes. In light of the model predictions, we propose possible deregulations of Fn14 leading to its overexpression in solid tumors and tissue injuries.

Multiscale Agent-Based and Hybrid Modeling of the Tumor Immune Microenvironment

Multiscale systems biology and systems pharmacology are powerful methodologies that are playing increasingly important roles in understanding the fundamental mechanisms of biological phenomena and in clinical applications. In this review, we summarize the state of the art in the applications of agent-based models (ABM) and hybrid modeling to the tumor immune microenvironment and cancer immune response, including immunotherapy. Heterogeneity is a hallmark of cancer; tumor heterogeneity at the molecular, cellular, and tissue scales is a major determinant of metastasis, drug resistance, and low response rate to molecular targeted therapies and immunotherapies. Agent-based modeling is an effective methodology to obtain and understand quantitative characteristics of these processes and to propose clinical solutions aimed at overcoming the current obstacles in cancer treatment. We review models focusing on intra-tumor heterogeneity, particularly on interactions between cancer cells and stromal cells, including immune cells, the role of tumor-associated vasculature in the immune response, immune-related tumor mechanobiology, and cancer immunotherapy. We discuss the role of digital pathology in parameterizing and validating spatial computational models and potential applications to therapeutics.

Multicellular Models Bridging Intracellular Signaling and Gene Transcription to Population Dynamics

Cell signaling and gene transcription occur at faster time scales compared to cellular death, division, and evolution. Bridging these multiscale events in a model is computationally challenging. We introduce a framework for the systematic development of multiscale cell population models. Using message passing interface (MPI) parallelism, the framework creates a population model from a single-cell biochemical network model. It launches parallel simulations on a single-cell model and treats each stand-alone parallel process as a cell object. MPI mediates cell-to-cell and cell-to-environment communications in a server-client fashion. In the framework, model-specific higher level rules link the intracellular molecular events to cellular functions, such as death, division, or phenotype change. Cell death is implemented by terminating a parallel process, while cell division is carried out by creating a new process (daughter cell) from an existing one (mother cell). We first demonstrate these capabilities by creating two simple example models. In one model, we consider a relatively simple scenario where cells can evolve independently. In the other model, we consider interdependency among the cells, where cellular communication determines their collective behavior and evolution under a temporally evolving growth condition. We then demonstrate the framework’s capability by simulating a full-scale model of bacterial quorum sensing, where the dynamics of a population of bacterial cells is dictated by the intercellular communications in a time-evolving growth environment.

Integrative analysis of oncogenic fusion genes and their functional impact in colorectal cancer

Background

Fusion genes are good candidates of molecular targets for cancer therapy. However, there is insufficient research on the clinical implications and functional characteristics of fusion genes in colorectal cancer (CRC).

Methods

In this study, we analysed RNA sequencing data of CRC patients (147 tumour and 47 matched normal tissues) to identify oncogenic fusion genes and evaluated their role in CRC.

Results

We validated 24 fusion genes, including novel fusions, by three algorithms and Sanger sequencing. Fusions from most patients were mutually exclusive CRC oncogenes and included tumour suppressor gene mutations. Eleven fusion genes from 13 patients (8.8%) were determined as oncogenic fusion genes by analysing their gene expression and function. To investigate their oncogenic impact, we performed proliferation and migration assays of CRC cell lines expressing fusion genes of GTF3A-CDK8, NAGLU- IKZF3, RNF121- FOLR2, and STRN-ALK. Overexpression of these fusion genes increased cell proliferation except GTF3A-CDK8. In addition, overexpression of NAGLU-IKZF3 enhanced migration of CRC cells. We demonstrated that NAGLU-IKZF3, RNF121-FOLR2, and STRN-ALK had tumourigenic effects in CRC.

Conclusion

In summary, we identified and characterised oncogenic fusion genes and their function in CRC, and implicated NAGLU-IKZF3 and RNF121-FOLR2 as novel molecular targets for personalised medicine development.

Timescale Separation of Positive and Negative Signaling Creates History-Dependent Responses to IgE Receptor Stimulation

The high-affinity receptor for IgE expressed on the surface of mast cells and basophils interacts with antigens, via bound IgE antibody, and triggers secretion of inflammatory mediators that contribute to allergic reactions. To understand how past inputs (memory) influence future inflammatory responses in mast cells, a microfluidic device was used to precisely control exposure of cells to alternating stimulatory and non-stimulatory inputs. We determined that the response to subsequent stimulation depends on the interval of signaling quiescence. For shorter intervals of signaling quiescence, the second response is blunted relative to the first response, whereas longer intervals of quiescence induce an enhanced second response. Through an iterative process of computational modeling and experimental tests, we found that these memory-like phenomena arise from a confluence of rapid, short-lived positive signals driven by the protein tyrosine kinase Syk; slow, long-lived negative signals driven by the lipid phosphatase Ship1; and slower degradation of Ship1 co-factors. This work advances our understanding of mast cell signaling and represents a generalizable approach for investigating the dynamics of signaling systems.

Extracellular inhibitors can attenuate tumorigenic Wnt pathway activity in adenomatous polyposis coli mutants: Predictions of a validated mathematical model

Background

Despite considerable investigational efforts, no method to overcome the pathogenesis caused by loss of function (LoF) mutations in tumor suppressor genes has been successfully translated to the clinic. The most frequent LoF mutation in human cancers is Adenomatous polyposis coli (APC), causing aberrant activation of the Wnt pathway. In nearly all colon cancer tumors, the APC protein is truncated, but still retains partial binding abilities.

Objective & methods

Here, we tested the hypothesis that extracellular inhibitors of the Wnt pathway, although acting upstream of the APC mutation, can restore normal levels of pathway activity in colon cancer cells. To this end, we developed and simulated a mathematical model for the Wnt pathway in different APC mutants, with or without the effects of the extracellular inhibitors, Secreted Frizzled-Related Protein1 (sFRP1) and Dickhopf1 (Dkk1). We compared our model predictions to experimental data in the literature.

Results

Our model accurately predicts T-cell factor (TCF) activity in mutant cells that vary in APC mutation. Model simulations suggest that both sFRP1 and DKK1 can reduce TCF activity in APC^1638N/1572T and Apc^min/min mutants, but restoration of normal activity levels is possible only in the former. When applied in combination, synergism between the two inhibitors can reduce their effective doses to one-fourth of the doses required under single inhibitor application. Overall, re-establishment of normal Wnt pathway activity is predicted for every APC mutant in whom TCF activity is increased by up to 11 fold.

Conclusions

Our work suggests that extracellular inhibitors can effectively restore normal Wnt pathway activity in APC-truncated cancer cells, even though these LoF mutations occur downstream of the inhibitory action. The insufficient activity of the truncated APC can be quantitatively balanced by the upstream intervention. This new concept of upstream intervention to control the effects of downstream mutations may be considered also for other partial LoF mutations in other signaling pathways.

Predictive Modeling of Drug Treatment in the Area of Personalized Medicine

Despite a growing body of knowledge on the mechanisms underlying the onset and progression of cancer, treatment success rates in oncology are at best modest. Current approaches use statistical methods that fail to embrace the inherent and expansive complexity of the tumor/patient/drug interaction. Computational modeling, in particular mechanistic modeling, has the power to resolve this complexity. Using fundamental knowledge on the interactions occurring between the components of a complex biological system, large-scale in silico models with predictive capabilities can be generated. Here, we describe how mechanistic virtual patient models, based on systematic molecular characterization of patients and their diseases, have the potential to shift the theranostic paradigm for oncology, both in the fields of personalized medicine and targeted drug development. In particular, we highlight the mechanistic modeling platform ModCell™ for individualized prediction of patient responses to treatment, emphasizing modeling techniques and avenues of application.

Network and systems biology: essential steps in virtualising drug discovery and development

The biological processes that keep us healthy or cause disease, as well as the mechanisms of action of possible drugs are inherently complex. In the face of this complexity, attempts at discovering new drugs to treat diseases have alternated between trial-and-error (typically on experimental systems) and grand simplification, usually based on much too little information. We now have the chance to combine these strategies through establishment of 'virtual patient' models, centred on a detailed molecular characterisation of thousands or even, in the future, millions of patients. In doing so, we lay the foundations for truly personalised therapy, as well as a far-reaching virtualisation of drug discovery and development in oncology and other areas of medicine.

Modeling for (physical) biologists: an introduction to the rule-based approach

Models that capture the chemical kinetics of cellular regulatory networks can be specified in terms of rules for biomolecular interactions. A rule defines a generalized reaction, meaning a reaction that permits multiple reactants, each capable of participating in a characteristic transformation and each possessing certain, specified properties, which may be local, such as the state of a particular site or domain of a protein. In other words, a rule defines a transformation and the properties that reactants must possess to participate in the transformation. A rule also provides a rate law. A rule-based approach to modeling enables consideration of mechanistic details at the level of functional sites of biomolecules and provides a facile and visual means for constructing computational models, which can be analyzed to study how system-level behaviors emerge from component interactions.

Mathematical modelling suggests a differential impact of β-transducin repeat-containing protein paralogues on Wnt/β-catenin signalling dynamics

The Wnt/β-catenin signalling pathway is involved in the regulation of a multitude of cellular processes by controlling the concentration of the transcriptional regulator β-catenin. Proteasomal degradation of β-catenin is mediated by two β-transducin repeat-containing protein paralogues, homologous to Slimb protein (HOS) and F-box/WD repeat-containing protein 1A (FWD1), which are functionally interchangeable and thereby considered to function redundantly in the pathway. HOS and FWD1 are both regulated by Wnt/β-catenin signalling, albeit in opposite directions, thus establishing interlocked negative and positive feedback loops. The functional relevance of the opposite regulation of HOS and FWD1 by Wnt/β-catenin signalling in conjunction with their redundant activities in proteasomal degradation of β-catenin remains unresolved. Using a detailed ordinary differential equation model, we investigated the specific influence of each individual feedback mechanism and their combination on Wnt/β-catenin signal transduction under wild-type and cancerous conditions. We found that, under wild-type conditions, the signalling dynamics are predominantly affected by the HOS feedback as a result of a higher concentration of HOS than FWD1. Transcriptional up-regulation of FWD1 by other signalling pathways reduced the impact of the HOS feedback. The opposite regulation of HOS and FWD1 expression by Wnt/β-catenin signalling allows the FWD1 feedback to be employed as a compensation mechanism against aberrant pathway activation as a result of a reduced HOS concentration. By contrast, the FWD1 feedback provides no protection against aberrant activation in adenomatous polyposis coli protein mutant cancer cells.

Article Commentary: Predictive Modeling of Drug Treatment in the Area of Personalized Medicine

Despite a growing body of knowledge on the mechanisms underlying the onset and progression of cancer, treatment success rates in oncology are at best modest. Current approaches use statistical methods that fail to embrace the inherent and expansive complexity of the tumor/patient/drug interaction. Computational modeling, in particular mechanistic modeling, has the power to resolve this complexity. Using fundamental knowledge on the interactions occurring between the components of a complex biological system, large-scale in silico models with predictive capabilities can be generated. Here, we describe how mechanistic virtual patient models, based on systematic molecular characterization of patients and their diseases, have the potential to shift the theranostic paradigm for oncology, both in the fields of personalized medicine and targeted drug development. In particular, we highlight the mechanistic modeling platform ModCell™ for individualized prediction of patient responses to treatment, emphasizing modeling techniques and avenues of application.

An Interaction Library for the FcεRI Signaling Network

Antigen receptors play a central role in adaptive immune responses. Although the molecular networks associated with these receptors have been extensively studied, we currently lack a systems-level understanding of how combinations of non-covalent interactions and post-translational modifications are regulated during signaling to impact cellular decision-making. To fill this knowledge gap, it will be necessary to formalize and piece together information about individual molecular mechanisms to form large-scale computational models of signaling networks. To this end, we have developed an interaction library for signaling by the high-affinity IgE receptor, FcεRI. The library consists of executable rules for protein–protein and protein–lipid interactions. This library extends earlier models for FcεRI signaling and introduces new interactions that have not previously been considered in a model. Thus, this interaction library is a toolkit with which existing models can be expanded and from which new models can be built. As an example, we present models of branching pathways from the adaptor protein Lat, which influence production of the phospholipid PIP₃ at the plasma membrane and the soluble second messenger IP₃. We find that inclusion of a positive feedback loop gives rise to a bistable switch, which may ensure robust responses to stimulation above a threshold level. In addition, the library is visualized to facilitate understanding of network circuitry and identification of network motifs.

Rule-based modeling: a computational approach for studying biomolecular site dynamics in cell signaling systems

Rule-based modeling was developed to address the limitations of traditional approaches for modeling chemical kinetics in cell signaling systems. These systems consist of multiple interacting biomolecules (e.g., proteins), which themselves consist of multiple parts (e.g., domains, linear motifs, and sites of phosphorylation). Consequently, biomolecules that mediate information processing generally have the potential to interact in multiple ways, with the number of possible complexes and posttranslational modification states tending to grow exponentially with the number of binary interactions considered. As a result, only large reaction networks capture all possible consequences of the molecular interactions that occur in a cell signaling system, which is problematic because traditional modeling approaches for chemical kinetics (e.g., ordinary differential equations) require explicit network specification. This problem is circumvented through representation of interactions in terms of local rules. With this approach, network specification is implicit and model specification is concise. Concise representation results in a coarse graining of chemical kinetics, which is introduced because all reactions implied by a rule inherit the rate law associated with that rule. Coarse graining can be appropriate if interactions are modular, and the coarseness of a model can be adjusted as needed. Rules can be specified using specialized model-specification languages, and recently developed tools designed for specification of rule-based models allow one to leverage powerful software engineering capabilities. A rule-based model comprises a set of rules, which can be processed by general-purpose simulation and analysis tools to achieve different objectives (e.g., to perform either a deterministic or stochastic simulation).