1
|
Szwabowski GL, Daigle BJ, Baker DL, Parrill AL. Structure-based pharmacophore modeling 2. Developing a novel framework for structure-based pharmacophore model generation and selection. J Mol Graph Model 2023; 122:108488. [PMID: 37121167 DOI: 10.1016/j.jmgm.2023.108488] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2022] [Accepted: 04/06/2023] [Indexed: 05/02/2023]
Abstract
Pharmacophore models are three-dimensional arrangements of molecular features required for biological activity that are used in ligand identification efforts for many biological targets, including G protein-coupled receptors (GPCR). Though GPCR are integral membrane proteins of considerable interest as targets for drug development, many of these receptors lack known ligands or experimentally determined structures necessary for ligand- or structure-based pharmacophore model generation, respectively. Thus, we here present a structure-based pharmacophore modeling approach that uses fragments placed with Multiple Copy Simultaneous Search (MCSS) to generate high-performing pharmacophore models in the context of experimentally determined, as well as modeled GPCR structures. Moreover, we have addressed the oft-neglected topic of pharmacophore model selection via development of a cluster-then-predict machine learning workflow. Herein score-based pharmacophore models were generated in experimentally determined and modeled structures of 13 class A GPCR and resulted in pharmacophore models exhibiting high enrichment factors when used to search a database containing 569 class A GPCR ligands. In addition, classification of pharmacophore models with the best performing cluster-then-predict logistic regression classifier resulted in positive predictive values (PPV) of 0.88 and 0.76 for selecting high enrichment pharmacophore models from among those generated in experimentally determined and modeled structures, respectively.
Collapse
Affiliation(s)
| | - Bernie J Daigle
- Departments of Biological Sciences and Computer Science, The University of Memphis, Memphis, TN, 38152, USA
| | - Daniel L Baker
- Department of Chemistry, The University of Memphis, Memphis, TN, 38152, USA
| | - Abby L Parrill
- Department of Chemistry, The University of Memphis, Memphis, TN, 38152, USA.
| |
Collapse
|
2
|
Szwabowski GL, Cole JA, Baker DL, Parrill AL. Structure-based pharmacophore modeling 1. Automated random pharmacophore model generation. J Mol Graph Model 2023; 121:108429. [PMID: 36804368 DOI: 10.1016/j.jmgm.2023.108429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2022] [Revised: 01/18/2023] [Accepted: 02/06/2023] [Indexed: 02/13/2023]
Abstract
Pharmacophores are three-dimensional arrangements of molecular features required for biological activity that are often used in virtual screening efforts to prioritize ligands for experimental testing. G protein-coupled receptors (GPCR) are integral membrane proteins of considerable interest as targets for ligand discovery and drug development. Ligand-based pharmacophore models can be constructed to identify structural commonalities between known bioactive ligands for targets including GPCR. However, structure-based pharmacophores (which only require an experimentally determined or modeled structure for a protein target) have gained more attention to aid in virtual screening efforts as the number of publicly available experimentally determined GPCR structures have increased (140 unique GPCR represented as of October 24, 2022). Thus, the goal of this study was to develop a method of structure-based pharmacophore model generation applicable to ligand discovery for GPCR that have few known ligands. Pharmacophore models were generated within the active sites of 8 class A GPCR crystal structures via automated annotation of 5 randomly selected functional group fragments to sample diverse combinations of pharmacophore features. Each of the 5000 generated pharmacophores was then used to search a database containing active and decoy/inactive compounds for 30 class A GPCR and scored using enrichment factor and goodness-of-hit metrics to assess performance. Application of this method to the set of 8 class A GPCR produced pharmacophore models possessing the theoretical maximum enrichment factor value in both resolved structures (8 of 8 cases) and homology models (7 of 8 cases), indicating that generated pharmacophore models can prove useful in the context of virtual screening.
Collapse
Affiliation(s)
| | - Judith A Cole
- Department of Biological Sciences, The University of Memphis, Memphis, TN, 38152, USA
| | - Daniel L Baker
- Department of Chemistry, The University of Memphis, Memphis, TN, 38152, USA
| | - Abby L Parrill
- Department of Chemistry, The University of Memphis, Memphis, TN, 38152, USA.
| |
Collapse
|
3
|
Jiang S, Feher M, Williams C, Cole B, Shaw DE. AutoPH4: An Automated Method for Generating Pharmacophore Models from Protein Binding Pockets. J Chem Inf Model 2020; 60:4326-4338. [PMID: 32639159 DOI: 10.1021/acs.jcim.0c00121] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Pharmacophore models are widely used in computational drug discovery (e.g., in the virtual screening of drug molecules) to capture essential information about interactions between ligands and a target protein. Generating pharmacophore models from protein structures is typically a manual process, but there has been growing interest in automated pharmacophore generation methods. Automation makes feasible the processing of large numbers of protein conformations, such as those generated by molecular dynamics (MD) simulations, and thus may help achieve the longstanding goal of incorporating protein flexibility into virtual screening workflows. Here, we present AutoPH4, a new automated method for generating pharmacophore models based on protein structures; we show that a virtual screening workflow incorporating AutoPH4 ranks compounds more accurately than any other pharmacophore-based virtual screening workflow for which results on a public benchmark have been reported. The strong performance of the virtual screening workflow indicates that the AutoPH4 component of the workflow generates high-quality pharmacophores, making AutoPH4 promising for use in future virtual screening workflows as well, such as ones that use conformations generated by MD simulations.
Collapse
Affiliation(s)
- Siduo Jiang
- D. E. Shaw Research, New York, New York 10036, United States
| | - Miklos Feher
- D. E. Shaw Research, New York, New York 10036, United States
| | - Chris Williams
- Chemical Computing Group, Montreal, Quebec H3A 2R7, Canada
| | - Brian Cole
- D. E. Shaw Research, New York, New York 10036, United States
| | - David E Shaw
- D. E. Shaw Research, New York, New York 10036, United States.,Department of Biochemistry and Molecular Biophysics, Columbia University, New York, New York 10032, United States
| |
Collapse
|
4
|
Arodola OA, Soliman MES. Quantum mechanics implementation in drug-design workflows: does it really help? Drug Des Devel Ther 2017; 11:2551-2564. [PMID: 28919707 PMCID: PMC5587087 DOI: 10.2147/dddt.s126344] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
The pharmaceutical industry is progressively operating in an era where development costs are constantly under pressure, higher percentages of drugs are demanded, and the drug-discovery process is a trial-and-error run. The profit that flows in with the discovery of new drugs has always been the motivation for the industry to keep up the pace and keep abreast with the endless demand for medicines. The process of finding a molecule that binds to the target protein using in silico tools has made computational chemistry a valuable tool in drug discovery in both academic research and pharmaceutical industry. However, the complexity of many protein-ligand interactions challenges the accuracy and efficiency of the commonly used empirical methods. The usefulness of quantum mechanics (QM) in drug-protein interaction cannot be overemphasized; however, this approach has little significance in some empirical methods. In this review, we discuss recent developments in, and application of, QM to medically relevant biomolecules. We critically discuss the different types of QM-based methods and their proposed application to incorporating them into drug-design and -discovery workflows while trying to answer a critical question: are QM-based methods of real help in drug-design and -discovery research and industry?
Collapse
Affiliation(s)
- Olayide A Arodola
- Department of Pharmaceutical Chemistry, University of KwaZulu-Natal, Durban, South Africa
| | - Mahmoud ES Soliman
- Department of Pharmaceutical Chemistry, University of KwaZulu-Natal, Durban, South Africa
- Department of Pharmaceutical Organic Chemistry, Faculty of Pharmacy, Zagazig University, Egypt
| |
Collapse
|
5
|
Abstract
Over the past two decades, solvent mapping has emerged as a useful tool for identifying hot spots within binding sites on proteins for drug-like molecules and suggesting properties of potential binders. While the experimental technique requires solving multiple crystal structures of a protein in different solvents, computational solvent mapping allows for fast analysis of a protein for potential binding sites and their druggability. Recent advances in genomics, systems biology and interactomics provide a multitude of potential targets for drug development and solvent mapping can provide useful information to help prioritize targets for drug discovery projects. Here, we review various approaches to computational solvent mapping, highlight some key advances and provide our opinion on future directions in the field.
Collapse
|
6
|
Yu W, Lakkaraju SK, Raman EP, Fang L, MacKerell AD. Pharmacophore modeling using site-identification by ligand competitive saturation (SILCS) with multiple probe molecules. J Chem Inf Model 2015; 55:407-20. [PMID: 25622696 DOI: 10.1021/ci500691p] [Citation(s) in RCA: 58] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Receptor-based pharmacophore modeling is an efficient computer-aided drug design technique that uses the structure of the target protein to identify novel leads. However, most methods consider protein flexibility and desolvation effects in a very approximate way, which may limit their use in practice. The Site-Identification by Ligand Competitive Saturation (SILCS) assisted pharmacophore modeling protocol (SILCS-Pharm) was introduced recently to address these issues, as SILCS naturally takes both protein flexibility and desolvation effects into account by using full molecular dynamics simulations to determine 3D maps of the functional group-affinity patterns on a target receptor. In the present work, the SILCS-Pharm protocol is extended to use a wider range of probe molecules including benzene, propane, methanol, formamide, acetaldehyde, methylammonium, acetate and water. This approach removes the previous ambiguity brought by using water as both the hydrogen-bond donor and acceptor probe molecule. The new SILCS-Pharm protocol is shown to yield improved screening results, as compared to the previous approach based on three target proteins. Further validation of the new protocol using five additional protein targets showed improved screening compared to those using common docking methods, further indicating improvements brought by the explicit inclusion of additional feature types associated with the wider collection of probe molecules in the SILCS simulations. The advantage of using complementary features and volume constraints, based on exclusion maps of the protein defined from the SILCS simulations, is presented. In addition, reranking using SILCS-based ligand grid free energies is shown to enhance the diversity of identified ligands for the majority of targets. These results suggest that the SILCS-Pharm protocol will be of utility in rational drug design.
Collapse
Affiliation(s)
- Wenbo Yu
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| | | | | | | | | |
Collapse
|
7
|
Abstract
Pharmacophore modeling incorporates geometric and chemical features of known inhibitors and/or targeted binding sites to rationally identify and design new drug leads. In this study, we have encoded a three-dimensional pharmacophore matching similarity (FMS) scoring function into the structure-based design program DOCK. Validation and characterization of the method are presented through pose reproduction, crossdocking, and enrichment studies. When used alone, FMS scoring dramatically improves pose reproduction success to 93.5% (∼20% increase) and reduces sampling failures to 3.7% (∼6% drop) compared to the standard energy score (SGE) across 1043 protein-ligand complexes. The combined FMS+SGE function further improves success to 98.3%. Crossdocking experiments using FMS and FMS+SGE scoring, for six diverse protein families, similarly showed improvements in success, provided proper pharmacophore references are employed. For enrichment, incorporating pharmacophores during sampling and scoring, in most cases, also yield improved outcomes when docking and rank-ordering libraries of known actives and decoys to 15 systems. Retrospective analyses of virtual screenings to three clinical drug targets (EGFR, IGF-1R, and HIVgp41) using X-ray structures of known inhibitors as pharmacophore references are also reported, including a customized FMS scoring protocol to bias on selected regions in the reference. Overall, the results and fundamental insights gained from this study should benefit the docking community in general, particularly researchers using the new FMS method to guide computational drug discovery with DOCK.
Collapse
Affiliation(s)
- Lingling Jiang
- Department of Applied Mathematics & Statistics, ‡Institute of Chemical Biology & Drug Discovery, §Laufer Center for Physical & Quantitative Biology, Stony Brook University , Stony Brook, New York 11794-3600, United States
| | | |
Collapse
|
8
|
Yu W, Lakkaraju SK, Raman EP, MacKerell AD. Site-Identification by Ligand Competitive Saturation (SILCS) assisted pharmacophore modeling. J Comput Aided Mol Des 2014; 28:491-507. [PMID: 24610239 DOI: 10.1007/s10822-014-9728-0] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2013] [Accepted: 02/04/2014] [Indexed: 12/14/2022]
Abstract
Database screening using receptor-based pharmacophores is a computer-aided drug design technique that uses the structure of the target molecule (i.e. protein) to identify novel ligands that may bind to the target. Typically receptor-based pharmacophore modeling methods only consider a single or limited number of receptor conformations and map out the favorable binding patterns in vacuum or with a limited representation of the aqueous solvent environment, such that they may suffer from neglect of protein flexibility and desolvation effects. Site-Identification by Ligand Competitive Saturation (SILCS) is an approach that takes into account these, as well as other, properties to determine 3-dimensional maps of the functional group-binding patterns on a target receptor (i.e. FragMaps). In this study, a method to use the FragMaps to automatically generate receptor-based pharmacophore models is presented. It converts the FragMaps into SILCS pharmacophore features including aromatic, aliphatic, hydrogen-bond donor and acceptor chemical functionalities. The method generates multiple pharmacophore hypotheses that are then quantitatively ranked using SILCS grid free energies. The pharmacophore model generation protocol is validated using three different protein targets, including using the resulting models in virtual screening. Improved performance and efficiency of the SILCS derived pharmacophore models as compared to published docking studies, as well as a recently developed receptor-based pharmacophore modeling method is shown, indicating the potential utility of the approach in rational drug design.
Collapse
Affiliation(s)
- Wenbo Yu
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD, 21201, USA
| | | | | | | |
Collapse
|
9
|
Joseph-McCarthy D, Campbell AJ, Kern G, Moustakas D. Fragment-Based Lead Discovery and Design. J Chem Inf Model 2014; 54:693-704. [DOI: 10.1021/ci400731w] [Citation(s) in RCA: 100] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Diane Joseph-McCarthy
- Infection Innovative Medicines Unit, AstraZeneca, R&D Boston, 35 Gatehouse Drive, Waltham, Massachusetts 02451, United States
| | - Arthur J. Campbell
- Infection Innovative Medicines Unit, AstraZeneca, R&D Boston, 35 Gatehouse Drive, Waltham, Massachusetts 02451, United States
| | - Gunther Kern
- Infection Innovative Medicines Unit, AstraZeneca, R&D Boston, 35 Gatehouse Drive, Waltham, Massachusetts 02451, United States
| | - Demetri Moustakas
- Infection Innovative Medicines Unit, AstraZeneca, R&D Boston, 35 Gatehouse Drive, Waltham, Massachusetts 02451, United States
| |
Collapse
|
10
|
Zhang X, Wong SE, Lightstone FC. Message passing interface and multithreading hybrid for parallel molecular docking of large databases on petascale high performance computing machines. J Comput Chem 2013; 34:915-27. [PMID: 23345155 DOI: 10.1002/jcc.23214] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2012] [Revised: 10/03/2012] [Accepted: 11/28/2012] [Indexed: 02/06/2023]
Abstract
A mixed parallel scheme that combines message passing interface (MPI) and multithreading was implemented in the AutoDock Vina molecular docking program. The resulting program, named VinaLC, was tested on the petascale high performance computing (HPC) machines at Lawrence Livermore National Laboratory. To exploit the typical cluster-type supercomputers, thousands of docking calculations were dispatched by the master process to run simultaneously on thousands of slave processes, where each docking calculation takes one slave process on one node, and within the node each docking calculation runs via multithreading on multiple CPU cores and shared memory. Input and output of the program and the data handling within the program were carefully designed to deal with large databases and ultimately achieve HPC on a large number of CPU cores. Parallel performance analysis of the VinaLC program shows that the code scales up to more than 15K CPUs with a very low overhead cost of 3.94%. One million flexible compound docking calculations took only 1.4 h to finish on about 15K CPUs. The docking accuracy of VinaLC has been validated against the DUD data set by the re-docking of X-ray ligands and an enrichment study, 64.4% of the top scoring poses have RMSD values under 2.0 Å. The program has been demonstrated to have good enrichment performance on 70% of the targets in the DUD data set. An analysis of the enrichment factors calculated at various percentages of the screening database indicates VinaLC has very good early recovery of actives.
Collapse
Affiliation(s)
- Xiaohua Zhang
- Biosciences & Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Lab, Livermore, CA 94550, USA
| | | | | |
Collapse
|
11
|
Kim ND, Lee YH, Han CK, Ahn SK. Discovery of Novel 11β-HSD1 Inhibitors by Pharmacophore-Based Virtual Screening. B KOREAN CHEM SOC 2012. [DOI: 10.5012/bkcs.2012.33.7.2365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
12
|
Pharmacophore based virtual screening, molecular docking studies to design potent heat shock protein 90 inhibitors. Eur J Med Chem 2011; 46:2937-47. [DOI: 10.1016/j.ejmech.2011.04.018] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2011] [Revised: 03/30/2011] [Accepted: 04/04/2011] [Indexed: 11/19/2022]
|
13
|
Wallach I. Pharmacophore inference and its application to computational drug discovery. Drug Dev Res 2010. [DOI: 10.1002/ddr.20398] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Izhar Wallach
- Department of Computer Science and the Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, Ontario, Canada
| |
Collapse
|
14
|
Cole DC, Olland AM, Jacob J, Brooks J, Bursavich MG, Czerwinski R, DeClercq C, Johnson M, Joseph-McCarthy D, Ellingboe JW, Lin L, Nowak P, Presman E, Strand J, Tam A, Williams CMM, Yao S, Tsao DHH, Fitz LJ. Identification and characterization of acidic mammalian chitinase inhibitors. J Med Chem 2010; 53:6122-8. [PMID: 20666458 DOI: 10.1021/jm100533p] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Acidic mammalian chitinase (AMCase) is a member of the glycosyl hydrolase 18 family (EC 3.2.1.14) that has been implicated in the pathophysiology of allergic airway disease such as asthma. Small molecule inhibitors of AMCase were identified using a combination of high-throughput screening, fragment screening, and virtual screening techniques and characterized by enzyme inhibition and NMR and Biacore binding experiments. X-ray structures of the inhibitors in complex with AMCase revealed that the larger more potent HTS hits, e.g. 5-(4-(2-(4-bromophenoxy)ethyl)piperazine-1-yl)-1H-1,2,4-triazol-3-amine 1, spanned from the active site pocket to a hydrophobic pocket. Smaller fragments identified by FBS occupy both these pockets independently and suggest potential strategies for linking fragments. Compound 1 is a 200 nM AMCase inhibitor which reduced AMCase enzymatic activity in the bronchoalveolar lavage fluid in allergen-challenged mice after oral dosing.
Collapse
Affiliation(s)
- Derek C Cole
- WorldWide Medicinal Chemistry: Inflammation & Immunology, Pfizer Global Research & Development, Cambridge, MA 01240, USA.
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
15
|
Cross JB, Thompson DC, Rai BK, Baber JC, Fan KY, Hu Y, Humblet C. Comparison of several molecular docking programs: pose prediction and virtual screening accuracy. J Chem Inf Model 2009; 49:1455-74. [PMID: 19476350 DOI: 10.1021/ci900056c] [Citation(s) in RCA: 316] [Impact Index Per Article: 21.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular docking programs are widely used modeling tools for predicting ligand binding modes and structure based virtual screening. In this study, six molecular docking programs (DOCK, FlexX, GLIDE, ICM, PhDOCK, and Surflex) were evaluated using metrics intended to assess docking pose and virtual screening accuracy. Cognate ligand docking to 68 diverse, high-resolution X-ray complexes revealed that ICM, GLIDE, and Surflex generated ligand poses close to the X-ray conformation more often than the other docking programs. GLIDE and Surflex also outperformed the other docking programs when used for virtual screening, based on mean ROC AUC and ROC enrichment values obtained for the 40 protein targets in the Directory of Useful Decoys (DUD). Further analysis uncovered general trends in accuracy that are specific for particular protein families. Modifying basic parameters in the software was shown to have a significant effect on docking and virtual screening results, suggesting that expert knowledge is critical for optimizing the accuracy of these methods.
Collapse
Affiliation(s)
- Jason B Cross
- Wyeth Research, Chemical Sciences, Collegeville, Pennsylvania 19426, USA.
| | | | | | | | | | | | | |
Collapse
|
16
|
Challenges of fragment screening. J Comput Aided Mol Des 2009; 23:449-51. [DOI: 10.1007/s10822-009-9293-0] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2009] [Accepted: 06/10/2009] [Indexed: 10/20/2022]
|
17
|
Schubert CR, Stultz CM. The multi-copy simultaneous search methodology: a fundamental tool for structure-based drug design. J Comput Aided Mol Des 2009; 23:475-89. [PMID: 19506805 DOI: 10.1007/s10822-009-9287-y] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2009] [Accepted: 05/20/2009] [Indexed: 10/20/2022]
Abstract
Fragment-based ligand design approaches, such as the multi-copy simultaneous search (MCSS) methodology, have proven to be useful tools in the search for novel therapeutic compounds that bind pre-specified targets of known structure. MCSS offers a variety of advantages over more traditional high-throughput screening methods, and has been applied successfully to challenging targets. The methodology is quite general and can be used to construct functionality maps for proteins, DNA, and RNA. In this review, we describe the main aspects of the MCSS method and outline the general use of the methodology as a fundamental tool to guide the design of de novo lead compounds. We focus our discussion on the evaluation of MCSS results and the incorporation of protein flexibility into the methodology. In addition, we demonstrate on several specific examples how the information arising from the MCSS functionality maps has been successfully used to predict ligand binding to protein targets and RNA.
Collapse
Affiliation(s)
- Christian R Schubert
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | | |
Collapse
|
18
|
Lerner MG, Meagher KL, Carlson HA. Automated clustering of probe molecules from solvent mapping of protein surfaces: new algorithms applied to hot-spot mapping and structure-based drug design. J Comput Aided Mol Des 2008; 22:727-36. [PMID: 18679808 PMCID: PMC2856601 DOI: 10.1007/s10822-008-9231-6] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2008] [Accepted: 07/21/2008] [Indexed: 10/21/2022]
Abstract
Use of solvent mapping, based on multiple-copy minimization (MCM) techniques, is common in structure-based drug discovery. The minima of small-molecule probes define locations for complementary interactions within a binding pocket. Here, we present improved methods for MCM. In particular, a Jarvis-Patrick (JP) method is outlined for grouping the final locations of minimized probes into physical clusters. This algorithm has been tested through a study of protein-protein interfaces, showing the process to be robust, deterministic, and fast in the mapping of protein "hot spots." Improvements in the initial placement of probe molecules are also described. A final application to HIV-1 protease shows how our automated technique can be used to partition data too complicated to analyze by hand. These new automated methods may be easily and quickly extended to other protein systems, and our clustering methodology may be readily incorporated into other clustering packages.
Collapse
Affiliation(s)
- Michael G. Lerner
- Department of Biophysics, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109-1055
| | - Kristin L. Meagher
- Department of Medicinal Chemistry, College of Pharmacy, 418 Church St., University of Michigan, Ann Arbor, Michigan 48109-1065
| | - Heather A. Carlson
- Department of Biophysics, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109-1055
- Department of Medicinal Chemistry, College of Pharmacy, 418 Church St., University of Michigan, Ann Arbor, Michigan 48109-1065
| |
Collapse
|
19
|
Thompson DC, Humblet C, Joseph-McCarthy D. Investigation of MM-PBSA Rescoring of Docking Poses. J Chem Inf Model 2008; 48:1081-91. [DOI: 10.1021/ci700470c] [Citation(s) in RCA: 103] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David C. Thompson
- Wyeth Research Chemical and Screening Sciences 200 Cambridge Park Drive, Cambridge, Massachusetts 02140, and Wyeth Research Chemical and Screening Sciences 865 Ridge Road Princeton, New Jersey 08543
| | - Christine Humblet
- Wyeth Research Chemical and Screening Sciences 200 Cambridge Park Drive, Cambridge, Massachusetts 02140, and Wyeth Research Chemical and Screening Sciences 865 Ridge Road Princeton, New Jersey 08543
| | - Diane Joseph-McCarthy
- Wyeth Research Chemical and Screening Sciences 200 Cambridge Park Drive, Cambridge, Massachusetts 02140, and Wyeth Research Chemical and Screening Sciences 865 Ridge Road Princeton, New Jersey 08543
| |
Collapse
|
20
|
Evensen E, Joseph-McCarthy D, Weiss GA, Schreiber SL, Karplus M. Ligand design by a combinatorial approach based on modeling and experiment: application to HLA-DR4. J Comput Aided Mol Des 2007; 21:395-418. [PMID: 17657565 DOI: 10.1007/s10822-007-9119-x] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2006] [Accepted: 04/19/2007] [Indexed: 01/02/2023]
Abstract
Combinatorial synthesis and large scale screening methods are being used increasingly in drug discovery, particularly for finding novel lead compounds. Although these "random" methods sample larger areas of chemical space than traditional synthetic approaches, only a relatively small percentage of all possible compounds are practically accessible. It is therefore helpful to select regions of chemical space that have greater likelihood of yielding useful leads. When three-dimensional structural data are available for the target molecule this can be achieved by applying structure-based computational design methods to focus the combinatorial library. This is advantageous over the standard usage of computational methods to design a small number of specific novel ligands, because here computation is employed as part of the combinatorial design process and so is required only to determine a propensity for binding of certain chemical moieties in regions of the target molecule. This paper describes the application of the Multiple Copy Simultaneous Search (MCSS) method, an active site mapping and de novo structure-based design tool, to design a focused combinatorial library for the class II MHC protein HLA-DR4. Methods for the synthesizing and screening the computationally designed library are presented; evidence is provided to show that binding was achieved. Although the structure of the protein-ligand complex could not be determined, experimental results including cross-exclusion of a known HLA-DR4 peptide ligand (HA) by a compound from the library. Computational model building suggest that at least one of the ligands designed and identified by the methods described binds in a mode similar to that of native peptides.
Collapse
Affiliation(s)
- Erik Evensen
- Committee on Higher Degrees in Biophysics, Harvard University, Cambridge, MA, USA
| | | | | | | | | |
Collapse
|
21
|
Damm KL, Carlson HA. Exploring Experimental Sources of Multiple Protein Conformations in Structure-Based Drug Design. J Am Chem Soc 2007; 129:8225-35. [PMID: 17555316 DOI: 10.1021/ja0709728] [Citation(s) in RCA: 119] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Receptor flexibility must be incorporated into structure-based drug design in order to portray a more accurate representation of a protein in solution. Our approach is to generate pharmacophore models based on multiple conformations of a protein and is very similar to solvent mapping of hot spots. Previously, we had success using computer-generated conformations of apo human immunodeficiency virus-1 protease (HIV-1p). Here, we examine the use of an NMR ensemble versus a collection of crystal structures, and we compare back to our previous study based on computer-generated conformations. To our knowledge, this is the first direct comparison of an NMR ensemble and a collection of crystal structures to incorporate protein flexibility in structure-based drug design. To provide an accurate comparison between the experimental sources, we used bound structures for our multiple protein structure (MPS) pharmacophore models. The models from an NMR ensemble and a collection of crystal structures were both able to discriminate known HIV-1p inhibitors from decoy molecules and displayed superior performance over models created from single conformations of the protein. Although the active-site conformations were already predefined by bound ligands, the use of MPS allows us to overcome the cross-docking problem and generate a model that does not simply reproduce the chemical characteristics of a specific ligand class. We show that there is more structural variation between 28 structures in an NMR ensemble than 90 crystal structures bound to a variety of ligands. MPS models from both sources performed well, but the model determined using the NMR ensemble appeared to be the most general yet accurate representation of the active site. This work encourages the use of NMR models in structure-based design.
Collapse
Affiliation(s)
- Kelly L Damm
- Department of Medicinal Chemistry, University of Michigan, Ann Arbor, MI 48109-1065, USA
| | | |
Collapse
|
22
|
Bowman AL, Lerner MG, Carlson HA. Protein flexibility and species specificity in structure-based drug discovery: dihydrofolate reductase as a test system. J Am Chem Soc 2007; 129:3634-40. [PMID: 17335207 DOI: 10.1021/ja068256d] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Anna L Bowman
- Department of Medicinal Chemistry and Biophysics Research Division, University of Michigan, Ann Arbor, Michigan 48109-1065, USA
| | | | | |
Collapse
|
23
|
McConnell O, Bach A, Balibar C, Byrne N, Cai Y, Carter G, Chlenov M, Di L, Fan K, Goljer I, He Y, Herold D, Kagan M, Kerns E, Koehn F, Kraml C, Marathias V, Marquez B, McDonald L, Nogle L, Petucci C, Schlingmann G, Tawa G, Tischler M, Williamson RT, Sutherland A, Watts W, Young M, Zhang MY, Zhang Y, Zhou D, Ho D. Enantiomeric separation and determination of absolute stereochemistry of asymmetric molecules in drug discovery—Building chiral technology toolboxes. Chirality 2007; 19:658-82. [PMID: 17390370 DOI: 10.1002/chir.20399] [Citation(s) in RCA: 82] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The application of Chiral Technology, or the (extensive) use of techniques or tools for the determination of absolute stereochemistry and the enantiomeric or chiral separation of racemic small molecule potential lead compounds, has been critical to successfully discovering and developing chiral drugs in the pharmaceutical industry. This has been due to the rapid increase over the past 10-15 years in potential drug candidates containing one or more asymmetric centers. Based on the experiences of one pharmaceutical company, a summary of the establishment of a Chiral Technology toolbox, including the implementation of known tools as well as the design, development, and implementation of new Chiral Technology tools, is provided.
Collapse
Affiliation(s)
- Oliver McConnell
- Wyeth Research, Chemical and Screening Sciences, Collegeville, PA 19426, USA.
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
24
|
Hoffmann M, Eitner K, von Grotthuss M, Rychlewski L, Banachowicz E, Grabarkiewicz T, Szkoda T, Kolinski A. Three dimensional model of severe acute respiratory syndrome coronavirus helicase ATPase catalytic domain and molecular design of severe acute respiratory syndrome coronavirus helicase inhibitors. J Comput Aided Mol Des 2006; 20:305-19. [PMID: 16972168 PMCID: PMC7088412 DOI: 10.1007/s10822-006-9057-z] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2005] [Accepted: 07/17/2006] [Indexed: 10/29/2022]
Abstract
The modeling of the severe acute respiratory syndrome coronavirus helicase ATPase catalytic domain was performed using the protein structure prediction Meta Server and the 3D Jury method for model selection, which resulted in the identification of 1JPR, 1UAA and 1W36 PDB structures as suitable templates for creating a full atom 3D model. This model was further utilized to design small molecules that are expected to block an ATPase catalytic pocket thus inhibit the enzymatic activity. Binding sites for various functional groups were identified in a series of molecular dynamics calculation. Their positions in the catalytic pocket were used as constraints in the Cambridge structural database search for molecules having the pharmacophores that interacted most strongly with the enzyme in a desired position. The subsequent MD simulations followed by calculations of binding energies of the designed molecules were compared to ATP identifying the most successful candidates, for likely inhibitors - molecules possessing two phosphonic acid moieties at distal ends of the molecule.
Collapse
Affiliation(s)
- Marcin Hoffmann
- BioInfoBank Institute, ul. Limanowskiego 24A, 60-744 Poznan, Poland.
| | | | | | | | | | | | | | | |
Collapse
|
25
|
Adcock SA, McCammon JA. Molecular dynamics: survey of methods for simulating the activity of proteins. Chem Rev 2006; 106:1589-615. [PMID: 16683746 PMCID: PMC2547409 DOI: 10.1021/cr040426m] [Citation(s) in RCA: 757] [Impact Index Per Article: 42.1] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Stewart A. Adcock
- NSF Center for Theoretical Biological Physics, Department of Chemistry and Biochemistry, University of California at San Diego, 9500 Gilman Drive, La Jolla, California 92093-0365
| | - J. Andrew McCammon
- NSF Center for Theoretical Biological Physics, Department of Chemistry and Biochemistry, University of California at San Diego, 9500 Gilman Drive, La Jolla, California 92093-0365
| |
Collapse
|
26
|
Milac AL, Avram S, Petrescu AJ. Evaluation of a neural networks QSAR method based on ligand representation using substituent descriptors. Application to HIV-1 protease inhibitors. J Mol Graph Model 2005; 25:37-45. [PMID: 16325439 DOI: 10.1016/j.jmgm.2005.09.014] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2005] [Revised: 06/17/2005] [Accepted: 09/29/2005] [Indexed: 11/18/2022]
Abstract
We present here a neural networks method designed to predict biological activity based on a local representation of the ligand. The compounds of the series are represented by a vector mapping for each of four substituent properties: volume, log P, dipole moment and a simple 'steric' parameter relating to its shape. This ligand representation was tested using neural networks on a set of 42 cyclic-urea derivatives, inhibiting HIV-1 protease. The leave-one-out cross-validation using all descriptors in the input gave a correlation factor between prediction and experiment of 0.76 for the overall set and 0.88 when three outliers were left out. To rank the significance of the four descriptors, we further tested all combinations of two and three parameters for each substituent, using two disjunctive testing sets of five inhibitors. In these sets, vectors with extreme descriptor values were used either in the training or the testing set (sets A and B, respectively). The method is a very good interpolator (set A, 95+/-2% accuracy) but a less effective extrapolator (set B, 85+/-2% accuracy). Generally, the combinations including the 'steric' parameter predict better than average, while those containing the volume are less effective. The best prediction, 98.8+/-1.2%, was obtained when log P, the dipole and the steric parameter were used on set A. At the opposite end, the lowest ranked descriptor set was obtained when replacing log P with the volume, giving 92.3+/-6.7% accuracy over the set A.
Collapse
Affiliation(s)
- Adina-Luminiţa Milac
- Institute of Biochemistry, Splaiul Independenţei 296, Sector 6, Bucharest, Romania
| | | | | |
Collapse
|
27
|
Arnold JR, Burdick KW, Pegg SCH, Toba S, Lamb ML, Kuntz ID. SitePrint: three-dimensional pharmacophore descriptors derived from protein binding sites for family based active site analysis, classification, and drug design. ACTA ACUST UNITED AC 2005; 44:2190-8. [PMID: 15554689 DOI: 10.1021/ci049814f] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Integrating biological and chemical information is one key task in drug discovery, and one approach to attaining this goal is via three-dimensional pharmacophore descriptors derived from protein binding sites. The SitePrint program generates, aligns, scores, and classifies three-dimensional pharmacophore descriptors, active site grids, and ligand surfaces. The descriptors are formed from molecular fragments that have been docked, minimized, filtered, and clustered in protein active sites. The descriptors have geometric coordinates derived from the fragment positions, and they capture the shape, electrostatics, locations, and angles of entry into pockets of the recognition sites: they also provide a direct link to databases of organic molecules. The descriptors have been shown to be robust with respect to small changes in protein structure observed when multiple compounds are cocrystallized in a protein. Five aligned thrombin cocrystals with an average core alpha-carbon RMSD of 0.7 A gave three-dimensional pharmacophore descriptors with an average RMSD of 1.1 A. On a larger test set, alignment and scoring of the descriptors using clique-based alignment, and a best first search strategy with an adapted forward-looking Ullmann heuristic was able to select the global minimum three-dimensional alignment in twenty-nine out of thirty cases in less than one CPU second on a workstation. A protein family based analysis was then performed to demonstrate the usefulness of the method in producing a correlation of active site pharmacophore descriptors to protein function. Each protein in a test set of thirty was assigned membership to a family based on computed active site similarity to the following families: kinases, nuclear receptors, the aspartyl, cysteine, serine, and metallo proteases. This method of classifying proteins is complementary to approaches based on sequence or fold homology. The values within protein families for correctly assigning membership of a protein to a family ranged from 25% to 80%.
Collapse
Affiliation(s)
- James R Arnold
- Department of Pharmaceutical Chemistry, School of Pharmacy, University of California San Francisco, Box 2240, N474-A Genentech Hall, San Francisco, California 94143-2240, USA
| | | | | | | | | | | |
Collapse
|
28
|
Joseph-McCarthy D. Chapter 12 Structure-Based Lead Optimization. ACTA ACUST UNITED AC 2005. [DOI: 10.1016/s1574-1400(05)01012-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
|
29
|
Abstract
Receptor-based virtual screening has become a viable source of novel leads in the pharmaceutical industry. The rapidly growing availability of structural information across protein families, the accessibility to increased computational power at affordable cost, as well as an improved understanding on how to effectively apply virtual screening technologies has contributed to their emergence. Nonetheless, continued improvement in the accuracy of scoring functions and a greater understanding of protein mobility is critical to advance the technology further.
Collapse
Affiliation(s)
- Juan C Alvarez
- Chemical and Screening Sciences, Wyeth Research, 200 Cambridge Park Drive, Cambridge Massachusetts 02140, USA.
| |
Collapse
|
30
|
Shao D, Berrodin TJ, Manas E, Hauze D, Powers R, Bapat A, Gonder D, Winneker RC, Frail DE. Identification of novel estrogen receptor alpha antagonists. J Steroid Biochem Mol Biol 2004; 88:351-60. [PMID: 15145444 DOI: 10.1016/j.jsbmb.2004.01.007] [Citation(s) in RCA: 51] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/24/2003] [Accepted: 01/16/2004] [Indexed: 11/26/2022]
Abstract
We have identified novel estrogen receptor alpha (ERalpha) antagonists using both cell-based and computer-based virtual screening strategies. A mammalian two-hybrid screen was used to select compounds that disrupt the interaction between the ERalpha ligand binding domain (LBD) and the coactivator SRC-3. A virtual screen was designed to select compounds that fit onto the LxxLL peptide-binding surface of the receptor, based on the X-ray crystal structure of the ERalpha LBD complexed with a LxxLL peptide. All selected compounds effectively inhibited 17-beta-estradiol induced coactivator recruitment with potency ranging from nano-molar to micromolar. However, in contrast to classical ER antagonists, these novel inhibitors poorly displace estradiol in the ER-ligand competition assay. Nuclear magnetic resonance (NMR) suggested direct binding of these compounds to the receptors pre-complexed with estradiol and further demonstrated that no estradiol displacement occurred. Partial proteolytic enzyme digestion revealed that, when compared with 17-beta-estradiol- and 4 hydroxy-tamoxifen (4-OHT) bound receptors, at least one of these compounds might induce a unique receptor conformation. These small molecules may represent new classes of ER antagonists, and may have the potential to provide an alternative for the current anti-estrogen therapy.
Collapse
Affiliation(s)
- Dalei Shao
- Women's Health and Bone, Wyeth Research, 500 Arcola Road, RN2294, Collegeville, PA 19426, USA.
| | | | | | | | | | | | | | | | | |
Collapse
|
31
|
Joseph-McCarthy D, Thomas BE, Belmarsh M, Moustakas D, Alvarez JC. Pharmacophore-based molecular docking to account for ligand flexibility. Proteins 2003; 51:172-88. [PMID: 12660987 DOI: 10.1002/prot.10266] [Citation(s) in RCA: 64] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Rapid computational mining of large 3D molecular databases is central to generating new drug leads. Accurate virtual screening of large 3D molecular databases requires consideration of the conformational flexibility of the ligand molecules. Ligand flexibility can be included without prohibitively increasing the search time by docking ensembles of precomputed conformers from a conformationally expanded database. A pharmacophore-based docking method whereby conformers of the same or different molecules are overlaid by their largest 3D pharmacophore and simultaneously docked by partial matches to that pharmacophore is presented. The method is implemented in DOCK 4.0.
Collapse
Affiliation(s)
- Diane Joseph-McCarthy
- Wyeth Research, Biological Chemistry Department, Cambridge, Massachusetts 02140, USA.
| | | | | | | | | |
Collapse
|