Adipocyte signaling, normally and in type 2 diabetes, is far from fully understood. We have earlier developed detailed dynamic mathematical models for several well-studied, partially overlapping, signaling pathways in adipocytes. Still, these models only cover a fraction of the total cellular response. For a broader coverage of the response, large-scale phosphoproteomic data and systems level knowledge on protein interactions are key. However, methods to combine detailed dynamic models with large-scale data, using information about the confidence of included interactions, are lacking. We have developed a method to first establish a core model by connecting existing models of adipocyte cellular signaling for: (1) lipolysis and fatty acid release, (2) glucose uptake, and (3) the release of adiponectin. Next, we use publicly available phosphoproteome data for the insulin response in adipocytes together with prior knowledge on protein interactions, to identify phosphosites downstream of the core model. In a parallel pairwise approach with low computation time, we test whether identified phosphosites can be added to the model. We iteratively collect accepted additions into layers and continue the search for phosphosites downstream of these added layers. For the first 30 layers with the highest confidence (311 added phosphosites), the model predicts independent data well (70-90% correct), and the predictive capability gradually decreases when we add layers of decreasing confidence. In total, 57 layers (3059 phosphosites) can be added to the model with predictive ability kept. Finally, our large-scale, layered model enables dynamic simulations of systems-wide alterations in adipocytes in type 2 diabetes.

In healthy cells, pro- and anti-apoptotic BCL2 family and BH3-only proteins are expressed in a delicate equilibrium. In contrast, this homeostasis is frequently perturbed in cancer cells due to the overexpression of anti-apoptotic BCL2 family proteins. Variability in the expression and sequestration of these proteins in Diffuse Large B cell Lymphoma (DLBCL) likely contributes to variability in response to BH3-mimetics. Successful deployment of BH3-mimetics in DLBCL requires reliable predictions of which lymphoma cells will respond. Here we show that a computational systems biology approach enables accurate prediction of the sensitivity of DLBCL cells to BH3-mimetics. We found that fractional killing of DLBCL, can be explained by cell-to-cell variability in the molecular abundances of signaling proteins. Importantly, by combining protein interaction data with a knowledge of genetic lesions in DLBCL cells, our in silico models accurately predict in vitro response to BH3-mimetics. Furthermore, through virtual DLBCL cells we predict synergistic combinations of BH3-mimetics, which we then experimentally validated. These results show that computational systems biology models of apoptotic signaling, when constrained by experimental data, can facilitate the rational assignment of efficacious targeted inhibitors in B cell malignancies, paving the way for development of more personalized approaches to treatment.

Pancreatic neuroendocrine tumors (PanNETs) are a rare tumor entity with largely unpredictable progression and increasing incidence in developed countries. Molecular pathways involved in PanNETs development are still not elucidated, and specific biomarkers are missing. Moreover, the heterogeneity of PanNETs makes their treatment challenging and most approved targeted therapeutic options for PanNETs lack objective responses. Here, we applied a systems biology approach integrating dynamic modeling strategies, foreign classifier tailored approaches, and patient expression profiles to predict PanNETs progression as well as resistance mechanisms to clinically approved treatments such as the mammalian target of rapamycin complex 1 (mTORC1) inhibitors. We set up a model able to represent frequently reported PanNETs drivers in patient cohorts, such as Menin-1 (MEN1), Death domain associated protein (DAXX), Tuberous Sclerosis (TSC), as well as wild-type tumors. Model-based simulations suggested drivers of cancer progression as both first and second hits after MEN1 loss. In addition, we could predict the benefit of mTORC1 inhibitors on differentially mutated cohorts and hypothesize resistance mechanisms. Our approach sheds light on a more personalized prediction and treatment of PanNET mutant phenotypes.

