Binding Hot Spots Research Articles

Structural basis of GM-CSF and IL-2 sequestration by the viral decoy receptor GIF.

Subversion of the host immune system by viruses is often mediated by molecular decoys that sequester host proteins pivotal to mounting effective immune responses. The widespread mammalian pathogen parapox Orf virus deploys GIF, a member of the poxvirus immune evasion superfamily, to antagonize GM-CSF (granulocyte macrophage colony-stimulating factor) and IL-2 (interleukin-2), two pleiotropic cytokines of the mammalian immune system. However, structural and mechanistic insights into the unprecedented functional duality of GIF have remained elusive. Here we reveal that GIF employs a dimeric binding platform that sequesters two copies of its target cytokines with high affinity and slow dissociation kinetics to yield distinct complexes featuring mutually exclusive interaction footprints. We illustrate how GIF serves as a competitive decoy receptor by leveraging binding hotspots underlying the cognate receptor interactions of GM-CSF and IL-2, without sharing any structural similarity with the cytokine receptors. Our findings contribute to the tracing of novel molecular mimicry mechanisms employed by pathogenic viruses.

Identification of an Allosteric Binding Site on Human Lysosomal Alpha-Galactosidase Opens the Way to New Pharmacological Chaperones for Fabry Disease

Personalized therapies are required for Fabry disease due to its large phenotypic spectrum and numerous different genotypes. In principle, missense mutations that do not affect the active site could be rescued with pharmacological chaperones. At present pharmacological chaperones for Fabry disease bind the active site and couple a stabilizing effect, which is required, to an inhibitory effect, which is deleterious. By in silico docking we identified an allosteric hot-spot for ligand binding where a drug-like compound, 2,6-dithiopurine, binds preferentially. 2,6-dithiopurine stabilizes lysosomal alpha-galactosidase in vitro and rescues a mutant that is not responsive to a mono-therapy with previously described pharmacological chaperones, 1-deoxygalactonojirimycin and galactose in a cell based assay.

Global mapping of antibody recognition of the hepatitis C virus E2 glycoprotein: Implications for vaccine design

The E2 envelope glycoprotein is the primary target of human neutralizing antibody response against hepatitis C virus (HCV), and is thus a major focus of vaccine and immunotherapeutics efforts. There is emerging evidence that E2 is a highly complex, dynamic protein with residues across the protein that are modulating antibody recognition, local and global E2 stability, and viral escape. To comprehensively map these determinants, we performed global E2 alanine scanning with a panel of 16 human monoclonal antibodies (hmAbs), resulting in an unprecedented dataset of the effects of individual alanine substitutions across the E2 protein (355 positions) on antibody recognition. Analysis of shared energetic effects across the antibody panel identified networks of E2 residues involved in antibody recognition and local and global E2 stability, as well as predicted contacts between residues across the entire E2 protein. Further analysis of antibody binding hotspot residues defined groups of residues essential for E2 conformation and recognition for all 14 conformationally dependent E2 antibodies and subsets thereof, as well as residues that enhance antibody recognition when mutated to alanine, providing a potential route to engineer E2 vaccine immunogens. By incorporating E2 sequence variability, we found a number of E2 polymorphic sites that are responsible for loss of neutralizing antibody binding. These data and analyses provide fundamental insights into antibody recognition of E2, highlighting the dynamic and complex nature of this viral envelope glycoprotein, and can serve as a reference for development and rational design of E2-targeting vaccines and immunotherapeutics.

Identification of Ligand Binding Hot Spots of the Histamine H1 Receptor following Structure-Based Fragment Optimization.

Developments in G protein-coupled receptor (GPCR) structural biology provide insights into GPCR-ligand binding. Compound 1 (4-(2-benzylphenoxy)piperidine) with high ligand efficiency for the histamine H1 receptor (H1R) was used to design derivatives to investigate the roles of (i) the amine-binding region, (ii) the upper and lower aromatic region, and (iii) binding site solvation. SAR analysis showed that the amine-binding region serves as the primary binding hot spot, preferably binding small tertiary amines. In silico prediction of water network energetics and mutagenesis studies indicated that the displacement of a water molecule from the amine-binding region is most likely responsible for the increased affinity of the N-methylated analog of 1. Deconstruction of 1 showed that the lower aromatic region serves as a secondary binding hot spot. This study demonstrates that an X-ray structure in combination with tool compounds, assessment of water energetics, and mutagenesis studies enables SAR exploration to map GPCR-ligand binding hot spots.

