Disentangling the Complexity of HGF Signaling by Combining Qualitative and Quantitative Modeling

Basal MET phosphorylation is an indicator of hepatocyte dysregulation in liver disease

Chronic liver diseases are worldwide on the rise. Due to the rapidly increasing incidence, in particular in Western countries, metabolic dysfunction-associated steatotic liver disease (MASLD) is gaining importance as the disease can develop into hepatocellular carcinoma. Lipid accumulation in hepatocytes has been identified as the characteristic structural change in MASLD development, but molecular mechanisms responsible for disease progression remained unresolved. Here, we uncover in primary hepatocytes from a preclinical model fed with a Western diet (WD) an increased basal MET phosphorylation and a strong downregulation of the PI3K-AKT pathway. Dynamic pathway modeling of hepatocyte growth factor (HGF) signal transduction combined with global proteomics identifies that an elevated basal MET phosphorylation rate is the main driver of altered signaling leading to increased proliferation of WD-hepatocytes. Model-adaptation to patient-derived hepatocytes reveal patient-specific variability in basal MET phosphorylation, which correlates with patient outcome after liver surgery. Thus, dysregulated basal MET phosphorylation could be an indicator for the health status of the liver and thereby inform on the risk of a patient to suffer from liver failure after surgery.

Deciphering signal transduction networks in the liver by mechanistic mathematical modelling

In health and disease, liver cells are continuously exposed to cytokines and growth factors. While individual signal transduction pathways induced by these factors were studied in great detail, the cellular responses induced by repeated or combined stimulations are complex and less understood. Growth factor receptors on the cell surface of hepatocytes were shown to be regulated by receptor interactions, receptor trafficking and feedback regulation. Here, we exemplify how mechanistic mathematical modelling based on quantitative data can be employed to disentangle these interactions at the molecular level. Crucial is the analysis at a mechanistic level based on quantitative longitudinal data within a mathematical framework. In such multi-layered information, step-wise mathematical modelling using submodules is of advantage, which is fostered by sharing of standardized experimental data and mathematical models. Integration of signal transduction with metabolic regulation in the liver and mechanistic links to translational approaches promise to provide predictive tools for biology and personalized medicine.

Systems biology of angiogenesis signaling: Computational models and omics

Angiogenesis is a highly regulated multiscale process that involves a plethora of cells, their cellular signal transduction, activation, proliferation, differentiation, as well as their intercellular communication. The coordinated execution and integration of such complex signaling programs is critical for physiological angiogenesis to take place in normal growth, development, exercise, and wound healing, while its dysregulation is critically linked to many major human diseases such as cancer, cardiovascular diseases, and ocular disorders; it is also crucial in regenerative medicine. Although huge efforts have been devoted to drug development for these diseases by investigation of angiogenesis‐targeted therapies, only a few therapeutics and targets have proved effective in humans due to the innate multiscale complexity and nonlinearity in the process of angiogenic signaling. As a promising approach that can help better address this challenge, systems biology modeling allows the integration of knowledge across studies and scales and provides a powerful means to mechanistically elucidate and connect the individual molecular and cellular signaling components that function in concert to regulate angiogenesis. In this review, we summarize and discuss how systems biology modeling studies, at the pathway‐, cell‐, tissue‐, and whole body‐levels, have advanced our understanding of signaling in angiogenesis and thereby delivered new translational insights for human diseases.

This article is categorized under:

Cardiovascular Diseases > Computational Models

Cancer > Computational Models

TopoFilter: a MATLAB package for mechanistic model identification in systems biology

Background

To develop mechanistic dynamic models in systems biology, one often needs to identify all (or minimal) representations of the biological processes that are consistent with experimental data, out of a potentially large set of hypothetical mechanisms. However, a simple enumeration of all alternatives becomes quickly intractable when the number of model parameters grows. Selecting appropriate dynamic models out of a large ensemble of models, taking the uncertainty in our biological knowledge and in the experimental data into account, is therefore a key current problem in systems biology.

Results

The TopoFilter package addresses this problem in a heuristic and automated fashion by implementing the previously described topological filtering method for Bayesian model selection. It includes a core heuristic for searching the space of submodels of a parametrized model, coupled with a sampling-based exploration of the parameter space. Recent developments of the method allow to balance exhaustiveness and speed of the model space search, to efficiently re-sample parameters, to parallelize the search, and to use custom scoring functions. We use a theoretical example to motivate these features and then demonstrate TopoFilter’s applicability for a yeast signaling network with more than 250’000 possible model structures.

Conclusions

TopoFilter is a flexible software framework that makes Bayesian model selection and reduction efficient and scalable to network models of a complexity that represents contemporary problems in, for example, cell signaling. TopoFilter is open-source, available under the GPL-3.0 license at https://gitlab.com/csb.ethz/TopoFilter. It includes installation instructions, a quickstart guide, a description of all package options, and multiple examples.

Mechanistically detailed systems biology modeling of the HGF/Met pathway in hepatocellular carcinoma

Hepatocyte growth factor (HGF) signaling through its receptor Met has been implicated in hepatocellular carcinoma tumorigenesis and progression. Met interaction with integrins is shown to modulate the downstream signaling to Akt and ERK (extracellular-regulated kinase). In this study, we developed a mechanistically detailed systems biology model of HGF/Met signaling pathway that incorporated specific interactions with integrins to investigate the efficacy of integrin-binding peptide, AXT050, as monotherapy and in combination with other therapeutics targeting this pathway. Here we report that the modeled dynamics of the response to AXT050 revealed that receptor trafficking is sufficient to explain the effect of Met-integrin interactions on HGF signaling. Furthermore, the model predicted patient-specific synergy and antagonism of efficacy and potency for combination of AXT050 with sorafenib, cabozantinib, and rilotumumab. Overall, the model provides a valuable framework for studying the efficacy of drugs targeting receptor tyrosine kinase interaction with integrins, and identification of synergistic drug combinations for the patients.

