Structure-based virtual screening and biological evaluation of novel small-molecule BTK inhibitors

Bruton tyrosine kinase (BTK) is linked to multiple signalling pathways that regulate cellular survival, activation, and proliferation. A covalent BTK inhibitor has shown favourable outcomes for treating B cell malignant leukaemia. However, covalent inhibitors require a high reactive warhead that may contribute to unexpected toxicity, poor selectivity, or reduced effectiveness in solid tumours. Herein, we report the identification of a novel noncovalent BTK inhibitor. The binding interactions (i.e. interactions from known BTK inhibitors) for the BTK binding site were identified and incorporated into a structure-based virtual screening (SBVS). Top-rank compounds were selected and testing revealed a BTK inhibitor with >50% inhibition at 10 µM concentration. Examining analogues revealed further BTK inhibitors. When tested across solid tumour cell lines, one inhibitor showed favourable inhibitory activity, suggesting its potential for targeting BTK malignant tumours. This inhibitor could serve as a basis for developing an effective BTK inhibitor targeting solid cancers.     

Identification of novel small-molecule histone deacetylase inhibitors by medium-throughput screening using a fluorigenic assay

Cationic peptides, known to disrupt bacterial membranes, are being developed as promising agents for therapeutic intervention against infectious disease. In the present study, we investigate structure-activity relationships in the bacterial membrane disruptor betapep-25, a peptide 33-mer. For insight into which amino acid residues are functionally important, we synthesized alanine-scanning variants of betapep-25 and assessed their ability to kill bacteria (Escherichia coli, Pseudomonas aeruginosa and Staphylococcus aureus) and to neutralize LPS (lipopolysaccharide). Activity profiles were found to vary with the bacterial strain examined. Specific cationic and smaller hydrophobic alkyl residues were crucial to optimal bactericidal activity against the Gram-negative bacteria, whereas larger hydrophobic and cationic residues mediated optimal activity against Gram-positive Staph. aureus. Lysine-substituted norleucine (n-butyl group) variants demonstrated that both charge and alkyl chain length mediate optimal activity. In terms of LPS neutralization, activity profiles were essentially the same against four species of LPS (E. coli 055 and 0111, Salmonella enterica serotype Typhimurium and Klebsiella pneumoniae), and different for two others (Ps. aeruginosa and Serratia marcescens), with specific hydrophobic, cationic and, surprisingly, anionic residues being functionally important. Furthermore, disulfide-bridged analogues demonstrated that an anti parallel beta-sheet structure is the bioactive conformation of betapep-25 in terms of its bactericidal, but not LPS endotoxin neutralizing, activity. Moreover, betapep-25 variants, like the parent peptide, do not lyse eukaryotic cells. This research contributes to the development and design of novel antibiotics.     

Identification of small-molecule urea derivatives as novel NAMPT inhibitors via pharmacophore-based virtual screening

Nicotinamide phosphoribosyltransferase (NAMPT) catalyzes the condensation of nicotinamide (NAM) with 5-phosphoribosyl-1-prophosphate (PRPP) to yield nicotinamide mononucleotide (NMN), a rate limiting enzyme in a mammalian salvage pathway of nicotinamide adenine dinucleotide (NAD⁺) synthesis. Recently, intracellular NAD⁺ has received substantial attention due to the recent discovery that several enzymes including poly(ADP-ribose) polymerases (PARPs), mono(ADP-ribose) transferases (ARTs), and sirtuins (SIRTs), use NAD⁺ as a substrate, suggesting that intracellular NAD⁺ level may regulate cytokine production, metabolism, and aging through these enzymes. NAMPT is found to be upregulated in various types of cancer, and given its importance in the NAD⁺ salvage pathway, NAMPT is considered as an attractive target for the development of new cancer therapies. In this study, the reported NAMPT inhibitors bearing amide, cyanoguanidine, and urea scaffolds were used to generate pharmacophore models and pharmacophore-based virtual screening studies were performed against ZINC database. Following the filtering steps, ten hits were identified and evaluated for their in vitro NAMPT inhibitory effects. Compounds GF4 (NAMPT IC₅₀ = 2.15 ± 0.22 μM) and GF8 (NAMPT IC₅₀ = 7.31 ± 1.59 μM) were identified as new urea-typed inhibitors of NAMPT which also displayed cytotoxic activities against human HepG2 hepatocellular carcinoma cell line with IC₅₀ values of 15.20 ± 1.28 and 24.28 ± 6.74 μM, respectively.     

