Detection of ligand binding hot spots on protein surfaces via fragment-based methods: application to DJ-1 and glucocerebrosidase

Expanding FTMap for Fragment-Based Identification of Pharmacophore Regions in Ligand Binding Sites

The knowledge of ligand binding hot spots and of the important interactions within such hot spots is crucial for the design of lead compounds in the early stages of structure-based drug discovery. The computational solvent mapping server FTMap can reliably identify binding hot spots as consensus clusters, free energy minima that bind a variety of organic probe molecules. However, in its current implementation, FTMap provides limited information on regions within the hot spots that tend to interact with specific pharmacophoric features of potential ligands. E-FTMap is a new server that expands on the original FTMap protocol. E-FTMap uses 119 organic probes, rather than the 16 in the original FTMap, to exhaustively map binding sites, and identifies pharmacophore features as atomic consensus sites where similar chemical groups bind. We validate E-FTMap against a set of 109 experimentally derived structures of fragment-lead pairs, finding that highly ranked pharmacophore features overlap with the corresponding atoms in both fragments and lead compounds. Additionally, comparisons of mapping results to ensembles of bound ligands reveal that pharmacophores generated with E-FTMap tend to sample highly conserved protein-ligand interactions. E-FTMap is available as a web server at https://eftmap.bu.edu.

Fragment library screening by X-ray crystallography and binding site analysis on thioredoxin glutathione reductase of Schistosoma mansoni

Schistosomiasis is caused by parasites of the genus Schistosoma, which infect more than 200 million people. Praziquantel (PZQ) has been the main drug for controlling schistosomiasis for over four decades, but despite that it is ineffective against juvenile worms and size and taste issues with its pharmaceutical forms impose challenges for treating school-aged children. It is also important to note that PZQ resistant strains can be generated in laboratory conditions and observed in the field, hence its extensive use in mass drug administration programs raises concerns about resistance, highlighting the need to search for new schistosomicidal drugs. Schistosomes survival relies on the redox enzyme thioredoxin glutathione reductase (TGR), a validated target for the development of new anti-schistosomal drugs. Here we report a high-throughput fragment screening campaign of 768 compounds against S. mansoni TGR (SmTGR) using X-ray crystallography. We observed 49 binding events involving 35 distinct molecular fragments which were found to be distributed across 16 binding sites. Most sites are described for the first time within SmTGR, a noteworthy exception being the "doorstop pocket" near the NADPH binding site. We have compared results from hotspots and pocket druggability analysis of SmTGR with the experimental binding sites found in this work, with our results indicating only limited coincidence between experimental and computational results. Finally, we discuss that binding sites at the doorstop/NADPH binding site and in the SmTGR dimer interface, should be prioritized for developing SmTGR inhibitors as new antischistosomal drugs.

In Silico Exploration of the Trypanothione Reductase (TryR) of L. mexicana

Human leishmaniasis is a neglected tropical disease which affects nearly 1.5 million people every year, with Mexico being an important endemic region. One of the major defense mechanisms of these parasites is based in the polyamine metabolic pathway, as it provides the necessary compounds for its survival. Among the enzymes in this route, trypanothione reductase (TryR), an oxidoreductase enzyme, is crucial for the Leishmania genus' survival against oxidative stress. Thus, it poses as an attractive drug target, yet due to the size and features of its catalytic pocket, modeling techniques such as molecular docking focusing on that region is not convenient. Herein, we present a computational study using several structure-based approaches to assess the druggability of TryR from L. mexicana, the predominant Leishmania species in Mexico, beyond its catalytic site. Using this consensus methodology, three relevant pockets were found, of which the one we call σ-site promises to be the most favorable one. These findings may help the design of new drugs of trypanothione-related diseases.

FTMove: A Web Server for Detection and Analysis of Cryptic and Allosteric Binding Sites by Mapping Multiple Protein Structures

Protein mapping distributes many copies of different molecular probes on the surface of a target protein in order to determine binding hot spots, regions that are highly preferable for ligand binding. While mapping of X-ray structures by the FTMap server is inherently static, this limitation can be overcome by the simultaneous analysis of multiple structures of the protein. FTMove is an automated web server that implements this approach. From the input of a target protein, by PDB code, the server identifies all structures of the protein available in the PDB, runs mapping on them, and combines the results to form binding hot spots and binding sites. The user may also upload their own protein structures, bypassing the PDB search for similar structures. Output of the server consists of the consensus binding sites and the individual mapping results for each structure - including the number of probes located in each binding site, for each structure. This level of detail allows the users to investigate how the strength of a binding site relates to the protein conformation, other binding sites, and the presence of ligands or mutations. In addition, the structures are clustered on the basis of their binding properties. The use of FTMove is demonstrated by application to 22 proteins with known allosteric binding sites; the orthosteric and allosteric binding sites were identified in all but one case, and the sites were typically ranked among the top five. The FTMove server is publicly available at https://ftmove.bu.edu.

Carnitine is a pharmacological allosteric chaperone of the human lysosomal α-glucosidase

Pompe disease is an inherited metabolic disorder due to the deficiency of the lysosomal acid α-glucosidase (GAA). The only approved treatment is enzyme replacement therapy with the recombinant enzyme (rhGAA). Further approaches like pharmacological chaperone therapy, based on the stabilising effect induced by small molecules on the target enzyme, could be a promising strategy. However, most known chaperones could be limited by their potential inhibitory effects on patient’s enzymes. Here we report on the discovery of novel chaperones for rhGAA, L- and D-carnitine, and the related compound acetyl-D-carnitine. These drugs stabilise the enzyme at pH and temperature without inhibiting the activity and acted synergistically with active-site directed pharmacological chaperones. Remarkably, they enhanced by 4-fold the acid α-glucosidase activity in fibroblasts from three Pompe patients with added rhGAA. This synergistic effect of L-carnitine and rhGAA has the potential to be translated into improved therapeutic efficacy of ERT in Pompe disease.

Conservation of binding properties in protein models

