Selective inhibitors of human mPGES-1 from structure-based computational screening

Human mPGES-1 is recognized as a promising target for next generation of anti-inflammatory drugs. Although various mPGES-1 inhibitors have been reported in literature, few have entered clinical trials and none has been proven clinically useful so far. It is highly desired for developing the next generation of therapeutics for inflammation-related diseases to design and discover novel inhibitors of mPGES-1 with new scaffolds. Here, we report the identification of a series of new, potent and selective inhibitors of human mPGES-1 with diverse scaffolds through combined computational and experimental studies. The computationally modeled binding structures of these new inhibitors of mPGES-1 provide some interesting clues for rational design of modified structures of the inhibitors to more favorably bind with mPGES-1.     

Fragment-based discovery of novel and selective mPGES-1 inhibitors Part 1: Identification of sulfonamido-1,2,3-triazole-4,5-dicarboxylic acid

Microsomal prostaglandin E synthase-1 (mPGES-1) is an inducible prostaglandin E synthase that catalyzes the conversion of prostaglandin PGH₂ to PGE₂ and represents a novel target for therapeutic treatment of inflammatory disorders. It is essential to identify mPGES-1 inhibitor with novel scaffold as new hit or lead compound for the purpose of the next-generation anti-inflammatory drugs. Herein we report the discovery of sulfonamido-1,2,3-triazole-4,5-dicarboxylic derivatives as a novel class of mPGES-1 inhibitors identified through fragment-based virtual screening and in vitro assays on the inhibitory activity of the actual compounds. 1-[2-(N-Phenylbenzenesulfonamido)ethyl]-1H-1,2,3-triazole-4,5-dicarboxylic acid (6f) inhibits human mPGES-1 (IC₅₀ of 1.1 μM) with high selectivity (ca.1000-fold) over both COX-1 and COX-2 in a cell-free assay. In addition, the activity of compound 6f was again tested at 10 μM concentration in presence of 0.1% Triton X-100 and found to be reduced to 1/4 of its original activity without this detergent. Compared to the complete loss of activity of nuisance inhibitor with the detergent, therefore, compound 6f would be regarded as a partial nuisance inhibitor of mPGES-1 with a novel scaffold for the optimal design of more potent mPGES-1 inhibitors.     

Discovery and development of heat shock protein 90 inhibitors

Heat shock protein 90 (Hsp90) is an important target in cancer because of its role in maintaining transformation and has recently become the focus of several drug discovery and development efforts. While compounds with different modes of action are known, the focus of this review is on those classes of compounds which inhibit Hsp90 by binding to the N-terminal ATP pocket. These include natural product inhibitors such as geldanamycin and radicicol and synthetic inhibitors comprised of purines, pyrazoles, isoxazoles and other scaffolds. The synthetic inhibitors have been discovered either by structure-based design, high throughput screening and more recently using fragment-based design and virtual screening techniques. This review will discuss the discovery of these different classes, as well as their development as potential clinical agents.     

Structure-based and ligand-based drug design for microsomal prostaglandin E synthase-1 inhibitors

