Design, quality and validation of the EU-OPENSCREEN fragment library poised to a high-throughput screening collection

The EU-OPENSCREEN (EU-OS) European Research Infrastructure Consortium (ERIC) is a multinational, not-for-profit initiative that integrates high-capacity screening platforms and chemistry groups across Europe to facilitate research in chemical biology and early drug discovery. Over the years, the EU-OS has assembled a high-throughput screening compound collection, the European Chemical Biology Library (ECBL), that contains approximately 100 000 commercially available small molecules and a growing number of thousands of academic compounds crowdsourced through our network of European and non-European chemists. As an extension of the ECBL, here we describe the computational design, quality control and use case screenings of the European Fragment Screening Library (EFSL) composed of 1056 mini and small chemical fragments selected from a substructure analysis of the ECBL. Access to the EFSL is open to researchers from both academia and industry. Using EFSL, eight fragment screening campaigns using different structural and biophysical methods have successfully identified fragment hits in the last two years. As one of the highlighted projects for antibiotics, we describe the screening by Bio-Layer Interferometry (BLI) of the EFSL, the identification of a 35 μM fragment hit targeting the beta-ketoacyl-ACP synthase 2 (FabF), its binding confirmation to the protein by X-ray crystallography (PDB 8PJ0), its subsequent rapid exploration of its surrounding chemical space through hit-picking of ECBL compounds that contain the fragment hit as a core substructure, and the final binding confirmation of two follow-up hits by X-ray crystallography (PDB 8R0I and 8R1V).

CIPSI: An open chemical intellectual property service for medicinal chemists

The availability of patent chemical data offers public access to a chemical space that is not well covered by other sources collecting small molecules from scholarly literature. However, open applications to facilitate the search and analysis of biologically-relevant molecular structures present in patents are still largely missing. We have developed CIPSI, an open Chemical Intellectual Property Service @ IMIM to assist medicinal chemists in searching and analysing molecules in SureChEMBL patents. The current version contains 6,240,500 molecules from 236,689 pharmacological patents, of which 5,949,214 are confidently assigned to core chemical structures reminiscent of the Markush structure in the patent claim. The platform includes some graphical tools to facilitate comparative patent analyses between drugs, chemical substructures, and company assignees. CIPSI is available at https://cipsi.org.

In silico approaches in carcinogenicity hazard assessment: case study of pregabalin, a nongenotoxic mouse carcinogen

In silico toxicology protocols are meant to support computationally-based assessments using principles that ensure that results can be generated, recorded, communicated, archived, and then evaluated in a uniform, consistent, and reproducible manner. We investigated the availability of in silico models to predict the carcinogenic potential of pregabalin using the ten key characteristics of carcinogens as a framework for organizing mechanistic studies. Pregabalin is a single-species carcinogen producing only one type of tumor, hemangiosarcomas in mice via a nongenotoxic mechanism. The overall goal of this exercise is to test the ability of in silico models to predict nongenotoxic carcinogenicity with pregabalin as a case study. The established mode of action (MOA) of pregabalin is triggered by tissue hypoxia, leading to oxidative stress (KC5), chronic inflammation (KC6), and increased cell proliferation (KC10) of endothelial cells. Of these KCs, in silico models are available only for selected endpoints in KC5, limiting the usefulness of computational tools in prediction of pregabalin carcinogenicity. KC1 (electrophilicity), KC2 (genotoxicity), and KC8 (receptor-mediated effects), for which predictive in silico models exist, do not play a role in this mode of action. Confidence in the overall assessments is considered to be medium to high for KCs 1, 2, 5, 6, 7 (immune system effects), 8, and 10 (cell proliferation), largely due to the high-quality experimental data. In order to move away from dependence on animal data, development of reliable in silico models for prediction of oxidative stress, chronic inflammation, immunosuppression, and cell proliferation will be critical for the ability to predict nongenotoxic compound carcinogenicity.

Impact of the COVID-19 pandemic on the spontaneous reporting and signal detection of adverse drug events