We study the models submitted to round 12 of the Critical Assessment of protein Structure Prediction (CASP) experiment to assess how well the binding properties are conserved when the X-ray structures of the target proteins are replaced by their models. To explore small molecule binding we generate distributions of molecular probes - which are fragment-sized organic molecules of varying size, shape, and polarity - around the protein, and count the number of interactions between each residue and the probes, resulting in a vector of interactions we call a binding fingerprint. The similarity between two fingerprints, one for the X-ray structure and the other for a model of the protein, is determined by calculating the correlation coefficient between the two vectors. The resulting correlation coefficients are shown to correlate with global measures of accuracy established in CASP, and the relationship yields an accuracy threshold that has to be reached for meaningful binding surface conservation. The clusters formed by the probe molecules reliably predict binding hot spots and ligand binding sites in both X-ray structures and reasonably accurate models of the target, but ensembles of models may be needed for assessing the availability of proper binding pockets. We explored ligand docking to the few targets that had bound ligands in the X-ray structure. More targets were available to assess the ability of the models to reproduce protein-protein interactions by docking both the X-ray structures and models to their interaction partners in complexes. It was shown that this application is more difficult than finding small ligand binding sites, and the success rates heavily depend on the local structure in the potential interface. In particular, predicted conformations of flexible loops are frequently incorrect in otherwise highly accurate models, and may prevent predicting correct protein-protein interactions.

Inhibition of intracellular Ca2+ mobilization and potassium channels activation are involved in the vasorelaxation induced by 7-hydroxycoumarin

Coumarins exhibit a wide variety of biological effects, including activities in the cardiovascular system and the aim of this study was to evaluate the vascular therapeutic potential of 7-Hydroxicoumarin (7-HC). The vascular effects induced by 7-HC (0.001 μM-300 μM), were investigated by in vitro approaches using isometric tension measurements in rat superior mesenteric arteries and by in silico assays using Ligand-based analysis. Our results suggest that the vasorelaxant effect of 7-HC seems to rely on potassium channels, notably through large conductance Ca²⁺-activated K⁺ (BK_Ca) channels activation. In fact, 7-HC (300 μM) significantly reduced CaCl₂-induced contraction as well as the reduction of intracellular calcium mobilization. However, the relaxation induced by 7-HC was independent of store-operated calcium entry (SOCE). Moreover, in silico analysis suggests that potassium channels have a common binding pocket, where 7-HC may bind and hint that its binding profile is more similar to quinine's than verapamil's. These results are compatible with the inhibition of Ca²⁺ release from intracellular stores, which is prompted by phenylephrine and caffeine. Taken together, these results demonstrate a therapeutic potential of 7-HC on the cardiovascular system, making it a promising lead compound for the development of drugs useful in the treatment of cardiovascular diseases.

Trapping conformational states of a flavin-dependent N-monooxygenase in crystallo reveals protein and flavin dynamics

The siderophore biosynthetic enzyme A (SidA) ornithine hydroxylase from Aspergillus fumigatus is a fungal disease drug target involved in the production of hydroxamate-containing siderophores, which are used by the pathogen to sequester iron. SidA is an N-monooxygenase that catalyzes the NADPH-dependent hydroxylation of l-ornithine through a multistep oxidative mechanism, utilizing a C4a-hydroperoxyflavin intermediate. Here we present four new crystal structures of SidA in various redox and ligation states, including the first structure of oxidized SidA without NADP(H) or l-ornithine bound (resting state). The resting state structure reveals a new out active site conformation characterized by large rotations of the FAD isoalloxazine around the C1-'C2' and N10-C1' bonds, coupled to a 10-Å movement of the Tyr-loop. Additional structures show that either flavin reduction or the binding of NADP(H) is sufficient to drive the FAD to the in conformation. The structures also reveal protein conformational changes associated with the binding of NADP(H) and l-ornithine. Some of these residues were probed using site-directed mutagenesis. Docking was used to explore the active site of the out conformation. These calculations identified two potential ligand-binding sites. Altogether, our results provide new information about conformational dynamics in flavin-dependent monooxygenases. Understanding the different active site conformations that appear during the catalytic cycle may allow fine-tuning of inhibitor discovery efforts.

A Novel Pyrazolopyrimidine Ligand of Human PGK1 and Stress Sensor DJ1 Modulates the Shelterin Complex and Telomere Length Regulation

Telomere signaling and metabolic dysfunction are hallmarks of cell aging. New agents targeting these processes might provide therapeutic opportunities, including chemoprevention strategies against cancer predisposition. We report identification and characterization of a pyrazolopyrimidine compound series identified from screens focused on cell immortality and whose targets are glycolytic kinase PGK1 and oxidative stress sensor DJ1. We performed structure–activity studies on the series to develop a photoaffinity probe to deconvolute the cellular targets. In vitro binding and structural analyses confirmed these targets, suggesting that PGK1/DJ1 interact, which we confirmed by immunoprecipitation. Glucose homeostasis and oxidative stress are linked to telomere signaling and exemplar compound CRT0063465 blocked hypoglycemic telomere shortening. Intriguingly, PGK1 and DJ1 bind to TRF2 and telomeric DNA. Compound treatment modulates these interactions and also affects Shelterin complex composition, while conferring cellular protection from cytotoxicity due to bleomycin and desferroxamine. These results demonstrate therapeutic potential of the compound series.

Why Some Targets Benefit from beyond Rule of Five Drugs

Beyond rule-of-five (bRo5) compounds are increasingly used in drug discovery. Here we analyze 37 target proteins that have bRo5 drugs or clinical candidates. Targets can benefit from bRo5 drugs if they have "complex" hot spot structure with four or more hots spots, including some strong ones. Complex I targets show positive correlation between binding affinity and molecular weight. These targets are conventionally druggable, but reaching additional hot spots enables improved pharmaceutical properties. Complex II targets, mostly protein kinases, also have strong hot spots but show no correlation between affinity and ligand molecular weight, and the primary motivation for creating larger drugs is to increase selectivity. Each target considered as complex III has some specific reason for requiring bRo5 drugs. Finally, targets with "simple" hot spot structure, i.e., three or fewer weak hot spots, must use larger compounds that interact with surfaces beyond the hot spot region to achieve acceptable affinity.

