Hotspots API: A Python Package for the Detection of Small Molecule Binding Hotspots and Application to Structure-Based Drug Design.

Methods that survey protein surfaces for binding hotspots can help to evaluate target tractability and guide exploration of potential ligand binding regions. Fragment Hotspot Maps builds upon interaction data mined from the CSD (Cambridge Structural Database) and exploits the idea of identifying hotspots using small chemical fragments, which is now widely used to design new drug leads. Prior to this publication, Fragment Hotspot Maps was only publicly available through a web application. To increase the accessibility of this algorithm we present the Hotspots API (application programming interface), a toolkit that offers programmatic access to the core Fragment Hotspot Maps algorithm, thereby facilitating the interpretation and application of the analysis. To demonstrate the package's utility, we present a workflow which automatically derives protein hydrogen-bond constraints for molecular docking with GOLD. The Hotspots API is available from https://github.com/prcurran/hotspots under the MIT license and is dependent upon the commercial CSD Python API.

Polyphony: superposition independent methods for ensemble-based drug discovery.

BACKGROUND: Structure-based drug design is an iterative process, following cycles of structural biology, computer-aided design, synthetic chemistry and bioassay. In favorable circumstances, this process can lead to the structures of hundreds of protein-ligand crystal structures. In addition, molecular dynamics simulations are increasingly being used to further explore the conformational landscape of these complexes. Currently, methods capable of the analysis of ensembles of crystal structures and MD trajectories are limited and usually rely upon least squares superposition of coordinates. RESULTS: Novel methodologies are described for the analysis of multiple structures of a protein. Statistical approaches that rely upon residue equivalence, but not superposition, are developed. Tasks that can be performed include the identification of hinge regions, allosteric conformational changes and transient binding sites. The approaches are tested on crystal structures of CDK2 and other CMGC protein kinases and a simulation of p38α. Known interaction - conformational change relationships are highlighted but also new ones are revealed. A transient but druggable allosteric pocket in CDK2 is predicted to occur under the CMGC insert. Furthermore, an evolutionarily-conserved conformational link from the location of this pocket, via the αEF-αF loop, to phosphorylation sites on the activation loop is discovered. CONCLUSIONS: New methodologies are described and validated for the superimposition independent conformational analysis of large collections of structures or simulation snapshots of the same protein. The methodologies are encoded in a Python package called Polyphony, which is released as open source to accompany this paper [http://wrpitt.bitbucket.org/polyphony/].

Annular tautomerism: experimental observations and quantum mechanics calculations.

The use of MP2 level quantum mechanical (QM) calculations on isolated heteroaromatic ring systems for the prediction of the tautomeric propensities of whole molecules in a crystalline environment was examined. A Polarisable Continuum Model was used in the calculations to account for environment effects on the tautomeric relative stabilities. The calculated relative energies of tautomers were compared to relative abundances within the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB). The work was focussed on 84 annular tautomeric forms of 34 common ring systems. Good agreement was found between the calculations and the experimental data even if the quantity of these data was limited in many cases. The QM results were compared to those produced by much faster semiempirical calculations. In a search for other sources of the useful experimental data, the relative numbers of known compounds in which prototropic positions were often substituted by heavy atoms were also analysed. A scheme which groups all annular tautomeric transformations into 10 classes was developed. The scheme was designed to encompass a comprehensive set of known and theoretically possible tautomeric ring systems generated as part of a previous study. General trends across analogous ring systems were detected as a result. The calculations and statistics collected on crystallographic data as well as the general trends observed should be useful for the better modelling of annular tautomerism in the applications such as computer-aided drug design, small molecule crystal structure prediction, the naming of compounds and the interpretation of protein-small molecule crystal structures.

Metabolism by Aldehyde Oxidase: Drug Design and Complementary Approaches to Challenges in Drug Discovery.

