Computational Multitarget Drug Design

Computational prediction of the molecular mechanism of statin group of drugs against SARS-CoV-2 pathogenesis

Recently published clinical data from COVID-19 patients indicated that statin therapy is associated with a better clinical outcome and a significant reduction in the risk of mortality. In this study by computational analysis, we have aimed to predict the possible mechanism of the statin group of drugs by which they can inhibit SARS-CoV-2 pathogenesis. Blind docking of the critical structural and functional proteins of SARS-CoV-2 like RNA-dependent RNA polymerase, M-protease of 3-CL-Pro, Helicase, and the Spike proteins ( wild type and mutants from different VOCs) were performed using the Schrodinger docking tool. We observed that fluvastatin and pitavastatin showed fair, binding affinities to RNA polymerase and 3-CL-Pro, whereas fluvastatin showed the strongest binding affinity to the helicase. Fluvastatin also showed the highest affinity for the SpikeDelta and a fair docking score for other spike variants. Additionally, molecular dynamics simulation confirmed the formation of a stable drug-protein complex between Fluvastatin and target proteins. Thus our study shows that of all the statins, fluvastatin can bind to multiple target proteins of SARS-CoV-2, including the spike-mutant proteins. This property might contribute to the potent antiviral efficacy of this drug.

Differentiating Inhibitors of Closely Related Protein Kinases with Single- or Multi-Target Activity via Explainable Machine Learning and Feature Analysis

Protein kinases are major drug targets. Most kinase inhibitors are directed against the adenosine triphosphate (ATP) cofactor binding site, which is largely conserved across the human kinome. Hence, such kinase inhibitors are often thought to be promiscuous. However, experimental evidence and activity data for publicly available kinase inhibitors indicate that this is not generally the case. We have investigated whether inhibitors of closely related human kinases with single- or multi-kinase activity can be differentiated on the basis of chemical structure. Therefore, a test system consisting of two distinct kinase triplets has been devised for which inhibitors with reported triple-kinase activities and corresponding single-kinase activities were assembled. Machine learning models derived on the basis of chemical structure distinguished between these multi- and single-kinase inhibitors with high accuracy. A model-independent explanatory approach was applied to identify structural features determining accurate predictions. For both kinase triplets, the analysis revealed decisive features contained in multi-kinase inhibitors. These features were found to be absent in corresponding single-kinase inhibitors, thus providing a rationale for successful machine learning. Mapping of features determining accurate predictions revealed that they formed coherent and chemically meaningful substructures that were characteristic of multi-kinase inhibitors compared with single-kinase inhibitors.

3β-Acetoxy-21α-H-hop-22(29)ene, a novel multitargeted phytocompound against SARS-CoV-2: in silico screening

The present pandemic disease COVID-19 demands an urgent need for more efficient antiviral drugs against SARS-CoV-2. Computational drug designing and discovery enable us to explore ethnomedicinal plants as a source of various lead molecules that can be used against present and future pathogens. Adiantum latifolium Lam., a common fern, is resistant to pathogens mainly due to the presence of various phytochemicals having antimicrobial properties. In our previous study, 3β-acetoxy-21α-H-hop-22(29)ene, a terpenoid has been characterized from the methanol extract of leaves of A. latifolium. The manuscript evaluates the antiviral potency of the compound against SARS-CoV-2 through molecular docking method. Proteins essential for SARS-CoV-2 multiplication in host cells are the target sites. The study revealed strong binding affinity of the compound for all the ten proteins selected, including seven nonstructural proteins, two structural proteins and one receptor protein, with a binding energy of -4.67 to -8.76 kcal/mol. MDS and MMPBSA analysis of the best ranked complex further confirmed the results. The multitargeted compound can be considered as a natural lead molecule in drug designing against COVID-19, but requires wet-lab experimentation and clinical trials.Communicated by Ramaswamy H. Sarma.

Development of Dual Inhibitors Targeting Epidermal Growth Factor Receptor in Cancer Therapy

Epidermal growth factor receptor (EGFR) is of great significance in mediating cell signaling transduction and tumor behaviors. Currently, third-generation inhibitors of EGFR, especially osimertinib, are at the clinical frontier for the treatment of EGFR-mutant non-small-cell lung cancer (NSCLC). Regrettably, the rapidly developing drug resistance caused by EGFR mutations and the compensatory mechanism have largely limited their clinical efficacy. Given the synergistic effect between EGFR and other compensatory targets during tumorigenesis and tumor development, EGFR dual-target inhibitors are promising for their reduced risk of drug resistance, higher efficacy, lower dosage, and fewer adverse events than those of single-target inhibitors. Hence, we present the synergistic mechanism underlying the role of EGFR dual-target inhibitors against drug resistance, their structure-activity relationships, and their therapeutic potential. Most importantly, we emphasize the optimal target combinations and design strategies for EGFR dual-target inhibitors and provide some perspectives on new challenges and future directions in this field.

Multi-Phase In Silico Discovery of Potential SARS-CoV-2 RNA-Dependent RNA Polymerase Inhibitors among 3009 Clinical and FDA-Approved Related Drugs

Proceeding our prior studies of SARS-CoV-2, the inhibitory potential against SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) has been investigated for a collection of 3009 clinical and FDA-approved drugs. A multi-phase in silico approach has been employed in this study. Initially, a molecular fingerprint experiment of Remdesivir (RTP), the co-crystallized ligand of the examined protein, revealed the most similar 150 compounds. Among them, 30 compounds were selected after a structure similarity experiment. Subsequently, the most similar 30 compounds were docked against SARS-CoV-2 RNA-dependent RNA polymerase (PDB ID: 7BV2). Aloin 359, Baicalin 456, Cefadroxil 1273, Sophoricoside 1459, Hyperoside 2109, and Vitexin 2286 exhibited the most precise binding modes, as well as the best binding energies. To confirm the obtained results, MD simulations experiments have been conducted for Hyperoside 2109, the natural flavonoid glycoside that exhibited the best docking scores, against RdRp (PDB ID: 7BV2) for 100 ns. The achieved results authenticated the correct binding of 2109, showing low energy and optimum dynamics. Our team presents these outcomes for scientists all over the world to advance in vitro and in vivo examinations against COVID-19 for the promising compounds.

Candidate Therapeutics by Screening for Multitargeting Ligands: Combining the CB2 Receptor With CB1, PPARγ and 5-HT4 Receptors

