BioNSi: A Discrete Biological Network Simulator Tool

Application of the PHENotype SIMulator for rapid identification of potential candidates in effective COVID-19 drug repurposing

The current, rapidly diversifying pandemic has accelerated the need for efficient and effective identification of potential drug candidates for COVID-19. Knowledge on host-immune response to SARS-CoV-2 infection, however, remains limited with few drugs approved to date. Viable strategies and tools are rapidly arising to address this, especially with repurposing of existing drugs offering significant promise. Here we introduce a systems biology tool, the PHENotype SIMulator, which -by leveraging available transcriptomic and proteomic databases-allows modeling of SARS-CoV-2 infection in host cells in silico to i) determine with high sensitivity and specificity (both>96%) the viral effects on cellular host-immune response, resulting in specific cellular SARS-CoV-2 signatures and ii) utilize these cell-specific signatures to identify promising repurposable therapeutics. Powered by this tool, coupled with domain expertise, we identify several potential COVID-19 drugs including methylprednisolone and metformin, and further discern key cellular SARS-CoV-2-affected pathways as potential druggable targets in COVID-19 pathogenesis.

In silico model of atherosclerosis with individual patient calibration to enable precision medicine for cardiovascular disease

OBJECTIVE

Guidance for preventing myocardial infarction and ischemic stroke by tailoring treatment for individual patients with atherosclerosis is an unmet need. Such development may be possible with computational modeling. Given the multifactorial biology of atherosclerosis, modeling must be based on complete biological networks that capture protein-protein interactions estimated to drive disease progression. Here, we aimed to develop a clinically relevant scale model of atherosclerosis, calibrate it with individual patient data, and use it to simulate optimized pharmacotherapy for individual patients.

APPROACH AND RESULTS

The study used a uniquely constituted plaque proteomic dataset to create a comprehensive systems biology disease model for simulating individualized responses to pharmacotherapy. Plaque tissue was collected from 18 patients with 6735 proteins at two locations per patient. 113 pathways were identified and included in the systems biology model of endothelial cells, vascular smooth muscle cells, macrophages, lymphocytes, and the integrated intima, altogether spanning 4411 proteins, demonstrating a range of 39-96% plaque instability. After calibrating the systems biology models for individual patients, we simulated intensive lipid-lowering, anti-inflammatory, and anti-diabetic drugs. We also simulated a combination therapy. Drug response was evaluated as the degree of change in plaque stability, where an improvement was defined as a reduction of plaque instability. In patients with initially unstable lesions, simulated responses varied from high (20%, on combination therapy) to marginal improvement, whereas patients with initially stable plaques showed generally less improvement.

CONCLUSION

In this pilot study, proteomics-based system biology modeling was shown to simulate drug response based on atherosclerotic plaque instability with a power of 90%, providing a potential strategy for improved personalized management of patients with cardiovascular disease.

Systems neuroimmunology: a review of multiomics methodologies to characterize neuroimmunological interactions in spinal and cranial diseases

Neuroimmunology plays a critical role in our understanding of the pathophysiological processes that underlie a variety of diseases treated by neurosurgeons, including degenerative disc disease (DDD), glioblastoma (GBM), aneurysmal subarachnoid hemorrhage (aSAH), and others. Compared with traditional methods in neuroimmunology, which study one pathway or gene at a time, emerging multiomics methodologies allow for holistic interrogation of multiple immune-signaling pathways to test hypotheses and the effects of therapeutics at a systems level. In this review, the authors summarize key concepts for gathering and analyzing multiomics data so that neurosurgeons can contribute to the emerging field of systems neuroimmunology. Additionally, they describe 3 use cases, based on original research published by their group and others, that utilize transcriptomic, metabolomic, and proteomic analyses to study immune-signaling pathways in DDD, aSAH, and GBM. Through these use cases, techniques for performing machine learning and network-based analyses to generate new clinical insights from multiomics data are shared. The authors hope that neurosurgeons might use this review as a summary of common tools and principles in systems immunology to better engage in creating the immunotherapies of tomorrow.

PHENSIM: Phenotype Simulator

