William R. Pitt

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.

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.

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.

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.

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. Proteins 2007.

AQUARIUS2: knowledge-based modeling of solvent sites around proteins

The program AQUARIUS2 calculates probable positions for water molecules within the first hydration shell of any protein for which atomic coordinates are known. Like its predecessor, AQUARIUS, it uses a knowledge of water molecules sites from crystallographically determined protein structures. Energy calculations are not employed. It differs substantially from the original program in that a 3‐D probability map (for solvent sites) is generated around the surface of the protein instead of the previously used discrete points. The accuracy of the program has been gauged by comparison with experimentally derived water molecule positions for proteins not used in the knowledge base of the program. It has also been tested by combining the probability density maps with crystallographically determined electron density maps for the protein porphobilinogen deaminase. This procedure filters the most likely solvent electron density peaks from the background noise and has been used in the determination of the solvent structure around the protein nerve growth factor.

New methods for the analysis of the protein-solvent interface

Abstract The protein-solvent interface is complex and may include solvent channels and cavities as well as the normal surface water molecules. We describe several algorithms for investigating the intra- and inter-molecular interactions of proteins in general but with the aim of developing methods to accurately and definitively characterise the interactions of water and other small ligands with proteins. Specifically, we present the methods which underlie three programs (AQUARIUS2, HOLE and PRO_ACT) which can be used to to look at different aspects of these interactions.