Network perturbation analysis of gene transcriptional profiles reveals protein targets and mechanism of action of drugs and influenza A viral infection

Decoding HiPSC-CM's Response to SARS-CoV-2: mapping the molecular landscape of cardiac injury

BACKGROUND

Acute cardiac injury caused by coronavirus disease 2019 (COVID-19) increases mortality. Acute cardiac injury caused by COVID-19 requires understanding how severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) directly infects cardiomyocytes. This study provides a solid foundation for related studies by using a model of SARS-CoV-2 infection in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) at the transcriptome level, highlighting the relevance of this study to related studies. SARS-CoV-2 infection in hiPSC-CMs has previously been studied by bioinformatics without presenting the full molecular biological process. We present a unique bioinformatics view of the complete molecular biological process of SARS-CoV-2 infection in hiPSC-CMs.

METHODS

To validate the RNA-seq datasets, we used GSE184715 and GSE150392 for the analytical studies, GSE193722 for validation at the cellular level, and GSE169241 for validation in heart tissue samples. GeneCards and MsigDB databases were used to find genes associated with the phenotype. In addition to differential expression analysis and principal component analysis (PCA), we also performed protein-protein interaction (PPI) analysis, functional enrichment analysis, hub gene analysis, upstream transcription factor prediction, and drug prediction.

RESULTS

Differentially expressed genes (DEGs) were classified into four categories: cardiomyocyte cytoskeletal protein inhibition, proto-oncogene activation and inflammation, mitochondrial dysfunction, and intracellular cytoplasmic physiological function. Each of the hub genes showed good diagnostic prediction, which was well validated in other datasets. Inhibited biological functions included cardiomyocyte cytoskeletal proteins, adenosine triphosphate (ATP) synthesis and electron transport chain (ETC), glucose metabolism, amino acid metabolism, fatty acid metabolism, pyruvate metabolism, citric acid cycle, nucleic acid metabolism, replication, transcription, translation, ubiquitination, autophagy, and cellular transport. Proto-oncogenes, inflammation, nuclear factor-kappaB (NF-κB) pathways, and interferon signaling were activated, as well as inflammatory factors. Viral infection activates multiple pathways, including the interferon pathway, proto-oncogenes and mitochondrial oxidative stress, while inhibiting cardiomyocyte backbone proteins and energy metabolism. Infection limits intracellular synthesis and metabolism, as well as the raw materials for mitochondrial energy synthesis. Mitochondrial dysfunction and energy abnormalities are ultimately caused by proto-oncogene activation and SARS-CoV-2 infection. Activation of the interferon pathway, proto-oncogene up-regulation, and mitochondrial oxidative stress cause the inflammatory response and lead to diminished cardiomyocyte contraction. Replication, transcription, translation, ubiquitination, autophagy, and cellular transport are among the functions that decline physiologically.

CONCLUSION

SARS-CoV-2 infection in hiPSC-CMs is fundamentally mediated via mitochondrial dysfunction. Therapeutic interventions targeting mitochondrial dysfunction may alleviate the cardiovascular complications associated with SARS-CoV-2 infection.

Drug target prediction through deep learning functional representation of gene signatures

Many machine learning applications in bioinformatics currently rely on matching gene identities when analyzing input gene signatures and fail to take advantage of preexisting knowledge about gene functions. To further enable comparative analysis of OMICS datasets, including target deconvolution and mechanism of action studies, we develop an approach that represents gene signatures projected onto their biological functions, instead of their identities, similar to how the word2vec technique works in natural language processing. We develop the Functional Representation of Gene Signatures (FRoGS) approach by training a deep learning model and demonstrate that its application to the Broad Institute's L1000 datasets results in more effective compound-target predictions than models based on gene identities alone. By integrating additional pharmacological activity data sources, FRoGS significantly increases the number of high-quality compound-target predictions relative to existing approaches, many of which are supported by in silico and/or experimental evidence. These results underscore the general utility of FRoGS in machine learning-based bioinformatics applications. Prediction networks pre-equipped with the knowledge of gene functions may help uncover new relationships among gene signatures acquired by large-scale OMICs studies on compounds, cell types, disease models, and patient cohorts.

An overview of recent advances and challenges in predicting compound-protein interaction (CPI)

