Global network alignment is the computational problem of determining the similarity between nodes of different networks to establish a one-to-one correspondence between them. It has important applications in the biological field, particularly for discovering similar roles between the elements of different systems or for transferring knowledge from a well-studied system to another. In this paper, we present SHELLEY, a tool that facilitates the development, testing, and combination of learning-based network alignment algorithms by providing a set of modules that allow for the recreation and combination of both representation learning methods (RLMs) and deep matching methods (DMMs). We then present a case study in which we apply this tool to a protein-protein interaction network (PPI), demonstrating how the representation phase of RLMs is crucial for model robustness against noise.The code of SHELLEY is available at: https://github.com/rickydeluca/shelley

Early detection of prediabetes is crucial to preventing its progression to diabetes. Providing individuals with a personalized sense of their risk could improve prevention efforts. While complex mathematical models that simulate metabolic and inflammatory processes offer detailed and patient-specific insights, their computational cost usually makes them impractical for real-time prediction on mobile platforms.This work introduces a long short-term memory (LSTM) surrogate for the MT2D model, that simulates the main metabolic and inflammatory processes undergoing the transition to prediabetes. The model is developed using a dataset of 43 669 simulated subjects, each with lifestyle inputs and biomarker outputs over six months. Using 8 time series inputs, the surrogate predicts the dynamics of 11 key metabolic and inflammatory outputs, closely replicating the behaviour of the MT2D model.After training, the proposed LSTM model reduces computational time from an average of 8.4 hours to 0.1 seconds per simulation, making it suitable for mobile device deployment. The model achieves root mean squared errors on the order of 10−2on scaled data, and shows promise for prediabetes risk assessment by capturing trends in inflammatory biomarkers.This surrogate model can provide real-time and patient-specific insights into the metabolic health, potentially improving the understanding of prediabetes risk.

Drug repositioning is a promising strategy to discover new therapeutic applications for existing drugs, significantly reducing the time and costs associated with traditional drug development. This study employs a network medicine approach to analyze successful cases of drug repositioning, focusing on the exploratory hypothesis that the efficacy of repositioning may be determined by functional similarity between between diseases for which the drug was originally designed and diseases for which the same drug is reused. Network medicine tools were employed to investigate the connections between disease-associated genes, proteins, and approved drugs. Biological networks, including protein-protein interactions and functional interactions networks, as well as gene- and drug-disease association data are analyzed to identify functional similarities and possible molecular connections between diseases and treatments. Using clustering techniques and topological analysis, the results reveal a suggestive overlap of involved genes and functional interactions, emphasizing the value of computational methods in accelerating drug repositioning efforts and improving understanding of drug repositioning dynamics for more efficient therapeutic interventions.

The importance of the immune system (IS) in tuberculosis (TB) drug development is often underestimated because of the intricate nature of experiments and the specialized knowledge needed. In vitro and animal studies fall short in replicating the intricate reactions of the human IS to drugs and infections. In this study, we present our initial efforts in employing an in silico approach to comprehend how an individual’s IS impacts the efficacy of therapy, particularly in managing mycobacterium tuberculosis (Mtb) infection and minimizing the risk of relapse. We employed a well-established agent-based IS simulator called C-IMMSIM. We conducted simulations to investigate the long-term outcomes of TB disease in a virtual cohort infected with Mtb over a 50-year period. Our simulations revealed that individuals with competent IS showed a high success rate in containing Mtb infection. Furthermore, to better understand the dynamic interactions between Mtb and the IS, we deliberately introduced specific IS deficiencies, thus successfully inducing short-term relapses and mortality. These results confirm the model’s ability to elucidate the mechanisms underlying the interactions between Mtb and the IS.

