Small and colorful stones make beautiful mosaics: fragment-based chemogenomics

An Automated, Open-Source Workflow for the Generation of (3D) Fragment Libraries

The recent success of fragment-based drug discovery (FBDD) is inextricably linked to adequate library design. To guide the design of our fragment libraries, we have constructed an automated workflow in the open-source KNIME software. The workflow considers chemical diversity and novelty of the fragments, and can also take into account the three-dimensional (3D) character. This design tool can be used to create large and diverse libraries but also to select a small number of representative compounds as a focused set of unique screening compounds to enrich existing fragment libraries. To illustrate the procedures, the design and synthesis of a 10-membered focused library is reported based on the cyclopropane scaffold, which is underrepresented in our existing fragment screening library. Analysis of the focused compound set indicates significant shape diversity and a favorable overall physicochemical profile. By virtue of its modular setup, the workflow can be readily adjusted to design libraries that focus on properties other than 3D shape.

Fragment screening yields a small-molecule stabilizer of 14-3-3 dimers that modulates client protein interactions

The development of protein‐protein interaction (PPI) inhibitors has been a successful strategy in drug development. However, the identification of PPI stabilizers has proven much more challenging. Here we report a fragment‐based drug screening approach using the regulatory hub‐protein 14‐3‐3 as a platform for identifying PPI stabilizers. A homogenous time‐resolved FRET assay was used to monitor stabilization of 14‐3‐3/peptide binding using the known interaction partner estrogen receptor alpha. Screening of an in‐house fragment library identified fragment 2 (VUF15640) as a putative PPI stabilizer capable of cooperatively stabilizing 14‐3‐3 PPIs in a cooperative fashion with Fusicoccin‐A. Mechanistically, fragment 2 appears to enhance 14‐3‐3 dimerization leading to increased client‐protein binding. Functionally, fragment 2 enhanced potency of 14‐3‐3 in a cell‐free system inhibiting the enzyme activity of the nitrate reductase. In conclusion, we identified a general PPI stabilizer targeting 14‐3‐3, which could be used as a tool compound for investigating 14‐3‐3 client protein interactions.

Bromo-Cyclobutenaminones as New Covalent UDP-N-Acetylglucosamine Enolpyruvyl Transferase (MurA) Inhibitors

Drug discovery programs against the antibacterial target UDP-N-acetylglucosamine enolpyruvyl transferase (MurA) have already resulted in covalent inhibitors having small three- and five-membered heterocyclic rings. In the current study, the reactivity of four-membered rings was carefully modulated to obtain a novel family of covalent MurA inhibitors. Screening a small library of cyclobutenone derivatives led to the identification of bromo-cyclobutenaminones as new electrophilic warheads. The electrophilic reactivity and cysteine specificity have been determined in a glutathione (GSH) and an oligopeptide assay, respectively. Investigating the structure-activity relationship for MurA suggests a crucial role for the bromine atom in the ligand. In addition, MS/MS experiments have proven the covalent labelling of MurA at Cys115 and the observed loss of the bromine atom suggests a net nucleophilic substitution as the covalent reaction. This new set of compounds might be considered as a viable chemical starting point for the discovery of new MurA inhibitors.

Structure-based exploration and pharmacological evaluation of N-substituted piperidin-4-yl-methanamine CXCR4 chemokine receptor antagonists

Using the available structural information of the chemokine receptor CXCR4, we present hit finding and hit exploration studies that make use of virtual fragment screening, design, synthesis and structure-activity relationship (SAR) studies. Fragment 2 was identified as virtual screening hit and used as a starting point for the exploration of 31 N-substituted piperidin-4-yl-methanamine derivatives to investigate and improve the interactions with the CXCR4 binding site. Additionally, subtle structural ligand changes lead to distinct interactions with CXCR4 resulting in a full to partial displacement of CXCL12 binding and competitive and/or non-competitive antagonism. Three-dimensional quantitative structure-activity relationship (3D-QSAR) and binding model studies were used to identify important hydrophobic interactions that determine binding affinity and indicate key ligand-receptor interactions.

