Structure-based drug design: exploring the proper filling of apolar pockets at enzyme active sites

The proper filling of apolar pockets at enzyme active sites is central for increasing binding activity and selectivity of hits and leads in medicinal chemistry. In our structure-based design approach toward the generation of potent enzyme inhibitors, we encountered a variety of challenges in gaining suitable binding affinity from the occupation of such pockets. We summarize them here for the first time. A fluorine scan of tricyclic thrombin inhibitors led to the discovery of favorable orthogonal dipolar C-F...CO interactions. Efficient cation-pi interactions were established in the S4 pocket of factor Xa, another serine protease from the blood coagulation cascade. Changing from mono- to bisubstrate inhibitors of catechol O-methyltransferase, a target in the L-Dopa-based treatment of Parkinson's disease, enabled the full exploitation of a previously unexplored hydrophobic pocket. Conformational preorganization of a pocket at an enzyme active site is crucial for harvesting binding affinity. This is demonstrated for two enzymes from the nonmevalonate pathway of isoprenoid biosynthesis, IspE and IspF, which are pursued as antimalarial targets. Disrupting crystallographically defined water networks on the way into a pocket might cost all of the binding free enthalpy gained from its occupation, as revealed in studies with tRNA-guanine transglycosylase, a target against shigellosis. Investigations of the active site of plasmepsin II, another antimalarial target, showed that principles for proper apolar cavity filling, originally developed for synthetic host-guest systems, are also applicable to enzyme environments.     

Structure-based combinatorial library design: discovery of non-peptidic inhibitors of caspases 3 and 8

Structure-based design of a combinatorial array was carried out in order to identify non-peptidic thiomethylketone inhibitors of caspases 3 and 8. Five compounds from the designed array were active against caspase 3, and two were active against caspase 8. A 2.5-Å resolution co-crystal structure of caspase 3 and a thiomethylketone array member is reported. The structure-based design strategy has proved useful for identifying caspase inhibitors.     

Defining Cdk5 ligand chemical space with small molecule inhibitors of tau phosphorylation

Cyclin-dependent kinase 5 (Cdk5) is widely viewed as a possible target for a wide variety of neurological disorders. One pathological role attributed to Cdk5 is the abnormal phosphorylation of tau that may lead to the neuronal inclusions known as neurofibrillary tangles. A high through-put screen for inhibitors of Cdk5-mediated phosphorylation of tau resulted in three compounds with distinct mechanisms of action. One compound is competitive with ATP and has a high affinity for the Cdk5 ATP binding pocket. The second compound also competes with ATP, is noncompetitive with tau, and (uniquely among this class of inhibitors) displaces adjacent amino acid residues to make room for the nitrophenyl group. A third compound did not compete with ATP, but did compete with tau at low concentrations of tau. The SAR and charge optimization derived from cocrystals of the two ATP competitors along with cocrystals of three other ATP competitors map out the importance of filling and properly charging different regions of the ATP binding pocket. Taken together, this analysis shows how the structure of Cdk5 constrains the space of potential inhibitors and reveals a pocket unfilled in all of the structures. These leads could be a starting point for structure-based drug design of more potent and selective inhibitors.     

In search of new lead compounds for trypanosomiasis drug design: A protein structure-based linked-fragment approach

Summary A modular method for pursuing structure-based inhibitor design in the framework of a design cycle is presented. The approach entails four stages: (1) a design pathway is defined in the three-dimensional structure of a target protein; (2) this pathway is divided into subregions; (3) complementary building blocks, also called fragments, are designed in each subregion; complementarity is defined in terms of shape, hydrophobicity, hydrogen bond properties and electrostatics; and (4) fragments from different subregions are linked into potential lead compounds. Stages (3) and (4) are qualitatively guided by force-field calculations. In addition, the designed fragments serve as entries for retrieving existing compounds from chemical databases. This linked-fragment approach has been applied in the design of potentially selective inhibitors of triosephosphate isomerase from Trypanosoma brucei, the causative agent of sleeping sickness.     

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.     

Structure-based design and in-parallel synthesis of inhibitors of AmpC beta-lactamase.