Aldehyde oxidase (AO) catalyzes oxidations of azaheterocycles and aldehydes, amide hydrolysis, and diverse reductions. AO substrates are rare among marketed drugs, and many candidates failed due to poor pharmacokinetics, interspecies differences, and adverse effects. As most issues arise from complex and poorly understood AO biology, an effective solution is to stop or decrease AO metabolism. This perspective focuses on rational drug design approaches to modulate AO-mediated metabolism in drug discovery. AO biological aspects are also covered, as they are complementary to chemical design and important when selecting the experimental system for risk assessment. The authors' recommendation is an early consideration of AO-mediated metabolism supported by computational and in vitro experimental methods but not an automatic avoidance of AO structural flags, many of which are versatile and valuable building blocks. Preferably, consideration of AO-mediated metabolism should be part of the multiparametric drug optimization process, with the goal to improve overall drug-like properties.

IRAK-4 inhibitors. Part II: a structure-based assessment of imidazo[1,2-a]pyridine binding.

A potent IRAK-4 inhibitor was identified through routine project cross screening. The binding mode was inferred using a combination of in silico docking into an IRAK-4 homology model, surrogate crystal structure analysis and chemical analogue SAR.

IRAK-4 inhibitors. Part 1: a series of amides.

The synthesis and profile of a series of amides are described. Some of these compounds were potent IRAK-4 inhibitors and two examples were evaluated in vivo.

IRAK-4 inhibitors. Part III: a series of imidazo[1,2-a]pyridines.

Following the identification of a potent IRAK inhibitor through routine project cross screening, a novel class of IRAK-4 inhibitor was established. The SAR of imidazo[1,2-a]pyridino-pyridines and benzimidazolo-pyridines was explored.

Identification of novel fragment compounds targeted against the pY pocket of v-Src SH2 by computational and NMR screening and thermodynamic evaluation.

Discovery of small molecule inhibitors of protein-protein interactions is a major challenge to pharmaceutical development. Fragment-based approaches have begun to be widely adopted as an effective way of exploring chemical space on a protein surface with reduced library size. On completion of a fragment screen, the subsequent selection of appropriate "hit" molecules for development is a key decision point. Thermodynamic parameters can be used in this decision process. In this work, a fragment identification protocol based on a virtual fragment analysis and selection followed by 19F NMR screening was directed at the phosphotyrosine binding site of the Src SH2 domain. Three new ligands were identified. Isothermal titration calorimetry was used to provide thermodynamic parameters for the physiologically relevant ligand and the selected fragments. One of these fragments possesses a highly favorable enthalpic contribution to complex formation compared to other fragments and to the physiologically relevant ligand suggesting that it would make a good candidate for compound development.

Structure-based virtual screening for novel ligands.

Computational searches for novel ligands for a given protein binding site have become ubiquitous in the pharmaceutical industry, and are potentially equally useful in helping identify small-molecule tools for biology. Here we describe the steps needed to carry out a high-throughput docking (HTD) or three-dimensional (3D) pharmacophore virtual screen starting with a model of the target protein's structure. The advice given is, in most cases, software independent but some tips are provided which apply only to certain popular programs. Useful work can be carried out using free programs on a modest workstation. Of course, any resultant "hits" remain in the virtual world until they are experimentally tested.

Heteroaromatic rings of the future.

Small aromatic ring systems are of central importance in the development of novel synthetic protein ligands. Here we generate a complete list of 24,847 such ring systems. We call this list and associated annotations VEHICLe, which stands for virtual exploratory heterocyclic library. Searches of literature and compound databases, using this list as substructure queries, identified only 1701 as synthesized. Using a carefully validated machine learning approach, we were able to estimate that the number of unpublished, but synthetically tractable, VEHICLe rings could be over 3000. However, analysis also shows that the rate of publication of novel examples to be as low as 5-10 per year. With this work, we aim to provide fresh stimulus to creative organic chemists by highlighting a small set of apparently simple ring systems that are predicted to be tractable but are, to the best of our knowledge, unconquered.

On the origins of enzyme inhibitor selectivity and promiscuity: a case study of protein kinase binding to staurosporine.