Compound-protein interactions (CPIs) are critical in drug discovery for identifying therapeutic targets, drug side effects, and repurposing existing drugs. Machine learning (ML) algorithms have emerged as powerful tools for CPI prediction, offering notable advantages in cost-effectiveness and efficiency. This review provides an overview of recent advances in both structure-based and non-structure-based CPI prediction ML models, highlighting their performance and achievements. It also offers insights into CPI prediction-related datasets and evaluation benchmarks. Lastly, the article presents a comprehensive assessment of the current landscape of CPI prediction, elucidating the challenges faced and outlining emerging trends to advance the field.

Gene network inference from single-cell omics data and domain knowledge for constructing COVID-19-specific ICAM1-associated pathways

Introduction: Intercellular adhesion molecule 1 (ICAM-1) is a critical molecule responsible for interactions between cells. Previous studies have suggested that ICAM-1 triggers cell-to-cell transmission of HIV-1 or HTLV-1, that SARS-CoV-2 shares several features with these viruses via interactions between cells, and that SARS-CoV-2 cell-to-cell transmission is associated with COVID-19 severity. From these previous arguments, it is assumed that ICAM-1 can be related to SARS-CoV-2 cell-to-cell transmission in COVID-19 patients. Indeed, the time-dependent change of the ICAM-1 expression level has been detected in COVID-19 patients. However, signaling pathways that consist of ICAM-1 and other molecules interacting with ICAM-1 are not identified in COVID-19. For example, the current COVID-19 Disease Map has no entry for those pathways. Therefore, discovering unknown ICAM1-associated pathways will be indispensable for clarifying the mechanism of COVID-19. Materials and methods: This study builds ICAM1-associated pathways by gene network inference from single-cell omics data and multiple knowledge bases. First, single-cell omics data analysis extracts coexpressed genes with significant differences in expression levels with spurious correlations removed. Second, knowledge bases validate the models. Finally, mapping the models onto existing pathways identifies new ICAM1-associated pathways. Results: Comparison of the obtained pathways between different cell types and time points reproduces the known pathways and indicates the following two unknown pathways: (1) upstream pathway that includes proteins in the non-canonical NF-κB pathway and (2) downstream pathway that contains integrins and cytoskeleton or motor proteins for cell transformation. Discussion: In this way, data-driven and knowledge-based approaches are integrated into gene network inference for ICAM1-associated pathway construction. The results can contribute to repairing and completing the COVID-19 Disease Map, thereby improving our understanding of the mechanism of COVID-19.

Benchmarking causal reasoning algorithms for gene expression-based compound mechanism of action analysis

BACKGROUND

Elucidating compound mechanism of action (MoA) is beneficial to drug discovery, but in practice often represents a significant challenge. Causal Reasoning approaches aim to address this situation by inferring dysregulated signalling proteins using transcriptomics data and biological networks; however, a comprehensive benchmarking of such approaches has not yet been reported. Here we benchmarked four causal reasoning algorithms (SigNet, CausalR, CausalR ScanR and CARNIVAL) with four networks (the smaller Omnipath network vs. 3 larger MetaBase™ networks), using LINCS L1000 and CMap microarray data, and assessed to what extent each factor dictated the successful recovery of direct targets and compound-associated signalling pathways in a benchmark dataset comprising 269 compounds. We additionally examined impact on performance in terms of the functions and roles of protein targets and their connectivity bias in the prior knowledge networks.

RESULTS

According to statistical analysis (negative binomial model), the combination of algorithm and network most significantly dictated the performance of causal reasoning algorithms, with the SigNet recovering the greatest number of direct targets. With respect to the recovery of signalling pathways, CARNIVAL with the Omnipath network was able to recover the most informative pathways containing compound targets, based on the Reactome pathway hierarchy. Additionally, CARNIVAL, SigNet and CausalR ScanR all outperformed baseline gene expression pathway enrichment results. We found no significant difference in performance between L1000 data or microarray data, even when limited to just 978 'landmark' genes. Notably, all causal reasoning algorithms also outperformed pathway recovery based on input DEGs, despite these often being used for pathway enrichment. Causal reasoning methods performance was somewhat correlated with connectivity and biological role of the targets.

CONCLUSIONS

