Statistical ensemble of gene regulatory networks of macrophage differentiation

Pancreatic cancer mutationscape: revealing the link between modular restructuring and intervention efficacy amidst common mutations

There is increasing evidence that biological systems are modular in both structure and function. Complex biological signaling networks such as gene regulatory networks (GRNs) are proving to be composed of subcategories that are interconnected and hierarchically ranked. These networks contain highly dynamic processes that ultimately dictate cellular function over time, as well as influence phenotypic fate transitions. In this work, we use a stochastic multicellular signaling network of pancreatic cancer (PC) to show that the variance in topological rankings of the most phenotypically influential modules implies a strong relationship between structure and function. We further show that induction of mutations alters the modular structure, which analogously influences the aggression and controllability of the disease in silico. We finally present evidence that the impact and location of mutations with respect to PC modular structure directly corresponds to the efficacy of single agent treatments in silico, because topologically deep mutations require deep targets for control.

Macrophage phenotype transitions in a stochastic gene-regulatory network model

Polarization is the process by which a macrophage cell commits to a phenotype based on external signal stimulation. To know how this process is affected by random fluctuations and events within a cell is of utmost importance to better understand the underlying dynamics and predict possible phenotype transitions. For this purpose, we develop a stochastic modeling approach for the macrophage polarization process. We classify phenotype states using the Robust Perron Cluster Analysis and quantify transition pathways and probabilities by applying Transition Path Theory. Depending on the model parameters, we identify four bistable and one tristable phenotype configuration. We find that bistable transitions are fast but their states less robust. In contrast, phenotype transitions in the tristable situation have a comparatively long time duration, which reflects the robustness of the states. The results indicate parallels in the overall transition behavior of macrophage cells with other heterogeneous and plastic cell types, such as cancer cells. Our approach allows for a probabilistic interpretation of macrophage phenotype transitions and biological inference on phenotype robustness. In general, the methodology can easily be adapted to other systems where random state switches are known to occur.

A multiscale hybrid model for exploring the effect of Resolvin D1 on macrophage polarization during acute inflammation

Dysregulated inflammation underlies various diseases. Specialized pro-resolving mediators (SPMs) like Resolvin D1 (RvD1) have been shown to resolve inflammation and halt disease progression. Macrophages, key immune cells that drive inflammation, respond to the presence of RvD1 by polarizing to an anti-inflammatory type (M₂). However, RvD1's mechanisms, roles, and utility are not fully understood. This paper introduces a gene-regulatory network (GRN) model that contains pathways for RvD1 and other SPMs and proinflammatory molecules like lipopolysaccharides. We couple this GRN model to a partial differential equation - agent-based hybrid model using a multiscale framework to simulate an acute inflammatory response with and without the presence of RvD1. We calibrate and validate the model using experimental data from two animal models. The model reproduces the dynamics of key immune components and the effects of RvD1 during acute inflammation. Our results suggest RvD1 can drive macrophage polarization through the G protein-coupled receptor 32 (GRP32) pathway. The presence of RvD1 leads to an earlier and increased M₂ polarization, reduced neutrophil recruitment, and faster apoptotic neutrophil clearance. These results support a body of literature that suggests that RvD1 is a promising candidate for promoting the resolution of acute inflammation. We conclude that once calibrated and validated on human data, the model can identify critical sources of uncertainty, which could be further elucidated in biological experiments and assessed for clinical use.

Expression Profile of Genes Related to the Th17 Pathway in Macrophages Infected by Leishmania major and Leishmania amazonensis: The Use of Gene Regulatory Networks in Modeling This Pathway

Leishmania amazonensis and Leishmania major are the causative agents of cutaneous and mucocutaneous diseases. The infections‘ outcome depends on host–parasite interactions and Th1/Th2 response, and in cutaneous form, regulation of Th17 cytokines has been reported to maintain inflammation in lesions. Despite that, the Th17 regulatory scenario remains unclear. With the aim to gain a better understanding of the transcription factors (TFs) and genes involved in Th17 induction, in this study, the role of inducing factors of the Th17 pathway in Leishmania–macrophage infection was addressed through computational modeling of gene regulatory networks (GRNs). The Th17 GRN modeling integrated experimentally validated data available in the literature and gene expression data from a time-series RNA-seq experiment (4, 24, 48, and 72 h post-infection). The generated model comprises a total of 10 TFs, 22 coding genes, and 16 cytokines related to the Th17 immune modulation. Addressing the Th17 induction in infected and uninfected macrophages, an increase of 2- to 3-fold in 4–24 h was observed in the former. However, there was a decrease in basal levels at 48–72 h for both groups. In order to evaluate the possible outcomes triggered by GRN component modulation in the Th17 pathway. The generated GRN models promoted an integrative and dynamic view of Leishmania–macrophage interaction over time that extends beyond the analysis of single-gene expression.