Development of a high-throughput screening assay for the discovery of small-molecule SecA inhibitors

A major pathway for bacterial preprotein translocation is provided by the Sec-dependent preprotein translocation pathway. Proteins destined for Sec-dependent translocation are synthesized as preproteins with an N-terminal signal peptide, which targets them to the SecYEG translocase channel. The driving force for the translocation reaction is provided by the peripheral membrane ATPase SecA, which couples the hydrolysis of ATP to the stepwise transport of unfolded preproteins across the bacterial membrane. Since SecA is essential, highly conserved among bacterial species, and has no close human homologues, it represents a promising target for antibacterial chemotherapy. However, high-throughput screening (HTS) campaigns to identify SecA inhibitors are hampered by the low intrinsic ATPase activity of SecA and the requirement of hydrophobic membranes for measuring the membrane or translocation ATPase activity of SecA. To address this issue, we have developed a colorimetric high-throughput screening assay in a 384-well format, employing an Escherichia coli (E. coli) SecA mutant with elevated intrinsic ATPase activity. The assay was applied for screening of a chemical library consisting of ∼27,000 compounds and proved to be highly reliable (average Z′ factor of 0.89). In conclusion, a robust HTS assay has been established that will facilitate the search for novel SecA inhibitors.     

A three-fluorophore FRET assay for high-throughput screening of small-molecule inhibitors of ribosome assembly

In one of the first steps of prokaryotic ribosome assembly, the ribosomal protein S15 binds to a three-way junction in the central domain of the 16S rRNA. Binding causes a conformational change that is required for subsequent binding events. Using a novel fluorescence resonance energy transfer assay with three fluorophores, two on the RNA and one on the S15 protein, small-molecule libraries can be screened for potential inhibitors of this initial step in ribosome assembly. The employment of three fluorophores allows both the conformational change of the RNA and the binding of S15 to be monitored in a single assay.     

Discovery of novel small-molecule Src kinase inhibitors via a kinase-focused druglikeness rule and structure-based virtual screening

Overexpression of the non-receptor tyrosine kinase Src is implicated in the development and progression of various human cancers. Blocking signalling pathways mediated by Src is believed to be a promising anticancer strategy. We report herein the discovery of novel small-molecule Src inhibitors by crystal structure-based virtual screening. A kinase-focused druglikeness rule was proposed and used in the design of compound library. Combination of large-scale docking with DOCK and rescoring with GOLD resulted in 6 hits with moderate to potent inhibitory activity against Src. Among them, compound 1 with an IC₅₀ of 1.2 μM shows the most potent inhibitory activity. By using molecular docking, binding models of the top 3 hits (ranked by potency and ligand efficiency) with Src were constructed to provide a rational strategy that simultaneously exploits hydrogen bonding interaction in the hinge region and hydrophobic stacking in the back pocket. This approach is instructive and meaningful to further structure-based drug design of Src inhibitors.     

Inhibition of protein-protein interactions: the discovery of druglike beta-catenin inhibitors by combining virtual and biophysical screening

The interaction between beta-catenin and Tcf family members is crucial for the Wnt signal transduction pathway, which is commonly mutated in cancer. This interaction extends over a very large surface area (4800 A(2)), and inhibiting such interactions using low molecular weight inhibitors is a challenge. However, protein surfaces frequently contain "hot spots," small patches that are the main mediators of binding affinity. By making tight interactions with a hot spot, a small molecule can compete with a protein. The Tcf3/Tcf4-binding surface on beta-catenin contains a well-defined hot spot around residues K435 and R469. A 17,700 compounds subset of the Pharmacia corporate collection was docked to this hot spot with the QXP program; 22 of the best scoring compounds were put into a biophysical (NMR and ITC) screening funnel, where specific binding to beta-catenin, competition with Tcf4 and finally binding constants were determined. This process led to the discovery of three druglike, low molecular weight Tcf4-competitive compounds with the tightest binder having a K(D) of 450 nM. Our approach can be used in several situations (e.g., when selecting compounds from external collections, when no biochemical functional assay is available, or when no HTS is envisioned), and it may be generally applicable to the identification of inhibitors of protein-protein interactions.     

