Predicting hiCE inhibitors based upon pharmacophore models derived from the receptor and its ligands

Human intestinal carboxyl esterase (hiCE) is a drug target for ameliorating irinotecan-induced diarrhea. By reducing irinotecan-induced diarrhea, hiCE inhibitors can improve the anti-cancer efficacy of irinotecan. To find effective hiCE inhibitors, a new virtual screening protocol that combines pharmacophore models derived from the hiCE structure and its ligands has been proposed. The hiCE structure has been constructed through homology techniques using hCES1’s crystal structure. The hiCE structure was optimized via molecular dynamics simulations with the most known active hiCE inhibitors docked into the structure. An optimized pharmacophore, derived from the receptor, was then generated. A ligand-based pharmacophore was also generated from a larger set of known hiCE inhibitors. The final hiCE inhibitor predictions were based upon the virtual screening hits from both ligand-based and receptor-based pharmacophore models. The hit rates from the ligand-based and receptor-based pharmacophore models are 88% and 86%, respectively. The final hit rate is 94%. The two models are highly consistent with one another (85%). This proves that both models are reliable.     

Development of predictive pharmacophore model for in silico screening,and 3D QSAR CoMFA and CoMSIA studies for lead optimization,for designing of potent tumor necrosis factor alpha converting enzyme inhibitors

A chemical feature-based pharmacophore model was developed for Tumor Necrosis Factor-α converting enzyme (TACE) inhibitors. A five point pharmacophore model having two hydrogen bond acceptors (A), one hydrogen bond donor (D) and two aromatic rings (R) with discrete geometries as pharmacophoric features was developed. The pharmacophore model so generated was then utilized for in silico screening of a database. The pharmacophore model so developed was validated by using four compounds having proven TACE inhibitory activity which were grafted into the database. These compounds mapped well onto the five listed pharmacophoric features. This validated pharmacophore model was also used for alignment of molecules in CoMFA and CoMSIA analysis. The contour maps of the CoMFA/CoMSIA models were utilized to provide structural insight for activity improvement of potential novel TACE inhibitors. The pharmacophore model so developed could be used for in silico screening of any commercial/in house database for identification of TACE inhibiting lead compounds, and the leads so identified could be optimized using the developed CoMSIA model. The present work highlights the tremendous potential of the two mutually complementary ligand-based drug designing techniques (i.e. pharmacophore mapping and 3D-QSAR analysis) using TACE inhibitors as prototype biologically active molecules.     

PhAST: pharmacophore alignment search tool

We present a ligand-based virtual screening technique (PhAST) for rapid hit and lead structure searching in large compound databases. Molecules are represented as strings encoding the distribution of pharmacophoric features on the molecular graph. In contrast to other text-based methods using SMILES strings, we introduce a new form of text representation that describes the pharmacophore of molecules. This string representation opens the opportunity for revealing functional similarity between molecules by sequence alignment techniques in analogy to homology searching in protein or nucleic acid sequence databases. We favorably compared PhAST with other current ligand-based virtual screening methods in a retrospective analysis using the BEDROC metric. In a prospective application, PhAST identified two novel inhibitors of 5-lipoxygenase product formation with minimal experimental effort. This outcome demonstrates the applicability of PhAST to drug discovery projects and provides an innovative concept of sequence-based compound screening with substantial scaffold hopping potential.     

Identification of potential Tpx inhibitors against pathogen-host interactions

Yersinia organisms cause many infectious diseases by invading human cells and delivering their virulence factors via the type three secretion system (T3SS). One alternative strategy in the fight against these pathogenic organisms is to interfere with their T3SS. Previous studies demonstrated that thiol peroxidase, Tpx is functional in the assembly of T3SS and its inhibition by salicylidene acylhydrazides prevents the secretion of pathogenic effectors. In this study, the aim was to identify potential inhibitors of Tpx using an integrated approach starting with high throughput virtual screening and ending with molecular dynamics simulations of selected ligands. Virtual screening of ZINC database of 500,000 compounds via ligand-based and structure-based pharmacophore models retrieved 10,000 hits. The structure-based pharmacophore model was validated using high-throughput virtual screening (HTVS). After multistep docking (SP and XP), common scaffolds were used to find common substructures and the ligand binding poses were optimized using induced fit docking. The stability of the protein–ligand complex was examined with molecular dynamics simulations and the binding free energy of the complex was calculated. As a final outcome eight compounds with different chemotypes were proposed as potential inhibitors for Tpx. The eight ligands identified by a detailed virtual screening protocol can serve as leads in future drug design efforts against the destructive actions of pathogenic bacteria.     