Chemical Diversity in the G Protein-Coupled Receptor Superfamily

G protein-coupled receptors (GPCRs) are the largest family of cell signaling transmembrane proteins that can be modulated by a plethora of chemical compounds. Systematic cheminformatics analysis of structurally and pharmacologically characterized GPCR ligands shows that cocrystallized GPCR ligands cover a significant part of chemical ligand space, despite their limited number. Many GPCR ligands and substructures interact with multiple receptors, providing a basis for polypharmacological ligand design. Experimentally determined GPCR structures represent a variety of binding sites and receptor-ligand interactions that can be translated to chemically similar ligands for which structural data are lacking. This integration of structural, pharmacological, and chemical information on GPCR-ligand interactions enables the extension of the structural GPCR-ligand interactome and the structure-based design of novel modulators of GPCR function.

Structural Analysis of Chemokine Receptor-Ligand Interactions

This review focuses on the construction and application of structural chemokine receptor models for the elucidation of molecular determinants of chemokine receptor modulation and the structure-based discovery and design of chemokine receptor ligands. A comparative analysis of ligand binding pockets in chemokine receptors is presented, including a detailed description of the CXCR4, CCR2, CCR5, CCR9, and US28 X-ray structures, and their implication for modeling molecular interactions of chemokine receptors with small-molecule ligands, peptide ligands, and large antibodies and chemokines. These studies demonstrate how the integration of new structural information on chemokine receptors with extensive structure–activity relationship and site-directed mutagenesis data facilitates the prediction of the structure of chemokine receptor–ligand complexes that have not been crystallized. Finally, a review of structure-based ligand discovery and design studies based on chemokine receptor crystal structures and homology models illustrates the possibilities and challenges to find novel ligands for chemokine receptors.

Molecular interaction fingerprint approaches for GPCR drug discovery

Protein-ligand interaction fingerprints (IFPs) are binary 1D representations of the 3D structure of protein-ligand complexes encoding the presence or absence of specific interactions between the binding pocket amino acids and the ligand. Various implementations of IFPs have been developed and successfully applied for post-processing molecular docking results for G Protein-Coupled Receptor (GPCR) ligand binding mode prediction and virtual ligand screening. Novel interaction fingerprint methods enable structural chemogenomics and polypharmacology predictions by complementing the increasing amount of GPCR structural data. Machine learning methods are increasingly used to derive relationships between bioactivity data and fingerprint descriptors of chemical and structural information of binding sites, ligands, and protein-ligand interactions. Factors that influence the application of IFPs include structure preparation, binding site definition, fingerprint similarity assessment, and data processing and these factors pose challenges as well possibilities to optimize interaction fingerprint methods for GPCR drug discovery.

Twenty years on: the impact of fragments on drug discovery

After 20 years of sometimes quiet growth, fragment-based drug discovery (FBDD) has become mainstream. More than 30 drug candidates derived from fragments have entered the clinic, with two approved and several more in advanced trials. FBDD has been widely applied in both academia and industry, as evidenced by the large number of papers from universities, non-profit research institutions, biotechnology companies and pharmaceutical companies. Moreover, FBDD draws on a diverse range of disciplines, from biochemistry and biophysics to computational and medicinal chemistry. As the promise of FBDD strategies becomes increasingly realized, now is an opportune time to draw lessons and point the way to the future. This Review briefly discusses how to design fragment libraries, how to select screening techniques and how to make the most of information gleaned from them. It also shows how concepts from FBDD have permeated and enhanced drug discovery efforts.

KLIFS: a structural kinase-ligand interaction database

