Steven Shave

Ligand discovery and virtual screening using the program LIDAEUS

This paper discusses advances in docking and scoring approaches with examples from the high‐throughput virtual screening program LIDAEUS. We describe the discovery of small molecule inhibitors for the immunophilin CypA, the cyclin‐dependent kinase CDK2 and the cyclapolin series of potent Polo‐like kinase inhibitors. These results are discussed in the context of advances in massively parallel computing and in the development of annotated databases.

Inhibition of the ERCC1-XPF structure-specific endonuclease to overcome cancer chemoresistance

ERCC1-XPF is a structure-specific endonuclease that is required for the repair of DNA lesions, generated by the widely used platinum-containing cancer chemotherapeutics such as cisplatin, through the Nucleotide Excision Repair and Interstrand Crosslink Repair pathways. Based on mouse xenograft experiments, where ERCC1-deficient melanomas were cured by cisplatin therapy, we proposed that inhibition of ERCC1-XPF could enhance the effectiveness of platinum-based chemotherapy. Here we report the identification and properties of inhibitors against two key targets on ERCC1-XPF. By targeting the ERCC1-XPF interaction domain we proposed that inhibition would disrupt the ERCC1-XPF heterodimer resulting in destabilisation of both proteins. Using in silico screening, we identified an inhibitor that bound to ERCC1-XPF in a biophysical assay, reduced the level of ERCC1-XPF complexes in ovarian cancer cells, inhibited Nucleotide Excision Repair and sensitised melanoma cells to cisplatin. We also utilised high throughput and in silico screening to identify the first reported inhibitors of the other key target, the XPF endonuclease domain. We demonstrate that two of these compounds display specificity in vitro for ERCC1-XPF over two other endonucleases, bind to ERCC1-XPF, inhibit Nucleotide Excision Repair in two independent assays and specifically sensitise Nucleotide Excision Repair-proficient, but not Nucleotide Excision Repair-deficient human and mouse cells to cisplatin.

EDULISS: a small-molecule database with data-mining and pharmacophore searching capabilities

