Ligand-based virtual-screening identified a novel CFTR ligand which improves the defective cell surface expression of misfolded ABC transporters.

Cystic fibrosis (CF) is a monogenetic disease caused by the mutation of CFTR, a cAMP-regulated Cl- channel expressing at the apical plasma membrane (PM) of epithelia. ∆F508-CFTR, the most common mutant in CF, fails to reach the PM due to its misfolding and premature degradation at the endoplasmic reticulum (ER). Recently, CFTR modulators have been developed to correct CFTR abnormalities, with some being used as therapeutic agents for CF treatment. One notable example is Trikafta, a triple combination of CFTR modulators (TEZ/ELX/IVA), which significantly enhances the functionality of ΔF508-CFTR on the PM. However, there's room for improvement in its therapeutic effectiveness since TEZ/ELX/IVA doesn't fully stabilize ΔF508-CFTR on the PM. To discover new CFTR modulators, we conducted a virtual screening of approximately 4.3 million compounds based on the chemical structures of existing CFTR modulators. This effort led us to identify a novel CFTR ligand named FR3. Unlike clinically available CFTR modulators, FR3 appears to operate through a distinct mechanism of action. FR3 enhances the functional expression of ΔF508-CFTR on the apical PM in airway epithelial cell lines by stabilizing NBD1. Notably, FR3 counteracted the degradation of mature ΔF508-CFTR, which still occurs despite the presence of TEZ/ELX/IVA. Furthermore, FR3 corrected the defective PM expression of a misfolded ABCB1 mutant. Therefore, FR3 may be a potential lead compound for addressing diseases resulting from the misfolding of ABC transporters.

Galectin-10 in serum extracellular vesicles reflects asthma pathophysiology.

BACKGROUND: Novel biomarkers (BMs) are urgently needed for bronchial asthma (BA) with various phenotypes and endotypes. OBJECTIVE: We sought to identify novel BMs reflecting tissue pathology from serum extracellular vesicles (EVs). METHODS: We performed data-independent acquisition of serum EVs from 4 healthy controls, 4 noneosinophilic asthma (NEA) patients, and 4 eosinophilic asthma (EA) patients to identify novel BMs for BA. We confirmed EA-specific BMs via data-independent acquisition validation in 61 BA patients and 23 controls. To further validate these findings, we performed data-independent acquisition for 6 patients with chronic rhinosinusitis without nasal polyps and 7 patients with chronic rhinosinusitis with nasal polyps. RESULTS: We identified 3032 proteins, 23 of which exhibited differential expression in EA. Ingenuity pathway analysis revealed that protein signatures from each phenotype reflected disease characteristics. Validation revealed 5 EA-specific BMs, including galectin-10 (Gal10), eosinophil peroxidase, major basic protein, eosinophil-derived neurotoxin, and arachidonate 15-lipoxygenase. The potential of Gal10 in EVs was superior to that of eosinophils in terms of diagnostic capability and detection of airway obstruction. In rhinosinusitis patients, 1752 and 8413 proteins were identified from EVs and tissues, respectively. Among 11 BMs identified in EVs and tissues from patients with chronic rhinosinusitis with nasal polyps, 5 (including Gal10 and eosinophil peroxidase) showed significant correlations between EVs and tissues. Gal10 release from EVs was implicated in eosinophil extracellular trapped cell death in vitro and in vivo. CONCLUSION: Novel BMs such as Gal10 from serum EVs reflect disease pathophysiology in BA and may represent a new target for liquid biopsy approaches.

Automatic generation of functional peptides with desired bioactivity and membrane permeability using Bayesian optimization.

Peptides are potentially useful modalities of drugs; however, cell membrane permeability is an obstacle in peptide drug discovery. The identification of bioactive peptides for a therapeutic target is also challenging because of the huge amino acid sequence patterns of peptides. In this study, we propose a novel computational method, PEptide generation system using Neural network Trained on Amino acid sequence data and Gaussian process-based optimizatiON (PENTAGON), to automatically generate new peptides with desired bioactivity and cell membrane permeability. In the algorithm, we mapped peptide amino acid sequences onto the latent space constructed using a variational autoencoder and searched for peptides with desired bioactivity and cell membrane permeability using Bayesian optimization. We used our proposed method to generate peptides with cell membrane permeability and bioactivity for each of the nine therapeutic targets, such as the estrogen receptor (ER). Our proposed method outperformed a previously developed peptide generator in terms of similarity to known active peptide sequences and the length of generated peptide sequences.

