Integrative network analysis suggests prioritised drugs for atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD) is a prevalent chronic inflammatory skin disease whose pathophysiology involves the interplay between genetic and environmental factors, ultimately leading to dysfunction of the epidermis. While several treatments are effective in symptom management, many existing therapies offer only temporary relief and often come with side effects. For this reason, the formulation of an effective therapeutic plan is challenging and there is a need for more effective and targeted treatments that address the root causes of the condition. Here, we hypothesise that modelling the complexity of the molecular buildup of the atopic dermatitis can be a concrete means to drive drug discovery. METHODS: We preprocessed, harmonised and integrated publicly available transcriptomics datasets of lesional and non-lesional skin from AD patients. We inferred co-expression network models of both AD lesional and non-lesional skin and exploited their interactional properties by integrating them with a priori knowledge in order to extrapolate a robust AD disease module. Pharmacophore-based virtual screening was then utilised to build a tailored library of compounds potentially active for AD. RESULTS: In this study, we identified a core disease module for AD, pinpointing known and unknown molecular determinants underlying the skin lesions. We identified skin- and immune-cell type signatures expressed by the disease module, and characterised the impaired cellular functions underlying the complex phenotype of atopic dermatitis. Therefore, by investigating the connectivity of genes belonging to the AD module, we prioritised novel putative biomarkers of the disease. Finally, we defined a tailored compound library by characterising the therapeutic potential of drugs targeting genes within the disease module to facilitate and tailor future drug discovery efforts towards novel pharmacological strategies for AD. CONCLUSIONS: Overall, our study reveals a core disease module providing unprecedented information about genetic, transcriptional and pharmacological relationships that foster drug discovery in atopic dermatitis.

Advancing chemical safety assessment through an omics-based characterization of the test system-chemical interaction.

Assessing chemical safety is essential to evaluate the potential risks of chemical exposure to human health and the environment. Traditional methods relying on animal testing are being replaced by 3R (reduction, refinement, and replacement) principle-based alternatives, mainly depending on in vitro test methods and the Adverse Outcome Pathway framework. However, these approaches often focus on the properties of the compound, missing the broader chemical-biological interaction perspective. Currently, the lack of comprehensive molecular characterization of the in vitro test system results in limited real-world representation and contextualization of the toxicological effect under study. Leveraging omics data strengthens the understanding of the responses of different biological systems, emphasizing holistic chemical-biological interactions when developing in vitro methods. Here, we discuss the relevance of meticulous test system characterization on two safety assessment relevant scenarios and how omics-based, data-driven approaches can improve the future generation of alternative methods.

Expanding adverse outcome pathways towards one health models for nanosafety.

The ever-growing production of nano-enabled products has generated the need for dedicated risk assessment strategies that ensure safety for humans and the environment. Transdisciplinary approaches are needed to support the development of new technologies while respecting environmental limits, as also highlighted by the EU Green Deal Chemicals Strategy for Sustainability and its safe and sustainable by design (SSbD) framework. The One Health concept offers a holistic multiscale approach for the assessment of nanosafety. However, toxicology is not yet capable of explaining the interaction between chemicals and biological systems at the multiscale level and in the context of the One Health framework. Furthermore, there is a disconnect between chemical safety assessment, epidemiology, and other fields of biology that, if unified, would enable the adoption of the One Health model. The development of mechanistic toxicology and the generation of omics data has provided important biological knowledge of the response of individual biological systems to nanomaterials (NMs). On the other hand, epigenetic data have the potential to inform on interspecies mechanisms of adaptation. These data types, however, need to be linked to concepts that support their intuitive interpretation. Adverse Outcome Pathways (AOPs) represent an evolving framework to anchor existing knowledge to chemical risk assessment. In this perspective, we discuss the possibility of integrating multi-level toxicogenomics data, including toxicoepigenetic insights, into the AOP framework. We anticipate that this new direction of toxicogenomics can support the development of One Health models applicable to groups of chemicals and to multiple species in the tree of life.

DREAM: an R package for druggability evaluation of human complex diseases.