Analysis of tractable allosteric sites in G protein-coupled receptors

Allosteric modulation of G protein-coupled receptors represent a promising mechanism of pharmacological intervention. Dramatic developments witnessed in the structural biology of membrane proteins continue to reveal that the binding sites of allosteric modulators are widely distributed, including along protein surfaces. Here we restrict consideration to intrahelical and intracellular sites together with allosteric conformational locks, and show that the protein mapping tools FTMap and FTSite identify 83% and 88% of such experimentally confirmed allosteric sites within the three strongest sites found. The methods were also able to find partially hidden allosteric sites that were not fully formed in X-ray structures crystallized in the absence of allosteric ligands. These results confirm that the intrahelical sites capable of binding druglike allosteric modulators are among the strongest ligand recognition sites in a large fraction of GPCRs and suggest that both FTMap and FTSite are useful tools for identifying allosteric sites and to aid in the design of such compounds in a range of GPCR targets.

In Silico Studies in Drug Research Against Neurodegenerative Diseases

Background

Neurodegenerative diseases such as Alzheimer's disease (AD), amyotrophic lateral sclerosis, Parkinson's disease (PD), spinal cerebellar ataxias, and spinal and bulbar muscular atrophy are described by slow and selective degeneration of neurons and axons in the central nervous system (CNS) and constitute one of the major challenges of modern medicine. Computer-aided or in silico drug design methods have matured into powerful tools for reducing the number of ligands that should be screened in experimental assays.

Methods

In the present review, the authors provide a basic background about neurodegenerative diseases and in silico techniques in the drug research. Furthermore, they review the various in silico studies reported against various targets in neurodegenerative diseases, including homology modeling, molecular docking, virtual high-throughput screening, quantitative structure activity relationship (QSAR), hologram quantitative structure activity relationship (HQSAR), 3D pharmacophore mapping, proteochemometrics modeling (PCM), fingerprints, fragment-based drug discovery, Monte Carlo simulation, molecular dynamic (MD) simulation, quantum-mechanical methods for drug design, support vector machines, and machine learning approaches.

Results

Detailed analysis of the recently reported case studies revealed that the majority of them use a sequential combination of ligand and structure-based virtual screening techniques, with particular focus on pharmacophore models and the docking approach.

Conclusion

Neurodegenerative diseases have a multifactorial pathoetiological origin, so scientists have become persuaded that a multi-target therapeutic strategy aimed at the simultaneous targeting of multiple proteins (and therefore etiologies) involved in the development of a disease is recommended in future.

Discovery and Optimization of Inhibitors of the Parkinson's Disease Associated Protein DJ-1

DJ-1 is a Parkinson's disease associated protein endowed with enzymatic, redox sensing, regulatory, chaperoning, and neuroprotective activities. Although DJ-1 has been vigorously studied for the past decade and a half, its exact role in the progression of the disease remains uncertain. In addition, little is known about the spatiotemporal regulation of DJ-1, or the biochemical basis explaining its numerous biological functions. Progress has been hampered by the lack of inhibitors with precisely known mechanisms of action. Herein, we have employed biophysical methodologies and X-ray crystallography to identify and to optimize a family of compounds inactivating the critical Cys106 residue of human DJ-1. We demonstrate these compounds are potent inhibitors of various activities of DJ-1 in vitro and in cell-based assays. This study reports a new family of DJ-1 inhibitors with a defined mechanism of action, and contributes toward the understanding of the biological function of DJ-1.

Stereodivergent synthesis of right- and left-handed iminoxylitol heterodimers and monomers. Study of their impact on β-glucocerebrosidase activity

A library of dimers and heterodimers of both enantiomers of 2-O-alkylated iminoxylitol derivatives has been synthesised and evaluated on β-glucocerebrosidase (GCase), the enzyme responsible for Gaucher disease (GD). Although the objective was to target simultaneously the active site and a secondary binding site of the glucosidase, the (-)-2-iminoxylitol moiety seemed detrimental for imiglucerase inhibition and no significant enhancement was obtained in G202R, N370S and L444P fibroblasts. However, all compounds having at least one (+)-2-O-alkyl iminoxylitol are GCase inhibitors in the nano molar range and are significant GCase activity enhancers in G202R fibroblats, as confirmed by a decrease of glucosylceramide levels and by co-localization studies.

Structural Biology of the DJ-1 Superfamily

The DJ-1 (also called the DJ-1/PfpI, ThiJ/PfpI, or DJ-1/ThiJ/PfpI) superfamily is a structural and functional diverse group of proteins that are present in most organisms. Many of these proteins remain poorly characterized at the biochemical level, but include some known chaperones, proteases, and various stress response proteins that remain mechanistically mysterious. This chapter outlines what is known from a structural perspective about the cellular and biochemical functions of many of these proteins from distinct clades of the superfamily in several organisms. In humans, DJ-1 appears to function primarily as a redox-responsive protein that may act as a sensor for imbalances in cellular redox state. Because mutations in human DJ-1 cause certain types of heritable Parkinson's disease, the role of oxidative posttranslational modifications and pathogenic mutations in human DJ-1 is emphasized in the latter sections of this chapter.

DJ-1 as a Therapeutic Target Against Cancer

DJ-1 is a gene involved in various cellular processes, including transcriptional regulation, oxidative stress response, fertilization, mitochondrial regulation, inflammatory and fibrogenic niche formation, and glycation damage prevention. Although a disease-associated genetic study within the past decade has demonstrated that the mutation of DJ-1 is associated with autosomal early-onset Parkinson's disease, increasing evidence suggests that DJ-1 also plays a critical role in tumor development and progression. In this review, we provide an overview of current knowledge concerning the role and the mechanism of DJ-1 in cancer and also discuss the possibility of DJ-1 as a therapeutic target against cancer.

Exploring the structural origins of cryptic sites on proteins