We present the relational database EDULISS (EDinburgh University Ligand Selection System), which stores structural, physicochemical and pharmacophoric properties of small molecules. The database comprises a collection of over 4 million commercially available compounds from 28 different suppliers. A user-friendly web-based interface for EDULISS (available at http://eduliss.bch.ed.ac.uk/) has been established providing a number of data-mining possibilities. For each compound a single 3D conformer is stored along with over 1600 calculated descriptor values (molecular properties). A very efficient method for unique compound recognition, especially for a large scale database, is demonstrated by making use of small subgroups of the descriptors. Many of the shape and distance descriptors are held as pre-calculated bit strings permitting fast and efficient similarity and pharmacophore searches which can be used to identify families of related compounds for biological testing. Two ligand searching applications are given to demonstrate how EDULISS can be used to extract families of molecules with selected structural and biophysical features.

Discovery of a new class of inhibitors for the protein arginine deiminase type 4 (PAD4) by structure-based virtual screening

BackgroundRheumatoid arthritis (RA) is an autoimmune disease with unknown etiology. Anticitrullinated protein autoantibody has been documented as a highly specific autoantibody associated with RA. Protein arginine deiminase type 4 (PAD4) is the enzyme responsible for catalyzing the conversion of peptidylarginine into peptidylcitrulline. PAD4 is a new therapeutic target for RA treatment. In order to search for inhibitors of PAD4, structure-based virtual screening was performed using LIDAEUS (Ligand discovery at Edinburgh university). Potential inhibitors were screened experimentally by inhibition assays.ResultsTwenty two of the top-ranked water-soluble compounds were selected for inhibitory screening against PAD4. Three compounds showed significant inhibition of PAD4 and their IC50 values were investigated. The structures of the three compounds show no resemblance with previously discovered PAD4 inhibitors, nor with existing drugs for RA treatment.ConclusionThree compounds were discovered as potential inhibitors of PAD4 by virtual screening. The compounds are commercially available and can be used as scaffolds to design more potent inhibitors against PAD4.

11β-Hydroxysteroid dehydrogenase type 1 contributes to the balance between 7-keto- and 7-hydroxy-oxysterols in vivo

Graphical abstract

High-Throughput Chemical Screening for Antivirulence Developmental Phenotypes in Trypanosoma brucei

ABSTRACT In the bloodstream of mammalian hosts, the sleeping sickness parasite, Trypanosoma brucei, exists as a proliferative slender form or a nonproliferative, transmissible, stumpy form. The transition between these developmental forms is controlled by a density-dependent mechanism that is important for the parasites infection dynamics, immune evasion via ordered antigenic variation, and disease transmissibility. However, stumpy formation has been lost in most laboratory-adapted trypanosome lines, generating monomorphic parasites that proliferate uncontrolled as slender forms in vitro and in vivo. Nonetheless, these forms are readily amenable to cell culture and high-throughput screening for trypanocidal lead compounds. Here, we have developed and exploited a high-throughput screen for developmental phenotypes using a transgenic monomorphic cell line expressing a reporter under the regulation of gene control signals from the stumpy-specific molecule PAD1. Using a whole-cell fluorescence-based assay to screen over 6,000 small molecules from a kinase-focused compound library, small molecules able to activate stumpy-specific gene expression and proliferation arrest were assayed in a rapid assay format. Independent follow-up validation identified one hit able to induce modest, yet specific, changes in mRNA expression indicative of a partial differentiation to stumpy forms in monomorphs. Further, in pleomorphs this compound induced a stumpy-like phenotype, entailing growth arrest, morphological changes, PAD1 expression, and enhanced differentiation to procyclic forms. This not only provides a potential tool compound for the further understanding of stumpy formation but also demonstrates the use of high-throughput screening in the identification of compounds able to induce specific phenotypes, such as differentiation, in African trypanosomes.

Identification and X-ray co-crystal structure of a small-molecule activator of LFA-1-ICAM-1 binding.

Stabilization of protein–protein interactions by small molecules is a concept with few examples reported to date. Herein we describe the identification and X-ray co-crystal structure determination of IBE-667, an ICAM-1 binding enhancer for LFA-1. IBE-667 was designed based on the SAR information obtained from an on-bead screen of tagged one-bead one-compound combinatorial libraries by confocal nanoscanning and bead picking (CONA). Cellular assays demonstrate the activity of IBE-667 in promoting the binding of LFA-1 on activated immune cells to ICAM-1.

UFSRAT: Ultra-fast Shape Recognition with Atom Types--the discovery of novel bioactive small molecular scaffolds for FKBP12 and 11βHSD1.

Motivation Using molecular similarity to discover bioactive small molecules with novel chemical scaffolds can be computationally demanding. We describe Ultra-fast Shape Recognition with Atom Types (UFSRAT), an efficient algorithm that considers both the 3D distribution (shape) and electrostatics of atoms to score and retrieve molecules capable of making similar interactions to those of the supplied query. Results Computational optimization and pre-calculation of molecular descriptors enables a query molecule to be run against a database containing 3.8 million molecules and results returned in under 10 seconds on modest hardware. UFSRAT has been used in pipelines to identify bioactive molecules for two clinically relevant drug targets; FK506-Binding Protein 12 and 11β-hydroxysteroid dehydrogenase type 1. In the case of FK506-Binding Protein 12, UFSRAT was used as the first step in a structure-based virtual screening pipeline, yielding many actives, of which the most active shows a KD, app of 281 µM and contains a substructure present in the query compound. Success was also achieved running solely the UFSRAT technique to identify new actives for 11β-hydroxysteroid dehydrogenase type 1, for which the most active displays an IC50 of 67 nM in a cell based assay and contains a substructure radically different to the query. This demonstrates the valuable ability of the UFSRAT algorithm to perform scaffold hops. Availability and Implementation A web-based implementation of the algorithm is freely available at http://opus.bch.ed.ac.uk/ufsrat/.

Ligand discovery on massively parallel systems

Virtual screening is an approach for identifying promising leads for drugs and is used in the pharmaceutical industry. We present the parallelization of LIDAEUS (LIgand Discovery At Edinburgh UniverSity), creating a massively parallel high-throughput virtual-screening code. This program is being used to predict the binding modes involved in the docking of small ligands to proteins. Parallelization efforts have focused on achieving maximum parallel efficiency and developing a memory-efficient parallel sorting routine. Using an IBM Blue Gene/L™ supercomputer, runtimes have been reduced from 8 days on a modest seven-node cluster to 62 minutes on 1,024 processors using a standard dataset of 1.67 million small molecules and FKBP12, a protein target of interest in immunosuppressive therapies. Using more-complex datasets, the code scales upward to make use for the full processor set of 2,048. The code has been successfully used for the task of gathering data on approximately 1.67 million small molecules binding to approximately 400 high-quality crystallographically determined ligand-bound protein structures, generating data on more than 646 million protein-ligand complexes. A number of novel ligands have already been discovered and validated experimentally.

PuLSE: Quality control and quantification of peptide sequences explored by phage display libraries

The design of highly diverse phage display libraries is based on assumption that DNA bases are incorporated at similar rates within the randomized sequence. As library complexity increases and expected copy numbers of unique sequences decrease, the exploration of library space becomes sparser and the presence of truly random sequences becomes critical. We present the program PuLSE (Phage Library Sequence Evaluation) as a tool for assessing randomness and therefore diversity of phage display libraries. PuLSE runs on a collection of sequence reads in the fastq file format and generates tables profiling the library in terms of unique DNA sequence counts and positions, translated peptide sequences, and normalized ‘expected’ occurrences from base to residue codon frequencies. The output allows at-a-glance quantitative quality control of a phage library in terms of sequence coverage both at the DNA base and translated protein residue level, which has been missing from toolsets and literature. The open source program PuLSE is available in two formats, a C++ source code package for compilation and integration into existing bioinformatics pipelines and precompiled binaries for ease of use.