Structural and Functional View of Polypharmacology

Predicting human and viral protein variants affecting COVID-19 susceptibility and repurposing therapeutics

The COVID-19 disease is an ongoing global health concern. Although vaccination provides some protection, people are still susceptible to re-infection. Ostensibly, certain populations or clinical groups may be more vulnerable. Factors causing these differences are unclear and whilst socioeconomic and cultural differences are likely to be important, human genetic factors could influence susceptibility. Experimental studies indicate SARS-CoV-2 uses innate immune suppression as a strategy to speed-up entry and replication into the host cell. Therefore, it is necessary to understand the impact of variants in immunity-associated human proteins on susceptibility to COVID-19. In this work, we analysed missense coding variants in several SARS-CoV-2 proteins and their human protein interactors that could enhance binding affinity to SARS-CoV-2. We curated a dataset of 19 SARS-CoV-2: human protein 3D-complexes, from the experimentally determined structures in the Protein Data Bank and models built using AlphaFold2-multimer, and analysed the impact of missense variants occurring in the protein-protein interface region. We analysed 468 missense variants from human proteins and 212 variants from SARS-CoV-2 proteins and computationally predicted their impacts on binding affinities for the human viral protein complexes. We predicted a total of 26 affinity-enhancing variants from 13 human proteins implicated in increased binding affinity to SARS-CoV-2. These include key-immunity associated genes (TOMM70, ISG15, IFIH1, IFIT2, RPS3, PALS1, NUP98, AXL, ARF6, TRIMM, TRIM25) as well as important spike receptors (KREMEN1, AXL and ACE2). We report both common (e.g., Y13N in IFIH1) and rare variants in these proteins and discuss their likely structural and functional impact, using information on known and predicted functional sites. Potential mechanisms associated with immune suppression implicated by these variants are discussed. Occurrence of certain predicted affinity-enhancing variants should be monitored as they could lead to increased susceptibility and reduced immune response to SARS-CoV-2 infection in individuals/populations carrying them. Our analyses aid in understanding the potential impact of genetic variation in immunity-associated proteins on COVID-19 susceptibility and help guide drug-repurposing strategies.

Systematic identification and repurposing of FDA-approved drugs as antibacterial agents against Streptococcus pyogenes: In silico and in vitro studies

Streptococcus pyogenes (Group A Streptococcus - GAS) is a human pathogen causing wide range of infections and toxin-mediated diseases in human beings of all age groups with fatality of 10-30 %. The limited success of antibiotics and the non-availability of vaccines makes GAS a global burden. The multi-subunit RNA polymerase (RNAP) is a validated bacterial therapeutic target as it is involved in transcription and can arrest growth. Of the five subunits of this enzyme complex, the β-subunit (RpoC) has attracted specific attention as a drug target, particularly in the switch region. Here we attempt to repurpose non-antimicrobial drugs to act as RpoC inhibitors against S. pyogenes. In this study, 1826 FDA approved drugs have been identified through high-throughput virtual screening. Free Energy Perturbation (FEP) based binding free energy calculations have been performed at the final step of the virtual screening funnel to ensure high accuracy in silico results. Three compounds identified have been tested for susceptibility of S. pyogenes MTCC 442 strain and two antibiotic-resistant clinical isolates of S. pyogenes using microdilution assay. Among the three, two drugs Amlodipine Besylate (Amd) and Ranitidine hydrochloride (Rnt) have shown inhibition against all the tested strains and its mechanism of interaction with RpoC has been studied. The docked complexes were analyzed to understand the binding mode of the drugs to the target. Classical Molecular Dynamics studies for RpoC-Rnt complex and the two stable conformations of RpoC-Amd complex was carried out. Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of Gyration (RoG) and Solvent Accessible Surface Area (SASA) of the complexes were plotted and studied. The thermodynamic parameters of protein-drug were experimentally determined using Isothermal Titration Calorimetry (ITC). Infrared spectroscopic studies and Fluorescence quenching studies provided insights into the secondary structural changes in RpoC on binding to the drugs.

Integrated untargeted and targeted proteomics to unveil plasma prognostic markers for patients with acute paraquat poisoning: A pilot study