Molecular dynamics (MD) simulations of proteins reveal the existence of many transient surface pockets; however, the factors determining what small subset of these represent druggable or functionally relevant ligand binding sites, called "cryptic sites," are not understood. Here, we examine multiple X-ray structures for a set of proteins with validated cryptic sites, using the computational hot spot identification tool FTMap. The results show that cryptic sites in ligand-free structures generally have a strong binding energy hot spot very close by. As expected, regions around cryptic sites exhibit above-average flexibility, and close to 50% of the proteins studied here have unbound structures that could accommodate the ligand without clashes. Nevertheless, the strong hot spot neighboring each cryptic site is almost always exploited by the bound ligand, suggesting that binding may frequently involve an induced fit component. We additionally evaluated the structural basis for cryptic site formation, by comparing unbound to bound structures. Cryptic sites are most frequently occluded in the unbound structure by intrusion of loops (22.5%), side chains (19.4%), or in some cases entire helices (5.4%), but motions that create sites that are too open can also eliminate pockets (19.4%). The flexibility of cryptic sites frequently leads to missing side chains or loops (12%) that are particularly evident in low resolution crystal structures. An interesting observation is that cryptic sites formed solely by the movement of side chains, or of backbone segments with fewer than five residues, result only in low affinity binding sites with limited use for drug discovery.

Testing inhomogeneous solvation theory in structure-based ligand discovery

Binding-site water is often displaced upon ligand recognition, but is commonly neglected in structure-based ligand discovery. Inhomogeneous solvation theory (IST) has become popular for treating this effect, but it has not been tested in controlled experiments at atomic resolution. To do so, we turned to a grid-based version of this method, GIST, readily implemented in molecular docking. Whereas the term only improves docking modestly in retrospective ligand enrichment, it could be added without disrupting performance. We thus turned to prospective docking of large libraries to investigate GIST's impact on ligand discovery, geometry, and water structure in a model cavity site well-suited to exploring these terms. Although top-ranked docked molecules with and without the GIST term often overlapped, many ligands were meaningfully prioritized or deprioritized; some of these were selected for testing. Experimentally, 13/14 molecules prioritized by GIST did bind, whereas none of the molecules that it deprioritized were observed to bind. Nine crystal complexes were determined. In six, the ligand geometry corresponded to that predicted by GIST, for one of these the pose without the GIST term was wrong, and three crystallographic poses differed from both predictions. Notably, in one structure, an ordered water molecule with a high GIST displacement penalty was observed to stay in place. Inclusion of this water-displacement term can substantially improve the hit rates and ligand geometries from docking screens, although the magnitude of its effects can be small and its impact in drug binding sites merits further controlled studies.

Do Fragments and Crystallization Additives Bind Similarly to Drug-like Ligands?

The success of fragment-based drug design (FBDD) hinges upon the optimization of low-molecular-weight compounds (MW < 300 Da) with weak binding affinities to lead compounds with high affinity and selectivity. Usually, structural information from fragment-protein complexes is used to develop ideas about the binding mode of similar but drug-like molecules. In this regard, crystallization additives such as cryoprotectants or buffer components, which are highly abundant in crystal structures, are frequently ignored. Thus, the aim of this study was to investigate the information present in protein complexes with fragments as well as those with additives and how they relate to the binding modes of their drug-like counterparts. We present a thorough analysis of the binding modes of crystallographic additives, fragments, and drug-like ligands bound to four diverse targets of wide interest in drug discovery and highly represented in the Protein Data Bank: cyclin-dependent kinase 2, β-secretase 1, carbonic anhydrase 2, and trypsin. We identified a total of 630 unique molecules bound to the catalytic binding sites, among them 31 additives, 222 fragments, and 377 drug-like ligands. In general, we observed that, independent of the target, protein-fragment interaction patterns are highly similar to those of drug-like ligands and mostly cover the residues crucial for binding. Crystallographic additives are also able to show conserved binding modes and recover the residues important for binding in some of the cases. Moreover, we show evidence that the information from fragments and drug-like ligands can be applied to rescore docking poses in order to improve the prediction of binding modes.

Hot spot mapping of protein surfaces with TEMPOL: Bovine pancreatic RNase A as a model system

TEMPOL spin-label has been used to identify surface exposure of protein nuclei from NMR analysis of the induced paramagnetic relaxation enhancements (PRE). The absence of linear dependence between atom depths and observed PRE reveals that specific mechanisms drive the approach of the paramagnet to the protein surface. RNase A represents a unique protein system to explore the fine details of the information offered by TEMPOL induced PRE, due to the abundance of previous results, obtained in solution and in the crystal, dealing with surface dynamics behavior of this protein. MD simulations in explicit solvent have been performed, also in the presence of TEMPOL, in order to delineate the role of intermolecular hydrogen bonds (HB) on PRE extents. Comparison of our results with the ones obtained from multiple solvent crystal structure (MSCS) studies yields information on the specificities that these two techniques have for characterizing protein-ligand interactions, a fundamental step in the development of reliable surface druggability predictors.

Driving Structure-Based Drug Discovery through Cosolvent Molecular Dynamics

Identifying binding hotspots on protein surfaces is of prime interest in structure-based drug discovery, either to assess the tractability of pursuing a protein target or to drive improved potency of lead compounds. Computational approaches to detect such regions have traditionally relied on energy minimization of probe molecules onto static protein conformations in the absence of the natural aqueous environment. Advances in high performance computing now allow us to assess hotspots using molecular dynamics (MD) simulations. MD simulations integrate protein flexibility and the complicated role of water, thereby providing a more realistic assessment of the complex kinetics and thermodynamics at play. In this review, we describe the evolution of various cosolvent-based MD techniques and highlight a myriad of potential applications for such technologies in computational drug development.

From structural biology to designing therapy for inborn errors of metabolism