Network-based identification of diagnosis-specific trans-omic biomarkers via integration of multiple omics data.

The integration of multiple omics data promises to reveal new insights into the pathogenic mechanisms of complex human diseases, with the potential to identify avenues for the development of targeted therapies for disease subtypes. However, the extraction of diagnostic/disease-specific biomarkers from multiple omics data with biological pathway knowledge is a challenging issue in precision medicine. In this paper, we present a novel computational method to identify diagnosis-specific trans-omic biomarkers from multiple omics data. In the algorithm, we integrated multi-class sparse canonical correlation analysis (MSCCA) and molecular pathway analysis in order to derive discriminative molecular features that are correlated across different omics layers. We applied our proposed method to analyzing proteome and metabolome data of heart failure (HF), and extracted trans-omic biomarkers for HF subtypes; specifically, ischemic cardiomyopathy (ICM) and dilated cardiomyopathy (DCM). We were able to detect not only individual proteins that were previously reported from single-omics studies but also correlated protein-metabolite pairs characteristic of HF disease subtypes. For example, we identified hexokinase1(HK1)-d-fructose-6-phosphate as a paired trans-omic biomarker for DCM, which could significantly perturb amino-sugar metabolism. Our proposed method is expected to be useful for various applications in precision medicine.

DIRECTEUR: transcriptome-based prediction of small molecules that replace transcription factors for direct cell conversion.

MOTIVATION: Direct reprogramming (DR) is a process that directly converts somatic cells to target cells. Although DR via small molecules is safer than using transcription factors (TFs) in terms of avoidance of tumorigenic risk, the determination of DR-inducing small molecules is challenging. RESULTS: Here we present a novel in silico method, DIRECTEUR, to predict small molecules that replace TFs for DR. We extracted DR-characteristic genes using transcriptome profiles of cells in which DR was induced by TFs, and performed a variant of simulated annealing to explore small molecule combinations with similar gene expression patterns with DR-inducing TFs. We applied DIRECTEUR to predicting combinations of small molecules that convert fibroblasts into neurons or cardiomyocytes, and were able to reproduce experimentally verified and functionally related molecules inducing the corresponding conversions. The proposed method is expected to be useful for practical applications in regenerative medicine. AVAILABILITY AND IMPLEMENTATION: The code and data are available at the following link: https://github.com/HamanoLaboratory/DIRECTEUR.git.

A community effort in SARS-CoV-2 drug discovery.

The COVID-19 pandemic continues to pose a substantial threat to human lives and is likely to do so for years to come. Despite the availability of vaccines, searching for efficient small-molecule drugs that are widely available, including in low- and middle-income countries, is an ongoing challenge. In this work, we report the results of an open science community effort, the "Billion molecules against COVID-19 challenge", to identify small-molecule inhibitors against SARS-CoV-2 or relevant human receptors. Participating teams used a wide variety of computational methods to screen a minimum of 1 billion virtual molecules against 6 protein targets. Overall, 31 teams participated, and they suggested a total of 639,024 molecules, which were subsequently ranked to find 'consensus compounds'. The organizing team coordinated with various contract research organizations (CROs) and collaborating institutions to synthesize and test 878 compounds for biological activity against proteases (Nsp5, Nsp3, TMPRSS2), nucleocapsid N, RdRP (only the Nsp12 domain), and (alpha) spike protein S. Overall, 27 compounds with weak inhibition/binding were experimentally identified by binding-, cleavage-, and/or viral suppression assays and are presented here. Open science approaches such as the one presented here contribute to the knowledge base of future drug discovery efforts in finding better SARS-CoV-2 treatments.

