New syntheses of tetrazolylmethylphenylalanine and O-malonyltyrosine as pTyr mimetics for the design of STAT3 dimerization inhibitors

Investigation within the pTyr-binding pocket of the STAT3 SH2 domain led us to develop a novel synthesis of two pTyr mimetics, l-tetrazolylmethylphenylalanine (l-Tmp) and l-O-malonyltyrosine (l-OMT), that were next incorporated in a high affinity ligand of STAT3 SH2 domain. Biological evaluation of peptidomimetics on STAT3 dimerization identified l-OMT as the first non-phosphorus pTyr mimetic so far reported against STAT3 SH2 domain, harboring an activity similar to that of the Pmp-containing reference peptidomimetic.     

Induced fit in liver X receptor beta: a molecular dynamics-based investigation

Ligand induced fit phenomenon occurring at the ligand binding domain of the liver X receptor beta (LXRbeta) was investigated by means of molecular dynamics. Reliability of a 4-ns trajectory was tested from two distinct LXRbeta crystal complexes 1PQ6B/GW and 1PQ9B/T09 characterized by an open and a closed state of the pocket, respectively. Crossed complexes 1PQ6B/T09 and 1PQ9B/GW were then submitted to the same molecular dynamic conditions, which were able to recover LXRbeta conformations similar to the original crystallography data. Analysis of "open to closed" and "closed to open" conformational transitions pointed out the dynamic role of critical residues lining the ligand binding pocket involved in the local remodeling upon ligand binding (e.g., Phe271, Phe329, Phe340, Arg319, Glu281). Altogether, the present study indicates that the molecular dynamic protocol is a consistent approach for managing LXRbeta-related induced fit process. This protocol could therefore be used for refining ligand docking solutions of a structure-based design strategy.     

Binding of antifusion peptides with HIVgp41 from molecular dynamics simulations: quantitative correlation with experiment

Peptides based on C-terminal regions of the human immunodeficiency virus (HIV) viral protein gp41 represent an important new class of antiviral therapeutics called peptide fusion inhibitors. In this study, computational methods were used to model the binding of six peptides that contain residues that pack into a conserved hydrophobic pocket on HIVgp41, an attractive target site for the development of small molecule inhibitors. Free energies of binding were computed using molecular mechanics Generalized Born surface area (MM-GBSA) methods from molecular dynamics (MD) simulations, which employed either explicit (TIP3P) or continuum Generalized Born (GB) water models and strong correlations between experimental and computational affinities were obtained in both cases. Energy decomposition of the TIP3P-MD results (r2 = 0.75) reveals that variation in experimental affinity is highly correlated with changes in intermolecular van der Waals energies (deltaE(vdw)) on both a local (residue-based, r2 = 0.94) and global (peptide-based, r2 = 0.84) scale. The results show that differential association of C-peptides with HIVgp41 is driven solely by changes within the conserved pocket supporting the hypothesis that this region is an important drug target site. Such strong agreement with experiment is notable given the large size of the ligands (34 amino-acids) relative to the small range of experimental affinities (2 kcal/mol) and demonstrates good sensitivity of this computational method for simulating peptide fusion inhibitors. Finally, inspection of simulation trajectories identified a highly populated pi-type hydrogen bond, which formed between Gln575 on the receptor and the aromatic ring of peptide ligand Phe631, which could have important implications for drug design.     

Binding hotspots on K‐ras: Consensus ligand binding sites and other reactive regions from probe‐based molecular dynamics analysis

We have used probe‐based molecular dynamics (pMD) simulations to search for interaction hotspots on the surface of the therapeutically highly relevant oncogenic K‐Ras G12D. Combining the probe‐based query with an ensemble‐based pocket identification scheme and an analysis of existing Ras‐ligand complexes, we show that (i) pMD is a robust and cost‐effective strategy for binding site identification, (ii) all four of the previously reported ligand binding sites are suitable for structure‐based ligand design, and (iii) in some cases probe binding and expanded sampling of configurational space enable pocket expansion and increase the likelihood of site identification. Furthermore, by comparing the distribution of hotspots in nonpocket‐like regions with known protein‐ and membrane‐interacting interfaces, we propose that pMD has the potential to predict surface patches responsible for protein‐biomolecule interactions. These observations have important implications for future drug design efforts and will facilitate the search for potential interfaces responsible for the proposed transient oligomerization or interaction of Ras with other biomolecules in the cellular milieu. Proteins 2015; 83:898–909. © 2015 Wiley Periodicals, Inc.     