MOTIVATION: De novo drug development is a long and expensive process that poses significant challenges from the design to the preclinical testing, making the introduction into the market slow and difficult. This limitation paved the way to the development of drug repurposing, which consists in the re-usage of already approved drugs, developed for other therapeutic indications. Although several efforts have been carried out in the last decade in order to achieve clinically relevant drug repurposing predictions, the amount of repurposed drugs that have been employed in actual pharmacological therapies is still limited. On one hand, mechanistic approaches, including profile-based and network-based methods, exploit the wealth of data about drug sensitivity and perturbational profiles as well as disease transcriptomics profiles. On the other hand, chemocentric approaches, including structure-based methods, take into consideration the intrinsic structural properties of the drugs and their molecular targets. The poor integration between mechanistic and chemocentric approaches is one of the main limiting factors behind the poor translatability of drug repurposing predictions into the clinics. RESULTS: In this work, we introduce DREAM, an R package aimed to integrate mechanistic and chemocentric approaches in a unified computational workflow. DREAM is devoted to the druggability evaluation of pathological conditions of interest, leveraging robust drug repurposing predictions. In addition, the user can derive optimized sets of drugs putatively suitable for combination therapy. In order to show the functionalities of the DREAM package, we report a case study on atopic dermatitis. AVAILABILITY AND IMPLEMENTATION: DREAM is freely available at https://github.com/fhaive/dream. The docker image of DREAM is available at: https://hub.docker.com/r/fhaive/dream.

ESPERANTO: a GLP-field sEmi-SuPERvised toxicogenomics metadAta curatioN TOol.

SUMMARY: Biological data repositories are an invaluable source of publicly available research evidence. Unfortunately, the lack of convergence of the scientific community on a common metadata annotation strategy has resulted in large amounts of data with low FAIRness (Findable, Accessible, Interoperable and Reusable). The possibility of generating high-quality insights from their integration relies on data curation, which is typically an error-prone process while also being expensive in terms of time and human labour. Here, we present ESPERANTO, an innovative framework that enables a standardized semi-supervised harmonization and integration of toxicogenomics metadata and increases their FAIRness in a Good Laboratory Practice-compliant fashion. The harmonization across metadata is guaranteed with the definition of an ad hoc vocabulary. The tool interface is designed to support the user in metadata harmonization in a user-friendly manner, regardless of the background and the type of expertise. AVAILABILITY AND IMPLEMENTATION: ESPERANTO and its user manual are freely available for academic purposes at https://github.com/fhaive/esperanto. The input and the results showcased in Supplementary File S1 are available at the same link.

KNeMAP: a network mapping approach for knowledge-driven comparison of transcriptomic profiles.

MOTIVATION: Transcriptomic data can be used to describe the mechanism of action (MOA) of a chemical compound. However, omics data tend to be complex and prone to noise, making the comparison of different datasets challenging. Often, transcriptomic profiles are compared at the level of individual gene expression values, or sets of differentially expressed genes. Such approaches can suffer from underlying technical and biological variance, such as the biological system exposed on or the machine/method used to measure gene expression data, technical errors and further neglect the relationships between the genes. We propose a network mapping approach for knowledge-driven comparison of transcriptomic profiles (KNeMAP), which combines genes into similarity groups based on multiple levels of prior information, hence adding a higher-level view onto the individual gene view. When comparing KNeMAP with fold change (expression) based and deregulated gene set-based methods, KNeMAP was able to group compounds with higher accuracy with respect to prior information as well as is less prone to noise corrupted data. RESULT: We applied KNeMAP to analyze the Connectivity Map dataset, where the gene expression changes of three cell lines were analyzed after treatment with 676 drugs as well as the Fortino et al. dataset where two cell lines with 31 nanomaterials were analyzed. Although the expression profiles across the biological systems are highly different, KNeMAP was able to identify sets of compounds that induce similar molecular responses when exposed on the same biological system. AVAILABILITY AND IMPLEMENTATION: Relevant data and the KNeMAP function is available at: https://github.com/fhaive/KNeMAP and 10.5281/zenodo.7334711.

Toxicogenomics Data for Chemical Safety Assessment and Development of New Approach Methodologies: An Adverse Outcome Pathway-Based Approach.