Crystallographic Structure of Truncated CCL21 and the Putative Sulfotyrosine-Binding Site.

CCL21 chemokine binds the G protein-coupled receptor CCR7, aiding not only in immune response but also in cancer metastasis. Compared with other chemokines, CCL21 has a unique extended unstructured C-terminus that is truncated in some naturally occurring variants. We have determined the X-ray crystallographic structure of a truncated CCL21 (residues 1-79) lacking the extended C-terminus and identified, via two-dimensional nuclear magnetic resonance (NMR), a putative sulfotyrosine-binding site that may recognize such post-translationally modified tyrosine residues on the receptor. Compared to the previously determined NMR structure of full-length CCL21, the crystal structure presents new druggable binding hot spots resulting from an alternative N-loop conformation. In addition, whereas the previous NMR structure did not provide any structural information after residue 70, the C-terminus of the truncated CCL21, ordered up to Ala77 in our crystal structure, is placed near the N-loop and sulfotyrosine-binding site, indicating that the extended C-terminus of full-length CCL21 can interact with this important region for receptor binding. These observations suggest a potential origin for the autoinhibition of CCL21 activity that was recently described. The new crystal structure and binding hot spot analysis have important implications for the function of the CCL21 C-terminus and drug discovery.

Docking-undocking combination applied to the D3R Grand Challenge 2015.

Novel methods for drug discovery are constantly under development and independent exercises to test and validate them for different goals are extremely useful. The drug discovery data resource (D3R) Grand Challenge 2015 offers an excellent opportunity as an external assessment and validation experiment for Computer-Aided Drug Discovery methods. The challenge comprises two protein targets and prediction tests: binding mode and ligand ranking. We have faced both of them with the same strategy: pharmacophore-guided docking followed by dynamic undocking (a new method tested experimentally here) and, where possible, critical assessment of the results based on pre-existing information. In spite of using methods that are qualitative in nature, our results for binding mode and ligand ranking were amongst the best on Hsp90. Results for MAP4K4 were less positive and we track the different performance across systems to the level of previous knowledge about accessible conformational states. We conclude that docking is quite effective if supplemented by dynamic undocking and empirical information (e.g. binding hot spots, productive protein conformations). This setup is well suited for virtual screening, a frequent application that was not explicitly tested in this edition of the D3R Grand Challenge 2015. Protein flexibility remains as the main cause for hard failures.

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.

Nucleosome eviction in mitosis assists condensin loading and chromosome condensation.

Condensins associate with DNA and shape mitotic chromosomes. Condensins are enriched nearby highly expressed genes during mitosis, but how this binding is achieved and what features associated with transcription attract condensins remain unclear. Here, we report that condensin accumulates at or in the immediate vicinity of nucleosome-depleted regions during fission yeast mitosis. Two transcriptional coactivators, the Gcn5 histone acetyltransferase and the RSC chromatin-remodelling complex, bind to promoters adjoining condensin-binding sites and locally evict nucleosomes to facilitate condensin binding and allow efficient mitotic chromosome condensation. The function of Gcn5 is closely linked to condensin positioning, since neither the localization of topoisomerase II nor that of the cohesin loader Mis4 is altered in gcn5 mutant cells. We propose that nucleosomes act as a barrier for the initial binding of condensin and that nucleosome-depleted regions formed at highly expressed genes by transcriptional coactivators constitute access points into chromosomes where condensin binds free genomic DNA.

A "hotspot" for autoimmune T cells in type 1 diabetes.

The ability of a single T cell antigen receptor (TCR) to cross-react with multiple antigens allows the finite number of T cells within an organism to respond to the compendium of pathogen challenges faced during a lifetime. Effective immune surveillance, however, comes at a price. TCR cross-reactivity can allow molecular mimics to spuriously activate autoimmune T cells; it also underlies T cell rejection of organ transplants and drives graft-versus-host disease. In this issue of the JCI, Cole and colleagues provide insight into how an insulin-reactive T cell cross-reacts with pathogen-derived antigens by focusing on a limited portion of the peptides to provide a hotspot for binding. These findings dovetail with recent studies of alloreactive and autoimmune TCRs and suggest that the biochemical principles that govern conventional protein-protein interactions may allow the specificity and cross-reactivity profiles of T cells to be predicted.

Identification of Novel Potential gp120 of HIV-1 Antagonist Using Per-Residue Energy Contribution-Based Pharmacophore modelling.

