Improving the odds of drug development success through human genomics: modelling study

Genetic Insights into Therapeutic Targets for Neuromyelitis Optica Spectrum Disorders: A Mendelian Randomization Study

Neuromyelitis optica spectrum disorder (NMOSD) is a severe central nervous system disease primarily characterized by optic neuritis and myelitis, which can result in vision loss and limb paralysis. Current treatment options are limited in their ability to prevent relapses and mitigate disease progression, underscoring the urgent need for new drug targets to develop more effective therapies. The objective of this study is to identify potential drug targets associated with a reduced risk of NMOSD attacks or relapses through Mendelian randomization (MR) analysis, thereby addressing the limitations of existing treatment methods and providing better clinical options for patients. To identify therapeutic targets for NMOSD, a MR analysis was conducted. The cis-expression quantitative trait loci (cis-eQTL, exposure) data were sourced from the eQTLGen consortium, which included a sample size of 31,684. NMOSD (outcome) summary data were obtained from two of the largest independent cohorts: one cohort consisted of 86 NMOSD cases and 460 controls derived from whole-genome sequencing data, while the other cohort included 129 NMOSD patients and 784 controls. We performed a two-sample MR analysis to evaluate the association between single nucleotide polymorphisms (SNPs) and copy number variations with NMOSD. The MR analysis utilized the inverse variance weighted (IVW) method, supplemented by MR-Egger, weighted median, simple mode, and weighted mode methods. Sensitivity analyses were conducted to assess the presence of horizontal pleiotropy and heterogeneity. Colocalization analysis was employed to test whether NMOSD risk and gene expression are driven by common SNPs. Additionally, a phenome-wide association study (PheWAS) was performed to detect disease outcomes associated with NEU1. Supplementary analyses included single-nucleus RNA sequencing (snRNA-seq) data analysis, protein-protein interaction (PPI) networks, and drug feasibility assessments to prioritize potential therapeutic targets. Two drug targets, COL4A1 and NEU1, demonstrated significant MR results in two independent datasets. Notably, NEU1 showed substantial evidence of colocalization with NMOSD. Additionally, apart from the association between NEU1 and NMOSD, no other associations were observed between gene-proxied NEU1 inhibition and the increased risk of other NMOSD-related diseases. This study supports the potential of targeting NEU1 for drug inhibition to reduce the risk of NMOSD. Further preclinical research and drug development are warranted to validate the efficacy and safety of NEU1 as a therapeutic target and to explore its potential in NMOSD treatment.

Domain adaptable language modeling of chemical compounds identifies potent pathoblockers for Pseudomonas aeruginosa

Computational techniques for predicting molecular properties are emerging as key components for streamlining drug development, optimizing time and financial investments. Here, we introduce ChemLM, a transformer language model for this task. ChemLM leverages self-supervised domain adaptation on chemical molecules to enhance its predictive performance. Within the framework of ChemLM, chemical compounds are conceptualized as sentences composed of distinct chemical 'words', which are employed for training a specialized chemical language model. On the standard benchmark datasets, ChemLM either matched or surpassed the performance of current state-of-the-art methods. Furthermore, we evaluated the effectiveness of ChemLM in identifying highly potent pathoblockers targeting Pseudomonas aeruginosa (PA), a pathogen that has shown an increased prevalence of multidrug-resistant strains and has been identified as a critical priority for the development of new medications. ChemLM demonstrated substantially higher accuracy in identifying highly potent pathoblockers against PA when compared to state-of-the-art approaches. An intrinsic evaluation demonstrated the consistency of the chemical language model's representation concerning chemical properties. The results from benchmarking, experimental data and intrinsic analysis of the ChemLM space confirm the wide applicability of ChemLM for enhancing molecular property prediction within the chemical domain.

Impact of Glucagon-like Peptide-1 Receptor Agonists on Mental Illness: Evidence from a Mendelian Randomization Study

Emerging evidence suggests that glucagon-like peptide-1 receptor (GLP1R) agonists may have potential benefits for mental illnesses. However, their exact effects remain unclear. This study investigated the causal relationship between glucagon-like peptide-1 receptor agonist (GLP1RA) and the risk of 10 common mental illnesses, including attention deficit and hyperactivity disorder, anorexia nervosa, anxiety disorder, autism spectrum disorder, bipolar disorder, major depressive disorder, post-traumatic stress disorder, schizophrenia, cannabis use disorder, and alcohol use disorder. We selected GLP1RA as the exposure and conducted a Mendelian randomization (MR) analysis. The cis-eQTLs of the drug target gene GLP1R, provided by eQTLGen, were used to simulate the pharmacological effects of GLP1RA. Type 2 diabetes and BMI were included as positive controls. Using data from both the Psychiatric Genomic Consortium and FinnGen, we conducted separate MR analyses for the same disease across these two independent databases. Meta-analysis was used to pool the results. We found genetic evidence suggesting a causal relationship between GLP1RA and a reduced risk of schizophrenia [OR (95% CI) = 0.84 (0.71-0.98), I2 = 0.0%, common effects model]. Further mediation analysis indicated that this effect might be unrelated to improvements in glycemic control but rather mediated by BMI. However, the findings of this study provide insufficient evidence to support a causal relationship between GLP1RA and other mental illnesses. Sensitivity analyses did not reveal any potential bias due to horizontal pleiotropy or heterogeneity in the above results (p > 0.05). This study suggests that genetically proxied activation of glucagon-like peptide-1 receptor is associated with a lower risk of schizophrenia. GLP1R is implicated in schizophrenia pathogenesis, and its agonists may exert potential benefits through weight management. Our study provides useful information for understanding the neuropsychiatric effects of GLP1RA, which may contribute to refining future research designs and guiding clinical management. Moreover, our findings could have significant implications for overweight individuals at high risk of schizophrenia when selecting weight-loss medications. Future research should further investigate the potential mechanisms underlying the relationship between GLP1RA and schizophrenia.

Drug target assessments: classifying target modulation and associated health effects using multi-level BERT-based classification models

Motivation

Drug target selection determines the success of the drug development pipeline. Therefore, novel drug targets need to be assessed for their therapeutic benefits/risks at the earliest stage possible. Where manual risk/benefit analyses are often user-biased and time-consuming, Large Language Models can offer a systematic and efficient approach to curating and analysing literature. Currently, publicly available Large Language Models are lacking for this task, while public platforms for target assessments are limited to co-occurrences.

Results

BERT-models for multi-level classification of drug target-health effect relationships described in PubMed were developed. Relationships were classified based on (i) causality; (ii) direction of target modulation; (iii) direction of the associated health effect. The models showed competitive performances with F1 scores between 0.86 and 0.92 and their applicability was demonstrated using ADAM33 and OSM as case study. The developed classification pipeline is the first to allow detailed classification of drug target-health effect relationships. The models provide mechanistic insight into how target modulation affects health and disease, both from an efficacy and safety perspective. The models, deployed on the whole of PubMed and available through the TargetTri platform, are expected to offer a significant advancement in artificial intelligence-assisted target identification and evaluation.

Availability and implementation

https://www.targettri.com.

Evolution, genetic diversity, and health

Human genetic diversity in today's world has been shaped by evolutionary history, demographic shifts and environmental exposures, influencing complex traits, disease susceptibility and drug responses. Capturing this diversity is essential for advancing precision medicine and promoting equitable healthcare. Despite the great progress achieved with initiatives such as the human Pangenome and large biobanks that aim for a better representation of human diversity, important challenges remain. In this Perspective, we discuss the importance of diversity in clinical genomics through an evolutionary lens. We highlight progress and challenges and outline key clinical applications of diverse genetic data. We argue that diversifying both datasets and methodologies-integrating ancestral and environmental factors-is crucial for fully understanding the genetic basis of human health and disease.

Accelerating antiviral drug discovery: early hazard detection with a dual zebrafish and cell culture screen of a 403 compound library

The constant emergence of new viral pathogens underscores the need for continually evolving, effective antiviral drugs. A key challenge is identifying compounds that are both efficacious and safe, as many candidates fail during development due to unforeseen toxicity. To address this, the embryonic zebrafish morphology, mortality, and behavior (ZBE) screen and the SYSTEMETRIC® Cell Health Screen (CHS) were employed to evaluate the safety of 403 compounds from the Cayman Antiviral Screening Library. Of these compounds, 114 were FDA-approved, 17 were discontinued, and 97 remained on the market. CHS identified 25% (104 compounds) as toxic, with a Cell Health Index™ (CHI) > 0.5. The embryonic zebrafish model identified an additional 20% as toxic (79), bringing the total to 183. ZBEscreen flagged 19 toxic hits among compounds still on the market, seven of which were also identified by CHS. The combined use of CHS and zebrafish models enhanced hazard detection. Together, CHS and ZBEscreen identified 45.5% of the library as potentially hazardous. Notably, the zebrafish non-hazardous compounds correlated strongly with over-the-counter or prescribed antiviral drugs, confirming their known safety profile. Over 130 hazard-associated compounds warranted further investigation. Using self-organizing maps, six distinct neighborhoods of compound similarity were identified. This dual approach streamlined the early detection of hazards associated with promising leads and is expected to facilitate faster, safer antiviral discovery.

Therapeutic Drug Distribution across the Mouse Brain Is Heterogeneous as Revealed by In Vivo, Spatially Resolved Aptamer-Based Sensing