In recent years, the cannabinoid type 2 receptor (CB2R) has become a major target for treating many disease conditions. The old therapeutic paradigm of “one disease-one target-one drug” is being transformed to “complex disease-many targets-one drug.” Multitargeting, therefore, attracts much attention as a promising approach. We thus focus on designing single multitargeting agents (MTAs), which have many advantages over combined therapies. Using our ligand-based approach, the “Iterative Stochastic Elimination” (ISE) algorithm, we produce activity models of agonists and antagonists for desired therapeutic targets and anti-targets. These models are used for sequential virtual screening and scoring large libraries of molecules in order to pick top-scored candidates for testing in vitro and in vivo. In this study, we built activity models for CB2R and other targets for combinations that could be used for several indications. Those additional targets are the cannabinoid 1 receptor (CB1R), peroxisome proliferator-activated receptor gamma (PPARγ), and 5-Hydroxytryptamine receptor 4 (5-HT4R). All these models have high statistical parameters and are reliable. Many more CB2R/CBIR agonists were found than combined CB2R agonists with CB1R antagonist activity (by 200 fold). CB2R agonism combined with PPARγ or 5-HT4R agonist activity may be used for treating Inflammatory Bowel Disease (IBD). Combining CB2R agonism with 5-HT4R generates more candidates (14,008) than combining CB2R agonism with agonists for the nuclear receptor PPARγ (374 candidates) from an initial set of ∼2.1 million molecules. Improved enrichment of true vs. false positives may be achieved by requiring a better ISE score cutoff or by performing docking. Those candidates can be purchased and tested experimentally to validate their activity. Further, we performed docking to CB2R structures and found lower statistical performance of the docking (“structure-based”) compared to ISE modeling (“ligand-based”). Therefore, ISE modeling may be a better starting point for molecular discovery than docking.

Structure-based pharmacophore mapping and virtual screening of natural products to identify polypharmacological inhibitor against c-MET/EGFR/VEGFR-2

Three receptor tyrosine kinases (RTKs), c-MET, EGFR, and VEGFR-2 have been identified as potential oncogenic targets involved in tumor development, metastasis, and invasion. Designing inhibitors that can simultaneously interact with multiple targets is a promising approach, therefore, inhibiting these three RTKs with a single chemical component might give an effective chemotherapeutic strategy for addressing the disease while limiting adverse effects. The in-silico methods have been developed to identify the polypharmacological inhibitors particularly for drug repurposing and multitarget drug design. Here, to find a viable inhibitor from natural source against these three RTKs, structure-based pharmacophore mapping and virtual screening of SN-II database were carried out. The filtered compound SN00020821, identified as Cedeodarin, from different computational approaches, demonstrated good interactions with all the three targets, c-MET/EGFR/VEGFR-2, with interaction energies of -42.35 kcal/mol, -49.32 kcal/mol and -44.83 kcal/mol, respectively. SN00020821displayed stable key interactions with critical amino acids of all the three receptors' kinase catalytic domains including "DFG motif" explored through the MD simulations. Furthermore, it also met the ADMET requirements and was determined to be drug-like as predicted from the Lipinski's rule of five and Veber's rule. Finally, SN00020821 provides a novel molecular scaffold that could be investigated further as a polypharmacological anticancer therapeutic candidate that targets the three RTKs.Communicated by Ramaswamy H. Sarma.

Natural product drug discovery in the artificial intelligence era

Natural products (NPs) are primarily recognized as privileged structures to interact with protein drug targets. Their unique characteristics and structural diversity continue to marvel scientists for developing NP-inspired medicines, even though the pharmaceutical industry has largely given up. High-performance computer hardware, extensive storage, accessible software and affordable online education have democratized the use of artificial intelligence (AI) in many sectors and research areas. The last decades have introduced natural language processing and machine learning algorithms, two subfields of AI, to tackle NP drug discovery challenges and open up opportunities. In this article, we review and discuss the rational applications of AI approaches developed to assist in discovering bioactive NPs and capturing the molecular "patterns" of these privileged structures for combinatorial design or target selectivity.

Polypharmacological Approaches for CNS Diseases: Focus on Endocannabinoid Degradation Inhibition

Polypharmacology breaks up the classical paradigm of “one-drug, one target, one disease” electing multitarget compounds as potential therapeutic tools suitable for the treatment of complex diseases, such as metabolic syndrome, psychiatric or degenerative central nervous system (CNS) disorders, and cancer. These diseases often require a combination therapy which may result in positive but also negative synergistic effects. The endocannabinoid system (ECS) is emerging as a particularly attractive therapeutic target in CNS disorders and neurodegenerative diseases including Parkinson’s disease (PD), Alzheimer’s disease (AD), Huntington’s disease (HD), multiple sclerosis (MS), amyotrophic lateral sclerosis (ALS), stroke, traumatic brain injury (TBI), pain, and epilepsy. ECS is an organized neuromodulatory network, composed by endogenous cannabinoids, cannabinoid receptors type 1 and type 2 (CB₁ and CB₂), and the main catabolic enzymes involved in the endocannabinoid inactivation such as fatty acid amide hydrolase (FAAH) and monoacylglycerol lipase (MAGL). The multiple connections of the ECS with other signaling pathways in the CNS allows the consideration of the ECS as an optimal source of inspiration in the development of innovative polypharmacological compounds. In this review, we focused our attention on the reported polypharmacological examples in which FAAH and MAGL inhibitors are involved.

Alvarez KL, Aguilar-Pineda JA, Vera-Lopez KJ, Lino Cardenas CL. In Silico Analysis of Metabolites from Peruvian Native Plants as Potential Therapeutics against Alzheimer's Disease

Background: Despite research on the molecular bases of Alzheimer’s disease (AD), effective therapies against its progression are still needed. Recent studies have shown direct links between AD progression and neurovascular dysfunction, highlighting it as a potential target for new therapeutics development. In this work, we screened and evaluated the inhibitory effect of natural compounds from native Peruvian plants against tau protein, amyloid beta, and angiotensin II type 1 receptor (AT1R) pathologic AD markers. Methods: We applied in silico analysis, such as virtual screening, molecular docking, molecular dynamics simulation (MD), and MM/GBSA estimation, to identify metabolites from Peruvian plants with inhibitory properties, and compared them to nicotinamide, telmisartan, and grapeseed extract drugs in clinical trials. Results: Our results demonstrated the increased bioactivity of three plants’ metabolites against tau protein, amyloid beta, and AT1R. The MD simulations indicated the stability of the AT1R:floribundic acid, amyloid beta:rutin, and tau:brassicasterol systems. A polypharmaceutical potential was observed for rutin due to its high affinity to AT1R, amyloid beta, and tau. The metabolite floribundic acid showed bioactivity against the AT1R and tau, and the metabolite brassicasterol showed bioactivity against the amyloid beta and tau. Conclusions: This study has identified molecules from native Peruvian plants that have the potential to bind three pathologic markers of AD.

An in silico pipeline for the discovery of multitarget ligands: A case study for epi-polypharmacology based on DNMT1/HDAC2 inhibition

The search for novel therapeutic compounds remains an overwhelming task owing to the time-consuming and expensive nature of the drug development process and low success rates. Traditional methodologies that rely on the one drug-one target paradigm have proven insufficient for the treatment of multifactorial diseases, leading to a shift to multitarget approaches. In this emerging paradigm, molecules with off-target and promiscuous interactions may result in preferred therapies. In this study, we developed a general pipeline combining machine learning algorithms and a deep generator network to train a dual inhibitor classifier capable of identifying putative pharmacophoric traits. As a case study, we focused on dual inhibitors targeting DNA methyltransferase 1 (DNMT) and histone deacetylase 2 (HDAC2), two enzymes that play a central role in epigenetic regulation. We used this approach to identify dual inhibitors from a novel large natural product database in the public domain. We used docking and atomistic simulations as complementary approaches to establish the ligand-interaction profiles between the best hits and DNMT1/HDAC2. By using the combined ligand- and structure-based approaches, we discovered two promising novel scaffolds that can be used to simultaneously target both DNMT1 and HDAC2. We conclude that the flexibility and adaptability of the proposed pipeline has predictive capabilities of similar or derivative methods and is readily applicable to the discovery of small molecules targeting many other therapeutically relevant proteins.