De Novo Generation of Chemical Structures of Inhibitor and Activator Candidates for Therapeutic Target Proteins by a Transformer-Based Variational Autoencoder and Bayesian Optimization.

Deep generative models for molecular generation have been gaining much attention as structure generators to accelerate drug discovery. However, most previously developed methods are chemistry-centric approaches, and comprehensive biological responses in the cell have not been taken into account. In this study, we propose a novel computational method, TRIOMPHE-BOA (transcriptome-based inference and generation of molecules with desired phenotypes using the Bayesian optimization algorithm), to generate new chemical structures of inhibitor or activator candidates for therapeutic target proteins by integrating chemically and genetically perturbed transcriptome profiles. In the algorithm, the substructures of multiple molecules that were selected based on the transcriptome analysis are fused in the design of new chemical structures by exploring the latent space of a Transformer-based variational autoencoder using Bayesian optimization. Our results demonstrate the usefulness of the proposed method in terms of having high reproducibility of existing ligands for 10 therapeutic target proteins when compared with previous methods. Moreover, this method can be applied to proteins without detailed 3D structures or known ligands and is expected to become a powerful tool for more efficient hit identification.

Myocardial DNA Damage Predicts Heart Failure Outcome in Various Underlying Diseases.

BACKGROUND: Reliable predictors of treatment efficacy in heart failure have been long awaited. DNA damage has been implicated as a cause of heart failure. OBJECTIVES: The purpose of this study was to investigate the association of DNA damage in myocardial tissue with treatment response and prognosis of heart failure. METHODS: The authors performed immunostaining of DNA damage markers poly(ADP-ribose) (PAR) and γ-H2A.X in endomyocardial biopsy specimens from 175 patients with heart failure with reduced ejection fraction (HFrEF) of various underlying etiologies. They calculated the percentage of nuclei positive for each DNA damage marker (%PAR and %γ-H2A.X). The primary outcome was left ventricular reverse remodeling (LVRR) at 1 year, and the secondary outcome was a composite of cardiovascular death, heart transplantation, and ventricular assist device implantation. RESULTS: Patients who did not achieve LVRR after the optimization of medical therapies presented with significantly higher %PAR and %γ-H2A.X. The ROC analysis demonstrated good performance of both %PAR and %γ-H2A.X for predicting LVRR (AUCs: 0.867 and 0.855, respectively). There was a negative correlation between the mean proportion of DNA damage marker-positive nuclei and the probability of LVRR across different underlying diseases. In addition, patients with higher %PAR or %γ-H2A.X had more long-term clinical events (PAR HR: 1.63 [95% CI: 1.31-2.01]; P < 0.001; γ-H2A.X HR: 1.48 [95% CI: 1.27-1.72]; P < 0.001). CONCLUSIONS: DNA damage determines the consequences of human heart failure. Assessment of DNA damage is useful to predict treatment efficacy and prognosis of heart failure patients with various underlying etiologies.

Elucidating disease-associated mechanisms triggered by pollutants via the epigenetic landscape using large-scale ChIP-Seq data.

