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.

Activator-blocker model of transcriptional regulation by pioneer-like factors

Zygotic genome activation (ZGA) in the development of flies, fish, frogs and mammals depends on pioneer-like transcription factors (TFs). Those TFs create open chromatin regions, promote histone acetylation on enhancers, and activate transcription. Here, we use the panel of single, double and triple mutants for zebrafish genome activators Pou5f3, Sox19b and Nanog, multi-omics and mathematical modeling to investigate the combinatorial mechanisms of genome activation. We show that Pou5f3 and Nanog act differently on synergistic and antagonistic enhancer types. Pou5f3 and Nanog both bind as pioneer-like TFs on synergistic enhancers, promote histone acetylation and activate transcription. Antagonistic enhancers are activated by binding of one of these factors. The other TF binds as non-pioneer-like TF, competes with the activator and blocks all its effects, partially or completely. This activator-blocker mechanism mutually restricts widespread transcriptional activation by Pou5f3 and Nanog and prevents premature expression of late developmental regulators in the early embryo.

Kinetics of mRNA nuclear export regulate innate immune response gene expression

The abundance and stimulus-responsiveness of mature mRNA is thought to be determined by nuclear synthesis, processing, and cytoplasmic decay. However, the rate and efficiency of moving mRNA to the cytoplasm almost certainly contributes, but has rarely been measured. Here, we investigated mRNA export rates for innate immune genes. We generated high spatio-temporal resolution RNA-seq data from endotoxin-stimulated macrophages and parameterized a mathematical model to infer kinetic parameters with confidence intervals. We find that the effective chromatin-to-cytoplasm export rate is gene-specific, varying 100-fold: for some genes, less than 5% of synthesized transcripts arrive in the cytoplasm as mature mRNAs, while others show high export efficiency. Interestingly, effective export rates do not determine temporal gene responsiveness, but complement the wide range of mRNA decay rates; this ensures similar abundances of short- and long-lived mRNAs, which form successive innate immune response expression waves.

Disentangling ERBB Signaling in Breast Cancer Subtypes-A Model-Based Analysis

Simple Summary

Breast cancer subtypes are characterized by the expression and activity of estrogen-, progesterone- and HER2-receptors and differ by the treatment as well as patient prognosis. Tumors of the HER2-subtype overexpress this receptor and are successfully targeted with anti-HER2 therapies. We wanted to know if the HER2-receptor and the downstream signaling network act similarly also in the other subtypes and if this network could potentially be a therapeutic target beyond the HER2-positive subtype. To this end, we quantitatively assessed the wiring of signaling events in the individual subtypes to unravel the characteristics of HER-signaling. Our data along with a model-based analysis suggest that major parts of the intracellular signal transduction network are unchanged between the different breast cancer subtypes and that the clinical differences mostly come from the different levels at which these receptors are present in tumor cells as well as from the particular mutations that are present in individual tumors.

Abstract

Targeted therapies have shown striking success in the treatment of cancer over the last years. However, their specific effects on an individual tumor appear to be varying and difficult to predict. Using an integrative modeling approach that combines mechanistic and regression modeling, we gained insights into the response mechanisms of breast cancer cells due to different ligand–drug combinations. The multi-pathway model, capturing ERBB receptor signaling as well as downstream MAPK and PI3K pathways was calibrated on time-resolved data of the luminal breast cancer cell lines MCF7 and T47D across an array of four ligands and five drugs. The same model was then successfully applied to triple negative and HER2-positive breast cancer cell lines, requiring adjustments mostly for the respective receptor compositions within these cell lines. The additional relevance of cell-line-specific mutations in the MAPK and PI3K pathway components was identified via L₁ regularization, where the impact of these mutations on pathway activation was uncovered. Finally, we predicted and experimentally validated the proliferation response of cells to drug co-treatments. We developed a unified mathematical model that can describe the ERBB receptor and downstream signaling in response to therapeutic drugs targeting this clinically relevant signaling network in cell line that represent three major subtypes of breast cancer. Our data and model suggest that alterations in this network could render anti-HER therapies relevant beyond the HER2-positive subtype.