Drug discovery for central nervous system (CNS) targets is a high stakes process with estimated success rates below ten percent. Dose scaling, penetration through the blood-brain-barrier (BBB), and potency are among the various challenges involved in developing drugs for CNS targets. The standard approach to evaluate some of these parameters is dosing lead therapeutic compounds via intravenous delivery and assessing their brain levels via tissue homogenization and ex vivo quantification. Although a cost and time effective approach, brain homogenization lacks pharmacokinetic spatial resolution and normalizes drug levels to the entire brain volume. The brain, however, is known to have regional differences in cellular composition, transporters, BBB permeability, and drug-metabolizing enzymes, factors that could significantly affect pharmacological assessments during drug discovery. In this study we employ electrochemical aptamer-based sensors, a technology that allows in situ, real-time molecular monitoring in live animals, to reveal significant differences in the pharmacokinetics of drug uptake and accumulation in the brain of mice. Using vancomycin in the context of penetrating brain injury (PBI), our results highlight that potency may be significantly affected by PBI location. Additionally, more accurate dose scaling and delivery for deep brain wounds could be achieved by adjusting route of administration based on real-time-measured pharmacokinetic profiles, for example by changing delivery from intravenous to intracerebroventricular dosing. We emphasize the issue of establishing accurate pharmacological parameters during preclinical drug discovery efforts and underline the value of aptamer-based sensors for precise estimations of drug pharmacokinetics, transport across the BBB, and effective dose delivery during preclinical trials.

The role and application of bioinformatics techniques and tools in drug discovery

The process of drug discovery and development is both lengthy and intricate, demanding a substantial investment of time and financial resources. Bioinformatics techniques and tools can not only accelerate the identification of drug targets and the screening and refinement of drug candidates, but also facilitate the characterization of side effects and the prediction of drug resistance. High-throughput data from genomics, transcriptomics, proteomics, and metabolomics make significant contributions to mechanics-based drug discovery and drug reuse. This paper summarizes bioinformatics technologies and tools in drug research and development and their roles and applications in drug research and development, aiming to provide references for the development of new drugs and the realization of precision medicine.

Human Genetic Evidence to Inform Clinical Development of IL-6 Signaling Inhibition for Abdominal Aortic Aneurysm

BACKGROUND

Abdominal aortic aneurysm (AAA) represents a significant cause of mortality, yet no medical therapies have proven efficacious. The aim of the current study was to leverage human genetic evidence to inform clinical development of IL-6 (interleukin-6) signaling inhibition for the treatment of AAA.

METHODS

Associations of rs2228145, a missense variant in the IL6R gene region, are expressed per additional copy of the C allele, corresponding to the genetically predicted effect of IL-6 signaling inhibition. We consider genetic associations with AAA risk in the AAAgen consortium (39 221 cases and 1 086 107 controls) and UK Biobank (1963 cases and 365 680 controls). To validate against known effects of IL-6 signaling inhibition, we present associations with rheumatoid arthritis, polymyalgia rheumatica, and severe COVID-19. To explore mechanism specificity, we present associations with thoracic aortic aneurysm, intracranial aneurysm, and coronary artery disease. We further explored genetic associations in clinically relevant subgroups of the population.

RESULTS

We observed strong genetic associations with AAA risk in the AAAgen consortium, UK Biobank, and FinnGen (odds ratios: 0.91 [95% CI, 0.90-0.92], P=4×10-30; 0.90 [95% CI, 0.84-0.96], P=0.001; and 0.86 [95% CI, 0.82-0.91], P=7×10-9, respectively). The association was similar for fatal AAA but with greater uncertainty due to the lower number of events. The association with AAA was of greater magnitude than associations with coronary artery disease and even rheumatological disorders for which IL-6 inhibitors have been approved. No strong associations were observed with thoracic aortic aneurysm or intracranial aneurysm. Associations attenuated toward the null in populations with concomitant rheumatological or connective tissue disease.

CONCLUSIONS

Inhibition of IL-6 signaling is a promising strategy for treating AAA but not other types of aneurysmal disease. These findings serve to help inform clinical development of IL-6 signaling inhibition for AAA treatment.

Advancing tumor microenvironment and lymphoid tissue research through 3D bioprinting and biofabrication

Cancer progression is significantly influenced by the complex interactions within the tumor microenvironment (TME). Immune cells, in particular, play a critical role by infiltrating tumors from the circulation and surrounding lymphoid tissues in an attempt to control their spread. However, they often fail in this task. Current in vivo and in vitro preclinical models struggle to fully capture these intricate interactions affecting our ability to understand immune evasion and predict drugs behaviour in the clinic. To address this challenge, biofabrication and particularly 3D bioprinting has emerged as a promising tool for modeling both tumors and the immune system. Its ability to incorporate multiple cell types into 3D matrices, enable tissue compartmentalization with high spatial accuracy, and integrate vasculature makes it a valuable approach. Nevertheless, limited research has focused on capturing the complex tumor-immune interplay in vitro. This review highlights the composition and significance of the TME, the architecture and function of lymphoid tissues, and innovative approaches to modeling their interactions in vitro, while proposing the concept of an extended TME.

Combining Single-Cell RNA Sequencing and Mendelian Randomization to Explore Novel Drug Targets for Parkinson's Disease

Neuroinflammation is a key pathological factor of PD, and T cells play a central role in neuroinflammatory progression. However, the causal effect of T cell-related genes on the risk of PD is still unclear. We explored single-cell RNA sequencing (scRNA-Seq) datasets of the peripheral blood T cells of PD patients and healthy controls, and screened the differentially expressed genes (DEGs) in the cytotoxic CD4 + T cells relative to the other T cell subsets. Pseudo-time series analysis, cell-cell communication analysis, and metabolic pathway analysis was performed for the cytotoxic CD4 + T cells. The DEGs were also functionally annotated through GO and KEGG pathway enrichment analyses. The MR approach was used to establish causal effects of the DEGs (exposure) on PD risk (outcome), and explore new drug targets for PD. The findings of MR analysis were further validated by Steiger filtering, bidirectional MR, Bayesian colocalization analysis, and phenotype scanning, and the GWAS data from an independent PD case-control cohort was used for external validation of the results. Finally, differences in gene expression between PD patients and healthy controls were further validated in scRNA-Seq and bulk transcriptome sequencing data. We found that increased expression of IL-32, GNLY, MT2A, and ARPC2 was significantly associated with a higher risk of PD. In contrast, the increase in ARRB2 was closely related to a lower risk of PD. IL32, GNLY, MT2A, ARRB2, and ARPC2 are the causal genes and potential drug targets of PD. Cytotoxic CD4 + T cells are likely the key effectors of PD-related neuroinflammation. These findings provide new insights into the pathogenesis and treatment options for PD, and further research and clinical trials based on the five potential drug targets and neuroinflammation are necessary.

Comprehensive evaluation of pure and hybrid collaborative filtering in drug repurposing

Drug development is known to be a costly and time-consuming process, which is prone to high failure rates. Drug repurposing allows drug discovery by reusing already approved compounds. The outcomes of past clinical trials can be used to predict novel drug-disease associations by leveraging drug- and disease-related similarities. To tackle this classification problem, collaborative filtering with implicit feedback (and potentially additional data on drugs and diseases) has become popular. It can handle large imbalances between negative and positive known associations and known and unknown associations. However, properly evaluating the improvement over the state of the art is challenging, as there is no consensus approach to compare models. We propose a reproducible methodology for comparing collaborative filtering-based drug repurposing. We illustrate this method by comparing 11 models from the literature on eight diverse drug repurposing datasets. Based on this benchmark, we derive guidelines to ensure a fair and comprehensive evaluation of the performance of those models. In particular, an uncontrolled bias on unknown associations might lead to severe data leakage and a misestimation of the model's true performance. Moreover, in drug repurposing, the ability of a model to extrapolate beyond its training distribution is crucial and should also be assessed. Finally, we identified a subcategory of collaborative filtering that seems efficient and robust to distribution shifts. Benchmarks constitute an essential step towards increased reproducibility and more accessible development of competitive drug repurposing methods.

Metabolite Associations with Childhood and Juvenile Absence Epilepsy: A Bidirectional Mendelian Randomization Study

Background

The precise involvement of metabolites in the pathogenesis of Childhood absence epilepsy (CAE) and juvenile absence epilepsy (JAE) remains elusive. Consequently, this investigation introduces bidirectional Mendelian randomization (MR) as a tool to explore causality and underlying mechanisms.

Methods

Bidirectional MR analysis was conducted employing a comprehensive set comprising 1091 human blood metabolites and 309 metabolite ratios, systematically probing potential causal associations with JAE and CAE. Genome-wide association study (GWAS) data pertaining to these epileptic conditions were meticulously obtained from the International League Against Epilepsy (ILAE) consortium. Sensitivity analyses were rigorously performed to evaluate for heterogeneity and pleiotropy. Reverse MR analysis was also conducted to verify the direction of causality, and no significant reverse causal relationships were identified.

Results