BACKGROUND: Despite well-documented effects on human health, the action modes of environmental pollutants are incompletely understood. Although transcriptome-based approaches are widely used to predict associations between chemicals and disorders, the molecular cues regulating pollutant-derived gene expression changes remain unclear. Therefore, we developed a data-mining approach, termed "DAR-ChIPEA," to identify transcription factors (TFs) playing pivotal roles in the action modes of pollutants. METHODS: Large-scale public ChIP-Seq data (human, n = 15,155; mouse, n = 13,156) were used to predict TFs that are enriched in the pollutant-induced differentially accessible genomic regions (DARs) obtained from epigenome analyses (ATAC-Seq). The resultant pollutant-TF matrices were then cross-referenced to a repository of TF-disorder associations to account for pollutant modes of action. We subsequently evaluated the performance of the proposed method using a chemical perturbation data set to compare the outputs of the DAR-ChIPEA and our previously developed differentially expressed gene (DEG)-ChIPEA methods using pollutant-induced DEGs as input. We then adopted the proposed method to predict disease-associated mechanisms triggered by pollutants. RESULTS: The proposed approach outperformed other methods using the area under the receiver operating characteristic curve score. The mean score of the proposed DAR-ChIPEA was significantly higher than that of our previously described DEG-ChIPEA (0.7287 vs. 0.7060; Q = 5.278 × 10-42; two-tailed Wilcoxon rank-sum test). The proposed approach further predicted TF-driven modes of action upon pollutant exposure, indicating that (1) TFs regulating Th1/2 cell homeostasis are integral in the pathophysiology of tributyltin-induced allergic disorders; (2) fine particulates (PM2.5) inhibit the binding of C/EBPs, Rela, and Spi1 to the genome, thereby perturbing normal blood cell differentiation and leading to immune dysfunction; and (3) lead induces fatty liver by disrupting the normal regulation of lipid metabolism by altering hepatic circadian rhythms. CONCLUSIONS: Highlighting genome-wide chromatin change upon pollutant exposure to elucidate the epigenetic landscape of pollutant responses outperformed our previously described method that focuses on gene-adjacent domains only. Our approach has the potential to reveal pivotal TFs that mediate deleterious effects of pollutants, thereby facilitating the development of strategies to mitigate damage from environmental pollution.

Design and structural optimization of thiadiazole derivatives with potent GLS1 inhibitory activity.

GLS1 is an attractive target not only as anticancer agents but also as candidates for various potential pharmaceutical applications such as anti-aging and anti-obesity treatments. We performed docking simulations based on the complex crystal structure of GLS1 and its inhibitor CB-839 and found that compound A bearing a thiadiazole skeleton exhibits GLS1 inhibition. Furthermore, we synthesized 27 thiadiazole derivatives in an effort to obtain a more potent GLS1 inhibitor. Among the synthesized derivatives, 4d showed more potent GLS1 inhibitory activity (IC50 of 46.7 µM) than known GLS1 inhibitor DON and A. Therefore, 4d is a very promising novel GLS1 inhibitor.

De novo drug design based on patient gene expression profiles via deep learning.

Computational de novo drug design is a challenging issue in medicine, and it is desirable to consider all of the relevant information of the biological systems in a disease state. Here, we propose a novel computational method to generate drug candidate molecular structures from patient gene expression profiles via deep learning, which we call DRAGONET. Our model can generate new molecules that are likely to counteract disease-specific gene expression patterns in patients, which is made possible by exploring the latent space constructed by a transformer-based variational autoencoder and integrating the substructures of disease-correlated molecules. We applied DRAGONET to generate drug candidate molecules for gastric cancer, atopic dermatitis, and Alzheimer's disease, and demonstrated that the newly generated molecules were chemically similar to registered drugs for each disease. This approach is applicable to diseases with unknown therapeutic target proteins and will make a significant contribution to the field of precision medicine.

Association between immune checkpoint inhibitor-induced myocarditis and concomitant use of thiazide diuretics.

Although an association has been reported between diuretics and myocarditis, it is unclear whether the risk of immune checkpoint inhibitor (ICI)-induced myocarditis is affected by concomitant diuretics. Thus, the aim of this work was to evaluate the impact of concomitant diuretics on ICI-induced myocarditis. This cross-sectional study used disproportionality analysis and a pharmacovigilance database to assess the risk of myocarditis with various diuretics in patients receiving ICIs via the analysis of data entered into the VigiBase database through December 2022. Multiple logistic regression analysis was performed to identify risk factors for myocarditis in patients who received ICIs. A total of 90 611 patients who received ICIs, including 975 cases of myocarditis, were included as the eligible dataset. A disproportionality in myocarditis was observed for loop diuretic use (reporting odds ratio 1.47, 95% confidence interval [CI] 1.02-2.04, P = .03) and thiazide use (reporting odds ratio 1.76, 95% CI 1.20-2.50, P < .01) in patients who received ICIs. The results of the multiple logistic regression analysis showed that the use of thiazides (odds ratio 1.67, 95% CI 1.15-2.34, P < .01) was associated with an increased risk of myocarditis in patients who received ICIs. Our findings may help to predict the risk of myocarditis in patients receiving ICIs.