Cancer is widely considered a genetic disease. Notably, recent works have highlighted that every human gene may possibly be associated with cancer. Thus, the distinction between genes that drive oncogenesis and those that are associated to the disease, but do not play a role, requires attention. Here we investigated single cells and bulk (cell-population) datasets of several cancer transcriptomes and proteomes in relation to their healthy counterparts. When analyzed by machine learning and statistical approaches in bulk datasets, both general and cancer-specific oncogenes, as defined by the Cancer Genes Census, show invariant behavior to randomly selected gene sets of the same size for all cancers. However, when protein-protein interaction analyses were performed, the oncogenes-derived networks show higher connectivity than those relative to random genes. Moreover, at single-cell scale, we observe variant behavior in a subset of oncogenes for each considered cancer type. Moving forward, we concur that the role of oncogenes needs to be further scrutinized by adopting protein causality and higher-resolution single-cell analyses.

Increasing evidence points to a role of the circadian clock in the regulation of cancer hallmarks with a strong impact on the understanding and treatment of this disease. Anti-cancer treatment can be personalized considering treatment timing. Here we present a new mathematical model based on data from three colorectal cancer cell lines and core-clock knock-outs, which couples the circadian and drug metabolism network, and that allows to determine toxicity profiles for a given drug and cell type. Moreover, this model integrates external Zeitgebers and thus may be used to fine-tune toxicity by using external factors, such as light, and therefore, to a certain extent, help fitting the endogenous rhythms of the patients to a defined clinic routine facilitating the implementation of time-dependent treatment in clinical practice.

Dimensionality reduction offers unique insights into high-dimensional microbiome dynamics by leveraging collective abundance fluctuations of multiple bacteria driven by similar ecological perturbations. However, methods providing lower-dimensional representations of microbiome dynamics both at the community and individual taxa levels are not currently available. To that end, we present EMBED: Essential MicroBiomE Dynamics, a probabilistic nonlinear tensor factorization approach. Like normal mode analysis in structural biophysics, EMBED infers ecological normal modes (ECNs), which represent the unique orthogonal modes capturing the collective behavior of microbial communities. Using multiple real and synthetic datasets, we show that a very small number of ECNs can accurately approximate microbiome dynamics. Inferred ECNs reflect specific ecological behaviors, providing natural templates along which the dynamics of individual bacteria may be partitioned. Moreover, the multi-subject treatment in EMBED systematically identifies subject-specific and universal abundance dynamics that are not detected by traditional approaches. Collectively, these results highlight the utility of EMBED as a versatile dimensionality reduction tool for studies of microbiome dynamics.

Microridges are evolutionarily conserved actin-rich protrusions present on the apical surface of squamous epithelial cells. In zebrafish epidermal cells, microridges form self-evolving patterns due to the underlying actomyosin network dynamics. However, their morphological and dynamic characteristics have remained poorly understood owing to a lack of computational methods. We achieved ~95% pixel-level accuracy with a deep learning microridge segmentation strategy enabling quantitative insights into their bio-physical-mechanical characteristics. From the segmented images, we estimated an effective microridge persistence length of ~6.1 μm. We discovered the presence of mechanical fluctuations and found relatively greater stresses stored within patterns of yolk than flank, indicating distinct regulation of their actomyosin networks. Furthermore, spontaneous formations and positional fluctuations of actin clusters within microridges were associated with pattern rearrangements over short length/time-scales. Our framework allows large-scale spatiotemporal analysis of microridges during epithelial development and probing of their responses to chemical and genetic perturbations to unravel the underlying patterning mechanisms.

Mathematical models based on homogenized representation of cardiac tissue have greatly improved our understanding of cardiac electrophysiology. However, these models are too coarse to investigate the dynamics at the level of the myocytes since the cells are not present in homogenized models. Recently, fine scale models have been proposed to allow for cell-level resolution of the dynamics, but these models are too computationally expensive to be used in applications like whole heart simulations of large animals. To address this issue, we propose a model that balances computational demands and physiological accuracy. The model is founded on Kirchhoff's current law, and represents every myocyte in the tissue. This allows specific properties to be assigned to individual cardiomyocytes, and other cell types like fibroblasts can be added to the model in an accurate manner while keeping the computing efforts reasonable.