Mechanistic toxicology provides a powerful approach to inform on the safety of chemicals and the development of safe-by-design compounds. Although toxicogenomics supports mechanistic evaluation of chemical exposures, its implementation into the regulatory framework is hindered by uncertainties in the analysis and interpretation of such data. The use of mechanistic evidence through the adverse outcome pathway (AOP) concept is promoted for the development of new approach methodologies (NAMs) that can reduce animal experimentation. However, to unleash the full potential of AOPs and build confidence into toxicogenomics, robust associations between AOPs and patterns of molecular alteration need to be established. Systematic curation of molecular events to AOPs will create the much-needed link between toxicogenomics and systemic mechanisms depicted by the AOPs. This, in turn, will introduce novel ways of benefitting from the AOPs, including predictive models and targeted assays, while also reducing the need for multiple testing strategies. Hence, a multi-step strategy to annotate AOPs is developed, and the resulting associations are applied to successfully highlight relevant adverse outcomes for chemical exposures with strong in vitro and in vivo convergence, supporting chemical grouping and other data-driven approaches. Finally, a panel of AOP-derived in vitro biomarkers for pulmonary fibrosis (PF) is identified and experimentally validated.

The integration of large-scale public data and network analysis uncovers molecular characteristics of psoriasis.

In recent years, a growing interest in the characterization of the molecular basis of psoriasis has been observed. However, despite the availability of a large amount of molecular data, many pathogenic mechanisms of psoriasis are still poorly understood. In this study, we performed an integrated analysis of 23 public transcriptomic datasets encompassing both lesional and uninvolved skin samples from psoriasis patients. We defined comprehensive gene co-expression network models of psoriatic lesions and uninvolved skin. Moreover, we curated and exploited a wide range of functional information from multiple public sources in order to systematically annotate the inferred networks. The integrated analysis of transcriptomics data and co-expression networks highlighted genes that are frequently dysregulated and show aberrant patterns of connectivity in the psoriatic lesion compared with the unaffected skin. Our approach allowed us to also identify plausible, previously unknown, actors in the expression of the psoriasis phenotype. Finally, we characterized communities of co-expressed genes associated with relevant molecular functions and expression signatures of specific immune cell types associated with the psoriasis lesion. Overall, integrating experimental driven results with curated functional information from public repositories represents an efficient approach to empower knowledge generation about psoriasis and may be applicable to other complex diseases.

A three component model for superdiffusive motion effectively describes migration of eukaryotic cells moving freely or under a directional stimulus.

Although the simple diffusion model can effectively describe the movement of eukaryotic cells on a culture surface observed at relatively low sampling frequency, at higher sampling rates more complex models are often necessary to better fit the experimental data. Currently available models can describe motion paths by involving additional parameters, such as linearity or directional persistence in time. However sometimes difficulties arise as it is not easy to effectively evaluate persistence in presence of a directional bias. Here we present a procedure which helps solve this problem, based on a model which describes displacement as the vectorial sum of three components: diffusion, persistence and directional bias. The described model has been tested by analysing the migratory behaviour of simulated cell populations and used to analyse a collection of experimental datasets, obtained by observing cell cultures in time lapse microscopy. Overall, the method produces a good description of migration behaviour as it appears to capture the expected increase in the directional bias in presence of wound without a large concomitant increase in the persistence module, allowing it to remain as a physically meaningful quantity in the presence of a directional stimulus.

Factors Modulating COVID-19: A Mechanistic Understanding Based on the Adverse Outcome Pathway Framework.

Addressing factors modulating COVID-19 is crucial since abundant clinical evidence shows that outcomes are markedly heterogeneous between patients. This requires identifying the factors and understanding how they mechanistically influence COVID-19. Here, we describe how eleven selected factors (age, sex, genetic factors, lipid disorders, heart failure, gut dysbiosis, diet, vitamin D deficiency, air pollution and exposure to chemicals) influence COVID-19 by applying the Adverse Outcome Pathway (AOP), which is well-established in regulatory toxicology. This framework aims to model the sequence of events leading to an adverse health outcome. Several linear AOPs depicting pathways from the binding of the virus to ACE2 up to clinical outcomes observed in COVID-19 have been developed and integrated into a network offering a unique overview of the mechanisms underlying the disease. As SARS-CoV-2 infectibility and ACE2 activity are the major starting points and inflammatory response is central in the development of COVID-19, we evaluated how those eleven intrinsic and extrinsic factors modulate those processes impacting clinical outcomes. Applying this AOP-aligned approach enables the identification of current knowledge gaps orientating for further research and allows to propose biomarkers to identify of high-risk patients. This approach also facilitates expertise synergy from different disciplines to address public health issues.

Characterization of ENM Dynamic Dose-Dependent MOA in Lung with Respect to Immune Cells Infiltration.