Following rigorous genetic variant selection, significant associations were identified based on PIVW < .05, PWM < .05, and PMR-Egger < .05 criteria in MR analysis. Only 1 metabolite, (2 or 3)-decaonate levels, exhibited an association with JAE (P = .005, OR=0.987, 95% CI=0.978-0.996). Childhood absence epilepsy was associated with 5 metabolites: X-23648 (P = .012, OR=0.982, 95% CI=0.968-0.996), X-21845 levels (P = .045, OR=1.018, 95% CI=1.001-1.035), 2'-o-methylcytidine (P = .008, OR=0.995, 95% CI=0.991-1.001), 2'-o-methyluridine (P = .007, OR=0.995, 95% CI=0.99-0.999), and spermidine-topyruvate ratio (P = .014, OR=0.973, 95% CI=0.954-0.992). No evidence of reverse causality was found between JAE and CAE and the aforementioned metabolites.

Conclusion

The study establishes causal relationships between the aforementioned 6 metabolites and CAE and JAE. This integration of genomics with metabolism offers novel insights into epilepsy mechanisms and has important implications for screening and prevention.

Mechanism-free repurposing of drugs for C9orf72-related ALS/FTD using large-scale genomic data

Repeat expansions in the C9orf72 gene are the most common genetic cause of (ALS) and frontotemporal dementia (FTD). Like other genetic forms of neurodegeneration, pinpointing the precise mechanism(s) by which this mutation leads to neuronal death remains elusive, and this lack of knowledge hampers the development of therapy for C9orf72-related disease. We used an agnostic approach based on genomic data (n = 41,273 ALS and healthy samples, and n = 1,516 C9orf72 carriers) to overcome these bottlenecks. Our drug-repurposing screen, based on gene- and expression-pattern matching and information about the genetic variants influencing onset age among C9orf72 carriers, identified acamprosate, a γ-aminobutyric acid analog, as a potentially repurposable treatment for patients carrying C9orf72 repeat expansions. We validated its neuroprotective effect in cell models and showed comparable efficacy to riluzole, the current standard of care. Our work highlights the potential value of genomics in repurposing drugs in situations where the underlying pathomechanisms are inherently complex. VIDEO ABSTRACT.

Drug risks associated with sarcopenia: a real-world and GWAS study

INTRODUCTION

Drug-induced sarcopenia has not received adequate attention. Meanwhile, there is growing recognition of the importance of effective pharmacovigilance in evaluating the benefits and risks of medications.

AIMS

The primary aim of this study is to investigate the potential association between drug use and sarcopenia through an analysis of adverse event reports from the Food and Drug Administration (FDA) Adverse Event Reporting System (FAERS) and to evaluate the genetic factors contributing to drug-induced sarcopenia using summary-data-based Mendelian randomization (SMR).

METHODS

We obtained reports of adverse drug reactions from FAERS. Primary outcomes included sarcopenia and potential sarcopenia. We calculated the Proportional reporting ratio (PRR) to assess the risk of specific adverse events associated with various drugs, applying chi-square tests for statistical significance. Additionally, we used SMR based on Genome-wide association study (GWAS) to evaluate the potential associations between drug target genes of some significant medications and sarcopenia outcomes. The outcome data for sarcopenia included metrics as hand grip strength and appendicular lean mass (ALM).

RESULTS

A total of 55 drugs were identified as inducing potential sarcopenia, and 3 drugs were identified as inducing sarcopenia. The top 5 drugs causing a potential risk of sarcopenia were levofloxacin (PRR = 9.96, χ2 = 1057), pregabalin (PRR = 7.20, χ2 = 1023), atorvastatin (PRR = 4.68, χ2 = 903), duloxetine (PRR = 4.76, χ2 = 527) and venlafaxine (PRR = 5.56, χ2 = 504), and the 3 drugs that had been proved to induced sarcopenia included metformin (PRR = 7.41, χ2 = 58), aspirin (PRR = 5.93, χ2 = 35), and acetaminophen (PRR = 4.73, χ2 = 25). We identified electron-transfer flavoprotein dehydrogenase (ETFDH) and protein Kinase AMP-Activated Non-Catalytic Subunit Beta 1 (PRKAB1) as the primary drug target genes for metformin, while Prostaglandin-endoperoxide Synthase 1 (PTGS1) and Prostaglandin-endoperoxide Synthase 2 (PTGS2) were considered the primary action target genes for aspirin and acetaminophen according to DrugBank database. SMR showed that the expression abundance of ETFDH was negatively correlated with right hand grip strength (blood: OR = 1.01, p-value = 1.27e-02; muscle: OR = 1.01, p-value = 1.42e-02) and negatively correlated with appendicular lean mass (blood: OR = 1.03, p-value = 7.73e-08; muscle: OR = 1.03, p-value = 1.67e-07).

CONCLUSIONS

We find that metformin, aspirin, and acetaminophen are specifically noted for their potential to induce sarcopenia based on the analyses conducted. We perform signal mining for drug-associated sarcopenia events based on real-world data and provides certain guidance for the safe use of medications to prevent sarcopenia.

Validation of TYK2 and exploration of PRSS36 as drug targets for psoriasis using Mendelian randomization

Psoriasis is a chronic inflammatory skin disorder with multiple causes, including genetic and environmental factors. Despite advances in treatment, there remains a need to identify novel therapeutic targets. A Mendelian randomization (MR) analysis was conducted to identify therapeutic targets for psoriasis. Data on cis-expression quantitative trait loci were obtained from the eQTLGen Consortium (n = 31,684). Summary statistics for psoriasis (outcome) were sourced from the GWAS Catalog with a sample size of 484,598, including 5,427 cases and 479,171 controls. Colocalization analysis was used to assess whether psoriasis risk and gene expression were driven by shared single nucleotide polymorphisms. Drug prediction and molecular docking were utilized to validate the pharmacological value of the drug targets. The MR analysis found that 81 drug targets were significantly associated, and two (TYK2 and PRSS36) were supported by colocalization analysis (PP.H4 > 0.80). Phenome-wide association studies did not show any associations with other traits at the gene level. Biologically, these genes were closely related to immune function. Molecular docking revealed strong binding with drugs and proteins, as supported by available structural data. This study validated TYK2 as a drug target for psoriasis, in line with its existing clinical use, including the development of decucravacitinib. PRSS36 is a potential novel target requiring further investigation.

Disease coverage of human genome-wide association studies and pharmaceutical research and development

BACKGROUND

Despite the growing interest in the use of human genomic data for drug target identification and validation, the extent to which the spectrum of human disease has been addressed by genome-wide association studies (GWAS), or by drug development, and the degree to which these efforts overlap remain unclear.

METHODS

In this study we harmonize and integrate different data sources to create a sample space of all the human drug targets and diseases and identify points of convergence or divergence of GWAS and drug development efforts.

RESULTS

We show that only 612 of 11,158 diseases listed in Human Disease Ontology have an approved drug treatment in at least one region of the world. Of the 1414 diseases that are the subject of preclinical or clinical phase drug development, only 666 have been investigated in GWAS. Conversely, of the 1914 human diseases that have been the subject of GWAS, 1121 have yet to be investigated in drug development.

CONCLUSIONS

We produce target-disease indication lists to help the pharmaceutical industry to prioritize future drug development efforts based on genetic evidence, academia to prioritize future GWAS for diseases without effective treatments, and both sectors to harness genetic evidence to expand the indications for licensed drugs or to identify repurposing opportunities for clinical candidates that failed in their originally intended indication.

Integrating text mining with network models for successful target identification: in vitro validation in MASH-induced liver fibrosis

An in silico target discovery pipeline was developed by including a directional and weighted molecular disease network for metabolic dysfunction-associated steatohepatitis (MASH)-induced liver fibrosis. This approach integrates text mining, network biology, and artificial intelligence/machine learning with clinical transcriptome data for optimal translational power. At the mechanistic level, the critical components influencing disease progression were identified from the disease network using in silico knockouts. The top-ranked genes were then subjected to a target efficacy analysis, following which the top-5 candidate targets were validated in vitro. Three targets, including EP300, were confirmed for their roles in liver fibrosis. EP300 gene-silencing was found to significantly reduce collagen by 37%; compound intervention studies performed in human primary hepatic stellate cells and the hepatic stellate cell line LX-2 showed significant inhibition of collagen to the extent of 81% compared to the TGFβ-stimulated control (1 μM inobrodib in LX-2 cells). The validated in silico pipeline presents a unique approach for the identification of human-disease-mechanism-relevant drug targets. The directionality of the network ensures adherence to physiologically relevant signaling cascades, while the inclusion of clinical data boosts its translational power and ensures identification of the most relevant disease pathways. In silico knockouts thus provide crucial molecular insights for successful target identification.

Protein Structure Inspired Discovery of a Novel Inducer of Anoikis in Human Melanoma

Drug discovery historically starts with an established function, either that of compounds or proteins. This can hamper discovery of novel therapeutics. As structure determines function, we hypothesized that unique 3D protein structures constitute primary data that can inform novel discovery. Using a computationally intensive physics-based analytical platform operating at supercomputing speeds, we probed a high-resolution protein X-ray crystallographic library developed by us. For each of the eight identified novel 3D structures, we analyzed binding of sixty million compounds. Top-ranking compounds were acquired and screened for efficacy against breast, prostate, colon, or lung cancer, and for toxicity on normal human bone marrow stem cells, both using eight-day colony formation assays. Effective and non-toxic compounds segregated to two pockets. One compound, Dxr2-017, exhibited selective anti-melanoma activity in the NCI-60 cell line screen. In eight-day assays, Dxr2-017 had an IC50 of 12 nM against melanoma cells, while concentrations over 2100-fold higher had minimal stem cell toxicity. Dxr2-017 induced anoikis, a unique form of programmed cell death in need of targeted therapeutics. Our findings demonstrate proof-of-concept that protein structures represent high-value primary data to support the discovery of novel acting therapeutics. This approach is widely applicable.

