Drug repositioning by structure-based virtual screening

Repurposing of USFDA-approved drugs to identify leads for inhibition of acetylcholinesterase enzyme: a plausible utility as an anti-Alzheimer agent

In the quest to identify new anti-Alzheimer agents, we employed drug repositioning or drug repositioning techniques on approved USFDA small molecules. Herein, we report the structure-based virtual screening (SBVS) of 1880 USFDA-approved drugs. The in silico-based identification was followed by calculating Prime MMGB-SA binding energy and molecular dynamics simulation studies. The cumulative analysis led to identifying domperidone as an identified hit. Domperidone was further corroborated in vitro using anticholinesterase-based assessment, keeping donepezil as a positive control. The analysis revealed that the identified lead (domperidone) could induce an inhibitory effect on AChE in a dose-dependent manner with an IC50 of 3.67 μM as compared to donepezil, which exhibited an IC50 of 1.37 μM. However, as domperidone is known to have poor BBB permeability, we rationally proposed new analogues utilizing the principles of bioisosterism. The bioisostere-clubbed analogues were found to have better BBB permeability, affinity, and stability within the catalytic domain of AChE via molecular docking and dynamics studies. The proposed bioisosteres may be synthesized in the future. They may plausibly be explored for their implication in the developmental progress of new anti-Alzheimer agent achieved via repurposing techniques in future.

Drug repurposing for diabetes mellitus: In silico and in vitro investigation of DrugBank database for α-glucosidase inhibitors

The process of developing novel compounds/drugs is arduous, time-intensive, and financially burdensome, characterized by a notably low success rate and relatively high attrition rates. To alleviate these challenges, compound/drug repositioning strategies are employed to predict potential therapeutic effects for DrugBank-approved compounds across various diseases. In this study, we devised a computational and enzyme inhibitory mechanistic approach to identify promising compounds from the pool of DrugBank-approved substances targeting Diabetes Mellitus (DM). Molecular docking analyses were employed to validate the binding interaction patterns and conformations of the screened compounds within the active site of α-glucosidase. Notably, Asp352 and Glu277 participated in interactions within the α-glucosidase-ligand complexes, mediated by conventional hydrogen bonding and van der Waals forces, respectively. The stability of the docked complexes (α-glucosidase-compounds) was scrutinized through Molecular Dynamics (MD) simulations. Subsequent in vitro analyses assessed the therapeutic potential of the repositioned compounds against α-glucosidase. Kinetic studies revealed that "Forodesine" exhibited a lower IC50 (0.24 ± 0.04 mM) compared to the control, and its inhibitory pattern corresponds to that of competitive inhibitors. In-depth in silico secondary structure content analysis detailed the interactions between Forodesine and α-glucosidase, unveiling significant alterations in enzyme conformation upon binding, impacting its catalytic activity. Overall, our findings underscore the potential of Forodesine as a promising candidate for DM treatment through α-glucosidase inhibition. Further validation through in vitro and in vivo studies is imperative to confirm the therapeutic benefits of Forodesine in conformational diseases such as DM.

Drug repurposing for neurodegenerative diseases

Neurodegenerative diseases (NDDs) are neuronal problems that include the brain and spinal cord and result in loss of sensory and motor dysfunction. Common NDDs include Alzheimer's disease (AD), Parkinson's disease (PD), Huntington's disease (HD), Multiple Sclerosis (MS), and Amyotrophic Lateral Sclerosis (ALS) etc. The occurrence of these diseases increases with age and is one of the challenging problems among elderly people. Though, several scientific research has demonstrated the key pathologies associated with NDDs still the underlying mechanisms and molecular details are not well understood and need to be explored and this poses a lack of effective treatments for NDDs. Several lines of evidence have shown that NDDs have a high prevalence and affect more than a billion individuals globally but still, researchers need to work forward in identifying the best therapeutic target for NDDs. Thus, several researchers are working in the directions to find potential therapeutic targets to alter the disease pathology and treat the diseases. Several steps have been taken to identify the early detection of the disease and drug repurposing for effective treatment of NDDs. Moreover, it is logical that current medications are being evaluated for their efficacy in treating such disorders; therefore, drug repurposing would be an efficient, safe, and cost-effective way in finding out better medication. In the current manuscript we discussed the utilization of drugs that have been repurposed for the treatment of AD, PD, HD, MS, and ALS.

Identifying novel aryl hydrocarbon receptor (AhR) modulators from clinically approved drugs: In silico screening and In vitro validation

The aryl hydrocarbon receptor (AhR) functions as a vital ligand-activated transcription factor, governing both physiological and pathophysiological processes. Notably, it responds to xenobiotics, leading to a diverse array of outcomes. In the context of drug repurposing, we present here a combined approach of utilizing structure-based virtual screening and molecular dynamics simulations. This approach aims to identify potential AhR modulators from Drugbank repository of clinically approved drugs. By focusing on the AhR PAS-B binding pocket, our screening protocol included binding affinities calculations, complex stability, and interactions within the binding site as a filtering method. Comprehensive evaluations of all DrugBank small molecule database revealed ten promising hits. This included flibanserin, butoconazole, luliconazole, naftifine, triclabendazole, rosiglitazone, empagliflozin, benperidol, nebivolol, and zucapsaicin. Each exhibiting diverse binding behaviors and remarkably very low binding free energy. Experimental studies further illuminated their modulation of AhR signaling, and showing that they are consistently reducing AhR activity, except for luliconazole, which intriguingly enhances the AhR activity. This work demonstrates the possibility of using computational modelling as a quick screening tool to predict new AhR modulators from extensive drug libraries. Importantly, these findings hold immense therapeutic potential for addressing AhR-associated disorders. Consequently, it offers compelling prospects for innovative interventions through drug repurposing.

Benchmarking ANI potentials as a rescoring function and screening FDA drugs for SARS-CoV-2 Mpro

Here, we introduce the use of ANI-ML potentials as a rescoring function in the host-guest interaction in molecular docking. Our results show that the "docking power" of ANI potentials can compete with the current scoring functions at the same level of computational cost. Benchmarking studies on CASF-2016 dataset showed that ANI is ranked in the top 5 scoring functions among the other 34 tested. In particular, the ANI predicted interaction energies when used in conjunction with GOLD-PLP scoring function can boost the top ranked solution to be the closest to the x-ray structure. Rapid and accurate calculation of interaction energies between ligand and protein also enables screening of millions of drug candidates/docking poses. Using a unique protocol in which docking by GOLD-PLP, rescoring by ANI-ML potentials and extensive MD simulations along with end state free energy methods are combined, we have screened FDA approved drugs against the SARS-CoV-2 main protease (Mpro). The top six drug molecules suggested by the consensus of these free energy methods have already been in clinical trials or proposed as potential drug molecules in previous theoretical and experimental studies, approving the validity and the power of accuracy in our screening method.

Predicting drug-target interactions using matrix factorization with self-paced learning and dual similarity information

BACKGROUND

Drug repositioning (DR) refers to a method used to find new targets for existing drugs. This method can effectively reduce the development cost of drugs, save time on drug development, and reduce the risks of drug design. The traditional experimental methods related to DR are time-consuming, expensive, and have a high failure rate. Several computational methods have been developed with the increase in data volume and computing power. In the last decade, matrix factorization (MF) methods have been widely used in DR issues. However, these methods still have some challenges. (1) The model easily falls into a bad local optimal solution due to the high noise and high missing rate in the data. (2) Single similarity information makes the learning power of the model insufficient in terms of identifying the potential associations accurately.

OBJECTIVE

We proposed self-paced learning with dual similarity information and MF (SPLDMF), which introduced the self-paced learning method and more information related to drugs and targets into the model to improve prediction performance.

METHODS

Combining self-paced learning first can effectively alleviate the model prone to fall into a bad local optimal solution because of the high noise and high data missing rate. Then, we incorporated more data into the model to improve the model's capacity for learning.