BACKGROUND: Group I beta-lactamases are a major cause of antibiotic resistance to beta-lactams such as penicillins and cephalosporins. These enzymes are only modestly affected by classic beta-lactam-based inhibitors, such as clavulanic acid. Conversely, small arylboronic acids inhibit these enzymes at sub-micromolar concentrations. Structural studies suggest these inhibitors bind to a well-defined cleft in the group I beta-lactamase AmpC; this cleft binds the ubiquitous R1 side chain of beta-lactams. Intriguingly, much of this cleft is left unoccupied by the small arylboronic acids. RESULTS: To investigate if larger boronic acids might take advantage of this cleft, structure-guided in-parallel synthesis was used to explore new inhibitors of AmpC. Twenty-eight derivatives of the lead compound, 3-aminophenylboronic acid, led to an inhibitor with 80-fold better binding (2; K(i) 83 nM). Molecular docking suggested orientations for this compound in the R1 cleft. Based on the docking results, 12 derivatives of 2 were synthesized, leading to inhibitors with K(i) values of 60 nM and with improved solubility. Several of these inhibitors reversed the resistance of nosocomial Gram-positive bacteria, though they showed little activity against Gram-negative bacteria. The X-ray crystal structure of compound 2 in complex with AmpC was subsequently determined to 2.1 A resolution. The placement of the proximal two-thirds of the inhibitor in the experimental structure corresponds with the docked structure, but a bond rotation leads to a distinctly different placement of the distal part of the inhibitor. In the experimental structure, the inhibitor interacts with conserved residues in the R1 cleft whose role in recognition has not been previously explored. CONCLUSIONS: Combining structure-based design with in-parallel synthesis allowed for the rapid exploration of inhibitor functionality in the R1 cleft of AmpC. The resulting inhibitors differ considerably from beta-lactams but nevertheless inhibit the enzyme well. The crystal structure of 2 (K(i) 83 nM) in complex with AmpC may guide exploration of a highly conserved, largely unexplored cleft, providing a template for further design against AmpC beta-lactamase.     

The Kynurenine Pathway and Kynurenine 3-Monooxygenase Inhibitors

Under normal physiological conditions, the kynurenine pathway (KP) plays a critical role in generating cellular energy and catabolizing tryptophan. Under inflammatory conditions, however, there is an upregulation of the KP enzymes, particularly kynurenine 3-monooxygenase (KMO). KMO has garnered much attention due to its production of toxic metabolites that have been implicated in many diseases and disorders. With many of these illnesses having an inadequate or modest treatment, there exists a need to develop KMO inhibitors that reduce the production of these toxic metabolites. Though prior efforts to find an appropriate KMO inhibitor were unpromising, the development of a KMO crystal structure has provided the opportunity for a rational structure-based design in the development of inhibitors. Therefore, the purpose of this review is to describe the kynurenine pathway, the kynurenine 3-monooxygenase enzyme, and KMO inhibitors and their potential candidacy for clinical use.     

Structure-guided fragment-based<Emphasis Type="Italic"> in silico</Emphasis> drug design of dengue protease inhibitors

An in silico fragment-based drug design approach was devised and applied towards the identification of small molecule inhibitors of the dengue virus (DENV) NS2B-NS3 protease. Currently, no DENV protease co-crystal structure with bound inhibitor and fully formed substrate binding site is available. Therefore a homology model of DENV NS2B-NS3 protease was generated employing a multiple template spatial restraints method and used for structure-based design. A library of molecular fragments was derived from the ZINC screening database with help of the retrosynthetic combinatorial analysis procedure (RECAP). 150,000 molecular fragments were docked to the DENV protease homology model and the docking poses were rescored using a target-specific scoring function. High scoring fragments were assembled to small molecule candidates by an implicit linking cascade. The cascade included substructure searching and structural filters focusing on interactions with the S1 and S2 pockets of the protease. The chemical space adjacent to the promising candidates was further explored by neighborhood searching. A total of 23 compounds were tested experimentally and two compounds were discovered to inhibit dengue protease (IC₅₀ = 7.7 μM and 37.9 μM, respectively) and the related West Nile virus protease (IC₅₀ = 6.3 μM and 39.0 μM, respectively). This study demonstrates the successful application of a structure-guided fragment-based in silico drug design approach for dengue protease inhibitors providing straightforward hit generation using a combination of homology modeling, fragment docking, chemical similarity and structural filters.     

