Rational methods for the selection of diverse screening compounds.

Traditionally a pursuit of large pharmaceutical companies, high-throughput screening assays are becoming increasingly common within academic and government laboratories. This shift has been instrumental in enabling projects that have not been commercially viable, such as chemical probe discovery and screening against high-risk targets. Once an assay has been prepared and validated, it must be fed with screening compounds. Crafting a successful collection of small molecules for screening poses a significant challenge. An optimized collection will minimize false positives while maximizing hit rates of compounds that are amenable to lead generation and optimization. Without due consideration of the relevant protein targets and the downstream screening assays, compound filtering and selection can fail to explore the great extent of chemical diversity and eschew valuable novelty. Herein, we discuss the different factors to be considered and methods that may be employed when assembling a structurally diverse compound collection for screening. Rational methods for selecting diverse chemical libraries are essential for their effective use in high-throughput screens.We are grateful for financial support from the MRC, Wellcome Trust, CRUK, EPSRC, BBSRC and Newman Trust.This is the author accepted manuscript. The final version is available from American Chemical Society via http://dx.doi.org/10.1021/cb100420

Identification of Small-Molecule Inhibitors against Meso-2, 6-Diaminopimelate Dehydrogenase from Porphyromonas gingivalis

Species-specific antimicrobial therapy has the potential to combat the increasing threat of antibiotic resistance and alteration of the human microbiome. We therefore set out to demonstrate the beginning of a pathogen-selective drug discovery method using the periodontal pathogen Porphyromonas gingivalis as a model. Through our knowledge of metabolic networks and essential genes we identified a “druggable” essential target, meso-diaminopimelate dehydrogenase, which is found in a limited number of species. We adopted a high-throughput virtual screen method on the ZINC chemical library to select a group of potential small-molecule inhibitors. Meso-diaminopimelate dehydrogenase from P. gingivaliswas first expressed and purified in Escherichia coli then characterized for enzymatic inhibitor screening studies. Several inhibitors with similar structural scaffolds containing a sulfonamide core and aromatic substituents showed dose-dependent inhibition. These compounds were further assayed showing reasonable whole-cell activity and the inhibition mechanism was determined. We conclude that the establishment of this target and screening strategy provides a model for the future development of new antimicrobials

Structure- and Ligand-Based Design of Novel Antimicrobial Agents

The use of computer based techniques in the design of novel therapeutic agents is a rapidly emerging field. Although the drug-design techniques utilized by Computational Medicinal Chemists vary greatly, they can roughly be classified into structure-based and ligand-based approaches. Structure-based methods utilize a solved structure of the design target, protein or DNA, usually obtained by X-ray or NMR methods to design or improve compounds with activity against the target. Ligand-based methods use active compounds with known affinity for a target that may yet be unresolved. These methods include Pharmacophore-based searching for novel active compounds or Quantitative Structure-Activity Relationship (QSAR) studies. The research presented here utilized both structure and ligand-based methods against two bacterial targets: Bacillus anthracis and Mycobacterium tuberculosis. The first part of this thesis details our efforts to design novel inhibitors of the enzyme dihydropteroate synthase from B. anthracis using crystal structures with known inhibitors bound. The second part describes a QSAR study that was performed using a series of novel nitrofuranyl compounds with known, whole-cell, inhibitory activity against M. tuberculosis. Dihydropteroate synthase (DHPS) catalyzes the addition of p-amino benzoic acid (pABA) to dihydropterin pyrophosphate (DHPP) to form pteroic acid as a key step in bacterial folate biosynthesis. It is the traditional target of the sulfonamide class of antibiotics. Unfortunately, bacterial resistance and adverse effects have limited the clinical utility of the sulfonamide antibiotics. Although six bacterial crystal structures are available, the flexible loop regions that enclose pABA during binding and contain key sulfonamide resistance sites have yet to be visualized in their functional conformation. To gain a new understanding of the structural basis of sulfonamide resistance, the molecular mechanism of DHPS action, and to generate a screening structure for high-throughput virtual screening, molecular dynamics simulations were applied to model the conformations of the unresolved loops in the active site. Several series of molecular dynamics simulations were designed and performed utilizing enzyme substrates and inhibitors, a transition state analog, and a pterin-sulfamethoxazole adduct. The positions of key mutation sites conserved across several bacterial species were closely monitored during these analyses. These residues were shown to interact closely with the sulfonamide binding site. The simulations helped us gain new understanding of the positions of the flexible loops during inhibitor binding that has allowed the development of a DHPS structural model that could be used for high-through put virtual screening (HTVS). Additionally, insights gained on the location and possible function of key mutation sites on the flexible loops will facilitate the design of new, potent inhibitors of DHPS that can bypass resistance mutations that render sulfonamides inactive. Prior to performing high-throughput virtual screening, the docking and scoring functions to be used were validated using established techniques against the B. anthracis DHPS target. In this validation study, five commonly used docking programs, FlexX, Surflex, Glide, GOLD, and DOCK, as well as nine scoring functions, were evaluated for their utility in virtual screening against the novel pterin binding site. Their performance in ligand docking and virtual screening against this target was examined by their ability to reproduce a known inhibitor conformation and to correctly detect known active compounds seeded into three separate decoy sets. Enrichment was demonstrated by calculated enrichment factors at 1% and Receiver Operating Characteristic (ROC) curves. The effectiveness of post-docking relaxation prior to rescoring and consensus scoring were also evaluated. Of the docking and scoring functions evaluated, Surflex with SurflexScore and Glide with GlideScore performed best overall for virtual screening against the DHPS target. The next phase of the DHPS structure-based drug design project involved high-throughput virtual screening against the DHPS structural model previously developed and docking methodology validated against this target. Two general virtual screening methods were employed. First, large, virtual libraries were pre-filtered by 3D pharmacophore and modified Rule-of-Three fragment constraints. Nearly 5 million compounds from the ZINC databases were screened generating 3,104 unique, fragment-like hits that were subsequently docked and ranked by score. Second, fragment docking without pharmacophore filtering was performed on almost 285,000 fragment-like compounds obtained from databases of commercial vendors. Hits from both virtual screens with high predicted affinity for the pterin binding pocket, as determined by docking score, were selected for in vitro testing. Activity and structure-activity relationship of the active fragment compounds have been developed. Several compounds with micromolar activity were identified and taken to crystallographic trials. Finally, in our ligand-based research into M. tuberculosis active agents, a series of nitrofuranylamide and related aromatic compounds displaying potent activity was investigated utilizing 3-Dimensional Quantitative Structure-Activity Relationship (3D-QSAR) techniques. Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) methods were used to produce 3D-QSAR models that correlated the Minimum Inhibitory Concentration (MIC) values against M. tuberculosis with the molecular structures of the active compounds. A training set of 95 active compounds was used to develop the models, which were then evaluated by a series of internal and external cross-validation techniques. A test set of 15 compounds was used for the external validation. Different alignment and ionization rules were investigated as well as the effect of global molecular descriptors including lipophilicity (cLogP, LogD), Polar Surface Area (PSA), and steric bulk (CMR), on model predictivity. Models with greater than 70% predictive ability, as determined by external validation and high internal validity (cross validated r2 \u3e .5) were developed. Incorporation of lipophilicity descriptors into the models had negligible effects on model predictivity. The models developed will be used to predict the activity of proposed new structures and advance the development of next generation nitrofuranyl and related nitroaromatic anti-tuberculosis agents