Discovery and optimization of novel small-molecule HIV-1 entry inhibitors using field-based virtual screening and bioisosteric replacement

With the emergence of drug-resistant strains and the cumulative toxicities associated with current therapies, demand remains for new inhibitors of HIV-1 replication. The inhibition of HIV-1 entry is an attractive, yet underexploited therapeutic approach with implications for salvage and preexposure prophylactic regimens, as well as topical microbicides. Using the combination of a field-derived bioactive conformation template to perform virtual screening and iterative bioisosteric replacements, coupled with in silico predictions of absorption, distribution, metabolism, and excretion, we have identified new leads for HIV-1 entry inhibitors.     

Virtual screening based identification of novel small-molecule inhibitors targeted to the HIV-1 capsid

The hydrophobic cavity of the C-terminal domain (CTD) of HIV-1 capsid has been recently validated as potential target for antiviral drugs by peptide-based inhibitors; however, there is no report yet of any small molecule compounds that target this hydrophobic cavity. In order to fill this gap and discover new classes of ant-HIV-1 inhibitors, we undertook a docking-based virtual screening and subsequent analog search, and medicinal chemistry approaches to identify small molecule inhibitors against this target. This article reports for the first time, to the best of our knowledge, identification of diverse classes of inhibitors that efficiently inhibited the formation of mature-like viral particles verified under electron microscope (EM) and showed potential as anti-HIV-1 agents in a viral infectivity assay against a wide range of laboratory-adapted as well as primary isolates in MT-2 cells and PBMC. In addition, the virions produced after the HIV-1 infected cells were treated with two of the most active compounds showed drastically reduced infectivity confirming the potential of these compounds as anti-HIV-1 agents. We have derived a comprehensive SAR from the antiviral data. The SAR analyses will be useful in further optimizing the leads to potential anti-HIV-1 agents.     

TR-FRET-based high-throughput screening assay for identification of UBC13 inhibitors

UBC13 is a noncanonical ubiquitin conjugating enzyme (E2) that has been implicated in a variety of cellular signaling processes due to its ability to catalyze formation of lysine 63-linked polyubiquitin chains on various substrates. In particular, UBC13 is required for signaling by a variety of receptors important in immune regulation, making it a candidate target for inflammatory diseases. UBC13 is also critical for double-strand DNA repair and thus a potential radiosensitizer and chemosensitizer target for oncology. The authors developed a high-throughput screening (HTS) assay for UBC13 based on the method of time-resolved fluorescence resonance energy transfer (TR-FRET). The TR-FRET assay combines fluorochrome (Fl)-conjugated ubiquitin (fluorescence acceptor) with terbium (Tb)-conjugated ubiquitin (fluorescence donor), such that the assembly of mixed chains of Fl- and Tb-ubiquitin creates a robust TR-FRET signal. The authors defined conditions for optimized performance of the TR-FRET assay in both 384- and 1536-well formats. Chemical library screens (total 456 865 compounds) were conducted in high-throughput mode using various compound collections, affording superb Z' scores (typically >0.7) and thus validating the performance of the assays. Altogether, the HTS assays described here are suitable for large-scale, automated screening of chemical libraries in search of compounds with inhibitory activity against UBC13.     

Identification of novel small-molecule inhibitors for human transketolase by high-throughput screening with fluorescent intensity (FLINT) assay