External factors severely affecting in a short period of time the spontaneous reporting of adverse events (AEs) can significantly impact drug safety signal detection. Coronavirus disease 2019 (COVID-19) represented an enormous challenge for health systems, with over 767 million cases and massive vaccination campaigns involving over 70% of the worldwide population. This study investigates the potential masking effect on certain AEs caused by the substantial increase in reports solely related to COVID-19 vaccines within various spontaneous reporting systems (SRSs). Three SRSs were used to monitor AEs reporting before and during the pandemic, namely, the World Health Organisation (WHO) global individual case safety reports database (VigiBase®), the United States Food and Drug Administration Adverse Event Reporting System (FAERS) and the Japanese Adverse Drug Event Report database (JADER). Findings revealed a sudden over-reporting of 35 AEs (≥ 200%) during the pandemic, with an increment of the RRF value in 2021 of at least double the RRF reported in 2020. This translates into a substantial reduction in signals of disproportionate reporting (SDR) due to the massive inclusion of COVID-19 vaccine reports. To mitigate the masking effect of COVID-19 vaccines in post-marketing SRS analyses, we recommend utilizing COVID-19-corrected versions for a more accurate assessment.

A multiscale orchestrated computational framework to reveal emergent phenomena in neuroblastoma

Neuroblastoma is a complex and aggressive type of cancer that affects children. Current treatments involve a combination of surgery, chemotherapy, radiotherapy, and stem cell transplantation. However, treatment outcomes vary due to the heterogeneous nature of the disease. Computational models have been used to analyse data, simulate biological processes, and predict disease progression and treatment outcomes. While continuum cancer models capture the overall behaviour of tumours, and agent-based models represent the complex behaviour of individual cells, multiscale models represent interactions at different organisational levels, providing a more comprehensive understanding of the system. In 2018, the PRIMAGE consortium was formed to build a cloud-based decision support system for neuroblastoma, including a multi-scale model for patient-specific simulations of disease progression. In this work we have developed this multi-scale model that includes data such as patient's tumour geometry, cellularity, vascularization, genetics and type of chemotherapy treatment, and integrated it into an online platform that runs the simulations on a high-performance computation cluster using Onedata and Kubernetes technologies. This infrastructure will allow clinicians to optimise treatment regimens and reduce the number of costly and time-consuming clinical trials. This manuscript outlines the challenging framework's model architecture, data workflow, hypothesis, and resources employed in its development.

Evaluation of QSAR models for predicting mutagenicity: outcome of the Second Ames/QSAR international challenge project

Quantitative structure-activity relationship (QSAR) models are powerful in silico tools for predicting the mutagenicity of unstable compounds, impurities and metabolites that are difficult to examine using the Ames test. Ideally, Ames/QSAR models for regulatory use should demonstrate high sensitivity, low false-negative rate and wide coverage of chemical space. To promote superior model development, the Division of Genetics and Mutagenesis, National Institute of Health Sciences, Japan (DGM/NIHS), conducted the Second Ames/QSAR International Challenge Project (2020-2022) as a successor to the First Project (2014-2017), with 21 teams from 11 countries participating. The DGM/NIHS provided a curated training dataset of approximately 12,000 chemicals and a trial dataset of approximately 1,600 chemicals, and each participating team predicted the Ames mutagenicity of each trial chemical using various Ames/QSAR models. The DGM/NIHS then provided the Ames test results for trial chemicals to assist in model improvement. Although overall model performance on the Second Project was not superior to that on the First, models from the eight teams participating in both projects achieved higher sensitivity than models from teams participating in only the Second Project. Thus, these evaluations have facilitated the development of QSAR models.

Illuminating the Chemical Space of Untargeted Proteins