Potential drug targets for osteoporosis identified: A Mendelian randomization study

Background

Osteoporosis is a prevalent global health condition, primarily affecting the aging population, and several therapies for osteoporosis have been widely used. However, available drugs for osteoporosis are far from satisfactory because they cannot alleviate disease progression. This study aimed to explore potential drug targets for osteoporosis through Mendelian randomization analysis.

Methods

Using cis-expression quantitative trait loci (cis-eQTL) data of druggable genes and two genome-wide association studies (GWAS) datasets related to osteoporosis (UK Biobank and FinnGen cohorts), we employed mendelian randomization (MR) analysis to identify the druggable genes with causal relationships with osteoporosis. Subsequently, a series of follow-up analyses were conducted, such as colocalization analysis, cell-type specificity analysis, and correlation analysis with risk factors. The association between potential drug targets and osteoporosis was validated by qRT-PCR.

Results

Six druggable genes with causal relationships with osteoporosis were identified and successfully replicated, including ACPP, DNASE1L3, IL32, PPOX, ST6GAL1, and TGM3. Cell-type specificity analysis revealed that PPOX and ST6GAL1 were expressed in all cell types in the bone samples, while IL32, ACPP, DNASE1L3, and TGM3 were expressed in specific cell types. The GWAS data showed there were seven risk factors for osteoporosis, including vitamin D deficiency, COPD, physical activity, BMI, MMP-9, ALP and PTH. Furthermore, ACPP was associated with vitamin D deficiency and COPD; DNASE1L3 was linked to physical activity; IL32 correlated with BMI and MMP-9; and ST6GAL1 was related to ALP, physical activity, and MMP-9. Among these risk factors, only MMP-9 had a high genetic correlation with osteoporosis. The results of qRT-PCR demonstrated that IL32 was upregulated while ST6GAL1 was downregulated in peripheral blood of osteoporosis patients.

Conclusion

Our findings suggested that those six druggable genes offer potential drug targets for osteoporosis and require further clinical investigation, especially IL32 and ST6GAL1.

Identification of potential therapeutic targets for skin cutaneous melanoma on the basic of transcriptomics

BACKGROUND

Advanced skin cutaneous melanoma (SKCM) is responsible for the majority of skin cancer-related deaths. Apart from the rare BRAF V600F mutation, which can be targeted with specific drugs, there are currently no other novel effective therapeutic targets.

METHODS

We used SMR analysis with cis-expressed quantitative trait locus (cis-eQTL) as the exposure variable and SKCM as the outcome variable to identify potential therapeutic targets for SKCM. Colocalization assays and HEIDI tests are used to test whether SKCM risk and gene expression are driven by common SNPs. Replication analysis further validated the findings, and we also constructed protein-protein interaction networks to explore the relationship between the identified genes and known SKCM targets. Drug prediction and molecular docking further validated the medicinal value of drug targets. Transcriptome differential analysis further validated that there were differences between normal tissues and SKCM for the selected targets.

RESULTS

We identified 13 genes significantly associated with the risk of SKCM, including five protective genes and eight harmful genes. The HEIDI test and co-localization analysis further indicates a causal association between genes (SOX4, MAFF) and SKCM, categorized as Class 1 evidence targets. The remaining 11 genes, except for HELZ2 show a moderately causal association with SKCM, categorized as Class 2 evidence targets. Target druggability predictions from DGIdb suggest that SOX4, MAFF, ACSF3, CDK10, SPG7, and TCF25 are likely to be future drug targets.

CONCLUSION

The study provides genetic evidence for targeting available drug genes for the treatment of SKCM.

Epidermal Growth Factor Receptor Targeting in Colorectal Carcinoma: Antibodies and Patient-Derived Organoids as a Smart Model to Study Therapy Resistance

Colorectal cancer (CRC) is the second leading cause of cancer-related death worldwide. Therefore, the need for new therapeutic strategies is still a challenge. Surgery and chemotherapy represent the first-line interventions; nevertheless, the prognosis for metastatic CRC (mCRC) patients remains unacceptable. An important step towards targeted therapy came from the inhibition of the epidermal growth factor receptor (EGFR) pathway, by the anti-EGFR antibody, Cetuximab, or by specific tyrosine kinase inhibitors (TKI). Cetuximab, a mouse-human chimeric monoclonal antibody (mAb), binds to the extracellular domain of EGFR thus impairing EGFR-mediated signaling and reducing cell proliferation. TKI can affect the EGFR biochemical pathway at different steps along the signaling cascade. Apart from Cetuximab, other anti-EGFR mAbs have been developed, such as Panitumumab. Both antibodies have been approved for the treatment of KRAS-NRAS wild type mCRC, alone or in combination with chemotherapy. These antibodies display strong differences in activating the host immune system against CRC, due to their different immunoglobulin isotypes. Although anti-EGFR antibodies are efficient, drug resistance occurs with high frequency. Resistant tumor cell populations can either already be present before therapy or develop later by biochemical adaptations or new genomic mutations in the EGFR pathway. Numerous efforts have been made to improve the efficacy of the anti-EGFR mAbs or to find new agents that are able to block downstream EGFR signaling cascade molecules. Indeed, we examined the importance of analyzing the anti-EGFR antibody-drug conjugates (ADC) developed to overcome resistance and/or stimulate the tumor host's immunity against CRC growth. Also, patient-derived CRC organoid cultures represent a useful and feasible in vitro model to study tumor behavior and therapy response. Organoids can reflect tumor genetic heterogeneity found in the tissue of origin, representing a unique tool for personalized medicine. Thus, CRC-derived organoid cultures are a smart model for studying the tumor microenvironment and for the preclinical assay of anti-EGFR drugs.

Renshen Yangrong decoction for secondary malaise and fatigue: network pharmacology and Mendelian randomization study

Background

Renshen Yangrong decoction (RSYRD) has been shown therapeutic effects on secondary malaise and fatigue (SMF). However, to date, its bioactive ingredients and potential targets remain unclear.

Purpose

The purpose of this study is to assess the potential ingredients and targets of RSYRD on SMF through a comprehensive strategy integrating network pharmacology, Mendelian randomization as well as molecular docking verification.

Methods

Search for potential active ingredients and corresponding protein targets of RSYRD on TCMSP and BATMAN-TCM for network pharmacology analysis. Mendelian randomization (MR) was performed to find therapeutic targets for SMF. The eQTLGen Consortium (sample sizes: 31,684) provided data on cis-expression quantitative trait loci (cis-eQTL, exposure). The summary data on SMF (outcome) from genome-wide association studies (GWAS) were gathered from the MRC-IEU Consortium (sample sizes: 463,010). We built a target interaction network between the probable active ingredient targets of RSYRD and the therapeutic targets of SMF. We next used drug prediction and molecular docking to confirm the therapeutic value of the therapeutic targets.

Results

In RSYRD, network pharmacology investigations revealed 193 possible active compounds and 234 associated protein targets. The genetically predicted amounts of 176 proteins were related to SMF risk in the MR analysis. Thirty-seven overlapping targets for RSYRD in treating SMF, among which six (NOS3, GAA, IMPA1, P4HTM, RB1, and SLC16A1) were prioritized with the most convincing evidence. Finally, the 14 active ingredients of RSYRD were identified as potential drug molecules. The strong affinity between active components and putative protein targets was established by molecular docking.

Conclusion

This study revealed several active components and possible RSYRD protein targets for the therapy of SMF and provided novel insights into the feasibility of using Mendelian randomization for causal inference between Chinese medical formula and disease.

Base editing strategies to convert CAG to CAA diminish the disease-causing mutation in Huntington's disease

An expanded CAG repeat in the huntingtin gene (HTT) causes Huntington's disease (HD). Since the length of uninterrupted CAG repeat, not polyglutamine, determines the age-at-onset in HD, base editing strategies to convert CAG to CAA are anticipated to delay onset by shortening the uninterrupted CAG repeat. Here, we developed base editing strategies to convert CAG in the repeat to CAA and determined their molecular outcomes and effects on relevant disease phenotypes. Base editing strategies employing combinations of cytosine base editors and guide RNAs (gRNAs) efficiently converted CAG to CAA at various sites in the CAG repeat without generating significant indels, off-target edits, or transcriptome alterations, demonstrating their feasibility and specificity. Candidate BE strategies converted CAG to CAA on both expanded and non-expanded CAG repeats without altering HTT mRNA and protein levels. In addition, somatic CAG repeat expansion, which is the major disease driver in HD, was significantly decreased in the liver by a candidate BE strategy treatment in HD knock-in mice carrying canonical CAG repeats. Notably, CAG repeat expansion was abolished entirely in HD knock-in mice carrying CAA-interrupted repeats, supporting the therapeutic potential of CAG-to-CAA conversion strategies in HD and potentially other repeat expansion disorders.

Identification of drug targets for Sjögren's syndrome: multi-omics Mendelian randomization and colocalization analyses