Paraquat (PQ) is a widely used but strongly toxic herbicide, which can induce multiple organ failure. The overall survival rate of the poisoned patients was only 54.4% due to lack of specific antidotes. Besides, the definite pathogenic mechanism of PQ is still not fully understood. In this pilot study, untargeted and targeted proteomics were integrated to explore the expression characteristics of plasma protein in PQ poisoned patients, and identify the differentially expressed proteins between survivors and non-survivors. A total of 494 plasma proteins were detected, and of which 47 were upregulated and 44 were downregulated in PQ poisoned patients compared to healthy controls. Among them, five differential plasma proteins (S100A9, S100A8, MB, ACTB and RAB11FIP3) were further validated by multiple reaction monitoring (MRM)-based targeted proteomic approach, and three of them (S100A9, S100A8 and ACTB) were confirmed to be correlated with PQ poisoning. Meanwhile, 84 dysregulated plasma proteins were identified in non-survivors compared with survivors. Moreover, targeted proteomic and ROC analysis suggested that ACTB had a good performance in predicting the prognosis of PQ poisoned patients. These findings highlighted the value of label-free and mass spectrometry-based proteomics in screening prognostic biomarkers of PQ poisoning and studying the mechanism of PQ toxicity.

Network analyses reveal new insights into the effect of multicomponent Tr14 compared to single-component diclofenac in an acute inflammation model

BACKGROUND

Modifying the acute inflammatory response has wide clinical benefits. Current options include non-steroidal anti-inflammatory drugs (NSAIDs) and therapies that may resolve inflammation. Acute inflammation involves multiple cell types and various processes. We, therefore, investigated whether an immunomodulatory drug that acts simultaneously at multiple sites shows greater potential to resolve acute inflammation more effectively and with fewer side effects than a common anti-inflammatory drug developed as a small molecule for a single target. In this work, we used time-series gene expression profiles from a wound healing mouse model to compare the effects of Traumeel (Tr14), a multicomponent natural product, to diclofenac, a single component NSAID on inflammation resolution.

RESULTS

We advance previous studies by mapping the data onto the "Atlas of Inflammation Resolution", followed by in silico simulations and network analysis. We found that Tr14 acts primarily on the late phase of acute inflammation (during resolution) compared to diclofenac, which suppresses acute inflammation immediately after injury.

CONCLUSIONS

Our results provide new insights how network pharmacology of multicomponent drugs may support inflammation resolution in inflammatory conditions.

KinFams: De-Novo Classification of Protein Kinases Using CATH Functional Units

Protein kinases are important targets for treating human disorders, and they are the second most targeted families after G-protein coupled receptors. Several resources provide classification of kinases into evolutionary families (based on sequence homology); however, very few systematically classify functional families (FunFams) comprising evolutionary relatives that share similar functional properties. We have developed the FunFam-MARC (Multidomain ARchitecture-based Clustering) protocol, which uses multi-domain architectures of protein kinases and specificity-determining residues for functional family classification. FunFam-MARC predicts 2210 kinase functional families (KinFams), which have increased functional coherence, in terms of EC annotations, compared to the widely used KinBase classification. Our protocol provides a comprehensive classification for kinase sequences from >10,000 organisms. We associate human KinFams with diseases and drugs and identify 28 druggable human KinFams, i.e., enriched in clinically approved drugs. Since relatives in the same druggable KinFam tend to be structurally conserved, including the drug-binding site, these KinFams may be valuable for shortlisting therapeutic targets. Information on the human KinFams and associated 3D structures from AlphaFold2 are provided via our CATH FTP website and Zenodo. This gives the domain structure representative of each KinFam together with information on any drug compounds available. For 32% of the KinFams, we provide information on highly conserved residue sites that may be associated with specificity.

Integrating human brain proteomic data with genome-wide association study findings identifies novel brain proteins in substance use traits

Despite the identification of a growing number of genetic risk loci for substance use traits (SUTs), the impact of these loci on protein abundance and the potential utility of relevant proteins as therapeutic targets are unknown. We conducted a proteome-wide association study (PWAS) in which we integrated human brain proteomes from discovery (Banner; N = 152) and validation (ROSMAP; N = 376) datasets with genome-wide association study (GWAS) summary statistics for 4 SUTs. The 4 samples comprised GWAS of European-ancestry individuals for smoking initiation [Smk] (N = 1,232,091), alcohol use disorder [AUD] (N = 313,959), cannabis use disorder [CUD] (N = 384,032), and opioid use disorder [OUD] (N = 302,585). We conducted transcriptome-wide association studies (TWAS) with human brain transcriptomic data to examine the overlap of genetic effects at the proteomic and transcriptomic levels and characterize significant genes through conditional, colocalization, and fine-mapping analyses. We identified 27 genes (Smk = 21, AUD = 3, CUD = 2, OUD = 1) that were significantly associated with cis-regulated brain protein abundance. Of these, 7 showed evidence for causality (Smk: NT5C2, GMPPB, NQO1, RHOT2, SRR and ACTR1B; and AUD: CTNND1). Cis-regulated transcript levels for 8 genes (Smk = 6, CUD = 1, OUD = 1) were associated with SUTs, indicating that genetic loci could confer risk for these SUTs by modulating both gene expression and proteomic abundance. Functional studies of the high-confidence risk proteins identified here are needed to determine whether they are modifiable targets and useful in developing medications and biomarkers for these SUTs.