According to the Illuminating the Druggable Genome (IDG) initiative, 90% of the proteins encoded by the human genome still lack an identified active ligand, that is, a small molecule with biologically relevant binding potency or functional activity in an in vitro assay. Under this scenario, there is an urgent need for new approaches to chemically address these yet untargeted proteins. It is widely recognized that the best starting point for generating novel small molecules for proteins is to exploit the expected polypharmacology of known active ligands across phylogenetically related proteins following the paradigm that similar proteins are likely to interact with similar ligands. Here, we introduce a computational strategy to identify privileged structures that, when chemically expanded, are highly probable to contain active small molecules for untargeted proteins. The protocol was first tested on a set of 576 currently targeted proteins having at least one protein family sibling the year before their first active ligand was reported. A privileged structure contained in active ligands that were identified in the following years was correctly anticipated for 214 (37%) of those targeted proteins, a lower-bound recall estimate when considering data completeness issues. When applied to a set of 1184 untargeted potential druggable genes in cancer, the identification of privileged structures from known bioactive ligands of protein family siblings allowed for extracting a priority list of diverse commercially available small molecules for 960 of them. Assuming a minimum success rate of 37%, the chemical library selections should be able to deliver active ligands for at least 355 currently untargeted proteins associated with cancer.

Primary and Community Care Transformation in Post-COVID Era: Nationwide General Practitioner Survey

Introduction: The health emergency caused by COVID-19 has led to substantial changes in the usual working system of primary healthcare centers and in relations with users. The Catalan Society of Family and Community Medicine designed a survey that aimed to collect the opinions and facilitate the participation of its partners on what the future work model of general practitioners (GPs) should look like post-COVID-19. Methodology: Online survey of Family and Community Medicine members consisting of filiation data, 22 Likert-type multiple-choice questions grouped in five thematic axes, and a free text question. Results: The number of respondents to the questionnaire was 1051 (22.6% of all members): 83.2% said they spent excessive time on bureaucratic tasks; 91.8% were against call center systems; 66% believed that home care is the responsibility of every family doctor; 77.5% supported continuity of care as a fundamental value of patient-centered care; and >90% defended the contracting of complementary tests and first hospital visits from primary healthcare (PHC). Conclusions: The survey responses describe a strong consensus on the identity and competencies of the GP and on the needs of and the threats to the PHC system. The demand for an increase in health resources, greater professional leadership, elimination of bureaucracy, an increase in the number of health professionals, and greater management autonomy, are the axes towards which a new era in PHC should be directed.

Chemoinformatics and artificial intelligence colloquium: progress and challenges in developing bioactive compounds

We report the main conclusions of the first Chemoinformatics and Artificial Intelligence Colloquium, Mexico City, June 15-17, 2022. Fifteen lectures were presented during a virtual public event with speakers from industry, academia, and non-for-profit organizations. Twelve hundred and ninety students and academics from more than 60 countries. During the meeting, applications, challenges, and opportunities in drug discovery, de novo drug design, ADME-Tox (absorption, distribution, metabolism, excretion and toxicity) property predictions, organic chemistry, peptides, and antibiotic resistance were discussed. The program along with the recordings of all sessions are freely available at https://www.difacquim.com/english/events/2022-colloquium/ .

Fragment-based virtual screening identifies a first-in-class preclinical drug candidate for Huntington's disease

Currently, there are no therapies available to modify the disease progression of Huntington's disease (HD). Recent clinical trial failures of antisense oligonucleotide candidates in HD have demonstrated the need for new therapeutic approaches. Here, we developed a novel in-silico fragment scanning approach across the surface of mutant huntingtin (mHTT) polyQ and predicted four hit compounds. Two rounds of compound analoging using a strategy of testing structurally similar compounds in an affinity assay rapidly identified GLYN122. In vitro, GLYN122 directly binds and reduces mHTT and induces autophagy in neurons. In vivo, our results confirm that GLYN122 can reduce mHTT in the cortex and striatum of the R/2 mouse model of Huntington's disease and subsequently improve motor symptoms. Thus, the in-vivo pharmacology profile of GLYN122 is a potential new preclinical candidate for the treatment of HD.

Environmental levels of carbaryl impair zebrafish larvae behaviour: The potential role of ADRA2B and HTR2B