RESULTS

Our model achieved the best results on each dataset tested. For example, the area under the receiver operating characteristic curve and the precision-recall curve of SPLDMF was 0.982 and 0.815, respectively, outperforming the state-of-the-art methods.

CONCLUSION

The experimental results on five benchmark datasets and two extended datasets demonstrated the effectiveness of our approach in predicting drug-target interactions.

MCL-DTI: using drug multimodal information and bi-directional cross-attention learning method for predicting drug-target interaction

BACKGROUND

Prediction of drug-target interaction (DTI) is an essential step for drug discovery and drug reposition. Traditional methods are mostly time-consuming and labor-intensive, and deep learning-based methods address these limitations and are applied to engineering. Most of the current deep learning methods employ representation learning of unimodal information such as SMILES sequences, molecular graphs, or molecular images of drugs. In addition, most methods focus on feature extraction from drug and target alone without fusion learning from drug-target interacting parties, which may lead to insufficient feature representation.

MOTIVATION

In order to capture more comprehensive drug features, we utilize both molecular image and chemical features of drugs. The image of the drug mainly has the structural information and spatial features of the drug, while the chemical information includes its functions and properties, which can complement each other, making drug representation more effective and complete. Meanwhile, to enhance the interactive feature learning of drug and target, we introduce a bidirectional multi-head attention mechanism to improve the performance of DTI.

RESULTS

To enhance feature learning between drugs and targets, we propose a novel model based on deep learning for DTI task called MCL-DTI which uses multimodal information of drug and learn the representation of drug-target interaction for drug-target prediction. In order to further explore a more comprehensive representation of drug features, this paper first exploits two multimodal information of drugs, molecular image and chemical text, to represent the drug. We also introduce to use bi-rectional multi-head corss attention (MCA) method to learn the interrelationships between drugs and targets. Thus, we build two decoders, which include an multi-head self attention (MSA) block and an MCA block, for cross-information learning. We use a decoder for the drug and target separately to obtain the interaction feature maps. Finally, we feed these feature maps generated by decoders into a fusion block for feature extraction and output the prediction results.

CONCLUSIONS

MCL-DTI achieves the best results in all the three datasets: Human, C. elegans and Davis, including the balanced datasets and an unbalanced dataset. The results on the drug-drug interaction (DDI) task show that MCL-DTI has a strong generalization capability and can be easily applied to other tasks.

Epithelial-mesenchymal transition-related gene signature for prognosis of lung squamous cell carcinoma

Epithelial-mesenchymal transition (EMT) is associated with tumor invasion and progression, and is regulated by DNA methylation. A prognostic signature of lung squamous cell carcinoma (LUSC) with EMT-related gene data has not yet been established. In our study, we constructed a co-expression network using differentially expressed genes (DEGs) obtained from The Cancer Genome Atlas (TCGA) to identify hub genes. We conducted a correlation analysis between the differentially methylated hub genes and differentially expressed EMT-related genes to screen EMT-related differentially methylated genes (ERDMGs). Functional enrichment was performed to annotate the ERDMGs. The least absolute shrinkage and selection operator (LASSO) and stepwise Cox regression analyses were performed to build a survival prognosis prediction model. Additionally, druggability analysis was performed to predict the potential drug targets of ERDMGs. We screened 11 ERDMGs that were enriched in cell adhesion molecules and other signaling pathways. Finally, we constructed a 4-ERDMG model, which showed good ability to predict survival prognosis in the training and validation sets. The model could serve as an independent predictive factor for patients with LUSC. Additionally, our druggability analysis predicted that CC chemokine ligand 23 (CCL23) and Hepatocyte nuclear factor 1b (HNF1B) may be the underlying drug targets of LUSC. We established a new risk score (RS) system as a prognostic indicator to predict the outcome of patients with LUSC, which will help in the improvement of treatment strategies.

LM-DTI: a tool of predicting drug-target interactions using the node2vec and network path score methods

Introduction: Drug-target interaction (DTI) prediction is a key step in drug function discovery and repositioning. The emergence of large-scale heterogeneous biological networks provides an opportunity to identify drug-related target genes, which led to the development of several computational methods for DTI prediction. Methods: Considering the limitations of conventional computational methods, a novel tool named LM-DTI based on integrated information related to lncRNAs and miRNAs was proposed, which adopted the graph embedding (node2vec) and the network path score methods. First, LM-DTI innovatively constructed a heterogeneous information network containing eight networks composed of four types of nodes (drug, target, lncRNA, and miRNA). Next, the node2vec method was used to obtain feature vectors of drug as well as target nodes, and the path score vector of each drug-target pair was calculated using the DASPfind method. Finally, the feature vectors and path score vectors were merged and input into the XGBoost classifier to predict potential drug-target interactions. Results and Discussion: The 10-fold cross validations evaluate the classification accuracies of the LM-DTI. The prediction performance of LM-DTI in AUPR reached 0.96, which showed a significant improvement compared with those of conventional tools. The validity of LM-DTI has also been verified by manually searching literature and various databases. LM-DTI is scalable and computing efficient; thus representing a powerful drug relocation tool that can be accessed for free at http://www.lirmed.com:5038/lm_dti.

Drug repositioning: doxazosin attenuates the virulence factors and biofilm formation in Gram-negative bacteria

The resistance development is an increasing global health risk that needs innovative solutions. Repurposing drugs to serve as anti-virulence agents is suggested as an advantageous strategy to diminish bacterial resistance development. Bacterial virulence is controlled by quorum sensing (QS) system that orchestrates the expression of biofilm formation, motility, and virulence factors production as enzymes and virulent pigments. Interfering with QS could lead to bacterial virulence mitigation without affecting bacterial growth that does not result in bacterial resistance development. This study investigated the probable anti-virulence and anti-QS activities of α-adrenoreceptor blocker doxazosin against Proteus mirabilis and Pseudomonas aeruginosa. Besides in silico study, in vitro and in vivo investigations were conducted to assess the doxazosin anti-virulence actions. Doxazosin significantly diminished the biofilm formation and release of QS-controlled Chromobacterium violaceum pigment and virulence factors in P. aeruginosa and P. mirabilis, and downregulated the QS encoding genes in P. aeruginosa. Virtually, doxazosin interfered with QS proteins, and in vivo protected mice against P. mirabilis and P. aeruginosa. The role of the membranal sensors as QseC and PmrA was recognized in enhancing the Gram-negative virulence. Doxazosin downregulated the membranal sensors PmR and QseC encoding genes and could in silico interfere with them. In conclusion, this study preliminary documents the probable anti-QS and anti-virulence activities of doxazosin, which indicate its possible application as an alternative or in addition to antibiotics. However, extended toxicological and pharmacological investigations are essential to approve the feasible clinical application of doxazosin as novel efficient anti-virulence agent. KEY POINTS: • Anti-hypertensive doxazosin acquires anti-quorum sensing activities • Doxazosin diminishes the virulence of Proteus mirabilis and Pseudomonas aeruginosa • Doxazosin could dimmish the bacterial espionage.

Pharmacological Chaperones and Protein Conformational Diseases: Approaches of Computational Structural Biology

Whenever a protein fails to fold into its native structure, a profound detrimental effect is likely to occur, and a disease is often developed. Protein conformational disorders arise when proteins adopt abnormal conformations due to a pathological gene variant that turns into gain/loss of function or improper localization/degradation. Pharmacological chaperones are small molecules restoring the correct folding of a protein suitable for treating conformational diseases. Small molecules like these bind poorly folded proteins similarly to physiological chaperones, bridging non-covalent interactions (hydrogen bonds, electrostatic interactions, and van der Waals contacts) loosened or lost due to mutations. Pharmacological chaperone development involves, among other things, structural biology investigation of the target protein and its misfolding and refolding. Such research can take advantage of computational methods at many stages. Here, we present an up-to-date review of the computational structural biology tools and approaches regarding protein stability evaluation, binding pocket discovery and druggability, drug repurposing, and virtual ligand screening. The tools are presented as organized in an ideal workflow oriented at pharmacological chaperones' rational design, also with the treatment of rare diseases in mind.