Identification of New Toxicity Mechanisms in Drug-Induced Liver Injury through Systems Pharmacology

Among adverse drug reactions, drug-induced liver injury presents particular challenges because of its complexity, and the underlying mechanisms are still not completely characterized. Our knowledge of the topic is limited and based on the assumption that a drug acts on one molecular target. We have leveraged drug polypharmacology, i.e., the ability of a drug to bind multiple targets and thus perturb several biological processes, to develop a systems pharmacology platform that integrates all drug–target interactions. Our analysis sheds light on the molecular mechanisms of drugs involved in drug-induced liver injury and provides new hypotheses to study this phenomenon.

Porphyromonas gingivalis resistance and virulence: An integrated functional network analysis

BACKGROUND

The gram-negative bacteria Porphyromonas gingivalis (PG) is the most prevalent cause of periodontal diseases and multidrug-resistant (MDR) infections. Periodontitis and MDR infections are severe due to PG's ability to efflux antimicrobial and virulence factors. This gives rise to colonisation, biofilm development, evasion, and modulation of the host defence system. Despite extensive studies on the MDR efflux pump in other pathogens, little is known about the efflux pump and its association with the virulence factor in PG. Prolonged infection of PG leads to complete loss of teeth and other systemic diseases. This necessitates the development of new therapeutic interventions to prevent and control MDR.

OBJECTIVE

The study aims to identify the most indispensable proteins that regulate both resistance and virulence in PG, which could therefore be used as a target to fight against the MDR threat to antibiotics.

METHODS

We have adopted a hierarchical network-based approach to construct a protein interaction network. Firstly, individual networks of four major efflux pump proteins and two virulence regulatory proteins were constructed, followed by integrating them into one. The relationship between proteins was investigated using a combination of centrality scores, k-core network decomposition, and functional annotation, to computationally identify the indispensable proteins.

RESULTS

Our study identified four topologically significant genes, PG_0538, PG_0539, PG_0285, and PG_1797, as potential pharmacological targets. PG_0539 and PG_1797 were identified to have significant associations between the efflux pump and virulence genes. This type of underpinning research may help in narrowing the drug spectrum used for treating periodontal diseases, and may also be exploited to look into antibiotic resistance and pathogenicity in bacteria other than PG.

Mass Spectrometric Behavior and Molecular Mechanisms of Fermented Deoxyanthocyanidins to Alleviate Ulcerative Colitis Based on Network Pharmacology

Aims. Ulcerative colitis (UC) is a type of chronic idiopathic inflammatory bowel disease with a multifactorial pathogenesis and limited treatment options. The aim of the present study is to investigate the hydrogen deuterium exchange mass spectrometry (HDX-MS) behaviors of fermented deoxyanthocyanidins and their molecular mechanisms to alleviate UC by using quantum chemistry and network pharmacology. Methods. Tandem MS indicated at least two fragmentation pathways through which deuterated vinylphenol-deoxyanthocyanidins could generate different product ions. Quantum calculations were conducted to determine the transition states of the relevant molecules and analyze their optimized configuration, vibrational characteristics, intrinsic reaction coordinates, and corresponding energies. The potential targets of deoxyanthocyanidins in UC were screened from a public database. The R package was used for Gene Ontology (GO) and KEGG pathway analyses, and the protein–protein interactions (PPIs) of the targets were assessed using Search Tool for the Retrieval of Interacting Genes (STRING). Finally, molecular docking was implemented to analyze the binding energies and action modes of the target compounds through the online tool CB-Dock. Results. Quantum calculations indicated two potential fragmentation pathways involving the six-membered ring and dihydrogen cooperative transfer reactions of the vinylphenol-deoxyanthocyanidins. A total of 146 and 57 intersecting targets of natural and fermented deoxyanthocyanidins were separately screened out from the UC database and significant overlaps in GO terms and KEGG pathways were noted. Three shared hub targets (i.e., PTGS2, ESR1, and EGFR) were selected from the two PPI networks by STRING. Molecular docking results showed that all deoxyanthocyanidins have a good binding potential with the hub target proteins and that fermented deoxyanthocyanidins have lower binding energies and more stable conformations compared with natural ones. Conclusions. Deoxyanthocyanidins may provide anti-inflammatory, antioxidative, and immune system regulatory effects to suppress UC progression. It is proposed for the first time that fermentation of deoxyanthocyanidins can help adjust the structure of the intestinal microbiota and increase the biological activity of the natural compounds against UC. Furthermore, HDX-MS is a helpful strategy to analyze deoxyanthocyanidin metabolites with unknown structures.

