Assessment of Drug-Induced Liver Injury through Cell Morphology and Gene Expression Analysis.

Drug-induced liver injury (DILI) is a significant concern in drug development, often leading to drug withdrawal. Although many studies aim to identify biomarkers and gene/pathway signatures related to liver toxicity and aim to predict DILI compounds, this remains a challenge in drug discovery. With a strong development of high-content screening/imaging (HCS/HCI) for phenotypic screening, we explored the morphological cell perturbations induced by DILI compounds. In the first step, cell morphological signatures were associated with two datasets of DILI chemicals (DILIRank and eTox). The mechanisms of action were then analyzed for chemicals having transcriptomics data and sharing similar morphological perturbations. Signaling pathways associated with liver toxicity (cell cycle, cell growth, apoptosis, ...) were then captured, and a hypothetical relation between cell morphological perturbations and gene deregulation was illustrated within our analysis. Finally, using the cell morphological signatures, machine learning approaches were developed to predict chemicals with a potential risk of DILI. Some models showed relevant performance with validation set balanced accuracies between 0.645 and 0.739. Overall, our findings demonstrate the utility of combining HCI with transcriptomics data to identify the morphological and gene expression signatures related to DILI chemicals. Moreover, our protocol could be extended to other toxicity end points, offering a promising avenue for comprehensive toxicity assessment in drug discovery.

Assessing network-based methods in the context of system toxicology.

Introduction: Network-based methods are promising approaches in systems toxicology because they can be used to predict the effects of drugs and chemicals on health, to elucidate the mode of action of compounds, and to identify biomarkers of toxicity. Over the years, the network biology community has developed a wide range of methods, and users are faced with the task of choosing the most appropriate method for their own application. Furthermore, the advantages and limitations of each method are difficult to determine without a proper standard and comparative evaluation of their performance. This study aims to evaluate different network-based methods that can be used to gain biological insight into the mechanisms of drug toxicity, using valproic acid (VPA)-induced liver steatosis as a benchmark. Methods: We provide a comprehensive analysis of the results produced by each method and highlight the fact that the experimental design (how the method is applied) is relevant in addition to the method specifications. We also contribute with a systematic methodology to analyse the results of the methods individually and in a comparative manner. Results: Our results show that the evaluated tools differ in their performance against the benchmark and in their ability to provide novel insights into the mechanism of adverse effects of the drug. We also suggest that aggregation of the results provided by different methods provides a more confident set of candidate genes and processes to further the knowledge of the drug's mechanism of action. Discussion: By providing a detailed and systematic analysis of the results of different network-based tools, we aim to assist users in making informed decisions about the most appropriate method for systems toxicology applications.

Linking chemicals, genes and morphological perturbations to diseases.

The progress in image-based high-content screening technology has facilitated high-throughput phenotypic profiling notably the quantification of cell morphology perturbation by chemicals. However, understanding the mechanism of action of a chemical and linking it to cell morphology and phenotypes remains a challenge in drug discovery. In this study, we intended to integrate molecules that induced transcriptomic perturbations and cellular morphological changes into a biological network in order to assess chemical-phenotypic relationships in humans. Such a network was enriched with existing disease information to suggest molecular and cellular profiles leading to phenotypes. Two datasets were used for this study. Firstly, we used the "Cell Painting morphological profiling assay" dataset, composed of 30,000 compounds tested on human osteosarcoma cells (named U2OS). Secondly, we used the "L1000 mRNA profiling assay" dataset, a collection of transcriptional expression data from cultured human cells treated with approximately 20,000 bioactive small molecules from the Library of Integrated Network-based Cellular Signatures (LINCS). Furthermore, pathways, gene ontology terms and disease enrichments were performed on the transcriptomics data. Overall, our study makes it possible to develop a biological network combining chemical-gene-pathway-morphological perturbation and disease relationships. It contains an ensemble of 9989 chemicals, 732 significant morphological features and 12,328 genes. Through diverse examples, we demonstrated that some drugs shared similar genes, pathways and morphological profiles that, taken together, could help in deciphering chemical-phenotype observations.

Development of a chemogenomics library for phenotypic screening.

With the development of advanced technologies in cell-based phenotypic screening, phenotypic drug discovery (PDD) strategies have re-emerged as promising approaches in the identification and development of novel and safe drugs. However, phenotypic screening does not rely on knowledge of specific drug targets and needs to be combined with chemical biology approaches to identify therapeutic targets and mechanisms of actions induced by drugs and associated with an observable phenotype. In this study, we developed a system pharmacology network integrating drug-target-pathway-disease relationships as well as morphological profile from an existing high content imaging-based high-throughput phenotypic profiling assay known as "Cell Painting". Furthermore, from this network, a chemogenomic library of 5000 small molecules that represent a large and diverse panel of drug targets involved in diverse biological effects and diseases has been developed. Such a platform and a chemogenomic library could assist in the target identification and mechanism deconvolution of some phenotypic assays. The usefulness of the platform is illustrated through examples.

Drug-target-ADR Network and Possible Implications of Structural Variants in Adverse Events.