Overall, we conclude that causal reasoning performs well at recovering signalling proteins related to compound MoA upstream from gene expression changes by leveraging prior knowledge networks, and that the choice of network and algorithm has a profound impact on the performance of causal reasoning algorithms. Based on the analyses presented here this is true for both microarray-based gene expression data as well as those based on the L1000 platform.

Reverse tracking from drug-induced transcriptomes through multilayer molecular networks reveals hidden drug targets

Identifying molecular targets of a drug is an essential process for drug discovery and development. The recent in-silico approaches are usually based on the structure information of chemicals and proteins. However, 3D structure information is hard to obtain and machine-learning methods using 2D structure suffer from data imbalance problem. Here, we present a reverse tracking method from genes to target proteins using drug-perturbed gene transcriptional profiles and multilayer molecular networks. We scored how well the protein explains gene expression changes perturbed by a drug. We validated the protein scores of our method in predicting known targets of drugs. Our method performs better than other methods using the gene transcriptional profiles and shows the ability to suggest the molecular mechanism of drugs. Furthermore, our method has the potential to predict targets for objects that do not have rigid structural information, such as coronavirus.

Roles of Skeletal Muscle in Development: A Bioinformatics and Systems Biology Overview

The ability to assess various cellular events consequent to perturbations, such as genetic mutations, disease states and therapies, has been recently revolutionized by technological advances in multiple "omics" fields. The resulting deluge of information has enabled and necessitated the development of tools required to both process and interpret the data. While of tremendous value to basic researchers, the amount and complexity of the data has made it extremely difficult to manually draw inference and identify factors key to the study objectives. The challenges of data reduction and interpretation are being met by the development of increasingly complex tools that integrate disparate knowledge bases and synthesize coherent models based on current biological understanding. This chapter presents an example of how genomics data can be integrated with biological network analyses to gain further insight into the developmental consequences of genetic perturbations. State of the art methods for conducting similar studies are discussed along with modern methods used to analyze and interpret the data.

Fuzzy and Rough Set Theory Based Computational Framework for Mining Genetic Interaction Triplets From Gene Expression Profiles for Lung Adenocarcinoma

Genetic interactions are very helpful in understanding different disease and discovering drugs for it. Compared to the gene pairs that represent the genetic interactions between two genes, the gene triplets are more informative and useful. However, existing works on genetic interactions among gene triplets have primarily focused on detecting gene triplets from time series gene expression profiles. Generating the time series gene expression profiles for humans is quite impracticable but the labeled gene expression profiles are available for different diseases in case of humans. In this paper, a computational framework has been proposed to detect gene triplets from labeled gene expression profiles. First, it employs Rough Set Theory for extracting the key genes and then designs a fuzzy inference system for generating possible gene triplets. Further, Root Mean Squared Error measure has been used to prune out the irrelevant gene triplets. In the present work, the proposed computational framework has been applied to labeled lung adenocarcinoma dataset and can be applied to any other labeled gene expression dataset. The extracted gene triplets and their functionalities have been verified with existing biological literature and benchmark databases and the results of verification signify that the proposed framework is promising in terms of finding useful genetic triplets. Further, the proposed framework has been found more efficient as compared to an existing mutual information-based technique in terms of detecting known genetic interactions.

FL-DTD: an integrated pipeline to predict the drug interacting targets by feedback loop-based network analysis

Drug target discovery is an essential step to reveal the mechanism of action (MoA) underlying drug therapeutic effects and/or side effects. Most of the approaches are usually labor-intensive while unable to identify the tissue-specific interacting targets, especially the targets with weaker drug binding affinity. In this work, we proposed an integrated pipeline, FL-DTD, to predict the drug interacting targets of novel compounds in a tissue-specific manner. This method was built based on a hypothesis that cells under a status of homeostasis would take responses to drug perturbation by activating feedback loops. Therefore, the drug interacting targets can be predicted by analyzing the network responses after drug perturbation. We evaluated this method using the expression data of estrogen stimulation, gene manipulation and drug perturbation and validated its good performance to identify the annotated drug targets. Using STAT3 as a target protein, we applied this method to drug perturbation data of 500 natural compounds and predicted five compounds with STAT3 interacting activities. Experimental assay validated the STAT3-interacting activities of four compounds. Overall, our evaluation suggests that FL-DTD predicts the drug interacting targets with good accuracy and can be used for drug target discovery.

Computational Identification of Potential Multitarget Inhibitors of Nipah Virus by Molecular Docking and Molecular Dynamics