Exploiting protein family and protein network data to identify novel drug targets for bladder cancer

Bladder cancer remains one of the most common forms of cancer and yet there are limited small molecule targeted therapies. Here, we present a computational platform to identify new potential targets for bladder cancer therapy. Our method initially exploited a set of known driver genes for bladder cancer combined with predicted bladder cancer genes from mutationally enriched protein domain families. We enriched this initial set of genes using protein network data to identify a comprehensive set of 323 putative bladder cancer targets. Pathway and cancer hallmarks analyses highlighted putative mechanisms in agreement with those previously reported for this cancer and revealed protein network modules highly enriched in potential drivers likely to be good targets for targeted therapies. 21 of our potential drug targets are targeted by FDA approved drugs for other diseases — some of them are known drivers or are already being targeted for bladder cancer (FGFR3, ERBB3, HDAC3, EGFR). A further 4 potential drug targets were identified by inheriting drug mappings across our in-house CATH domain functional families (FunFams). Our FunFam data also allowed us to identify drug targets in families that are less prone to side effects i.e., where structurally similar protein domain relatives are less dispersed across the human protein network. We provide information on our novel potential cancer driver genes, together with information on pathways, network modules and hallmarks associated with the predicted and known bladder cancer drivers and we highlight those drivers we predict to be likely drug targets.

Polypharmacology: The science of multi-targeting molecules

Polypharmacology is a concept where a molecule can interact with two or more targets simultaneously. It offers many advantages as compared to the conventional single-targeting molecules. A multi-targeting drug is much more efficacious due to its cumulative efficacy at all of its individual targets making it much more effective in complex and multifactorial diseases like cancer, where multiple proteins and pathways are involved in the onset and development of the disease. For a molecule to be polypharmacologic in nature, it needs to possess promiscuity which is the ability to interact with multiple targets; and at the same time avoid binding to antitargets which would otherwise result in off-target adverse effects. There are certain structural features and physicochemical properties which when present would help researchers to predict if the designed molecule would possess promiscuity or not. Promiscuity can also be identified via advanced state-of-the-art computational methods. In this review, we also elaborate on the methods by which one can intentionally incorporate promiscuity in their molecules and make them polypharmacologic. The polypharmacology paradigm of "one drug-multiple targets" has numerous applications especially in drug repurposing where an already established drug is redeveloped for a new indication. Though designing a polypharmacological drug is much more difficult than designing a single-targeting drug, with the current technologies and information regarding different diseases and chemical functional groups, it is plausible for researchers to intentionally design a polypharmacological drug and unlock its advantages.

Centrality of drug targets in protein networks

Background

In the pharmaceutical industry, competing for few validated drug targets there is a drive to identify new ways of therapeutic intervention. Here, we attempted to define guidelines to evaluate a target’s ‘fitness’ based on its node characteristics within annotated protein functional networks to complement contingent therapeutic hypotheses.

Results

We observed that targets of approved, selective small molecule drugs exhibit high node centrality within protein networks relative to a broader set of investigational targets spanning various development stages. Targets of approved drugs also exhibit higher centrality than other proteins within their respective functional class. These findings expand on previous reports of drug targets’ network centrality by suggesting some centrality metrics such as low topological coefficient as inherent characteristics of a ‘good’ target, relative to other exploratory targets and regardless of its functional class. These centrality metrics could thus be indicators of an individual protein’s ‘fitness’ as potential drug target. Correlations between protein nodes’ network centrality and number of associated publications underscored the possibility of knowledge bias as an inherent limitation to such predictions.

Conclusions

Despite some entanglement with knowledge bias, like structure-oriented ‘druggability’ assessments of new protein targets, centrality metrics could assist early pharmaceutical discovery teams in evaluating potential targets with limited experimental proof of concept and help allocate resources for an effective drug discovery pipeline.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04342-x.

Zebrafish disease models in drug discovery: from preclinical modelling to clinical trials

Numerous drug treatments that have recently entered the clinic or clinical trials have their genesis in zebrafish. Zebrafish are well established for their contribution to developmental biology and have now emerged as a powerful preclinical model for human disease, as their disease characteristics, aetiology and progression, and molecular mechanisms are clinically relevant and highly conserved. Zebrafish respond to small molecules and drug treatments at physiologically relevant dose ranges and, when combined with cell-specific or tissue-specific reporters and gene editing technologies, drug activity can be studied at single-cell resolution within the complexity of a whole animal, across tissues and over an extended timescale. These features enable high-throughput and high-content phenotypic drug screening, repurposing of available drugs for personalized and compassionate use, and even the development of new drug classes. Often, drugs and drug leads explored in zebrafish have an inter-organ mechanism of action and would otherwise not be identified through targeted screening approaches. Here, we discuss how zebrafish is an important model for drug discovery, the process of how these discoveries emerge and future opportunities for maximizing zebrafish potential in medical discoveries.

