Reproducing crystal binding modes of ligand functional groups using Site-Identification by Ligand Competitive Saturation (SILCS) simulations

The applicability of a computational method, Site Identification by Ligand Competitive Saturation (SILCS), to identify regions on a protein surface with which different types of functional groups on low-molecular weight inhibitors interact is demonstrated. The method involves molecular dynamics (MD) simulations of a protein in an aqueous solution of chemically diverse small molecules from which probability distributions of fragments types, termed FragMaps, are obtained. In the present application, SILCS simulations are performed with an aqueous solution of 1 M benzene and propane to map the affinity pattern of the protein for aromatic and aliphatic functional groups. In addition, water hydrogen and oxygen atoms serve as probes for hydrogen-bond donor and acceptor affinity, respectively. The method is tested using a set of 7 proteins for which crystal structures of complexes with several high affinity inhibitors are known. Good agreement is obtained between FragMaps and the positions of chemically similar functional groups in inhibitors as observed in the X-ray crystallographic structures. Quantitative capabilities of the SILCS approach are demonstrated by converting FragMaps to free energies, termed Grid Free Energies (GFE), and showing correlation between the GFE values and experimental binding affinities. For proteins for which ligand decoy sets are available, GFE values are shown to typically score the crystal conformation and conformations similar to it more favorable than decoys. Additionally, SILCS is tested for its ability to capture the subtle differences in ligand affinity across homologous proteins, information which may be of utility toward specificity-guided drug design. Taken together, our results show that SILCS can recapitulate the known location of functional groups of bound inhibitors for a number of proteins, suggesting that the method may be of utility for rational drug design.     

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.     

人类腺苷受体A3亚型拮抗剂的构效关系分析

构建人类腺苷受体A3亚型药效团模型和三维蛋白结构模型用于作用模式研究.以18个来源于文献具有腺苷受体A3亚型拮抗活性的化合物作为训练集,使用HypoGen方法构建药效团模型.通过同源模建和分子动力学模拟构建了人类腺苷受体A3亚型的三维蛋白模型,并利用PROCHECK方法评估该模型的合理性,对所得的结构使用分子对接程序进行作用模式分析,药效团模型和同源模建结果相互匹配较好.使用新药效团模型对MDL药物数据库(MDDR)中包含的约120000个化合物进行虚拟筛选,得到了8个候选化合物,用于进一步的生物学评价和活性测定.本工作对于人类腺苷受体A3亚型拮抗剂的设计和抗哮喘药物的研发具有一定的理论指导和应用价值.     

Development of biologically active compounds by combining 3D QSAR and structure-based design methods

One of the major challenges in computational approaches to drug design is the accurate prediction of the binding affinity of novel biomolecules. In the present study an automated procedure which combines docking and 3D-QSAR methods was applied to several drug targets. The developed receptor-based 3D-QSAR methodology was tested on several sets of ligands for which the three-dimensional structure of the target protein has been solved – namely estrogen receptor, acetylcholine esterase and protein-tyrosine-phosphatase 1B. The molecular alignments of the studied ligands were determined using the docking program AutoDock and were compared with the X-ray structures of the corresponding protein-ligand complexes. The automatically generated protein-based ligand alignment obtained was subsequently taken as basis for a comparative field analysis applying the GRID/GOLPE approach. Using GRID interaction fields and applying variable selection procedures, highly predictive models were obtained. It is expected that concepts from receptor-based 3D QSAR will be valuable tools for the analysis of high-throughput screening as well as virtual screening data     

Combined receptor and ligand-based approach to the universal pharmacophore model development for studies of drug blockade to the hERG1 pore domain

Long QT syndrome, LQTS, results in serious cardiovascular disorders, such as tachyarrhythmia and sudden cardiac death. A promiscuous binding of different drugs to the intracavitary binding site in the pore domain (PD) of human ether-a-go-go related gene (hERG) channels leads to a similar dysfunction, known as a drug-induced LQTS. Therefore, an assessment of the blocking ability for potent drugs is of great pragmatic value for molecular pharmacology and medicinal chemistry of hERGs. Thus, we attempted to create an in silico model aimed at blinded drug screening for their blocking ability to the hERG1 PD. Two distinct approaches to the drug blockage, ligand-based QSAR and receptor-based molecular docking methods, are combined for development of a universal pharmacophore model, which provides rapid assessment of drug blocking ability to the hERG1 channel. The best 3D-QSAR model (AAADR.7) from PHASE modeling was selected from a pool consisting of 44 initial candidates. The constructed model using 31 hERG blockers was validated with 9 test set compounds. The resulting model correctly predicted the pIC(50) values of test set compounds as true unknowns. To further evaluate the pharmacophore model, 14 hERG blockers with diverse hERG blocking potencies were selected from literature and they were used as additional external blind test sets. The resulting average deviation between in vitro and predicted pIC(50) values of external test set blockers is found as 0.29 suggesting that the model is able to accuretely predict the pIC(50) values as true unknowns. These pharmacophore models were merged with a previously developed atomistic receptor model for the hERG1 PD and exhibited a high consistency between ligand-based and receptor-based models. Therefore, the developed 3D-QSAR model provides a predictive tool for profiling candidate compounds before their synthesis. This model also indicated the key functional groups determining a high-affinity blockade of the hERG1 channel. To cross-validate consistency between the constructed hERG1 pore domain and the pharmacophore models, we performed docking studies using the homology model of hERG1. To understand how polar or nonpolar moieties of inhibitors stimulate channel inhibition, critical amino acid replacement (i.e., T623, S624, S649, Y652 and F656) at the hERG cavity was examined by in silico mutagenesis. The average docking score differences between wild type and mutated hERG channels was found to have the following order: F656A > Y652A > S624A > T623A > S649A. These results are in agreement with experimental data.     