MOTIVATION: Gene-disease associations are fundamental for understanding disease etiology and developing effective interventions and treatments. Identifying genes not yet associated with a disease due to a lack of studies is a challenging task in which prioritization based on prior knowledge is an important element. The computational search for new candidate disease genes may be eased by positive-unlabeled learning, the machine learning (ML) setting in which only a subset of instances are labeled as positive while the rest of the dataset is unlabeled. In this work, we propose a set of effective network-based features to be used in a novel Markov diffusion-based multi-class labeling strategy for putative disease gene discovery. RESULTS: The performances of the new labeling algorithm and the effectiveness of the proposed features have been tested on 10 different disease datasets using three ML algorithms. The new features have been compared against classical topological and functional/ontological features and a set of network- and biological-derived features already used in gene discovery tasks. The predictive power of the integrated methodology in searching for new disease genes has been found to be competitive against state-of-the-art algorithms. AVAILABILITY AND IMPLEMENTATION: The source code of NIAPU can be accessed at https://github.com/AndMastro/NIAPU. The source data used in this study are available online on the respective websites. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Drug repurposing is a highly active research area, aiming at finding novel uses for drugs that have been previously developed for other therapeutic purposes. Despite the flourishing of methodologies, success is still partial, and different approaches offer, each, peculiar advantages. In this composite landscape, we present a novel methodology focusing on an efficient mathematical procedure based on gene similarity scores and biased random walks which rely on robust drug-gene-disease association data sets. The recommendation mechanism is further unveiled by means of the Markov chain underlying the random walk process, hence providing explainability about how findings are suggested. Performances evaluation and the analysis of a case study on rheumatoid arthritis show that our approach is accurate in providing useful recommendations and is computationally efficient, compared to the state of the art of drug repurposing approaches.

The importance of faster drug development has never been more evident than in present time when the whole world is struggling to cope up with the COVID-19 pandemic. At times when timely development of effective drugs and treatment plans could potentially save millions of lives, drug repurposing is one area of medicine that has garnered much of research interest. Apart from experimental drug repurposing studies that happen within wet labs, lot many new quantitative methods have been proposed in the literature. In this paper, one such quantitative methods for drug repurposing is implemented and evaluated. DruSiLa (DRUg in-SIlico LAboratory) is an in-silico drug repurposing method that leverages disease similarity measures to quantitatively rank existing drugs for their potential therapeutic efficacy against novel diseases. The proposed method makes use of available, manually curated, and open datasets on diseases, their genetic origins, and disease-related patho-phenotypes. DruSiLa evaluates pairwise disease similarity scores of any given target disease to each known disease in our dataset. Such similarity scores are then propagated through disease-drug associations, and aggregated at drug nodes to rank them for their predicted effectiveness against the target disease.

The possibility to computationally prioritize candidate disease genes capitalizing on existing information has led to a speedup in the discovery of new methods. Many gene discovery techniques exploit network data, like protein-protein interactions (PPIs), in order to extract knowledge from the network structure relying on several network metrics. We here present PROCONSUL, a method that builds on top of the concept of connectivity significance (CS) and exploits the idea of probabilistic exploration of the space of putative disease genes. We show that our methodology is able to outperform the state-of-the-art tool based on CS in several settings, and propose different, effective gene discovery strategies according to specific disease network properties.

Systems medicine (SM) has emerged as a powerful tool for studying the human body at the systems level with the aim of improving our understanding, prevention and treatment of complex diseases. Being able to automatically extract relevant features needed for a given task from high-dimensional, heterogeneous data, deep learning (DL) holds great promise in this endeavour. This review paper addresses the main developments of DL algorithms and a set of general topics where DL is decisive, namely, within the SM landscape. It discusses how DL can be applied to SM with an emphasis on the applications to predictive, preventive and precision medicine. Several key challenges have been highlighted including delivering clinical impact and improving interpretability. We used some prototypical examples to highlight the relevance and significance of the adoption of DL in SM, one of them is involving the creation of a model for personalized Parkinson's disease. The review offers valuable insights and informs the research in DL and SM.