Brain proteome-wide association study implicates novel proteins in depression pathogenesis

Depression is a common condition, but current treatments are only effective in a subset of individuals. To identify new treatment targets, we integrated depression genome-wide association study (GWAS) results (N = 500,199) with human brain proteomes (N = 376) to perform a proteome-wide association study of depression followed by Mendelian randomization. We identified 19 genes that were consistent with being causal in depression, acting via their respective cis-regulated brain protein abundance. We replicated nine of these genes using an independent depression GWAS (N = 307,353) and another human brain proteomic dataset (N = 152). Eleven of the 19 genes also had cis-regulated mRNA levels that were associated with depression, based on integration of the depression GWAS with human brain transcriptomes (N = 888). Meta-analysis of the discovery and replication proteome-wide association study analyses identified 25 brain proteins consistent with being causal in depression, 20 of which were not previously implicated in depression by GWAS. Together, these findings provide promising brain protein targets for further mechanistic and therapeutic studies.

Design, synthesis and evaluation of novel indole-2-carboxamides for growth inhibition of Mycobacterium tuberculosis and paediatric brain tumour cells

The omnipresent threat of tuberculosis (TB) and the scant treatment options thereof necessitate the development of new antitubercular agents, preferably working via a novel mechanism of action distinct from the current drugs. Various studies identified the mycobacterial membrane protein large 3 transporter (MmpL3) as the target of several classes of compounds, including the indole-2-caboxamides. Herein, several indoleamide analogues were rationally designed, synthesised, and evaluated for their antitubercular and antitumour activities. Compound 8g displayed the highest activity (MIC = 0.32 μM) against the drug-sensitive (DS) Mycobacterium tuberculosis (M. tb) H37Rv strain. This compound also exhibited high selective activity towards M. tb over mammalian cells [IC₅₀ (Vero cells) = 40.9 μM, SI = 128], suggesting its minimal cytotoxicity. In addition, when docked into the MmpL3 active site, 8g adopted a binding profile similar to the indoleamide ligand ICA38. A related compound 8f showed dual antitubercular (MIC = 0.62 μM) and cytotoxic activities against paediatric glioblastoma multiforme (GBM) cell line KNS42 [IC₅₀ (viability) = 0.84 μM]. Compound 8f also showed poor cytotoxic activity against healthy Vero cells (IC₅₀ = 39.9 μM). Compounds 9a and 15, which were inactive against M. tb, showed potent cytotoxic (IC₅₀ = 8.25 and 5.04 μM, respectively) and antiproliferative activities (IC₅₀ = 9.85 and 6.62 μM, respectively) against KNS42 cells. Transcriptional analysis of KNS42 cells treated with compound 15 revealed a significant downregulation in the expression of the carbonic anhydrase 9 (CA9) and the spleen tyrosine kinase (SYK) genes. The expression levels of these genes in GBM tumours were previously shown to contribute to tumour progression, suggesting their involvement in our observed antitumour activities. Compounds 9a and 15 were selected for further evaluations against three different paediatric brain tumour cell lines (BT12, BT16 and DAOY) and non-neoplastic human fibroblast cells HFF1. Compound 9a showed remarkable cytotoxic (IC₅₀ = 0.89 and 1.81 μM, respectively) and antiproliferative activities (IC₅₀ = 7.44 and 6.06 μM, respectively) against the two tested atypical teratoid/rhabdoid tumour (AT/RT) cells BT12 and BT16. Interestingly, compound 9a was not cytotoxic when tested against non-neoplastic HFF1 cells [IC₅₀ (viability) = 119 μM]. This suggests that an indoleamide scaffold can be fine-tuned to confer a set of derivatives with selective antitubercular and/or antitumour activities.

In this study, we demonstrated that an indoleamide scaffold can be fine-tuned to confer a set of derivatives with selective antitubercular and/or antitumour activities.

Assessing Protein Function Through Structural Similarities with CATH

The functional diversity of proteins is closely related to their differences in sequence and structure. Despite variations in functional sites, global structural similarity is a valuable source of information when assessing potential functional similarities between proteins. The CATH database contains a well-established hierarchical classification of more than 430,000 protein domain structures and nearly 95 million protein domain sequences, with integrated functional annotations for each represented family. The present chapter provides an overview of the main features of CATH with emphasis on exploiting structural similarities to obtain functional information for proteins.