The molecular effects of exposures to engineered nanomaterials (ENMs) are still largely unknown. In classical inhalation toxicology, cell composition of bronchoalveolar lavage (BAL) is a toxicity indicator at the lung tissue level that can aid in interpreting pulmonary histological changes. Toxicogenomic approaches help characterize the mechanism of action (MOA) of ENMs by investigating the differentially expressed genes (DEG). However, dissecting which molecular mechanisms and events are directly induced by the exposure is not straightforward. It is now generally accepted that direct effects follow a monotonic dose-dependent pattern. Here, we applied an integrated modeling approach to study the MOA of four ENMs by retrieving the DEGs that also show a dynamic dose-dependent profile (dddtMOA). We further combined the information of the dddtMOA with the dose dependency of four immune cell populations derived from BAL counts. The dddtMOA analysis highlighted the specific adaptation pattern to each ENM. Furthermore, it revealed the distinct effect of the ENM physicochemical properties on the induced immune response. Finally, we report three genes dose-dependent in all the exposures and correlated with immune deregulation in the lung. The characterization of dddtMOA for ENM exposures, both for apical endpoints and molecular responses, can further promote toxicogenomic approaches in a regulatory context.

Integrated Network Pharmacology Approach for Drug Combination Discovery: A Multi-Cancer Case Study.

Despite remarkable efforts of computational and predictive pharmacology to improve therapeutic strategies for complex diseases, only in a few cases have the predictions been eventually employed in the clinics. One of the reasons behind this drawback is that current predictive approaches are based only on the integration of molecular perturbation of a certain disease with drug sensitivity signatures, neglecting intrinsic properties of the drugs. Here we integrate mechanistic and chemocentric approaches to drug repositioning by developing an innovative network pharmacology strategy. We developed a multilayer network-based computational framework integrating perturbational signatures of the disease as well as intrinsic characteristics of the drugs, such as their mechanism of action and chemical structure. We present five case studies carried out on public data from The Cancer Genome Atlas, including invasive breast cancer, colon adenocarcinoma, lung squamous cell carcinoma, hepatocellular carcinoma and prostate adenocarcinoma. Our results highlight paclitaxel as a suitable drug for combination therapy for many of the considered cancer types. In addition, several non-cancer-related genes representing unusual drug targets were identified as potential candidates for pharmacological treatment of cancer.

Nextcast: A software suite to analyse and model toxicogenomics data.