Revisiting the phosphotyrosine binding pocket of Fyn SH2 domain led to the identification of novel SH2 superbinders.

Protein engineering through directed evolution is an effective way to obtain proteins with novel functions with the potential applications as tools for diagnosis or therapeutics. Many natural proteins have undergone directed evolution in vitro in the test tubes in the laboratories worldwide, resulting in the numerous protein variants with novel or enhanced functions. we constructed here an SH2 variant library by randomizing 8 variable residues in its phosphotyrosine (pTyr) binding pocket. Selection of this library by a pTyr peptide led to the identification of SH2 variants with enhanced affinities measured by EC50. Fluorescent polarization was then applied to quantify the binding affinities of the newly identified SH2 variants. As a result, three SH2 variants, named V3, V13 and V24, have comparable binding affinities with the previously identified SH2 triple‐mutant superbinder. Biolayer Interferometry assay was employed to disclose the kinetics of the binding of these SH2 superbinders to the phosphotyrosine peptide. The results indicated that all the SH2 superbinders have two‐orders increase of the dissociation rate when binding the pTyr peptide while there was no significant change in their associate rates. Intriguingly, though binding the pTyr peptide with comparable affinity with other SH2 superbinders, the V3 does not bind to the sTyr peptide. However, variant V13 and V24 have cross‐reactivity with both pTyr and sTyr peptides. The newly identified superbinders could be utilized as tools for the identification of pTyr‐containing proteins from tissues under different physiological or pathophysiological conditions and may have the potential in the therapeutics.     

Characterizing tyrosine phosphorylation signaling in lung cancer using SH2 profiling

Background

Tyrosine kinases drive the proliferation and survival of many human cancers. Thus profiling the global state of tyrosine phosphorylation of a tumor is likely to provide a wealth of information that can be used to classify tumors for prognosis and prediction. However, the comprehensive analysis of tyrosine phosphorylation of large numbers of human cancer specimens is technically challenging using current methods.

Methodology/Principal Findings

We used a phosphoproteomic method termed SH2 profiling to characterize the global state of phosphotyrosine (pTyr) signaling in human lung cancer cell lines. This method quantifies the phosphorylated binding sites for SH2 domains, which are used by cells to respond to changes in pTyr during signaling. Cells could be grouped based on SH2 binding patterns, with some clusters correlated with EGF receptor (EGFR) or K-RAS mutation status. Binding of specific SH2 domains, most prominently RAS pathway activators Grb2 and ShcA, correlated with EGFR mutation and sensitivity to the EGFR inhibitor erlotinib. SH2 binding patterns also reflected MET activation and could identify cells driven by multiple kinases. The pTyr responses of cells treated with kinase inhibitors provided evidence of distinct mechanisms of inhibition.

Conclusions/Significance

This study illustrates the potential of modular protein domains and their proteomic binding profiles as powerful molecular diagnostic tools for tumor classification and biomarker identification.     

Structure-based design of novel nonpeptide inhibitors of the Src SH2 domain:phosphotyrosine mimetics exploiting multifunctional group replacement chemistry

A series of novel nonpeptide inhibitors of the pp60(c-Src) (Src) SH2 domain is described that exploit multifunctional group replacement of the phenylphosphate moiety of phosphotyrosine (pTyr). Relative to an x-ray structure of citrate complexed to the pTyr binding site of the Src SH2 domain, these nonpeptide ligands illustrate the systematic replacement of the phosphate group by multiple nonhydrolyzable, mono- or dianionic functionalities. Specifically, several phenylalanine (Phe) analogs incorporating key 4' and 3' substituents were synthesized and incorporated into a bicyclic benzamide template previously reported (W. C. Shakespeare et al., Proceedings of the National Academy of Science USA, 2000, Vol. 97, pp. 9373-9378). These pTyr mimetics included 4',3'-diphosphono-Phe (Dpp), 4',3'-dicarboxymethyloxy-Phe (Dcp), and 4'-phosphono-3'-carboxymethyloxy-Phe (Cpp). Noteworthy were nonpeptide inhibitors 8-11 that were 5- to 10-fold more potent than the cognate tetrapeptide ligand Ac-pTyr-Glu-Glu-Ile-NH(2) in binding to the Src SH2 domain.     