In silico screening of inhibitors against human dihydrofolate reductase to identify potential anticancer compounds

In all species, dihydrofolate reductase (DHFR) is an essential enzyme that regulates the cellular amount of tetrahydrofolate. Human DHFR (hDHFR) activity inhibition results in tetrahydrofolate depletion and cell death. This property has made hDHFR a therapeutic target for cancer. Methotrexate is a well-known hDHFR inhibitor, but its administration has shown some light to severe adverse effects. Therefore, we aimed to find new potential hDHFR inhibitors using structure-based virtual screening, ADMET prediction, molecular docking, and molecular dynamics simulations. Here, we used the PubChem database to find all compounds with at least 90% structural similarity to known natural DHFR inhibitors. To explore their interaction pattern and estimate their binding affinities, the screened compounds (2023) were subjected to structure-based molecular docking against hDHFR. The fifteen compounds that showed higher binding affinity to the hDHFR than the reference compound (methotrexate) displayed important molecular orientation and interactions with key residues in the enzyme's active site. These compounds were subjected to Lipinski and ADMET prediction. PubChem CIDs: 46886812 and 638190 were identified as putative inhibitors. In addition, molecular dynamics simulations revealed that the binding of compounds (CIDs: 46886812 and 63819) stabilized the hDHFR structure and caused minor conformational changes. Our findings suggest that two compounds (CIDs: 46886812 and 63819) could be promising potential inhibitors of hDHFR in cancer therapy.Communicated by Ramaswamy H. Sarma.

Structure based Drug Designing Approaches in SARS-CoV-2 Spike Inhibitor Design

The COVID-19 outbreak and the pandemic situation have hastened the research community to design a novel drug and vaccine against its causative organism, the SARS-CoV-2. The spike glycoprotein present on the surface of this pathogenic organism plays an immense role in viral entry and antigenicity. Hence, it is considered an important drug target in COVID-19 drug design. Several three-dimensional crystal structures of this SARS-CoV-2 spike protein have been identified and deposited in the Protein DataBank during the pandemic period. This accelerated the research in computer- aided drug designing, especially in the field of structure-based drug designing. This review summarizes various structure-based drug design approaches applied to this SARS-CoV-2 spike protein and its findings. Specifically, it is focused on different structure-based approaches such as molecular docking, high-throughput virtual screening, molecular dynamics simulation, drug repurposing, and target-based pharmacophore modelling and screening. These structural approaches have been applied to different ligands and datasets such as FDA-approved drugs, small molecular chemical compounds, chemical libraries, chemical databases, structural analogs, and natural compounds, which resulted in the prediction of spike inhibitors, spike-ACE-2 interface inhibitors, and allosteric inhibitors.

A Hybrid Docking and Machine Learning Approach to Enhance the Performance of Virtual Screening Carried out on Protein-Protein Interfaces

The modulation of protein-protein interactions (PPIs) by small chemical compounds is challenging. PPIs play a critical role in most cellular processes and are involved in numerous disease pathways. As such, novel strategies that assist the design of PPI inhibitors are of major importance. We previously reported that the knowledge-based DLIGAND2 scoring tool was the best-rescoring function for improving receptor-based virtual screening (VS) performed with the Surflex docking engine applied to several PPI targets with experimentally known active and inactive compounds. Here, we extend our investigation by assessing the vs. potential of other types of scoring functions with an emphasis on docking-pose derived solvent accessible surface area (SASA) descriptors, with or without the use of machine learning (ML) classifiers. First, we explored rescoring strategies of Surflex-generated docking poses with five GOLD scoring functions (GoldScore, ChemScore, ASP, ChemPLP, ChemScore with Receptor Depth Scaling) and with consensus scoring. The top-ranked poses were post-processed to derive a set of protein and ligand SASA descriptors in the bound and unbound states, which were combined to derive descriptors of the docked protein-ligand complexes. Further, eight ML models (tree, bagged forest, random forest, Bayesian, support vector machine, logistic regression, neural network, and neural network with bagging) were trained using the derivatized SASA descriptors and validated on test sets. The results show that many SASA descriptors are better than Surflex and GOLD scoring functions in terms of overall performance and early recovery success on the used dataset. The ML models were superior to all scoring functions and rescoring approaches for most targets yielding up to a seven-fold increase in enrichment factors at 1% of the screened collections. In particular, the neural networks and random forest-based ML emerged as the best techniques for this PPI dataset, making them robust and attractive vs. tools for hit-finding efforts. The presented results suggest that exploring further docking-pose derived SASA descriptors could be valuable for structure-based virtual screening projects, and in the present case, to assist the rational design of small-molecule PPI inhibitors.

Ibrahim T, Khafagy ES, Lopes BS, Ali MAM, Hegazy WAH, Elfaky MA. Characterization of the Anti-Biofilm and Anti-Quorum Sensing Activities of the β-Adrenoreceptor Antagonist Atenolol against Gram-Negative Bacterial Pathogens

The development of bacterial resistance to antibiotics is an increasing public health issue that worsens with the formation of biofilms. Quorum sensing (QS) orchestrates the bacterial virulence and controls the formation of biofilm. Targeting bacterial virulence is promising approach to overcome the resistance increment to antibiotics. In a previous detailed in silico study, the anti-QS activities of twenty-two β-adrenoreceptor blockers were screened supposing atenolol as a promising candidate. The current study aims to evaluate the anti-QS, anti-biofilm and anti-virulence activities of the β-adrenoreceptor blocker atenolol against Gram-negative bacteria Serratia marcescens, Pseudomonas aeruginosa, and Proteus mirabilis. An in silico study was conducted to evaluate the binding affinity of atenolol to S. marcescens SmaR QS receptor, P. aeruginosa QscR QS receptor, and P. mirabilis MrpH adhesin. The atenolol anti-virulence activity was evaluated against the tested strains in vitro and in vivo. The present finding shows considerable ability of atenolol to compete with QS proteins and significantly downregulated the expression of QS- and virulence-encoding genes. Atenolol showed significant reduction in the tested bacterial biofilm formation, virulence enzyme production, and motility. Furthermore, atenolol significantly diminished the bacterial capacity for killing and protected mice. In conclusion, atenolol has potential anti-QS and anti-virulence activities against S. marcescens, P. aeruginosa, and P. mirabilis and can be used as an adjuvant in treatment of aggressive bacterial infections.

Muting Bacterial Communication: Evaluation of Prazosin Anti-Quorum Sensing Activities against Gram-Negative Bacteria Pseudomonas aeruginosa, Proteus mirabilis, and Serratia marcescens

Simple Summary

Bacterial infections are considered one of the main challenges to global health. Bacterial virulence is controlled by interplayed systems to regulate bacterial invasion and infection in host tissues. Quorum sensing (QS) plays a crucial role in regulating virulence factor production, thus could be considered as the bacterial communication system in the bacterial population. The current study aimed to assess the anti-QS and anti-virulence activities of α-adrenoreceptor prazosin against three virulent Gram-negative bacteria. It was demonstrated that prazosin significantly downregulates the expression of QS-encoding genes and shows considered ability to compete on QS proteins in tested strains. Prazosin can significantly diminish biofilm formation and production of virulent enzymes and mitigate the virulence factors of tested strains. However, more testing is required alongside pharmacological and toxicological studies to assure the potential clinical use of prazosin as an adjuvant anti-QS and anti-virulence agent.

Abstract