Background

Targeted therapy for Sjögren's syndrome (SS) has become an important focus for clinicians. Multi-omics-wide Mendelian randomization (MR) analyses have provided new ideas for identifying potential drug targets.

Methods

We conducted summary-data-based Mendelian randomization (SMR) analysis to evaluate therapeutic targets associated with SS by integrating DNA methylation, gene expression and protein quantitative trait loci (mQTL, eQTL, and pQTL, respectively). Genetic associations with SS were derived from the FinnGen study (discovery) and the GWAS catalog (replication). Colocalization analyses were employed to determine whether two potentially relevant phenotypes share the same genetic factors in a given region. Moreover, to delve deeper into potential regulation among DNA methylation, gene expression, and protein abundance, we conducted MR analysis to explore the causal relationship between candidate gene methylation and expression, as well as between gene expression and protein abundance. Drug prediction and molecular docking were further employed to validate the pharmacological activity of the candidate drug targets.

Results

Upon integrating the multi-omics data, we identified three genes associated with SS risk: TNFAIP3, BTN3A1, and PLAU. The methylation of cg22068371 in BTN3A1 was positively associated with protein levels, consistent with the negative effect of cg22068371 methylation on the risk of SS. Additionally, positive correlations were observed between the gene methylation of PLAU (cg04939496) and expression, as well as between expression and protein levels. This consistency elucidates the promotional effects of PLAU on SS risk at the DNA methylation, gene expression, and protein levels. At the protein level, genetically predicted TNFAIP3 (OR 2.47, 95% CI 1.56-3.92) was positively associated with SS risk, while BTN3A1 (OR 2.96E-03, 95% CI 2.63E-04-3.33E-02) was negatively associated with SS risk. Molecular docking showed stable binding for candidate drugs and target proteins.

Conclusion

Our study reveals promising therapeutic targets for the treatment of SS, providing valuable insights into targeted therapy for SS. However, further validation through future experiments is warranted.

Protein Structure Inspired Drug Discovery

Drug discovery starts with known function, either of a compound or a protein, in-turn prompting investigations to probe 3D structure of the compound-protein interface. As protein structure determines function, we hypothesized that unique 3D structural motifs represent primary information denoting unique function that can drive discovery of novel agents. Using a physics-based protein structure analysis platform developed by us, designed to conduct computationally intensive analysis at supercomputing speeds, we probed a high-resolution protein x-ray crystallographic library developed by us. We selected 3D structural motifs whose function was not otherwise established, that offered environments supporting binding of drug-like chemicals and were present on proteins that were not established therapeutic targets. For each of eight potential binding pockets on six different proteins we accessed a 60 million compound library and used our analysis platform to evaluate binding. Using eight-day colony formation assays acquired compounds were screened for efficacy against human breast, prostate, colon and lung cancer cells and toxicity against human bone marrow stem cells. Compounds selectively inhibiting cancer growth segregated to two pockets on separate proteins. The compound, Dxr2-017, exhibited selective activity against human melanoma cells in the NCI-60 cell line screen, had an IC50 of 19 nM against human melanoma M14 cells in our eight-day assay, while over 2100-fold higher concentrations inhibited stem cells by less than 30%. We show that Dxr2-017 induces anoikis, a unique form of programmed cell death in need of targeted therapeutics. The predicted target protein for Dxr2-017 is expressed in bacteria, not in humans. This supports our strategy of focusing on unique 3D structural motifs. It is known that functionally important 3D structures are evolutionarily conserved. Here we demonstrate proof-of-concept that protein structure represents high value primary data to support discovery of novel therapeutics. This approach is widely applicable.

Refining the impact of genetic evidence on clinical success

The cost of drug discovery and development is driven primarily by failure1, with only about 10% of clinical programmes eventually receiving approval2-4. We previously estimated that human genetic evidence doubles the success rate from clinical development to approval5. In this study we leverage the growth in genetic evidence over the past decade to better understand the characteristics that distinguish clinical success and failure. We estimate the probability of success for drug mechanisms with genetic support is 2.6 times greater than those without. This relative success varies among therapy areas and development phases, and improves with increasing confidence in the causal gene, but is largely unaffected by genetic effect size, minor allele frequency or year of discovery. These results indicate we are far from reaching peak genetic insights to aid the discovery of targets for more effective drugs.

Virtual Screening of a Chemically Diverse "Superscaffold" Library Enables Ligand Discovery for a Key GPCR Target

The advent of ultra-large libraries of drug-like compounds has significantly broadened the possibilities in structure-based virtual screening, accelerating the discovery and optimization of high-quality lead chemotypes for diverse clinical targets. Compared to traditional high-throughput screening, which is constrained to libraries of approximately one million compounds, the ultra-large virtual screening approach offers substantial advantages in both cost and time efficiency. By expanding the chemical space with compounds synthesized from easily accessible and reproducible reactions and utilizing a large, diverse set of building blocks, we can enhance both the diversity and quality of the discovered lead chemotypes. In this study, we explore new chemical spaces using reactions of sulfur(VI) fluorides to create a combinatorial library consisting of several hundred million compounds. We screened this virtual library for cannabinoid type II receptor (CB2) antagonists using the high-resolution structure in conjunction with a rationally designed antagonist, AM10257. The top-predicted compounds were then synthesized and tested in vitro for CB2 binding and functional antagonism, achieving an experimentally validated hit rate of 55%. Our findings demonstrate the effectiveness of reliable reactions, such as sulfur fluoride exchange, in diversifying ultra-large chemical spaces and facilitate the discovery of new lead compounds for important biological targets.

Unveiling the Potential of Ent-Kaurane Diterpenoids: Multifaceted Natural Products for Drug Discovery

Natural products hold immense potential for drug discovery, yet many remain unexplored in vast libraries and databases. In an attempt to fill this gap and meet the growing demand for effective drugs, this study delves into the promising world of ent-kaurane diterpenoids, a class of natural products with huge therapeutic potential. With a dataset of 570 ent-kaurane diterpenoids obtained from the literature, we conducted an in silico analysis, evaluating their physicochemical, pharmacokinetic, and toxicological properties with a focus on their therapeutic implications. Notably, these natural compounds exhibit drug-like properties, aligning closely with those of FDA-approved drugs, indicating a high potential for drug development. The ranges of the physicochemical parameters were as follows: molecular weights-288.47 to 626.82 g/mol; number of heavy atoms-21 to 44; the number of hydrogen bond donors and acceptors-0 to 8 and 1 to 11, respectively; the number of rotatable bonds-0 to 11; fraction Csp3-0.65 to 1; and TPSA-20.23 to 189.53 Ų. Additionally, the majority of these molecules display favorable safety profiles, with only 0.70%, 1.40%, 0.70%, and 46.49% exhibiting mutagenic, tumorigenic, reproduction-enhancing, and irritant properties, respectively. Importantly, ent-kaurane diterpenoids exhibit promising biopharmaceutical properties. Their average lipophilicity is optimal for drug absorption, while over 99% are water-soluble, facilitating delivery. Further, 96.5% and 28.20% of these molecules exhibited intestinal and brain bioavailability, expanding their therapeutic reach. The predicted pharmacological activities of these compounds encompass a diverse range, including anticancer, immunosuppressant, chemoprotective, anti-hepatic, hepatoprotectant, anti-inflammation, antihyperthyroidism, and anti-hepatitis activities. This multi-targeted profile highlights ent-kaurane diterpenoids as highly promising candidates for further drug discovery endeavors.

Using Mendelian randomization provides genetic insights into potential targets for sepsis treatment

Sepsis is recognized as a major contributor to the global disease burden, but there is a lack of specific and effective therapeutic agents. Utilizing Mendelian randomization (MR) methods alongside evidence of causal genetics presents a chance to discover novel targets for therapeutic intervention. MR approach was employed to investigate potential drug targets for sepsis. Pooled statistics from IEU-B-4980 comprising 11,643 cases and 474,841 controls were initially utilized, and the findings were subsequently replicated in the IEU-B-69 (10,154 cases and 454,764 controls). Causal associations were then validated through colocalization. Furthermore, a range of sensitivity analyses, including MR-Egger intercept tests and Cochran's Q tests, were conducted to evaluate the outcomes of the MR analyses. Three drug targets (PSMA4, IFNAR2, and LY9) exhibited noteworthy MR outcomes in two separate datasets. Notably, PSMA4 demonstrated not only an elevated susceptibility to sepsis (OR 1.32, 95% CI 1.20-1.45, p = 1.66E-08) but also exhibited a robust colocalization with sepsis (PPH4 = 0.74). According to the present MR analysis, PSMA4 emerges as a highly encouraging pharmaceutical target for addressing sepsis. Suppression of PSMA4 could potentially decrease the likelihood of sepsis.

Predicting Protein Functions Based on Heterogeneous Graph Attention Technique