Chemoenzymatic synthesis of 4-diphosphocytidyl-2-C-methyl-d-erythritol: a substrate for IspE

Enantiomerically pure 2-C-methyl-d-erythritol 4-phosphate 1 (MEP) is synthesized from 1,2-O-isopropylidene-α-d-xylofuranose via facile benzylation in good yield. Subsequently, 1 is used for enzymatic synthesis of 4-diphosphocytidyl-2-C-methyl-d-erythritol 2 (CDP-ME) using 4-diphosphocytidyl-2-C-methyl-d-erythritol synthase (IspD). The chemoenzymatically synthesized 2 can be used as substrate for assay of IspE and for high throughput screening to identify IspE inhibitors.     

Structure-based design of potent non-peptide MDM2 inhibitors

A successful structure-based design of a class of non-peptide small-molecule MDM2 inhibitors targeting the p53-MDM2 protein-protein interaction is reported. The most potent compound 1d binds to MDM2 protein with a Ki value of 86 nM and is 18 times more potent than a natural p53 peptide (residues 16-27). Compound 1d is potent in inhibition of cell growth in LNCaP prostate cancer cells with wild-type p53 and shows only a weak activity in PC-3 prostate cancer cells with a deleted p53. Importantly, 1d has a minimal toxicity to normal prostate epithelial cells. Our studies provide a convincing example that structure-based strategy can be employed to design highly potent, non-peptide, cell-permeable, small-molecule inhibitors to target protein-protein interaction, which remains a very challenging area in chemical biology and drug design.     

The structural basis for the specificity of pyridinylimidazole inhibitors of p38 MAP kinase

Background: The p38 mitogen-activated protein (MAP) kinase regulates signal transduction in response to environmental stress. Pyridinylimidazole compounds are specific inhibitors of p38 MAP kinase that block the production of the cytokines interleukin-1 β and tumor necrosis factor α, and they are effective in animal models of arthritis, bone resorption and endotoxin shock. These compounds have been useful probes for studying the physiological functions of the p38-mediated MAP kinase pathway.Results: We report the crystal structure of a novel pyridinylimidazole compound complexed with p38 MAP kinase, and we demonstrate that this compound binds to the same site on the kinase as does ATP. Mutagenesis showed that a single residue difference between p38 MAP kinase and other MAP kinases is sufficient to confer selectivity among pyridinylimidazole compounds.Conclusions: Our results reveal how pyridinylimidazole compounds are potent and selective inhibitors of p38 MAP kinase but not other MAP kinases. It should now be possible to design other specific inhibitors of activated p38 MAP kinase using the structure of the nonphosphorylated enzyme.     

Structure-based discovery and in-parallel optimization of novel competitive inhibitors of thymidylate synthase.

BACKGROUND: The substrate sites of enzymes are attractive targets for structure-based inhibitor design. Two difficulties hinder efforts to discover and elaborate new (nonsubstrate-like) inhibitors for these sites. First, novel inhibitors often bind at nonsubstrate sites. Second, a novel scaffold introduces chemistry that is frequently unfamiliar, making synthetic elaboration challenging. RESULTS: In an effort to discover and elaborate a novel scaffold for a substrate site, we combined structure-based screening with in-parallel synthetic elaboration. These techniques were used to find new inhibitors that bound to the folate site of Lactobacillus casei thymidylate synthase (LcTS), an enzyme that is a potential target for proliferative diseases, and is highly studied. The available chemicals directory was screened, using a molecular-docking computer program, for molecules that complemented the three-dimensional structure of this site. Five high-ranking compounds were selected for testing. Activity and docking studies led to a derivative of one of these, dansyltyrosine (Ki 65 microM). Using solid-phase in-parallel techniques 33 derivatives of this lead were synthesized and tested. These analogs are dissimilar to the substrate but bind competitively with it. The most active analog had a Ki of 1.3 microM. The tighter binding inhibitors were also the most specific for LcTS versus related enzymes. CONCLUSIONS: TS can recognize inhibitors that are dissimilar to, but that bind competitively with, the folate substrate. Combining structure-based discovery with in-parallel synthetic techniques allowed the rapid elaboration of this series of compounds. More automated versions of this approach can be envisaged.     