TEAD1 trapping by the Q353R-Lamin A/C causes dilated cardiomyopathy.

Mutations in the LMNA gene encoding Lamin A and C (Lamin A/C), major components of the nuclear lamina, cause laminopathies including dilated cardiomyopathy (DCM), but the underlying molecular mechanisms have not been fully elucidated. Here, by leveraging single-cell RNA sequencing (RNA-seq), assay for transposase-accessible chromatin using sequencing (ATAC-seq), protein array, and electron microscopy analysis, we show that insufficient structural maturation of cardiomyocytes owing to trapping of transcription factor TEA domain transcription factor 1 (TEAD1) by mutant Lamin A/C at the nuclear membrane underlies the pathogenesis of Q353R-LMNA-related DCM. Inhibition of the Hippo pathway rescued the dysregulation of cardiac developmental genes by TEAD1 in LMNA mutant cardiomyocytes. Single-cell RNA-seq of cardiac tissues from patients with DCM with the LMNA mutation confirmed the dysregulated expression of TEAD1 target genes. Our results propose an intervention for transcriptional dysregulation as a potential treatment of LMNA-related DCM.

A trial of topiramate for patients with hereditary spinocerebellar ataxia.

In an open pilot trial, six patients with various hereditary forms of spinocerebellar ataxia (SCA) were assigned to topiramate (50 mg/day) for 24 weeks. Four patients completed the protocol without adverse events. Of these four patients, topiramate was effective for three patients. Some patients with SCA could respond to treatment with topiramate.

Characterization of Immune Checkpoint Inhibitor-Induced Myasthenia Gravis Using the US Food and Drug Administration Adverse Event Reporting System.

Myasthenia gravis (MG) is a rare but fatal adverse event of immune checkpoint inhibitors (ICIs). We assessed whether patient characteristics differed between those with ICI-related myasthenia gravis and those with idiopathic myasthenia gravis. Reports from the US Food and Drug Administration Adverse Event Reporting System were analyzed. Multivariate analyses were conducted to evaluate the associations between age, sex, and ICI treatment and the reporting rate of myasthenia gravis. Among 5 464 099 cases between 2011 and 2019, 53 447 were treated with ICIs. Myasthenia gravis was reported more often in ICI users. Multiple logistic regression analyses showed that the reporting rate of ICI-related myasthenia gravis did not differ significantly between men and women; however, it was higher in older people than in younger people (adjusted odds ratio, 2.4 [95%CI, 1.84-3.13]). We also investigated useful signs for the early detection of myositis and myocarditis, which are fatal when overlapping with ICI-related myasthenia gravis. Patients with elevated serum creatine kinase or troponin levels were more likely to have concurrent myositis and myocarditis. Unlike idiopathic myasthenia gravis, there was no sex difference in the development of ICI-related myasthenia gravis, which may be more common in older people. Considering the physiological muscle weakness that occurs in the elderly, it may be necessary to monitor ICI-related myasthenia gravis more closely in older people.

Differential effects of proton pump inhibitors and vonoprazan on vascular endothelial growth factor expression in cancer cells.

Proton pump inhibitors (PPIs) are potent inhibitors of gastric acid secretion, used as first-line agents in treating peptic ulcers. However, we have previously reported that PPIs may diminish the therapeutic effect of anti-vascular endothelial growth factor (VEGF) drugs in patients with cancer. In this study, we explored the effects of vonoprazan, a novel gastric acid secretion inhibitor used for the treatment of peptic ulcers, on the secretion of VEGF in cancer cells and attempted to propose it as an alternative PPI for cancer chemotherapy. The effects of PPI and vonoprazan on VEGF expression in cancer cells were compared by real-time reverse transcription-polymerase chain reaction and ELISA. The interaction of vonoprazan and PPIs with transcriptional regulators by docking simulation analysis. In various cancer cell lines, including the human colorectal cancer cell line (LS174T), PPI increased VEGF messenger RNA expression and VEGF protein secretion, while this effect was not observed with vonoprazan. Molecular docking simulation analysis showed that vonoprazan had a lower binding affinity for estrogen receptor alpha (ER-α), one of the transcriptional regulators of VEGF, compared to PPI. Although the PPI-induced increase in VEGF expression was counteracted by pharmacological ER-α inhibition, the effect of vonoprazan on VEGF expression was unchanged. Vonoprazan does not affect VEGF expression in cancer cells, which suggests that vonoprazan might be an alternative to PPIs, with no interference with the therapeutic effects of anti-VEGF cancer chemotherapy.