In bioinformatics, protein function prediction stands as a fundamental area of research and plays a crucial role in addressing various biological challenges, such as the identification of potential targets for drug discovery and the elucidation of disease mechanisms. However, known functional annotation databases usually provide positive experimental annotations that proteins carry out a given function, and rarely record negative experimental annotations that proteins do not carry out a given function. Therefore, existing computational methods based on deep learning models focus on these positive annotations for prediction and ignore these scarce but informative negative annotations, leading to an underestimation of precision. To address this issue, we introduce a deep learning method that utilizes a heterogeneous graph attention technique. The method first constructs a heterogeneous graph that covers the protein-protein interaction network, ontology structure, and positive and negative annotation information. Then, it learns embedding representations of proteins and ontology terms by using the heterogeneous graph attention technique. Finally, it leverages these learned representations to reconstruct the positive protein-term associations and score unobserved functional annotations. It can enhance the predictive performance by incorporating these known limited negative annotations into the constructed heterogeneous graph. Experimental results on three species (i.e., Human, Mouse, and Arabidopsis) demonstrate that our method can achieve better performance in predicting new protein annotations than state-of-the-art methods.

SGCLDGA: unveiling drug-gene associations through simple graph contrastive learning

Drug repurposing offers a viable strategy for discovering new drugs and therapeutic targets through the analysis of drug-gene interactions. However, traditional experimental methods are plagued by their costliness and inefficiency. Despite graph convolutional network (GCN)-based models' state-of-the-art performance in prediction, their reliance on supervised learning makes them vulnerable to data sparsity, a common challenge in drug discovery, further complicating model development. In this study, we propose SGCLDGA, a novel computational model leveraging graph neural networks and contrastive learning to predict unknown drug-gene associations. SGCLDGA employs GCNs to extract vector representations of drugs and genes from the original bipartite graph. Subsequently, singular value decomposition (SVD) is employed to enhance the graph and generate multiple views. The model performs contrastive learning across these views, optimizing vector representations through a contrastive loss function to better distinguish positive and negative samples. The final step involves utilizing inner product calculations to determine association scores between drugs and genes. Experimental results on the DGIdb4.0 dataset demonstrate SGCLDGA's superior performance compared with six state-of-the-art methods. Ablation studies and case analyses validate the significance of contrastive learning and SVD, highlighting SGCLDGA's potential in discovering new drug-gene associations. The code and dataset for SGCLDGA are freely available at https://github.com/one-melon/SGCLDGA.

Drug development advances in human genetics-based targets

Drug development is a long and costly process, with a high degree of uncertainty from the identification of a drug target to its market launch. Targeted drugs supported by human genetic evidence are expected to enter phase II/III clinical trials or be approved for marketing more quickly, speeding up the drug development process. Currently, genetic data and technologies such as genome-wide association studies (GWAS), whole-exome sequencing (WES), and whole-genome sequencing (WGS) have identified and validated many potential molecular targets associated with diseases. This review describes the structure, molecular biology, and drug development of human genetics-based validated beneficial loss-of-function (LOF) mutation targets (target mutations that reduce disease incidence) over the past decade. The feasibility of eight beneficial LOF mutation targets (PCSK9, ANGPTL3, ASGR1, HSD17B13, KHK, CIDEB, GPR75, and INHBE) as targets for drug discovery is mainly emphasized, and their research prospects and challenges are discussed. In conclusion, we expect that this review will inspire more researchers to use human genetics and genomics to support the discovery of novel therapeutic drugs and the direction of clinical development, which will contribute to the development of new drug discovery and drug repurposing.

Health Digital Twins in Life Science and Health Care Innovation

Health digital twins (HDTs) are virtual representations of real individuals that can be used to simulate human physiology, disease, and drug effects. HDTs can be used to improve drug discovery and development by providing a data-driven approach to inform target selection, drug delivery, and design of clinical trials. HDTs also offer new applications into precision therapies and clinical decision making. The deployment of HDTs at scale could bring a precision approach to public health monitoring and intervention. Next steps include challenges such as addressing socioeconomic barriers and ensuring the representativeness of the technology based on the training and validation data sets. Governance and regulation of HDT technology are still in the early stages.

Genetic interactions of schizophrenia using gene-based statistical epistasis exclusively identify nervous system-related pathways and key hub genes

Background: The relationship between genotype and phenotype is governed by numerous genetic interactions (GIs), and the mapping of GI networks is of interest for two main reasons: 1) By modelling biological robustness, GIs provide a powerful opportunity to infer compensatory biological mechanisms via the identification of functional relationships between genes, which is of interest for biological discovery and translational research. Biological systems have evolved to compensate for genetic (i.e., variations and mutations) and environmental (i.e., drug efficacy) perturbations by exploiting compensatory relationships between genes, pathways and biological processes; 2) GI facilitates the identification of the direction (alleviating or aggravating interactions) and magnitude of epistatic interactions that influence the phenotypic outcome. The generation of GIs for human diseases is impossible using experimental biology approaches such as systematic deletion analysis. Moreover, the generation of disease-specific GIs has never been undertaken in humans. Methods: We used our Indian schizophrenia case-control (case-816, controls-900) genetic dataset to implement the workflow. Standard GWAS sample quality control procedure was followed. We used the imputed genetic data to increase the SNP coverage to analyse epistatic interactions across the genome comprehensively. Using the odds ratio (OR), we identified the GIs that increase or decrease the risk of a disease phenotype. The SNP-based epistatic results were transformed into gene-based epistatic results. Results: We have developed a novel approach by conducting gene-based statistical epistatic analysis using an Indian schizophrenia case-control genetic dataset and transforming these results to infer GIs that increase the risk of schizophrenia. There were ∼9.5 million GIs with a p-value ≤ 1 × 10-5. Approximately 4.8 million GIs showed an increased risk (OR > 1.0), while ∼4.75 million GIs had a decreased risk (OR <1.0) for schizophrenia. Conclusion: Unlike model organisms, this approach is specifically viable in humans due to the availability of abundant disease-specific genome-wide genotype datasets. The study exclusively identified brain/nervous system-related processes, affirming the findings. This computational approach fills a critical gap by generating practically non-existent heritable disease-specific human GIs from human genetic data. These novel datasets can train innovative deep-learning models, potentially surpassing the limitations of conventional GWAS.

Development of a human genetics-guided priority score for 19,365 genes and 399 drug indications

Studies have shown that drug targets with human genetic support are more likely to succeed in clinical trials. Hence, a tool integrating genetic evidence to prioritize drug target genes is beneficial for drug discovery. We built a genetic priority score (GPS) by integrating eight genetic features with drug indications from the Open Targets and SIDER databases. The top 0.83%, 0.28% and 0.19% of the GPS conferred a 5.3-, 9.9- and 11.0-fold increased effect of having an indication, respectively. In addition, we observed that targets in the top 0.28% of the score were 1.7-, 3.7- and 8.8-fold more likely to advance from phase I to phases II, III and IV, respectively. Complementary to the GPS, we incorporated the direction of genetic effect and drug mechanism into a directional version of the score called the GPS with direction of effect. We applied our method to 19,365 protein-coding genes and 399 drug indications and made all results available through a web portal.

The citrate transporter SLC13A5 as a therapeutic target for kidney disease: evidence from Mendelian randomization to inform drug development

BACKGROUND

Solute carrier family 13 member 5 (SLC13A5) is a Na+-coupled citrate co-transporter that mediates entry of extracellular citrate into the cytosol. SLC13A5 inhibition has been proposed as a target for reducing progression of kidney disease. The aim of this study was to leverage the Mendelian randomization paradigm to gain insight into the effects of SLC13A5 inhibition in humans, towards prioritizing and informing clinical development efforts.

METHODS

The primary Mendelian randomization analyses investigated the effect of SLC13A5 inhibition on measures of kidney function, including creatinine and cystatin C-based measures of estimated glomerular filtration rate (creatinine-eGFR and cystatin C-eGFR), blood urea nitrogen (BUN), urine albumin-creatinine ratio (uACR), and risk of chronic kidney disease and microalbuminuria. Secondary analyses included a paired plasma and urine metabolome-wide association study, investigation of secondary traits related to SLC13A5 biology, a phenome-wide association study (PheWAS), and a proteome-wide association study. All analyses were compared to the effect of genetically predicted plasma citrate levels using variants selected from across the genome, and statistical sensitivity analyses robust to the inclusion of pleiotropic variants were also performed. Data were obtained from large-scale genetic consortia and biobanks, with sample sizes ranging from 5023 to 1,320,016 individuals.

RESULTS

We found evidence of associations between genetically proxied SLC13A5 inhibition and higher creatinine-eGFR (p = 0.002), cystatin C-eGFR (p = 0.005), and lower BUN (p = 3 × 10-4). Statistical sensitivity analyses robust to the inclusion of pleiotropic variants suggested that these effects may be a consequence of higher plasma citrate levels. There was no strong evidence of associations of genetically proxied SLC13A5 inhibition with uACR or risk of CKD or microalbuminuria. Secondary analyses identified evidence of associations with higher plasma calcium levels (p = 6 × 10-13) and lower fasting glucose (p = 0.02). PheWAS did not identify any safety concerns.

CONCLUSIONS

This Mendelian randomization analysis provides human-centric insight to guide clinical development of an SLC13A5 inhibitor. We identify plasma calcium and citrate as biologically plausible biomarkers of target engagement, and plasma citrate as a potential biomarker of mechanism of action. Our human genetic evidence corroborates evidence from various animal models to support effects of SLC13A5 inhibition on improving kidney function.

TIE1 and TEK signalling, intraocular pressure, and primary open-angle glaucoma: a Mendelian randomization study

BACKGROUND