Investigation of phosphotyrosine recognition by the SH2 domain of the Src kinase.

The binding of tyrosine phosphorylated targets by SH2 domains is required for propagation of many cellular signals in higher eukaryotes; however, the determinants of phosphotyrosine (pTyr) recognition by SH2 domains are not well understood. In order to identify the attributes of pTyr required for high affinity interaction with SH2 domains, the binding of the SH2 domain of the Src kinase (Src SH2 domain) to a dephosphorylated peptide, a phosphoserine-containing peptide, and the amino acid pTyr was studied using titration calorimetry and compared with the binding of a high affinity tyrosyl phosphopeptide. The dephosphorylated peptide and the phosphoserine containing peptide both bind extremely weakly to the Src SH2 domain (DeltaGo (dephosphorylated)=-3.6 kcal/mol, DeltaGo (phosphoserine) >-3.7 kcal/mol); however, the DeltaGo value of pTyr binding is more favorable (-4.7 kcal/mol, or 50 % of the entire binding free energy of a high affinity tyrosyl phosphopeptide). These results indicate that both the phosphate and the tyrosine ring of the pTyr are critical determinants of high affinity binding. Alanine mutagenesis was also used to evaluate the energetic contribution to binding of ten residues located in the pTyr-binding site. Mutation of the strictly conserved Arg betaB5 resulted in a large increase in DeltaGo (DeltaDeltaGo=3.2 kcal/mol) while elimination of the other examined residues each resulted in a significantly smaller (DeltaDeltaGo<1.4 kcal/mol) reduction in affinity, indicating that Arg betaB5 is the single most important determinant of pTyr recognition. However, mutation of Cys betaC3, a residue unique to the Src SH2 domain, surprisingly increased affinity by eightfold (DeltaDeltaGo=-1.1 kcal/mol). Using a double mutant cycle analysis, it was revealed that residues of the pTyr-binding pocket are not coupled to the peptide residues C-terminal to the pTyr. In addition, comparison of each residue's DeltaDeltaGo value upon mutation with that residue's sequence conservation among SH2 domains revealed only a modest correlation between a residue's energetic contribution to pTyr recognition and its conservation throughout evolution. The results of this investigation highlight the importance of a single critical interaction, the buried ionic bond between the phosphate of the pTyr and Arg betaB5 of the SH2 domain, driving the binding of SH2 domains to tyrosine phosphorylated targets.     

Direct interaction of STAT4 with the IL-12 receptor.

Signal transduction by interleukin-12 (IL-12) requires phosphorylation and activation of STAT4. Direct interaction of the SH2 domain of STAT4 with a phosphotyrosine residue in the IL-12 receptor has been proposed to be required for the subsequent STAT4 phosphorylation. The IL-12 receptor beta2 subunit contains three tyrosine residues in its cytoplasmic domain. To test the hypothesis that one of these tyrosines is involved in binding STAT4, phosphopeptides were synthesized according to the amino acid sequences surrounding each of these tyrosine residues. Only the phosphopeptide containing pTyr800 strongly bound to STAT4 in a cell-free binding assay. When this phosphopeptide was introduced into TALL-104 cells, it blocked IL-12-induced STAT4 phosphorylation by competing with the IL-12 receptor for binding to STAT4. A series of alanine replacements was performed in this phosphopeptide to elucidate which amino acids surrounding the pTyr800 residue are critical for STAT4 binding. To summarize, the site on the IL-12 receptor which binds STAT4 can be described as -T-X-X-G-pY(800)-L-, where the core G-pY(800)-L motif is critical for the binding; the threonine at the pY-4 position has only a minor contribution and X represents amino acids not critical for the binding. These results demonstrate that only a small region of the IL-12 receptor is critically involved in binding STAT4 and suggest the feasibility that small molecule inhibitors could be identified which interfere with IL-12 signal transduction for treatment of autoimmune diseases.     