Molecular Modeling Techniques Applied to the Design of Multitarget Drugs: Methods and Applications

Multifactorial diseases, such as cancer and diabetes present a challenge for the traditional "one-target, one disease" paradigm due to their complex pathogenic mechanisms. Although a combination of drugs can be used, a multitarget drug may be a better choice face of its efficacy, lower adverse effects and lower chance of resistance development. The computer-based design of these multitarget drugs can explore the same techniques used for single-target drug design, but the difficulties associated to the obtention of drugs that are capable of modulating two or more targets with similar efficacy impose new challenges, whose solutions involve the adaptation of known techniques and also to the development of new ones, including machine-learning approaches. In this review, some SBDD and LBDD techniques for the multitarget drug design are discussed, together with some cases where the application of such techniques led to effective multitarget ligands.

Explainable machine learning predictions of dual-target compounds reveal characteristic structural features

Compounds with defined multi-target activity play an increasingly important role in drug discovery. Structural features that might be signatures of such compounds have mostly remained elusive thus far. We have explored the potential of explainable machine learning to uncover structural motifs that are characteristic of dual-target compounds. For a pharmacologically relevant target pair-based test system designed for our study, accurate prediction models were derived and the influence of molecular representation features of test compounds was quantified to explain the predictions. The analysis revealed small numbers of specific features whose presence in dual-target and absence in single-target compounds determined accurate predictions. These features formed coherent substructures in dual-target compounds. From computational analysis of specific feature contributions, structural motifs emerged that were confirmed to be signatures of different dual-target activities. Our findings demonstrate the ability of explainable machine learning to bridge between predictions and intuitive chemical analysis and reveal characteristic substructures of dual-target compounds.

In Silico Screening of Semi-Synthesized Compounds as Potential Inhibitors for SARS-CoV-2 Papain-like Protease: Pharmacophoric Features, Molecular Docking, ADMET, Toxicity and DFT Studies