Quorum sensing (QS) controls the production of several bacterial virulence factors. There is accumulative evidence to support that targeting QS can ensure a significant diminishing of bacterial virulence. Lessening bacterial virulence has been approved as an efficient strategy to overcome the development of antimicrobial resistance. The current study aimed to assess the anti-QS and anti-virulence activities of α-adrenoreceptor prazosin against three virulent Gram-negative bacteria Pseudomonades aeruginosa, Proteus mirabilis, and Serratia marcescens. The evaluation of anti-QS was carried out on a series of in vitro experiments, while the anti-virulence activities of prazosin were tested in an in vivo animal model. The prazosin anti-QS activity was assessed on the production of QS-controlled Chromobacterium violaceum pigment violacein and the expression of QS-encoding genes in P. aeruginosa. In vitro tests were performed to evaluate the prazosin effects on biofilm formation and production of extracellular enzymes by P. aeruginosa, P. mirabilis, and S. marcescens. A protective assay was conducted to evaluate the in vivo anti-virulence activity of prazosin against P. aeruginosa, P. mirabilis, and S. marcescens. Moreover, precise in silico molecular docking was performed to test the prazosin affinity to different QS receptors. The results revealed that prazosin significantly decreased the production of violacein and the virulent enzymes, protease and hemolysins, in the tested strains. Prazosin significantly diminished biofilm formation in vitro and bacterial virulence in vivo. The prazosin anti-QS activity was proven by its downregulation of QS-encoding genes and its obvious binding affinity to QS receptors. In conclusion, prazosin could be considered an efficient anti-virulence agent to be used as an adjuvant to antibiotics, however, it requires further pharmacological evaluations prior to clinical application.

Protocol for high-throughput screening of ACE2 enzymatic activators to treat COVID-19-induced metabolic complications

Drug repositioning represents a cost- and time-efficient strategy for drug development. Here, we present a workflow of in silico screening of ACE2 enzymatic activators to treat COVID-19-induced metabolic complications. By using structure-based virtual screening and signature-based off-target effect identification via the Connectivity Map database, we provide a ranked list of the repositioning candidates as potential ACE2 enzymatic activators to ameliorate COVID-19-induced metabolic complications. The workflow can also be applied to other diseases with ACE2 as a potential target.

For complete details on the use and execution of this protocol, please refer to Li et al. (2022).

•

Structure-based high-throughput virtual screening for ACE2 enzymatic activators

•

Signature-based drug repositioning for COVID-19-induced metabolic complications

•

Workflow applicable for other diseases with ACE2 as a potential target

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Recent advances of IDH1 mutant inhibitor in cancer therapy

Isocitrate dehydrogenase (IDH) is the key metabolic enzyme that catalyzes the conversion of isocitrate to α-ketoglutarate (α-KG). Two main types of IDH1 and IDH2 are present in humans. In recent years, mutations in IDH have been observed in several tumors, including glioma, acute myeloid leukemia, and chondrosarcoma. Among them, the frequency of IDH1 mutations is higher than IDH2. IDH1 mutations have been shown to increase the conversion of α-KG to 2-hydroxyglutarate (2-HG). IDH1 mutation-mediated accumulation of 2-HG leads to epigenetic dysregulation, altering gene expression, and impairing cell differentiation. A rapidly emerging therapeutic approach is through the development of small molecule inhibitors targeting mutant IDH1 (mIDH1), as evidenced by the recently approved of the first selective IDH1 mutant inhibitor AG-120 (ivosidenib) for the treatment of IDH1-mutated AML. This review will focus on mIDH1 as a therapeutic target and provide an update on IDH1 mutant inhibitors in development and clinical trials.

Protein-Ligand Docking in the Machine-Learning Era

Molecular docking plays a significant role in early-stage drug discovery, from structure-based virtual screening (VS) to hit-to-lead optimization, and its capability and predictive power is critically dependent on the protein-ligand scoring function. In this review, we give a broad overview of recent scoring function development, as well as the docking-based applications in drug discovery. We outline the strategies and resources available for structure-based VS and discuss the assessment and development of classical and machine learning protein-ligand scoring functions. In particular, we highlight the recent progress of machine learning scoring function ranging from descriptor-based models to deep learning approaches. We also discuss the general workflow and docking protocols of structure-based VS, such as structure preparation, binding site detection, docking strategies, and post-docking filter/re-scoring, as well as a case study on the large-scale docking-based VS test on the LIT-PCBA data set.

Insights into the Cardiotoxic Effects of Veratrum Lobelianum Alkaloids: Pilot Study

Jervine, protoveratrine A (proA), and protoveratrine B (proB) are Veratrum alkaloids that are presented in some remedies obtained from Veratrum lobelianum, such as Veratrum aqua. This paper reports on a single-center pilot cardiotoxic mechanism study of jervine, proA, and proB in case series. The molecular aspects were studied via molecular dynamic simulation, molecular docking with cardiac sodium channel Na_V1.5, and machine learning-based structure–activity relationship modeling. HPLC-MS/MS method in combination with clinical events were used to analyze Veratrum alkaloid cardiotoxicity in patients. Jervine demonstrates the highest docking score (−10.8 kcal/mol), logP value (4.188), and pK_a value (9.64) compared with proA and proB. Also, this compound is characterized by the lowest calculated IC₅₀. In general, all three analyzed alkaloids show the affinity to Na_V1.5 that highly likely results in cardiotoxic action. The clinical data of seven cases of intoxication by Veratrum aqua confirms the results of molecular modeling. Patients exhibited nausea, muscle weakness, bradycardia, and arterial hypotension. The association between alkaloid concentrations in blood and urine and severity of patient condition is described. These experiments, while primary, confirmed that jervine, proA, and proB contribute to cardiotoxicity by Na_V1.5 inhibition.

Screening of potent STAT3-SH2 domain inhibitors from JAK/STAT compound library through molecular dynamics simulation

The Signal Transducer and Activator of Transcription 3 (STAT3) protein is activated consistently in the tumor cells and thus studied as a potent target for cancer prevention. The TYR705-phosphorylated (pTyr) STAT3 forms a homo-dimer by binding to its recognition site in the Src Homology 2 (SH2) domain of another STAT3 monomer, causing cellular survival, proliferation, inflammation, and tumor invasion. Many inhibitors of STAT3-SH2 have recently been identified using both computational and experimental approaches. In this study, we used molecular docking, Absorption, Distribution, Metabolism, and Excretion/Toxicological (ADME/tox) and molecular dynamics modeling to examine binding affinities and specificities of 191 inhibitor drugs from the SELLECKCHEM database. The binding free energies of the inhibitors were calculated by Induced Fit Docking (IFD) prime energy. The binding hotspots of STAT3-SH2 were evaluated via binding energy decomposition and hydrogen bond distribution analysis, and the inhibitor compound's stability was assessed through MD simulation. (-)-Epigallocatechin gallate, Kaempferol-3-O-rutinoside, Picroside I, Saikosaponin D, and Ginsenoside Rk1 were found to be the top hit inhibitor compounds. They exhibited an exceptional docking score, a low binding free energy, interacted with the key amino acid residue, and showed significant ADME/tox moderation. These compounds were further proved to be favorable by their stability in an MD simulation run for 100 ns using GROMACS software. The inhibitors (-)-Epigallocatechin gallate, Kaempferol-3-O-rutinoside, and Saikosaponin D show improved stability in molecular dynamic modeling and are expected to have a significant STAT3-SH2 inhibitory effect against cancer.

Structure-Based Virtual Screening Identifies Novobiocin, Montelukast, and Cinnarizine as TRPV1 Modulators with Anticonvulsant Activity In Vivo