Protein flexibility and species specificity in structure-based drug discovery: dihydrofolate reductase as a test system

In structure-based drug discovery, researchers would like to identify all possible scaffolds for a given target. However, techniques that push the boundaries of chemical space could lead to many false positives or inhibitors that lack specificity for the target. Is it possible to broadly identify the appropriate chemical space for the inhibitors and yet maintain target specificity? To address this question, we have turned to dihydrofolate reductase (DHFR), a well-studied metabolic enzyme of pharmacological relevance. We have extended our multiple protein structure (MPS) method for receptor-based pharmacophore models to use multiple X-ray crystallographic structures. Models were created for DHFR from human and Pneumocystis carinii. These models incorporate a fair degree of protein flexibility and are highly selective for known DHFR inhibitors over drug-like non-inhibitors. Despite sharing a highly conserved active site, the pharmacophore models reflect subtle differences between the human and P. carinii forms, which identify species-specific, high-affinity inhibitors. We also use structures of DHFR from Candida albicans as a counter example. The available crystal structures show little flexibility, and the resulting models give poorer performance in identifying species-specific inhibitors. Therapeutic success for this system may depend on achieving species specificity between the related human host and these key fungal targets. The MPS technique is a promising advance for structure-based drug discovery for DHFR and other proteins of biomedical interest.     

Incorporating protein flexibility in structure-based drug discovery: using HIV-1 protease as a test case

We have developed a receptor-based pharmacophore method which utilizes a collection of protein structures to account for inherent protein flexibility in structure-based drug design. Several procedures were systematically evaluated to derive the most general protocol for using multiple protein structures. Most notably, incorporating more protein flexibility improved the performance of the method. The pharmacophore models successfully discriminate known inhibitors from drug-like non-inhibitors. Furthermore, the models correctly identify the bound conformations of some ligands. We used unliganded HIV-1 protease to develop and validate this method. Drug design is always initiated with a protein-ligand structure, and such success with unbound protein structures is remarkable - particularly in the case of HIV-1 protease, which has a large conformational change upon binding. This technique holds the promise of successful computer-based drug design before bound crystal structures are even discovered, which can mean a jump-start of 1-3 years in tackling some medically relevant systems with computational methods.     

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.     

Comparison of ligand-based and receptor-based virtual screening of HIV entry inhibitors for the CXCR4 and CCR5 receptors using 3D ligand shape matching and ligand-receptor docking

HIV infection is initiated by fusion of the virus with the target cell through binding of the viral gp120 protein with the CD4 cell surface receptor protein and the CXCR4 or CCR5 co-receptors. There is currently considerable interest in developing novel ligands that can modulate the conformations of these co-receptors and, hence, ultimately block virus-cell fusion. This article describes a detailed comparison of the performance of receptor-based and ligand-based virtual screening approaches to find CXCR4 and CCR5 antagonists that could potentially serve as HIV entry inhibitors. Because no crystal structures for these proteins are available, homology models of CXCR4 and CCR5 have been built, using bovine rhodopsin as the template. For ligand-based virtual screening, several shape-based and property-based molecular comparison approaches have been compared, using high-affinity ligands as query molecules. These methods were compared by virtually screening a library assembled by us, consisting of 602 known CXCR4 and CCR5 inhibitors and some 4700 similar presumed inactive molecules. For each receptor, the library was queried using known binders, and the enrichment factors and diversity of the resulting virtual hit lists were analyzed. Overall, ligand-based shape-matching searches yielded higher enrichments than receptor-based docking, especially for CXCR4. The results obtained for CCR5 suggest the possibility that different active scaffolds bind in different ways within the CCR5 pocket.     