Papain-like protease is an essential enzyme in the proteolytic processing required for the replication of SARS-CoV-2. Accordingly, such an enzyme is an important target for the development of anti-SARS-CoV-2 agents which may reduce the mortality associated with outbreaks of SARS-CoV-2. A set of 69 semi-synthesized molecules that exhibited the structural features of SARS-CoV-2 papain-like protease inhibitors (PLPI) were docked against the coronavirus papain-like protease (PLpro) enzyme (PDB ID: (4OW0). Docking studies showed that derivatives 34 and 58 were better than the co-crystallized ligand while derivatives 17, 28, 31, 40, 41, 43, 47, 54, and 65 exhibited good binding modes and binding free energies. The pharmacokinetic profiling study was conducted according to the four principles of the Lipinski rules and excluded derivative 31. Furthermore, ADMET and toxicity studies showed that derivatives 28, 34, and 47 have the potential to be drugs and have been demonstrated as safe when assessed via seven toxicity models. Finally, comparing the molecular orbital energies and the molecular electrostatic potential maps of 28, 34, and 47 against the co-crystallized ligand in a DFT study indicated that 28 is the most promising candidate to interact with the target receptor (PLpro).

In Silico Exploration of Potential Natural Inhibitors against SARS-Cov-2 nsp10

In continuation of our previous effort, different in silico selection methods were applied to 310 naturally isolated metabolites that exhibited antiviral potentialities before. The applied selection methods aimed to pick the most relevant inhibitor of SARS-CoV-2 nsp10. At first, a structural similarity study against the co-crystallized ligand, S-Adenosyl Methionine (SAM), of SARS-CoV-2 nonstructural protein (nsp10) (PDB ID: 6W4H) was carried out. The similarity analysis culled 30 candidates. Secondly, a fingerprint study against SAM preferred compounds 44, 48, 85, 102, 105, 182, 220, 221, 282, 284, 285, 301, and 302. The docking studies picked 48, 182, 220, 221, and 284. While the ADMET analysis expected the likeness of the five candidates to be drugs, the toxicity study preferred compounds 48 and 182. Finally, a density-functional theory (DFT) study suggested vidarabine (182) to be the most relevant SARS-Cov-2 nsp10 inhibitor.

22-Hydroxyhopane, a novel multitargeted phytocompound against SARS-CoV-2 from Adiantum latifolium Lam

The present pandemic disease COVID-19 demands an urgent need for more efficient antiviral drugs against SARS-CoV-2. 22-Hydroxyhopane is a bioactive triterpenoid compound with antibacterial activity, present in the leaves of Adiantum latifolium. In this study, molecular docking method revealed strong binding affinity of the compound for ten proteins essential for SARS-CoV-2 multiplication in host cells, including seven nonstructural proteins, two structural proteins and one receptor protein, with a binding energy of -7.61 to -9.82 kcal/mol and inhibition constant <1 μM. MDS and MM-PBSA analysis of the best ranked complex further confirmed the results. The targets selected include six enzymes, RNA binding protein, spike protein, membrane protein and ACE2 receptor of SARS-CoV-2. It is the first report of a natural compound from A. latifolium having multitargeted activity against SARS-CoV-2. We conclude that 22-hydroxyhopane may be used as a best source for the development of novel therapeutic drugs for COVID-19, but requires further evaluations.

Exploration of potential inhibitors for tuberculosis via structure-based drug design, molecular docking, and molecular dynamics simulation studies

Drug resistance in tuberculosis is major threat to human population. In the present investigation, we aimed to identify novel and potent benzimidazole molecules to overcome the resistance management. A series of 20 benzimidazole derivatives were examined for its activity as selective antitubercular agents. Initially, AutodockVina algorithm was performed to assess the efficacy of the molecules. The results are further enriched by redocking by means of Glide algorithm. The binding free energies of the compounds were then calculated by MM-generalized-born surface area method. Molecular docking studies elucidated that benzimidazole derivatives has revealed formation of hydrogen bond and strong binding affinity in the active site of Mycobacterium tuberculosis protein. Note that ARG308, GLY189, VAL312, LEU403, and LEU190 amino acid residues of Mycobacterium tuberculosis protein PrpR are involved in binding with ligands of benzimidazoles. Interestingly, the ligands exhibited same binding potential to the active site of protein complex PrpR in both the docking programs. In essence, the result portrays that benzimidazole derivatives such as 1p, 1q, and 1 t could be potent and selective antitubercular agents than the standard drug isoniazid. These compounds were then subjected to molecular dynamics simulation to validate the dynamics activity of the compounds against PrpR. Finally, the inhibitory behavior of compounds was predicted using a machine learning algorithm trained on a data collection of 15,000 compounds utilizing graph-based signatures. Overall, the study concludes that designed benzimidazoles can be employed as antitubercular agents. Indeed, the results are helpful for the experimental biologists to develop safe and non-toxic drugs against tuberculosis.

In Silico-Based Design and In Vivo Evaluation of an Anthranilic Acid Derivative as a Multitarget Drug in a Diet-Induced Metabolic Syndrome Model

Metabolic syndrome (MetS) is a complex disease that affects almost a quarter of the world's adult population. In MetS, diabetes, obesity, hyperglycemia, high cholesterol, and high blood pressure are the most common disorders. Polypharmacy is the most used strategy for managing conditions related to MetS, but it has drawbacks such as low medication adherence. Multitarget ligands have been proposed as an interesting approach to developing drugs to treat complex diseases. However, suitable preclinical models that allow their evaluation in a context closer to a clinical situation of a complex disease are needed. From molecular docking studies, compound 1b, a 5-aminoanthranilic acid derivative substituted with 4'-trifluoromethylbenzylamino and 3',4'-dimethoxybenzamide moieties, was identified as a potential multitarget drug, as it showed high in silico affinity against targets related to MetS, including PPAR-α, PPAR-γ, and HMG-CoA reductase. It was evaluated in a diet-induced MetS rat model and simultaneously lowered blood pressure, glucose, total cholesterol, and triglyceride levels after a 14-day treatment. No toxicity events were observed during an acute lethal dose evaluation test at 1500 mg/kg. Hence, the diet-induced MetS model is suitable for evaluating treatments for MetS, and compound 1b is an attractive starting point for developing multitarget drugs.

Artificial intelligence to deep learning: machine intelligence approach for drug discovery

Drug designing and development is an important area of research for pharmaceutical companies and chemical scientists. However, low efficacy, off-target delivery, time consumption, and high cost impose a hurdle and challenges that impact drug design and discovery. Further, complex and big data from genomics, proteomics, microarray data, and clinical trials also impose an obstacle in the drug discovery pipeline. Artificial intelligence and machine learning technology play a crucial role in drug discovery and development. In other words, artificial neural networks and deep learning algorithms have modernized the area. Machine learning and deep learning algorithms have been implemented in several drug discovery processes such as peptide synthesis, structure-based virtual screening, ligand-based virtual screening, toxicity prediction, drug monitoring and release, pharmacophore modeling, quantitative structure-activity relationship, drug repositioning, polypharmacology, and physiochemical activity. Evidence from the past strengthens the implementation of artificial intelligence and deep learning in this field. Moreover, novel data mining, curation, and management techniques provided critical support to recently developed modeling algorithms. In summary, artificial intelligence and deep learning advancements provide an excellent opportunity for rational drug design and discovery process, which will eventually impact mankind. The primary concern associated with drug design and development is time consumption and production cost. Further, inefficiency, inaccurate target delivery, and inappropriate dosage are other hurdles that inhibit the process of drug delivery and development. With advancements in technology, computer-aided drug design integrating artificial intelligence algorithms can eliminate the challenges and hurdles of traditional drug design and development. Artificial intelligence is referred to as superset comprising machine learning, whereas machine learning comprises supervised learning, unsupervised learning, and reinforcement learning. Further, deep learning, a subset of machine learning, has been extensively implemented in drug design and development. The artificial neural network, deep neural network, support vector machines, classification and regression, generative adversarial networks, symbolic learning, and meta-learning are examples of the algorithms applied to the drug design and discovery process. Artificial intelligence has been applied to different areas of drug design and development process, such as from peptide synthesis to molecule design, virtual screening to molecular docking, quantitative structure-activity relationship to drug repositioning, protein misfolding to protein-protein interactions, and molecular pathway identification to polypharmacology. Artificial intelligence principles have been applied to the classification of active and inactive, monitoring drug release, pre-clinical and clinical development, primary and secondary drug screening, biomarker development, pharmaceutical manufacturing, bioactivity identification and physiochemical properties, prediction of toxicity, and identification of mode of action.

Recent Progress on Tubulin Inhibitors with Dual Targeting Capabilities for Cancer Therapy

Microtubules play a crucial role in multiple cellular functions including mitosis, cell signaling, and organelle trafficking, which makes the microtubule an important target for cancer therapy. Despite the great successes of microtubule-targeting agents in the clinic, the development of drug resistance and dose-limiting toxicity restrict their clinical efficacy. In recent years, multitarget therapy has been considered an effective strategy to achieve higher therapeutic efficacy, in particular dual-target drugs. In terms of the synergetic effect of tubulin and other antitumor agents such as receptor tyrosine kinases inhibitors, histone deacetylases inhibitors, DNA-damaging agents, and topoisomerase inhibitors in combination therapy, designing dual-target tubulin inhibitors is regarded as a promising approach to overcome these limitations and improve therapeutic efficacy. In this Perspective, we discussed rational target combinations, design strategies, structure-activity relationships, and future directions of dual-target tubulin inhibitors.

Fine-tuning of a generative neural network for designing multi-target compounds

Exploring the origin of multi-target activity of small molecules and designing new multi-target compounds are highly topical issues in pharmaceutical research. We have investigated the ability of a generative neural network to create multi-target compounds. Data sets of experimentally confirmed multi-target, single-target, and consistently inactive compounds were extracted from public screening data considering positive and negative assay results. These data sets were used to fine-tune the REINVENT generative model via transfer learning to systematically recognize multi-target compounds, distinguish them from single-target or inactive compounds, and construct new multi-target compounds. During fine-tuning, the model showed a clear tendency to increasingly generate multi-target compounds and structural analogs. Our findings indicate that generative models can be adopted for de novo multi-target compound design.

Compound dataset and custom code for deep generative multi-target compound design

AIM

Generating a data and software infrastructure for evaluating multi-target compound (MT-CPD) design via deep generative modeling.

METHODOLOGY

The REINVENT 2.0 approach for generative modeling was extended for MT-CPD design and a large benchmark data set was curated.

EXEMPLARY RESULTS & DATA

Proof-of-concept for deep generative MT-CPD design was established. Custom code and the benchmark set comprising 2809 MT-CPDs, 61,928 single-target and 295,395 inactive compounds from biological screens are made freely available.

LIMITATIONS & NEXT STEPS

MT-CPD design via deep learning is still at its conceptual stages. It will be required to demonstrate experimental impact. The data and software we provide enable further investigation of MT-CPD design and generation of candidate molecules for experimental programs.

Introducing the metacore concept for multi-target ligand design

In this work, we introduce the concept of "metacores" (MCs) for the organization of analog series (ASs) and multi-target (MT) ligand design. Generating compounds that are active against distantly related or unrelated targets is a central task in polypharmacology-oriented drug discovery. MCs are obtained by two-stage extraction of structural cores from ASs. The methodology is chemically intuitive and generally applicable. Each MC represents a set of related ASs and a template for the generation of new structures. We have systematically identified ASs that exclusively consisted of analogs with MT activity and determined their target profiles. From these ASs, a large set of 317 structurally diverse MCs was extracted, 127 of which were associated with different target families. The newly generated MCs were characterized and further prioritized on the basis of AS, compound, and target coverage. The analysis indicated that 260 MCs were pharmaceutically relevant. These MCs and the compound and target information they capture are made freely available for medicinal chemistry applications.

Multiple Target Drug Design Using LigBuilder 3

Designing drugs that directly interact with multiple targets is a promising approach for treating complicated diseases. In order to successfully bind to multiple targets of different families and achieve the desired ligand efficiency, multi-target-directed ligands (MTDLs) require a higher level of diversity and complexity. De novo design strategies for creating more diverse chemical entities with desired properties may present an improved approach for developing MTDLs. In this chapter, we describe a computational protocol for developing MTDLs using the first reported multi-target de novo program, LigBuilder 3, which combines a binding site prediction module with de novo drug design and optimization modules. As an illustration of each detailed procedure, we design dual-functional compounds of two well-characterized virus enzymes, HIV protease and reverse transcriptase (PR and RT, respectively), using fragments extracted from known inhibitors. LigBuilder 3 is accessible at http://www.pkumdl.cn/ligbuilder3/ .

Machine learning reveals that structural features distinguishing promiscuous and non-promiscuous compounds depend on target combinations

Compounds with defined multi-target activity (promiscuity) play an increasingly important role in drug discovery. However, the molecular basis of multi-target activity is currently only little understood. In particular, it remains unclear whether structural features exist that generally characterize promiscuous compounds and set them apart from compounds with single-target activity. We have devised a test system using machine learning to systematically examine structural features that might characterize compounds with multi-target activity. Using this system, more than 860,000 diagnostic predictions were carried out. The analysis provided compelling evidence for the presence of structural characteristics of promiscuous compounds that were dependent on given target combinations, but not generalizable. Feature weighting and mapping identified characteristic substructures in test compounds. Taken together, these findings are relevant for the design of compounds with desired multi-target activity.

Structured data sets of compounds with multi-target and corresponding single-target activity from biological assays

Aim:

Providing compound data sets for promiscuity analysis with single-target (ST) and multi-target (MT) activity, taking confirmed inactivity against targets into account.

Methodology:

Compounds and target annotations are extracted from screening assays. For a given combination of targets, MT and ST compounds are identified, ensuring test data completeness.

Exemplary results & data:

A total of 1242 MT compounds active against five or more targets and 6629 corresponding ST compounds are characterized, organized and made freely available.

Limitations & next steps:

Screening campaigns typically cover a smaller target space than compounds from the medicinal chemistry literature and their activity annotations might be of lesser quality. Reported compound groups will be subjected to target set-based promiscuity analysis and predictions.

The ability of a compound to bind to multiple biological targets by defined mechanisms is termed promiscuity. Analyzing compound promiscuity helps to better understand how drugs function that are capable of interacting with multiple therapeutic targets. In drug discovery, this phenomenon is referred to as polypharmacology. Machine learning using data sets of compounds with multi-target and corresponding single-target activity aids in identifying structural features that distinguish these compounds.

Drug design targeting active posttranslational modification protein isoforms

Modern drug design aims to discover novel lead compounds with attractable chemical profiles to enable further exploration of the intersection of chemical space and biological space. Identification of small molecules with good ligand efficiency, high activity, and selectivity is crucial toward developing effective and safe drugs. However, the intersection is one of the most challenging tasks in the pharmaceutical industry, as chemical space is almost infinity and continuous, whereas the biological space is very limited and discrete. This bottleneck potentially limits the discovery of molecules with desirable properties for lead optimization. Herein, we present a new direction leveraging posttranslational modification (PTM) protein isoforms target space to inspire drug design termed as "Post-translational Modification Inspired Drug Design (PTMI-DD)." PTMI-DD aims to extend the intersections of chemical space and biological space. We further rationalized and highlighted the importance of PTM protein isoforms and their roles in various diseases and biological functions. We then laid out a few directions to elaborate the PTMI-DD in drug design including discovering covalent binding inhibitors mimicking PTMs, targeting PTM protein isoforms with distinctive binding sites from that of wild-type counterpart, targeting protein-protein interactions involving PTMs, and hijacking protein degeneration by ubiquitination for PTM protein isoforms. These directions will lead to a significant expansion of the biological space and/or increase the tractability of compounds, primarily due to precisely targeting PTM protein isoforms or complexes which are highly relevant to biological functions. Importantly, this new avenue will further enrich the personalized treatment opportunity through precision medicine targeting PTM isoforms.

Analysis of Biological Screening Compounds with Single- or Multi-Target Activity via Diagnostic Machine Learning

Predicting compounds with single- and multi-target activity and exploring origins of compound specificity and promiscuity is of high interest for chemical biology and drug discovery. We present a large-scale analysis of compound promiscuity including two major components. First, high-confidence datasets of compounds with multi- and corresponding single-target activity were extracted from biological screening data. Positive and negative assay results were taken into account and data completeness was ensured. Second, these datasets were investigated using diagnostic machine learning to systematically distinguish between compounds with multi- and single-target activity. Models built on the basis of chemical structure consistently produced meaningful predictions. These findings provided evidence for the presence of structural features differentiating promiscuous and non-promiscuous compounds. Machine learning under varying conditions using modified datasets revealed a strong influence of nearest neighbor relationship on the predictions. Many multi-target compounds were found to be more similar to other multi-target compounds than single-target compounds and vice versa, which resulted in consistently accurate predictions. The results of our study confirm the presence of structural relationships that differentiate promiscuous and non-promiscuous compounds.

Discovery of novel 1,3,5-triazine as adenosine A2A receptor antagonist for benefit in Parkinson's disease

Parkinson's disease (PD) is a chronic neuro-degenerative ailment characterized by impairment in various motor and nonmotor functions of the body. In the past few years, adenosine A₂ A receptor (A₂ AR) antagonists have attracted much attention due to significant relief in PD. Therefore, in the current study, we intend to disclose the development of novel 1,3,5-triazines as A₂ AR antagonist. The radioligand binding and selectivity of analogs were tested in HEK293 (human embryonic kidney) and the cells were transfected with pcDNA 3.1(+) containing full-length human A₂ AR cDNA and pcDNA 3.1(+) containing full-length human A₁ R cDNA, where they exhibit selective affinity for A₂ AR. Molecular docking analysis was also conducted to rationalize the probable mode of action, binding affinity, and orientation of the most potent molecule (7c) at the active site of A₂ AR. It has been shown that compound 7c form numerous nonbonded interactions in the active site of A₂ AR by interacting with Ala59, Ala63, Ile80, Val84 Glu169, Phe168, Met270, and Ile274. The study revealed 1,3,5-triazines as a novel class of A₂ AR antagonists.

Searching for target-specific and multi-targeting organics for Covid-19 in the Drugbank database with a double scoring approach

The current outbreak of Covid-19 infection due to SARS-CoV-2, a virus from the coronavirus family, has become a major threat to human healthcare. The virus has already infected more than 44 M people and the number of deaths reported has reached more than 1.1 M which may be attributed to lack of medicine. The traditional drug discovery approach involves many years of rigorous research and development and demands for a huge investment which cannot be adopted for the ongoing pandemic infection. Rather we need a swift and cost-effective approach to inhibit and control the viral infection. With the help of computational screening approaches and by choosing appropriate chemical space, it is possible to identify lead drug-like compounds for Covid-19. In this study, we have used the Drugbank database to screen compounds against the most important viral targets namely 3C-like protease (3CLpro), papain-like protease (PLpro), RNA-dependent RNA polymerase (RdRp) and the spike (S) protein. These targets play a major role in the replication/transcription and host cell recognition, therefore, are vital for the viral reproduction and spread of infection. As the structure based computational screening approaches are more reliable, we used the crystal structures for 3C-like main protease and spike protein. For the remaining targets, we used the structures based on homology modeling. Further, we employed two scoring methods based on binding free energies implemented in AutoDock Vina and molecular mechanics-generalized Born surface area approach. Based on these results, we propose drug cocktails active against the three viral targets namely 3CLpro, PLpro and RdRp. Interestingly, one of the identified compounds in this study i.e. Baloxavir marboxil has been under clinical trial for the treatment of Covid-19 infection. In addition, we have identified a few compounds such as Phthalocyanine, Tadalafil, Lonafarnib, Nilotinib, Dihydroergotamine, R-428 which can bind to all three targets simultaneously and can serve as multi-targeting drugs. Our study also included calculation of binding energies for various compounds currently under drug trials. Among these compounds, it is found that Remdesivir binds to targets, 3CLpro and RdRp with high binding affinity. Moreover, Baricitinib and Umifenovir were found to have superior target-specific binding while Darunavir is found to be a potential multi-targeting drug. As far as we know this is the first study where the compounds from the Drugbank database are screened against four vital targets of SARS-CoV-2 and illustrates that the computational screening using a double scoring approach can yield potential drug-like compounds against Covid-19 infection.

Systematic Data Analysis and Diagnostic Machine Learning Reveal Differences between Compounds with Single- and Multitarget Activity

Small molecules with multitarget activity are capable of triggering polypharmacological effects and are of high interest in drug discovery. Compared to single-target compounds, promiscuity also affects drug distribution and pharmacodynamics and alters ADMET characteristics. Features distinguishing between compounds with single- and multitarget activity are currently only little understood. On the basis of systematic data analysis, we have assembled large sets of promiscuous compounds with activity against related or functionally distinct targets and the corresponding compounds with single-target activity. Machine learning predicted promiscuous compounds with surprisingly high accuracy. Molecular similarity analysis combined with control calculations under varying conditions revealed that accurate predictions were largely determined by structural nearest-neighbor relationships between compounds from different classes. We also found that large proportions of promiscuous compounds with activity against related or unrelated targets and corresponding single-target compounds formed analog series with distinct chemical space coverage, which further rationalized the predictions. Moreover, compounds with activity against proteins from functionally distinct classes were often active against unique targets that were not covered by other promiscuous compounds. The results of our analysis revealed that nearest-neighbor effects determined the prediction of promiscuous compounds and that preferential partitioning of compounds with single- and multitarget activity into structurally distinct analog series was responsible for such effects, hence providing a rationale for the presence of different structure-promiscuity relationships.

Cannabinoid Receptor Subtype 2 (CB2R) in a Multitarget Approach: Perspective of an Innovative Strategy in Cancer and Neurodegeneration

The cannabinoid receptor subtype 2 (CB2R) represents an interesting and new therapeutic target for its involvement in the first steps of neurodegeneration as well as in cancer onset and progression. Several studies, focused on different types of tumors, report a promising anticancer activity induced by CB2R agonists due to their ability to reduce inflammation and cell proliferation. Moreover, in neuroinflammation, the stimulation of CB2R, overexpressed in microglial cells, exerts beneficial effects in neurodegenerative disorders. With the aim to overcome current treatment limitations, new drugs can be developed by specifically modulating, together with CB2R, other targets involved in such multifactorial disorders. Building on successful case studies of already developed multitarget strategies involving CB2R, in this Perspective we aim at prompting the scientific community to consider new promising target associations involving HDACs (histone deacetylases) and σ receptors by employing modern approaches based on molecular hybridization, computational polypharmacology, and machine learning algorithms.

Molecular dynamics simulation, free energy landscape and binding free energy computations in exploration the anti-invasive activity of amygdalin against metastasis

BACKGROUND AND OBJECTIVE

Historically, amygdalin has been used as alternative medicine or in vitro and in vivo studies, but no single study exists which discusses the structural mechanism of amygdalin at a molecular level. This paper inquiries into the inhibitory actions of amygdalin on the selected targets: AKT1, FAK, and ILK, which are regulators for various mediated signaling pathways, and are associated with cell adhesion, migration, and differentiation. In order to get details at the molecular level of amygdalin's inhibitory activities against chosen proteins, molecular modeling and simulation techniques including double docking, molecular dynamics simulation, free energy landscape analysis, and binding free energy calculation were exerted.

METHODS

To get molecular level details of amygdalin inhibitory effects against the relevant proteins; here the utilized tools are the following: the double docking, molecular dynamics simulation, free energy landscape analysis, g_mmpbsa, and interaction entropy were used to evaluate the inhibitory activity against targeted proteins.

RESULTS

The computational calculations revealed that amygdalin inhibits the selected targets via block the ATP-binding pocket of AKT1, FAK, and ILK by forming stable hydrogen bonds. Moreover, free energy landscape, FEL exposed that amygdalin stabilized the global conformations of both FAK and ILK proteins to the minimum global energy besides it reduced the essential dynamics of FAK and ILK proteins. MMPBSA computations provided further evidence for amygdalin's stability inside the ATP-binding pocket of AKT1, FAK, and ILK with a binding free energy of 45.067, -13.033, 13.109 kJ/mol, respectively. The binding free energies are lastly consistent with the hydrogen bonding and pairs within 0.35 nm results. The decomposition of binding energy shows the pivotal amino acid residues responsible for the stability of amygdalin's interactions inside the ATP-binding sites by forming hydrogen bonds.

CONCLUSIONS

Before this work, it was enigmatic to make predictions about how amygdalin inhibits metastasis of cancer. But the computational results contribute in several ways to our understanding of amygdalin activity and provide a basic insight into the activity of amygdalin as a multi-target drug in the metastasis and invasion of cancer.

Recent progress on cheminformatics approaches to epigenetic drug discovery

The ability of epigenetic markers to affect genome function has enabled transformative changes in drug discovery, especially in cancer and other emerging therapeutic areas. Concordant with the introduction of the term 'epi-informatics', the size of the epigenetically relevant chemical space has grown substantially and so did the number of applications of cheminformatic methods to epigenetics. Recent progress in epi-informatics has improved our understanding of the structure-epigenetic activity relationships and boosted the development of models predicting novel epigenetic agents. Herein, we review the advances in computational approaches to drug discovery of small molecules with epigenetic modulation profiles, summarize the current chemogenomics data available for epigenetic targets, and provide a perspective on the greater utility of biomedical knowledge mining as a means to advance the epigenetic drug discovery.

Computational/in silico methods in drug target and lead prediction

Drug-like compounds are most of the time denied approval and use owing to the unexpected clinical side effects and cross-reactivity observed during clinical trials. These unexpected outcomes resulting in significant increase in attrition rate centralizes on the selected drug targets. These targets may be disease candidate proteins or genes, biological pathways, disease-associated microRNAs, disease-related biomarkers, abnormal molecular phenotypes, crucial nodes of biological network or molecular functions. This is generally linked to several factors, including incomplete knowledge on the drug targets and unpredicted pharmacokinetic expressions upon target interaction or off-target effects. A method used to identify targets, especially for polygenic diseases, is essential and constitutes a major bottleneck in drug development with the fundamental stage being the identification and validation of drug targets of interest for further downstream processes. Thus, various computational methods have been developed to complement experimental approaches in drug discovery. Here, we present an overview of various computational methods and tools applied in predicting or validating drug targets and drug-like molecules. We provide an overview on their advantages and compare these methods to identify effective methods which likely lead to optimal results. We also explore major sources of drug failure considering the challenges and opportunities involved. This review might guide researchers on selecting the most efficient approach or technique during the computational drug discovery process.

Multi-targeted drug design strategies for the treatment of schizophrenia

INTRODUCTION

Schizophrenia is a complex psychiatric disease (or a conglomeration of disorders) manifesting with positive, negative and cognitive symptoms. The pathophysiology of schizophrenia is not completely known; however, it involves many neurotransmitters and their receptors. In order to treat schizophrenia, drugs need to be multi-target drugs. Indeed, the action of second and third generation antipsychotics involves interactions with many receptors, belonging mainly to aminergic GPCRs.

AREAS COVERED

In this review, the authors summarize current concepts of schizophrenia with the emphasis on the modern dopaminergic, serotoninergic, and glutamatergic hypotheses. Next, they discuss treatments of the disease, stressing multi-target antipsychotics. They cover different aspects of design of multi-target ligands, including the application of molecular modeling approaches for the design and benefits and limitations of multifunctional compounds. Finally, they present successful case studies of multi-target drug design against schizophrenia.

EXPERT OPINION

Treatment of schizophrenia requires the application of multi-target drugs. While designing single target drugs is relatively easy, designing multifunctional compounds is a challenge due to the necessity to balance the affinity to many targets, while avoiding promiscuity and the problems with drug-likeness. Multi-target drugs bring many benefits: better efficiency, fewer adverse effects, and drug-drug interactions and better patient compliance to drug regime.

Computational methods-guided design of modulators targeting protein-protein interactions (PPIs)

Protein-protein interactions (PPIs) play a pivotal role in extensive biological processes and are thus crucial to human health and the development of disease states. Due to their critical implications, PPIs have been spotlighted as promising drug targets of broad-spectrum therapeutic interests. However, owing to the general properties of PPIs, such as flat surfaces, featureless conformations, difficult topologies, and shallow pockets, previous attempts were faced with serious obstacles when targeting PPIs and almost portrayed them as "intractable" for decades. To date, rapid progress in computational chemistry and structural biology methods has promoted the exploitation of PPIs in drug discovery. These techniques boost their cost-effective and high-throughput traits, and enable the study of dynamic PPI interfaces. Thus, computational methods represent an alternative strategy to target "undruggable" PPI interfaces and have attracted intense pharmaceutical interest in recent years, as exemplified by the accumulating number of successful cases. In this review, we first introduce a diverse set of computational methods used to design PPI modulators. Herein, we focus on the recent progress in computational strategies and provide a comprehensive overview covering various methodologies. Importantly, a list of recently-reported successful examples is highlighted to verify the feasibility of these computational approaches. Finally, we conclude the general role of computational methods in targeting PPIs, and also discuss future perspectives on the development of such aids.

Zinc-mediated conformational preselection mechanism in the allosteric control of DNA binding to the zinc transcriptional regulator (ZitR)

The zinc transcriptional regulator (ZitR) functions as a metalloregulator that fine tunes transcriptional regulation through zinc-dependent DNA binding. However, the molecular mechanism of zinc-driven allosteric control of the DNA binding to ZitR remains elusive. Here, we performed enhanced sampling accelerated molecular dynamics simulations to figure out the mechanism, revealing the role of protein dynamics in the zinc-induced allosteric control of DNA binding to ZitR. The results suggest that zinc-free ZitR samples distinct conformational states, only a handful of which are compatible with DNA binding. Remarkably, zinc binding reduces the conformational plasticity of the DNA-binding domain of ZitR, promoting the population shift in the ZitR conformational ensemble towards the DNA binding-competent conformation. Further co-binding of DNA to the zinc–ZitR complex stabilizes this competent conformation. These findings suggest that ZitR–DNA interactions are allosterically regulated in a zinc-mediated conformational preselection manner, highlighting the importance of conformational dynamics in the regulation of transcription factor family.

Nanonetworks in Biomedical Applications

By interconnecting nanomachines and forming nanonetworks, the capacities of single nanomachines are expected to be enhanced, as the ensuing information exchange will allow them to cooperate towards a common goal. Nowadays, systems normally use electromagnetic signals to encode, send and receive information, however, in a novel communication paradigm, molecular transceivers, channel models or protocols use molecules. This article presents the current developments in nanomachines along with their future architecture to better understand nanonetwork scenarios in biomedical applications. Furthermore, to highlight the communication needs between nanomachines, two applications for nanonetworks are also presented: i) a new networking paradigm, called the Internet of NanoThings, that allows nanoscale devices to interconnect with existing communication networks, and ii) Molecular Communication, where the propagation of chemical compounds like drug particles, carry out the information exchange.