Nipah virus (NiV) is a recently emerged paramyxovirus that causes severe encephalitis and respiratory diseases in humans. Despite the severe pathogenicity of this virus and its pandemic potential, not even a single type of molecular therapeutics has been approved for human use. Considering the role of NiV attachment glycoprotein G (NiV-G), fusion glycoprotein (NiV-F), and nucleoprotein (NiV-N) in virus replication and spread, these are the most attractive targets for anti-NiV drug discovery. Therefore, to prospect for potential multitarget chemical/phytochemical inhibitor(s) against NiV, a sequential molecular docking and molecular-dynamics-based approach was implemented by simultaneously targeting NiV-G, NiV-F, and NiV-N. Information on potential NiV inhibitors was compiled from the literature, and their 3D structures were drawn manually, while the information and 3D structures of phytochemicals were retrieved from the established structural databases. Molecules were docked against NiV-G (PDB ID:2VSM), NiV-F (PDB ID:5EVM), and NiV-N (PDB ID:4CO6) and then prioritized based on (1) strong protein-binding affinity, (2) interactions with critically important binding-site residues, (3) ADME and pharmacokinetic properties, and (4) structural stability within the binding site. The molecules that bind to all the three viral proteins (NiV-G ∩ NiV-F ∩ NiV-N) were considered multitarget inhibitors. This study identified phytochemical molecules RASE0125 (17-O-Acetyl-nortetraphyllicine) and CARS0358 (NA) as distinct multitarget inhibitors of all three viral proteins, and chemical molecule ND_nw_193 (RSV604) as an inhibitor of NiV-G and NiV-N. We expect the identified compounds to be potential candidates for in vitro and in vivo antiviral studies, followed by clinical treatment of NiV.

Expanding the search for small-molecule antibacterials by multidimensional profiling

New techniques for systematic profiling of small-molecule effects can enhance traditional growth inhibition screens for antibiotic discovery and change how we search for new antibacterial agents. Computational models that integrate physicochemical compound properties with their phenotypic and molecular downstream effects can not only predict efficacy of molecules yet to be tested, but also reveal unprecedented insights on compound modes of action (MoAs). The unbiased characterization of compounds that themselves are not growth inhibitory but exhibit diverse MoAs, can expand antibacterial strategies beyond direct inhibition of core essential functions. Early and systematic functional annotation of compound libraries thus paves the way to new models in the selection of lead antimicrobial compounds. In this Review, we discuss how multidimensional small-molecule profiling and the ever-increasing computing power are accelerating the discovery of unconventional antibacterials capable of bypassing resistance and exploiting synergies with established antibacterial treatments and with protective host mechanisms.

Computational analyses of mechanism of action (MoA): data, methods and integration

The elucidation of a compound's Mechanism of Action (MoA) is a challenging task in the drug discovery process, but it is important in order to rationalise phenotypic findings and to anticipate potential side-effects. Bioinformatic approaches, advances in machine learning techniques and the increasing deposition of high-throughput data in public databases have significantly contributed to recent advances in the field, but it is not straightforward to decide which data and methods are most suitable to use in a given case. In this review, we focus on these methods and data and their applications in generating MoA hypotheses for subsequent experimental validation. We discuss compound-specific data such as -omics, cell morphology and bioactivity data, as well as commonly used supplementary prior knowledge such as network and pathway data, and provide information on databases where this data can be accessed. In terms of methodologies, we discuss both well-established methods (connectivity mapping, pathway enrichment) as well as more developing methods (neural networks and multi-omics integration). Finally, we review case studies where the MoA of a compound was successfully suggested from computational analysis by incorporating multiple data modalities and/or methodologies. Our aim for this review is to provide researchers with insights into the benefits and drawbacks of both the data and methods in terms of level of understanding, biases and interpretation - and to highlight future avenues of investigation which we foresee will improve the field of MoA elucidation, including greater public access to -omics data and methodologies which are capable of data integration.

Drug target inference by mining transcriptional data using a novel graph convolutional network framework