From Infection to Immunity: Understanding the Response to SARS-CoV2 Through In-Silico Modeling

Background

Immune system conditions of the patient is a key factor in COVID-19 infection survival. A growing number of studies have focused on immunological determinants to develop better biomarkers for therapies.

Aim

Studies of the insurgence of immunity is at the core of both SARS-CoV-2 vaccine development and therapies. This paper attempts to describe the insurgence (and the span) of immunity in COVID-19 at the population level by developing an in-silico model. We simulate the immune response to SARS-CoV-2 and analyze the impact of infecting viral load, affinity to the ACE2 receptor, and age in an artificially infected population on the course of the disease.

Methods

We use a stochastic agent-based immune simulation platform to construct a virtual cohort of infected individuals with age-dependent varying degrees of immune competence. We use a parameter set to reproduce known inter-patient variability and general epidemiological statistics.

Results

By assuming the viremia at day 30 of the infection to be the proxy for lethality, we reproduce in-silico several clinical observations and identify critical factors in the statistical evolution of the infection. In particular, we evidence the importance of the humoral response over the cytotoxic response and find that the antibody titers measured after day 25 from the infection are a prognostic factor for determining the clinical outcome of the infection. Our modeling framework uses COVID-19 infection to demonstrate the actionable effectiveness of modeling the immune response at individual and population levels. The model developed can explain and interpret observed patterns of infection and makes verifiable temporal predictions. Within the limitations imposed by the simulated environment, this work proposes quantitatively that the great variability observed in the patient outcomes in real life can be the mere result of subtle variability in the infecting viral load and immune competence in the population. In this work, we exemplify how computational modeling of immune response provides an important view to discuss hypothesis and design new experiments, in particular paving the way to further investigations about the duration of vaccine-elicited immunity especially in the view of the blundering effect of immunosenescence.

Potential predictors of type-2 diabetes risk: machine learning, synthetic data and wearable health devices

BACKGROUND

The aim of a recent research project was the investigation of the mechanisms involved in the onset of type 2 diabetes in the absence of familiarity. This has led to the development of a computational model that recapitulates the aetiology of the disease and simulates the immunological and metabolic alterations linked to type-2 diabetes subjected to clinical, physiological, and behavioural features of prototypical human individuals.

RESULTS

We analysed the time course of 46,170 virtual subjects, experiencing different lifestyle conditions. We then set up a statistical model able to recapitulate the simulated outcomes.

CONCLUSIONS

The resulting machine learning model adequately predicts the synthetic dataset and can, therefore, be used as a computationally-cheaper version of the detailed mathematical model, ready to be implemented on mobile devices to allow self-assessment by informed and aware individuals. The computational model used to generate the dataset of this work is available as a web-service at the following address: http://kraken.iac.rm.cnr.it/T2DM .

Bifurcation and sensitivity analysis reveal key drivers of multistability in a model of macrophage polarization

In this paper, we present and analyze a mathematical model for polarization of a single macrophage which, despite its simplicity, exhibits complex dynamics in terms of multistability. In particular, we demonstrate that an asymmetry in the regulatory mechanisms and parameter values is important for observing multiple phenotypes. Bifurcation and sensitivity analyses show that external signaling cues are necessary for macrophage commitment and emergence to a phenotype, but that the intrinsic macrophage pathways are equally important. Based on our numerical results, we formulate hypotheses that could be further investigated by laboratory experiments to deepen our understanding of macrophage polarization.

A generalizable data-driven multicellular model of pancreatic ductal adenocarcinoma

BACKGROUND

Mechanistic models, when combined with pertinent data, can improve our knowledge regarding important molecular and cellular mechanisms found in cancer. These models make the prediction of tissue-level response to drug treatment possible, which can lead to new therapies and improved patient outcomes. Here we present a data-driven multiscale modeling framework to study molecular interactions between cancer, stromal, and immune cells found in the tumor microenvironment. We also develop methods to use molecular data available in The Cancer Genome Atlas to generate sample-specific models of cancer.