Predicting drug properties with parameter-free machine learning: pareto-optimal embedded modeling (POEM)

The prediction of absorption, distribution, metabolism, excretion, and toxicity (ADMET) of small molecules from their molecular structure is a central problem in medicinal chemistry with great practical importance in drug discovery. Creating predictive models conventionally requires substantial trial-and-error for the selection of molecular representations, machine learning (ML) algorithms, and hyperparameter tuning. A generally applicable method that performs well on all datasets without tuning would be of great value but is currently lacking. Here, we describe pareto-optimal embedded modeling (POEM), a similarity-based method for predicting molecular properties. POEM is a non-parametric, supervised ML algorithm developed to generate reliable predictive models without need for optimization. POEM’s predictive strength is obtained by combining multiple different representations of molecular structures in a context-specific manner, while maintaining low dimensionality. We benchmark POEM relative to industry-standard ML algorithms and published results across 17 classifications tasks. POEM performs well in all cases and reduces the risk of overfitting.

Epigenetic Modulators as Potential Multi-targeted Drugs Against Hedgehog Pathway for Treatment of Cancer

The Sonic hedgehog signalling is known to play a crucial role in regulating embryonic development, cancer stem cell maintenance and tissue patterning. Dysregulated hedgehog signalling has been reported to affect tumorigenesis and drug response in various human malignancies. Epigenetic therapy relying on DNA methyltransferase and Histone deacetylase inhibitors are being proposed as potential drug candidates considering their efficiency in preventing development of cancer progenitor cells, killing drug resistant cells and also dictating "on/off" switch of tumor suppressor genes and oncogenes. In this docking approach, epigenetic modulators were virtually screened for their efficiency in inhibiting key regulators of SHH pathway viz., sonic hedgehog, Smoothened and Gli using polypharmacological approach. The control drugs and epigenetic modulators were docked with PDB protein structures using AutoDock vina and further checked for their drug-likeness properties. Further molecular dynamics simulation using VMD and NAMD, and MMP/GBSA energy calculation were employed for verifying the stability and entropy of the ligand-receptor complex. EPZ-6438 and GSK 343 (EZH2 inhibitors), CHR 3996 and Mocetinostat (HDAC inhibitors), GSK 126 (HKMT inhibitor) and UNC 1215 (L3MBTL3 antagonist) exhibited multiple-targeted approach in modulating HH signalling. This is the first study to report these epigenetic drugs as potential multi-targeted hedgehog pathway inhibitors. Thus, epigenetic polypharmacology approach can be explored as a better alternative to challenges of acute long term toxicity and drug resistance occurring due to traditional single targeted chemotherapy in the future.