At the SSIEM Symposium in Istanbul 2010, I presented an overview of protein structural approaches in the study of inborn errors of metabolism (Yue and Oppermann 2011). Five years on, the field is going strong with new protein structures, uncovered catalytic functions and novel chemical matters for metabolic enzymes, setting the stage for the next generation of drug discovery. This article aims to update on recent advances and lessons learnt on inborn errors of metabolism via the protein-centric approach, citing examples of work from my group, collaborators and co-workers that cover diverse pathways of transsulfuration, cobalamin and glycogen metabolism. Taking into consideration that many inborn errors of metabolism result in the loss of enzyme function, this presentation aims to outline three key principles that guide the design of small molecule therapy in this technically challenging field: (1) integrating structural, biochemical and cell-based data to evaluate the wide spectrum of mutation-driven enzyme defects in stability, catalysis and protein-protein interaction; (2) studying multi-domain proteins and multi-protein complexes as examples from nature, to learn how enzymes are activated by small molecules; (3) surveying different regions of the enzyme, away from its active site, that can be targeted for the design of allosteric activators and inhibitors.

Pharmacological Chaperones: Design and Development of New Therapeutic Strategies for the Treatment of Conformational Diseases

Errors in protein folding may result in premature clearance of structurally aberrant proteins, or in the accumulation of toxic misfolded species or protein aggregates. These pathological events lead to a large range of conditions known as conformational diseases. Several research groups have presented possible therapeutic solutions for their treatment by developing novel compounds, known as pharmacological chaperones. These cell-permeable molecules selectively provide a molecular scaffold around which misfolded proteins can recover their native folding and, thus, their biological activities. Here, we review therapeutic strategies, clinical potentials, and cost-benefit impacts of several classes of pharmacological chaperones for the treatment of a series of conformational diseases.

Focused grid-based resampling for protein docking and mapping

The fast Fourier transform (FFT) sampling algorithm has been used with success in application to protein-protein docking and for protein mapping, the latter docking a variety of small organic molecules for the identification of binding hot spots on the target protein. Here we explore the local rather than global usage of the FFT sampling approach in docking applications. If the global FFT based search yields a near-native cluster of docked structures for a protein complex, then focused resampling of the cluster generally leads to a substantial increase in the number of conformations close to the native structure. In protein mapping, focused resampling of the selected hot spot regions generally reveals further hot spots that, while not as strong as the primary hot spots, also contribute to ligand binding. The detection of additional ligand binding regions is shown by the improved overlap between hot spots and bound ligands.

Lessons from Hot Spot Analysis for Fragment-Based Drug Discovery

Analysis of binding energy hot spots at protein surfaces can provide crucial insights into the prospects for successful application of fragment-based drug discovery (FBDD), and whether a fragment hit can be advanced into a high-affinity, drug-like ligand. The key factor is the strength of the top ranking hot spot, and how well a given fragment complements it. We show that published data are sufficient to provide a sophisticated and quantitative understanding of how hot spots derive from a protein 3D structure, and how their strength, number, and spatial arrangement govern the potential for a surface site to bind to fragment-sized and larger ligands. This improved understanding provides important guidance for the effective application of FBDD in drug discovery.

New Frontiers in Druggability

A powerful early approach to evaluating the druggability of proteins involved determining the hit rate in NMR-based screening of a library of small compounds. Here, we show that a computational analog of this method, based on mapping proteins using small molecules as probes, can reliably reproduce druggability results from NMR-based screening and can provide a more meaningful assessment in cases where the two approaches disagree. We apply the method to a large set of proteins. The results show that, because the method is based on the biophysics of binding rather than on empirical parametrization, meaningful information can be gained about classes of proteins and classes of compounds beyond those resembling validated targets and conventionally druglike ligands. In particular, the method identifies targets that, while not druggable by druglike compounds, may become druggable using compound classes such as macrocycles or other large molecules beyond the rule-of-five limit.

Ligand deconstruction: Why some fragment binding positions are conserved and others are not

Fragment-based drug discovery (FBDD) relies on the premise that the fragment binding mode will be conserved on subsequent expansion to a larger ligand. However, no general condition has been established to explain when fragment binding modes will be conserved. We show that a remarkably simple condition can be developed in terms of how fragments coincide with binding energy hot spots--regions of the protein where interactions with a ligand contribute substantial binding free energy--the locations of which can easily be determined computationally. Because a substantial fraction of the free energy of ligand binding comes from interacting with the residues in the energetically most important hot spot, a ligand moiety that sufficiently overlaps with this region will retain its location even when other parts of the ligand are removed. This hypothesis is supported by eight case studies. The condition helps identify whether a protein is suitable for FBDD, predicts the size of fragments required for screening, and determines whether a fragment hit can be extended into a higher affinity ligand. Our results show that ligand binding sites can usefully be thought of in terms of an anchor site, which is the top-ranked hot spot and dominates the free energy of binding, surrounded by a number of weaker satellite sites that confer improved affinity and selectivity for a particular ligand and that it is the intrinsic binding potential of the protein surface that determines whether it can serve as a robust binding site for a suitably optimized ligand.

Pharmacological Chaperone Therapy: Preclinical Development, Clinical Translation, and Prospects for the Treatment of Lysosomal Storage Disorders

Lysosomal storage disorders (LSDs) are a group of inborn metabolic diseases caused by mutations in genes that encode proteins involved in different lysosomal functions, in most instances acidic hydrolases. Different therapeutic approaches have been developed to treat these disorders. Pharmacological chaperone therapy (PCT) is an emerging approach based on small-molecule ligands that selectively bind and stabilize mutant enzymes, increase their cellular levels, and improve lysosomal trafficking and activity. Compared to other approaches, PCT shows advantages, particularly in terms of oral administration, broad biodistribution, and positive impact on patients' quality of life. After preclinical in vitro and in vivo studies, PCT is now being translated in the first clinical trials, either as monotherapy or in combination with enzyme replacement therapy, for some of the most prevalent LSDs. For some LSDs, the results of the first clinical trials are encouraging and warrant further development. Future research in the field of PCT will be directed toward the identification of novel chaperones, including new allosteric drugs, and the exploitation of synergies between chaperone treatment and other therapeutic approaches.