RESULTS

By combining published models of different cells relevant to pancreatic ductal adenocarcinoma (PDAC), we built an agent-based model of the multicellular pancreatic tumor microenvironment, formally describing cell type-specific molecular interactions and cytokine-mediated cell-cell communications. We used an ensemble-based modeling approach to systematically explore how variations in the tumor microenvironment affect the viability of cancer cells. The results suggest that the autocrine loop involving EGF signaling is a key interaction modulator between pancreatic cancer and stellate cells. EGF is also found to be associated with previously described subtypes of PDAC. Moreover, the model allows a systematic exploration of the effect of possible therapeutic perturbations; our simulations suggest that reducing bFGF secretion by stellate cells will have, on average, a positive impact on cancer apoptosis.

CONCLUSIONS

The developed framework allows model-driven hypotheses to be generated regarding therapeutically relevant PDAC states with potential molecular and cellular drivers indicating specific intervention strategies.

The greater inflammatory pathway-high clinical potential by innovative predictive, preventive, and personalized medical approach

BACKGROUND AND LIMITATIONS

Impaired wound healing (WH) and chronic inflammation are hallmarks of non-communicable diseases (NCDs). However, despite WH being a recognized player in NCDs, mainstream therapies focus on (un)targeted damping of the inflammatory response, leaving WH largely unaddressed, owing to three main factors. The first is the complexity of the pathway that links inflammation and wound healing; the second is the dual nature, local and systemic, of WH; and the third is the limited acknowledgement of genetic and contingent causes that disrupt physiologic progression of WH.

PROPOSED APPROACH

Here, in the frame of Predictive, Preventive, and Personalized Medicine (PPPM), we integrate and revisit current literature to offer a novel systemic view on the cues that can impact on the fate (acute or chronic inflammation) of WH, beyond the compartmentalization of medical disciplines and with the support of advanced computational biology.

CONCLUSIONS

This shall open to a broader understanding of the causes for WH going awry, offering new operational criteria for patients' stratification (prediction and personalization). While this may also offer improved options for targeted prevention, we will envisage new therapeutic strategies to reboot and/or boost WH, to enable its progression across its physiological phases, the first of which is a transient acute inflammatory response versus the chronic low-grade inflammation characteristic of NCDs.

A mechanistic integrative computational model of macrophage polarization: Implications in human pathophysiology

Macrophages respond to signals in the microenvironment by changing their functional phenotypes, a process known as polarization. Depending on the context, they acquire different patterns of transcriptional activation, cytokine expression and cellular metabolism which collectively constitute a continuous spectrum of phenotypes, of which the two extremes are denoted as classical (M1) and alternative (M2) activation. To quantitatively decode the underlying principles governing macrophage phenotypic polarization and thereby harness its therapeutic potential in human diseases, a systems-level approach is needed given the multitude of signaling pathways and intracellular regulation involved. Here we develop the first mechanism-based, multi-pathway computational model that describes the integrated signal transduction and macrophage programming under M1 (IFN-γ), M2 (IL-4) and cell stress (hypoxia) stimulation. Our model was calibrated extensively against experimental data, and we mechanistically elucidated several signature feedbacks behind the M1-M2 antagonism and investigated the dynamical shaping of macrophage phenotypes within the M1-M2 spectrum. Model sensitivity analysis also revealed key molecular nodes and interactions as targets with potential therapeutic values for the pathophysiology of peripheral arterial disease and cancer. Through simulations that dynamically capture the signal integration and phenotypic marker expression in the differential macrophage polarization responses, our model provides an important computational basis toward a more quantitative and network-centric understanding of the complex physiology and versatile functions of macrophages in human diseases.

Modeling the Effect of High Calorie Diet on the Interplay between Adipose Tissue, Inflammation, and Diabetes

Background

Type 2 diabetes (T2D) is a chronic metabolic disease potentially leading to serious widespread tissue damage. Human organism develops T2D when the glucose-insulin control is broken for reasons that are not fully understood but have been demonstrated to be linked to the emergence of a chronic inflammation. Indeed such low-level chronic inflammation affects the pancreatic production of insulin and triggers the development of insulin resistance, eventually leading to an impaired control of the blood glucose concentration. On the contrary, it is well-known that obesity and inflammation are strongly correlated.

Aim