A fundamental challenge that arises in biomedicine is the need to characterize compounds in a relevant cellular context in order to reveal potential on-target or off-target effects. Recently, the fast accumulation of gene transcriptional profiling data provides us an unprecedented opportunity to explore the protein targets of chemical compounds from the perspective of cell transcriptomics and RNA biology. Here, we propose a novel Siamese spectral-based graph convolutional network (SSGCN) model for inferring the protein targets of chemical compounds from gene transcriptional profiles. Although the gene signature of a compound perturbation only provides indirect clues of the interacting targets, and the biological networks under different experiment conditions further complicate the situation, the SSGCN model was successfully trained to learn from known compound-target pairs by uncovering the hidden correlations between compound perturbation profiles and gene knockdown profiles. On a benchmark set and a large time-split validation dataset, the model achieved higher target inference accuracy as compared to previous methods such as Connectivity Map. Further experimental validations of prediction results highlight the practical usefulness of SSGCN in either inferring the interacting targets of compound, or reversely, in finding novel inhibitors of a given target of interest.

Pathway-Based Drug Repurposing with DPNetinfer: A Method to Predict Drug-Pathway Associations via Network-Based Approaches

Identification of drug-pathway associations plays an important role in pathway-based drug repurposing. However, it is time-consuming and costly to uncover new drug-pathway associations experimentally. The drug-induced transcriptomics data provide a global view of cellular pathways and tell how these pathways change under different treatments. These data enable computational approaches for large-scale prediction of drug-pathway associations. Here we introduced DPNetinfer, a novel computational method to predict potential drug-pathway associations based on substructure-drug-pathway networks via network-based approaches. The results demonstrated that DPNetinfer performed well in a pan-cancer network with an AUC (area under curve) = 0.9358. Meanwhile, DPNetinfer was shown to have a good capability of generalization on two external validation sets (AUC = 0.8519 and 0.7494, respectively). As a case study, DPNetinfer was used in pathway-based drug repurposing for cancer therapy. Unexpected anticancer activities of some nononcology drugs were then identified on the PI3K-Akt pathway. Considering tumor heterogeneity, seven primary site-based models were constructed by DPNetinfer in different drug-pathway networks. In a word, DPNetinfer provides a powerful tool for large-scale prediction of drug-pathway associations in pathway-based drug repurposing. A web tool for DPNetinfer is freely available at http://lmmd.ecust.edu.cn/netinfer/.

TREAP: A New Topological Approach to Drug Target Inference

Over 50% of drugs fail in stage 3 clinical trials, many because of a poor understanding of the drug's mechanisms of action (MoA). A better comprehension of drug MoA will significantly improve research and development (R&D). Current proposed algorithms, such as ProTINA and DeMAND, can be overly complex. Additionally, they are unable to predict whether the drug-induced gene expression or the topology of the networks used to model gene regulation primarily impacts accurate drug target inference. In this work, we evaluate how network and gene expression data affect ProTINA's accuracy. We find that network topology predominantly determines the accuracy of ProTINA's predictions. We further show that the size of an interaction network and/or selecting cell-specific networks has a limited effect on accuracy. We then demonstrate that a specific network topology measure, betweenness, can be used to improve drug target prediction. Based on these results, we create a new algorithm, TREAP, that combines betweenness values and adjusted p-values for target inference. TREAP offers an alternative approach to drug target inference and is advantageous because it is not computationally demanding, provides easy-to-interpret results, and is often more accurate at predicting drug targets than current state-of-the-art approaches.

LiPLike: towards gene regulatory network predictions of high certainty

MOTIVATION

High correlation in expression between regulatory elements is a persistent obstacle for the reverse-engineering of gene regulatory networks. If two potential regulators have matching expression patterns, it becomes challenging to differentiate between them, thus increasing the risk of false positive identifications.

RESULTS

To allow for gene regulation predictions of high confidence, we propose a novel method, the Linear Profile Likelihood (LiPLike), that assumes a regression model and iteratively searches for interactions that cannot be replaced by a linear combination of other predictors. To compare the performance of LiPLike with other available inference methods, we benchmarked LiPLike using three independent datasets from the Dialogue on Reverse Engineering Assessment and Methods 5 (DREAM5) network inference challenge. We found that LiPLike could be used to stratify predictions of other inference tools, and when applied to the predictions of DREAM5 participants, we observed an average improvement in accuracy of >140% compared to individual methods. Furthermore, LiPLike was able to independently predict networks better than all DREAM5 participants when applied to biological data. When predicting the Escherichia coli network, LiPLike had an accuracy of 0.38 for the top-ranked 100 interactions, whereas the corresponding DREAM5 consensus model yielded an accuracy of 0.11.

