Boolean modeling reveals that cyclic attractors in macrophage polarization serve as reservoirs of states to balance external perturbations from the tumor microenvironment.

Cyclic attractors generated from Boolean models may explain the adaptability of a cell in response to a dynamical complex tumor microenvironment. In contrast to this idea, we postulate that cyclic attractors in certain cases could be a systemic mechanism to face the perturbations coming from the environment. To justify our conjecture, we present a dynamic analysis of a highly curated transcriptional regulatory network of macrophages constrained into a cancer microenvironment. We observed that when M1-associated transcription factors (STAT1 or NF-κB) are perturbed and the microenvironment balances to a hyper-inflammation condition, cycle attractors activate genes whose signals counteract this effect implicated in tissue damage. The same behavior happens when the M2-associated transcription factors are disturbed (STAT3 or STAT6); cycle attractors will prevent a hyper-regulation scenario implicated in providing a suitable environment for tumor growth. Therefore, here we propose that cyclic macrophage phenotypes can serve as a reservoir for balancing the phenotypes when a specific phenotype-based transcription factor is perturbed in the regulatory network of macrophages. We consider that cyclic attractors should not be simply ignored, but it is necessary to carefully evaluate their biological importance. In this work, we suggest one conjecture: the cyclic attractors can serve as a reservoir to balance the inflammatory/regulatory response of the network under external perturbations.

Transcriptional and Microenvironmental Landscape of Macrophage Transition in Cancer: A Boolean Analysis.

The balance between pro- and anti-inflammatory immune system responses is crucial to face and counteract complex diseases such as cancer. Macrophages are an essential population that contributes to this balance in collusion with the local tumor microenvironment. Cancer cells evade the attack of macrophages by liberating cytokines and enhancing the transition to the M2 phenotype with pro-tumoral functions. Despite this pernicious effect on immune systems, the M1 phenotype still exists in the environment and can eliminate tumor cells by liberating cytokines that recruit and activate the cytotoxic actions of TH1 effector cells. Here, we used a Boolean modeling approach to understand how the tumor microenvironment shapes macrophage behavior to enhance pro-tumoral functions. Our network reconstruction integrates experimental data and public information that let us study the polarization from monocytes to M1, M2a, M2b, M2c, and M2d subphenotypes. To analyze the dynamics of our model, we modeled macrophage polarization in different conditions and perturbations. Notably, our study identified new hybrid cell populations, undescribed before. Based on the in vivo macrophage behavior, we explained the hybrid macrophages' role in the tumor microenvironment. The in silico model allowed us to postulate transcriptional factors that maintain the balance between macrophages with anti- and pro-tumoral functions. In our pursuit to maintain the balance of macrophage phenotypes to eliminate malignant tumor cells, we emulated a theoretical genetically modified macrophage by modifying the activation of NFκB and a loss of function in HIF1-α and discussed their phenotype implications. Overall, our theoretical approach is as a guide to design new experiments for unraveling the principles of the dual host-protective or -harmful antagonistic roles of transitional macrophages in tumor immunoediting and cancer cell fate decisions.

Some perspectives on network modeling in therapeutic target prediction.

Drug target identification is of significant commercial interest to pharmaceutical companies, and there is a vast amount of research done related to the topic of therapeutic target identification. Interdisciplinary research in this area involves both the biological network community and the graph algorithms community. Key steps of a typical therapeutic target identification problem include synthesizing or inferring the complex network of interactions relevant to the disease, connecting this network to the disease-specific behavior, and predicting which components are key mediators of the behavior. All of these steps involve graph theoretical or graph algorithmic aspects. In this perspective, we provide modelling and algorithmic perspectives for therapeutic target identification and highlight a number of algorithmic advances, which have gotten relatively little attention so far, with the hope of strengthening the ties between these two research communities.