Dynamic modeling of Nrf2 pathway activation in liver cells after toxicant exposure

Cells are exposed to oxidative stress and reactive metabolites every day. The Nrf2 signaling pathway responds to oxidative stress by upregulation of antioxidants like glutathione (GSH) to compensate the stress insult and re-establish homeostasis. Although mechanisms describing the interaction between the key pathway constituents Nrf2, Keap1 and p62 are widely reviewed and discussed in literature, quantitative dynamic models bringing together these mechanisms with time-resolved data are limited. Here, we present an ordinary differential equation (ODE) based dynamic model to describe the dynamic response of Nrf2, Keap1, Srxn1 and GSH to oxidative stress caused by the soft-electrophile diethyl maleate (DEM). The time-resolved data obtained by single-cell confocal microscopy of green fluorescent protein (GFP) reporters and qPCR of the Nrf2 pathway components complemented with siRNA knock down experiments, is accurately described by the calibrated mathematical model. We show that the quantitative model can describe the activation of the Nrf2 pathway by compounds with a different mechanism of activation, including drugs which are known for their ability to cause drug induced liver-injury (DILI) i.e., diclofenac (DCF) and omeprazole (OMZ). Finally, we show that our model can reveal differences in the processes leading to altered activation dynamics amongst DILI inducing drugs.

Clustering Based Parameter Estimation of Thyroid Hormone Pathway

An ordinary differential equation (ODE) model of the working of the thyroid system for euthyroidism has been presented. As clinical data for thyroid hormones is relatively scarce, such modelling offers potential benefits over wet lab procedures. Genetic algorithms developed for determining of parameters of the ODE system using the available data have been presented and evaluated. This approach enables subject specific parameter estimation towards characterisation of individual thyroid operation. Initially, a simple steady state model was used. Later a cosinor model for the circadian variation of thyroid hormones was used to obtain more reliable results, as indicated through sensitivity analysis in conjunction with other statistical methods. Our parameter determination method has been tested on groups of patients with similar observed values of thyroid stimulating hormone (TSH), free T ₃ and free T ₄ (identified through clustering) to determine their parameter values jointly. This approach appears to produce parameter sets with lower variation than parameters determined independently, thus leading to better parameter determination.

Evaluation framework for systems models

As decisions in drug development increasingly rely on predictions from mechanistic systems models, assessing the predictive capability of such models is becoming more important. Several frameworks for the development of quantitative systems pharmacology (QSP) models have been proposed. In this paper, we add to this body of work with a framework that focuses on the appropriate use of qualitative and quantitative model evaluation methods. We provide details and references for those wishing to apply these methods, which include sensitivity and identifiability analyses, as well as concepts such as validation and uncertainty quantification. Many of these methods have been used successfully in other fields, but are not as common in QSP modeling. We illustrate how to apply these methods to evaluate QSP models, and propose methods to use in two case studies. We also share examples of misleading results when inappropriate analyses are used.

Magnitude of Ubiquitination Determines the Fate of Epidermal Growth Factor Receptor Upon Ligand Stimulation

Receptor tyrosine kinases (RTK) bind growth factors and are critical for cell proliferation and differentiation. Their dysregulation leads to a loss of growth control, often resulting in cancer. Epidermal growth factor receptor (EGFR) is the prototypic RTK and can bind several ligands exhibiting distinct mitogenic potentials. Whereas the phosphorylation on individual EGFR sites and their roles for downstream signaling have been extensively studied, less is known about ligand-specific ubiquitination events on EGFR, which are crucial for signal attenuation and termination. We used a proteomics-based workflow for absolute quantitation combined with mathematical modeling to unveil potentially decisive ubiquitination events on EGFR from the first 30 seconds to 15 minutes of stimulation. Four ligands were used for stimulation: epidermal growth factor (EGF), heparin-binding-EGF like growth factor, transforming growth factor-α and epiregulin. Whereas only little differences in the order of individual ubiquitination sites were observed, the overall amount of modified receptor differed depending on the used ligand, indicating that absolute magnitude of EGFR ubiquitination, and not distinctly regulated ubiquitination sites, is a major determinant for signal attenuation and the subsequent cellular outcomes.

