Cell reprogramming design by transfer learning of functional transcriptional networks.

Recent developments in synthetic biology, next-generation sequencing, and machine learning provide an unprecedented opportunity to rationally design new disease treatments based on measured responses to gene perturbations and drugs to reprogram cells. The main challenges to seizing this opportunity are the incomplete knowledge of the cellular network and the combinatorial explosion of possible interventions, both of which are insurmountable by experiments. To address these challenges, we develop a transfer learning approach to control cell behavior that is pre-trained on transcriptomic data associated with human cell fates, thereby generating a model of the network dynamics that can be transferred to specific reprogramming goals. The approach combines transcriptional responses to gene perturbations to minimize the difference between a given pair of initial and target transcriptional states. We demonstrate our approach's versatility by applying it to a microarray dataset comprising >9,000 microarrays across 54 cell types and 227 unique perturbations, and an RNASeq dataset consisting of >10,000 sequencing runs across 36 cell types and 138 perturbations. Our approach reproduces known reprogramming protocols with an AUROC of 0.91 while innovating over existing methods by pre-training an adaptable model that can be tailored to specific reprogramming transitions. We show that the number of gene perturbations required to steer from one fate to another increases with decreasing developmental relatedness and that fewer genes are needed to progress along developmental paths than to regress. These findings establish a proof-of-concept for our approach to computationally design control strategies and provide insights into how gene regulatory networks govern phenotype.

BPP: a platform for automatic biochemical pathway prediction.

A biochemical pathway consists of a series of interconnected biochemical reactions to accomplish specific life activities. The participating reactants and resultant products of a pathway, including gene fragments, proteins, and small molecules, coalesce to form a complex reaction network. Biochemical pathways play a critical role in the biochemical domain as they can reveal the flow of biochemical reactions in living organisms, making them essential for understanding life processes. Existing studies of biochemical pathway networks are mainly based on experimentation and pathway database analysis methods, which are plagued by substantial cost constraints. Inspired by the success of representation learning approaches in biomedicine, we develop the biochemical pathway prediction (BPP) platform, which is an automatic BPP platform to identify potential links or attributes within biochemical pathway networks. Our BPP platform incorporates a variety of representation learning models, including the latest hypergraph neural networks technology to model biochemical reactions in pathways. In particular, BPP contains the latest biochemical pathway-based datasets and enables the prediction of potential participants or products of biochemical reactions in biochemical pathways. Additionally, BPP is equipped with an SHAP explainer to explain the predicted results and to calculate the contributions of each participating element. We conduct extensive experiments on our collected biochemical pathway dataset to benchmark the effectiveness of all models available on BPP. Furthermore, our detailed case studies based on the chronological pattern of our dataset demonstrate the effectiveness of our platform. Our BPP web portal, source code and datasets are freely accessible at https://github.com/Glasgow-AI4BioMed/BPP.

Leveraging developmental landscapes for model selection in Boolean gene regulatory networks.

Boolean models are a well-established framework to model developmental gene regulatory networks (DGRNs) for acquisition of cellular identities. During the reconstruction of Boolean DGRNs, even if the network structure is given, there is generally a large number of combinations of Boolean functions that will reproduce the different cell fates (biological attractors). Here we leverage the developmental landscape to enable model selection on such ensembles using the relative stability of the attractors. First we show that previously proposed measures of relative stability are strongly correlated and we stress the usefulness of the one that captures best the cell state transitions via the mean first passage time (MFPT) as it also allows the construction of a cellular lineage tree. A property of great computational importance is the insensitivity of the different stability measures to changes in noise intensities. That allows us to use stochastic approaches to estimate the MFPT and thereby scale up the computations to large networks. Given this methodology, we revisit different Boolean models of Arabidopsis thaliana root development, showing that a most recent one does not respect the biologically expected hierarchy of cell states based on relative stabilities. We therefore developed an iterative greedy algorithm that searches for models which satisfy the expected hierarchy of cell states and found that its application to the root development model yields many models that meet this expectation. Our methodology thus provides new tools that can enable reconstruction of more realistic and accurate Boolean models of DGRNs.