Ligand design by a combinatorial approach based on modeling and experiment: application to HLA-DR4

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.     

Dihydrofolate Reductase Inhibitors: The Pharmacophore as a Guide for Co-Crystal Screening

In this work, co-crystal screening was carried out for two important dihydrofolate reductase (DHFR) inhibitors, trimethoprim (TMP) and pyrimethamine (PMA), and for 2,4-diaminopyrimidine (DAP), which is the pharmacophore of these active pharmaceutical ingredients (API). The isomeric pyridinecarboxamides and two xanthines, theophylline (THEO) and caffeine (CAF), were used as co-formers in the same experimental conditions, in order to evaluate the potential for the pharmacophore to be used as a guide in the screening process. In silico co-crystal screening was carried out using BIOVIA COSMOquick and experimental screening was performed by mechanochemistry and supported by (solid + liquid) binary phase diagrams, infrared spectroscopy (FTIR) and X-ray powder diffraction (XRPD). The in silico prediction of low propensities for DAP, TMP and PMA to co-crystallize with pyridinecarboxamides was confirmed: a successful outcome was only observed for DAP + nicotinamide. Successful synthesis of multicomponent solid forms was achieved for all three target molecules with theophylline, with DAP co-crystals revealing a greater variety of stoichiometries. The crystalline structures of a (1:2) TMP:THEO co-crystal and of a (1:2:1) DAP:THEO:ethyl acetate solvate were solved. This work demonstrated the possible use of the pharmacophore of DHFR inhibitors as a guide for co-crystal screening, recognizing some similar trends in the outcome of association in the solid state and in the molecular aggregation in the co-crystals, characterized by the same supramolecular synthons.     

人类丝氨酸消旋酶的同源模建及其与多肽类抑制剂的分子对接

利用同源模建和分子动力学模拟方法构建了人类丝氨酸消旋酶(hSR)的三维结构, 并利用profile-3D和procheck方法评估了模型的可靠性. 在此基础上用分子对接程序(affinity)将多肽类抑制剂A和B分别与hSR进行对接, 获得了其复合物结构的理论模型. 通过配体与受体之间相互作用能和结构分析给出了此类抑制剂与hSR的具体结合方式, 明确了hSR与此类抑制剂结合时起重要作用的氨基酸残基, 为基于人类丝氨酸消旋酶三维结构的药物设计提供重要的参考信息.     

Exploring experimental sources of multiple protein conformations in structure-based drug design

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.     

Efficient identification of inhibitors targeting the closed active site conformation of the HPRT from Trypanosoma cruzi

BACKGROUND: Currently, only two drugs are recommended for treatment of infection with Trypanosoma cruzi, the etiologic agent of Chagas' disease. These compounds kill the trypomastigote forms of the parasite circulating in the bloodstream, but are relatively ineffective against the intracellular stage of the parasite life cycle. Neither drug is approved by the FDA for use in the US. The hypoxanthine phosphoribosyltransferase (HPRT) from T. cruzi is a possible new target for antiparasite chemotherapy. The crystal structure of the HPRT in a conformation approximating the transition state reveals a closed active site that provides a well-defined target for computational structure-based drug discovery. RESULTS: A flexible ligand docking program incorporating a desolvation correction was used to screen the Available Chemicals Directory for inhibitors targeted to the closed conformation of the trypanosomal HPRT. Of 22 potential inhibitors identified, acquired and tested, 16 yielded K(i)'s between 0.5 and 17 microM versus the substrate phosphoribosylpyrophosphate. Surprisingly, three of eight compounds tested were effective in inhibiting the growth of parasites in infected mammalian cells. CONCLUSIONS: This structure-based docking method provided a remarkably efficient path for the identification of inhibitors targeting the closed conformation of the trypanosomal HPRT. The inhibition constants of the lead inhibitors identified are unusually favorable, and the trypanostatic activity of three of the compounds in cell culture suggests that they may provide useful starting points for drug design for the treatment of Chagas' disease.     