PEtab-Interoperable specification of parameter estimation problems in systems biology

Reproducibility and reusability of the results of data-based modeling studies are essential. Yet, there has been-so far-no broadly supported format for the specification of parameter estimation problems in systems biology. Here, we introduce PEtab, a format which facilitates the specification of parameter estimation problems using Systems Biology Markup Language (SBML) models and a set of tab-separated value files describing the observation model and experimental data as well as parameters to be estimated. We already implemented PEtab support into eight well-established model simulation and parameter estimation toolboxes with hundreds of users in total. We provide a Python library for validation and modification of a PEtab problem and currently 20 example parameter estimation problems based on recent studies.

Confidence intervals by constrained optimization-An algorithm and software package for practical identifiability analysis in systems biology

Practical identifiability of Systems Biology models has received a lot of attention in recent scientific research. It addresses the crucial question for models’ predictability: how accurately can the models’ parameters be recovered from available experimental data. The methods based on profile likelihood are among the most reliable methods of practical identification. However, these methods are often computationally demanding or lead to inaccurate estimations of parameters’ confidence intervals. Development of methods, which can accurately produce parameters’ confidence intervals in reasonable computational time, is of utmost importance for Systems Biology and QSP modeling.

We propose an algorithm Confidence Intervals by Constraint Optimization (CICO) based on profile likelihood, designed to speed-up confidence intervals estimation and reduce computational cost. The numerical implementation of the algorithm includes settings to control the accuracy of confidence intervals estimates. The algorithm was tested on a number of Systems Biology models, including Taxol treatment model and STAT5 Dimerization model, discussed in the current article.