Coupled feedback regulation of nuclear factor of activated T-cells (NFAT) modulates activation-induced cell death of T cells

A properly functioning immune system is vital for an organism’s wellbeing. Immune tolerance is a critical feature of the immune system that allows immune cells to mount effective responses against exogenous pathogens such as viruses and bacteria, while preventing attack to self-tissues. Activation-induced cell death (AICD) in T lymphocytes, in which repeated stimulations of the T-cell receptor (TCR) lead to activation and then apoptosis of T cells, is a major mechanism for T cell homeostasis and helps maintain peripheral immune tolerance. Defects in AICD can lead to development of autoimmune diseases. Despite its importance, the regulatory mechanisms that underlie AICD remain poorly understood, particularly at an integrative network level. Here, we develop a dynamic multi-pathway model of the integrated TCR signalling network and perform model-based analysis to characterize the network-level properties of AICD. Model simulation and analysis show that amplified activation of the transcriptional factor NFAT in response to repeated TCR stimulations, a phenomenon central to AICD, is tightly modulated by a coupled positive-negative feedback mechanism. NFAT amplification is predominantly enabled by a positive feedback self-regulated by NFAT, while opposed by a NFAT-induced negative feedback via Carabin. Furthermore, model analysis predicts an optimal therapeutic window for drugs that help minimize proliferation while maximize AICD of T cells. Overall, our study provides a comprehensive mathematical model of TCR signalling and model-based analysis offers new network-level insights into the regulation of activation-induced cell death in T cells.

Computational Modeling in Liver Surgery

The need for extended liver resection is increasing due to the growing incidence of liver tumors in aging societies. Individualized surgical planning is the key for identifying the optimal resection strategy and to minimize the risk of postoperative liver failure and tumor recurrence. Current computational tools provide virtual planning of liver resection by taking into account the spatial relationship between the tumor and the hepatic vascular trees, as well as the size of the future liver remnant. However, size and function of the liver are not necessarily equivalent. Hence, determining the future liver volume might misestimate the future liver function, especially in cases of hepatic comorbidities such as hepatic steatosis. A systems medicine approach could be applied, including biological, medical, and surgical aspects, by integrating all available anatomical and functional information of the individual patient. Such an approach holds promise for better prediction of postoperative liver function and hence improved risk assessment. This review provides an overview of mathematical models related to the liver and its function and explores their potential relevance for computational liver surgery. We first summarize key facts of hepatic anatomy, physiology, and pathology relevant for hepatic surgery, followed by a description of the computational tools currently used in liver surgical planning. Then we present selected state-of-the-art computational liver models potentially useful to support liver surgery. Finally, we discuss the main challenges that will need to be addressed when developing advanced computational planning tools in the context of liver surgery.

Protein abundance of AKT and ERK pathway components governs cell type-specific regulation of proliferation

Signaling through the AKT and ERK pathways controls cell proliferation. However, the integrated regulation of this multistep process, involving signal processing, cell growth and cell cycle progression, is poorly understood. Here, we study different hematopoietic cell types, in which AKT and ERK signaling is triggered by erythropoietin (Epo). Although these cell types share the molecular network topology for pro‐proliferative Epo signaling, they exhibit distinct proliferative responses. Iterating quantitative experiments and mathematical modeling, we identify two molecular sources for cell type‐specific proliferation. First, cell type‐specific protein abundance patterns cause differential signal flow along the AKT and ERK pathways. Second, downstream regulators of both pathways have differential effects on proliferation, suggesting that protein synthesis is rate‐limiting for faster cycling cells while slower cell cycles are controlled at the G1‐S progression. The integrated mathematical model of Epo‐driven proliferation explains cell type‐specific effects of targeted AKT and ERK inhibitors and faithfully predicts, based on the protein abundance, anti‐proliferative effects of inhibitors in primary human erythroid progenitor cells. Our findings suggest that the effectiveness of targeted cancer therapy might become predictable from protein abundance.

Clustering and Network Analysis of Reverse Phase Protein Array Data

Molecular profiling of proteins and phosphoproteins using a reverse phase protein array (RPPA) platform, with a panel of target-specific antibodies, enables the parallel, quantitative proteomic analysis of many biological samples in a microarray format. Hence, RPPA analysis can generate a high volume of multidimensional data that must be effectively interrogated and interpreted. A range of computational techniques for data mining can be applied to detect and explore data structure and to form functional predictions from large datasets. Here, two approaches for the computational analysis of RPPA data are detailed: the identification of similar patterns of protein expression by hierarchical cluster analysis and the modeling of protein interactions and signaling relationships by network analysis. The protocols use freely available, cross-platform software, are easy to implement, and do not require any programming expertise. Serving as data-driven starting points for further in-depth analysis, validation, and biological experimentation, these and related bioinformatic approaches can accelerate the functional interpretation of RPPA data.

In silico model-based inference: an emerging approach for inverse problems in engineering better medicines

Identifying the network of biochemical interactions that underpin disease pathophysiology is a key hurdle in drug discovery. While many components involved in these biological processes are identified, how components organize differently in health and disease remains unclear. In chemical engineering, mechanistic modeling provides a quantitative framework to capture our understanding of a reactive system and test this knowledge against data. Here, we describe an emerging approach to test this knowledge against data that leverages concepts from probability, Bayesian statistics, and chemical kinetics by focusing on two related inverse problems. The first problem is to identify the causal structure of the reaction network, given uncertainty as to how the reactive components interact. The second problem is to identify the values of the model parameters, when a network is known a priori.