PanomiR: a systems biology framework for analysis of multi-pathway targeting by miRNAs.

Charting microRNA (miRNA) regulation across pathways is key to characterizing their function. Yet, no method currently exists that can quantify how miRNAs regulate multiple interconnected pathways or prioritize them for their ability to regulate coordinate transcriptional programs. Existing methods primarily infer one-to-one relationships between miRNAs and pathways using differentially expressed genes. We introduce PanomiR, an in silico framework for studying the interplay of miRNAs and disease functions. PanomiR integrates gene expression, mRNA-miRNA interactions and known biological pathways to reveal coordinated multi-pathway targeting by miRNAs. PanomiR utilizes pathway-activity profiling approaches, a pathway co-expression network and network clustering algorithms to prioritize miRNAs that target broad-scale transcriptional disease phenotypes. It directly resolves differential regulation of pathways, irrespective of their differential gene expression, and captures co-activity to establish functional pathway groupings and the miRNAs that may regulate them. PanomiR uses a systems biology approach to provide broad but precise insights into miRNA-regulated functional programs. It is available at https://bioconductor.org/packages/PanomiR.

Improving therapeutic synergy score predictions with adverse effects using multi-task heterogeneous network learning.

Drug combinations could trigger pharmacological therapeutic effects (TEs) and adverse effects (AEs). Many computational methods have been developed to predict TEs, e.g. the therapeutic synergy scores of anti-cancer drug combinations, or AEs from drug-drug interactions. However, most of the methods treated the AEs and TEs predictions as two separate tasks, ignoring the potential mechanistic commonalities shared between them. Based on previous clinical observations, we hypothesized that by learning the shared mechanistic commonalities between AEs and TEs, we could learn the underlying MoAs (mechanisms of actions) and ultimately improve the accuracy of TE predictions. To test our hypothesis, we formulated the TE prediction problem as a multi-task heterogeneous network learning problem that performed TE and AE learning tasks simultaneously. To solve this problem, we proposed Muthene (multi-task heterogeneous network embedding) and evaluated it on our collected drug-drug interaction dataset with both TEs and AEs indications. Our experimental results showed that, by including the AE prediction as an auxiliary task, Muthene generated more accurate TE predictions than standard single-task learning methods, which supports our hypothesis. Using a drug pair Vincristine-Dasatinib as a case study, we demonstrated that our method not only provides a novel way of TE predictions but also helps us gain a deeper understanding of the MoAs of drug combinations.

Network integration of thermal proteome profiling with multi-omics data decodes PARP inhibition.

Complex disease phenotypes often span multiple molecular processes. Functional characterization of these processes can shed light on disease mechanisms and drug effects. Thermal Proteome Profiling (TPP) is a mass-spectrometry (MS) based technique assessing changes in thermal protein stability that can serve as proxies of functional protein changes. These unique insights of TPP can complement those obtained by other omics technologies. Here, we show how TPP can be integrated with phosphoproteomics and transcriptomics in a network-based approach using COSMOS, a multi-omics integration framework, to provide an integrated view of transcription factors, kinases and proteins with altered thermal stability. This allowed us to recover consequences of Poly (ADP-ribose) polymerase (PARP) inhibition in ovarian cancer cells on cell cycle and DNA damage response as well as interferon and hippo signaling. We found that TPP offers a complementary perspective to other omics data modalities, and that its integration allowed us to obtain a more complete molecular overview of PARP inhibition. We anticipate that this strategy can be used to integrate functional proteomics with other omics to study molecular processes.

Computational modeling of protracted HCMV replication using genome substrates and protein temporal profiles.