Scoring functions and enrichment: a case study on Hsp90.

BACKGROUND: The need for fast and accurate scoring functions has been driven by the increased use of in silico virtual screening twinned with high-throughput screening as a method to rapidly identify potential candidates in the early stages of drug development. We examine the ability of some the most common scoring functions (GOLD, ChemScore, DOCK, PMF, BLEEP and Consensus) to discriminate correctly and efficiently between active and non-active compounds among a library of approximately 3,600 diverse decoy compounds in a virtual screening experiment against heat shock protein 90 (Hsp90). RESULTS: Firstly, we investigated two ranking methodologies, GOLDrank and BestScorerank. GOLDrank is based on ranks generated using GOLD. The various scoring functions, GOLD, ChemScore, DOCK, PMF, BLEEP and Consensus, are applied to the pose ranked number one by GOLD for that ligand. BestScorerank uses multiple poses for each ligand and independently chooses the best ranked pose of the ligand according to each different scoring function. Secondly, we considered the effect of introducing the Thr184 hydrogen bond tether to guide the docking process towards a particular solution, and its effect on enrichment. Thirdly, we considered normalisation to account for the known bias of scoring functions to select larger molecules. All the scoring functions gave fairly similar enrichments, with the exception of PMF which was consistently the poorest performer. In most cases, GOLD was marginally the best performing individual function; the Consensus score usually performed similarly to the best single scoring function. Our best results were obtained using the Thr184 tether in combination with the BestScorerank protocol and normalisation for molecular weight. For that particular combination, DOCK was the best individual function; DOCK recovered 90% of the actives in the top 10% of the ranked list; Consensus similarly recovered 89% of the actives in its top 10%. CONCLUSION: Overall, we demonstrate the validity of virtual screening as a method for identifying new leads from a pool of ligands with similar physicochemical properties and we believe that the outcome of this study provides useful insight into the setting up of a suitable docking and scoring protocol, resulting in enrichment of 'target active' compounds

Low molecular weight compounds from mushrooms as potential Bcl-2 inhibitors: docking and virtual screening studies