Ensemble docking of multiple protein structures: considering protein structural variations in molecular docking

One approach to incorporate protein flexibility in molecular docking is the use of an ensemble consisting of multiple protein structures. Sequentially docking each ligand into a large number of protein structures is computationally too expensive to allow large-scale database screening. It is challenging to achieve a good balance between docking accuracy and computational efficiency. In this work, we have developed a fast, novel docking algorithm utilizing multiple protein structures, referred to as ensemble docking, to account for protein structural variations. The algorithm can simultaneously dock a ligand into an ensemble of protein structures and automatically select an optimal protein structure that best fits the ligand by optimizing both ligand coordinates and the conformational variable m, where m represents the m-th structure in the protein ensemble. The docking algorithm was validated on 10 protein ensembles containing 105 crystal structures and 87 ligands in terms of binding mode and energy score predictions. A success rate of 93% was obtained with the criterion of root-mean-square deviation <2.5 A if the top five orientations for each ligand were considered, comparable to that of sequential docking in which scores for individual docking are merged into one list by re-ranking, and significantly better than that of single rigid-receptor docking (75% on average). Similar trends were also observed in binding score predictions and enrichment tests of virtual database screening. The ensemble docking algorithm is computationally efficient, with a computational time comparable to that for docking a ligand into a single protein structure. In contrast, the computational time for the sequential docking method increases linearly with the number of protein structures in the ensemble. The algorithm was further evaluated using a more realistic ensemble in which the corresponding bound protein structures of inhibitors were excluded. The results show that ensemble docking successfully predicts the binding modes of the inhibitors, and discriminates the inhibitors from a set of noninhibitors with similar chemical properties. Although multiple experimental structures were used in the present work, our algorithm can be easily applied to multiple protein conformations generated by computational methods, and helps improve the efficiency of other existing multiple protein structure(MPS)-based methods to accommodate protein flexibility.     

Selective inhibition of Src SH2 by a novel thiol-targeting tricarbonyl-modified inhibitor and mechanistic analysis by (1)H/(13)C NMR spectroscopy.

Detailed analysis of Src SH2 binding by peptides containing a novel tricarbonyl-modified pTyr moiety is described. We envisaged that Src SH2 selectivity might be obtained by exploiting the thiol group of Cys188 present in the pTyr binding pocket of the protein at the betaC3 position. Peptidyl as well as non-peptidyl compounds 1-4 possessing a 4-alpha,beta-diketoester-modified pTyr mimic exhibited micromolar affinity to Src SH2. Furthermore, these tricarbonyl compounds were selective for Src SH2 to the extent they showed no significant affinity for either Cys188Ser or Cys188Ala Src SH2 mutants. Upon closer examination of the binding of these tricarbonyls to Src SH2 using NMR of 13C-labeled compounds (6a, 6b, and 6c), we found that after the initial binding event the molecule disproportionated in a 'retro-Claisen' fashion to provide benzoic acid 16 and, following hydrolysis of the methyl ester 17, the hemiketal adduct of glyoxalic acid 18.     

Structural basis for specificity switching of the Src SH2 domain

The Src SH2 domain binds pYEEI-containing phosphopeptides in an extended conformation with a hydrophobic pocket, which includes ThrEF1, binding Ile(pY +3). Mutating ThrEF1 to tryptophan switches specificity to an Asn(pY +2) requirement, yielding a biological mimic of the Grb2 SH2 domain. Here we show that the Src ThrEF1Trp SH2 domain mutant binds pYVNV phosphopeptides in a beta turn conformation, which, despite differing conformations of the interacting tryptophan, closely resembles the native Grb2/pYVNV cognate peptide binding mode. The ThrEF1Trp substitution therefore switches specificity by physically occluding the pTyr +3 binding pocket and by providing additional interaction surface area for Asn(pY +2). This demonstrates structurally how novel SH2 domain specificities may rapidly evolve through single amino acid substitutions and suggests how new signaling pathways may develop.     

Structure, dynamics, and binding thermodynamics of the v-Src SH2 domain: implications for drug design