The transient receptor potential vanilloid 1 (TRPV1) receptor is a nonselective cation channel, known to be involved in the regulation of many important physiological and pathological processes. In the last few years, it has been proposed as a promising target to develop novel anticonvulsant compounds. However, thermoregulatory effects associated with the channel inhibition have hampered the path for TRPV1 antagonists to become marketed drugs. In this regard, we conducted a structure-based virtual screening campaign to find potential TRPV1 modulators among approved drugs, which are known to be safe and thermally neutral. To this end, different docking models were developed and validated by assessing their pose and score prediction powers. Novobiocin, montelukast, and cinnarizine were selected from the screening as promising candidates for experimental testing and all of them exhibited nanomolar inhibitory activity. Moreover, the in vivo profiles showed promising results in at least one of the three models of seizures tested.

Natural Compound ZINC12899676 Reduces Porcine Epidemic Diarrhea Virus Replication by Inhibiting the Viral NTPase Activity

Porcine epidemic diarrhea virus (PEDV) is an alphacoronavirus (α-CoV) that causes high mortality in suckling piglets, leading to severe economic losses worldwide. No effective vaccine or commercial antiviral drug is readily available. Several replicative enzymes are responsible for coronavirus replication. In this study, the potential candidates targeting replicative enzymes (PLP2, 3CLpro, RdRp, NTPase, and NendoU) were screened from 187,119 compounds in ZINC natural products library, and seven compounds had high binding potential to NTPase and showed drug-like property. Among them, ZINC12899676 was identified to significantly inhibit the NTPase activity of PEDV by targeting its active pocket and causing its conformational change, and ZINC12899676 significantly inhibited PEDV replication in IPEC-J2 cells. It first demonstrated that ZINC12899676 inhibits PEDV replication by targeting NTPase, and then, NTPase may serve as a novel target for anti-PEDV.

Structure-based drug repurposing: traditional and advanced AI/ML-aided methods

The current global health emergency in the form of the Coronavirus 2019 (COVID-19) pandemic has highlighted the need for fast, accurate, and efficient drug discovery pipelines. Traditional drug discovery projects relying on in vitro high-throughput screening (HTS) involve large investments and sophisticated experimental set-ups, affordable only to big biopharmaceutical companies. In this scenario, application of efficient state-of-the-art computational methods and modern artificial intelligence (AI)-based algorithms for rapid screening of repurposable chemical space [approved drugs and natural products (NPs) with proven pharmacokinetic profiles] to identify the initial leads is a powerful option to save resources and time. Structure-based drug repurposing is a popular in silico repurposing approach. In this review, we discuss traditional and modern AI-based computational methods and tools applied at various stages for structure-based drug discovery (SBDD) pipelines. Additionally, we highlight the role of generative models in generating molecules with scaffolds from repurposable chemical space. Teaser: This review highlights the importance of repurposable chemical space, and the contributions of conventional in silico approaches and modern machine-learning algorithms for rapid structure-based drug repurposing.

Data Mining Meets Machine Learning: A Novel ANN-based Multi-Body Interaction Docking Scoring Function (MBI-Score) based on Utilizing Frequent Geometric and Chemical Patterns of Interfacial Atoms in Native Protein-Ligand Complexes

Accurate prediction of binding poses is crucial to structure-based drug design. We employ two powerful artificial intelligence (AI) approaches, data-mining and machine-learning, to design artificial neural network (ANN) based pose-scoring function. It is a simple machine-learning-based statistical function that employs frequent geometric and chemical patterns of interacting atoms at protein-ligand interfaces. The patterns are derived by mining interfaces of "native" protein-ligand complexes. Each interface is represented by a graph where nodes are atoms and edges connect protein-ligand interfacial atoms located within certain cutoff distance of each other. Applying frequent subgraph mining to these interfaces provides "native" frequent patterns of interacting atoms. Subsequently, given a pose for a protein-ligand complex of interest, the pose-scoring function (the information-processing unit or neuron) calculates the degree of matching between the interaction patterns present at the pose's interface and the native frequent patterns. The pose-scoring function takes into account the frequency of occurrence of the matching native patterns, the size of the match, and the degree of geometrical similarity between pose-specific and matching native frequent patterns. This novel "multi-body interaction" pose-scoring function (MBI-Score) was validated using two databases, PDBbind and Astex-85, and it outperformed seven commonly used commercial scoring functions. MBI-Score is available at www.khashanlab.org/mbi-score.

Screening of Clinically Approved and Investigation Drugs as Potential Inhibitors of SARS-CoV-2: A Combined in silico and in vitro Study

In the current study, we used 7922 FDA approved small molecule drugs as well as compounds in clinical investigation from NIH's NPC database in our drug repurposing study. SARS-CoV-2 main protease as well as Spike protein/ACE2 targets were used in virtual screening and top-100 compounds from each docking simulations were considered initially in short molecular dynamics (MD) simulations and their average binding energies were calculated by MM/GBSA method. Promising hit compounds selected based on average MM/GBSA scores were then used in long MD simulations. Based on these numerical calculations following compounds were found as hit inhibitors for the SARS-CoV-2 main protease: Pinokalant, terlakiren, ritonavir, cefotiam, telinavir, rotigaptide, and cefpiramide. In addition, following 3 compounds were identified as inhibitors for Spike/ACE2: Denopamine, bometolol, and rotigaptide. In order to verify the predictions of in silico analyses, 4 compounds (ritonavir, rotigaptide, cefotiam, and cefpiramide) for the main protease and 2 compounds (rotigaptide and denopamine) for the Spike/ACE2 interactions were tested by in vitro experiments. While the concentration-dependent inhibition of the ritonavir, rotigaptide, and cefotiam was observed for the main protease; denopamine was effective at the inhibition of Spike/ACE2 binding.

Pharmacological inhibition of KDM5A for cancer treatment

Lysine-specific demethylase 5A (KDM5A, also named RBP2 or JARID1A) is a demethylase that can remove methyl groups from histones H3K4me1/2/3. It is aberrantly expressed in many cancers, where it impedes differentiation and contributes to cancer cell proliferation, cell metastasis and invasiveness, drug resistance, and is associated with poor prognosis. Pharmacological inhibition of KDM5A has been reported to significantly attenuate tumor progression in vitro and in vivo in a range of solid tumors and acute myeloid leukemia. This review will present the structural aspects of KDM5A, its role in carcinogenesis, a comparison of currently available approaches for screening KDM5A inhibitors, a classification of KDM5A inhibitors, and its potential as a drug target in cancer therapy.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

Comparative analysis of network-based approaches and machine learning algorithms for predicting drug-target interactions

Computational prediction of drug-target interactions (DTIs) is of particular importance in the process of drug repositioning because of its efficiency in selecting potential candidates for DTIs. A variety of computational methods for predicting DTIs have been proposed over the past decade. Our interest is which methods or techniques are the most advantageous for increasing prediction accuracy. This article provides a comprehensive overview of network-based, machine learning, and integrated DTI prediction methods. The network-based methods handle a DTI network along with drug and target similarities in a matrix form and apply graph-theoretic algorithms to identify new DTIs. Machine learning methods use known DTIs and the features of drugs and target proteins as training data to build a predictive model. Integrated methods combine these two techniques. We assessed the prediction performance of the selected state-of-the-art methods using two different benchmark datasets. Our experimental results demonstrate that the integrated methods outperform the others in general. Some previous methods showed low accuracy on predicting interactions of unknown drugs which do not exist in the training dataset. Combining similarity matrices from multiple features by data fusion was not beneficial in increasing prediction accuracy. Finally, we analyzed future directions for further improvements in DTI predictions.

WaterMap and Molecular Dynamic Simulation-Guided Discovery of Potential PAK1 Inhibitors Using Repurposing Approaches