SHAFTS: a hybrid approach for 3D molecular similarity calculation. 1. Method and assessment of virtual screening

We developed a novel approach called SHAFTS (SHApe-FeaTure Similarity) for 3D molecular similarity calculation and ligand-based virtual screening. SHAFTS adopts a hybrid similarity metric combined with molecular shape and colored (labeled) chemistry groups annotated by pharmacophore features for 3D similarity calculation and ranking, which is designed to integrate the strength of pharmacophore matching and volumetric overlay approaches. A feature triplet hashing method is used for fast molecular alignment poses enumeration, and the optimal superposition between the target and the query molecules can be prioritized by calculating corresponding "hybrid similarities". SHAFTS is suitable for large-scale virtual screening with single or multiple bioactive compounds as the query "templates" regardless of whether corresponding experimentally determined conformations are available. Two public test sets (DUD and Jain's sets) including active and decoy molecules from a panel of useful drug targets were adopted to evaluate the virtual screening performance. SHAFTS outperformed several other widely used virtual screening methods in terms of enrichment of known active compounds as well as novel chemotypes, thereby indicating its robustness in hit compounds identification and potential of scaffold hopping in virtual screening.     

On the value of homology models for virtual screening: discovering hCXCR3 antagonists by pharmacophore-based and structure-based approaches

Human chemokine receptor CXCR3 (hCXCR3) antagonists have potential therapeutic applications as antivirus, antitumor, and anti-inflammatory agents. A novel virtual screening protocol, which combines pharmacophore-based and structure-based approaches, was proposed. A three-dimensional QSAR pharmacophore model and a structure-based docking model were built to virtually screen for hCXCR3 antagonists. The hCXCR3 antagonist binding site was constructed by homology modeling and molecular dynamics (MD) simulation. By combining the structure-based and ligand-based screenings results, 95% of the compounds satisfied either pharmacophore or docking score criteria and would be chosen as hits if the union of the two searches was taken. The false negative rates were 15% for the pharmacophore model, 14% for the homology model, and 5% for the combined model. Therefore, the consistency of the pharmacophore model and the structural binding model is 219/273 = 80%. The hit rate for the virtual screening protocol is 273/286 = 95%. This work demonstrated that the quality of both the pharmacophore model and homology model can be measured by the consistency of the two models, and the false negatives in virtual screening can be reduced by combining two virtual screening approaches.     

New serotonin 5-HT(6) ligands from common feature pharmacophore hypotheses

Serotonin 5-HT6 receptor antagonists are thought to play an important role in the treatment of psychiatry, Alzheimer's disease, and probably obesity. To find novel and potent 5-HT6 antagonists and to provide a new idea for drug design, we used a ligand-based pharmacophore to perform the virtual screening of a commercially available database. A three-dimensional common feature pharmacophore model was developed by using the HipHop program provided in Catalyst software and was used as a query for screening the database. A recursive partitioning (RP) model which can separate active and inactive compounds was used as a filtering system. Finally a sequential virtual screening procedure (SQSP) was conducted, wherein both the common feature pharmacophore and the RP model were used in succession to improve the results. Some of the hits were selected based on druglikeness, ADME properties, structural diversity, and synthetic accessibility for real biological evaluation. The best hit compound showed a significant IC50 value of 9.6 nM and can be used as a lead for further drug development.     

Pharmacophore identification, in silico screening, and virtual library design for inhibitors of the human factor Xa

Factor Xa inhibitors are innovative anticoagulant agents that provide a better safety/efficacy profile compared to other anticoagulative drugs. A chemical feature-based modeling approach was applied to identify crucial pharmacophore patterns from 3D crystal structures of inhibitors bound to human factor Xa (Pdb entries 1fjs, 1kns, 1eqz) using the software LIGANDSCOUT and CATALYST. The complex structures were selected regarding the criteria of high inhibitory potency (i.e. all ligands show K(i) values against factor Xa in the subnanomolar range) and good resolution (i.e. at least 2.2 A) in order to generate selective and high quality pharmacophore models. The resulting chemical-feature based hypotheses were used for virtual screening of commercial molecular databases such as the WDI database. Furthermore, a ligand-based molecular modeling approach was performed to obtain common-feature hypotheses that represent the relevant chemical interactions between 10 bioactive factor Xa inhibitors and the protein, respectively. In a next step a virtual combinatorial library was designed in order to generate new compounds with similar chemical and spatial properties as known inhibitors. The software tool ILIB DIVERSE was used for this procedure in order to provide new scaffolds of this group of anticoagulants. Finally we present the combination of these two techniques, hence virtual screening was performed with selective pharmacophore models in a focused virtual combinatorial database. De novo derived molecular scaffolds that were able to adequately satisfy the pharmacophore criteria are revealed and are promising templates for candidates for further development.     

2D conformationally sampled pharmacophore: a ligand-based pharmacophore to differentiate delta opioid agonists from antagonists

Pharmacophores are widely used for rational drug design and include those based on receptor binding sites or on known ligands. To date, ligand-based pharmacophores have typically used one or a small number of conformers of known receptor ligands. However, this method does not take into account the inherent dynamic nature of molecules, which sample a wide range of conformations, any of which could be the bound form. In the present study, molecular dynamics (MD) simulations were used as a means to sample the conformational space of ligands to include all accessible conformers at room temperature in the development of a pharmacophore. On the basis of these conformers, probability distributions of selected distances and angles in a series of delta specific opioid ligands were obtained and correlated with agonist versus antagonist activities. Individually, the distributions did not allow for unique agonist and antagonist pharmacophores to be identified. However, by extending the conformational analysis to two dimensions, a 2D conformationally sampled pharmacophore (CSP) for distinguishing delta receptor agonists and antagonists was developed. Application of this model to the compound DPI2505 suggests that it may have agonist activity. It is anticipated that the CSP method, which does not require alignment of compounds during pharmacophore development, will be a useful tool for obtaining structure-function relationships of ligands particularly in systems where the receptor 3D structure is not known.     

The discovery of Kv1.5 blockers as a case study for the application of virtual screening approaches

Different virtual screening techniques are available as alternatives to high throughput screening. These different techniques have been rarely used together on the same target. We had the opportunity to do so in order to discover novel blockers of the voltage-dependent potassium channel Kv1.5, a potential target for the treatment of atrial fibrillation. Our corporate database was searched, using a protein-based pharmacophore, derived from a homology model, as query. As a result, 244 molecules were screened in vitro, 19 of them (7.8%) were found to be active. Five of them, belonging to five different chemical classes, exhibited IC50 values under 10 microM. The performance of this structure-based virtual screening protocol has been compared with those of similarity and ligand-based pharmacophore searches. The analysis of the results supports the conventional wisdom of using as many virtual screening techniques as possible in order to maximize the chance of finding as many chemotypes as possible.     

Effective virtual screening protocol for CYP2C9 ligands using a screening site constructed from flurbiprofen and S-warfarin pockets

An effective virtual screening protocol was developed against an extended active site of CYP2C9, which was derived from X-ray structures complexed with flubiprofen and S-warfarin. Virtual screening has been effectively supported by our structure-based pharmacophore model. Importance of hot residues identified by mutation data and structural analysis was first estimated in an enrichment study. Key role of Arg108 and Phe114 in ligand binding was also underlined. Our screening protocol successfully identified 76% of known CYP2C9 ligands in the top 1% of the ranked database resulting 76-fold enrichment relative to random situation. Relevance of the protocol was further confirmed in selectivity studies, when 89% of CYP2C9 ligands were retrieved from a mixture of CYP2C9 and CYP2C8 ligands, while only 22% of CYP2C8 ligands were found applying the structure-based pharmacophore constraints. Moderate discrimination of CYP2C9 ligands from CYP2C18 and CYP2C19 ligands could also be achieved extending the application domain of our virtual screening protocol for the entire CYP2C family. Our findings further demonstrate the existence of an active site comprising of at least two binding pockets and strengthens the need of involvement of protein flexibility in virtual screening.     

Toward the discovery of functional transthyretin amyloid inhibitors: application of virtual screening methods

Inhibition of amyloid fibril formation by stabilization of the native form of the protein transthyretin (TTR) is a viable approach for the treatment of familial amyloid polyneuropathy that has been gaining momentum in the field of amyloid research. The TTR stabilizer molecules discovered to date have shown efficacy at inhibiting fibrilization in vitro but display impairing issues of solubility, affinity for TTR in the blood plasma and/or adverse effects. In this study we present a benchmark of four protein- and ligand-based virtual screening (VS) methods for identifying novel TTR stabilizers: (i) two-dimensional (2D) similarity searches with chemical hashed, pharmacophore, and UNITY fingerprints, (ii) 3D searches based on shape, chemical, and electrostatic similarity, (iii) LigMatch, a new ligand-based method which uses multiple templates and combines 3D geometric hashing with a 2D preselection process, and (iv) molecular docking to consensus X-ray crystal structures of TTR. We illustrate the potential of the best-performing VS protocols to retrieve promising new leads by ranking a tailored library of 2.3 million commercially available compounds. Our predictions show that the top-scoring molecules possess distinctive features from the known TTR binders, holding better solubility, fraction of halogen atoms, and binding affinity profiles. To the best of our knowledge, this is the first attempt to rationalize the utilization of a large battery of in silico screening techniques toward the identification of a new generation of TTR amyloid inhibitors.     

A pseudo-ligand approach to virtual screening

A virtual screening method is presented that is grounded on a receptor-derived pharmacophore model termed "virtual ligand" or "pseudo-ligand". The model represents an idealized constellation of potential ligand sites that interact with residues of the binding pocket. For rapid virtual screening of compound libraries the potential pharmacophore points of the virtual ligand are encoded as an alignment-free correlation vector, avoiding spatial alignment of pharmacophore features between the pharmacophore query (i.e., the virtual ligand) and the candidate molecule. The method was successfully applied to retrieving factor Xa inhibitors from a Ugi three-component combinatorial library, and yielded high enrichment of actives in a retrospective search for cyclooxygenase-2 (COX-2) inhibitors. The approach provides a concept for "de-orphanizing" potential drug targets and identifying ligands for hitherto unexplored or allosteric binding pockets.     

Identification of potential HIV-1 integrase strand transfer inhibitors: In silico virtual screening and QM/MM docking studies

HIV-1 integrase (IN) is a retroviral enzyme that catalyses integration of the reverse-transcribed viral DNA into the host genome, which is necessary for efficient viral replication. In this study, we have performed an in silico virtual screening for the identification of potential HIV-1 IN strand transfer (ST) inhibitors. Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to 3-Hydroxypyrimidine-2,4-diones. Based on the ligand-based pharmacophore model, we obtained a five-point pharmacophore with two hydrogen bond acceptors (A), one hydrogen bond donor (D), one hydrophobic group (H) and one aromatic ring (R) as pharmacophoric features. The pharmacophore hypothesis AADHR was used as a 3D query in a sequential virtual screening study to filter small molecule databases Maybridge, ChemBridge and Asinex. Hits matching with pharmacophore hypothesis AADHR were retrieved and passed progressively through Lipinski’s rule of five filtering, molecular docking and hierarchical clustering. The five compounds with best hits with novel and diverse chemotypes were subjected to QM/MM docking, which showed improved docking accuracy. We further performed molecular dynamics simulation and found three compounds that form stable interactions with key residues. These compounds could be used as a leads for further drug development and rational design of HIV-1 IN inhibitors.