The insecticide carbaryl is commonly found in indirectly exposed freshwater ecosystems at low concentrations considered safe for fish communities. In this study, we showed that after only 24 h of exposure to environmental concentrations of carbaryl (0.066-660 ng/L), zebrafish larvae exhibit impairments in essential behaviours. Interestingly, the observed behavioural effects induced by carbaryl were acetylcholinesterase-independent. To elucidate the molecular initiating event that resulted in the observed behavioural effects, in silico predictions were followed by in vitro validation. We identified two target proteins that potentially interacted with carbaryl, the α2B adrenoceptor (ADRA2B) and the serotonin 2B receptor (HTR2B). Using a pharmacological approach, we then tested the hypothesis that carbaryl had antagonistic interactions with both receptors. Similar to yohimbine and SB204741, which are prototypic antagonists of ADRA2B and HTR2B, respectively, carbaryl increased the heart rate of zebrafish larvae. When we compared the behavioural effects of a 24-h exposure to these pharmacological antagonists with those of carbaryl, a high degree of similarity was found. These results strongly suggest that antagonism of both ADRA2B and HTR2B is the molecular initiating event that leads to adverse outcomes in zebrafish larvae that have undergone 24 h of exposure to environmentally relevant levels of carbaryl.

A Workflow of Integrated Resources to Catalyze Network Pharmacology Driven COVID-19 Research