In primary open-angle glaucoma (POAG), lowering intraocular pressure (IOP) is the only proven way of slowing vision loss. Schlemm's canal (SC) is a hybrid vascular and lymphatic vessel that mediates aqueous humour drainage from the anterior ocular chamber. Animal studies support the importance of SC endothelial angiopoietin-TEK signalling, and more recently TIE1 signalling, in maintaining normal IOP. However, human genetic support for a causal role of TIE1 and TEK signalling in lowering IOP is currently lacking.

METHODS

GWAS summary statistics were obtained for plasma soluble TIE1 (sTIE1) protein levels (N = 35,559), soluble TEK (sTEK) protein levels (N = 35,559), IOP (N = 139,555) and POAG (Ncases = 16,677, Ncontrols = 199,580). Mendelian randomization (MR) was performed to estimate the association of genetically proxied TIE1 and TEK protein levels with IOP and POAG liability. Where significant MR estimates were obtained, genetic colocalization was performed to assess the probability of a shared causal variant (PPshared) versus distinct (PPdistinct) causal variants underlying TIE1/TEK signalling and the outcome. Publicly available single-nucleus RNA-sequencing data were leveraged to investigate differential expression of TIE1 and TEK in the human ocular anterior segment.

RESULTS

Increased genetically proxied TIE1 signalling and TEK signalling associated with a reduction in IOP (- 0.21 mmHg per SD increase in sTIE1, 95% CI = - 0.09 to - 0.33 mmHg, P = 6.57 × 10-4, and - 0.14 mmHg per SD decrease in sTEK, 95% CI = - 0.03 to - 0.25 mmHg, P = 0.011), but not with POAG liability. Colocalization analysis found that the probability of a shared causal variant was greater for TIE1 and IOP than for TEK and IOP (PPshared/(PPdistinct + PPshared) = 0.98 for TIE1 and 0.30 for TEK). In the anterior segment, TIE1 and TEK were preferentially expressed in SC, lymphatic, and vascular endothelium.

CONCLUSIONS

This study provides novel human genetic support for a causal role of both TIE1 and TEK signalling in regulating IOP. Here, combined evidence from cis-MR and colocalization analyses provide stronger support for TIE1 than TEK as a potential IOP-lowering therapeutic target.

Optical Screening and Classification of Drug Binding to Proteins in Human Blood Serum

Protein-drug interactions in the human bloodstream are important factors in applications ranging from drug design, where protein binding influences efficacy and dose delivery, to biomedical diagnostics, where rapid, quantitative measurements could guide optimized treatment regimes. Current measurement approaches use multistep assays, which probe the protein-bound drug fraction indirectly and do not provide fundamental structural or dynamic information about the in vivo protein-drug interaction. We demonstrate that ultrafast 2D-IR spectroscopy can overcome these issues by providing a direct, label-free optical measurement of protein-drug binding in blood serum samples. Four commonly prescribed drugs, known to bind to human serum albumin (HSA), were added to pooled human serum at physiologically relevant concentrations. In each case, spectral changes to the amide I band of the serum sample were observed, consistent with binding to HSA, but were distinct for each of the four drugs. A machine-learning-based classification of the serum samples achieved a total cross-validation prediction accuracy of 92% when differentiating serum-only samples from those with a drug present. Identification on a per-drug basis achieved correct drug identification in 75% of cases. These unique spectroscopic signatures of the drug-protein interaction thus enable the detection and differentiation of drug containing samples and give structural insight into the binding process as well as quantitative information on protein-drug binding. Using currently available instrumentation, the 2D-IR data acquisition required just 1 min and 10 μL of serum per sample, and so these results pave the way to fast, specific, and quantitative measurements of protein-drug binding in vivo with potentially invaluable applications for the development of novel therapies and personalized medicine.

Impact of the Physical Cellular Microenvironment on the Structure and Function of a Model Hepatocyte Cell Line for Drug Toxicity Applications

It is widely recognised that cells respond to their microenvironment, which has implications for cell culture practices. Growth cues provided by 2D cell culture substrates are far removed from native 3D tissue structure in vivo. Geometry is one of many factors that differs between in vitro culture and in vivo cellular environments. Cultured cells are far removed from their native counterparts and lose some of their predictive capability and reliability. In this study, we examine the cellular processes that occur when a cell is cultured on 2D or 3D surfaces for a short period of 8 days prior to its use in functional assays, which we term: "priming". We follow the process of mechanotransduction from cytoskeletal alterations, to changes to nuclear structure, leading to alterations in gene expression, protein expression and improved functional capabilities. In this study, we utilise HepG2 cells as a hepatocyte model cell line, due to their robustness for drug toxicity screening. Here, we demonstrate enhanced functionality and improved drug toxicity profiles that better reflect the in vivo clinical response. However, findings more broadly reflect in vitro cell culture practises across many areas of cell biology, demonstrating the fundamental impact of mechanotransduction in bioengineering and cell biology.

Effects of bumetanide on neonatal seizures: A systematic review of animal and human studies

BACKGROUND

Bumetanide, an inhibitor of the sodium-potassium-chloride cotransporter-1, has been suggested as an adjunct to phenobarbital for treating neonatal seizures.

METHODS

A systematic review of animal and human studies was conducted to evaluate the efficacy and safety of bumetanide for neonatal seizures. PubMed, Embase, CINAHL and Cochrane databases were searched in March 2023.

RESULTS

26 animal (rat or mice) studies describing 38 experiments (28 in-vivo and ten in-vitro) and two human studies (one RCT and one open-label dose-finding) were included. The study designs, methods to induce seizures, bumetanide dose, and outcome measures were heterogeneous, with only 4/38 experiments being in animal hypoxia/ischaemia models. Among 38 animal experiments, bumetanide was reported to have antiseizure effects in 21, pro-seizure in six and ineffective in 11. The two human studies (n = 57) did not show the benefits of bumetanide as an add-on agent to phenobarbital in their primary analyses, but one study reported benefit on post-hoc analysis. Overall, hearing impairment was detected in 5/37 surviving infants in the bumetanide group vs. 0/13 in controls. Four of the five infants with hearing impairment had received aminoglycosides concurrently. Other adverse effects reported were diuresis, mild-to-moderate dehydration, hypotension, and electrolyte disturbances. The studies did not report on long-term neurodevelopment. The certainty of the evidence was very low.

CONCLUSION

Animal data suggest that bumetanide has inconsistent effects as an antiseizure medication in neonates. Data from human studies are scarce and raise some concerns regarding ototoxicity when given with aminoglycosides. Well conducted studies in animal models of hypoxic-ischaemic encephalopathy are urgently needed. Future RCTs, if conducted in human neonates, should have an adequate sample size, assess neurodevelopment, minimize using aminoglycosides, be transparent about the potential ototoxicity in the parent information sheet, conduct early hearing tests and have trial-stopping rules that include hearing impairment as an outcome.

Identification of potential drug targets for rheumatoid arthritis from genetic insights: a Mendelian randomization study

INTRODUCTION

Rheumatoid arthritis (RA) is a chronic inflammatory illness that mostly affects the joints of the hands and feet and can reduce life expectancy by an average of 3 to 10 years. Although tremendous progress has been achieved in the treatment of RA, a large minority of patients continue to respond poorly to existing medications, owing in part to a lack of appropriate therapeutic targets.

METHODS

To find therapeutic targets for RA, a Mendelian randomization (MR) was performed. Cis-expression quantitative trait loci (cis-eQTL, exposure) data were obtained from the eQTLGen Consortium (sample size 31,684). Summary statistics for RA (outcome) were obtained from two largest independent cohorts: sample sizes of 97,173 (22,350 cases and 74,823 controls) and 269,377 (8279 cases and 261,098), respectively. Colocalisation analysis was used to test whether RA risk and gene expression were driven by common SNPs. Drug prediction and molecular docking was further used to validate the medicinal value of drug targets.

RESULTS

Seven drug targets were significant in both cohorts in MR analysis and supported by localization. PheWAS at the gene level showed only ATP2A1 associated with other traits. These genes are strongly associated with immune function in terms of biological significance. Molecular docking showed excellent binding for drugs and proteins with available structural data.

CONCLUSION

This study identifies seven potential drug targets for RA. Drugs designed to target these genes have a higher chance of success in clinical trials and is expected to help prioritise RA drug development and save on drug development costs.

A drug repurposing method based on inhibition effect on gene regulatory network

Numerous computational drug repurposing methods have emerged as efficient alternatives to costly and time-consuming traditional drug discovery approaches. Some of these methods are based on the assumption that the candidate drug should have a reversal effect on disease-associated genes. However, such methods are not applicable in the case that there is limited overlap between disease-related genes and drug-perturbed genes. In this study, we proposed a novel Drug Repurposing method based on the Inhibition Effect on gene regulatory network (DRIE) to identify potential drugs for cancer treatment. DRIE integrated gene expression profile and gene regulatory network to calculate inhibition score by using the shortest path in the disease-specific network. The results on eleven datasets indicated the superior performance of DRIE when compared to other state-of-the-art methods. Case studies showed that our method effectively discovered novel drug-disease associations. Our findings demonstrated that the top-ranked drug candidates had been already validated by CTD database. Additionally, it clearly identified potential agents for three cancers (colorectal, breast, and lung cancer), which was beneficial when annotating drug-disease relationships in the CTD. This study proposed a novel framework for drug repurposing, which would be helpful for drug discovery and development.