The FTMap family of web servers for determining and characterizing ligand-binding hot spots of proteins

FTMap is a computational mapping server that identifies binding hot spots of macromolecules-i.e., regions of the surface with major contributions to the ligand-binding free energy. To use FTMap, users submit a protein, DNA or RNA structure in PDB (Protein Data Bank) format. FTMap samples billions of positions of small organic molecules used as probes, and it scores the probe poses using a detailed energy expression. Regions that bind clusters of multiple probe types identify the binding hot spots in good agreement with experimental data. FTMap serves as the basis for other servers, namely FTSite, which is used to predict ligand-binding sites, FTFlex, which is used to account for side chain flexibility, FTMap/param, used to parameterize additional probes and FTDyn, for mapping ensembles of protein structures. Applications include determining the druggability of proteins, identifying ligand moieties that are most important for binding, finding the most bound-like conformation in ensembles of unliganded protein structures and providing input for fragment-based drug design. FTMap is more accurate than classical mapping methods such as GRID and MCSS, and it is much faster than the more-recent approaches to protein mapping based on mixed molecular dynamics. By using 16 probe molecules, the FTMap server finds the hot spots of an average-size protein in <1 h. As FTFlex performs mapping for all low-energy conformers of side chains in the binding site, its completion time is proportionately longer.

Pharmacological chaperone therapy for lysosomal storage diseases

Pharmacological chaperone therapy is an emerging approach to treat lysosomal storage diseases. Small-molecule chaperones interact with mutant enzymes, favor their correct conformation and enhance their stability. This approach shows significant advantages when compared with existing therapies, particularly in terms of the bioavailability of drugs, oral administration and positive impact on the quality of patients' lives. On the other hand, future research in this field must confront important challenges. The identification of novel chaperones is indispensable to expanding the number of patients amenable to this treatment and to optimize therapeutic efficacy. It is important to develop new allosteric drugs, to address the risk of inhibiting target enzymes. Future research must also be directed towards the exploitation of synergies between chaperone treatment and other therapeutic approaches.

DJ-1 as a human oncogene and potential therapeutic target

DJ-1 is a cancer- and Parkinson's disease-associated protein that participates in different intracellular signaling pathways to protect cells from toxic stresses. DJ-1 expression, oxidation, localization, and phosphorylation are often altered in human tumors, and DJ-1 has been implicated in various aspects of transformation, including uncontrolled proliferation, invasion, metastasis, and resistance to chemotherapy and apoptosis. Despite the strong relationship between DJ-1 and cancer, which made it a particularly attractive therapeutic target for cancer treatment, the detailed mechanisms of how this oncogene coordinates altered signaling with cell survival remains elusive. In this commentary, we discuss the role of DJ-1 in transformation, highlight some of the significant aspects of and prospects for therapeutically targeting the DJ-1 signaling in cancer, and describe what the future may hold.

Potential pharmacological chaperones targeting cancer-associated MCL-1 and Parkinson disease-associated α-synuclein

Pharmacological chaperones are small molecules that bind to proteins and stabilize them against thermal denaturation or proteolytic degradation, as well as assist or prevent certain protein-protein assemblies. These activities are being exploited for the development of treatments for diseases caused by protein instability and/or aberrant protein-protein interactions, such as those found in certain forms of cancers and neurodegenerative diseases. However, designing or discovering pharmacological chaperones for specific targets is challenging because of the relatively featureless protein target surfaces, the lack of suitable chemical libraries, and the shortage of efficient high-throughput screening methods. In this study, we attempted to address all these challenges by synthesizing a diverse library of small molecules that mimic protein α-helical secondary structures commonly found in protein-protein interaction surfaces. This was accompanied by establishing a facile "on-bead" high-throughput screening method that allows for rapid and efficient discovery of potential pharmacological chaperones and for identifying novel chaperones/inhibitors against a cancer-associated protein, myeloid cell leukemia 1 (MCL-1), and a Parkinson disease-associated protein, α-synuclein. Our data suggest that the compounds and methods described here will be useful tools for the development of pharmaceuticals for complex-disease targets that are traditionally deemed "undruggable."

Identifying ligand-binding hot spots in proteins using brominated fragments

High-quality crystals of Thermus thermophilus EF-Tu in the GTP-bound conformation at 1.7-2.7 Å resolution were used to test 18 small organic molecules, all brominated for confident identification in the anomalous difference maps. From this relatively small collection, it was possible to identify a small molecule bound in the functionally important tRNA CCA-end binding pocket. The antibiotic GE2270 A is known to interact with the same pocket in EF-Tu and to disrupt the association with tRNA. Bromide could be located from peaks in the anomalous map in data truncated to very low resolution without refining the structure. Considering the speed with which diffraction data can be collected today, it is proposed that it is worthwhile to collect the extra data from fragment screens while crystals are at hand to increase the knowledge of biological function and drug binding in an experimental structural context.

Comprehensive experimental and computational analysis of binding energy hot spots at the NF-κB essential modulator/IKKβ protein-protein interface

We report a comprehensive analysis of binding energy hot spots at the protein-protein interaction (PPI) interface between nuclear factor kappa B (NF-κB) essential modulator (NEMO) and IκB kinase subunit β (IKKβ), an interaction that is critical for NF-κB pathway signaling, using experimental alanine scanning mutagenesis and also the FTMap method for computational fragment screening. The experimental results confirm that the previously identified NEMO binding domain (NBD) region of IKKβ contains the highest concentration of hot-spot residues, the strongest of which are W739, W741, and L742 (ΔΔG = 4.3, 3.5, and 3.2 kcal/mol, respectively). The region occupied by these residues defines a potentially druggable binding site on NEMO that extends for ~16 Å to additionally include the regions that bind IKKβ L737 and F734. NBD residues D738 and S740 are also important for binding but do not make direct contact with NEMO, instead likely acting to stabilize the active conformation of surrounding residues. We additionally found two previously unknown hot-spot regions centered on IKKβ residues L708/V709 and L719/I723. The computational approach successfully identified all three hot-spot regions on IKKβ. Moreover, the method was able to accurately quantify the energetic importance of all hot-spot residues involving direct contact with NEMO. Our results provide new information to guide the discovery of small-molecule inhibitors that target the NEMO/IKKβ interaction. They additionally clarify the structural and energetic complementarity between "pocket-forming" and "pocket-occupying" hot-spot residues, and further validate computational fragment mapping as a method for identifying hot spots at PPI interfaces.