Protein kinases play a crucial role in cell signaling and are important drug targets in several therapeutic areas. The KLIFS database contains detailed structural kinase-ligand interaction information derived from all (>2900) structures of catalytic domains of human and mouse protein kinases deposited in the Protein Data Bank in order to provide insights into the structural determinants of kinase-ligand binding and selectivity. The kinase structures have been processed in a consistent manner by systematically analyzing the structural features and molecular interaction fingerprints (IFPs) of a predefined set of 85 binding site residues with bound ligands. KLIFS has been completely rebuilt and extended (>65% more structures) since its first release as a data set, including: novel automated annotation methods for (i) the assessment of ligand-targeted subpockets and the analysis of (ii) DFG and (iii) αC-helix conformations; improved and automated protocols for (iv) the generation of sequence/structure alignments, (v) the curation of ligand atom and bond typing for accurate IFP analysis and (vi) weekly database updates. KLIFS is now accessible via a website (http://klifs.vu-compmedchem.nl) that provides a comprehensive visual presentation of different types of chemical, biological and structural chemogenomics data, and allows the user to easily access, compare, search and download the data.

A fluorescent approach for identifying P2X1 ligands

There are no commercially available, small, receptor-specific P2X1 ligands. There are several synthetic derivatives of the natural agonist ATP and some structurally-complex antagonists including compounds such as PPADS, NTP-ATP, suramin and its derivatives (e.g. NF279, NF449). NF449 is the most potent and selective ligand, but potencies of many others are not particularly high and they can also act at other P2X, P2Y and non-purinergic receptors. While there is clearly scope for further work on P2X1 receptor pharmacology, screening can be difficult owing to rapid receptor desensitisation. To reduce desensitisation substitutions can be made within the N-terminus of the P2X1 receptor, but these could also affect ligand properties. An alternative is the use of fluorescent voltage-sensitive dyes that respond to membrane potential changes resulting from channel opening. Here we utilised this approach in conjunction with fragment-based drug-discovery. Using a single concentration (300 μM) we identified 46 novel leads from a library of 1443 fragments (hit rate = 3.2%). These hits were independently validated by measuring concentration-dependence with the same voltage-sensitive dye, and by visualising the competition of hits with an Alexa-647-ATP fluorophore using confocal microscopy; confocal yielded k_on (1.142 × 10⁶ M⁻¹ s⁻¹) and k_off (0.136 s⁻¹) for Alexa-647-ATP (K_d = 119 nM). The identified hit fragments had promising structural diversity. In summary, the measurement of functional responses using voltage-sensitive dyes was flexible and cost-effective because labelled competitors were not needed, effects were independent of a specific binding site, and both agonist and antagonist actions were probed in a single assay. The method is widely applicable and could be applied to all P2X family members, as well as other voltage-gated and ligand-gated ion channels.

This article is part of the Special Issue entitled ‘Fluorescent Tools in Neuropharmacology’.

•

A novel fluorescence-based screening approach for identifying P2X1 receptor ligand candidates.

•

Fragment-based drug discovery applied to ligand-gated ion channels.

•

The use of confocal microscopy to determine the kinetics and affinity of Alexa-647-ATP binding to P2X1 receptors.

•

Alexa-647-ATP for imaging P2X1 receptors on live cells.

Combinatorial Consensus Scoring for Ligand-Based Virtual Fragment Screening: A Comparative Case Study for Serotonin 5-HT(3)A, Histamine H(1), and Histamine H(4) Receptors

In the current study we have evaluated the applicability of ligand-based virtual screening (LBVS) methods for the identification of small fragment-like biologically active molecules using different similarity descriptors and different consensus scoring approaches. For this purpose, we have evaluated the performance of 14 chemical similarity descriptors in retrospective virtual screening studies to discriminate fragment-like ligands of three membrane-bound receptors from fragments that are experimentally determined to have no affinity for these proteins (true inactives). We used a complete fragment affinity data set of experimentally determined ligands and inactives for two G protein-coupled receptors (GPCRs), the histamine H1 receptor (H1R) and the histamine H4 receptor (H4R), and one ligand-gated ion channel (LGIC), the serotonin receptor (5-HT3AR), to validate our retrospective virtual screening studies. We have exhaustively tested consensus scoring strategies that combine the results of multiple actives (group fusion) or combine different similarity descriptors (similarity fusion), and for the first time systematically evaluated different combinations of group fusion and similarity fusion approaches. Our studies show that for these three case study protein targets both consensus scoring approaches can increase virtual screening enrichments compared to single chemical similarity search methods. Our cheminformatics analyses recommend to use a combination of both group fusion and similarity fusion for prospective ligand-based virtual fragment screening.

Computational methods for fragment-based ligand design: growing and linking

Fragment-based drug design has proved itself as a powerful technique for increasing the sampling and diversity of chemical space and enabling the design of novel leads and compounds. Computational techniques for identifying fragments, binding sites and particularly for linking, growing, and evolving fragments play a significant role in the process. Information from ADME studies and clustering property information in the form of toxicophores and chemotypes can play a significant role in aiding the design of novel, selective fragments with good activity profiles.

Structure-based virtual screening for fragment-like ligands of the G protein-coupled histamine H4 receptor

Structure-based virtual screening using H₁R- and β₂R-based histamine H₄R homology models identified 9 fragments with an affinity ranging from 0.14 to 6.3 μm for H₄R.

Fragment-Based Screening in Tandem with Phenotypic Screening Provides Novel Antiparasitic Hits

Methods to discover biologically active small molecules include target-based and phenotypic screening approaches. One of the main difficulties in drug discovery is elucidating and exploiting the relationship between drug activity at the protein target and disease modification, a phenotypic endpoint. Fragment-based drug discovery is a target-based approach that typically involves the screening of a relatively small number of fragment-like (molecular weight <300) molecules that efficiently cover chemical space. Here, we report a fragment screening on TbrPDEB1, an essential cyclic nucleotide phosphodiesterase (PDE) from Trypanosoma brucei, and human PDE4D, an off-target, in a workflow in which fragment hits and a series of close analogs are subsequently screened for antiparasitic activity in a phenotypic panel. The phenotypic panel contained T. brucei, Trypanosoma cruzi, Leishmania infantum, and Plasmodium falciparum, the causative agents of human African trypanosomiasis (sleeping sickness), Chagas disease, leishmaniasis, and malaria, respectively, as well as MRC-5 human lung cells. This hybrid screening workflow has resulted in the discovery of various benzhydryl ethers with antiprotozoal activity and low toxicity, representing interesting starting points for further antiparasitic optimization.

A structural chemogenomics analysis of aminergic GPCRs: lessons for histamine receptor ligand design

BACKGROUND AND PURPOSE

Chemogenomics focuses on the discovery of new connections between chemical and biological space leading to the discovery of new protein targets and biologically active molecules. G-protein coupled receptors (GPCRs) are a particularly interesting protein family for chemogenomics studies because there is an overwhelming amount of ligand binding affinity data available. The increasing number of aminergic GPCR crystal structures now for the first time allows the integration of chemogenomics studies with high-resolution structural analyses of GPCR-ligand complexes.

EXPERIMENTAL APPROACH

In this study, we have combined ligand affinity data, receptor mutagenesis studies, and amino acid sequence analyses to high-resolution structural analyses of (hist)aminergic GPCR-ligand interactions. This integrated structural chemogenomics analysis is used to more accurately describe the molecular and structural determinants of ligand affinity and selectivity in different key binding regions of the crystallized aminergic GPCRs, and histamine receptors in particular.

KEY RESULTS

Our investigations highlight interesting correlations and differences between ligand similarity and ligand binding site similarity of different aminergic receptors. Apparent discrepancies can be explained by combining detailed analysis of crystallized or predicted protein-ligand binding modes, receptor mutation studies, and ligand structure-selectivity relationships that identify local differences in essential pharmacophore features in the ligand binding sites of different receptors.

CONCLUSIONS AND IMPLICATIONS

We have performed structural chemogenomics studies that identify links between (hist)aminergic receptor ligands and their binding sites and binding modes. This knowledge can be used to identify structure-selectivity relationships that increase our understanding of ligand binding to (hist)aminergic receptors and hence can be used in future GPCR ligand discovery and design.

The receptor concept in 3D: from hypothesis and metaphor to GPCR-ligand structures

The first mentioning of the word "receptor" for the structure with which a bioactive compound should react for obtaining its specific influence on a physiological system goes back to the years around 1900. The receptor concept was adapted from the lock and key theory for the enzyme substrate and blockers interactions. Through the years the concept, in the beginning rather being a metaphor, not a model, was refined and became reality in recent years. Not only the structures of receptors were elucidated, also the receptor machineries were unraveled. Following a brief historical review we will describe how the recent breakthroughs in the experimental determination of G protein-coupled receptor (GPCR) crystal structures can be complemented by computational modeling, medicinal chemistry, biochemical, and molecular pharmacological studies to obtain new insights into the molecular determinants of GPCR-ligand binding and activation. We will furthermore discuss how this information can be used for structure-based discovery of novel GPCR ligands that bind specific (allosteric) binding sites with desired effects on GPCR functional activity.

Biofragments: an approach towards predicting protein function using biologically related fragments and its application to Mycobacterium tuberculosis CYP126

We present a novel fragment-based approach that tackles some of the challenges for chemical biology of predicting protein function. The general approach, which we have termed biofragments, comprises two key stages. First, a biologically relevant fragment library (biofragment library) can be designed and constructed from known sets of substrate-like ligands for a protein class of interest. Second, the library can be screened for binding to a novel putative ligand-binding protein from the same or similar class, and the characterization of hits provides insight into the basis of ligand recognition, selectivity, and function at the substrate level. As a proof-of-concept, we applied the biofragments approach to the functionally uncharacterized Mycobacterium tuberculosis (Mtb) cytochrome P450 isoform, CYP126. This led to the development of a tailored CYP biofragment library with notable 3D characteristics and a significantly higher screening hit rate (14%) than standard drug-like fragment libraries screened previously against Mtb CYP121 and 125 (4% and 1%, respectively). Biofragment hits were identified that make both substrate-like type-I and inhibitor-like type-II interactions with CYP126. A chemical-fingerprint-based substrate model was built from the hits and used to search a virtual TB metabolome, which led to the discovery that CYP126 has a strong preference for the recognition of aromatics and substrate-like type-I binding of chlorophenol moieties within the active site near the heme. Future catalytic analyses will be focused on assessing CYP126 for potential substrate oxidative dehalogenation.

ATP synthase in mycobacteria: special features and implications for a function as drug target

ATP synthase is a ubiquitous enzyme that is largely conserved across the kingdoms of life. This conservation is in accordance with its central role in chemiosmotic energy conversion, a pathway utilized by far by most living cells. On the other hand, in particular pathogenic bacteria whilst employing ATP synthase have to deal with energetically unfavorable conditions such as low oxygen tensions in the human host, e.g. Mycobacterium tuberculosis can survive in human macrophages for an extended time. It is well conceivable that such ATP synthases may carry idiosyncratic features that contribute to efficient ATP production. In this review genetic and biochemical data on mycobacterial ATP synthase are discussed in terms of rotary catalysis, stator composition, and regulation of activity. ATP synthase in mycobacteria is of particular interest as this enzyme has been validated as a target for promising new antibacterial drugs. A deeper understanding of the working of mycobacterial ATP synthase and its atypical features can provide insight in adaptations of bacterial energy metabolism. Moreover, pinpointing and understanding critical differences as compared with human ATP synthase may provide input for the design and development of selective ATP synthase inhibitors as antibacterials. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

Benchmarking of protein descriptor sets in proteochemometric modeling (part 2): modeling performance of 13 amino acid descriptor sets

Background

While a large body of work exists on comparing and benchmarking descriptors of molecular structures, a similar comparison of protein descriptor sets is lacking. Hence, in the current work a total of 13 amino acid descriptor sets have been benchmarked with respect to their ability of establishing bioactivity models. The descriptor sets included in the study are Z-scales (3 variants), VHSE, T-scales, ST-scales, MS-WHIM, FASGAI, BLOSUM, a novel protein descriptor set (termed ProtFP (4 variants)), and in addition we created and benchmarked three pairs of descriptor combinations. Prediction performance was evaluated in seven structure-activity benchmarks which comprise Angiotensin Converting Enzyme (ACE) dipeptidic inhibitor data, and three proteochemometric data sets, namely (1) GPCR ligands modeled against a GPCR panel, (2) enzyme inhibitors (NNRTIs) with associated bioactivities against a set of HIV enzyme mutants, and (3) enzyme inhibitors (PIs) with associated bioactivities on a large set of HIV enzyme mutants.

Results

The amino acid descriptor sets compared here show similar performance (<0.1 log units RMSE difference and <0.1 difference in MCC), while errors for individual proteins were in some cases found to be larger than those resulting from descriptor set differences ( > 0.3 log units RMSE difference and >0.7 difference in MCC). Combining different descriptor sets generally leads to better modeling performance than utilizing individual sets. The best performers were Z-scales (3) combined with ProtFP (Feature), or Z-Scales (3) combined with an average Z-Scale value for each target, while ProtFP (PCA8), ST-Scales, and ProtFP (Feature) rank last.

Conclusions

While amino acid descriptor sets capture different aspects of amino acids their ability to be used for bioactivity modeling is still – on average – surprisingly similar. Still, combining sets describing complementary information consistently leads to small but consistent improvement in modeling performance (average MCC 0.01 better, average RMSE 0.01 log units lower). Finally, performance differences exist between the targets compared thereby underlining that choosing an appropriate descriptor set is of fundamental for bioactivity modeling, both from the ligand- as well as the protein side.

Fragment-based lead discovery on G-protein-coupled receptors

INTRODUCTION

G-protein-coupled receptors (GPCRs) form one of the largest groups of potential targets for novel medications. Low druggability of many GPCR targets and inefficient sampling of chemical space in high-throughput screening expertise however often hinder discovery of drug discovery leads for GPCRs. Fragment-based drug discovery is an alternative approach to the conventional strategy and has proven its efficiency on several enzyme targets. Based on developments in biophysical screening techniques, receptor stabilization and in vitro assays, virtual and experimental fragment screening and fragment-based lead discovery recently became applicable for GPCR targets.

AREAS COVERED

This article provides a review of the biophysical as well as biological detection techniques suitable to study GPCRs together with their applications to screen fragment libraries and identify fragment-size ligands of cell surface receptors. The article presents several recent examples including both virtual and experimental protocols for fragment hit discovery and early hit to lead progress.

EXPERT OPINION

With the recent progress in biophysical detection techniques, the advantages of fragment-based drug discovery could be exploited for GPCR targets. Structural information on GPCRs will be more abundantly available for early stages of drug discovery projects, providing information on the binding process and efficiently supporting the progression of fragment hit to lead. In silico approaches in combination with biological assays can be used to address structurally challenging GPCRs and confirm biological relevance of interaction early in the drug discovery project.

Design and Evaluation of the Performance of an NMR Screening Fragment Library

The design of a suitable library is an essential prerequisite to establish a fragment-based screening capability. Several pharmaceutical companies have described their approaches to establishing fragment libraries; however there are few detailed reports of both design and analysis of performance for a fragment library maintained in an academic setting. Here we report our efforts towards the design of a fragment library for nuclear magnetic resonance spectroscopy-based screening, demonstrate the performance of the library through analysis of 14 screens, and present a comparison to previously reported fragment libraries.