Prospects of multitarget drug designing strategies by linking molecular docking and molecular dynamics to explore the protein-ligand recognition process

The designing of drugs that can simultaneously affect different protein targets is one novel and promising way to treat complex diseases. Multitarget drugs act on multiple protein receptors each implicated in the same disease state, and may be considered to be more beneficial than conventional drug therapies. For example, these drugs can have improved therapeutic potency due to synergistic effects on multiple targets, as well as improved safety and resistance profiles due to the combined regulation of potential primary therapeutic targets and compensatory elements and lower dosage typically required. This review analyzes in-silico methods that facilitate multitarget drug design that facilitate the discovery and development of novel therapeutic agents. Here presented is a summary of the progress in structure-based drug discovery techniques that study the process of molecular recognition of targets and ligands, moving from static molecular docking to improved molecular dynamics approaches in multitarget drug design, and the advantages and limitations of each.

A Definition of "Multitargeticity": Identifying Potential Multitarget and Selective Ligands Through a Vector Analysis

The design of multitarget drugs is an essential area of research in Medicinal Chemistry since they have been proposed as potential therapeutics for the management of complex diseases. However, defining a multitarget drug is not an easy task. In this work, we propose a vector analysis for measuring and defining "multitargeticity." We developed terms, such as order and force of a ligand, to finally reach two parameters: multitarget indexes 1 and 2. The combination of these two indexes allows discrimination of multitarget drugs. Several training sets were constructed to test the usefulness of the indexes: an experimental training set, with real affinities, a docking training set, within theoretical values, and an extensive database training set. The indexes proved to be useful, as they were used independently in silico and experimental data, identifying actual multitarget compounds and even selective ligands in most of the training sets. We then applied these indexes to evaluate a virtual library of potential ligands for targets related to multiple sclerosis, identifying 10 compounds that are likely leads for the development of multitarget drugs based on their in silico behavior. With this work, a new milestone is made in the way of defining multitargeticity and in drug design.