Inhibition of HIV-1 target cell entry, by targeting gp120, has been identified as a promising approach for the identification and development of prophylactic and salvage HIV infection inhibitors. A small molecule compound 18A is an important chemotype in the development of novel and diverse viral cell entry inhibitors, as it inhibits a wide variety of HIV strains by disrupting allosteric structuring on gp120. This study combines residue energy contribution (REC) pharmacophore mapping of 18A and in silico molecular docking in a virtual screening campaign to identify novel and diverse antagonists of gp120. The binding free energy of a validated docked complex of gp120-18A and the quantitative contribution of interacting residues were obtained with a more accurate molecular mechanics/generalised born surface area (MM/GBSA) method followed by mapping the energetically favourable residue contributions onto atom centres in 18A to obtain a pharmacophore model. The generated pharmacophore hypothesis was used to search the ZINC database for 3D structures that match the pharmacophore. Further, molecular docking, molecular dynamics simulations and binding free energy analysis were performed on retrieved hits in order to rank hits based on their affinity and interactions in the CD4 binding cavity of a gp120. Interestingly, the top scoring compound designated with ZINC database ID as ZINC64700951 (docking score=-8.8kcal/mol, ∆G=-43.77kcal/mol) showed higher affinity compared to compound 18A docking score=-7.3kcal/mol, ∆G=-31.97kcal/mol) and interaction of ZN64700951 with validated allosteric hot spot residues, Asp368 and Met426, and binding hot spot residues, Asn425, Glu370, Gly473, Trp427 and Met475 in gp120, suggest that ZN64700951 is a promising antagonist of gp120. Thus, ZN64700951 could serve as an additional prototype for further optimisation as an HIV target cell viral entry inhibitor.

Mechanism of pKID/KIX Association Studied by Molecular Dynamics Free Energy Simulations.

The phosphorylated kinase-inducible domain (pKID) associates with the kinase interacting domain (KIX) via a coupled folding and binding mechanism. The pKID domain is intrinsically disordered when unbound and upon phosphorylation at Ser133 binds to the KIX domain adopting a well-defined kinked two-helix structure. In order to identify putative hot spot residues of binding that could serve as an initial stable anchor, we performed in silico alanine scanning free energy simulations. The simulations indicate that charged residues including the phosphorylated central Ser133 of pKID make significant contributions to binding. However, these are of slightly smaller magnitude compared to several hydrophobic side chains not defining a single dominant binding hot spot. Both continuous molecular dynamics (MD) simulations and free energy analysis demonstrate that phosphorylation significantly stabilizes the central kinked motif around Ser133 of pKID and shifts the conformational equilibrium toward the bound conformation already in the absence of KIX. This result supports a view that pKID/KIX association follows in part a conformational selection process. During a 1.5 μs explicit solvent MD simulation, folding of pKID on the surface of KIX was observed after an initial contact at the bound position of the phosphorylation site was enforced following a sequential process of αA helix association and a stepwise association and folding of the second αB helix compatible with available experimental results.

Structure-Based Identification of Novel Ligands Targeting Multiple Sites within a Chemokine-G-Protein-Coupled-Receptor Interface.

CXCL12 is a human chemokine that recognizes the CXCR4 receptor and is involved in immune responses and metastatic cancer. Interactions between CXCL12 and CXCR4 are an important drug target but, like other elongated protein-protein interfaces, present challenges for small molecule ligand discovery due to the relatively shallow and featureless binding surfaces. Calculations using an NMR complex structure revealed a binding hot spot on CXCL12 that normally interacts with the I4/I6 residues from CXCR4. Virtual screening was performed against the NMR model, and subsequent testing has verified the specific binding of multiple docking hits to this site. Together with our previous results targeting two other binding pockets that recognize sulfotyrosine residues (sY12 and sY21) of CXCR4, including a new analog against the sY12 binding site reported herein, we demonstrate that protein-protein interfaces can often possess multiple sites for engineering specific small molecule ligands that provide lead compounds for subsequent optimization by fragment based approaches.

Fragment-Based Approaches to the Development of Mycobacterium tuberculosis CYP121 Inhibitors.