In the event of an outbreak due to an emerging pathogen, time is of the essence to contain or to mitigate the spread of the disease. Drug repositioning is one of the strategies that has the potential to deliver therapeutics relatively quickly. The SARS-CoV-2 pandemic has shown that integrating critical data resources to drive drug-repositioning studies, involving host-host, host-pathogen, and drug-target interactions, remains a time-consuming effort that translates to a delay in the development and delivery of a life-saving therapy. Here, we describe a workflow we designed for a semiautomated integration of rapidly emerging data sets that can be generally adopted in a broad network pharmacology research setting. The workflow was used to construct a COVID-19 focused multimodal network that integrates 487 host-pathogen, 63 278 host-host protein, and 1221 drug-target interactions. The resultant Neo4j graph database named "Neo4COVID19" is made publicly accessible via a web interface and via API calls based on the Bolt protocol. Details for accessing the database are provided on a landing page (https://neo4covid19.ncats.io/). We believe that our Neo4COVID19 database will be a valuable asset to the research community and will catalyze the discovery of therapeutics to fight COVID-19.

Rare Variation in Drug Metabolism and Long QT Genes and the Genetic Susceptibility to Acquired Long QT Syndrome

BACKGROUND

Acquired long QT syndrome (aLQTS) is a serious unpredictable adverse drug reaction. Pharmacogenomic markers may predict risk.

METHODS

Among 153 aLQTS patients (mean age 58 years [range, 14-88], 98.7% White, 85.6% symptomatic), computational methods identified proteins interacting most significantly with 216 QT-prolonging drugs. All cases underwent sequencing of 31 candidate genes arising from this analysis or associating with congenital LQTS. Variants were filtered using a minor allele frequency <1% and classified for susceptibility for aLQTS. Gene-burden analyses were then performed comparing the primary cohort to control exomes (n=452) and an independent replication aLQTS exome sequencing cohort.

RESULTS

In 25.5% of cases, at least one rare variant was identified: 22.2% of cases carried a rare variant in a gene associated with congenital LQTS, and in 4% of cases that variant was known to be pathogenic or likely pathogenic for congenital LQTS; 7.8% cases carried a cytochrome-P450 (CYP) gene variant. Of 12 identified CYP variants, 11 (92%) were in an enzyme known to metabolize at least one culprit drug to which the subject had been exposed. Drug-drug interactions that affected culprit drug metabolism were found in 19% of cases. More than one congenital LQTS variant, CYP gene variant, or drug interaction was present in 7.8% of cases. Gene-burden analyses of the primary cohort compared to control exomes (n=452), and an independent replication aLQTS exome sequencing cohort (n=67) and drug-tolerant controls (n=148) demonstrated an increased burden of rare (minor allele frequency<0.01) variants in CYP genes but not LQTS genes.

CONCLUSIONS

Rare susceptibility variants in CYP genes are emerging as potentially important pharmacogenomic risk markers for aLQTS and could form part of personalized medicine approaches in the future.

In Silico Approaches In Carcinogenicity Hazard Assessment: Current Status and Future Needs

Historically, identifying carcinogens has relied primarily on tumor studies in rodents, which require enormous resources in both money and time. In silico models have been developed for predicting rodent carcinogens but have not yet found general regulatory acceptance, in part due to the lack of a generally accepted protocol for performing such an assessment as well as limitations in predictive performance and scope. There remains a need for additional, improved in silico carcinogenicity models, especially ones that are more human-relevant, for use in research and regulatory decision-making. As part of an international effort to develop in silico toxicological protocols, a consortium of toxicologists, computational scientists, and regulatory scientists across several industries and governmental agencies evaluated the extent to which in silico models exist for each of the recently defined 10 key characteristics (KCs) of carcinogens. This position paper summarizes the current status of in silico tools for the assessment of each KC and identifies the data gaps that need to be addressed before a comprehensive in silico carcinogenicity protocol can be developed for regulatory use.

Extended connectivity interaction features: improving binding affinity prediction through chemical description

MOTIVATION

Machine-learning scoring functions (SFs) have been found to outperform standard SFs for binding affinity prediction of protein-ligand complexes. A plethora of reports focus on the implementation of increasingly complex algorithms, while the chemical description of the system has not been fully exploited.

RESULTS

Herein, we introduce Extended Connectivity Interaction Features (ECIF) to describe protein-ligand complexes and build machine-learning SFs with improved predictions of binding affinity. ECIF are a set of protein-ligand atom-type pair counts that take into account each atom's connectivity to describe it and thus define the pair types. ECIF were used to build different machine-learning models to predict protein-ligand affinities (pKd/pKi). The models were evaluated in terms of 'scoring power' on the Comparative Assessment of Scoring Functions 2016. The best models built on ECIF achieved Pearson correlation coefficients of 0.857 when used on its own, and 0.866 when used in combination with ligand descriptors, demonstrating ECIF descriptive power.

AVAILABILITY AND IMPLEMENTATION

Data and code to reproduce all the results are freely available at https://github.com/DIFACQUIM/ECIF.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Androgenic activation, impairment of the monoaminergic system and altered behavior in zebrafish larvae exposed to environmental concentrations of fenitrothion

Fenitrothion is an organophosphorus insecticide usually found in aquatic ecosystems at concentrations in the range of low ng/L. In this manuscript we show that 24 h exposure to environmental concentrations of fenitrothion, from ng/L to low μg/L, altered basal locomotor activity, visual-motor response and acoustic/vibrational escape response of zebrafish larvae. Furthermore, fenitrothion and expression of gap43a, gfap, atp2b1a, and mbp exhibited a significant non-monotonic concentration-response relationship. Once determined that environmental concentrations of fenitrothion were neurotoxic for zebrafish larvae, a computational analysis identified potential protein targets of this compound. Some of the predictions, including interactions with acetylcholinesterase, monoamine-oxidases and androgen receptor (AR), were experimentally validated. Binding to AR was the most suitable candidate for molecular initiating event, as indicated by both the up-regulation of cyp19a1b and sult2st3 and the non-monotonic relationship found between fenitrothion and the observed responses. Finally, when the integrity of the monoaminergic system was evaluated, altered levels of L-DOPA, DOPAC, HVA and 5-HIAA were found, as well as a significant up-regulation of slc18a2 expression at the lowest concentrations of fenitrothion. These data strongly suggest that concentrations of fenitrothion commonly found in aquatic ecosystems present a significant environmental risk for fish communities.

Identification of the Core Chemical Structure in SureChEMBL Patents

The SureChEMBL database provides open access to 17 million chemical entities mentioned in 14 million patents published since 1970. However, alongside with molecules covered by patent claims, the database is full of starting materials and intermediate products of little pharmacological relevance. Herein, we introduce a new filtering protocol to automatically select the core chemical structures best representing a congeneric series of pharmacologically relevant molecules in patents. The protocol is first validated against a selection of 890 SureChEMBL patents for which a total of 51,738 manually curated molecules are deposited in ChEMBL. Our protocol was able to select 92.5% of the molecules in ChEMBL from all 270,968 molecules in SureChEMBL for those patents. Subsequently, the protocol was applied to all 240,988 US pharmacological patents for which 9,111,706 molecules are available in SureChEMBL. The unsupervised filtering process selected 5,949,214 molecules (65.3% of the total number of molecules) that form highly congeneric chemical series in 188,795 of those patents (78.3% of the total number of patents). A SureChEMBL version enriched with molecules of pharmacological relevance is available for download at https://ftp.ebi.ac.uk/pub/databases/chembl/SureChEMBLccs.

Long Covid-19: Proposed Primary Care Clinical Guidelines for Diagnosis and Disease Management

Long COVID-19 may be defined as patients who, four weeks after the diagnosis of SARS-Cov-2 infection, continue to have signs and symptoms not explainable by other causes. The estimated frequency is around 10% and signs and symptoms may last for months. The main long-term manifestations observed in other coronaviruses (Severe Acute Respiratory Syndrome (SARS), Middle East respiratory syndrome (MERS)) are very similar to and have clear clinical parallels with SARS-CoV-2: mainly respiratory, musculoskeletal, and neuropsychiatric. The growing number of patients worldwide will have an impact on health systems. Therefore, the main objective of these clinical practice guidelines is to identify patients with signs and symptoms of long COVID-19 in primary care through a protocolized diagnostic process that studies possible etiologies and establishes an accurate differential diagnosis. The guidelines have been developed pragmatically by compiling the few studies published so far on long COVID-19, editorials and expert opinions, press releases, and the authors' clinical experience. Patients with long COVID-19 should be managed using structured primary care visits based on the time from diagnosis of SARS-CoV-2 infection. Based on the current limited evidence, disease management of long COVID-19 signs and symptoms will require a holistic, longitudinal follow up in primary care, multidisciplinary rehabilitation services, and the empowerment of affected patient groups.

A Workflow of Integrated Resources to Catalyze Network Pharmacology Driven COVID-19 Research

MOTIVATION

In the event of an outbreak due to an emerging pathogen, time is of the essence to contain or to mitigate the spread of the disease. Drug repositioning is one of the strategies that has the potential to deliver therapeutics relatively quickly. The SARS-CoV-2 pandemic has shown that integrating critical data resources to drive drug-repositioning studies, involving host-host, hostpathogen and drug-target interactions, remains a time-consuming effort that translates to a delay in the development and delivery of a life-saving therapy.

RESULTS

Here, we describe a workflow we designed for a semi-automated integration of rapidly emerging datasets that can be generally adopted in a broad network pharmacology research setting. The workflow was used to construct a COVID-19 focused multimodal network that integrates 487 host-pathogen, 74,805 host-host protein and 1,265 drug-target interactions. The resultant Neo4j graph database named "Neo4COVID19" is accessible via a web interface and via API calls based on the Bolt protocol. We believe that our Neo4COVID19 database will be a valuable asset to the research community and will catalyze the discovery of therapeutics to fight COVID-19.

AVAILABILITY

https://neo4covid19.ncats.io.

PRIMAGE project: predictive in silico multiscale analytics to support childhood cancer personalised evaluation empowered by imaging biomarkers

PRIMAGE is one of the largest and more ambitious research projects dealing with medical imaging, artificial intelligence and cancer treatment in children. It is a 4-year European Commission-financed project that has 16 European partners in the consortium, including the European Society for Paediatric Oncology, two imaging biobanks, and three prominent European paediatric oncology units. The project is constructed as an observational in silico study involving high-quality anonymised datasets (imaging, clinical, molecular, and genetics) for the training and validation of machine learning and multiscale algorithms. The open cloud-based platform will offer precise clinical assistance for phenotyping (diagnosis), treatment allocation (prediction), and patient endpoints (prognosis), based on the use of imaging biomarkers, tumour growth simulation, advanced visualisation of confidence scores, and machine-learning approaches. The decision support prototype will be constructed and validated on two paediatric cancers: neuroblastoma and diffuse intrinsic pontine glioma. External validation will be performed on data recruited from independent collaborative centres. Final results will be available for the scientific community at the end of the project, and ready for translation to other malignant solid tumours.

Closing the Gap Between Therapeutic Use and Mode of Action in Remedial Herbs

The ancient tradition of taking parts of a plant or preparing plant extracts for treating certain discomforts and maladies has long been lacking a scientific rationale to support its preparation and still widespread use in several parts of the world. In an attempt to address this challenge, we collected and integrated data connecting metabolites, plants, diseases, and proteins. A mechanistic hypothesis is generated when a metabolite is known to be present in a given plant, that plant is known to be used to treat a certain disease, that disease is known to be linked to the function of a given protein, and that protein is finally known or predicted to interact with the original metabolite. The construction of plant–protein networks from mutually connected metabolites and diseases facilitated the identification of plausible mechanisms of action for plants being used to treat analgesia, hypercholesterolemia, diarrhea, catarrh, and cough. Additional concrete examples using both experimentally known and computationally predicted, and subsequently experimentally confirmed, metabolite–protein interactions to close the connection circle between metabolites, plants, diseases, and proteins offered further proof of concept for the validity and scope of the approach to generate mode of action hypotheses for some of the therapeutic uses of remedial herbs.

Comparison of Zebrafish Larvae and hiPSC Cardiomyocytes for Predicting Drug-Induced Cardiotoxicity in Humans

Cardiovascular drug toxicity is responsible for 17% of drug withdrawals in clinical phases, half of post-marketed drug withdrawals and remains an important adverse effect of several marketed drugs. Early assessment of drug-induced cardiovascular toxicity is mandatory and typically done in cellular systems and mammals. Current in vitro screening methods allow high-throughput but are biologically reductionist. The use of mammal models, which allow a better translatability for predicting clinical outputs, is low-throughput, highly expensive, and ethically controversial. Given the analogies between the human and the zebrafish cardiovascular systems, we propose the use of zebrafish larvae during early drug discovery phases as a balanced model between biological translatability and screening throughput for addressing potential liabilities. To this end, we have developed a high-throughput screening platform that enables fully automatized in vivo image acquisition and analysis to extract a plethora of relevant cardiovascular parameters: heart rate, arrhythmia, AV blockage, ejection fraction, and blood flow, among others. We have used this platform to address the predictive power of zebrafish larvae for detecting potential cardiovascular liabilities in humans. We tested a chemical library of 92 compounds with known clinical cardiotoxicity profiles. The cross-comparison with clinical data and data acquired from human induced pluripotent stem cell cardiomyocytes calcium imaging showed that zebrafish larvae allow a more reliable prediction of cardiotoxicity than cellular systems. Interestingly, our analysis with zebrafish yields similar predictive performance as previous validation meta-studies performed with dogs, the standard regulatory preclinical model for predicting cardiotoxic liabilities prior to clinical phases.

The human endogenous metabolome as a pharmacology baseline for drug discovery

We have limited understanding of the variation in in vitro affinities of drugs for their targets. An analysis of a highly curated set of 815 interactions between 566 drugs and 129 primary targets reveals that 71% of drug-target affinities have values above that of the corresponding endogenous ligand, 96% of them fitting within a range of two orders of magnitude. Our findings suggest that the evolutionary optimised affinity of endogenous ligands for their native proteins can serve as a baseline for the primary pharmacology of drugs. We show that the degree of off-target selectivity and safety risks of drugs derived from their secondary pharmacology depend very much on that baseline. Thus, we propose a new approach for estimating safety margins.

Information Loss in Network Pharmacology

With the advent of increasing computational power and large-scale data acquisition, network analysis has become an attractive tool to study the organisation of complex systems and the interrelation of their constituent entities in various scientific domains. In many cases, relations only occur between entities of two different subsets, thereby forming a bipartite network. Often, the analysis of such bipartite networks involves the consideration of its two monopartite projections in order to focus on each entity subset individually as a means to deduce properties of the underlying original network. Although it is broadly acknowledged that this type of projection is not lossless, the inherent limitations of their interpretability are rarely discussed. In this work, we introduce two approaches for measuring the information loss associated with bipartite network projection. Application to two structurally distinct cases in network pharmacology, namely, drug-target and disease-gene bipartite networks, confirms that the major determinant of information loss is the degree of vertices omitted during the monopartite projection.