Targeting Beta-Blocker Drug-Drug Interactions with Fibrinogen Blood Plasma Protein: A Computational and Experimental Study

In this work, one of the most prevalent polypharmacology drug–drug interaction events that occurs between two widely used beta-blocker drugs—i.e., acebutolol and propranolol—with the most abundant blood plasma fibrinogen protein was evaluated. Towards that end, molecular docking and Density Functional Theory (DFT) calculations were used as complementary tools. A fibrinogen crystallographic validation for the three best ranked binding-sites shows 100% of conformationally favored residues with total absence of restricted flexibility. From those three sites, results on both the binding-site druggability and ligand transport analysis-based free energy trajectories pointed out the most preferred biophysical environment site for drug–drug interactions. Furthermore, the total affinity for the stabilization of the drug–drug complexes was mostly influenced by steric energy contributions, based mainly on multiple hydrophobic contacts with critical residues (THR22: P and SER50: Q) in such best-ranked site. Additionally, the DFT calculations revealed that the beta-blocker drug–drug complexes have a spontaneous thermodynamic stabilization following the same affinity order obtained in the docking simulations, without covalent-bond formation between both interacting beta-blockers in the best-ranked site. Lastly, experimental ultrasound density and velocity measurements were performed and allowed us to validate and corroborate the computational obtained results.

Barbiturate derivatives for managing multifaceted oncogenic pathways: A mini review

Development and progression of metastasis comprises synchronized erroneous expressions of several composite pathways, which are difficult to manage simultaneously with the representative anticancer molecules. The emergence of the drug resistance and the complex interplay between these pathways further potentiates cancer related complexities. Barbiturates and their derivatives present a commendable anticancer profile by attenuating the cancer manifesting metabolic and enzymatic pathways including, but not limited to matrix metalloproteinases, xanthine oxidase, amino peptidases, histone deacetylases, and Ras/mitogen-activated protein kinase. The derivatization and conjugation of barbiturates with pharmacophores delivers a suitable hybrid profile in containing the anomalous expression of these pathways. The present report presents a succinct collation of the barbiturates and their derivatives in managing the various cancer causing pathways.

Exploration in the Mechanism of Kaempferol for the Treatment of Gastric Cancer Based on Network Pharmacology

Background

Kaempferol is a natural polyphenol in lots of Chinese herbs, which has shown promising treatment for gastric cancer (GC). However, the molecular mechanisms of its action have not been systematically revealed yet. In this work, a network pharmacology approach was used to elucidate the potential mechanisms of kaempferol in the treatment of GC.

Methods

The kaempferol was input into the PharmMapper and SwissTargetPrediction database to get its targets, and the targets of GC were obtained by retrieving the Online Mendelian Inheritance in Man (OMIM) database, MalaCards database, Therapeutic Target Database (TTD), and Coolgen database. The molecular docking was performed to assess the interactions between kaempferol and these targets. Next, the overlap targets of kaempferol and GC were identified for GO and KEGG enrichment analyses. Afterward, a protein-protein interaction (PPI) network was constructed to get the hub targets, and the expression and overall survival analysis of the hub target were investigated. Finally, the overall survival (OS) analysis of hub targets was performed using the Kaplan-Meier Plotter online tool.

Results

A total of 990 genes related to GC and 10 overlapping genes were determined through matching the 24 potential targets of kaempferol with disease-associated genes. The result of molecular docking indicated that kaempferol can bind with these hub targets with good binding scores. These targets were further mapped to 140 GO biological process terms and 11 remarkable pathways. In the PPI network analysis, 3 key targets were identified, including ESR1, EGFR, and SRC. The mRNA and protein expression levels of EGFR and SRC were obviously higher in GC tissues. High expression of these targets was related to poor OS in GC patients.

Conclusions

This study provided a novel approach to reveal the therapeutic mechanisms of kaempferol on GC, which will ease the future clinical application of kaempferol in the treatment of GC.

Intracellular Binding/Unbinding Kinetics of Approved Drugs to Carbonic Anhydrase II Observed by in-Cell NMR

Candidate drugs rationally designed in vitro often fail due to low efficacy in vivo caused by low tissue availability or because of unwanted side effects. To overcome the limitations of in vitro rational drug design, the binding of candidate drugs to their target needs to be evaluated in the cellular context. Here, we applied in-cell NMR to investigate the binding of a set of approved drugs to the isoform II of carbonic anhydrase (CA) in living human cells. Some compounds were originally developed toward other targets and were later found to inhibit CAs. We observed strikingly different dose- and time-dependent binding, wherein some drugs exhibited a more complex behavior than others. Specifically, some compounds were shown to gradually unbind from intracellular CA II, even in the presence of free compound in the external medium, therefore preventing the quantitative formation of a stable protein–ligand complex. Such observations could be correlated to the known off-target binding activity of these compounds, suggesting that this approach could provide information on the pharmacokinetic profiles of lead candidates at the early stages of multitarget drug design.