p21-Activated kinase 1 (PAK1) is positioned at the nexus of several oncogenic signaling pathways. Currently, there are no approved inhibitors for disabling the transfer of phosphate in the active site directly, as they are limited by lower affinity, and poor kinase selectivity. In this work, a repurposing study utilizing FDA-approved drugs from the DrugBank database was pursued with an initial selection of 27 molecules out of ∼2162 drug molecules, based on their docking energies and molecular interaction patterns. From the molecules that were considered for WaterMap analysis, seven molecules, namely, Mitoxantrone, Labetalol, Acalabrutinib, Sacubitril, Flubendazole, Trazodone, and Niraparib, ascertained the ability to overlap with high-energy hydration sites. Considering many other displaced unfavorable water molecules, only Acalabrutinib, Flubendazole, and Trazodone molecules highlighted their prominence in terms of binding affinity gains through ΔΔG that ranges between 6.44 and 2.59 kcal/mol. Even if Mitoxantrone exhibited the highest docking score and greater interaction strength, it did not comply with the WaterMap and molecular dynamics simulation results. Moreover, detailed MD simulation trajectory analyses suggested that the drug molecules Flubendazole, Niraparib, and Acalabrutinib were highly stable, observed from their RMSD values and consistent interaction pattern with Glu315, Glu345, Leu347, and Asp407 including the hydrophobic interactions maintained in the three replicates. However, the drug molecule Trazodone displayed a loss of crucial interaction with Leu347, which was essential to inhibit the kinase activity of PAK1. The molecular orbital and electrostatic potential analyses elucidated the reactivity and strong complementarity potentials of the drug molecules in the binding pocket of PAK1. Therefore, the CADD-based reposition efforts, reported in this work, helped in the successful identification of new PAK1 inhibitors that requires further investigation by in vitro analysis.

Current status of structure-based drug repurposing against COVID-19 by targeting SARS-CoV-2 proteins

More than one and half years have passed, as of August 2021, since the COVID-19 caused by the novel coronavirus named SARS-CoV-2 emerged in 2019. While the recent success of vaccine developments likely reduces the severe cases, there is still a strong requirement of safety and effective therapeutic drugs for overcoming the unprecedented situation. Here we review the recent progress and the status of the drug discovery against COVID-19 with emphasizing a structure-based perspective. Structural data regarding the SARS-CoV-2 proteome has been rapidly accumulated in the Protein Data Bank, and up to 68% of the total amino acid residues encoded in the genome were covered by the structural data. Despite a global effort of in silico and in vitro screenings for drug repurposing, there is only a limited number of drugs had been successfully authorized by drug regulation organizations. Although many approved drugs and natural compounds, which exhibited antiviral activity in vitro, were considered potential drugs against COVID-19, a further multidisciplinary investigation is required for understanding the mechanisms underlying the antiviral effects of the drugs.

DTi2Vec: Drug-target interaction prediction using network embedding and ensemble learning

Drug-target interaction (DTI) prediction is a crucial step in drug discovery and repositioning as it reduces experimental validation costs if done right. Thus, developing in-silico methods to predict potential DTI has become a competitive research niche, with one of its main focuses being improving the prediction accuracy. Using machine learning (ML) models for this task, specifically network-based approaches, is effective and has shown great advantages over the other computational methods. However, ML model development involves upstream hand-crafted feature extraction and other processes that impact prediction accuracy. Thus, network-based representation learning techniques that provide automated feature extraction combined with traditional ML classifiers dealing with downstream link prediction tasks may be better-suited paradigms. Here, we present such a method, DTi2Vec, which identifies DTIs using network representation learning and ensemble learning techniques. DTi2Vec constructs the heterogeneous network, and then it automatically generates features for each drug and target using the nodes embedding technique. DTi2Vec demonstrated its ability in drug-target link prediction compared to several state-of-the-art network-based methods, using four benchmark datasets and large-scale data compiled from DrugBank. DTi2Vec showed a statistically significant increase in the prediction performances in terms of AUPR. We verified the "novel" predicted DTIs using several databases and scientific literature. DTi2Vec is a simple yet effective method that provides high DTI prediction performance while being scalable and efficient in computation, translating into a powerful drug repositioning tool.

Prediction of Drug-Target Interactions by Combining Dual-Tree Complex Wavelet Transform with Ensemble Learning Method

Identification of drug–target interactions (DTIs) is vital for drug discovery. However, traditional biological approaches have some unavoidable shortcomings, such as being time consuming and expensive. Therefore, there is an urgent need to develop novel and effective computational methods to predict DTIs in order to shorten the development cycles of new drugs. In this study, we present a novel computational approach to identify DTIs, which uses protein sequence information and the dual-tree complex wavelet transform (DTCWT). More specifically, a position-specific scoring matrix (PSSM) was performed on the target protein sequence to obtain its evolutionary information. Then, DTCWT was used to extract representative features from the PSSM, which were then combined with the drug fingerprint features to form the feature descriptors. Finally, these descriptors were sent to the Rotation Forest (RoF) model for classification. A 5-fold cross validation (CV) was adopted on four datasets (Enzyme, Ion Channel, GPCRs (G-protein-coupled receptors), and NRs (Nuclear Receptors)) to validate the proposed model; our method yielded high average accuracies of 89.21%, 85.49%, 81.02%, and 74.44%, respectively. To further verify the performance of our model, we compared the RoF classifier with two state-of-the-art algorithms: the support vector machine (SVM) and the k-nearest neighbor (KNN) classifier. We also compared it with some other published methods. Moreover, the prediction results for the independent dataset further indicated that our method is effective for predicting potential DTIs. Thus, we believe that our method is suitable for facilitating drug discovery and development.

A data-driven methodology towards evaluating the potential of drug repurposing hypotheses

Drug repurposing has become a widely used strategy to accelerate the process of finding treatments. While classical de novo drug development involves high costs, risks, and time-consuming paths, drug repurposing allows to reuse already-existing and approved drugs for new indications. Numerous research has been carried out in this field, both in vitro and in silico. Computational drug repurposing methods make use of modern heterogeneous biomedical data to identify and prioritize new indications for old drugs. In the current paper, we present a new complete methodology to evaluate new potentially repurposable drugs based on disease-gene and disease-phenotype associations, identifying significant differences between repurposing and non-repurposing data. We have collected a set of known successful drug repurposing case studies from the literature and we have analysed their dissimilarities with other biomedical data not necessarily participating in repurposing processes. The information used has been obtained from the DISNET platform. We have performed three analyses (at the genetical, phenotypical, and categorization levels), to conclude that there is a statistically significant difference between actual repurposing-related information and non-repurposing data. The insights obtained could be relevant when suggesting new potential drug repurposing hypotheses.

Targeting SARS-COV-2 non-structural protein 16: a virtual drug repurposing study

Non-Structural Protein 16 (nsp-16), a viral RNA methyltransferase (MTase), is one of the highly viable targets for drug discovery of coronaviruses including SARS-CoV-2. In this study, drug discovery of SARS-CoV-2 nsp-16 has been performed by a virtual drug repurposing approach. First, drug shape-based screening (among FDA approved drugs) with a known template of MTase inhibitor, sinefungin was done and best compounds with high similarity scores were selected. In addition to the selected compounds, 4 nucleoside analogs of anti-viral (Raltgravir, Maraviroc and Favipiravir) and anti-inflammatory (Prednisolone) drugs were selected for further investigations. Then, binding energies and interaction modes were found by molecular docking approaches and compouds with lower energy were selected for further investigation. After that, Molecular dynamics (MD) simulation was carried to test the potential selected compounds in a realistic environment. The results showed that Raltegravir and Maraviroc among other compounds can bind strongly to the active site of the protein compared to sinefungin, and can be potential candidates to inhibit NSP-16. Also, the MD simulation results suggested that the Maraviroc and Raltegravir are more effective drug candidates than Sinefungin for inhibiting the enzyme. It is concluded that Raltegravir and Maraviroc which may be used in the treatment of COVID-19 after Invitro and invivo studies and clinical trial for final confirmation of drug effectiveness. Communicated by Ramaswamy H. Sarma.