Regulome-based characterization of drug activity across the human diseasome.

Drugs are expected to recover the cell system away from the impaired state to normalcy through disease treatment. However, the understanding of gene regulatory machinery underlying drug activity or disease pathogenesis is far from complete. Here, we perform large-scale regulome analysis for various diseases in terms of gene regulatory machinery. Transcriptome signatures were converted into regulome signatures of transcription factors by integrating publicly available ChIP-seq data. Regulome-based correlations between diseases and their approved drugs were much clearer than the transcriptome-based correlations. For example, an inverse correlation was observed for cancers, whereas a positive correlation was observed for immune system diseases. After demonstrating the usefulness of the regulome-based drug discovery method in terms of accuracy and applicability, we predicted new drugs for nonsmall cell lung cancer and validated the anticancer activity in vitro. The proposed method is useful for understanding disease-disease relationships and drug discovery.

Small compound-based direct cell conversion with combinatorial optimization of pathway regulations.

MOTIVATION: Direct cell conversion, direct reprogramming (DR), is an innovative technology that directly converts source cells to target cells without bypassing induced pluripotent stem cells. The use of small compounds (e.g. drugs) for DR can help avoid carcinogenic risk induced by gene transfection; however, experimentally identifying small compounds remains challenging because of combinatorial explosion. RESULTS: In this article, we present a new computational method, COMPRENDRE (combinatorial optimization of pathway regulations for direct reprograming), to elucidate the mechanism of small compound-based DR and predict new combinations of small compounds for DR. We estimated the potential target proteins of DR-inducing small compounds and identified a set of target pathways involving DR. We identified multiple DR-related pathways that have not previously been reported to induce neurons or cardiomyocytes from fibroblasts. To overcome the problem of combinatorial explosion, we developed a variant of a simulated annealing algorithm to identify the best set of compounds that can regulate DR-related pathways. Consequently, the proposed method enabled to predict new DR-inducing candidate combinations with fewer compounds and to successfully reproduce experimentally verified compounds inducing the direct conversion from fibroblasts to neurons or cardiomyocytes. The proposed method is expected to be useful for practical applications in regenerative medicine. AVAILABILITY AND IMPLEMENTATION: The code supporting the current study is available at the http://labo.bio.kyutech.ac.jp/~yamani/comprendre. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

From drug repositioning to target repositioning: prediction of therapeutic targets using genetically perturbed transcriptomic signatures.

MOTIVATION: A critical element of drug development is the identification of therapeutic targets for diseases. However, the depletion of therapeutic targets is a serious problem. RESULTS: In this study, we propose the novel concept of target repositioning, an extension of the concept of drug repositioning, to predict new therapeutic targets for various diseases. Predictions were performed by a trans-disease analysis which integrated genetically perturbed transcriptomic signatures (knockdown of 4345 genes and overexpression of 3114 genes) and disease-specific gene transcriptomic signatures of 79 diseases. The trans-disease method, which takes into account similarities among diseases, enabled us to distinguish the inhibitory from activatory targets and to predict the therapeutic targetability of not only proteins with known target-disease associations but also orphan proteins without known associations. Our proposed method is expected to be useful for understanding the commonality of mechanisms among diseases and for therapeutic target identification in drug discovery. AVAILABILITY AND IMPLEMENTATION: Supplemental information and software are available at the following website [http://labo.bio.kyutech.ac.jp/~yamani/target_repositioning/]. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Corrigendum: A Web Server for GPCR-GPCR Interaction Pair Prediction.

[This corrects the article DOI: 10.3389/fendo.2022.825195.].