Adverse drug reactions (ADRs) are of major concern in drug safety. However, due to the biological complexity of human systems, understanding the underlying mechanisms involved in development of ADRs remains a challenging task. Here, we applied network sciences to analyze a tripartite network between 1000 drugs, 1407 targets, and 6164 ADRs. It allowed us to suggest drug targets susceptible to be associated to ADRs and organs, based on the system organ class (SOC). Furthermore, a score was developed to determine the contribution of a set of proteins to ADRs. Finally, we identified proteins that might increase the susceptibility of genes to ADRs, on the basis of knowledge about genomic structural variation in genes encoding proteins targeted by drugs. Such analysis should pave the way to individualize drug therapy and precision medicine.

Impacts of drug resistance mutations on the structural asymmetry of the HIV-2 protease.

BACKGROUND: Drug resistance is a severe problem in HIV treatment. HIV protease is a common target for the design of new drugs for treating HIV infection. Previous studies have shown that the crystallographic structures of the HIV-2 protease (PR2) in bound and unbound forms exhibit structural asymmetry that is important for ligand recognition and binding. Here, we investigated the effects of resistance mutations on the structural asymmetry of PR2. Due to the lack of structural data on PR2 mutants, the 3D structures of 30 PR2 mutants of interest have been modeled using an in silico protocol. Structural asymmetry analysis was carried out with an in-house structural-alphabet-based approach. RESULTS: The systematic comparison of the asymmetry of the wild-type structure and a large number of mutants highlighted crucial residues for PR2 structure and function. In addition, our results revealed structural changes induced by PR2 flexibility or resistance mutations. The analysis of the highlighted structural changes showed that some mutations alter protein stability or inhibitor binding. CONCLUSIONS: This work consists of a structural analysis of the impact of a large number of PR2 resistant mutants based on modeled structures. It suggests three possible resistance mechanisms of PR2, in which structural changes induced by resistance mutations lead to modifications in the dimerization interface, ligand recognition or inhibitor binding.

High Impact: The Role of Promiscuous Binding Sites in Polypharmacology.

The literature focuses on drug promiscuity, which is a drug's ability to bind to several targets, because it plays an essential role in polypharmacology. However, little work has been completed regarding binding site promiscuity, even though its properties are now recognized among the key factors that impact drug promiscuity. Here, we quantified and characterized the promiscuity of druggable binding sites from protein-ligand complexes in the high quality Mother Of All Databases while using statistical methods. Most of the sites (80%) exhibited promiscuity, irrespective of the protein class. Nearly half were highly promiscuous and able to interact with various types of ligands. The corresponding pockets were rather large and hydrophobic, with high sulfur atom and aliphatic residue frequencies, but few side chain atoms. Consequently, their interacting ligands can be large, rigid, and weakly hydrophilic. The selective sites that interacted with one ligand type presented less favorable pocket properties for establishing ligand contacts. Thus, their ligands were highly adaptable, small, and hydrophilic. In the dataset, the promiscuity of the site rather than the drug mainly explains the multiple interactions between the drug and target, as most ligand types are dedicated to one site. This underlines the essential contribution of binding site promiscuity to drug promiscuity between different protein classes.

Statistical Profiling of One Promiscuous Protein Binding Site: Illustrated by Urokinase Catalytic Domain.

While recent literature focuses on drug promiscuity, the characterization of promiscuous binding sites (ability to bind several ligands) remains to be explored. Here, we present a proteochemometric modeling approach to analyze diverse ligands and corresponding multiple binding sub-pockets associated with one promiscuous binding site to characterize protein-ligand recognition. We analyze both geometrical and physicochemical profile correspondences. This approach was applied to examine the well-studied druggable urokinase catalytic domain inhibitor binding site, which results in a large number of complex structures bound to various ligands. This approach emphasizes the importance of jointly characterizing pocket and ligand spaces to explore the impact of ligand diversity on sub-pocket properties and to establish their main profile correspondences. This work supports an interest in mining available 3D holo structures associated with a promiscuous binding site to explore its main protein-ligand recognition tendency.

Cavity Versus Ligand Shape Descriptors: Application to Urokinase Binding Pockets.

We analyzed 78 binding pockets of the human urokinase plasminogen activator (uPA) catalytic domain extracted from a data set of crystallized uPA-ligand complexes. These binding pockets were computed with an original geometric method that does NOT involve any arbitrary parameter, such as cutoff distances, angles, and so on. We measured the deviation from convexity of each pocket shape with the pocket convexity index (PCI). We defined a new pocket descriptor called distributional sphericity coefficient (DISC), which indicates to which extent the protein atoms of a given pocket lie on the surface of a sphere. The DISC values were computed with the freeware PCI. The pocket descriptors and their high correspondences with ligand descriptors are crucial for polypharmacology prediction. We found that the protein heavy atoms lining the urokinases binding pockets are either located on the surface of their convex hull or lie close to this surface. We also found that the radii of the urokinases binding pockets and the radii of their ligands are highly correlated (r = 0.9).