In this study, we investigate in silico the effect of overfeeding on the adipose tissue and the consequent set up of an inflammatory state. We model the emergence of the inflammation as the result of adipose mass increase which, in turn, is a direct consequence of a prolonged excess of high calorie intake.

Results

The model reproduces the fat accumulation due to excessive caloric intake observed in two clinical studies. Moreover, while showing consistent weight gains over long periods of time, it reveals a drift of the macrophage population toward the proinflammatory phenotype, thus confirming its association with fatness.

Gene Regulatory Network Modeling of Macrophage Differentiation Corroborates the Continuum Hypothesis of Polarization States

Macrophages derived from monocyte precursors undergo specific polarization processes which are influenced by the local tissue environment: classically activated (M1) macrophages, with a pro-inflammatory activity and a role of effector cells in Th1 cellular immune responses, and alternatively activated (M2) macrophages, with anti-inflammatory functions and involved in immunosuppression and tissue repair. At least three different subsets of M2 macrophages, namely, M2a, M2b, and M2c, are characterized in the literature based on their eliciting signals. The activation and polarization of macrophages is achieved through many, often intertwined, signaling pathways. To describe the logical relationships among the genes involved in macrophage polarization, we used a computational modeling methodology, namely, logical (Boolean) modeling of gene regulation. We integrated experimental data and knowledge available in the literature to construct a logical network model for the gene regulation driving macrophage polarization to the M1, M2a, M2b, and M2c phenotypes. Using the software GINsim and BoolNet, we analyzed the network dynamics under different conditions and perturbations to understand how they affect cell polarization. Dynamic simulations of the network model, enacting the most relevant biological conditions, showed coherence with the observed behavior of in vivo macrophages. The model could correctly reproduce the polarization toward the four main phenotypes as well as to several hybrid phenotypes, which are known to be experimentally associated to physiological and pathological conditions. We surmise that shifts among different phenotypes in the model mimic the hypothetical continuum of macrophage polarization, with M1 and M2 being the extremes of an uninterrupted sequence of states. Furthermore, model simulations suggest that anti-inflammatory macrophages are resilient to shift back to the pro-inflammatory phenotype.

Recent advances in the treatment of skin involvement in systemic sclerosis

Skin fibrosis is a devastating clinical condition commonly seen in skin-restricted and systemic disorders. The goal of skin fibrosis treatment is the restoration of abnormally activated dermal fibroblasts producing the excessive amount of extracellular matrix, which is generally a final consequence of the complex disease process including the activation of vascular and immune systems. Among various skin fibrotic conditions, the molecular mechanisms underlying dermal fibroblast activation have been mostly well studied in systemic sclerosis (SSc). SSc is a multisystem autoimmune and vascular disease resulting in extensive fibrosis of the skin and various internal organs. Since SSc pathogenesis is believed to include all the critical components regulating tissue fibrosis, the studies on anti-fibrotic drugs against SSc provide us much useful information regarding the strategy for the treatment of various skin fibrotic conditions. In the recent decade, as is the case with other autoimmune and inflammatory diseases, the molecular targeting therapy with monoclonal antibody has been clinically well examined in SSc. Promising clinical outcomes are so far reported in tocilizumab (an anti-IL-6 receptor antibody), rituximab (an anti-CD20 antibody), and fresolimumab (an anti-TGF-β antibody). The analysis of gene expression profiles in skin lesions of SSc patients treated with tocilizumab or fresolimumab revealed a critical role of monocyte-macrophage lineage cells in the development of skin fibrosis and the involvement of IL-6 and TGF-β in the activation of those cells. Considering that B cells modulate the differentiation and activation of macrophages, favorable clinical outcomes of rituximab treatment imply the central role of B cell/monocyte-macrophage lineage cell axis in the pathogenesis of SSc. This scenario may be applicable at least partly to other skin fibrotic conditions. In this review article, the currently available data on these drugs are summarized and the future directions are discussed.

Bioinformatics and systems biology research update from the 15th International Conference on Bioinformatics (InCoB2016)

The International Conference on Bioinformatics (InCoB) has been publishing peer-reviewed conference papers in BMC Bioinformatics since 2006. Of the 44 articles accepted for publication in supplement issues of BMC Bioinformatics, BMC Genomics, BMC Medical Genomics and BMC Systems Biology, 24 articles with a bioinformatics or systems biology focus are reviewed in this editorial. InCoB2017 is scheduled to be held in Shenzen, China, September 20-22, 2017.