MCN-CPI: Multiscale Convolutional Network for Compound-Protein Interaction Prediction

In the process of drug discovery, identifying the interaction between the protein and the novel compound plays an important role. With the development of technology, deep learning methods have shown excellent performance in various situations. However, the compound-protein interaction is complicated and the features extracted by most deep models are not comprehensive, which limits the performance to a certain extent. In this paper, we proposed a multiscale convolutional network that extracted the local and global features of the protein and the topological feature of the compound using different types of convolutional networks. The results showed that our model obtained the best performance compared with the existing deep learning methods.

Ubiquitination Regulators Discovered by Virtual Screening for the Treatment of Cancer

The ubiquitin-proteasome system oversees cellular protein degradation in order to regulate various critical processes, such as cell cycle control and DNA repair. Ubiquitination can serve as a marker for mutation, chemical damage, transcriptional or translational errors, and heat-induced denaturation. However, aberrant ubiquitination and degradation of tumor suppressor proteins may result in the growth and metastasis of cancer. Hence, targeting the ubiquitination cascade reaction has become a potential strategy for treating malignant diseases. Meanwhile, computer-aided methods have become widely accepted as fast and efficient techniques for early stage drug discovery. This review summarizes ubiquitination regulators that have been discovered via virtual screening and their applications for cancer treatment.

Identification of TNIK as a novel potential drug target in thyroid cancer based on protein druggability prediction

Thyroid cancer is a common endocrine malignancy; however, surgery remains its primary treatment option. A novel targeted drug for the development and application of targeted therapy in thyroid cancer treatment remain underexplored.We obtained RNA sequence data of thyroid cancer from The Cancer Genome Atlas database and identified differentially expressed genes (DEGs). Then, we constructed co-expression network with DEGs and combined it with differentially methylation analysis to screen the key genes in thyroid cancer. PockDrug-Server, an online tool, was applied to predict the druggability of the key genes. Finally, we constructed protein-protein interaction (PPI) network to observe potential targeted drugs for thyroid cancer.We identified 3 genes correlated with altered DNA methylation level and oncogenesis of thyroid cancer. According to the druggable analysis and PPI network, we predicted TRAF2 and NCK-interacting protein kinase (TNIK) sever as the drug targeted for thyroid cancer. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis indicated that genes in protein-protein interaction network of TNIK enriched in mitogen-activated protein kinase signaling pathway. For drug repositioning, we identified a targeted drug of genes in PPI network.Our study provides a bioinformatics method for screening drug targets and provides a theoretical basis for thyroid cancer targeted therapy.

Multi-PLI: interpretable multi-task deep learning model for unifying protein-ligand interaction datasets

The assessment of protein–ligand interactions is critical at early stage of drug discovery. Computational approaches for efficiently predicting such interactions facilitate drug development. Recently, methods based on deep learning, including structure- and sequence-based models, have achieved impressive performance on several different datasets. However, their application still suffers from a generalizability issue because of insufficient data, especially for structure based models, as well as a heterogeneity problem because of different label measurements and varying proteins across datasets. Here, we present an interpretable multi-task model to evaluate protein–ligand interaction (Multi-PLI). The model can run classification (binding or not) and regression (binding affinity) tasks concurrently by unifying different datasets. The model outperforms traditional docking and machine learning on both binary classification and regression tasks and achieves competitive results compared with some structure-based deep learning methods, even with the same training set size. Furthermore, combined with the proposed occlusion algorithm, the model can predict the important amino acids of proteins that are crucial for binding, thus providing a biological interpretation.

Evaluation of the effects of chlorhexidine and several flavonoids as antiviral purposes on SARS-CoV-2 main protease: molecular docking, molecular dynamics simulation studies

The recent outbreak of COVID-19 caused by a new human coronavirus called SARS-CoV-2, is continually causing worldwide human infections and deaths.The main protease (3CLpro), which plays a critical role in the life cycle of the virus, makes it an attractive target for the development of antiviral agents effective against coronaviruses (CoVs).Currently, there is no specific viral protein targeted therapeutics.Therefore, there is a need to investigate an alternative therapy which will prevent the spread of the infection, by focusing on the transmission of the virus.Chlorhexidine (CHX) and flavonoids agents have shown that they have a viral inactivation effect against enveloped viruses, and thus facilitate the struggle against oral transmission.Especially, some flavonoids have very strong antiviral activity in SARS-CoV and MERS-CoV main protease.This study was conducted to evaluate the CHX and flavonoids compounds potential antiviral effects on SARS-CoV-2 main protease through virtual screening for the COVID-19 treatment by molecular docking method.According to the results of this study, CHX, Kaempferol-3-rutinoside, Rutin, Quercetin 3-beta-D-glucoside and Isobavachalcone exhibited the best binding affinity against this enzyme, and also these compounds showed significant inhibitory effects compared to the SARS-CoV-2 main protease crystal structure inhibitor (N3).Especially, these compounds mainly interact with His41, Cys145, His163, Met165, Glu166 and Thr190 in SARS-CoV-2 main protease binding site. Further, MD simulation analysis also confirmed that stability of these interactions between the enzyme and these five compounds.The current study provides to guide clinical trials for broad-spectrum CHX and bioactive flavonoids to reduce the viral load of the infection and possibly disease progression.Communicated by Ramaswamy H. Sarma.

Novel Development of Predictive Feature Fingerprints to Identify Chemistry-Based Features for the Effective Drug Design of SARS-CoV-2 Target Antagonists and Inhibitors Using Machine Learning

A unique approach to bioactivity and chemical data curation coupled with random forest analyses has led to a series of target-specific and cross-validated predictive feature fingerprints (PFF) that have high predictability across multiple therapeutic targets and disease stages involved in the severe acute respiratory syndrome due to coronavirus 2 (SARS-CoV-2)-induced COVID-19 pandemic, which include plasma kallikrein, human immunodeficiency virus (HIV)-protease, nonstructural protein (NSP)5, NSP12, Janus kinase (JAK) family, and AT-1. The approach was highly accurate in determining the matched target for the different compound sets and suggests that the models could be used for virtual screening of target-specific compound libraries. The curation-modeling process was successfully applied to a SARS-CoV-2 phenotypic screen and could be used for predictive bioactivity estimation and prioritization for clinical trial selection; virtual screening of drug libraries for the repurposing of drug molecules; and analysis and direction of proprietary data sets.

Potential repurposing of four FDA approved compounds with antiplasmodial activity identified through proteome scale computational drug discovery and in vitro assay

Malaria elimination can benefit from time and cost-efficient approaches for antimalarials such as drug repurposing. In this work, 796 DrugBank compounds were screened against 36 Plasmodium falciparum targets using QuickVina-W. Hits were selected after rescoring using GRaph Interaction Matching (GRIM) and ligand efficiency metrics: surface efficiency index (SEI), binding efficiency index (BEI) and lipophilic efficiency (LipE). They were further evaluated in Molecular dynamics (MD). Twenty-five protein-ligand complexes were finally retained from the 28,656 (36 × 796) dockings. Hit GRIM scores (0.58 to 0.78) showed their molecular interaction similarity to co-crystallized ligands. Minimum LipE (3), SEI (23) and BEI (7) were in at least acceptable thresholds for hits. Binding energies ranged from -6 to -11 kcal/mol. Ligands showed stability in MD simulation with good hydrogen bonding and favorable protein-ligand interactions energy (the poorest being -140.12 kcal/mol). In vitro testing showed 4 active compounds with two having IC50 values in the single-digit μM range.

Engineering of fluorescent or photoactive Trojan probes for detection and eradication of β-Amyloids