An up-to-date overview of computational polypharmacology in modern drug discovery

INTRODUCTION

In recent years, computational polypharmacology has gained significant attention to study the promiscuous nature of drugs. Despite tremendous challenges, community-wide efforts have led to a variety of novel approaches for predicting drug polypharmacology. In particular, some rapid advances using machine learning and artificial intelligence have been reported with great success.

AREAS COVERED

In this article, the authors provide a comprehensive update on the current state-of-the-art polypharmacology approaches and their applications, focusing on those reports published after our 2017 review article. The authors particularly discuss some novel, groundbreaking concepts, and methods that have been developed recently and applied to drug polypharmacology studies.

EXPERT OPINION

Polypharmacology is evolving and novel concepts are being introduced to counter the current challenges in the field. However, major hurdles remain including incompleteness of high-quality experimental data, lack of in vitro and in vivo assays to characterize multi-targeting agents, shortage of robust computational methods, and challenges to identify the best target combinations and design effective multi-targeting agents. Fortunately, numerous national/international efforts including multi-omics and artificial intelligence initiatives as well as most recent collaborations on addressing the COVID-19 pandemic have shown significant promise to propel the field of polypharmacology forward.

Understanding Membrane Protein Drug Targets in Computational Perspective

Membrane proteins play crucial physiological roles in vivo and are the major category of drug targets for pharmaceuticals. The research on membrane protein is a significant part in the drug discovery. The biological process is a cycled network, and the membrane protein is a vital hub in the network since most drugs achieve the therapeutic effect via interacting with the membrane protein. In this review, typical membrane protein targets are described, including GPCRs, transporters and ion channels. Also, we conclude network servers and databases that are referring to the drug, drug-target information and their relevant data. Furthermore, we chiefly introduce the development and practice of modern medicines, particularly demonstrating a series of state-of-the-art computational models for the prediction of drug-target interaction containing network-based approach and machine-learningbased approach as well as showing current achievements. Finally, we discuss the prospective orientation of drug repurposing and drug discovery as well as propose some improved framework in bioactivity data, created or improved predicted approaches, alternative understanding approaches of drugs bioactivity and their biological processes.

X-ray Structure-Based Chemoinformatic Analysis Identifies Promiscuous Ligands Binding to Proteins from Different Classes with Varying Shapes

(1) Background: Compounds with multitarget activity are of interest in basic research to explore molecular foundations of promiscuous binding and in drug discovery as agents eliciting polypharmacological effects. Our study has aimed to systematically identify compounds that form complexes with proteins from distinct classes and compare their bioactive conformations and molecular properties. (2) Methods: A large-scale computational investigation was carried out that combined the analysis of complex X-ray structures, ligand binding modes, compound activity data, and various molecular properties. (3) Results: A total of 515 ligands with multitarget activity were identified that included 70 organic compounds binding to proteins from different classes. These multiclass ligands (MCLs) were often flexible and surprisingly hydrophilic. Moreover, they displayed a wide spectrum of binding modes. In different target structure environments, binding shapes of MCLs were often similar, but also distinct. (4) Conclusions: Combined structural and activity data analysis identified compounds with activity against proteins with distinct structures and functions. MCLs were found to have greatly varying shape similarity when binding to different protein classes. Hence, there were no apparent canonical binding shapes indicating multitarget activity. Rather, conformational versatility characterized MCL binding.

The Genome3D Consortium for Structural Annotations of Selected Model Organisms

Genome3D consortium is a collaborative project involving protein structure prediction and annotation resources developed by six world-leading structural bioinformatics groups, based in the United Kingdom (namely Blundell, Murzin, Gough, Sternberg, Orengo, and Jones). The main objective of Genome3D serves as a common portal to provide both predicted models and annotations of proteins in model organisms, using several resources developed by these labs such as CATH-Gene3D, DOMSERF, pDomTHREADER, PHYRE, SUPERFAMILY, FUGUE/TOCATTA, and VIVACE. These resources primarily use SCOP- and/or CATH-based protein domain assignments. Another objective of Genome3D is to compare structural classifications of protein domains in CATH and SCOP databases and to provide a consensus mapping of CATH and SCOP protein superfamilies. CATH/SCOP mapping analyses led to the identification of total of 1429 consensus superfamilies.Currently, Genome3D provides structural annotations for ten model organisms, including Homo sapiens, Arabidopsis thaliana, Mus musculus, Escherichia coli, Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, Plasmodium falciparum, Staphylococcus aureus, and Schizosaccharomyces pombe. Thus, Genome3D serves as a common gateway to each structure prediction/annotation resource and allows users to perform comparative assessment of the predictions. It, thus, assists researchers to broaden their perspective on structure/function predictions of their query protein of interest in selected model organisms.