The metabolic enzyme transketolase (TK) plays a crucial role in tumor cell nucleic acid synthesis, using glucose through the elevated nonoxidative pentose phosphate pathway (PPP). Identification of inhibitors specifically targeting TK and preventing the nonoxidative PPP from generating the RNA ribose precursor, ribose-5-phosphate, provides a novel approach for developing effective anticancer therapeutic agents. The full-length human transketolase gene was cloned and expressed in Escherichia coli and the recombinant human transketolase protein purified to homogeneity. A fluorescent intensity (FLINT) assay was developed and optimized. Library compounds were screened in a high-throughput screening (HTS) campaign using the FLINT assay. Fifty-four initial hits were identified. Among them, 2 scaffolds with high selectivity, ideal physiochemical properties, and low molecular weight were selected for lead optimization studies. These compounds specifically inhibited in vitro TK enzyme activity and suppressed tumor cell proliferation in at least 3 cancer cell lines: SW620, LS174T, and MIA PaCa-2. Identification of these active scaffolds represents a good starting point for development of drugs specifically targeting TK and the nonoxidative PPP for cancer therapy.     

A live-cell high-throughput screening assay for identification of fatty acid uptake inhibitors

We developed a live-cell high-throughput assay system using the baker's yeast Saccharomyces cerevisiae to screen for chemical compounds that will inhibit fatty acid uptake. The target for the inhibitors is a mammalian fatty acid transport protein (mmFATP2), which is involved in the fatty acid transport and activation pathway. The mmFATP2 was expressed in a S. cerevisiae mutant strain deficient in Fat1p-dependent fatty acid uptake and reduced in long-chain fatty acid activation, fat1Deltafaa1Delta. To detect fatty acid import, a fluorescent fatty acid analog, 4,4-difluoro-5-methyl-4-bora-3a,4a-diaza-s-indacene-3-dodecanoic acid (C1-BODIPY-C12), was incubated with cells expressing FATP2 in a 96-well plate. The mmFATP2-dependent C1-BODIPY-C12 uptake was monitored by measuring intracellular C1-BODIPY-C12 fluorescence on a microtiter plate reader, whereas extracellular fluorescence was quenched by a cell viability dye, trypan blue. Using this high-throughput screening method, we demonstrate that the uptake of the fluorescent fatty acid ligand was effectively competed by the natural fatty acid oleate. Inhibition of uptake was also demonstrated to occur when cells were pretreated with sodium azide or Triacsin C. This yeast live-cell-based assay is rapid to execute, inexpensive to implement, and has adequate sensitivity for high-throughput screening. The assay basis and limitations are discussed.     

Identification of small-molecule compounds targeting the dishevelled PDZ domain by virtual screening and binding studies

The Dishevelled (Dvl) protein, which conveys signals from receptors to the downstream effectors, is a critical constituent of the Wnt/β-catenin signaling pathway. Because the PDZ domain of Dvl protein functions through associations with a wide range of protein partners, Dvl protein involved in the Wnt signaling pathway has been considered to be therapeutic targets in cancers. In this study, we performed structure-based pharmacophore model of the Dvl PDZ domain to discover novel small-molecule binders and identified eight compounds with micromolar affinity. The most potent compound identified, BMD4702, efficiently bound to the Dvl PDZ domain with 11.2 μM affinity and had a 0.186 μM K_D value according to surface plasmon resonance and fluorescence spectroscopy, respectively. Combining both structural–kinetic relationship analyses and docking studies, we fourmulated that the ligand-binding site is composed of three H-bonds and three hydrophobic features. Thus, our approach led to the identification of potent binders of the Dvl PDZ domain and the findings provide novel insights into structure-based approaches to design high-affinity binders for the Dvl PDZ domain.     

High-throughput luciferase reporter assay for small-molecule inhibitors of microRNA function

MicroRNAs (miRNAs) are endogenous, single-stranded, noncoding RNAs of 21 to 23 nucleotides that regulate gene expression, typically by binding the 3' untranslated regions of target messenger RNAs. It is estimated that miRNAs are involved in the regulation of 30% of all genes and almost every genetic pathway. Recently, the misregulation of miRNAs has been linked to various human diseases including cancer and viral infections, identifying miRNAs as potential targets for drug discovery. Thus, small-molecule modifiers of miRNAs could serve as lead structures for the development of new therapeutic agents and be useful tools in the elucidation of detailed mechanisms of miRNA function. As a result, we have developed a high-throughput screen for potential small-molecule regulators of the liver-specific microRNA miR-122, which is involved in hepatocellular carcinoma development and hepatitis C virus infection. Our small-molecule screen employs a Huh7 human hepatoma cell line stably transfected with a Renilla luciferase sensor for endogenous miR-122. The assay was optimized and validated using an miR-122 antisense agent and a previously identified small-molecule miR-122 inhibitor. The described reporter assay will enable the high-throughput screening of small-molecule miR-122 inhibitors and can be readily extended to other miRNAs.     