Misfolding of galactose 1-phosphate uridylyltransferase can result in type I galactosemia

Type I galactosemia is a genetic disorder that is caused by the impairment of galactose-1-phosphate uridylyltransferase (GALT; EC 2.7.7.12). Although a large number of mutations have been detected through genetic screening of the human GALT (hGALT) locus, for many it is not known how they cause their effects. The majority of these mutations are missense, with predicted substitutions scattered throughout the enzyme structure and thus causing impairment by other means rather than direct alterations to the active site. To clarify the fundamental, molecular basis of hGALT impairment we studied five disease-associated variants p.D28Y, p.L74P, p.F171S, p.F194L and p.R333G using both a yeast model and purified, recombinant proteins. In a yeast expression system there was a correlation between lysate activity and the ability to rescue growth in the presence of galactose, except for p.R333G. Kinetic analysis of the purified proteins quantified each variant's level of enzymatic impairment and demonstrated that this was largely due to altered substrate binding. Increased surface hydrophobicity, altered thermal stability and changes in proteolytic sensitivity were also detected. Our results demonstrate that hGALT requires a level of flexibility to function optimally and that altered folding is the underlying reason of impairment in all the variants tested here. This indicates that misfolding is a common, molecular basis of hGALT deficiency and suggests the potential of pharmacological chaperones and proteostasis regulators as novel therapeutic approaches for type I galactosemia.

Fragment-guided design of subnanomolar β-lactamase inhibitors active in vivo

Fragment-based design was used to guide derivatization of a lead series of β-lactamase inhibitors that had heretofore resisted optimization for in vivo activity. X-ray structures of fragments overlaid with the lead suggested new, unanticipated functionality and points of attachment. Synthesis of three derivatives improved affinity over 20-fold and improved efficacy in cell culture. Crystal structures were consistent with the fragment-based design, enabling further optimization to a K(i) of 50 pM, a 500-fold improvement that required the synthesis of only six derivatives. One of these, compound 5, was tested in mice. Whereas cefotaxime alone failed to cure mice infected with β-lactamase-expressing Escherichia coli, 65% were cleared of infection when treated with a cefotaxime:5 combination. Fragment complexes offer a path around design hurdles, even for advanced molecules; the series described here may provide leads to overcome β-lactamase-based resistance, a key clinical challenge.

High-Throughput parallel blind Virtual Screening using BINDSURF

Background

Virtual Screening (VS) methods can considerably aid clinical research, predicting how ligands interact with drug targets. Most VS methods suppose a unique binding site for the target, usually derived from the interpretation of the protein crystal structure. However, it has been demonstrated that in many cases, diverse ligands interact with unrelated parts of the target and many VS methods do not take into account this relevant fact.

Results

We present BINDSURF, a novel VS methodology that scans the whole protein surface in order to find new hotspots, where ligands might potentially interact with, and which is implemented in last generation massively parallel GPU hardware, allowing fast processing of large ligand databases.

Conclusions

BINDSURF is an efficient and fast blind methodology for the determination of protein binding sites depending on the ligand, that uses the massively parallel architecture of GPUs for fast pre-screening of large ligand databases. Its results can also guide posterior application of more detailed VS methods in concrete binding sites of proteins, and its utilization can aid in drug discovery, design, repurposing and therefore help considerably in clinical research.

Relationship between hot spot residues and ligand binding hot spots in protein-protein interfaces

In the context of protein-protein interactions, the term "hot spot" refers to a residue or cluster of residues that makes a major contribution to the binding free energy, as determined by alanine scanning mutagenesis. In contrast, in pharmaceutical research, a hot spot is a site on a target protein that has high propensity for ligand binding and hence is potentially important for drug discovery. Here we examine the relationship between these two hot spot concepts by comparing alanine scanning data for a set of 15 proteins with results from mapping the protein surfaces for sites that can bind fragment-sized small molecules. We find the two types of hot spots are largely complementary; the residues protruding into hot spot regions identified by computational mapping or experimental fragment screening are almost always themselves hot spot residues as defined by alanine scanning experiments. Conversely, a residue that is found by alanine scanning to contribute little to binding rarely interacts with hot spot regions on the partner protein identified by fragment mapping. In spite of the strong correlation between the two hot spot concepts, they fundamentally differ, however. In particular, while identification of a hot spot by alanine scanning establishes the potential to generate substantial interaction energy with a binding partner, there are additional topological requirements to be a hot spot for small molecule binding. Hence, only a minority of hot spots identified by alanine scanning represent sites that are potentially useful for small inhibitor binding, and it is this subset that is identified by experimental or computational fragment screening.

Variable ligand- and receptor-binding hot spots in key strains of influenza neuraminidase

Influenza A continues to be a major public health concern due to its ability to cause epidemic and pandemic disease outbreaks in humans. Computational investigations of structural dynamics of the major influenza glycoproteins, especially the neuraminidase (NA) enzyme, are able to provide key insights beyond what is currently accessible with standard experimental techniques. In particular, all-atom molecular dynamics simulations reveal the varying degrees of flexibility for such enzymes. Here we present an analysis of the relative flexibility of the ligand- and receptor-binding area of three key strains of influenza A: highly pathogenic H5N1, the 2009 pandemic H1N1, and a human N2 strain. Through computational solvent mapping, we investigate the various ligand- and receptor-binding “hot spots” that exist on the surface of NA which interacts with both sialic acid receptors on the host cells and antiviral drugs. This analysis suggests that the variable cavities found in the different strains and their corresponding capacities to bind ligand functional groups may play an important role in the ability of NA to form competent reaction encounter complexes with other species of interest, including antiviral drugs, sialic acid receptors on the host cell surface, and the hemagglutinin protein. Such considerations may be especially useful for the prediction of how such complexes form and with what binding capacity.