Despite the unprecedented growth in our understanding of cell biology, it still remains challenging to connect it to experimental data obtained with cells and tissues’ physiopathological status under precise circumstances. This knowledge gap often results in difficulties in designing validation experiments, which are usually labor-intensive, expensive to perform, and hard to interpret. Here we propose PHENSIM, a computational tool using a systems biology approach to simulate how cell phenotypes are affected by the activation/inhibition of one or multiple biomolecules, and it does so by exploiting signaling pathways. Our tool’s applications include predicting the outcome of drug administration, knockdown experiments, gene transduction, and exposure to exosomal cargo. Importantly, PHENSIM enables the user to make inferences on well-defined cell lines and includes pathway maps from three different model organisms. To assess our approach’s reliability, we built a benchmark from transcriptomics data gathered from NCBI GEO and performed four case studies on known biological experiments. Our results show high prediction accuracy, thus highlighting the capabilities of this methodology. PHENSIM standalone Java application is available at https://github.com/alaimos/phensim, along with all data and source codes for benchmarking. A web-based user interface is accessible at https://phensim.tech/.

Despite the unprecedented growth in our understanding of cell biology, it still remains challenging to connect it to experimental data obtained with cells and tissues’ physiopathological status under precise circumstances. This knowledge gap often results in difficulties in designing validation experiments, which are usually labor-intensive, expensive to perform, and hard to interpret. In this context, ’in silico’ simulations can be extensively applied in massive scales, testing thousands of hypotheses under various conditions, which is usually experimentally infeasible. At present, many simulation models have become available. However, complex biological networks might pose challenges to their performance.

We propose PHENSIM, a computational tool using a systems biology approach to simulate how cell phenotypes are affected by the activation/inhibition of one or multiple biomolecules, and it does so by exploiting signaling pathways. We implemented our tool as a freely accessible web application, hoping to allow ’in silico’ simulations to play a more central role in the modeling and understanding of biological phenomena.

Rapid Identification of Druggable Targets and the Power of the PHENotype SIMulator for Effective Drug Repurposing in COVID-19

The current, rapidly diversifying pandemic has accelerated the need for efficient and effective identification of potential drug candidates for COVID-19. Knowledge on host-immune response to SARS-CoV-2 infection, however, remains limited with very few drugs approved to date. Viable strategies and tools are rapidly arising to address this, especially with repurposing of existing drugs offering significant promise. Here we introduce a systems biology tool, the PHENotype SIMulator, which - by leveraging available transcriptomic and proteomic databases - allows modeling of SARS-CoV-2 infection in host cells in silico to i) determine with high sensitivity and specificity (both > 96%) the viral effects on cellular host-immune response, resulting in a specific cellular SARS-CoV-2 signature and ii) utilize this specific signature to narrow down promising repurposable therapeutic strategies. Powered by this tool, coupled with domain expertise, we have identified several potential COVID-19 drugs including methylprednisolone and metformin, and further discern key cellular SARS-CoV-2-affected pathways as potential new druggable targets in COVID-19 pathogenesis.

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.

Systems Analysis of Insulin and IGF1 Receptors Networks in Breast Cancer Cells Identifies Commonalities and Divergences in Expression Patterns

Insulin and insulin-like growth factor-1 (IGF1), acting respectively via the insulin (INSR) and IGF1 (IGF1R) receptors, play key developmental and metabolic roles throughout life. In addition, both signaling pathways fulfill important roles in cancer initiation and progression. The present study was aimed at identifying mechanistic differences between INSR and IGF1R using a recently developed bioinformatics tool, the Biological Network Simulator (BioNSi). This application allows to import and merge multiple pathways and interaction information from the KEGG database into a single network representation. The BioNsi network simulation tool allowed us to exploit the availability of gene expression data derived from breast cancer cell lines with specific disruptions of the INSR or IGF1R genes in order to investigate potential differences in protein expression that might be linked to biological attributes of the specific receptor networks. Modeling-generated information was corroborated by experimental and biological assays. BioNSi analyses revealed that the expression of 75 and 71 genes changed during simulation of IGF1R-KD and INSR-KD, compared to control cells, respectively. Out of 16 proteins that BioNSi analysis was based on, validated by Western blotting, nine were shown to be involved in DNA repair, eight in cell cycle checkpoints, six in proliferation, eight in apoptosis, seven in oxidative stress, six in cell migration, two in energy homeostasis, and three in senescence. Taken together, analyses identified a number of commonalities and, most importantly, dissimilarities between the IGF1R and INSR pathways that might help explain the basis for the biological differences between these networks.