AVAILABILITY AND IMPLEMENTATION

We made LiPLike available to the community as a Python toolbox, available at https://gitlab.com/Gustafsson-lab/liplike. We believe that LiPLike will be used for high confidence predictions in studies where individual model interactions are of high importance, and to remove false positive predictions made by other state-of-the-art gene-gene regulation prediction tools.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

An integrative method to predict signalling perturbations for cellular transitions

Induction of specific cellular transitions is of clinical importance, as it allows to revert disease cellular phenotype, or induce cellular reprogramming and differentiation for regenerative medicine. Signalling is a convenient way to accomplish such transitions without transfer of genetic material. Here we present the first general computational method that systematically predicts signalling molecules, whose perturbations induce desired cellular transitions. This probabilistic method integrates gene regulatory networks (GRNs) with manually-curated signalling pathways obtained from MetaCore from Clarivate Analytics, to model how signalling cues are received and processed in the GRN. The method was applied to 219 cellular transition examples, including cell type transitions, and overall correctly predicted experimentally validated signalling molecules, consistently outperforming other well-established approaches, such as differential gene expression and pathway enrichment analyses. Further, we validated our method predictions in the case of rat cirrhotic liver, and identified the activation of angiopoietins receptor Tie2 as a potential target for reverting the disease phenotype. Experimental results indicated that this perturbation induced desired changes in the gene expression of key TFs involved in fibrosis and angiogenesis. Importantly, this method only requires gene expression data of the initial and desired cell states, and therefore is suited for the discovery of signalling interventions for disease treatments and cellular therapies.

Inferring the transcriptional regulatory mechanism of signal-dependent gene expression via an integrative computational approach

Eukaryotic transcription factors (TFs) coordinate different upstream signals to regulate the expression of their target genes. To unveil this regulatory network in B-cell receptor signaling, we developed a computational pipeline to systematically analyze the extracellular signal-regulated kinase (ERK)- and IκB kinase (IKK)-dependent transcriptome responses. We combined a bilinear regression method and kinetic modeling to identify the signal-to-TF and TF-to-gene dynamics, respectively. We input a set of time-course experimental data for B cells and concentrated on transcriptional activators. The results show that the combination of TFs differentially controlled by ERK and IKK could contribute divergent expression dynamics in orchestrating the B-cell response. Our findings provide insights into the regulatory mechanisms underlying signal-dependent gene expression in eukaryotic cells.

SourceSet: A graphical model approach to identify primary genes in perturbed biological pathways

Topological gene-set analysis has emerged as a powerful means for omic data interpretation. Although numerous methods for identifying dysregulated genes have been proposed, few of them aim to distinguish genes that are the real source of perturbation from those that merely respond to the signal dysregulation. Here, we propose a new method, called SourceSet, able to distinguish between the primary and the secondary dysregulation within a Gaussian graphical model context. The proposed method compares gene expression profiles in the control and in the perturbed condition and detects the differences in both the mean and the covariance parameters with a series of likelihood ratio tests. The resulting evidence is used to infer the primary and the secondary set, i.e. the genes responsible for the primary dysregulation, and the genes affected by the perturbation through network propagation. The proposed method demonstrates high specificity and sensitivity in different simulated scenarios and on several real biological case studies. In order to fit into the more traditional pathway analysis framework, SourceSet R package also extends the analysis from a single to multiple pathways and provides several graphical outputs, including Cytoscape visualization to browse the results.

The rapid increase in omic studies has created a need to understand the biological implications of their results. Gene-set analysis has emerged as a powerful means for gaining such understanding, evolving in the last decade from the classical enrichment analysis to the more powerful topological approaches. Although numerous methods for identifying dysregulated genes have been proposed, few of them aim to distinguish genes that are the real source of perturbation from those that merely respond to the signal dysregulation. This distinction is crucial for network medicine, where the prioritization of the effect of biological perturbations may help in the molecular understanding of drug treatments and diseases. Here we propose a new method, called SourceSet, able to distinguish between primary and secondary dysregulation within a graphical model context, demonstrating a high specificity and sensitivity in different simulated scenarios and on real biological case studies.

