The endocannabinoid anandamide (AEA) is an antinociceptive lipid that is inactivated through cellular uptake and subsequent catabolism by fatty acid amide hydrolase (FAAH). Fatty acid binding ...proteins (FABPs) are intracellular carriers that deliver AEA and related N-acylethanolamines (NAEs) to FAAH for hydrolysis. The mammalian brain expresses three FABP subtypes: FABP3, FABP5, and FABP7. Recent work from our group has revealed that pharmacological inhibition of FABPs reduces inflammatory pain in mice. The goal of the current work was to explore the effects of FABP inhibition upon nociception in diverse models of pain. We developed inhibitors with differential affinities for FABPs to elucidate the subtype(s) that contributes to the antinociceptive effects of FABP inhibitors. Inhibition of FABPs reduced nociception associated with inflammatory, visceral, and neuropathic pain. The antinociceptive effects of FABP inhibitors mirrored their affinities for FABP5, while binding to FABP3 and FABP7 was not a predictor of in vivo efficacy. The antinociceptive effects of FABP inhibitors were mediated by cannabinoid receptor 1 (CB1) and peroxisome proliferator-activated receptor alpha (PPARα) and FABP inhibition elevated brain levels of AEA, providing the first direct evidence that FABPs regulate brain endocannabinoid tone. These results highlight FABPs as novel targets for the development of analgesic and anti-inflammatory therapeutics.
Prostate cancer is a leading cause of cancer-related deaths in men in the United States. Although treatable when detected early, prostate cancer commonly transitions to an aggressive ...castration-resistant metastatic state. While taxane chemotherapeutics such as docetaxel are mainstay treatment options for prostate cancer, taxane resistance often develops. Fatty acid binding protein 5 (FABP5) is an intracellular lipid chaperone that is upregulated in advanced prostate cancer and is implicated as a key driver of its progression. The recent demonstration that FABP5 inhibitors produce synergistic inhibition of tumor growth when combined with taxane chemotherapeutics highlights the possibility that FABP5 may regulate other features of taxane function, including resistance. Employing taxane-resistant DU145-TXR cells and a combination of cytotoxicity, apoptosis, and cell cycle assays, our findings demonstrate that FABP5 knockdown sensitizes the cells to docetaxel. In contrast, docetaxel potency was unaffected by FABP5 knockdown in taxane-sensitive DU145 cells. Taxane-resistance in DU145-TXR cells stems from upregulation of the P-glycoprotein ATP binding cassette subfamily B member 1 (ABCB1). Expression analyses and functional assays confirmed that FABP5 knockdown in DU145-TXR cells markedly reduced ABCB1 expression and activity, respectively. Our study demonstrates a potential new function for FABP5 in regulating taxane sensitivity and the expression of a major P-glycoprotein efflux pump in prostate cancer cells.
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.
With an increasing interest in RNA therapeutics and for targeting RNA to treat disease, there is a need for the tools used in protein-based drug design, particularly DOCKing algorithms, to be ...extended or adapted for nucleic acids. Here, we have compiled a test set of RNA-ligand complexes to validate the ability of the DOCK suite of programs to successfully recreate experimentally determined binding poses. With the optimized parameters and a minimal scoring function, 70% of the test set with less than seven rotatable ligand bonds and 26% of the test set with less than 13 rotatable bonds can be successfully recreated within 2 A heavy-atom RMSD. When DOCKed conformations are rescored with the implicit solvent models AMBER generalized Born with solvent-accessible surface area (GB/SA) and Poisson-Boltzmann with solvent-accessible surface area (PB/SA) in combination with explicit water molecules and sodium counterions, the success rate increases to 80% with PB/SA for less than seven rotatable bonds and 58% with AMBER GB/SA and 47% with PB/SA for less than 13 rotatable bonds. These results indicate that DOCK can indeed be useful for structure-based drug design aimed at RNA. Our studies also suggest that RNA-directed ligands often differ from typical protein-ligand complexes in their electrostatic properties, but these differences can be accommodated through the choice of potential function. In addition, in the course of the study, we explore a variety of newly added DOCK functions, demonstrating the ease with which new functions can be added to address new scientific questions.
A database consisting of 780 ligand−receptor complexes, termed SB2010, has been derived from the Protein Databank to evaluate the accuracy of docking protocols for regenerating bound ligand ...conformations. The goal is to provide easily accessible community resources for development of improved procedures to aid virtual screening for ligands with a wide range of flexibilities. Three core experiments using the program DOCK, which employ rigid (RGD), fixed anchor (FAD), and flexible (FLX) protocols, were used to gauge performance by several different metrics: (1) global results, (2) ligand flexibility, (3) protein family, and (4) cross-docking. Global spectrum plots of successes and failures vs rmsd reveal well-defined inflection regions, which suggest the commonly used 2 Å criteria is a reasonable choice for defining success. Across all 780 systems, success tracks with the relative difficulty of the calculations: RGD (82.3%) > FAD (78.1%) > FLX (63.8%). In general, failures due to scoring strongly outweigh those due to sampling. Subsets of SB2010 grouped by ligand flexibility (7-or-less, 8-to-15, and 15-plus rotatable bonds) reveal that success degrades linearly for FAD and FLX protocols, in contrast to RGD, which remains constant. Despite the challenges associated with FLX anchor orientation and on-the-fly flexible growth, success rates for the 7-or-less (74.5%) and, in particular, the 8-to-15 (55.2%) subset are encouraging. Poorer results for the very flexible 15-plus set (39.3%) indicate substantial room for improvement. Family-based success appears largely independent of ligand flexibility, suggesting a strong dependence on the binding site environment. For example, zinc-containing proteins are generally problematic, despite moderately flexible ligands. Finally, representative cross-docking examples, for carbonic anhydrase, thermolysin, and neuraminidase families, show the utility of family-based analysis for rapid identification of particularly good or bad docking trends, and the type of failures involved (scoring/sampling), which will likely be of interest to researchers making specific receptor choices for virtual screening. SB2010 is available for download at http://rizzolab.org.
We report unrestrained, all-atom molecular dynamics simulations of HIV-1 protease that sample large conformational changes of the active site flaps. In particular, the unliganded protease undergoes ...multiple conversions between the "closed" and "semiopen" forms observed in crystal structures of inhibitor-bound and unliganded protease, respectively, including reversal of flap "handedness." Simulations in the presence of a cyclic urea inhibitor yield stable closed flaps. Furthermore, we observe several events in which the flaps of the unliganded protease open to a much greater degree than observed in crystal structures and subsequently return to the semiopen state. Our data strongly support the hypothesis that the unliganded protease predominantly populates the semiopen conformation, with closed and fully open structures being a minor component of the overall ensemble. The results also provide a model for the flap opening and closing that is considered to be essential to enzyme function.
Fatty acid binding proteins (FABPs), in particular FABP5 and FABP7, have recently been identified by us as intracellular transporters for the endocannabinoid anandamide (AEA). Furthermore, animal ...studies by others have shown that elevated levels of endocannabinoids resulted in beneficial pharmacological effects on stress, pain and inflammation and also ameliorate the effects of drug withdrawal. Based on these observations, we hypothesized that FABP5 and FABP7 would provide excellent pharmacological targets. Thus, we performed a virtual screening of over one million compounds using DOCK and employed a novel footprint similarity scoring function to identify lead compounds with binding profiles similar to oleic acid, a natural FABP substrate. Forty-eight compounds were purchased based on their footprint similarity scores (FPS) and assayed for biological activity against purified human FABP5 employing a fluorescent displacement-binding assay. Four compounds were found to exhibit approximately 50% inhibition or greater at 10 µM, as good as or better inhibitors of FABP5 than BMS309403, a commercially available inhibitor. The most potent inhibitor, γ-truxillic acid 1-naphthyl ester (ChemDiv 8009-2334), was determined to have K(i) value of 1.19±0.01 µM. Accordingly a novel α-truxillic acid 1-naphthyl mono-ester (SB-FI-26) was synthesized and assayed for its inhibitory activity against FABP5, wherein SB-FI-26 exhibited strong binding (K(i) 0.93±0.08 µM). Additionally, we found SB-FI-26 to act as a potent anti-nociceptive agent with mild anti-inflammatory activity in mice, which strongly supports our hypothesis that the inhibition of FABPs and subsequent elevation of anandamide is a promising new approach to drug discovery. Truxillic acids and their derivatives were also shown by others to have anti-inflammatory and anti-nociceptive effects in mice and to be the active component of Chinese a herbal medicine (Incarvillea sinensis) used to treat rheumatism and pain in humans. Our results provide a likely mechanism by which these compounds exert their effects.
Structure-based methods that employ principles of de novo design can be used to construct small organic molecules from scratch using pre-existing fragment libraries to sample chemical space and are ...an important class of computational algorithms for drug-lead discovery. Here, we present a powerful new design method for DOCK6 that employs a Descriptor-Driven De Novo strategy (termed D3N) in which user-defined cheminformatics descriptors (and their target ranges) are calculated at each layer of growth using the open-source toolkit RDKit. The objective is to tailor ligand growth toward desirable regions of chemical space. The approach was extensively validated through: (1) comparison of cheminformatics descriptors computed using the new DOCK6/RDKit interface versus the standard Python/RDKit installation, (2) examination of descriptor distributions generated using D3N growth under different conditions (target ranges and environments), and (3) construction of ligands with very tight (pinpoint) descriptor ranges using clinically relevant compounds as a reference. Our testing confirms that the new DOCK6/RDKit integration is robust, showcases how the new D3N routines can be used to direct sampling around user-defined chemical spaces, and highlights the utility of on-the-fly descriptor calculations for ligand design to important drug targets.
False negative docking outcomes for highly symmetric molecules are a barrier to the accurate evaluation of docking programs, scoring functions, and protocols. This work describes an implementation of ...a symmetry-corrected root-mean-square deviation (RMSD) method into the program DOCK based on the Hungarian algorithm for solving the minimum assignment problem, which dynamically assigns atom correspondence in molecules with symmetry. The algorithm adds only a trivial amount of computation time to the RMSD calculations and is shown to increase the reported overall docking success rate by approximately 5% when tested over 1043 receptor–ligand systems. For some families of protein systems the results are even more dramatic, with success rate increases up to 16.7%. Several additional applications of the method are also presented including as a pairwise similarity metric to compare molecules during de novo design, as a scoring function to rank-order virtual screening results, and for the analysis of trajectories from molecular dynamics simulation. The new method, including source code, is available to registered users of DOCK6 (http://dock.compbio.ucsf.edu).