Predictions from First-Principles of Membrane Permeability to Small Molecules: How Useful Are They in Practice?

Predicting from first-principles the rate of passive permeation of small molecules across the biological membrane represents a promising strategy for screening lead compounds upstream in the drug-discovery and development pipeline. One popular avenue for the estimation of permeation rates rests on computer simulations in conjunction with the inhomogeneous solubility-diffusion model, which requires the determination of the free-energy change and position-dependent diffusivity of the substrate along the translocation pathway through the lipid bilayer. In this Perspective, we will clarify the physical meaning of the membrane permeability inferred from such computer simulations, and how theoretical predictions actually relate to what is commonly measured experimentally. We will also examine why these calculations remain both technically challenging and overly computationally expensive, which has hitherto precluded their routine use in nonacademic settings. We finally synopsize possible research directions to meet these challenges, increase the predictive power of physics-based rates of passive permeation, and, by ricochet, improve their practical usefulness.

Biological and functional multimorbidity-from mechanisms to management

Globally, the number of people with multiple co-occurring diseases will increase substantially over the coming decades, with important consequences for patients, carers, healthcare systems and society. Addressing this challenge requires a shift in the prevailing clinical, educational and scientific thinking and organization-with a strong emphasis on the maintenance of generalist skills to balance the specialization trends of medical education and research. Multimorbidity is not a single entity but differs quantitively and qualitatively across life stages, ethnicities, sexes, socioeconomic groups and geographies. Data-driven science that quantifies the impact of disease co-occurrence-beyond the small number of currently well-studied long-term conditions (such as cardiometabolic diseases)-can help illuminate the pathological diversity of multimorbidity and identify common, mechanistically related, and prognostically relevant clusters. Broader access to data opportunities across modalities and disciplines will catalyze vertical and horizontal integration of multimorbidity research, to enable reconfiguring of medical services, clinical trials, guidelines and research in a way that accounts for the complexity of multimorbidity-and provides efficient, joined-up services for patients.

Renoprotective effects of genetically proxied fibroblast growth factor 21: Mendelian randomization, proteome-wide and metabolome-wide association study

BACKGROUND

Fibroblast growth factor 21 (FGF21) has demonstrated efficacy for reducing liver fat and reversing non-alcoholic steatohepatitis in phase 2 clinical trials. It is also postulated to have anti-fibrotic effects and therefore may be amenable to repurposing for the prevention and treatment of chronic kidney disease (CKD).

METHODS

We leverage a missense genetic variant, rs739320 in the FGF21 gene, that associates with magnetic resonance imaging-derived liver fat as a clinically validated and biologically plausible instrumental variable for studying the effects of FGF21 analogs. Performing Mendelian randomization, we ascertain associations between instrumented FGF21 and kidney phenotypes, cardiometabolic disease risk factors, as well as the circulating proteome (Somalogic, 4907 aptamers) and metabolome (Nightingale platform, 249 metabolites).

RESULTS

We report consistent renoprotective associations of genetically proxied FGF21 effect, including higher glomerular filtration rates (p = 1.9 × 10^-4), higher urinary sodium excretion (p = 5.1 × 10^-11), and lower urine albumin-creatinine ratio (p = 3.6 × 10^-5). These favorable effects translated to lower CKD risk (odds ratio per rs739320 C-allele, 0.96; 95%CI, 0.94-0.98; p = 3.2 × 10^-4). Genetically proxied FGF21 effect was also associated with lower fasting insulin, waist-to-hip ratio, blood pressure (systolic and diastolic BP, p < 1.0 × 10^-07) and blood lipid (low-density lipoprotein cholesterol, triglycerides and apolipoprotein B, p < 6.5 × 10^-24) profiles. The latter associations are replicated in our metabolome-wide association study. Proteomic perturbations associated with genetically predicted FGF21 effect were consistent with fibrosis reduction.

CONCLUSION

This study highlights the pleiotropic effects of genetically proxied FGF21 and supports a re-purposing opportunity for the treatment and prevention of kidney disease specifically. Further work is required to triangulate these findings, towards possible clinical development of FGF21 towards the treatment and prevention of kidney disease.

Exploring the Advantages of Quantum Generative Adversarial Networks in Generative Chemistry

De novo drug design with desired biological activities is crucial for developing novel therapeutics for patients. The drug development process is time- and resource-consuming, and it has a low probability of success. Recent advances in machine learning and deep learning technology have reduced the time and cost of the discovery process and therefore, improved pharmaceutical research and development. In this paper, we explore the combination of two rapidly developing fields with lead candidate discovery in the drug development process. First, artificial intelligence has already been demonstrated to successfully accelerate conventional drug design approaches. Second, quantum computing has demonstrated promising potential in different applications, such as quantum chemistry, combinatorial optimizations, and machine learning. This article explores hybrid quantum-classical generative adversarial networks (GAN) for small molecule discovery. We substituted each element of GAN with a variational quantum circuit (VQC) and demonstrated the quantum advantages in the small drug discovery. Utilizing a VQC in the noise generator of a GAN to generate small molecules achieves better physicochemical properties and performance in the goal-directed benchmark than the classical counterpart. Moreover, we demonstrate the potential of a VQC with only tens of learnable parameters in the generator of GAN to generate small molecules. We also demonstrate the quantum advantage of a VQC in the discriminator of GAN. In this hybrid model, the number of learnable parameters is significantly less than the classical ones, and it can still generate valid molecules. The hybrid model with only tens of training parameters in the quantum discriminator outperforms the MLP-based one in terms of both generated molecule properties and the achieved KL divergence. However, the hybrid quantum-classical GANs still face challenges in generating unique and valid molecules compared to their classical counterparts.

Genetic evidence implicating natriuretic peptide receptor-3 in cardiovascular disease risk: a Mendelian randomization study

BACKGROUND

C-type natriuretic peptide (CNP) is a known target for promoting growth and has been implicated as a therapeutic opportunity for the prevention and treatment of cardiovascular disease (CVD). This study aimed to explore the effect of CNP on CVD risk using the Mendelian randomization (MR) framework.

METHODS

Instrumental variables mimicking the effects of pharmacological intervention on CNP were identified as uncorrelated genetic variants located in the genes coding for its primary receptors, natriuretic peptide receptors-2 and 3 (NPR2 and NPR3), that associated with height. We performed MR and colocalization analyses to investigate the effects of NPR2 signalling and NPR3 function on CVD outcomes and risk factors. MR estimates were compared to those obtained when considering height variants from throughout the genome.

RESULTS

Genetically-proxied reduced NPR3 function was associated with a lower risk of CVD, with odds ratio (OR) 0.74 per standard deviation (SD) higher NPR3-predicted height, and 95% confidence interval (95% CI) 0.64-0.86. This effect was greater in magnitude than observed when considering height variants from throughout the genome. For CVD subtypes, similar MR associations for NPR3-predicted height were observed when considering the outcomes of coronary artery disease (0.75, 95% CI 0.60-0.92), stroke (0.69, 95% CI 0.50-0.95) and heart failure (0.77, 95% CI 0.58-1.02). Consideration of CVD risk factors identified systolic blood pressure (SBP) as a potential mediator of the NPR3-related CVD risk lowering. For stroke, we found that the MR estimate for NPR3 was greater in magnitude than could be explained by a genetically predicted SBP effect alone. Colocalization results largely supported the MR findings, with no evidence of results being driven by effects due to variants in linkage disequilibrium. There was no MR evidence supporting effects of NPR2 on CVD risk, although this null finding could be attributable to fewer genetic variants being identified to instrument this target.

CONCLUSIONS

This genetic analysis supports the cardioprotective effects of pharmacologically inhibiting NPR3 receptor function, which is only partly mediated by an effect on blood pressure. There was unlikely sufficient statistical power to investigate the cardioprotective effects of NPR2 signalling.

Explainable Drug Repurposing Approach From Biased Random Walks

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

Plug-and-Play Lymph Node-on-Chip: Secondary Tumor Modeling by the Combination of Cell Spheroid, Collagen Sponge and T-Cells

Towards the improvement of the efficient study of drugs and contrast agents, the 3D microfluidic platforms are currently being actively developed for testing these substances and particles in vitro. Here, we have elaborated a microfluidic lymph node-on-chip (LNOC) as a tissue engineered model of a secondary tumor in lymph node (LN) formed due to the metastasis process. The developed chip has a collagen sponge with a 3D spheroid of 4T1 cells located inside, simulating secondary tumor in the lymphoid tissue. This collagen sponge has a morphology and porosity comparable to that of a native human LN. To demonstrate the suitability of the obtained chip for pharmacological applications, we used it to evaluate the effect of contrast agent/drug carrier size, on the penetration and accumulation of particles in 3D spheroids modeling secondary tumor. For this, the 0.3, 0.5 and 4 μm bovine serum albumin (BSA)/tannic acid (TA) capsules were mixed with lymphocytes and pumped through the developed chip. The capsule penetration was examined by scanning with fluorescence microscopy followed by quantitative image analysis. The results show that capsules with a size of 0.3 μm passed more easily to the tumor spheroid and penetrated inside. We hope that the device will represent a reliable alternative to in vivo early secondary tumor models and decrease the amount of in vivo experiments in the frame of preclinical study.