Elements of nucleotide specificity in the Trypanosoma brucei mitochondrial RNA editing enzyme RET2

The causative agent of African sleeping sickness, Trypanosoma brucei , undergoes an unusual mitochondrial RNA editing process that is essential for its survival. RNA editing terminal uridylyl transferase 2 of T. brucei (TbRET2) is an indispensable component of the editosome machinery that performs this editing. TbRET2 is required to maintain the vitality of both the insect and bloodstream forms of the parasite, and with its high-resolution crystal structure, it poses as a promising pharmaceutical target. Neither the exclusive requirement of uridine 5'-triphosphate (UTP) for catalysis, nor the RNA primer preference of TbRET2 is well-understood. Using all-atom explicitly solvated molecular dynamics (MD) simulations, we investigated the effect of UTP binding on TbRET2 structure and dynamics, as well as the determinants governing TbRET2's exclusive UTP preference. Through our investigations of various nucleoside triphosphate substrates (NTPs), we show that UTP preorganizes the binding site through an extensive water-mediated H-bonding network, bringing Glu424 and Arg144 side chains to an optimum position for RNA primer binding. In contrast, cytosine 5'-triphosphate (CTP) and adenosine 5'-triphosphate (ATP) cannot achieve this preorganization and thus preclude productive RNA primer binding. Additionally, we have located ligand-binding "hot spots" of TbRET2 based on the MD conformational ensembles and computational fragment mapping. TbRET2 reveals different binding pockets in the apo and UTP-bound MD simulations, which could be targeted for inhibitor design.

Hydrophobicity of proteins and interfaces: insights from density fluctuations

Macroscopic characterizations of hydrophobicity (e.g., contact angle measurements) do not extend to the surfaces of proteins and nanoparticles. Molecular measures of hydrophobicity of such surfaces need to account for the behavior of hydration water. Theory and state-of-the-art simulations suggest that water density fluctuations provide such a measure; fluctuations are enhanced near hydrophobic surfaces and quenched with increasing surface hydrophilicity. Fluctuations affect conformational equilibria and dynamics of molecules at interfaces. Enhanced fluctuations are reflected in enhanced cavity formation, more favorable binding of hydrophobic solutes, increased compressibility of hydration water, and enhanced water-water correlations at hydrophobic surfaces. These density fluctuation-based measures can be used to develop practical methods to map the hydrophobicity/philicity of heterogeneous surfaces including those of proteins. They highlight that the hydrophobicity of a group is context dependent and is significantly affected by its environment (e.g., chemistry and topography) and especially by confinement. The ability to include information about hydration water in mapping hydrophobicity is expected to significantly impact our understanding of protein-protein interactions as well as improve drug design and discovery methods and bioseparation processes.

A molecular dynamics ensemble-based approach for the mapping of druggable binding sites

An expanding repertoire of "allosteric" drugs is revealing that structure-based drug design (SBDD) is not restricted to the "active site" of the target protein. Such compounds have been shown to bind distant regions of the protein topography, potentially providing higher levels of target specificity, reduced toxicity and access to new regions of chemical space. Unfortunately, the location of such allosteric pockets is not obvious in the absence of a bound crystal structure and the ability to predict their presence would be useful in the discovery of novel therapies. Here, we describe a method for the prediction of "druggable" binding sites that takes protein flexibility into account through the use of molecular dynamics (MD) simulation. By using a dynamic representation of the target, we are able to sample multiple protein conformations that may expose new drug-binding surfaces. We perform a fragment-based mapping analysis of individual structures in the MD ensemble using the FTMAP algorithm and then rank the most prolific probe-binding protein residues to determine potential "hot-spots" for further examination. This approach has recently been applied to a pair of human G-protein-coupled receptors (GPCRs), resulting in the detection of five potential allosteric sites.

Analysis of protein binding sites by computational solvent mapping

Computational solvent mapping globally samples the surface of target proteins using molecular probes-small molecules or functional groups-to identify potentially favorable binding positions. The method is based on X-ray and NMR screening studies showing that the binding sites of proteins also bind a large variety of fragment-sized molecules. We have developed the multistage mapping algorithm FTMap (available as a server at http://ftmap.bu.edu/ ) based on the fast Fourier transform (FFT) correlation approach. Identifying regions of low free energy rather than individual low energy conformations, FTMap reproduces the available experimental mapping results. Applications to a variety of proteins show that the probes always cluster in important subsites of the binding site, and the amino acid residues that interact with many probes also bind the specific ligands of the protein. The "consensus" sites at which a number of different probes cluster are likely to be "druggable" sites, capable of binding drug-size ligands with high affinity. Due to its sensitivity to conformational changes, the method can also be used for comparing the binding sites in different structures of a protein.

Allosteric regulation of human liver pyruvate kinase by peptides that mimic the phosphorylated/dephosphorylated N-terminus

An advantage of studying allosteric regulation over covalent modification is that allostery allows the experimentalist to vary the concentration of effector, thereby allowing independent quantification of effector binding and allosteric coupling. In turn, this capacity allows the use of effector analogues to determine which regions of the effector contribute to effector binding and which contribute to allosteric regulation. Like many other proteins, human liver pyruvate kinase (hL-PYK) is regulated by phosphorylation. The phosphorylation of hL-PYK occurs on Ser12 of the N-terminus. Phosphorylation appears to interrupt an interaction (distant from the active site) between the N-terminus and the main body of the protein. Since this interaction increases the affinity of hL-PYK for the substrate (phosphoenolpyruvate, PEP), phosphorylation-dependent interruption of the N-terminus/main-body interaction results in an antagonism of PEP binding. Due to the advantages of studying an allosteric system, we detail a protocol to express and purify N-terminal peptides of hL-PYK using a SUMO-fusion system. We further demonstrate that these peptides act as allosteric regulators that modulate the affinity of hL-PYK for PEP.