Sphingosine kinase 1 (SK1) allosteric inhibitors that target the dimerization site

The sphingosine kinase 1 (SK1)/sphingosine-1-phosphate (S1P) signaling pathway is a crucial target for numerous human diseases from cancer to cardiovascular diseases. However, available SK1 inhibitors that target the active site suffer from poor potency, selectivity and pharmacokinetic properties. The selectivity issue of the kinases, which share a highly-conserved ATP-pocket, can be overcome by targeting the less-conserved allosteric sites. SK1 is known to function minimally as a dimer; however, the crystal structure of the SK1 dimer has not been determined. In this study, a template-based algorithm implemented in PRISM was used to predict the SK1 dimer structure and then the possible allosteric sites at the dimer interface were determined via SiteMap. These sites were used in a virtual screening campaign that includes an integrated workflow of structure-based pharmacophore modeling, virtual screening, molecular docking, re-screening of common scaffolds to propose a series of compounds with different scaffolds as potential allosteric SK1 inhibitors. Finally, the stability of the SK1-ligand complexes was analyzed by molecular dynamics simulations. As a final outcome, ligand 7 having a 4,9-dihydro-1H-purine scaffold and ligand 12 having a 2,3,4,9-tetrahydro-1H-β-carboline scaffold were found to be potential selective inhibitors for SK1.     

Approaches to solving the rigid receptor problem by identifying a minimal set of flexible residues during ligand docking

BACKGROUND: Using fixed receptor sites derived from high-resolution crystal structures in structure-based drug design does not properly account for ligand-induced enzyme conformational change and imparts a bias into the discovery and design of novel ligands. We sought to facilitate the design of improved drug leads by defining residues most likely to change conformation, and then defining a minimal manifold of possible conformations of a target site for drug design based on a small number of identified flexible residues. RESULTS: The crystal structure of thymidylate synthase from an important pathogenic target Pneumocystis carinii (PcTS) bound to its substrate and the inhibitor, BW1843U89, is reported here and reveals a new conformation with respect to the structure of PcTS bound to substrate and the more conventional antifolate inhibitor, CB3717. We developed an algorithm for determining which residues provide 'soft spots' in the protein, regions where conformational adaptation suggests possible modifications for a drug lead that may yield higher affinity. Remodeling the active site of thymidylate synthase with new conformations for only three residues that were identified with this algorithm yields scores for ligands that are compatible with experimental kinetic data. CONCLUSIONS: Based on the examination of many protein/ligand complexes, we develop an algorithm (SOFTSPOTS) for identifying regions of a protein target that are more likely to accommodate plastically to regions of a drug molecule. Using these indicators we develop a second algorithm (PLASTIC) that provides a minimal manifold of possible conformations of a protein target for drug design, reducing the bias in structure-based drug design imparted by structures of enzymes co-crystallized with inhibitors.     

Targeting a Large Active Site: Structure-Based Design of Nanomolar Inhibitors of Trypanosoma brucei Trypanothione Reductase

Trypanothione reductase (TR) plays a key role in the unique redox metabolism of trypanosomatids, the causative agents of human African trypanosomiasis (HAT), Chagas’ disease, and leishmaniases. Introduction of a new, lean propargylic vector to a known class of TR inhibitors resulted in the strongest reported competitive inhibitor of Trypanosoma (T.) brucei TR, with an inhibition constant K_i of 73 nm , which is fully selective against human glutathione reductase (hGR). The best ligands exhibited in vitro IC₅₀ values (half-maximal inhibitory concentration) against the HAT pathogen, T. brucei rhodesiense, in the mid-nanomolar range, reaching down to 50 nm. X-Ray co-crystal structures confirmed the binding mode of the ligands and revealed the presence of a HEPES buffer molecule in the large active site. Extension of the propargylic vector, guided by structure-based design, to replace the HEPES buffer molecule should give inhibitors with low nanomolar K_i and IC₅₀ values for in vivo studies.