Combinatorial enzymatic assay for the screening of a new class of bacterial cell wall inhibitors

We have developed a screening assay by thin-layer chromatography (TLC) to identify inhibitors for the bacterial essential enzymes MurA, -B, and -C. Libraries of compounds were synthesized using the mix-and-split combinatorial chemistry approach. Screening of the pooled compounds using the developed assay revealed the presence of many pools active in vitro. Pools of interest were tested for antibacterial activity. Individual molecules in the active pools were synthesized and retested with the TLC assay and with bacteria. We focused on the best five compounds for further analysis. They were tested for inhibition on each of the three enzymes separately, and showed no inhibition of MurA or MurB activity but were all inhibitors of MurC enzyme. This approach yielded interesting lead compounds for the development of novel antibacterial agents.     

Ligand-based pharmacophore mapping and virtual screening for identification of potential discoidin domain receptor 1 inhibitors

Abbreviations ADME absorption, distribution, metabolism, and excretion

MMGB/SA molecular mechanics generalized born surface area

IFD induced fit docking

RTK receptor tyrosine kinase

NSCLC non-small-cell lung cancer

ATP adenosine triphosphate

OPLS optimized potential for liquid stimulation

RMSD root mean square deviation

HTVS high-throughput virtual screening

SP standard precision

XP extra precision

OPLS-AA optimized potential for liquid stimulation-all atom

MD molecular simulation

MME molecular mechanics energies

SGB surface generalized born

POPC membrane 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine membrane

PDB Protein Data Bank

DDR1 discoidin domain receptor 1

DDR2 discoidin domain receptor 2

DDRs discoidin domain receptors

ECM extracellular matrix

TIP4P transferable intermolecular potential 4 point

NPT constant particle number, pressure and temperature

RMSF root mean square fluctuation

Communicated by Ramaswamy H. Sarma     

3D QSAR pharmacophore-based virtual screening for the identification of potential inhibitors of tyrosinase

Abstract

Tyrosinase plays an important role in melanin biosynthesis and protects skin against ultraviolet radiations. Functional deficiency of tyrosinase results in serious dermatological diseases. Tyrosinase also participates in neuromelanin formation in the human brain, which leads to neurodegeneration resulting in Parkinson’s disease. In fruits and vegetables, tyrosinase plays a critical role in senescence, causing undesired browning that results in faster deterioration and shorter shelf lines. The only commercially available tyrosinase is mushroom tyrosinase and it shows the highest homology to the mammalian tyrosinase. Although kojic acid is currently used as a tyrosinase inhibitor, they have serious side effects such as dermatitis, carcinogenesis and hepatotoxicity. Therefore, in order to develop a more active and safer tyrosinase inhibitor, 3D QSAR pharmacophore models were generated based on experimentally known inhibitors. The pharmacophore model, Hypo1, was developed with a large cost difference, high correlation coefficient and low RMS deviation. Hypo1 showed a good spatial arrangement; consisting of five-point features including two hydrogen bond acceptor, one hydrogen bond donor and two hydrophobic features. Hypo1 was further validated by cost analysis, test set and Fisher’s randomisation method. Hypo1 was used as a 3D query for screening the in-house drug-like databases, and the hits were further selected by applying ADMET, Lipinski’s rule of five and fit value criteria. To identify binding conformations, the obtained hits were subjected to molecular docking. Finally, molecular dynamics simulations revealed the appropriate binding modes of hit compounds. To conclude, we propose the final three hit compounds with new structural scaffolds as a virtual candidate as tyrosinase inhibitors.

Communicated by Ramaswamy H. Sarma