LigBuilder V3: A Multi-Target de novo Drug Design Approach

With the rapid development of systems-based pharmacology and poly-pharmacology, method development for rational design of multi-target drugs has becoming urgent. In this paper, we present the first de novo multi-target drug design program LigBuilder V3, which can be used to design ligands to target multiple receptors, multiple binding sites of one receptor, or various conformations of one receptor. LigBuilder V3 is generally applicable in de novo multi-target drug design and optimization, especially for the design of concise ligands for protein targets with large difference in binding sites. To demonstrate the utility of LigBuilder V3, we have used it to design dual-functional inhibitors targeting HIV protease and HIV reverse transcriptase with three different strategy, including multi-target de novo design, multi-target growing, and multi-target linking. The designed compounds were computational validated by MM/GBSA binding free energy estimation as highly potential multi-target inhibitors for both HIV protease and HIV reverse transcriptase. The LigBuilder V3 program can be downloaded at “http://www.pkumdl.cn/ligbuilder3/”.

Discovery of Multitarget-Directed Ligands Against Influenza A Virus From Compound Yizhihao Through a Predictive System for Compound-Protein Interactions

Influenza A virus (IAV) is a threat to public health due to its high mutation rate and resistance to existing drugs. In this investigation, 15 targets selected from an influenza virus–host interaction network were successfully constructed as a multitarget virtual screening system for new drug discovery against IAV using Naïve Bayesian, recursive partitioning, and CDOCKER methods. The predictive accuracies of the models were evaluated using training sets and test sets. The system was then used to predict active constituents of Compound Yizhihao (CYZH), a Chinese medicinal compound used to treat influenza. Twenty-eight compounds with multitarget activities were selected for subsequent in vitro evaluation. Of the four compounds predicted to be active on neuraminidase (NA), chlorogenic acid, and orientin showed inhibitory activity in vitro. Linarin, sinensetin, cedar acid, isoliquiritigenin, sinigrin, luteolin, chlorogenic acid, orientin, epigoitrin, and rupestonic acid exhibited significant effects on TNF-α expression, which is almost consistent with predicted results. Results from a cytopathic effect (CPE) reduction assay revealed acacetin, indirubin, tryptanthrin, quercetin, luteolin, emodin, and apigenin had protective effects against wild-type strains of IAV. Quercetin, luteolin, and apigenin had good efficacy against resistant IAV strains in CPE reduction assays. Finally, with the aid of Gene Ontology biological process analysis, the potential mechanisms of CYZH action were revealed. In conclusion, a compound-protein interaction-prediction system was an efficient tool for the discovery of novel compounds against influenza, and the findings from CYZH provide important information for its usage and development.

Fingerprinting CANDO: Increased Accuracy with Structure- and Ligand-Based Shotgun Drug Repurposing

We have upgraded our Computational Analysis of Novel Drug Opportunities (CANDO) platform for shotgun drug repurposing by including ligand-based, data fusion, and decision tree pipelines. The goal of shotgun drug repurposing is to screen and rank every existing human use drug or compound for every disease/indication. The first version of CANDO implemented a structure-based pipeline that modeled interactions between compounds and proteins on a large scale, generating compound-proteome interaction signatures used to infer the similarity of drug behavior; the new pipelines accomplish this by incorporating molecular fingerprints and the Tanimoto coefficient. We obtain improved benchmarking performance with the new pipelines across all three evaluation metrics used: average indication accuracy, pairwise accuracy, and coverage. The best performing pipeline achieves an average indication accuracy of 19.0% at the top10 cutoff, compared to 11.7% for v1, and 2.2% for a random control. Our results demonstrate that the CANDO drug recovery accuracy is substantially improved by integrating multiple pipelines, thereby enhancing our ability to generate putative therapeutic repurposing candidates, and increasing drug discovery efficiency.