The recent advancements in toxicogenomics have led to the availability of large omics data sets, representing the starting point for studying the exposure mechanism of action and identifying candidate biomarkers for toxicity prediction. The current lack of standard methods in data generation and analysis hampers the full exploitation of toxicogenomics-based evidence in regulatory risk assessment. Moreover, the pipelines for the preprocessing and downstream analyses of toxicogenomic data sets can be quite challenging to implement. During the years, we have developed a number of software packages to address specific questions related to multiple steps of toxicogenomics data analysis and modelling. In this review we present the Nextcast software collection and discuss how its individual tools can be combined into efficient pipelines to answer specific biological questions. Nextcast components are of great support to the scientific community for analysing and interpreting large data sets for the toxicity evaluation of compounds in an unbiased, straightforward, and reliable manner. The Nextcast software suite is available at: ( https://github.com/fhaive/nextcast).

Computationally prioritized drugs inhibit SARS-CoV-2 infection and syncytia formation.

The pharmacological arsenal against the COVID-19 pandemic is largely based on generic anti-inflammatory strategies or poorly scalable solutions. Moreover, as the ongoing vaccination campaign is rolling slower than wished, affordable and effective therapeutics are needed. To this end, there is increasing attention toward computational methods for drug repositioning and de novo drug design. Here, multiple data-driven computational approaches are systematically integrated to perform a virtual screening and prioritize candidate drugs for the treatment of COVID-19. From the list of prioritized drugs, a subset of representative candidates to test in human cells is selected. Two compounds, 7-hydroxystaurosporine and bafetinib, show synergistic antiviral effects in vitro and strongly inhibit viral-induced syncytia formation. Moreover, since existing drug repositioning methods provide limited usable information for de novo drug design, the relevant chemical substructures of the identified drugs are extracted to provide a chemical vocabulary that may help to design new effective drugs.

Microarray Data Preprocessing: From Experimental Design to Differential Analysis.

DNA microarray data preprocessing is of utmost importance in the analytical path starting from the experimental design and leading to a reliable biological interpretation. In fact, when all relevant aspects regarding the experimental plan have been considered, the following steps from data quality check to differential analysis will lead to robust, trustworthy results. In this chapter, all the relevant aspects and considerations about microarray preprocessing will be discussed. Preprocessing steps are organized in an orderly manner, from experimental design to quality check and batch effect removal, including the most common visualization methods. Furthermore, we will discuss data representation and differential testing methods with a focus on the most common microarray technologies, such as gene expression and DNA methylation.

VOLTA: adVanced mOLecular neTwork Analysis.

MOTIVATION: Network analysis is a powerful approach to investigate biological systems. It is often applied to study gene co-expression patterns derived from transcriptomics experiments. Even though co-expression analysis is widely used, there is still a lack of tools that are open and customizable on the basis of different network types and analysis scenarios (e.g. through function accessibility), but are also suitable for novice users by providing complete analysis pipelines. RESULTS: We developed VOLTA, a Python package suited for complex co-expression network analysis. VOLTA is designed to allow users direct access to the individual functions, while they are also provided with complete analysis pipelines. Moreover, VOLTA offers when possible multiple algorithms applicable to each analytical step (e.g. multiple community detection or clustering algorithms are provided), hence providing the user with the possibility to perform analysis tailored to their needs. This makes VOLTA highly suitable for experienced users who wish to build their own analysis pipelines for a wide range of networks as well as for novice users for which a 'plug and play' system is provided. AVAILABILITY AND IMPLEMENTATION: The package and used data are available at GitHub: https://github.com/fhaive/VOLTA and 10.5281/zenodo.5171719. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Integrated network analysis reveals new genes suggesting COVID-19 chronic effects and treatment.

The COVID-19 disease led to an unprecedented health emergency, still ongoing worldwide. Given the lack of a vaccine or a clear therapeutic strategy to counteract the infection as well as its secondary effects, there is currently a pressing need to generate new insights into the SARS-CoV-2 induced host response. Biomedical data can help to investigate new aspects of the COVID-19 pathogenesis, but source heterogeneity represents a major drawback and limitation. In this work, we applied data integration methods to develop a Unified Knowledge Space (UKS) and used it to identify a new set of genes associated with SARS-CoV-2 host response, both in vitro and in vivo. Functional analysis of these genes reveals possible long-term systemic effects of the infection, such as vascular remodelling and fibrosis. Finally, we identified a set of potentially relevant drugs targeting proteins involved in multiple steps of the host response to the virus.

Covid-19 acute responses and possible long term consequences: What nanotoxicology can teach us.

Long-term effects of Covid-19 disease are still poorly understood. However, similarities between the responses to SARS-CoV-2 and certain nanomaterials suggest fibrotic pulmonary disease as a concern for public health in the next future. Cross-talk between nanotoxicology and other relevant disciplines can help us to deploy more effective Covid-19 therapies and management strategies.

TinderMIX: Time-dose integrated modelling of toxicogenomics data.

BACKGROUND: Omics technologies have been widely applied in toxicology studies to investigate the effects of different substances on exposed biological systems. A classical toxicogenomic study consists in testing the effects of a compound at different dose levels and different time points. The main challenge consists in identifying the gene alteration patterns that are correlated to doses and time points. The majority of existing methods for toxicogenomics data analysis allow the study of the molecular alteration after the exposure (or treatment) at each time point individually. However, this kind of analysis cannot identify dynamic (time-dependent) events of dose responsiveness. RESULTS: We propose TinderMIX, an approach that simultaneously models the effects of time and dose on the transcriptome to investigate the course of molecular alterations exerted in response to the exposure. Starting from gene log fold-change, TinderMIX fits different integrated time and dose models to each gene, selects the optimal one, and computes its time and dose effect map; then a user-selected threshold is applied to identify the responsive area on each map and verify whether the gene shows a dynamic (time-dependent) and dose-dependent response; eventually, responsive genes are labelled according to the integrated time and dose point of departure. CONCLUSIONS: To showcase the TinderMIX method, we analysed 2 drugs from the Open TG-GATEs dataset, namely, cyclosporin A and thioacetamide. We first identified the dynamic dose-dependent mechanism of action of each drug and compared them. Our analysis highlights that different time- and dose-integrated point of departure recapitulates the toxicity potential of the compounds as well as their dynamic dose-dependent mechanism of action.