Trojan horse technology institutes a potentially promising strategy to bring together a diagnostic or cell-based drug design and a delivery platform. It provides the opportunity to re-engineer a novel multimodal, neurovascular detection probe, or medicine to fuse with blood-brain barrier (BBB) molecular Trojan horse. In Alzheimer's disease (AD) this could allow the targeted delivery of detection or therapeutic probes across the BBB to the sites of plaques and tangles development to image or decrease amyloid load, enhance perivascular Aβ clearance, and improve cerebral blood flow, owing principally to the significantly improved cerebral permeation. A Trojan horse can also be equipped with photosensitizers, nanoparticles, quantum dots, or fluorescent molecules to function as multiple targeting theranostic compounds that could be activated following changes in disease-specific processes of the diseased tissue such as pH and protease activity, or exogenous stimuli such as, light. This concept review theorizes the use of receptor-mediated transport-based platforms to transform such novel ideas to engineer systemic and smart Trojan detection or therapeutic probes to advance the neurodegenerative field.

Harnessing endophenotypes and network medicine for Alzheimer's drug repurposing

Following two decades of more than 400 clinical trials centered on the "one drug, one target, one disease" paradigm, there is still no effective disease-modifying therapy for Alzheimer's disease (AD). The inherent complexity of AD may challenge this reductionist strategy. Recent observations and advances in network medicine further indicate that AD likely shares common underlying mechanisms and intermediate pathophenotypes, or endophenotypes, with other diseases. In this review, we consider AD pathobiology, disease comorbidity, pleiotropy, and therapeutic development, and construct relevant endophenotype networks to guide future therapeutic development. Specifically, we discuss six main endophenotype hypotheses in AD: amyloidosis, tauopathy, neuroinflammation, mitochondrial dysfunction, vascular dysfunction, and lysosomal dysfunction. We further consider how this endophenotype network framework can provide advances in computational and experimental strategies for drug-repurposing and identification of new candidate therapeutic strategies for patients suffering from or at risk for AD. We highlight new opportunities for endophenotype-informed, drug discovery in AD, by exploiting multi-omics data. Integration of genomics, transcriptomics, radiomics, pharmacogenomics, and interactomics (protein-protein interactions) are essential for successful drug discovery. We describe experimental technologies for AD drug discovery including human induced pluripotent stem cells, transgenic mouse/rat models, and population-based retrospective case-control studies that may be integrated with multi-omics in a network medicine methodology. In summary, endophenotype-based network medicine methodologies will promote AD therapeutic development that will optimize the usefulness of available data and support deep phenotyping of the patient heterogeneity for personalized medicine in AD.

Fenoterol and dobutamine as SARS-CoV-2 main protease inhibitors: A virtual screening study

Global health is under heavy threat by a worldwide pandemic caused by a new type of coronavirus (COVID-19) since its rapid spread in China in 2019 [1]. Currently, there are no approved specific drugs and effective treatment for COVID-19 infection, but several available drugs are known to facilitate tentative treatment. Since drug design, development and testing procedures are time-consuming [2], [1], [2], [3], virtual screening studies with the aid of available drug databases take the initiative at this point and save the time. Besides, drug repurposing strategies promises to identify new agents for the novel diseases in a time-critical fashion. In this study, we used structure based virtual screening method on FDA approved drugs and compounds in clinical trials. As a result of this study we choose three most prominent compounds for further studies. Here we show that these three compounds (dobutamine and its two derivatives) can be considered as promising inhibitors for SARS-CoV-2 main protease and results also demonstrate the possible interactions of dobutamine and its derivatives with SARS-CoV-2 main protease (6W63) [6]. Our efforts in this work directly address current urgency of a new drug discovery against COVID-19.

Ensemble Learning Prediction of Drug-Target Interactions Using GIST Descriptor Extracted from PSSM-Based Evolutionary Information

Identifying the drug-target interactions (DTIs) plays an essential role in new drug development. However, there still has the limited knowledge of DTIs and a significant number of unknown DTI pairs. Moreover, the traditional experimental methods have inevitable disadvantages such as high cost and time-consuming. Therefore, developing computational methods for predicting DTIs is attracting more and more attention. In this study, we report a novel computational approach for predicting DTI using GIST feature, position-specific scoring matrix (PSSM), and rotation forest (RF). Specifically, each target protein is first converted into a PSSM for retaining evolutionary information. Then, the GIST feature is extracted from PSSM and substructure fingerprint information is adopted to extract the feature of the drug. Finally, combining each protein and drug features to form a new drug-target pair, which is employed as input feature for RF classifier. In the experiment, the proposed method achieves high average accuracies of 89.25%, 85.93%, 82.36%, and 73.89% on enzyme, ion channel, G protein-coupled receptors (GPCRs), and nuclear receptor, respectively. For further evaluating the prediction performance of the proposed method, we compare it with the state-of-the-art support vector machine (SVM) classifier on the same golden standard dataset. These promising results illustrate that the proposed method is more effective and stable than other methods. We expect the proposed method to be a useful tool for predicting large-scale DTIs.

Fast Rescoring Protocols to Improve the Performance of Structure-Based Virtual Screening Performed on Protein-Protein Interfaces

Protein-protein interactions (PPIs) are attractive targets for drug design because of their essential role in numerous cellular processes and disease pathways. However, in general, PPIs display exposed binding pockets at the interface, and as such, have been largely unexploited for therapeutic interventions with low-molecular weight compounds. Here, we used docking and various rescoring strategies in an attempt to recover PPI inhibitors from a set of active and inactive molecules for 11 targets collected in ChEMBL and PubChem. Our focus is on the screening power of the various developed protocols and on using fast approaches so as to be able to apply such a strategy to the screening of ultralarge libraries in the future. First, we docked compounds into each target using the fast "pscreen" mode of the structure-based virtual screening (VS) package Surflex. Subsequently, the docking poses were postprocessed to derive a set of 3D topological descriptors: (i) shape similarity and (ii) interaction fingerprint similarity with a co-crystallized inhibitor, (iii) solvent-accessible surface area, and (iv) extent of deviation from the geometric center of a reference inhibitor. The derivatized descriptors, together with descriptor-scaled scoring functions, were utilized to investigate possible impacts on VS performance metrics. Moreover, four standalone scoring functions, RF-Score-VS (machine-learning), DLIGAND2 (knowledge-based), Vinardo (empirical), and X-SCORE (empirical), were employed to rescore the PPI compounds. Collectively, the results indicate that the topological scoring algorithms could be valuable both at a global level, with up to 79% increase in areas under the receiver operating characteristic curve for some targets, and in early stages, with up to a 4-fold increase in enrichment factors at 1% of the screened collections. Outstandingly, DLIGAND2 emerged as the best scoring function on this data set, outperforming all rescoring techniques in terms of VS metrics. The described methodology could help in the rational design of small-molecule PPI inhibitors and has direct applications in many therapeutic areas, including cancer, CNS, and infectious diseases such as COVID-19.

Overcoming cancer therapeutic bottleneck by drug repurposing

Ever present hurdles for the discovery of new drugs for cancer therapy have necessitated the development of the alternative strategy of drug repurposing, the development of old drugs for new therapeutic purposes. This strategy with a cost-effective way offers a rare opportunity for the treatment of human neoplastic disease, facilitating rapid clinical translation. With an increased understanding of the hallmarks of cancer and the development of various data-driven approaches, drug repurposing further promotes the holistic productivity of drug discovery and reasonably focuses on target-defined antineoplastic compounds. The "treasure trove" of non-oncology drugs should not be ignored since they could target not only known but also hitherto unknown vulnerabilities of cancer. Indeed, different from targeted drugs, these old generic drugs, usually used in a multi-target strategy may bring benefit to patients. In this review, aiming to demonstrate the full potential of drug repurposing, we present various promising repurposed non-oncology drugs for clinical cancer management and classify these candidates into their proposed administration for either mono- or drug combination therapy. We also summarize approaches used for drug repurposing and discuss the main barriers to its uptake.