MOTIVATION: Unsupervised learning approaches are frequently used to stratify patients into clinically relevant subgroups and to identify biomarkers such as disease-associated genes. However, clustering and biclustering techniques are oblivious to the functional relationship of genes and are thus not ideally suited to pinpoint molecular mechanisms along with patient subgroups. RESULTS: We developed the network-constrained biclustering approach Biclustering Constrained by Networks (BiCoN) which (i) restricts biclusters to functionally related genes connected in molecular interaction networks and (ii) maximizes the difference in gene expression between two subgroups of patients. This allows BiCoN to simultaneously pinpoint molecular mechanisms responsible for the patient grouping. Network-constrained clustering of genes makes BiCoN more robust to noise and batch effects than typical clustering and biclustering methods. BiCoN can faithfully reproduce known disease subtypes as well as novel, clinically relevant patient subgroups, as we could demonstrate using breast and lung cancer datasets. In summary, BiCoN is a novel systems medicine tool that combines several heuristic optimization strategies for robust disease mechanism extraction. BiCoN is well-documented and freely available as a python package or a web interface. AVAILABILITY AND IMPLEMENTATION: PyPI package: https://pypi.org/project/bicon. WEB INTERFACE: https://exbio.wzw.tum.de/bicon. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SUMMARY: Simulated data are crucial for evaluating epistasis detection tools in genome-wide association studies. Existing simulators are limited, as they do not account for linkage disequilibrium (LD), support limited interaction models of single nucleotide polymorphisms (SNPs) and only dichotomous phenotypes or depend on proprietary software. In contrast, EpiGEN supports SNP interactions of arbitrary order, produces realistic LD patterns and generates both categorical and quantitative phenotypes. AVAILABILITY AND IMPLEMENTATION: EpiGEN is implemented in Python 3 and is freely available at https://github.com/baumbachlab/epigen. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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 .

Network-based ranking methods (e.g., centrality analysis) have found extensive use in systems biology and network medicine for the prediction of essential proteins, for the prioritization of drug targets candidates in the treatment of several pathologies and in biomarker discovery, and for human disease genes identification. We here studied the connectivity of the human protein-protein interaction network (i.e., the interactome) to find the nodes whose removal has the heaviest impact on the network, i.e., maximizes its fragmentation. Such nodes are known as Critical Nodes (CNs). Specifically, we implemented a Critical Node Heuristic (CNH) and compared its performance against other four heuristics based on well known centrality measures. To better understand the structure of the interactome, the CNs' role played in the network, and the different heuristics' capabilities to grasp biologically relevant nodes, we compared the sets of nodes identified as CNs by each heuristic with two experimentally validated sets of essential genes, i.e., the genes whose removal impact on a given organism's ability to survive. Our results show that classical centrality measures (i.e., closeness centrality, degree) found more essential genes with respect to CNH on the current version of the human interactome, however the removal of such nodes does not have the greatest impact on interactome connectivity, while, interestingly, the genes identified by CNH show peculiar characteristics both from the topological and the biological point of view. Finally, even if a relevant fraction of essential genes is found via the classical centrality measures, the same measures seem to fail in identifying the whole set of essential genes, suggesting once again that some of them are not central in the network, that there may be biases in the current interaction data, and that different, combined graph theoretical and other techniques should be applied for their discovery.

Investigation about the mechanisms involved in the onset of type 2 diabetes in absence of familiarity is the focus of a research project which has led to the development of a computational model that recapitulates the aetiology of the disease. The model simulates the metabolic and immunological alterations related to type-2 diabetes associated to several clinical, physiological and behavioural characteristics of representative virtual patients. In this study, the results of 46170 simulations corresponding to the same number of virtual subjects, experiencing different lifestyle conditions, are analysed for the construction of a statistical model able to recapitulate the simulated dynamics. The resulting machine learning model adequately predicts the synthetic data 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.

Systems medicine holds many promises, but has so far provided only a limited number of proofs of principle. To address this road block, possible barriers and challenges of translating systems medicine into clinical practice need to be identified and addressed. The members of the European Cooperation in Science and Technology (COST) Action CA15120 Open Multiscale Systems Medicine (OpenMultiMed) wish to engage the scientific community of systems medicine and multiscale modelling, data science and computing, to provide their feedback in a structured manner. This will result in follow-up white papers and open access resources to accelerate the clinical translation of systems medicine.

