RESUMO
We here introduce the Aquamarine (AQM) dataset, an extensive quantum-mechanical (QM) dataset that contains the structural and electronic information of 59,783 low-and high-energy conformers of 1,653 molecules with a total number of atoms ranging from 2 to 92 (mean: 50.9), and containing up to 54 (mean: 28.2) non-hydrogen atoms. To gain insights into the solvent effects as well as collective dispersion interactions for drug-like molecules, we have performed QM calculations supplemented with a treatment of many-body dispersion (MBD) interactions of structures and properties in the gas phase and implicit water. Thus, AQM contains over 40 global and local physicochemical properties (including ground-state and response properties) per conformer computed at the tightly converged PBE0+MBD level of theory for gas-phase molecules, whereas PBE0+MBD with the modified Poisson-Boltzmann (MPB) model of water was used for solvated molecules. By addressing both molecule-solvent and dispersion interactions, AQM dataset can serve as a challenging benchmark for state-of-the-art machine learning methods for property modeling and de novo generation of large (solvated) molecules with pharmaceutical and biological relevance.
Assuntos
Teoria Quântica , Solventes , Solventes/química , Preparações Farmacêuticas/química , Água/química , Conformação MolecularRESUMO
Federated multipartner machine learning has been touted as an appealing and efficient method to increase the effective training data volume and thereby the predictivity of models, particularly when the generation of training data is resource-intensive. In the landmark MELLODDY project, indeed, each of ten pharmaceutical companies realized aggregated improvements on its own classification or regression models through federated learning. To this end, they leveraged a novel implementation extending multitask learning across partners, on a platform audited for privacy and security. The experiments involved an unprecedented cross-pharma data set of 2.6+ billion confidential experimental activity data points, documenting 21+ million physical small molecules and 40+ thousand assays in on-target and secondary pharmacodynamics and pharmacokinetics. Appropriate complementary metrics were developed to evaluate the predictive performance in the federated setting. In addition to predictive performance increases in labeled space, the results point toward an extended applicability domain in federated learning. Increases in collective training data volume, including by means of auxiliary data resulting from single concentration high-throughput and imaging assays, continued to boost predictive performance, albeit with a saturating return. Markedly higher improvements were observed for the pharmacokinetics and safety panel assay-based task subsets.
Assuntos
Benchmarking , Relação Quantitativa Estrutura-Atividade , Bioensaio , Aprendizado de MáquinaRESUMO
Drug-induced liver injury (DILI) is an important safety concern and a major reason to remove a drug from the market. Advancements in recent machine learning methods have led to a wide range of in silico models for DILI predictive methods based on molecule chemical structures (fingerprints). Existing publicly available DILI data sets used for model building are based on the interpretation of drug labels or patient case reports, resulting in a typical binary clinical DILI annotation. We developed a novel phenotype-based annotation to process hepatotoxicity information extracted from repeated dose in vivo preclinical toxicology studies using INHAND annotation to provide a more informative and reliable data set for machine learning algorithms. This work resulted in a data set of 430 unique compounds covering diverse liver pathology findings which were utilized to develop multiple DILI prediction models trained on the publicly available data (TG-GATEs) using the compound's fingerprint. We demonstrate that the TG-GATEs compounds DILI labels can be predicted well and how the differences between TG-GATEs and the external test compounds (Johnson & Johnson) impact the model generalization performance.
Assuntos
Doença Hepática Induzida por Substâncias e Drogas , Efeitos Colaterais e Reações Adversas Relacionados a Medicamentos , Humanos , Algoritmos , Aprendizado de Máquina , Simulação por ComputadorRESUMO
The search for chemical hit material is a lengthy and increasingly expensive drug discovery process. To improve it, ligand-based quantitative structure-activity relationship models have been broadly applied to optimize primary and secondary compound properties. Although these models can be deployed as early as the stage of molecule design, they have a limited applicability domainâif the structures of interest differ substantially from the chemical space on which the model was trained, a reliable prediction will not be possible. Image-informed ligand-based models partly solve this shortcoming by focusing on the phenotype of a cell caused by small molecules, rather than on their structure. While this enables chemical diversity expansion, it limits the application to compounds physically available and imaged. Here, we employ an active learning approach to capitalize on both of these methods' strengths and boost the model performance of a mitochondrial toxicity assay (Glu/Gal). Specifically, we used a phenotypic Cell Painting screen to build a chemistry-independent model and adopted the results as the main factor in selecting compounds for experimental testing. With the additional Glu/Gal annotation for selected compounds we were able to dramatically improve the chemistry-informed ligand-based model with respect to the increased recognition of compounds from a 10% broader chemical space.
Assuntos
Aprendizado Profundo , Relação Quantitativa Estrutura-Atividade , Ligantes , Descoberta de Drogas/métodosRESUMO
In this paper, we compare the most popular Atom-to-Atom Mapping (AAM) tools: ChemAxon,[1] Indigo,[2] RDTool,[3] NameRXN (NextMove),[4] and RXNMapper[5] which implement different AAM algorithms. An open-source RDTool program was optimized, and its modified version ("new RDTool") was considered together with several consensus mapping strategies. The Condensed Graph of Reaction approach was used to calculate chemical distances and develop the "AAM fixer" algorithm for an automatized correction of erroneous mapping. The benchmarking calculations were performed on a Golden dataset containing 1851 manually mapped and curated reactions. The best performing RXNMapper program together with the AMM Fixer was applied to map the USPTO database. The Golden dataset, mapped USPTO and optimized RDTool are available in the GitHub repository https://github.com/Laboratoire-de-Chemoinformatique.
Assuntos
Benchmarking , Fenômenos Bioquímicos , Algoritmos , Bases de Dados FactuaisRESUMO
Background: Accurate prediction of absorption, distribution, metabolism and excretion (ADME) properties can facilitate the identification of promising drug candidates. Methodology & Results: The authors present the Janssen generic Target Product Profile (gTPP) model, which predicts 18 early ADME properties, employs a graph convolutional neural network algorithm and was trained on between 1000-10,000 internal data points per predicted parameter. gTPP demonstrated stronger predictive power than pretrained commercial ADME models and automatic model builders. Through a novel logging method, the authors report gTPP usage for more than 200 Janssen drug discovery scientists. Conclusion: The investigators successfully enabled the rapid and systematic implementation of predictive ML tools across a drug discovery pipeline in all therapeutic areas. This experience provides useful guidance for other large-scale AI/ML deployment efforts.
Assuntos
Inibidores das Enzimas do Citocromo P-450/farmacologia , Sistema Enzimático do Citocromo P-450/metabolismo , Desenvolvimento de Medicamentos , Inibidores das Enzimas do Citocromo P-450/química , Humanos , Modelos MolecularesRESUMO
The quality of experimental data for chemical reactions is a critical consideration for any reaction-driven study. However, the curation of reaction data has not been extensively discussed in the literature so far. Here, we suggest a 4 steps protocol that includes the curation of individual structures (reactants and products), chemical transformations, reaction conditions and endpoints. Its implementation in Python3 using CGRTools toolkit has been used to clean three popular reaction databases Reaxys, USPTO and Pistachio. The curated USPTO database is available in the GitHub repository (Laboratoire-de-Chemoinformatique/Reaction_Data_Cleaning).
Assuntos
Curadoria de Dados , Bases de Dados Factuais , Padrões de ReferênciaRESUMO
CD38-targeted antibody, daratumumab, is approved for the treatment of multiple myeloma (MM). Phase 1/2 studies GEN501/SIRIUS revealed a novel immunomodulatory mechanism of action (MOA) of daratumumab that enhanced the immune response, reducing natural killer (NK) cells without affecting efficacy or safety. We further evaluated daratumumab's effects on immune cells in whole blood samples of relapsed/refractory MM patients from both treatment arms of the phase 3 POLLUX study (lenalidomide/dexamethasone [Rd] or daratumumab plus Rd [D-Rd]) at baseline (D-Rd, 40; Rd, 45) and after 2 months on treatment (D-Rd, 31; Rd, 33) using cytometry by time-of-flight. We confirmed previous reports of NK cell reduction with D-Rd. Persisting NK cells were phenotypically distinct, with increased expression of HLA-DR, CD69, CD127, and CD27. The proportion of T cells increased preferentially in deep responders to D-Rd, with a higher proportion of CD8+ versus CD4+ T cells. The expansion of CD8+ T cells correlated with clonality, indicating generation of adaptive immune response with D-Rd. D-Rd resulted in a higher proportion of effector memory T cells versus Rd. D-Rd reduced immunosuppressive CD38+ regulatory T cells. This study confirms daratumumab's immunomodulatory MOA in combination with immunomodulatory drugs and provides further insight into immune cell changes and activation status following daratumumab-based therapy.
Assuntos
Protocolos de Quimioterapia Combinada Antineoplásica/uso terapêutico , Biomarcadores/análise , Células Matadoras Naturais/imunologia , Mieloma Múltiplo/imunologia , Linfócitos T Reguladores/imunologia , Linfócitos T/imunologia , Anticorpos Monoclonais/administração & dosagem , Dexametasona/administração & dosagem , Humanos , Células Matadoras Naturais/efeitos dos fármacos , Lenalidomida/administração & dosagem , Mieloma Múltiplo/tratamento farmacológico , Mieloma Múltiplo/patologia , Linfócitos T/efeitos dos fármacos , Linfócitos T Reguladores/efeitos dos fármacosRESUMO
Image-based profiling is a maturing strategy by which the rich information present in biological images is reduced to a multidimensional profile, a collection of extracted image-based features. These profiles can be mined for relevant patterns, revealing unexpected biological activity that is useful for many steps in the drug discovery process. Such applications include identifying disease-associated screenable phenotypes, understanding disease mechanisms and predicting a drug's activity, toxicity or mechanism of action. Several of these applications have been recently validated and have moved into production mode within academia and the pharmaceutical industry. Some of these have yielded disappointing results in practice but are now of renewed interest due to improved machine-learning strategies that better leverage image-based information. Although challenges remain, novel computational technologies such as deep learning and single-cell methods that better capture the biological information in images hold promise for accelerating drug discovery.
Assuntos
Descoberta de Drogas/métodos , Indústria Farmacêutica/métodos , Processamento de Imagem Assistida por Computador/métodos , Aprendizado de Máquina , Animais , Biologia Computacional/métodos , Biologia Computacional/tendências , Descoberta de Drogas/tendências , Indústria Farmacêutica/tendências , Ensaios de Triagem em Larga Escala/métodos , Ensaios de Triagem em Larga Escala/tendências , Humanos , Processamento de Imagem Assistida por Computador/tendências , Aprendizado de Máquina/tendênciasRESUMO
Phenomic profiles are high-dimensional sets of readouts that can comprehensively capture the biological impact of chemical and genetic perturbations in cellular assay systems. Phenomic profiling of compound libraries can be used for compound target identification or mechanism of action (MoA) prediction and other applications in drug discovery. To devise an economical set of phenomic profiling assays, we assembled a library of 1,008 approved drugs and well-characterized tool compounds manually annotated to 218 unique MoAs, and we profiled each compound at four concentrations in live-cell, high-content imaging screens against a panel of 15 reporter cell lines, which expressed a diverse set of fluorescent organelle and pathway markers in three distinct cell lineages. For 41 of 83 testable MoAs, phenomic profiles accurately ranked the reference compounds (AUC-ROC ≥ 0.9). MoAs could be better resolved by screening compounds at multiple concentrations than by including replicates at a single concentration. Screening additional cell lineages and fluorescent markers increased the number of distinguishable MoAs but this effect quickly plateaued. There remains a substantial number of MoAs that were hard to distinguish from others under the current study's conditions. We discuss ways to close this gap, which will inform the design of future phenomic profiling efforts.
Assuntos
Produtos Biológicos/farmacologia , Proteínas Luminescentes/genética , Fenômica/métodos , Bibliotecas de Moléculas Pequenas/farmacologia , Células A549 , Linhagem Celular , Descoberta de Drogas , Regulação da Expressão Gênica/efeitos dos fármacos , Células Hep G2 , Humanos , Proteínas Luminescentes/metabolismoRESUMO
Artificial intelligence (AI) is undergoing a revolution thanks to the breakthroughs of machine learning algorithms in computer vision, speech recognition, natural language processing and generative modelling. Recent works on publicly available pharmaceutical data showed that AI methods are highly promising for Drug Target prediction. However, the quality of public data might be different than that of industry data due to different labs reporting measurements, different measurement techniques, fewer samples and less diverse and specialized assays. As part of a European funded project (ExCAPE), that brought together expertise from pharmaceutical industry, machine learning, and high-performance computing, we investigated how well machine learning models obtained from public data can be transferred to internal pharmaceutical industry data. Our results show that machine learning models trained on public data can indeed maintain their predictive power to a large degree when applied to industry data. Moreover, we observed that deep learning derived machine learning models outperformed comparable models, which were trained by other machine learning algorithms, when applied to internal pharmaceutical company datasets. To our knowledge, this is the first large-scale study evaluating the potential of machine learning and especially deep learning directly at the level of industry-scale settings and moreover investigating the transferability of publicly learned target prediction models towards industrial bioactivity prediction pipelines.
RESUMO
Daratumumab is a CD38-targeted human monoclonal antibody with direct anti-myeloma cell mechanisms of action. Flow cytometry in relapsed and/or refractory multiple myeloma (RRMM) patients treated with daratumumab revealed cytotoxic T-cell expansion and reduction of immune-suppressive populations, suggesting immune modulation as an additional mechanism of action. Here, we performed an in-depth analysis of the effects of daratumumab on immune-cell subpopulations using high-dimensional mass cytometry. Whole-blood and bone-marrow baseline and on-treatment samples from RRMM patients who participated in daratumumab monotherapy studies (SIRIUS and GEN501) were evaluated with high-throughput immunophenotyping. In daratumumab-treated patients, the intensity of CD38 marker expression decreased on many immune cells in SIRIUS whole-blood samples. Natural killer (NK) cells were depleted with daratumumab, with remaining NK cells showing increased CD69 and CD127, decreased CD45RA, and trends for increased CD25, CD27, and CD137 and decreased granzyme B. Immune-suppressive population depletion paralleled previous findings, and a newly observed reduction in CD38+ basophils was seen in patients who received monotherapy. After 2 months of daratumumab, the T-cell population in whole-blood samples from responders shifted to a CD8 prevalence with higher granzyme B positivity (P = 0.017), suggesting increased killing capacity and supporting monotherapy-induced CD8+ T-cell activation. High-throughput cytometry immune profiling confirms and builds upon previous flow cytometry data, including comparable CD38 marker intensity on plasma cells, NK cells, monocytes, and B/T cells. Interestingly, a shift toward cytolytic granzyme B+ T cells was also observed and supports adaptive responses in patients that may contribute to depth of response. © 2018 The Authors. Cytometry Part A published by Wiley Periodicals, Inc. on behalf of International Society for Advancement of Cytometry.
Assuntos
ADP-Ribosil Ciclase 1/imunologia , Anticorpos Monoclonais/uso terapêutico , Antineoplásicos/uso terapêutico , Células Matadoras Naturais/efeitos dos fármacos , Células Matadoras Naturais/imunologia , Mieloma Múltiplo/tratamento farmacológico , Mieloma Múltiplo/imunologia , Antígenos de Diferenciação de Linfócitos T/metabolismo , Basófilos/citologia , Basófilos/efeitos dos fármacos , Basófilos/imunologia , Células da Medula Óssea/citologia , Células da Medula Óssea/imunologia , Linfócitos T CD4-Positivos/citologia , Linfócitos T CD4-Positivos/imunologia , Linfócitos T CD8-Positivos/citologia , Linfócitos T CD8-Positivos/imunologia , Citometria de Fluxo , Granzimas/metabolismo , Humanos , Imunofenotipagem , Células Matadoras Naturais/citologia , Mieloma Múltiplo/sangue , Mieloma Múltiplo/metabolismo , RecidivaRESUMO
Deep learning is currently the most successful machine learning technique in a wide range of application areas and has recently been applied successfully in drug discovery research to predict potential drug targets and to screen for active molecules. However, due to (1) the lack of large-scale studies, (2) the compound series bias that is characteristic of drug discovery datasets and (3) the hyperparameter selection bias that comes with the high number of potential deep learning architectures, it remains unclear whether deep learning can indeed outperform existing computational methods in drug discovery tasks. We therefore assessed the performance of several deep learning methods on a large-scale drug discovery dataset and compared the results with those of other machine learning and target prediction methods. To avoid potential biases from hyperparameter selection or compound series, we used a nested cluster-cross-validation strategy. We found (1) that deep learning methods significantly outperform all competing methods and (2) that the predictive performance of deep learning is in many cases comparable to that of tests performed in wet labs (i.e., in vitro assays).
RESUMO
By adding biological information, beyond the chemical properties and desired effect of a compound, uncharted compound areas and connections can be explored. In this study, we add transcriptional information for 31K compounds of Janssen's primary screening deck, using the HT L1000 platform and assess (a) the transcriptional connection score for generating compound similarities, (b) machine learning algorithms for generating target activity predictions, and (c) the scaffold hopping potential of the resulting hits. We demonstrate that the transcriptional connection score is best computed from the significant genes only and should be interpreted within its confidence interval for which we provide the stats. These guidelines help to reduce noise, increase reproducibility, and enable the separation of specific and promiscuous compounds. The added value of machine learning is demonstrated for the NR3C1 and HSP90 targets. Support Vector Machine models yielded balanced accuracy values ≥80% when the expression values from DDIT4 & SERPINE1 and TMEM97 & SPR were used to predict the NR3C1 and HSP90 activity, respectively. Combining both models resulted in 22 new and confirmed HSP90-independent NR3C1 inhibitors, providing two scaffolds (i.e., pyrimidine and pyrazolo-pyrimidine), which could potentially be of interest in the treatment of depression (i.e., inhibiting the glucocorticoid receptor (i.e., NR3C1), while leaving its chaperone, HSP90, unaffected). As such, the initial hit rate increased by a factor 300, as less, but more specific chemistry could be screened, based on the upfront computed activity predictions.
Assuntos
Proteínas de Choque Térmico HSP90/genética , Ensaios de Triagem em Larga Escala , Pirazóis/farmacologia , Pirimidinas/farmacologia , Receptores de Glucocorticoides/genética , Transcriptoma , Proteínas de Choque Térmico HSP90/metabolismo , Humanos , Receptores de Glucocorticoides/metabolismo , Máquina de Vetores de SuporteRESUMO
In both academia and the pharmaceutical industry, large-scale assays for drug discovery are expensive and often impractical, particularly for the increasingly important physiologically relevant model systems that require primary cells, organoids, whole organisms, or expensive or rare reagents. We hypothesized that data from a single high-throughput imaging assay can be repurposed to predict the biological activity of compounds in other assays, even those targeting alternate pathways or biological processes. Indeed, quantitative information extracted from a three-channel microscopy-based screen for glucocorticoid receptor translocation was able to predict assay-specific biological activity in two ongoing drug discovery projects. In these projects, repurposing increased hit rates by 50- to 250-fold over that of the initial project assays while increasing the chemical structure diversity of the hits. Our results suggest that data from high-content screens are a rich source of information that can be used to predict and replace customized biological assays.
Assuntos
Reposicionamento de Medicamentos/métodos , Processamento de Imagem Assistida por Computador/métodos , Aprendizado de Máquina , Redes Neurais de Computação , Antineoplásicos/farmacologia , Linhagem Celular Tumoral , Ensaios de Triagem em Larga Escala/métodos , Humanos , Neoplasias/tratamento farmacológicoRESUMO
Clinical efficacy regularly requires the combination of drugs. For an early estimation of the clinical value of (potentially many) combinations of pharmacologic compounds during discovery, the observed combination effect is typically compared to that expected under a null model. Mechanistic accuracy of that null model is not aspired to; to the contrary, combinations that deviate favorably from the model (and thereby disprove its accuracy) are prioritized. Arguably the most popular null model is the Loewe Additivity model, which conceptually maps any assay under study to a (virtual) single-step enzymatic reaction. It is easy-to-interpret and requires no other information than the concentration-response curves of the individual compounds. However, the original Loewe model cannot accommodate concentration-response curves with different maximal responses and, by consequence, combinations of an agonist with a partial or inverse agonist. We propose an extension, named Biochemically Intuitive Generalized Loewe (BIGL), that can address different maximal responses, while preserving the biochemical underpinning and interpretability of the original Loewe model. In addition, we formulate statistical tests for detecting synergy and antagonism, which allow for detecting statistically significant greater/lesser observed combined effects than expected from the null model. Finally, we demonstrate the novel method through application to several publicly available datasets.
Assuntos
Agonismo de Drogas , Antagonismo de Drogas , Quimioterapia Combinada , Modelos Teóricos , Relação Dose-Resposta a Droga , Humanos , Resultado do TratamentoRESUMO
The identification of functional non-coding mutations is a key challenge in the field of genomics. Here we introduce µ-cisTarget to filter, annotate and prioritize cis-regulatory mutations based on their putative effect on the underlying "personal" gene regulatory network. We validated µ-cisTarget by re-analyzing the TAL1 and LMO1 enhancer mutations in T-ALL, and the TERT promoter mutation in melanoma. Next, we re-sequenced the full genomes of ten cancer cell lines and used matched transcriptome data and motif discovery to identify master regulators with de novo binding sites that result in the up-regulation of nearby oncogenic drivers. µ-cisTarget is available from http://mucistarget.aertslab.org .
Assuntos
Análise Mutacional de DNA/métodos , Redes Reguladoras de Genes , Genes Neoplásicos , Mutação , Neoplasias/genética , Sequências Reguladoras de Ácido Nucleico , Algoritmos , Sítios de Ligação , Linhagem Celular Tumoral , Feminino , Perfilação da Expressão Gênica , Genômica/métodos , Humanos , Masculino , Neoplasias/metabolismo , Medicina de Precisão/métodos , Fatores de Transcrição/metabolismoRESUMO
Chemogenomics data generally refers to the activity data of chemical compounds on an array of protein targets and represents an important source of information for building in silico target prediction models. The increasing volume of chemogenomics data offers exciting opportunities to build models based on Big Data. Preparing a high quality data set is a vital step in realizing this goal and this work aims to compile such a comprehensive chemogenomics dataset. This dataset comprises over 70 million SAR data points from publicly available databases (PubChem and ChEMBL) including structure, target information and activity annotations. Our aspiration is to create a useful chemogenomics resource reflecting industry-scale data not only for building predictive models of in silico polypharmacology and off-target effects but also for the validation of cheminformatics approaches in general.
RESUMO
BACKGROUND & AIMS: Simeprevir is an oral hepatitis C virus (HCV) NS3/4A protease inhibitor approved for treatment of chronic HCV infection. Baseline NS3 polymorphisms in all patients and emerging mutations in patients who failed to achieve sustained virologic response (SVR) with simeprevir plus peginterferon/ribavirin (PR) in Phase IIb/III studies are described. METHODS: Baseline sequencing data were available for 2007 genotype 1 (GT1)-infected patients. Post-baseline data were available for 197/245 simeprevir-treated patients who did not achieve SVR. In vitro simeprevir susceptibility was assessed in a transient replicon assay as site-directed mutants or in chimeric replicons with patient-derived NS3 protease sequences. RESULTS: Baseline NS3 polymorphisms at positions associated with reduced in vitro susceptibility to simeprevir (43, 80, 122, 155, 156, and/or 168; EC50 fold change >2.0) were uncommon (1.3% [26/2007]), with the exception of Q80K, which confers â¼10-fold reduction in simeprevir activity in vitro (13.7% [274/2007]; GT1a 29.5% [269/911], GT1b 0.5% [5/1096]). Baseline Q80K had minor effect on initial response to simeprevir/PR, but resulted in lower SVR rates. Overall, 91.4% of simeprevir-treated patients [180/197] without SVR had emerging mutations at NS3 positions 80, 122, 155, and/or 168 at failure (mainly R155K in GT1a with and without Q80K, and D168V in GT1b), conferring high-level resistance in vitro (EC50 fold change >50). Emerging mutations were no longer detectable by population sequencing at study end in 50% [90/180] of patients (median follow-up 28.4weeks). CONCLUSIONS: Simeprevir treatment failure was usually associated with emerging high-level resistance mutations, which became undetectable over time in half of the patients.