Human cytomegalovirus (HCMV) is a major cause of illness in immunocompromised individuals. The HCMV lytic cycle contributes to the clinical manifestations of infection. The lytic cycle occurs over ∼96 h in diverse cell types and consists of viral DNA (vDNA) genome replication and temporally distinct expression of hundreds of viral proteins. Given its complexity, understanding this elaborate system can be facilitated by the introduction of mechanistic computational modeling of temporal relationships. Therefore, we developed a multiplicity of infection (MOI)-dependent mechanistic computational model that simulates vDNA kinetics and late lytic replication based on in-house experimental data. The predictive capabilities were established by comparison to post hoc experimental data. Computational analysis of combinatorial regulatory mechanisms suggests increasing rates of protein degradation in association with increasing vDNA levels. The model framework also allows expansion to account for additional mechanisms regulating the processes. Simulating vDNA kinetics and the late lytic cycle for a wide range of MOIs yielded several unique observations. These include the presence of saturation behavior at high MOIs, inefficient replication at low MOIs, and a precise range of MOIs in which virus is maximized within a cell type, being 0.382 IU to 0.688 IU per fibroblast. The predicted saturation kinetics at high MOIs are likely related to the physical limitations of cellular machinery, while inefficient replication at low MOIs may indicate a minimum input material required to facilitate infection. In summary, we have developed and demonstrated the utility of a data-driven and expandable computational model simulating lytic HCMV infection.

A computational model of circRNA-associated diseases based on a graph neural network: prediction and case studies for follow-up experimental validation.