Microsomal prostaglandin E synthase-1 (mPGES-1) has been regarded as an attractive drug for inflammation-related diseases. In search of new mPGES-1 inhibitors, we performed virtual screening using our traditional Chinese medicine and natural products database (http://tcm.cmu.edu.tw/) and constructed comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) using a training set of 30 experimentally tested mPGES-1 inhibitors. The CoMFA and CoMSIA models derived were statistically significant with cross-validated coefficient values of 0.808 for CoMFA and 0.829 for CoMSIA and non-cross-validated coefficient values of 0.829 for CoMFA and 0.980 for CoMSIA. Docking and de novo evolution design gave three top derivatives, 2-O-caffeoyl tartaric acid-Evo_2, glucogallin-Evo_1 and 3-O-feruloylquinic acid-Evo_7 that have higher binding affinities than the control, glutathione. These three derivatives have interactions with Arg70, Arg73, Arg110, Arg126 and Arg38, which all are mPGES-1 key active site residues. In addition, these derivatives fit well into the CoMFA and CoMSIA models, with hydrophobic, hydrophilic and electropositive substructures mapped onto corresponding contour plots. Hence, we suggest that these three de novo compounds could be a starting basis for new mPGES-1 inhibitors.     

Identification of PPARgamma Partial Agonists of Natural Origin (I): Development of a Virtual Screening Procedure and In Vitro Validation

Background

Although there are successful examples of the discovery of new PPARγ agonists, it has recently been of great interest to identify new PPARγ partial agonists that do not present the adverse side effects caused by PPARγ full agonists. Consequently, the goal of this work was to design, apply and validate a virtual screening workflow to identify novel PPARγ partial agonists among natural products.

Methodology/Principal Findings

We have developed a virtual screening procedure based on structure-based pharmacophore construction, protein-ligand docking and electrostatic/shape similarity to discover novel scaffolds of PPARγ partial agonists. From an initial set of 89,165 natural products and natural product derivatives, 135 compounds were identified as potential PPARγ partial agonists with good ADME properties. Ten compounds that represent ten new chemical scaffolds for PPARγ partial agonists were selected for in vitro biological testing, but two of them were not assayed due to solubility problems. Five out of the remaining eight compounds were confirmed as PPARγ partial agonists: they bind to PPARγ, do not or only moderately stimulate the transactivation activity of PPARγ, do not induce adipogenesis of preadipocyte cells and stimulate the insulin-induced glucose uptake of adipocytes.

Conclusions/Significance

We have demonstrated that our virtual screening protocol was successful in identifying novel scaffolds for PPARγ partial agonists.     

An in silico approach for the discovery of CDK5/p25 interaction inhibitors

The lack of selectivity of all existing ATP competitive inhibitors for a single cyclin-dependent kinase (CDK) has led us to redirect the structure-based molecule design from targeting the classic ATP-binding pocket in CDK5 toward the CDK5/p25 interface. The aim was to seek novel inhibition mechanisms to interrupt protein-protein interactions. A combined strategy of alanine-scanning calculations for locating binding sites, virtual screening for small molecules, molecular dynamics simulations for examining the binding stability of virtual screening hits and bio-assays for testing the level of inhibition was set up and used to explore novel inhibitors capable of interrupting the interactions between the proteins, and consequently of inhibiting the kinase activity. Two compounds were shown to inhibit the complex formation between CDK5 and p25 through p25 binding. They could open avenues for the discovery of new types of structures that prevent interactions between CDK5 and p25 or other CDK and activator proteins, and, more importantly, provide leads in the development of selective inhibitors among CDKs.     

A novel class of Hsp90 inhibitors isolated by structure-based virtual screening

A novel class of 3-phenyl-2-styryl-3H-quinazolin-4-one Hsp90 inhibitors with in vitro anti-tumor activity are identified by structure-based virtual screening of a chemical database with docking simulations in the N-terminal ATP-binding site, in vitro ATPase assay using yeast Hsp90, and cell-based Her2 degradation assay in a consecutive fashion. These results exemplify the usefulness of the structure-based virtual screening with molecular docking in drug discovery. The structural features responsible for a tight binding of the inhibitors in the active site of Hsp90 are discussed in detail.     

Discovering novel ligands for macromolecules using X-ray crystallographic screening

The need to decrease the time scale for clinical compound discovery has led to innovations at several stages in the process, including genomics/proteomics for target identification, ultrahigh-throughput screening for lead identification, and structure-based drug design and combinatorial chemistry for lead optimization. A critical juncture in the process is the identification of a proper lead compound, because a poor choice may generate costly difficulties at later stages. Lead compounds are commonly identified from high-throughput screens of large compound libraries, derived from known substrates/inhibitors, or identified in computational prescreeusing X-ray crystal structures. Structural information is often consulted to efficiently optimize leads, but under the current paradigm, such data require preidentification and confirmation of compound binding. Here, we describe a new X-ray crystallography-driven screening technique that combines the steps of lead identification, structural assessment, and optimization. The method is rapid, efficient, and high-throughput, and it results in detailed crystallographic structure information. The utility of the method is demonstrated in the discovery and optimization of a new orally available class of urokinase inhibitors for the treatment of cancer.     

Joint effects of acetochlor and urea on germinating characteristics of crop seeds

Neuraminidase (NA) is one of the most important targets to screen the drugs of anti-influenza virus A and B. After virtual screening approaches were applied to a compound database which possesses more than 10000 compound structures, 160 compounds were selected for bioactivity assay, then a High Throughput Screening (HTS) model established for influenza virus NA inhibitors was applied to detect these compounds. Finally, three compounds among them displayed higher inhibitory activities, the range of their IC5o was from 0.1 μmol/L to 3 μmol/L. Their structural scaffolds are novel and different from those of NA inhibitors approved for influenza treatment, and will be useful for the design and research of new NA inhibitors. The result indicated that the combination of virtual screening with HTS was very significant to drug screening and drug discovery.     

Selective immunoproteasome inhibitors with non-peptide scaffolds identified from structure-based virtual screening

As a major component of the crucial nonlysosomal protein degradation pathway in the cells, the proteasome has been implicated in many diseases such as Alzheimer’s disease, Huntington’s disease, inflammatory bowel diseases, autoimmune diseases, multiple myeloma (MM) and other cancers. There are two main proteasome subtypes: the constitutive proteasome which is expressed in all eukaryotic cells and the immunoproteasome which is expressed in immune cells and can be induced in other cell types. Majority of currently available proteasome inhibitors are peptide backbone-based, having short half-lives in the body. It is highly desirable to identify novel, immunoproteasome-selective inhibitors with non-peptide scaffolds for development of novel therapeutics. Through combined virtual screening and experimental studies targeting the immunoproteasome, we have identified a set of novel immunoproteasome inhibitors with diverse non-peptide scaffolds. Some of the identified inhibitors have significant selectivity for the immunoproteasome over the constitutive proteasome. Unlike most of the currently available proteasome inhibitors, these new inhibitors lacking electrophilic pharmacophores are not expected to form a covalent bond with proteasome after the binding. These non-peptide scaffolds may provide a new platform for future rational drug design and discovery targeting the immunoproteasome.     

Drug screening for influenza neuraminidase inhibitors

Neuraminidase (NA) is one of the most important targets to screen the drugs of anti-influenza virus A and B. After virtual screening approaches were applied to a compound database which possesses more than 10000 compound structures, 160 compounds were selected for bioactivity assay, then a High Throughput Screening (HTS) model established for influenza virus NA inhibitors was applied to detect these compounds. Finally, three compounds among them displayed higher inhibitory activities, the range of their IC50 was from 0.1 μmol/L to 3μmol/L. Their structural scaffolds are novel and different from those of NA inhibitors approved for influenza treatment, and will be useful for the design and research of new NA inhibitors. The resuit indicated that the combination of virtual screening with HTS was very significant to drug screening and drug discovery.     

Structure-based virtual screening against SARS-3CL(pro) to identify novel non-peptidic hits

Severe acute respiratory syndrome is a highly infectious upper respiratory tract disease caused by SARS-CoV, a previously unidentified human coronavirus. SARS-3CL(pro) is a viral cysteine protease critical to the pathogen's life cycle and hence a therapeutic target of importance. The recently elucidated crystal structures of this enzyme provide an opportunity for the discovery of inhibitors through rational drug design. In the current study, Gold docking program was utilized to conduct extensive docking studies against the target crystal structure to develop a robust and predictive docking protocol. The validated docking protocol was used to conduct a structure-based virtual screening of the Asinex Platinum collection. Biological evaluation of a screened selection of compounds was carried out to identify novel inhibitors of the viral protease.     

Biarylimidazoles as inhibitors of microsomal prostaglandin E2 synthase-1

Microsomal prostaglandin E(2) synthase (mPGES-1) represents a potential target for novel analgesic and anti-inflammatory agents. High-throughput screening identified several leads of mPGES-1 inhibitors which were further optimized for potency and selectivity. A series of inhibitors bearing a biaryl imidazole scaffold exhibits excellent inhibition of PGE(2) production in enzymatic and cell-based assays. The synthesis of these molecules and their activities will be discussed.     

Structure-based and shape-complemented pharmacophore modeling for the discovery of novel checkpoint kinase 1 inhibitors

Checkpoint kinase 1 (Chk1), a member of the serine/threonine kinase family, is an attractive therapeutic target for anticancer combination therapy. A structure-based modeling approach complemented with shape components was pursued to develop a reliable pharmacophore model for ATP-competitive Chk1 inhibitors. Common chemical features of the pharmacophore model were derived by clustering multiple structure-based pharmacophore features from different Chk1-ligand complexes in comparable binding modes. The final model consisted of one hydrogen bond acceptor (HBA), one hydrogen bond donor (HBD), two hydrophobic (HY) features, several excluded volumes and shape constraints. In the validation study, this feature-shape query yielded an enrichment factor of 9.196 and performed fairly well at distinguishing active from inactive compounds, suggesting that the pharmacophore model can serve as a reliable tool for virtual screening to facilitate the discovery of novel Chk1 inhibitors. Besides, these pharmacophore features were assumed to be essential for Chk1 inhibitors, which might be useful for the identification of potential Chk1 inhibitors.     

Discovery of thienoquinolone derivatives as selective and ATP non-competitive CDK5/p25 inhibitors by structure-based virtual screening

Calpain mediated cleavage of CDK5 natural precursor p35 causes a stable complex formation of CDK5/p25, which leads to hyperphosphorylation of tau. Thus inhibition of this complex is a viable target for numerous acute and chronic neurodegenerative diseases involving tau protein, including Alzheimer’s disease. Since CDK5 has the highest sequence homology with its mitotic counterpart CDK2, our primary goal was to design selective CDK5/p25 inhibitors targeting neurodegeneration. A novel structure-based virtual screening protocol comprised of e-pharmacophore models and virtual screening workflow was used to identify nine compounds from a commercial database containing 2.84 million compounds. An ATP non-competitive and selective thieno[3,2-c]quinolin-4(5H)-one inhibitor (10) with ligand efficiency (LE) of 0.3 was identified as the lead molecule. Further SAR optimization led to the discovery of several low micromolar inhibitors with good selectivity. The research represents a new class of potent ATP non-competitive CDK5/p25 inhibitors with good CDK2/E selectivity.     

On scaffold hopping: Challenges in the discovery of sulfated small molecules as mimetics of glycosaminoglycans

The design of sulfated, small, nonsaccharide molecules as modulators of proteins is still in its infancy as standard drug discovery tools such as library of diverse sulfated molecules and in silico docking and scoring protocol have not been firmly established. Databases, such as ZINC, contain too few sulfate-containing nonsaccharide molecules, which severely limits the identification of new hits. Lack of a generally applicable protocol for scaffold hopping limits the development of sulfated small molecules as synthetic mimetics of the highly sulfated glycosaminoglycans. We explored a sequential ligand-based (LBVS) and structure-based virtual screening (SBVS) approach starting from our initial discovery of monosulfated benzofurans to discover alternative scaffolds as allosteric modulators of thrombin, a key coagulation enzyme. Screening the ZINC database containing nearly 1 million nonsulfated small molecules using a pharmacophore developed from the parent sulfated benzofurans followed by a genetic algorithm-based dual-filter docking and scoring screening identified a group of 10 promising hits, of which three top-scoring hits were synthesized. Each was found to selectively inhibit human alpha-thrombin suggesting the possibility of this approach for scaffold hopping. Michaelis–Menten kinetics showed allosteric inhibition mechanism for the best molecule and human plasma studies confirmed good anticoagulation potential as expected. Our simple sequential LBVS and SBVS approach is likely to be useful as a general strategy for identification of sulfated small molecules hits as modulators of glycosaminoglycan–protein interactions.     

Virtual screening using ligand-based pharmacophores for inhibitors of human tyrosyl-DNA phospodiesterase (hTdp1)

Human tyrosyl-DNA phosphodiesterase (hTdp1) inhibitors have become a major area of drug research and structure-based design since they have been shown to work synergistically with camptothecin (CPT) and selectively in cancer cells. The pharmacophore features of 14 hTdp1 inhibitors were used as a filter to screen the ChemNavigator iResearch Library of about 27 million purchasable samples. Docking of the inhibitors and hits obtained from virtual screening was performed into a structural model of hTdp1 based on a high resolution X-ray crystal structure of human Tdp1 in complex with vanadate, DNA and a human topoisomerase I (Top1)-derived peptide (PDB code: 1NOP). We present and discuss in some detail 46 compounds matching the three-dimensional arrangement of the pharmacophoric features. The presented novel chemotypes may provide new scaffolds for developing inhibitors of Tdp1.     

Identification of new γ-hydroxybutenolides that preferentially inhibit the activity of mPGES-1

Microsomal prostaglandin E(2) synthase-1 (mPGES-1) has been recognized as novel, promising drug target for anti-inflammatory and anticancer drugs. mPGES-1 catalyzes the synthesis of the inducible prostaglandin E(2) in response to pro-inflammatory stimuli, rendering this enzyme extremely interesting in drug discovery process owing to the drastic reduction of the severe side effects typical for traditional non-steroidal anti-inflammatory drugs. In the course of our investigations focused on this topic, we identified two interesting molecules bearing the γ-hydroxybutenolide scaffold which potently inhibit the activity of mPGES-1. Notably, the lead compound 2c that inhibited mPGES-1 with IC(50)=0.9μM, did not affect other related enzymes within the arachidonic acid cascade.     

Docking, virtual high throughput screening and in silico fragment-based drug design

The drug discovery process has been profoundly changed recently by the adoption of computational methods helping the design of new drug candidates more rapidly and at lower costs. In silico drug design consists of a collection of tools helping to make rational decisions at the different steps of the drug discovery process, such as the identification of a biomolecular target of therapeutical interest, the selection or the design of new lead compounds and their modification to obtain better affinities, as well as pharmacokinetic and pharmacodynamic properties. Among the different tools available, a particular emphasis is placed in this review on molecular docking, virtual high-throughput screening and fragment-based ligand design.     

Novel DNA Topoisomerase IIα Inhibitors from Combined Ligand- and Structure-Based Virtual Screening

DNA topoisomerases are enzymes responsible for the relaxation of DNA torsional strain, as well as for the untangling of DNA duplexes after replication, and are important cancer drug targets. One class of topoisomerase inhibitors, “poisons”, binds to the transient enzyme-DNA complex which occurs during the mechanism of action, and inhibits the religation of DNA. This ultimately leads to the accumulation of DNA double strand breaks and cell death. Different types of topoisomerases occur in human cells and several poisons of topoisomerase I and II are widely used clinically. However, their use is compromised by a variety of side effects. Recent studies confirm that the inhibition of the α-isoform of topoisomerase II is responsible for the cytotoxic effect, whereas the inhibition of the β-isoform leads to development of adverse drug reactions. Thus, the discovery of agents selective for topoisomerase IIα is an important strategy for the development of topoisomerase II poisons with improved clinical profiles. Here, we present a computer-aided drug design study leading to the identification of structurally novel topoisomerase IIα poisons. The study combines ligand- and structure-based drug design methods including pharmacophore models, homology modelling, docking, and virtual screening of the National Cancer Institute compound database. From the 8 compounds identified from the computational work, 6 were tested for their capacity to poison topoisomerase II in vitro: 4 showed selective inhibitory activity for the α- over the β-isoform and 3 of these exhibited cytotoxic activity. Thus, our study confirms the applicability of computer-aided methods for the discovery of novel topoisomerase II poisons, and presents compounds which could be investigated further as selective topoisomerase IIα inhibitors.