The CICO algorithm is implemented in a software package freely available in Julia (https://github.com/insysbio/LikelihoodProfiler.jl) and Python (https://github.com/insysbio/LikelihoodProfiler.py).

Differential equations-based models are widely used in Systems Biology and Quantitative Systems Pharmacology and play a significant role in the discovery of new disease-directed drugs. Complexity of models is a trade off from their employment to crucial fields of biology and medicine. These areas of application require large non-linear models with many unknown parameters. How accurately can the parameters of a model be recovered from experimental data? What is the identifiable subset of parameters? Can the model be reduced or reparameterized to become identifiable? All those questions of identifiability analysis are essential for model’s predictability and reliability. That explains why the topic of identifiability of Systems Biology models has received a lot of attention in recent scientific research. However, existing numerical methods of identifiability analysis are computationally demanding or often lead to inaccurate estimations. Development of methods, which can accurately produce parameters’ confidence intervals in reasonable computational time, is of utmost importance for Systems Biology and QSP modeling. We propose an algorithm and a software package to test identifiability of Systems Biology models, designed to speed-up confidence intervals estimation and reduce computational cost. The software package was tested on a number of Systems Biology models, including Taxol treatment model and STAT5 Dimerization model, discussed in the current article.

Should a viral genome stay in the host cell or leave? A quantitative dynamics study of how hepatitis C virus deals with this dilemma

Virus proliferation involves gene replication inside infected cells and transmission to new target cells. Once positive-strand RNA virus has infected a cell, the viral genome serves as a template for copying (“stay-strategy”) or is packaged into a progeny virion that will be released extracellularly (“leave-strategy”). The balance between genome replication and virion release determines virus production and transmission efficacy. The ensuing trade-off has not yet been well characterized. In this study, we use hepatitis C virus (HCV) as a model system to study the balance of the two strategies. Combining viral infection cell culture assays with mathematical modeling, we characterize the dynamics of two different HCV strains (JFH-1, a clinical isolate, and Jc1-n, a laboratory strain), which have different viral release characteristics. We found that 0.63% and 1.70% of JFH-1 and Jc1-n intracellular viral RNAs, respectively, are used for producing and releasing progeny virions. Analysis of the Malthusian parameter of the HCV genome (i.e., initial proliferation rate) and the number of de novo infections (i.e., initial transmissibility) suggests that the leave-strategy provides a higher level of initial transmission for Jc1-n, whereas, in contrast, the stay-strategy provides a higher initial proliferation rate for JFH-1. Thus, theoretical-experimental analysis of viral dynamics enables us to better understand the proliferation strategies of viruses, which contributes to the efficient control of virus transmission. Ours is the first study to analyze the stay-leave trade-off during the viral life cycle and the significance of the replication-release switching mechanism for viral proliferation.

A theoretical-experimental analysis of viral dynamics reveals a stay/leave trade-off during the viral life cycle and demonstrates the significance for viral proliferation of a replication-release switching mechanism.

Disentangling molecular mechanisms regulating sensitization of interferon alpha signal transduction

Tightly interlinked feedback regulators control the dynamics of intracellular responses elicited by the activation of signal transduction pathways. Interferon alpha (IFNα) orchestrates antiviral responses in hepatocytes, yet mechanisms that define pathway sensitization in response to prestimulation with different IFNα doses remained unresolved. We establish, based on quantitative measurements obtained for the hepatoma cell line Huh7.5, an ordinary differential equation model for IFNα signal transduction that comprises the feedback regulators STAT1, STAT2, IRF9, USP18, SOCS1, SOCS3, and IRF2. The model-based analysis shows that, mediated by the signaling proteins STAT2 and IRF9, prestimulation with a low IFNα dose hypersensitizes the pathway. In contrast, prestimulation with a high dose of IFNα leads to a dose-dependent desensitization, mediated by the negative regulators USP18 and SOCS1 that act at the receptor. The analysis of basal protein abundance in primary human hepatocytes reveals high heterogeneity in patient-specific amounts of STAT1, STAT2, IRF9, and USP18. The mathematical modeling approach shows that the basal amount of USP18 determines patient-specific pathway desensitization, while the abundance of STAT2 predicts the patient-specific IFNα signal response.

Comparative quantitative systems pharmacology modeling of anti-PCSK9 therapeutic modalities in hypercholesterolemia

Since the discovery of proprotein convertase subtilisin/kexin type 9 (PCSK9) as an attractive target in the treatment of hypercholesterolemia, multiple anti-PCSK9 therapeutic modalities have been pursued in drug development. The objective of this research is to set the stage for the quantitative benchmarking of two anti-PCSK9 pharmacological modality classes, monoclonal antibodies (mAbs) and small interfering RNA (siRNA). To this end, we developed an integrative mathematical model of lipoprotein homeostasis describing the dynamic interplay between PCSK9, LDL-cholesterol (LDL-C), VLDL-cholesterol, HDL-cholesterol (HDL-C), apoB, lipoprotein a [Lp(a)], and triglycerides (TGs). We demonstrate that LDL-C decreased proportionally to PCSK9 reduction for both mAb and siRNA modalities. At marketed doses, however, treatment with mAbs resulted in an additional ∼20% LDL-C reduction compared with siRNA. We further used the model as an evaluation tool and determined that no quantitative differences were observed in HDL-C, Lp(a), TG, or apoB responses, suggesting that the disruption of PCSK9 synthesis would provide no additional effects on lipoprotein-related biomarkers in the patient segment investigated. Predictive model simulations further indicate that siRNA therapies may reach reductions in LDL-C levels comparable to those achieved with mAbs if the current threshold of 80% PCSK9 inhibition via siRNA could be overcome.