DNA mismatch repair is required for the host innate response and controls cellular fate after influenza virus infection

Despite the cytopathic nature of influenza A virus (IAV) replication, we recently reported that a subset of lung epithelial club cells is able to intrinsically clear virus and survive infection. However, the mechanisms that drive cell survival during a normally lytic infection remained unclear. Using a loss-of-function screening approach, we discovered that the DNA mismatch repair (MMR) pathway is essential for club cell survival of IAV infection. Repair of virally-induced oxidative damage by the DNA MMR pathway not only allowed cell survival of infection but also facilitated host gene transcription, including the expression of antiviral and stress response genes. Enhanced viral suppression of the DNA MMR pathway prevented club cell survival and increased the severity of viral disease in vivo. Altogether, these results identify previously unappreciated roles for DNA MMR as a central modulator of cellular fate and a contributor to the innate antiviral response, which together, control influenza viral disease severity.

A dual controllability analysis of influenza virus-host protein-protein interaction networks for antiviral drug target discovery

Background

Host factors of influenza virus replication are often found in key topological positions within protein-protein interaction networks. This work explores how protein states can be manipulated through controllability analysis: the determination of the minimum manipulation needed to drive the cell system to any desired state. Here, we complete a two-part controllability analysis of two protein networks: a host network representing the healthy cell state and an influenza A virus-host network representing the infected cell state. In this context, controllability analyses aim to identify key regulating host factors of the infected cell’s progression. This knowledge can be utilized in further biological analysis to understand disease dynamics and isolate proteins for study as drug target candidates.

Results

Both topological and controllability analyses provide evidence of wide-reaching network effects stemming from the addition of viral-host protein interactions. Virus interacting and driver host proteins are significant both topologically and in controllability, therefore playing important roles in cell behavior during infection. Functional analysis finds overlap of results with previous siRNA studies of host factors involved in influenza replication, NF-kB pathway and infection relevance, and roles as interferon regulating genes. 24 proteins are identified as holding regulatory roles specific to the infected cell by measures of topology, controllability, and functional role. These proteins are recommended for further study as potential antiviral drug targets.

Conclusions

Seasonal outbreaks of influenza A virus are a major cause of illness and death around the world each year with a constant threat of pandemic infection. This research aims to increase the efficiency of antiviral drug target discovery using existing protein-protein interaction data and network analysis methods. These results are beneficial to future studies of influenza virus, both experimental and computational, and provide evidence that the combination of topology and controllability analyses may be valuable for future efforts in drug target discovery.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2917-z) contains supplementary material, which is available to authorized users.

Modeling gene-regulatory networks to describe cell fate transitions and predict master regulators

Complex organisms originate from and are maintained by the information encoded in the genome. A major challenge of systems biology is to develop algorithms that describe the dynamic regulation of genome functions from large omics datasets. Here, we describe TETRAMER, which reconstructs gene-regulatory networks from temporal transcriptome data during cell fate transitions to predict “master” regulators by simulating cascades of temporal transcription-regulatory events.

Boolean network modeling in systems pharmacology

Quantitative systems pharmacology (QSP) is an emerging discipline that aims to discover how drugs modulate the dynamics of biological components in molecular and cellular networks and the impact of those perturbations on human pathophysiology. The integration of systems-based experimental and computational approaches is required to facilitate the advancement of this field. QSP models typically consist of a series of ordinary differential equations (ODE). However, this mathematical framework requires extensive knowledge of parameters pertaining to biological processes, which is often unavailable. An alternative framework that does not require knowledge of system-specific parameters, such as Boolean network modeling, could serve as an initial foundation prior to the development of an ODE-based model. Boolean network models have been shown to efficiently describe, in a qualitative manner, the complex behavior of signal transduction and gene/protein regulatory processes. In addition to providing a starting point prior to quantitative modeling, Boolean network models can also be utilized to discover novel therapeutic targets and combinatorial treatment strategies. Identifying drug targets using a network-based approach could supplement current drug discovery methodologies and help to fill the innovation gap across the pharmaceutical industry. In this review, we discuss the process of developing Boolean network models and the various analyses that can be performed to identify novel drug targets and combinatorial approaches. An example for each of these analyses is provided using a previously developed Boolean network of signaling pathways in multiple myeloma. Selected examples of Boolean network models of human (patho-)physiological systems are also reviewed in brief.