Applications of molecular networks in biomedicine

Due to the large interdependence between the molecular components of living systems, many phenomena, including those related to pathologies, cannot be explained in terms of a single gene or a small number of genes. Molecular networks, representing different types of relationships between molecular entities, embody these large sets of interdependences in a framework that allow their mining from a systemic point of view to obtain information. These networks, often generated from high-throughput omics datasets, are used to study the complex phenomena of human pathologies from a systemic point of view. Complementing the reductionist approach of molecular biology, based on the detailed study of a small number of genes, systemic approaches to human diseases consider that these are better reflected in large and intricate networks of relationships between genes. These networks, and not the single genes, provide both better markers for diagnosing diseases and targets for treating them. Network approaches are being used to gain insight into the molecular basis of complex diseases and interpret the large datasets associated with them, such as genomic variants. Network formalism is also suitable for integrating large, heterogeneous and multilevel datasets associated with diseases from the molecular level to organismal and epidemiological scales. Many of these approaches are available to nonexpert users through standard software packages.

Peptide Conjugates with Small Molecules Designed to Enhance Efficacy and Safety

Peptides constitute molecular diversity with unique molecular mechanisms of action that are proven indispensable in the management of many human diseases, but of only a mere fraction relative to more traditional small molecule-based medicines. The integration of these two therapeutic modalities offers the potential to enhance and broaden pharmacology while minimizing dose-dependent toxicology. This review summarizes numerous advances in drug design, synthesis and development that provide direction for next-generation research endeavors in this field. Medicinal studies in this area have largely focused upon the application of peptides to selectively enhance small molecule cytotoxicity to more effectively treat multiple oncologic diseases. To a lesser and steadily emerging extent peptides are being therapeutically employed to complement and diversify the pharmacology of small molecule drugs in diseases other than just cancer. No matter the disease, the purpose of the molecular integration remains constant and it is to achieve superior therapeutic outcomes with diminished adverse effects. We review linker technology and conjugation chemistries that have enabled integrated and targeted pharmacology with controlled release. Finally, we offer our perspective on opportunities and obstacles in the field.

Genie in a bottle: controlled release helps tame natural polypharmacology?

Ability to faithfully report drug-target interactions constitutes a major critical parameter in preclinical/clinical settings. Yet the assessment of target engagement remains challenging, particularly for promiscuous and/or polypharmacologic ligands. Drawing from our improved insights into native electrophile signaling and emerging technologies that profile and interrogate these non-enzyme-assisted signaling subsystems, we posit that 'trained' polypharmocologic covalent inhibitors can be designed. Accumulating evidence indicates that electrophile-modified states at fractional occupancy can alter cell fate. Thus, by understanding sensing preferences and ligandable regions favored by the natural electrophilic signals at individual protein-ligand resolution, we can better evaluate target engagement and develop a function-guided understanding of polypharmacology.

Promiscuous Ligands from Experimentally Determined Structures, Binding Conformations, and Protein Family-Dependent Interaction Hotspots

Compound promiscuity is often attributed to nonspecific binding or assay artifacts. On the other hand, it is well-known that many pharmaceutically relevant compounds are capable of engaging multiple targets in vivo, giving rise to polypharmacology. To explore and better understand promiscuous binding characteristics of small molecules, we have searched X-ray structures (and very few qualifying solution structures) for ligands that bind to multiple distantly related or unrelated target proteins. Experimental structures of a given ligand bound to different targets represent high-confidence data for exploring promiscuous binding events. A total of 192 ligands were identified that formed crystallographic complexes with proteins from different families and for which activity data were available. These "multifamily" compounds included endogenous ligands and were often more polar than other bound compounds and active in the submicromolar range. Unexpectedly, many promiscuous ligands displayed conserved or similar binding conformations in different active sites. Others were found to conformationally adjust to binding sites of different architectures. A comprehensive analysis of ligand-target interactions revealed that multifamily ligands frequently formed different interaction hotspots in binding sites, even if their bound conformations were similar, thus providing a rationale for promiscuous binding events at the molecular level of detail. As a part of this work, all multifamily ligands we have identified and associated activity data are made freely available.