Hot-spots-guided receptor-based pharmacophores (HS-Pharm): a knowledge-based approach to identify ligand-anchoring atoms in protein cavities and prioritize structure-based pharmacophores

The design of biologically active compounds from ligand-free protein structures using a structure-based approach is still a major challenge. In this paper, we present a fast knowledge-based approach (HS-Pharm) that allows the prioritization of cavity atoms that should be targeted for ligand binding, by training machine learning algorithms with atom-based fingerprints of known ligand-binding pockets. The knowledge of hot spots for ligand binding is here used for focusing structure-based pharmacophore models. Three targets of pharmacological interest (neuraminidase, beta2 adrenergic receptor, and cyclooxygenase-2) were used to test the evaluated methodology, and the derived structure-based pharmacophores were used in retrospective virtual screening studies. The current study shows that structure-based pharmacophore screening is a powerful technique for the fast identification of potential hits in a chemical library, and that it is a valid alternative to virtual screening by molecular docking.     

Pharmacophore modeling and parallel screening for PPAR ligands

We describe the generation and validation of pharmacophore models for PPARs, as well as a large scale validation of the parallel screening approach by screening PPAR ligands against a large database of structure-based models. A large test set of 357 PPAR ligands was screened against 48 PPAR models to determine the best models for agonists of PPAR-alpha, PPAR-delta, and PPAR-gamma. Afterwards, a parallel screen was performed using the 357 PPAR ligands and 47 structure-based models for PPARs, which were integrated into a 1537 models comprising in-house pharmacophore database, to assess the enrichment of PPAR ligands within the PPAR hypotheses. For these purposes, we categorized the 1537 database models into 181 protein targets and developed a score that ranks the retrieved targets for each ligand. Thus, we tried to find out if the concept of parallel screening is able to predict the correct pharmacological target for a set of compounds. The PPAR target was ranked first more often than any other target. This confirms the ability of parallel screening to forecast the pharmacological active target for a set of compounds.     

A structure-based design approach for the identification of novel inhibitors: application to an alanine racemase