The essential enzyme CYP121 is a target for drug development against antibiotic resistant strains of Mycobacterium tuberculosis. A triazol-1-yl phenol fragment 1 was identified to bind to CYP121 using a cascade of biophysical assays. Synthetic merging and optimization of 1 produced a 100-fold improvement in binding affinity, yielding lead compound 2 (KD = 15 μM). Deconstruction of 2 into its component retrofragments allowed the group efficiency of structural motifs to be assessed, the identification of more LE scaffolds for optimization and highlighted binding affinity hotspots. Structure-guided addition of a metal-binding pharmacophore onto LE retrofragment scaffolds produced low nanomolar (KD = 15 nM) CYP121 ligands. Elaboration of these compounds to target binding hotspots in the distal active site afforded compounds with excellent selectivity against human drug-metabolizing P450s. Analysis of the factors governing ligand potency and selectivity using X-ray crystallography, UV–vis spectroscopy, and native mass spectrometry provides insight for subsequent drug development.

Binding mode of chitin and TLR2 via molecular docking and dynamics simulation

Innate immunity is an important part of immune system, providing immediate defence for the host against various infections through phagocytes. Toll-like receptors (TLRs) are major proteins expressed on the cell membrane known as pattern recognition receptors (PRR) that recognise non-self molecules (pathogen-associated molecular patterns (PAMPs)). Because TLRs have been implicated in many inflammatory diseases and cancer, TLRs targeted therapeutics have drawn great attention in clinical application in wide range of conditions. Many of them are undergoing evaluation in clinical trials. Chitin is the second most abundant polysaccharide detected in many insects and fungi. Studies have shown that chitin, as major PAMPs in host-infection, can activate TLR2-dependent innate immunity pathway. Therefore, chitin has potential use as an important agonist or antagonist to control key processes in innate immunity. However, no direct evidence has shown that chitin is the direct target of TLR2. This study first demonstrates a binding model of chitin and TLR2 and then confirmed its stability by molecular dynamic simulation and MM/PBSA (molecular mechanics/Poisson−Boltzmann surface area) calculations. The binding between chitin and TLR2 was taken place inside the binding pocket. Two hydrogen bonds were formed between chitin and TLR2, including Ser320 and Lys321. The van der Waals interaction has the major contribution in stabilising the binding of the chitin molecule with the protein. This study also suggests six hot-spots for specific binding of chitin in the binding site of TLR2, namely, Phe296, Phe299, Leu302, Thr309, Ser320 and Val322. Molecular dynamics simulation demonstrates that the complex of chitin and TLR2 is very stable with a total binding affinity of −27.2 kcal/mol from MM/PBSA calculation.

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.

The 340-cavity in neuraminidase provides new opportunities for influenza drug development: A molecular dynamics simulation study

Influenza neuraminidase (NA) is a pivotal target for viral infection control. However, the accumulating of mutations compromise the efficacy of NA inhibitors. Thus, it is critical to design new drugs targeted to different motifs of NA. Recently, a new motif called 340-cavity was discovered in NA subtypes close to the calcium binding site. The presence of calcium is known to influence NA activity and thermostability. Therefore, the 340-cavity is a putative ligand-binding site for affecting the normal function of NA. In this study, we performed molecular dynamics simulations of different NA subtypes to explore the mechanism of 340-loop formation. Ligand-binding site prediction and fragment library screening were also carried out to provide evidence for the 340-cavity as a druggable pocket. We found that residues G342 and P/R344 in the 340-loop determine the size of the 340-cavity, and the calcium ion plays an important role in maintaining the conformation of the 340-loop through contacts with G345 and Q347. In addition, the 340-cavity is predicted to be a ligand-binding site by metaPocket, and a sequence analysis method is proposed to predict the existence of the 340-cavity. Our study shows that the 340-cavity is not an occasional or atypical domain in NA subtypes, and it has potential to function as a new hotspot for influenza drug binding.

Co-Occurring Atomic Contacts for the Characterization of Protein Binding Hot Spots.

A binding hot spot is a small area at a protein-protein interface that can make significant contribution to binding free energy. This work investigates the substantial contribution made by some special co-occurring atomic contacts at a binding hot spot. A co-occurring atomic contact is a pair of atomic contacts that are close to each other with no more than three covalent-bond steps. We found that two kinds of co-occurring atomic contacts can play an important part in the accurate prediction of binding hot spot residues. One is the co-occurrence of two nearby hydrogen bonds. For example, mutations of any residue in a hydrogen bond network consisting of multiple co-occurring hydrogen bonds could disrupt the interaction considerably. The other kind of co-occurring atomic contact is the co-occurrence of a hydrophobic carbon contact and a contact between a hydrophobic carbon atom and a π ring. In fact, this co-occurrence signifies the collective effect of hydrophobic contacts. We also found that the B-factor measurements of several specific groups of amino acids are useful for the prediction of hot spots. Taking the B-factor, individual atomic contacts and the co-occurring contacts as features, we developed a new prediction method and thoroughly assessed its performance via cross-validation and independent dataset test. The results show that our method achieves higher prediction performance than well-known methods such as Robetta, FoldX and Hotpoint. We conclude that these contact descriptors, in particular the novel co-occurring atomic contacts, can be used to facilitate accurate and interpretable characterization of protein binding hot spots.