Mestrado com dupla diplomação com o Institut Superieur de Biotechnologie de MonastirMushrooms have the ability to promote apoptosis in tumor cell lines, but the mechanism of action is not quite well understood. Inhibition of the interaction between Bcl-2 and pro-apoptotic proteins could be an important step that leads to apoptosis. Therefore, the discovery of compounds with the ability to inhibit Bcl-2 is an ongoing research topic in drug discovery. In this study, we started by analyzing Bcl-2 experimental structures that are currently available in Protein Data Bank database. After analysis of the more relevant Bcl-2 structures, 4 were finally selected. An analysis of the best docking methodology was then performed using a cross-docking and re-docking approach while testing 2 docking softwares: AutoDock 4 and AutoDock Vina. Autodock4 provided the best docking results and was selected to perform a virtual screening study applied to a dataset of 40 Low Molecular Weight (LMW) compounds present in mushrooms, using the selected Bcl-2 structures as target. Results suggest that steroid are the more promising family, among the analyzed compounds, and may have the ability to interact with Bcl-2 and this way promoting tumor apoptosis. The steroids that presented lowest estimated binding energy (ΔG) were: Ganodermanondiol, Cerevisterol, Ganoderic Acid X and Lucidenic Lactone; with estimated ΔG values between -8,45 and -8,23 Kcal/mol. A detailed analysis of the docked conformation of these 4 top ranked LMW compounds was also performed and illustrates a plausible interaction between the 4 top raked steroids and Bcl-2, thus substantiating the accuracy of the predicted docked poses. Therefore, tumoral apoptosis promoted by mushroom might be related to Bcl-2 inhibition mediated by steroid family of compounds.Os cogumelos apresentam a capacidade de promover a apoptose em linhas células tumorais, No entanto o seu mecanismo de ação não é completamente conhecido. A inibição da interação entre Bcl-2 e proteínas pro-apoptóticas pode ser um passo importante na iniciação do processo de apoptose tumoral. Por essa razão, a descoberta de compostos que inibam a proteína Bcl-2 é uma área importante na descoberta de novos fármacos antitumorais. Neste estudo, começou-se por analisar as estruturas experimentais de Bcl-2 atualmente presentes na base de estruturas Protein Data Bank. Após análise das estruturas de Bcl-2 mais relevantes, 4 foram escolhidas. Um estudo de “cross-docking” e “re-docking” foi então realizado para escolher a metodologia de “docking” mais adequada. Testaram-se 2 softwares, o AutoDock 4 e o AutoDock Vina, e verificou-se que o AutoDock 4 apresentava melhores resultados, tendo sido o selecionado para realizar os ensaios de “screening” virtual dos 40 compostos de baixo peso molecular presentes em cogumelos, utilizando as 4 estruturas selecionadas. Os resultados obtidos sugerem que os esteroides são a família de compostos mais prometedores de entre as famílias de compostos estudadas. Os esteroides que apresentaram valores de energia de ligação (ΔG) mais baixos foram: Ganodermanondiol, Cerevisterol, Ácido Ganoderico X and Lactona Lucidénica, com valores de ΔG estimado entre -8,45 e -8,23 Kcal/mol. Uma análise detalhada da conformação de ligação foi também realizada dos 4 melhores compostos de baixo peso molecular melhor classificados. Esta análise demonstra um modo de interação plausível entre os compostos e a estrutura da Bcl-2, consubstanciando a eficácia dos resultados obtidos por “docking”. Conclui-se que o processo inibição de apoptose tumoral observada em cogumelos pode estar relacionado com a inibição da Bcl-2 por esteroides presentes nos cogumelos

Ligand-based virtual screening using binary kernel discrimination

This paper discusses the use of a machine-learning technique called binary kernel discrimination (BKD) for virtual screening in drug- and pesticide-discovery programmes. BKD is compared with several other ligand-based tools for virtual screening in databases of 2D structures represented by fragment bit-strings, and is shown to provide an effective, and reasonably efficient, way of prioritising compounds for biological screening

In silico assessment of potential druggable pockets on the surface of α1-Antitrypsin conformers

The search for druggable pockets on the surface of a protein is often performed on a single conformer, treated as a rigid body. Transient druggable pockets may be missed in this approach. Here, we describe a methodology for systematic in silico analysis of surface clefts across multiple conformers of the metastable protein α1-antitrypsin (A1AT). Pathological mutations disturb the conformational landscape of A1AT, triggering polymerisation that leads to emphysema and hepatic cirrhosis. Computational screens for small molecule inhibitors of polymerisation have generally focused on one major druggable site visible in all crystal structures of native A1AT. In an alternative approach, we scan all surface clefts observed in crystal structures of A1AT and in 100 computationally produced conformers, mimicking the native solution ensemble. We assess the persistence, variability and druggability of these pockets. Finally, we employ molecular docking using publicly available libraries of small molecules to explore scaffold preferences for each site. Our approach identifies a number of novel target sites for drug design. In particular one transient site shows favourable characteristics for druggability due to high enclosure and hydrophobicity. Hits against this and other druggable sites achieve docking scores corresponding to a Kd in the µM–nM range, comparing favourably with a recently identified promising lead. Preliminary ThermoFluor studies support the docking predictions. In conclusion, our strategy shows considerable promise compared with the conventional single pocket/single conformer approach to in silico screening. Our best-scoring ligands warrant further experimental investigation