The beneficial effects of physical activity for the prevention and management of several chronic diseases are widely recognized.Mathematical modeling of the effects of physical exercise in body metabolism and in particular its influence on the control of glucose homeostasis is of primary importance in the development of eHealth monitoring devices for a personalized medicine.Nonetheless, to date only a few mathematical models have been aiming at this specific purpose.We have developed a whole-body computational model of the effects on metabolic homeostasis of a bout of physical exercise.Built upon an existing model, it allows to detail better both subjects' characteristics and physical exercise, thus determining to a greater extent the dynamics of the hormones and the metabolites considered. Author summaryExercise has a great impact on human metabolism and the lack of physical activity represents one of the main causes of the metabolic disorders.The effectiveness of regular physical activity in the prevention and management of several chronic diseases is widely recognized.In the study of the metabolism and related disorders, mathematical models have proven useful in describing and quantifying physiological processes often not easily measurable in vivo.Formulating a model describing the metabolic responses to a physical exercise session is a challenging task since the effects vary depending on its intensity, duration, modality and are also dependent on the subjects' physical characteristics (e.g.age, gender, body weight, fitness status).To date, none of the existing computational models is able to provide this level of "personalization".Thus, starting from an existing model of fuel homeostasis during exercise, we have formulated a novel computational system that is more detailed in describing both the physical exercise and the subjects' characteristics.

One of the greatest challenges in biomedicine is to get a unified view of observations made from the molecular up to the organism scale. Towards this goal, multiscale models have been highly instrumental in contexts such as the cardiovascular field, angiogenesis, neurosciences and tumour biology. More recently, such models are becoming an increasingly important resource to address immunological questions as well. Systematic mining of the literature in multiscale modelling led us to identify three main fields of immunological applications: host-virus interactions, inflammatory diseases and their treatment and development of multiscale simulation platforms for immunological research and for educational purposes. Here, we review the current developments in these directions, which illustrate that multiscale models can consistently integrate immunological data generated at several scales, and can be used to describe and optimize therapeutic treatments of complex immune diseases.

BACKGROUND: Macrophages cover a major role in the immune system, being the most plastic cell yielding several key immune functions. METHODS: Here we derived a minimalistic gene regulatory network model for the differentiation of macrophages into the two phenotypes M1 (pro-) and M2 (anti-inflammatory). RESULTS: To test the model, we simulated a large number of such networks as in a statistical ensemble. In other words, to enable the inter-cellular crosstalk required to obtain an immune activation in which the macrophage plays its role, the simulated networks are not taken in isolation but combined with other cellular agents, thus setting up a discrete minimalistic model of the immune system at the microscopic/intracellular (i.e., genetic regulation) and mesoscopic/intercellular scale. CONCLUSIONS: We show that within the mesoscopic level description of cellular interaction and cooperation, the gene regulatory logic is coherent and contributes to the overall dynamics of the ensembles that shows, statistically, the expected behaviour.

MOTIVATION: Immune cells coordinate their efforts for the correct and efficient functioning of the immune system (IS). Each cell type plays a distinct role and communicates with other cell types through mediators such as cytokines, chemokines and hormones, among others, that are crucial for the functioning of the IS and its fine tuning. Nevertheless, a quantitative analysis of the topological properties of an immunological network involving this complex interchange of mediators among immune cells is still lacking. RESULTS: Here we present a method for quantifying the relevance of different mediators in the immune network, which exploits a definition of centrality based on the concept of efficient communication. The analysis, applied to the human IS, indicates that its mediators differ significantly in their network relevance. We found that cytokines involved in innate immunity and inflammation and some hormones rank highest in the network, revealing that the most prominent mediators of the IS are molecules involved in these ancestral types of defence mechanisms which are highly integrated with the adaptive immune response, and at the interplay among the nervous, the endocrine and the immune systems. CONTACT: [email protected].