We report a new structure-based strategy for the identification of novel inhibitors. This approach has been applied to Bacillus stearothermophilus alanine racemase (AlaR), an enzyme implicated in the biosynthesis of the bacterial cell wall. The enzyme catalyzes the racemization of l- and d-alanine using pyridoxal 5-phosphate (PLP) as a cofactor. The restriction of AlaR to bacteria and some fungi and the absolute requirement for d-alanine in peptidoglycan biosynthesis make alanine racemase a suitable target for drug design. Unfortunately, known inhibitors of alanine racemase are not specific and inhibit the activity of other PLP-dependent enzymes, leading to neurological and other side effects.This article describes the development of a receptor-based pharmacophore model for AlaR, taking into account receptor flexibility (i.e. a `dynamic' pharmacophore model). In order to accomplish this, molecular dynamics (MD) simulations were performed on the full AlaR dimer from Bacillus stearothermophilus (PDB entry, 1sft) with a d-alanine molecule in one active site and the non-covalent inhibitor, propionate, in the second active site of this homodimer. The basic strategy followed in this study was to utilize conformations of the protein obtained during MD simulations to generate a dynamic pharmacophore model using the property mapping capability of the LigBuilder program. Compounds from the Available Chemicals Directory that fit the pharmacophore model were identified and have been submitted for experimental testing.The approach described here can be used as a valuable tool for the design of novel inhibitors of other biomolecular targets.     

Incorporating dynamics in E. coli dihydrofolate reductase enhances structure-based drug discovery

Escherichia coli dihydrofolate reductase (DHFR) is a long-standing target for enzyme studies. The influence of protein motion on its catalytic cycle is significant, and the conformation of the M20 loop is of particular interest. We present receptor-based pharmacophore models-an equivalent of solvent-mapping of binding hotspots-based on ensembles of protein conformations from molecular dynamics simulations of DHFR.NADPH in both the closed and open conformation of the M20 loop. The optimal models identify DHFR inhibitors over druglike non-inhibitors; furthermore, high-affinity inhibitors of E. coli DHFR are preferentially identified over general DHFR inhibitors. As expected, models resulting from simulations with DHFR in the productive conformation with a closed M20 loop have better performance than those from the open-loop simulations. Model performance improves with increased dynamic sampling, indicating that including a greater degree of protein flexibility can enhance the quest for potent inhibitors. This was compared to the limited conformational sampling seen in crystal structures, which were suboptimal for this application.     

Homology model-based virtual screening for GPCR ligands using docking and target-biased scoring

The current study investigates the combination of two recently reported techniques for the improvement of homology model-based virtual screening for G-protein coupled receptor (GPCR) ligands. First, ligand-supported homology modeling was used to generate receptor models that were in agreement with mutagenesis data and structure-activity relationship information of the ligands. Second, interaction patterns from known ligands to the receptor were applied for scoring and rank ordering compounds from a virtual library using ligand-receptor interaction fingerprint-based similarity (IFS). Our approach was evaluated in retrospective virtual screening experiments for antagonists of the metabotropic glutamate receptor (mGluR) subtype 5. The results of our approach were compared to the results obtained by conventional scoring functions (Dock-Score, PMF-Score, Gold-Score, ChemScore, and FlexX-Score). The IFS lead to significantly higher enrichment rates, relative to the competing scoring functions. Though using a target-biased scoring approach, the results were not biased toward the chemical classes of the reference structures. Our results indicate that the presented approach has the potential to serve as a general setup for successful structure-based GPCR virtual screening.     

Using a pharmacophore representation concept to elucidate molecular similarity of dopamine antagonists

The pharmacophoric concept plays an important role in ligand-based drug design methods to describe the similarity and diversity of molecules, and could also be exploited as a molecular representation scheme. A three-point pharmacophore method was used as a molecular representation perception. This procedure was implemented for dopamine antagonists of the D(2) receptor subtype. The molecular structures of the antagonists included in this analysis were categorized into two structurally distinct classes. Using structural superposition with internal energy minimization, two pharmacophore models were deduced. Based on these two models other D(2) antagonists that fulfil them were derived and studied. This procedure aided the identification of the common 3D patterns present in diverse molecules that act at the same biological target and the extraction of a common molecular framework for the two structural classes. The pharmacophoric information was found to be suitable for guiding superposition of structurally diverse molecules, using a more biologically meaningful selection of the targeting points.     

Targets looking for drugs: a multistep computational protocol for the development of structure-based pharmacophores and their applications for hit discovery

Pharmacophoresthree-dimensional (3D) arrangements of essential features enabling a molecule to exert a particular biological effectconstitute a very useful tool in drug design both in hit discovery and hit-to-lead optimization process. Two basic approaches for pharmacophoric model generation can be used by chemists, depending on the availability or not of the target 3D structure. In view of the rapidly growing number of protein structures that are now available, receptor-based pharmacophore generation methods are becoming more and more used. Since most of them require the knowledge of the 3D structure of the ligand-target complex, they cannot be applied when no compounds targeting the binding site of interest are known. Here, a GRID-based procedure for the generation of receptor-based pharmacophores starting from the knowledge of the sole protein structure is described and successfully applied to address three different tasks in the field of medicinal chemistry.     

Comparing pharmacophore models derived from crystallography and NMR ensembles

NMR and X-ray crystallography are the two most widely used methods for determining protein structures. Our previous study examining NMR versus X-Ray sources of protein conformations showed improved performance with NMR structures when used in our Multiple Protein Structures (MPS) method for receptor-based pharmacophores (Damm, Carlson, J Am Chem Soc 129:8225–8235, 2007). However, that work was based on a single test case, HIV-1 protease, because of the rich data available for that system. New data for more systems are available now, which calls for further examination of the effect of different sources of protein conformations. The MPS technique was applied to Growth factor receptor bound protein 2 (Grb2), Src SH2 homology domain (Src-SH2), FK506-binding protein 1A (FKBP12), and Peroxisome proliferator-activated receptor-γ (PPAR-γ). Pharmacophore models from both crystal and NMR ensembles were able to discriminate between high-affinity, low-affinity, and decoy molecules. As we found in our original study, NMR models showed optimal performance when all elements were used. The crystal models had more pharmacophore elements compared to their NMR counterparts. The crystal-based models exhibited optimum performance only when pharmacophore elements were dropped. This supports our assertion that the higher flexibility in NMR ensembles helps focus the models on the most essential interactions with the protein. Our studies suggest that the “extra” pharmacophore elements seen at the periphery in X-ray models arise as a result of decreased protein flexibility and make very little contribution to model performance.     

Approximating protein flexibility through dynamic pharmacophore models: application to fatty acid amide hydrolase (FAAH)

A structure-based drug discovery method is described that incorporates target flexibility through the use of an ensemble of protein conformations. The approach was applied to fatty acid amide hydrolase (FAAH), a key deactivating enzyme in the endocannabinoid system. The resultant dynamic pharmacophore models are rapidly able to identify known FAAH inhibitors over drug-like decoys. Different sources of FAAH conformational ensembles were explored, with both snapshots from molecular dynamics simulations and a group of X-ray structures performing well. Results were compared to those from docking and pharmacophore models generated from a single X-ray structure. Increasing conformational sampling consistently improved the pharmacophore models, emphasizing the importance of incorporating target flexibility in structure-based drug design.