BACKGROUND: Circular RNAs (circRNAs) have been confirmed to play a vital role in the occurrence and development of diseases. Exploring the relationship between circRNAs and diseases is of far-reaching significance for studying etiopathogenesis and treating diseases. To this end, based on the graph Markov neural network algorithm (GMNN) constructed in our previous work GMNN2CD, we further considered the multisource biological data that affects the association between circRNA and disease and developed an updated web server CircDA and based on the human hepatocellular carcinoma (HCC) tissue data to verify the prediction results of CircDA. RESULTS: CircDA is built on a Tumarkov-based deep learning framework. The algorithm regards biomolecules as nodes and the interactions between molecules as edges, reasonably abstracts multiomics data, and models them as a heterogeneous biomolecular association network, which can reflect the complex relationship between different biomolecules. Case studies using literature data from HCC, cervical, and gastric cancers demonstrate that the CircDA predictor can identify missing associations between known circRNAs and diseases, and using the quantitative real-time PCR (RT-qPCR) experiment of HCC in human tissue samples, it was found that five circRNAs were significantly differentially expressed, which proved that CircDA can predict diseases related to new circRNAs. CONCLUSIONS: This efficient computational prediction and case analysis with sufficient feedback allows us to identify circRNA-associated diseases and disease-associated circRNAs. Our work provides a method to predict circRNA-associated diseases and can provide guidance for the association of diseases with certain circRNAs. For ease of use, an online prediction server ( http://server.malab.cn/CircDA ) is provided, and the code is open-sourced ( https://github.com/nmt315320/CircDA.git ) for the convenience of algorithm improvement.

Mechanistic interactions as the origin of modularity in biological networks.

Biological networks are often modular. Explanations for this peculiarity either assume an adaptive advantage of a modular design such as higher robustness, or attribute it to neutral factors such as constraints underlying network assembly. Interestingly, most insights on the origin of modularity stem from models in which interactions are either determined by highly simplistic mechanisms, or have no mechanistic basis at all. Yet, empirical knowledge suggests that biological interactions are often mediated by complex structural or behavioural traits. Here, we investigate the origins of modularity using a model in which interactions are determined by potentially complex traits. Specifically, we model system elements-such as the species in an ecosystem-as finite-state machines (FSMs), and determine their interactions by means of communication between the corresponding FSMs. Using this model, we show that modularity probably emerges for free. We further find that the more modular an interaction network is, the less complex are the traits that mediate the interactions. Altogether, our results suggest that the conditions for modularity to evolve may be much broader than previously thought.

Detection of pan-cancer surface protein biomarkers via a network-based approach on transcriptomics data.

Cell surface proteins have been used as diagnostic and prognostic markers in cancer research and as targets for the development of anticancer agents. Many of these proteins lie at the top of signaling cascades regulating cell responses and gene expression, therefore acting as 'signaling hubs'. It has been previously demonstrated that the integrated network analysis on transcriptomic data is able to infer cell surface protein activity in breast cancer. Such an approach has been implemented in a publicly available method called 'SURFACER'. SURFACER implements a network-based analysis of transcriptomic data focusing on the overall activity of curated surface proteins, with the final aim to identify those proteins driving major phenotypic changes at a network level, named surface signaling hubs. Here, we show the ability of SURFACER to discover relevant knowledge within and across cancer datasets. We also show how different cancers can be stratified in surface-activity-specific groups. Our strategy may identify cancer-wide markers to design targeted therapies and biomarker-based diagnostic approaches.

Prospects and challenges of cancer systems medicine: from genes to disease networks.

It is becoming evident that holistic perspectives toward cancer are crucial in deciphering the overwhelming complexity of tumors. Single-layer analysis of genome-wide data has greatly contributed to our understanding of cellular systems and their perturbations. However, fundamental gaps in our knowledge persist and hamper the design of effective interventions. It is becoming more apparent than ever, that cancer should not only be viewed as a disease of the genome but as a disease of the cellular system. Integrative multilayer approaches are emerging as vigorous assets in our endeavors to achieve systemic views on cancer biology. Herein, we provide a comprehensive review of the approaches, methods and technologies that can serve to achieve systemic perspectives of cancer. We start with genome-wide single-layer approaches of omics analyses of cellular systems and move on to multilayer integrative approaches in which in-depth descriptions of proteogenomics and network-based data analysis are provided. Proteogenomics is a remarkable example of how the integration of multiple levels of information can reduce our blind spots and increase the accuracy and reliability of our interpretations and network-based data analysis is a major approach for data interpretation and a robust scaffold for data integration and modeling. Overall, this review aims to increase cross-field awareness of the approaches and challenges regarding the omics-based study of cancer and to facilitate the necessary shift toward holistic approaches.

Network approaches for modeling the effect of drugs and diseases.

The network approach is quickly becoming a fundamental building block of computational methods aiming at elucidating the mechanism of action (MoA) and therapeutic effect of drugs. By modeling the effect of drugs and diseases on different biological networks, it is possible to better explain the interplay between disease perturbations and drug targets as well as how drug compounds induce favorable biological responses and/or adverse effects. Omics technologies have been extensively used to generate the data needed to study the mechanisms of action of drugs and diseases. These data are often exploited to define condition-specific networks and to study whether drugs can reverse disease perturbations. In this review, we describe network data mining algorithms that are commonly used to study drug's MoA and to improve our understanding of the basis of chronic diseases. These methods can support fundamental stages of the drug development process, including the identification of putative drug targets, the in silico screening of drug compounds and drug combinations for the treatment of diseases. We also discuss recent studies using biological and omics-driven networks to search for possible repurposed FDA-approved drug treatments for SARS-CoV-2 infections (COVID-19).

Topological ranks reveal functional knowledge encoded in biological networks: a comparative analysis.

MOTIVATION: Biological networks topology yields important insights into biological function, occurrence of diseases and drug design. In the last few years, different types of topological measures have been introduced and applied to infer the biological relevance of network components/interactions, according to their position within the network structure. Although comparisons of such measures have been previously proposed, to what extent the topology per se may lead to the extraction of novel biological knowledge has never been critically examined nor formalized in the literature. RESULTS: We present a comparative analysis of nine outstanding topological measures, based on compact views obtained from the rank they induce on a given input biological network. The goal is to understand their ability in correctly positioning nodes/edges in the rank, according to the functional knowledge implicitly encoded in biological networks. To this aim, both internal and external (gold standard) validation criteria are taken into account, and six networks involving three different organisms (yeast, worm and human) are included in the comparison. The results show that a distinct handful of best-performing measures can be identified for each of the considered organisms, independently from the reference gold standard. AVAILABILITY: Input files and code for the computation of the considered topological measures and K-haus distance are available at https://gitlab.com/MaryBonomo/ranking. CONTACT: simona.rombo@unipa.it. SUPPLEMENTARY INFORMATION: Supplementary data are available at Briefings in Bioinformatics online.

HyMM: hybrid method for disease-gene prediction by integrating multiscale module structure.

MOTIVATION: Identifying disease-related genes is an important issue in computational biology. Module structure widely exists in biomolecule networks, and complex diseases are usually thought to be caused by perturbations of local neighborhoods in the networks, which can provide useful insights for the study of disease-related genes. However, the mining and effective utilization of the module structure is still challenging in such issues as a disease gene prediction. RESULTS: We propose a hybrid disease-gene prediction method integrating multiscale module structure (HyMM), which can utilize multiscale information from local to global structure to more effectively predict disease-related genes. HyMM extracts module partitions from local to global scales by multiscale modularity optimization with exponential sampling, and estimates the disease relatedness of genes in partitions by the abundance of disease-related genes within modules. Then, a probabilistic model for integration of gene rankings is designed in order to integrate multiple predictions derived from multiscale module partitions and network propagation, and a parameter estimation strategy based on functional information is proposed to further enhance HyMM's predictive power. By a series of experiments, we reveal the importance of module partitions at different scales, and verify the stable and good performance of HyMM compared with eight other state-of-the-arts and its further performance improvement derived from the parameter estimation. CONCLUSIONS: The results confirm that HyMM is an effective framework for integrating multiscale module structure to enhance the ability to predict disease-related genes, which may provide useful insights for the study of the multiscale module structure and its application in such issues as a disease-gene prediction.

Phenotypic plasticity for improved light harvesting, in tandem with methylome repatterning in reef-building corals.

Acclimatization through phenotypic plasticity represents a more rapid response to environmental change than adaptation and is vital to optimize organisms' performance in different conditions. Generally, animals are less phenotypically plastic than plants, but reef-building corals exhibit plant-like properties. They are light dependent with a sessile and modular construction that facilitates rapid morphological changes within their lifetime. We induced phenotypic changes by altering light exposure in a reciprocal transplant experiment and found that coral plasticity is a colony trait emerging from comprehensive morphological and physiological changes within the colony. Plasticity in skeletal features optimized coral light harvesting and utilization and paralleled significant methylome and transcriptome modifications. Network-associated responses resulted in the identification of hub genes and clusters associated to the change in phenotype: inter-partner recognition and phagocytosis, soft tissue growth and biomineralization. Furthermore, we identified hub genes putatively involved in animal photoreception-phototransduction. These findings fundamentally advance our understanding of how reef-building corals repattern the methylome and adjust a phenotype, revealing an important role of light sensing by the coral animal to optimize photosynthetic performance of the symbionts.

Dynamic regulation of RAS and RAS signaling.

RAS proteins regulate most aspects of cellular physiology. They are mutated in 30% of human cancers and 4% of developmental disorders termed Rasopathies. They cycle between active GTP-bound and inactive GDP-bound states. When active, they can interact with a wide range of effectors that control fundamental biochemical and biological processes. Emerging evidence suggests that RAS proteins are not simple on/off switches but sophisticated information processing devices that compute cell fate decisions by integrating external and internal cues. A critical component of this compute function is the dynamic regulation of RAS activation and downstream signaling that allows RAS to produce a rich and nuanced spectrum of biological outputs. We discuss recent findings how the dynamics of RAS and its downstream signaling is regulated. Starting from the structural and biochemical properties of wild-type and mutant RAS proteins and their activation cycle, we examine higher molecular assemblies, effector interactions and downstream signaling outputs, all under the aspect of dynamic regulation. We also consider how computational and mathematical modeling approaches contribute to analyze and understand the pleiotropic functions of RAS in health and disease.

A guide to ERK dynamics, part 1: mechanisms and models.

Extracellular signal-regulated kinase (ERK) has long been studied as a key driver of both essential cellular processes and disease. A persistent question has been how this single pathway is able to direct multiple cell behaviors, including growth, proliferation, and death. Modern biosensor studies have revealed that the temporal pattern of ERK activity is highly variable and heterogeneous, and critically, that these dynamic differences modulate cell fate. This two-part review discusses the current understanding of dynamic activity in the ERK pathway, how it regulates cellular decisions, and how these cell fates lead to tissue regulation and pathology. In part 1, we cover the optogenetic and live-cell imaging technologies that first revealed the dynamic nature of ERK, as well as current challenges in biosensor data analysis. We also discuss advances in mathematical models for the mechanisms of ERK dynamics, including receptor-level regulation, negative feedback, cooperativity, and paracrine signaling. While hurdles still remain, it is clear that higher temporal and spatial resolution provide mechanistic insights into pathway circuitry. Exciting new algorithms and advanced computational tools enable quantitative measurements of single-cell ERK activation, which in turn inform better models of pathway behavior. However, the fact that current models still cannot fully recapitulate the diversity of ERK responses calls for a deeper understanding of network structure and signal transduction in general.

The narrowing of dendrite branches across nodes follows a well-defined scaling law.

The systematic variation of diameters in branched networks has tantalized biologists since the discovery of da Vinci's rule for trees. Da Vinci's rule can be formulated as a power law with exponent two: The square of the mother branch's diameter is equal to the sum of the squares of those of the daughters. Power laws, with different exponents, have been proposed for branching in circulatory systems (Murray's law with exponent 3) and in neurons (Rall's law with exponent 3/2). The laws have been derived theoretically, based on optimality arguments, but, for the most part, have not been tested rigorously. Using superresolution methods to measure the diameters of dendrites in highly branched Drosophila class IV sensory neurons, we have found that these types of power laws do not hold. In their place, we have discovered a different diameter-scaling law: The cross-sectional area is proportional to the number of dendrite tips supported by the branch plus a constant, corresponding to a minimum diameter of the terminal dendrites. The area proportionality accords with a requirement for microtubules to transport materials and nutrients for dendrite tip growth. The minimum diameter may be set by the force, on the order of a few piconewtons, required to bend membrane into the highly curved surfaces of terminal dendrites. Because the observed scaling differs from Rall's law, we propose that cell biological constraints, such as intracellular transport and protrusive forces generated by the cytoskeleton, are important in determining the branched morphology of these cells.

Possible New Candidates Involved to Thalidomide-Related Limbs and Cardiac Defects: A Systems Biology Approach.

Thalidomide is a known teratogen that causes malformations especially in heart and limbs. Its mechanism of teratogenicity is still not fully elucidated. Recently, a new target of thalidomide was described, TBX5, and was observed a new interaction between HAND2 and TBX5 that is disrupted in the presence of thalidomide. Therefore, our study aimed to raise potential candidates for thalidomide teratogenesis, through systems biology, evaluating HAND2 and TBX5 interaction and heart and limbs malformations of thalidomide. Genes and proteins related to TBX5 and HAND2 were selected through TF2DNA, REACTOME, Human Phenotype Ontology, and InterPro databases. Networks were assembled using STRING © database. Network analysis were performed in Cytoscape © and R v3.6.2. Differential gene expression (DGE) analysis was performed through gene expression omnibus. We constructed a network for HAND2 and TBX5 interaction; a network for heart and limbs malformations of TE; and the two joined networks. We observed that EP300 protein seemed to be important in all networks. We also looked for proteins containing C2H2 domain in the assembled networks. ZIC3, GLI1, GLI3, ZNF148, and PRDM16 were the ones present in both heart and limbs malformations of TE networks. Furthermore, in the DGE analysis after treatment with thalidomide, we observed that FANCB, ESCO2, and XRCC2 were downregulated and present both in heart and limbs networks. Through systems biology, we were able to point to different new proteins and genes, and selected specially EP300, which was important in all the analyzed networks, to be further evaluated in the TE teratogenicity.

Gene network modeling and pathway analysis of maize transcriptomes in response to Maize Iranian mosaic virus.

Maize Iranian mosaic virus (MIMV, family Rhabdoviridae) is one of the factors limiting cereal production in Iran. In the present study, we sought to find critical genes and key pathways involved in MIMV infection and analyzed gene networks, pathways and promoters using transcriptome data. We determined the hub genes involved in pathways related to the proteasome and ubiquitin. The results showed the important role of the cellular endoplasmic reticulum in MIMV infection. Network cluster analysis confirmed the result of GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The discovered miRNAs belonged to miR166, miR167, miR169, miR395, miR399, miR408 and miR482 families, which are involved in various pathogenicity or resistance processes against MIMV or other viruses. The results of this study provide a list of hub genes, important pathways and new insights for the future development of virus-resistant transgenic crops and clarify the basic mechanism of plant response.