A Note on GRegNetSim: A Tool for the Discrete Simulation and Analysis of Genetic Regulatory Networks

Discrete simulations of genetic regulatory networks were used to study subsystems of yeast successfully. However, implementations of existing models underlying these simulations do not support a graphic interface, and require computations necessary to analyze their results to be done manually. Furthermore, differences between existing models suggest that an enriched model, encompassing both existing models, is needed. We developed a software tool, GRegNetSim, that allows the end-user to describe genetic regulatory networks graphically. The user can specify various transition functions at different nodes of the network, supporting, for example, threshold and gradient effects, and then apply the network to a variety of inputs. GRegNetSim displays the relationship between the inputs and the mode of behavior of the network in a graphic form that is easy to interpret. Furthermore, it can automatically extract statistical data necessary to analyze the simulations. The discrete simulations performed by GRegNetSim can be used to elucidate and predict the behavior, structure and properties of genetic regulatory networks in a unified manner. GRegNetSim is implemented as a Cytoscape App. Installation files, examples and source code, along with a detailed user guide, are freely available at https://sites.google.com/site/gregnetsim/.

Interactions between the circadian clock and TGF-β signaling pathway in zebrafish

Background

TGF-β signaling is a cellular pathway that functions in most cells and has been shown to play a role in multiple processes, such as the immune response, cell differentiation and proliferation. Recent evidence suggests a possible interaction between TGF-β signaling and the molecular circadian oscillator. The current study aims to characterize this interaction in the zebrafish at the molecular and behavioral levels, taking advantage of the early development of a functional circadian clock and the availability of light-entrainable clock-containing cell lines.

Results

Smad3a, a TGF-β signaling-related gene, exhibited a circadian expression pattern throughout the brain of zebrafish larvae. Both pharmacological inhibition and indirect activation of TGF-β signaling in zebrafish Pac-2 cells caused a concentration dependent disruption of rhythmic promoter activity of the core clock gene Per1b. Inhibition of TGF-β signaling in intact zebrafish larvae caused a phase delay in the rhythmic expression of Per1b mRNA. TGF-β inhibition also reversibly disrupted, phase delayed and increased the period of circadian rhythms of locomotor activity in zebrafish larvae.

Conclusions

The current research provides evidence for an interaction between the TGF-β signaling pathway and the circadian clock system at the molecular and behavioral levels, and points to the importance of TGF-β signaling for normal circadian clock function. Future examination of this interaction should contribute to a better understanding of its underlying mechanisms and its influence on a variety of cellular processes including the cell cycle, with possible implications for cancer development and progression.

Simulation and visualization of multiple KEGG pathways using BioNSi

Motivation: Many biologists are discouraged from using network simulation tools because these require manual, often tedious network construction. This situation calls for building new tools or extending existing ones with the ability to import biological pathways previously deposited in databases and analyze them, in order to produce novel biological insights at the pathway level.

Results: We have extended a network simulation tool (BioNSi), which now allows merging of multiple pathways from the KEGG pathway database into a single, coherent network, and visualizing its properties. Furthermore, the enhanced tool enables loading experimental expression data into the network and simulating its dynamics under various biological conditions or perturbations. As a proof of concept, we tested two sets of published experimental data, one related to inflammatory bowel disease condition and the other to breast cancer treatment. We predict some of the major observations obtained following these laboratory experiments, and provide new insights that may shed additional light on these results.

Tool requirements: Cytoscape 3.x, JAVA 8

Availability: The tool is freely available at http://bionsi.wix.com/bionsi, where a complete user guide and a step-by-step manual can also be found.