Relationships between ligand binding and the shapes of the binding sites in families of homologous enzymes are investigated by comparing matrices of distances between key binding site atoms. Multiple linear regression is used to help identify key distances that influence ligand binding affinity. In order to illustrate the utility of this generic approach, we study protein kinase binding sites for ATP and the promiscuous competitive inhibitor, staurosporine. We show that the size of the gatekeeper residue and the closure between the first glycine of the GXGXXG motif and the aspartate of the DFG loop act together to promote tight binding. Our web-based tool, 'mapping analogous hetero-atoms onto residue interactions' (MAHORI), indicates that the greater the number of hydrogen bonds made by the kinase around the methylamine group of staurosporine, the tighter the binding. The conservation of surrounding atoms identified using our novel grid-based method clearly demonstrates that the most structurally conserved part of the binding site for staurosporine is the main chain of the hinge region. The critical role of interactions that are not dependent on side-chain identities is consistent with the promiscuous nature of this inhibitor.

The thermodynamics of protein-ligand interaction and solvation: insights for ligand design.

Isothermal titration calorimetry is able to provide accurate information on the thermodynamic contributions of enthalpy and entropy changes to free energies of binding. The Structure/Calorimetry of Reported Protein Interactions Online database of published isothermal titration calorimetry studies and structural information on the interactions between proteins and small-molecule ligands is used here to reveal general thermodynamic properties of protein-ligand interactions and to investigate correlations with changes in solvation. The overwhelming majority of interactions are found to be enthalpically favoured. Synthetic inhibitors and biological ligands form two distinct subpopulations in the data, with the former having greater average affinity due to more favourable entropy changes on binding. The greatest correlation is found between the binding free energy and apolar surface burial upon complex formation. However, the free-energy contribution per unit area buried is only 30-50% of that expected from earlier studies of transfer free energies of small molecules. A simple probability-based estimator for the maximal affinity of a binding site in terms of its apolar surface area is proposed. Polar surface area burial also contributes substantially to affinity but is difficult to express in terms of unit area due to the small variation in the amount of polar surface buried and a tendency for cancellation of its enthalpic and entropic contributions. Conventionally, the contribution of apolar desolvation to affinity is attributed to gain of entropy due to solvent release. Although data presented here are supportive of this notion, because the correlation of entropy change with apolar surface burial is relatively weak, it cannot, on present evidence, be confidently considered to be correct. Further, thermodynamic changes arising from small differences between ligands binding to individual proteins are relatively large and, in general, uncorrelated with changes in solvation, suggesting that trends identified across widely differing proteins are of limited use in explaining or predicting the effects of ligand modifications.

Synthesis and SAR of aminopyrimidines as novel c-Jun N-terminal kinase (JNK) inhibitors.

The development of a series of novel aminopyrimidines as inhibitors of c-Jun N-terminal kinases is described. The synthesis, in vitro inhibitory values for JNK1, JNK2 and CDK2, and the in vitro inhibitory value for a c-Jun cellular assay are discussed.

Low molecular weight indole fragments as IMPDH inhibitors.

The study of non-oxazole containing indole fragments as inhibitors of inosine monophosphate dehydrogenase (IMPDH) is described. The synthesis and in vitro inhibitory values for IMPDH II are discussed.

Novel indole inhibitors of IMPDH from fragments: synthesis and initial structure-activity relationships.

The elaboration of previously reported indole fragments as inhibitors of inosine monophosphate dehydrogenase (IMPDH) is described. The synthesis, in vitro inhibitory values for IMPDH II, PBMC proliferation and physicochemical properties are discussed.

Novel 7-methoxy-6-oxazol-5-yl-2,3-dihydro-1H-quinazolin-4-ones as IMPDH inhibitors.

The synthesis and biological activity of a novel series of 7-methoxy-6-oxazol-5-yl-2,3-dihydro-1H-quinazolin-4-ones are described. Some of these compounds were found to be potent inhibitors of inosine 5'-monophosphate dehydrogenase type II (IMPDH II).

Quinazolinethiones and quinazolinediones, novel inhibitors of inosine monophosphate dehydrogenase: synthesis and initial structure-activity relationships.

The development of a series of novel quinazolinethiones and quinazolinediones as inhibitors of inosine monophosphate dehydrogenase (IMPDH) is described. The synthesis, in vitro inhibitory values for IMPDH II and in vitro inhibitory value for PBMC proliferation are discussed.