Identifying binding hot spots on protein surfaces by mixed-solvent molecular dynamics: HIV-1 protease as a test case.

Mixed-solvent molecular dynamics (MixMD) simulations use full protein flexibility and competition between water and small organic probes to achieve accurate hot-spot mapping on protein surfaces. In this study, we improved MixMD using human immunodeficiency virus type-1 protease (HIVp) as the test case. We used three probe-water solutions (acetonitrile-water, isopropanol-water, and pyrimidine-water), first at 50% w/w concentration and later at 5% v/v. Paradoxically, better mapping was achieved by using fewer probes; 5% simulations gave a superior signal-to-noise ratio and far fewer spurious hot spots than 50% MixMD. Furthermore, very intense and well-defined probe occupancies were observed in the catalytic site and potential allosteric sites that have been confirmed experimentally. The Eye site, an allosteric site underneath the flap of HIVp, has been confirmed by the presence of a 5-nitroindole fragment in a crystal structure. MixMD also mapped two additional hot spots: the Exo site (between the Gly16-Gly17 and Cys67-Gly68 loops) and the Face site (between Glu21-Ala22 and Val84-Ile85 loops). The Exo site was observed to overlap with crystallographic additives such as acetate and dimethyl sulfoxide that are present in different crystal forms of the protein. Analysis of crystal structures of HIVp in different symmetry groups has shown that some surface sites are common interfaces for crystal contacts, which means that they are surfaces that are relatively easy to desolvate and complement with organic molecules. MixMD should identify these sites; in fact, their occupancy values help establish a solid cut-off where "druggable" sites are required to have higher occupancies than the crystal-packing faces.

Molecular dynamics simulation of the interactions between EHD1 EH domain and multiple peptides.

To provide essential information for peptide inhibitor design, the interactions of Eps15 homology domain of Eps15 homology domain-containing protein 1 (EHD1 EH domain) with three peptides containing NPF (asparagine-proline-phenylalanine), DPF (aspartic acid-proline-phenylalanine), and GPF (glycine-proline-phenylalanine) motifs were deciphered at the atomic level. The binding affinities and the underlying structure basis were investigated. Molecular dynamics (MD) simulations were performed on EHD1 EH domain/peptide complexes for 60 ns using the GROMACS package. The binding free energies were calculated and decomposed by molecular mechanics/generalized Born surface area (MM/GBSA) method using the AMBER package. The alanine scanning was performed to evaluate the binding hot spot residues using FoldX software. The different binding affinities for the three peptides were affected dominantly by van der Waals interactions. Intermolecular hydrogen bonds provide the structural basis of contributions of van der Waals interactions of the flanking residues to the binding. van der Waals interactions should be the main consideration when we design peptide inhibitors of EHD1 EH domain with high affinities. The ability to form intermolecular hydrogen bonds with protein residues can be used as the factor for choosing the flanking residues.

Molecular basis of death effector domain chain assembly and its role in caspase-8 activation.

Assembly of a death-inducing signaling complex is a key event in the extrinsic apoptotic pathway, enabling activation of the caspase cascade and subsequent cell death. However, the molecular events governing DISC assembly have remained largely elusive because of the lack of information on mechanism and specificity regulating the death effector domain (DED)-DED interaction network. Using molecular modeling, mutagenesis, and biochemical and ex vivo experiments, we identified the precise binding interface and hot spots crucial for intermolecular DED chain assembly. Mutation of key interface residues (Leu42/Phe45) in procaspase-8 DED-A completely abrogated DED chain formation in HEK293 cells and prevented its association with FADD. A significant 2.6-3.6-fold reduction in procaspase-8 activation was observed in functional cell-death assays after substitution of the interfacial residues. Based on our results we propose a new model for DISC formation that refines the current understanding of the activation mechanism. Upon stimulation, FADD self-associates weakly via reciprocal interaction between helices α1/α4 and α2/α3 of the DED to form an oligomeric signaling platform that provides a stage for the initial recruitment of procaspase-8 through direct interaction with α1/α4 of DED-A, followed by sequential interaction mediated by helices α2/α5 of DED-B, to form the procaspase-8 DED chain that is crucial for its activation and subsequent cell death.