SH2 domains provide fundamental recognition sites in tyrosine kinase-mediated signaling pathways which, when aberrant, give rise to disease states such as cancer, diabetes, and immune deficiency. Designing specific inhibitors that target the SH2 domain-binding site, however, have presented a major challenge. Despite well over a decade of intensive research, clinically useful SH2 domain inhibitors have yet to become available. A better understanding of the structural, dynamic, and thermodynamic contributions to ligand binding of individual SH2 domains will provide some insight as to whether inhibitor development is possible. We report the first high resolution solution structure of the apo-v-Src SH2 domain. This is accompanied by the analysis of backbone dynamics and pK(a) values within the apo- and peptide-bound states. Our results indicate that the phosphotyrosine (pY) pocket is tightly structured and hence not adaptable to exogenous ligands. On the other hand, the pocket which accommodates residues proximal and C-terminal of the pY (pY + 3) or so-called specificity determining region, is a large dynamic-binding surface. This appears to allow a high level of promiscuity in binding. Binding of a series of synthetic, phosphotyrosyl, peptidomimetic compounds designed to explore interactions in the pY + 3 pocket further demonstrates the ability of the SH2 domain to accommodate diverse ligands. The thermodynamic parameters of these interactions show dramatic enthalpy/entropy compensation. These data suggest that the v-Src SH2 domain does not have a highly specific secondary-binding site, which clearly presents a major hurdle to design selective inhibitors.     

Small ligands interacting with the phosphotyrosine binding pocket of the Src SH2 protein

Various small fragments bearing phosphate, phosphonate or phosphonic acid moieties have been prepared through parallel synthesis and their binding potencies evaluated on the Src SH2 protein using a BIAcore assay. This provided us insight into the requirement of the Src SH2 pTyr binding pocket and some promising small ligands have been characterised.     

Structure-based design of novel human Pin1 inhibitors (III): Optimizing affinity beyond the phosphate recognition pocket

The design of potent Pin1 inhibitors has been challenging because its active site specifically recognizes a phospho-protein epitope. The de novo design of phosphate-based Pin1 inhibitors focusing on the phosphate recognition pocket and the successful replacement of the phosphate group with a carboxylate have been previously reported. The potency of the carboxylate series is now further improved through structure-based optimization of ligand–protein interactions in the proline binding site which exploits the H-bond interactions necessary for Pin1 catalytic function. Further optimization using a focused library approach led to the discovery of low nanomolar non-phosphate small molecular Pin1 inhibitors. Structural modifications designed to improve cell permeability resulted in Pin1 inhibitors with low micromolar anti-proliferative activities against cancer cells.     

Reconstructing the binding site of factor Xa in trypsin reveals ligand-induced structural plasticity

In order to investigate issues of selectivity and specificity in protein-ligand interactions, we have undertaken the reconstruction of the binding pocket of human factor Xa in the structurally related rat trypsin by site-directed mutagenesis. Three sequential regions (the "99"-, the "175"- and the "190"- loops) were selected as representing the major structural differences between the ligand binding sites of the two enzymes. Wild-type rat trypsin and variants X99rT and X(99/175/190)rT were expressed in yeast, and analysed for their interaction with factor Xa and trypsin inhibitors. For most of the inhibitors studied, progressive loop replacement at the trypsin surface resulted in inhibitory profiles akin to factor Xa. Crystals of the variants were obtained in the presence of benzamidine (3), and could be soaked with the highly specific factor Xa inhibitor (1). Binding of the latter to X99rT results in a series of structural adaptations to the ligand, including the establishment of an "aromatic box" characteristic of factor Xa. In X(99/175/190)rT, introduction of the 175-loop results in a surprising re-orientation of the "intermediate helix", otherwise common to trypsin and factor Xa. The re-orientation is accompanied by an isomerisation of the Cys168-Cys182 disulphide bond, and burial of the critical Phe174 side-chain. In the presence of (1), a major re-organisation of the binding site takes place to yield a geometry identical to that of factor Xa. In all, binding of (1) to trypsin and its variants results in significant structural rearrangements, inducing a binding surface strongly reminiscent of factor Xa, against which the inhibitor was optimised. The structural data reveal a plasticity of the intermediate helix, which has been implicated in the functional cofactor dependency of many trypsin-like serine proteinases. This approach of grafting loops onto scaffolds of known related structures may serve to bridge the gap between structural genomics and drug design.