Assessing the structural conservation of protein pockets to study functional and allosteric sites: implications for drug discovery

Background

Nonribosomal peptide synthetases (NRPSs) are multienzymatic, multidomain megasynthases involved in the biosynthesis of pharmaceutically important nonribosomal peptides. The peptaibol synthetase from Trichoderma virens (TPS) is an important member of the NRPS family that exhibits antifungal properties. The majority of the NRPSs terminate peptide synthesis with the thioesterase (TE) domain, which either hydrolyzes the thioester linkage, releasing the free peptic acid, or catalyzes the intramolecular macrocyclization to produce a macrolactone product. TPS is an important NRPS that does not encompass a TE domain, but rather a reductase domain (R domain) to release the mature peptide product reductively with the aid of a NADPH cofactor. However, the catalytic mechanism of the reductase domain has not yet been elucidated.

Results

We present here a three-dimensional (3D) model of the reductase domain based on the crystal structure of vestitone reductase (VR). VR belongs to the short-chain dehydrogenase/reductase (SDR) superfamily and is responsible for the nicotinamide dinucleotide phosphate (NADPH)-dependent reduction of the substrate into its corresponding secondary alcohol product. The binding sites of the probable linear substrates, alamethicin, trichotoxin, antiamoebin I, chrysopermin C and gramicidin, were identified within the modeled R domain using multiple docking approaches. The docking results of the ligand in the active site of the R domain showed that reductase side chains have a high affinity towards ligand binding, while the thioester oxygen of each substrate forms a hydrogen bond with the OH group of Tyr176 and the thiol group of the substrate is closer to the Glu220. The modeling and docking studies revealed the reaction mechanism of reduction of thioester into a primary alcohol.

Conclusion

Peptaibol biosynthesis incorporates a single R domain, which appears to catalyze the four-electron reduction reaction of a peptidyl carrier protein (PCP)-bound peptide to its corresponding primary alcohol. Analysis of R domains present in the non-redundant (nr) database of the NCBI showed that the R domain always resides in the last NRPS module and is involved in either a two or four-electron reduction reaction.     

Are different stoichiometries feasible for complexes between lymphotoxin-alpha and tumor necrosis factor receptor 1?

Background

Tumor necrosis factors, TNF and lymphotoxin-?? (LT), are cytokines that bind to two receptors, TNFR1 and TNFR2 (TNF-receptor 1 and 2) to trigger their signaling cascades. The exact mechanism of ligand-induced receptor activation is still unclear. It is generally assumed that three receptors bind to the homotrimeric ligand to trigger a signaling event. Recent evidence, though, has raised doubts if the ligand:receptor stoichiometry should indeed be 3:3 for ligand-induced cellular response. We used molecular dynamics simulations, elastic network models, as well as MM/PBSA to analyze this question.

Results

Applying MM/PBSA methodology to different stoichiometric complexes of human LT-(TNFR1)_n=1,2,3 the free energy of binding in these complexes has been estimated by single-trajectory and separate-trajectory methods. Simulation studies rationalized the favorable binding energy in the LT-(TNFR1)₁ complex, as evaluated from single-trajectory analysis to be an outcome of the interaction of cysteine-rich domain 4 (CRD4) and the ligand. Elastic network models (ENMs) help to associate the difference in the global fluctuation of the receptors in these complexes. Functionally relevant transformation associated with these complexes reveal the difference in the dynamics of the receptor when free and in complex with LT.

Conclusions

MM/PBSA predicts complexes with a ligand-receptor molar ratio of 3:1 and 3:2 to be energetically favorable. The high affinity associated with LT-(TNFR1)₁ is due to the interaction between the CRD4 domain with LT. The global dynamics ascertained from ENMs have highlighted the differential dynamics of the receptor in different states.     

The DynDom database of protein domain motions

A relational database has been developed based on the results from the application of the DynDom program to a number of proteins for which multiple X-ray conformers are available. The database is populated via a web-based tool that allows visitors to the website to run the DynDom program server-side by selecting pairs of X-ray conformers by Protein Data Bank code and chain identifier. AVAILABILITY: The website can be found at: http://www.sys.uea.ac.uk/dyndom.     

Ligand-induced conformational changes in a thermophilic ribose-binding protein

Background

Members of the periplasmic binding protein (PBP) superfamily are involved in transport and signaling processes in both prokaryotes and eukaryotes. Biological responses are typically mediated by ligand-induced conformational changes in which the binding event is coupled to a hinge-bending motion that brings together two domains in a closed form. In all PBP-mediated biological processes, downstream partners recognize the closed form of the protein. This motion has also been exploited in protein engineering experiments to construct biosensors that transduce ligand binding to a variety of physical signals. Understanding the mechanistic details of PBP conformational changes, both global (hinge bending, twisting, shear movements) and local (rotamer changes, backbone motion), therefore is not only important for understanding their biological function but also for protein engineering experiments.

Results

Here we present biochemical characterization and crystal structure determination of the periplasmic ribose-binding protein (RBP) from the hyperthermophile Thermotoga maritima in its ribose-bound and unliganded state. The T. maritima RBP (tmRBP) has 39% sequence identity and is considerably more resistant to thermal denaturation ( ^app T _m value is 108°C) than the mesophilic Escherichia coli homolog (ecRBP) ( ^app T _m value is 56°C). Polar ligand interactions and ligand-induced global conformational changes are conserved among ecRBP and tmRBP; however local structural rearrangements involving side-chain motions in the ligand-binding site are not conserved.

Conclusion

Although the large-scale ligand-induced changes are mediated through similar regions, and are produced by similar backbone movements in tmRBP and ecRBP, the small-scale ligand-induced structural rearrangements differentiate the mesophile and thermophile. This suggests there are mechanistic differences in the manner by which these two proteins bind their ligands and are an example of how two structurally similar proteins utilize different mechanisms to form a ligand-bound state.     

Structure-Function Analysis of Diacylglycerol Acyltransferase Sequences from 70 Organisms

Background

Diacylglycerol acyltransferase families (DGATs) catalyze the final and rate-limiting step of triacylglycerol (TAG) biosynthesis in eukaryotic organisms. Understanding the roles of DGATs will help to create transgenic plants with value-added properties and provide clues for therapeutic intervention for obesity and related diseases. The objective of this analysis was to identify conserved sequence motifs and amino acid residues for better understanding of the structure-function relationship of these important enzymes.

Results

117 DGAT sequences from 70 organisms including plants, animals, fungi and human are obtained from database search using tung tree DGATs. Phylogenetic analysis separates these proteins into DGAT1 and DGAT2 subfamilies. These DGATs are integral membrane proteins with more than 40% of the total amino acid residues being hydrophobic. They have similar properties and amino acid composition except that DGAT1s are approximately 20 kDa larger than DGAT2s. DGAT1s and DGAT2s have 41 and 16 completely conserved amino acid residues, respectively, although only two of them are shared by all DGATs. These residues are distributed in 7 and 6 sequence blocks for DGAT1s and DGAT2s, respectively, and located at the carboxyl termini, suggesting the location of the catalytic domains. These conserved sequence blocks do not contain the putative neutral lipid-binding domain, mitochondrial targeting signal, or ER retrieval motif. The importance of conserved residues has been demonstrated by site-directed and natural mutants.

Conclusions

This study has identified conserved sequence motifs and amino acid residues in all 117 DGATs and the two subfamilies. None of the completely conserved residues in DGAT1s and DGAT2s is present in recently reported isoforms in the multiple sequences alignment, raising an important question how proteins with completely different amino acid sequences could perform the same biochemical reaction. The sequence analysis should facilitate studying the structure-function relationship of DGATs with the ultimate goal to identify critical amino acid residues for engineering superb enzymes in metabolic engineering and selecting enzyme inhibitors in therapeutic application for obesity and related diseases.     

In silico characterization of immunogenic epitopes presented by HLA-Cw*0401

Background

Arenaviruses are a family of rodent-borne viruses that cause several hemorrhagic fevers. These diseases can be devastating and are often lethal. Herein, to aid in the design and development of diagnostics, treatments and vaccines for arenavirus infections, we have developed a database containing protein sequences from the seven pathogenic arenaviruses (Junin, Guanarito, Sabia, Machupo, Whitewater Arroyo, Lassa and LCMV).

Results

The database currently contains a non-redundant set of 333 protein sequences which were manually annotated. All entries were linked to NCBI and cited PubMed references. The database has a convenient query interface including BLAST search. Sequence variability analyses were also performed and the results are hosted in the database.

Conclusion

The database is available at http://epitope.liai.org:8080/projects/arena and can be used to aid in studies that require proteomic information from pathogenic arenaviruses.     

Absence of carious lesions at margins of glass-ionomer cement and amalgam restorations: An update of systematic review evidence

Background

A recent analysis of protein sequences deposited in the NCBI RefSeq database indicates that ~8.5 million protein sequences are encoded in prokaryotic and eukaryotic genomes, where ~30% are explicitly annotated as "hypothetical" or "uncharacterized" protein. Our Comparison of Protein Active-Site Structures (CPASS v.2) database and software compares the sequence and structural characteristics of experimentally determined ligand binding sites to infer a functional relationship in the absence of global sequence or structure similarity. CPASS is an important component of our Functional Annotation Screening Technology by NMR (FAST-NMR) protocol and has been successfully applied to aid the annotation of a number of proteins of unknown function.

Findings

We report a major upgrade to our CPASS software and database that significantly improves its broad utility. CPASS v.2 is designed with a layered architecture to increase flexibility and portability that also enables job distribution over the Open Science Grid (OSG) to increase speed. Similarly, the CPASS interface was enhanced to provide more user flexibility in submitting a CPASS query. CPASS v.2 now allows for both automatic and manual definition of ligand-binding sites and permits pair-wise, one versus all, one versus list, or list versus list comparisons. Solvent accessible surface area, ligand root-mean square difference, and Cβ distances have been incorporated into the CPASS similarity function to improve the quality of the results. The CPASS database has also been updated.

Conclusions

CPASS v.2 is more than an order of magnitude faster than the original implementation, and allows for multiple simultaneous job submissions. Similarly, the CPASS database of ligand-defined binding sites has increased in size by ~ 38%, dramatically increasing the likelihood of a positive search result. The modification to the CPASS similarity function is effective in reducing CPASS similarity scores for false positives by ~30%, while leaving true positives unaffected. Importantly, receiver operating characteristics (ROC) curves demonstrate the high correlation between CPASS similarity scores and an accurate functional assignment. As indicated by distribution curves, scores ≥ 30% infer a functional similarity. Software URL: http://cpass.unl.edu.     

Novel and unexpected bacterial diversity in an arsenic-rich ecosystem revealed by culture-dependent approaches

Background

BLAST is a commonly-used software package for comparing a query sequence to a database of known sequences; in this study, we focus on protein sequences. Position-specific-iterated BLAST (PSI-BLAST) iteratively searches a protein sequence database, using the matches in round i to construct a position-specific score matrix (PSSM) for searching the database in round i?+?1. Biegert and S?ding developed Context-sensitive BLAST (CS-BLAST), which combines information from searching the sequence database with information derived from a library of short protein profiles to achieve better homology detection than PSI-BLAST, which builds its PSSMs from scratch.

Results

We describe a new method, called domain enhanced lookup time accelerated BLAST (DELTA-BLAST), which searches a database of pre-constructed PSSMs before searching a protein-sequence database, to yield better homology detection. For its PSSMs, DELTA-BLAST employs a subset of NCBI??s Conserved Domain Database (CDD). On a test set derived from ASTRAL, with one round of searching, DELTA-BLAST achieves a ROC₅₀₀₀ of 0.270 vs. 0.116 for CS-BLAST. The performance advantage diminishes in iterated searches, but DELTA-BLAST continues to achieve better ROC scores than CS-BLAST.

Conclusions

DELTA-BLAST is a useful program for the detection of remote protein homologs. It is available under the ??Protein BLAST?? link at http://blast.ncbi.nlm.nih.gov.

Reviewers

This article was reviewed by Arcady Mushegian, Nick V. Grishin, and Frank Eisenhaber.     

Genes encoding hub and bottleneck enzymes of the Arabidopsis metabolic network preferentially retain homeologs through whole genome duplication

Background

Whole genome duplication (WGD) occurs widely in angiosperm evolution. It raises the intriguing question of how interacting networks of genes cope with this dramatic evolutionary event.

Results

In study of the Arabidopsis metabolic network, we assigned each enzyme (node) with topological centralities (in-degree, out-degree and between-ness) to measure quantitatively their centralities in the network. The Arabidopsis metabolic network is highly modular and separated into 11 interconnected modules, which correspond well to the functional metabolic pathways. The enzymes with higher in-out degree and between-ness (defined as hub and bottleneck enzymes, respectively) tend to be more conserved and preferentially retain homeologs after WGD. Moreover, the simultaneous retention of homeologs encoding enzymes which catalyze consecutive steps in a pathway is highly favored and easily achieved, and enzyme-enzyme interactions contribute to the retention of one-third of WGD enzymes.

Conclusions

Our analyses indicate that the hub and bottleneck enzymes of metabolic network obtain great benefits from WGD, and this event grants clear evolutionary advantages in adaptation to different environments.     

Prevalence and risk factors of vertebral compression fractures in female SLE patients

Introduction

The present study objective was to evaluate the incidence of methotrexate (MTX)-specific liver lesions from the analysis of a liver biopsy of inflammatory arthritis patients with elevated liver enzymes.

Methods

A case-control study was performed with 1,571 arthritis patients on long-term low-dose MTX therapy. Results of liver biopsy were analyzed in 41 patients with elevated liver enzymes. The expression of autoimmune markers was also assessed. This population was compared with 41 disease control subjects obtained from the same database, also on MTX but without elevated liver enzymes, matched for age, sex and rheumatic disease.

Results

Compared with the disease controls, patients with liver biopsy showed lower disease duration and lower MTX exposure, weekly and cumulative doses, reflecting shorter treatment duration due to liver abnormalities. Liver biopsies showed 17 autoimmune hepatitis-like (AIH-like) lesions, 13 nonalcoholic steatohepatitis-like lesions, seven limited liver lesions, and two primary biliary cirrhoses. However, MTX-specific lesions with dystrophic nuclei in hepatocytes were seen in only two cases. Liver biopsy lesions were associated with autoimmune markers (P = 0.007); notably, AIH-like lesions were associated with rheumatoid arthritis and with the presence of the HLA-DR shared epitope.

Conclusions

MTX-specific liver lesions are rarely observed in arthritis patients under long-term MTX therapy and elevated liver enzymes.     

The impact of either 4-R-hydroxyproline or 4-R-fluoroproline on the conformation and SH3m-cort binding of HPK1 proline-rich peptide

SH3 domains are probably the most abundant molecular-recognition modules of the proteome. A common feature of these domains is their interaction with ligand proteins containing Pro-rich sequences. Crystal and NMR structures of SH3 domains complexes with Pro-rich peptides show that the peptide ligands are bound over a range of up to seven residues in a PPII helix conformation. Short proline-rich peptides usually adopt little or no ordered secondary structure before binding interactions, and consequently their association with the SH3 domain is characterized by unfavorable binding entropy due to a loss of rotational freedom on forming the PPII helix. With the aim to stabilize the PPII helix conformation into the proline-rich decapeptide PPPLPPKPKF (P2), we replaced some proline residues either with the 4(R)-4-fluoro-l-proline (FPro) or the 4(R)-4-hydroxy-l-proline (Hyp). The interactions of P2 analogues with the SH3 domain of cortactin (SH3_m-cort) were analyzed by circular dichroism spectroscopy, while CD thermal transition experiments have been used to determine their propensity to adopt a PPII helix conformation. Results show that the introduction of three residues of Hyp efficiently stabilizes the PPII helix conformation, while it does not improve the affinity towards the SH3 domain, suggesting that additional forces, e.g., electrostatic interactions, are involved in the SH3_m-cort substrate recognition.     

Multiple sclerosis and pregnancy: what does the patient think? a questionnaire study

Background

The first objective of a DNA microarray experiment is typically to generate a list of genes or probes that are found to be differentially expressed or represented (in the case of comparative genomic hybridizations and/or copy number variation) between two conditions or strains. Rank Products analysis comprises a robust algorithm for deriving such lists from microarray experiments that comprise small numbers of replicates, for example, less than the number required for the commonly used t-test. Currently, users wishing to apply Rank Products analysis to their own microarray data sets have been restricted to the use of command line-based software which can limit its usage within the biological community.

Findings

Here we have developed a web interface to existing Rank Products analysis tools allowing users to quickly process their data in an intuitive and step-wise manner to obtain the respective Rank Product or Rank Sum, probability of false prediction and p-values in a downloadable file.

Conclusions

The online interactive Rank Products analysis tool RankProdIt, for analysis of any data set containing measurements for multiple replicated conditions, is available at: http://strep-microarray.sbs.surrey.ac.uk/RankProducts     

Identification of conformational B-cell Epitopes in an antigen from its primary sequence

Background

One of the major challenges in the field of vaccine design is to predict conformational B-cell epitopes in an antigen. In the past, several methods have been developed for predicting conformational B-cell epitopes in an antigen from its tertiary structure. This is the first attempt in this area to predict conformational B-cell epitope in an antigen from its amino acid sequence.

Results

All Support vector machine (SVM) models were trained and tested on 187 non-redundant protein chains consisting of 2261 antibody interacting residues of B-cell epitopes. Models have been developed using binary profile of pattern (BPP) and physiochemical profile of patterns (PPP) and achieved a maximum MCC of 0.22 and 0.17 respectively. In this study, for the first time SVM model has been developed using composition profile of patterns (CPP) and achieved a maximum MCC of 0.73 with accuracy 86.59%. We compare our CPP based model with existing structure based methods and observed that our sequence based model is as good as structure based methods.

Conclusion

This study demonstrates that prediction of conformational B-cell epitope in an antigen is possible from is primary sequence. This study will be very useful in predicting conformational B-cell epitopes in antigens whose tertiary structures are not available. A web server CBTOPE has been developed for predicting B-cell epitope http://www.imtech.res.in/raghava/cbtope/.     

Brain-derived neurotrophic factor receptor TrkB exists as a preformed dimer in living cells

Background

Neurotrophins (NTs) and their receptors play crucial roles in the development, functions and maintenance of nervous systems. It is widely believed that NT-induced dimerization of the receptors initiates the transmembrane signaling. However, it is still controversial whether the receptor molecule has a monomeric or dimeric structure on the cell surface before its ligand binding.

Findings

Using chemical cross-linking, bimolecular fluorescence complementation (BiFC) and luciferase fragment complementation (LFC) assays, in this study, we show the brain-derived neurotrophic factor (BDNF) receptor TrkB exists as a homodimer before ligand binding. We have also found by using BiFC and LFC that the dimer forms in the endoplasmic reticulum (ER), and that the receptor lacking its intracellular domain cannot form the dimeric structure.

Conclusions

Most, if not all, of the TrkB receptor has a preformed, yet inactive, homodimeric structure before BDNF binding. The intracellular domain of TrkB plays a crucial role in the spontaneous dimerization of the newly synthesized receptors, which occurs in ER. These findings provide new insight into an understanding of a molecular mechanism underlying transmembrane signaling mediated by NT receptors.     

The crystal structure of the catalytic domain of a eukaryotic guanylate cyclase

Background

Soluble guanylate cyclases generate cyclic GMP when bound to nitric oxide, thereby linking nitric oxide levels to the control of processes such as vascular homeostasis and neurotransmission. The guanylate cyclase catalytic module, for which no structure has been determined at present, is a class III nucleotide cyclase domain that is also found in mammalian membrane-bound guanylate and adenylate cyclases.

Results

We have determined the crystal structure of the catalytic domain of a soluble guanylate cyclase from the green algae Chlamydomonas reinhardtii at 2.55 Å resolution, and show that it is a dimeric molecule.

Conclusion

Comparison of the structure of the guanylate cyclase domain with the known structures of adenylate cyclases confirms the close similarity in architecture between these two enzymes, as expected from their sequence similarity. The comparison also suggests that the crystallized guanylate cyclase is in an inactive conformation, and the structure provides indications as to how activation might occur. We demonstrate that the two active sites in the dimer exhibit positive cooperativity, with a Hill coefficient of ~1.5. Positive cooperativity has also been observed in the homodimeric mammalian membrane-bound guanylate cyclases. The structure described here provides a reliable model for functional analysis of mammalian guanylate cyclases, which are closely related in sequence.     

Apollo2Go: a web service adapter for the Apollo genome viewer to enable distributed genome annotation

Background

Apollo, a genome annotation viewer and editor, has become a widely used genome annotation and visualization tool for distributed genome annotation projects. When using Apollo for annotation, database updates are carried out by uploading intermediate annotation files into the respective database. This non-direct database upload is laborious and evokes problems of data synchronicity.

Results

To overcome these limitations we extended the Apollo data adapter with a generic, configurable web service client that is able to retrieve annotation data in a GAME-XML-formatted string and pass it on to Apollo's internal input routine.

Conclusion

This Apollo web service adapter, Apollo2Go, simplifies the data exchange in distributed projects and aims to render the annotation process more comfortable. The Apollo2Go software is freely available from ftp://ftpmips.gsf.de/plants/apollo_webservice.     

HDAPD: a web tool for searching the disease-associated protein structures

Background

PCR-restriction fragment length polymorphism (RFLP) assay is a cost-effective method for SNP genotyping and mutation detection, but the manual mining for restriction enzyme sites is challenging and cumbersome. Three years after we constructed SNP-RFLPing, a freely accessible database and analysis tool for restriction enzyme mining of SNPs, significant improvements over the 2006 version have been made and incorporated into the latest version, SNP-RFLPing 2.

Results

The primary aim of SNP-RFLPing 2 is to provide comprehensive PCR-RFLP information with multiple functionality about SNPs, such as SNP retrieval to multiple species, different polymorphism types (bi-allelic, tri-allelic, tetra-allelic or indels), gene-centric searching, HapMap tagSNPs, gene ontology-based searching, miRNAs, and SNP500Cancer. The RFLP restriction enzymes and the corresponding PCR primers for the natural and mutagenic types of each SNP are simultaneously analyzed. All the RFLP restriction enzyme prices are also provided to aid selection. Furthermore, the previously encountered updating problems for most SNP related databases are resolved by an on-line retrieval system.

Conclusions

The user interfaces for functional SNP analyses have been substantially improved and integrated. SNP-RFLPing 2 offers a new and user-friendly interface for RFLP genotyping that can be used in association studies and is freely available at http://bio.kuas.edu.tw/snp-rflping2.     

DFT study of the Lewis acid mediated synthesis of 3-acyltetramic acids

The recent emergence and re-emergence of alphaviruses, in particular the chikungunya virus (CHIKV), in numerous countries has invoked a worldwide threat to human health, while simultaneously generating an economic burden on affected countries. There are currently no vaccines or effective drugs available for the treatment of the CHIKV, and with few lead compounds reported, the vital medicinal chemistry is significantly more challenging. This study reports on the discovery of potential inhibitors for the nsP3 macro domain of CHIKV using molecular docking, virtual screening, and molecular dynamics simulations, as well as work done to evaluate and confirm the active site of nsP3. Virtual screening was carried out based on blind docking as well as focused docking, using the database of 1541 compounds from NCI Diversity Set II, to identify hit compounds for nsP3. The top hit compounds were further subjected to molecular dynamic simulations, yielding a greater understanding of the dynamic behavior of nsP3 and its complexes with various ligands, concurrently confirming the outcomes of docking, and establishing in silico lead compounds which target the CHIKV nsP3 enzyme.

Grb7 SH2 domain structure and interactions with a cyclic peptide inhibitor of cancer cell migration and proliferation

Background

Human growth factor receptor bound protein 7 (Grb7) is an adapter protein that mediates the coupling of tyrosine kinases with their downstream signaling pathways. Grb7 is frequently overexpressed in invasive and metastatic human cancers and is implicated in cancer progression via its interaction with the ErbB2 receptor and focal adhesion kinase (FAK) that play critical roles in cell proliferation and migration. It is thus a prime target for the development of novel anti-cancer therapies. Recently, an inhibitory peptide (G7-18NATE) has been developed which binds specifically to the Grb7 SH2 domain and is able to attenuate cancer cell proliferation and migration in various cancer cell lines.

Results

As a first step towards understanding how Grb7 may be inhibited by G7-18NATE, we solved the crystal structure of the Grb7 SH2 domain to 2.1 Å resolution. We describe the details of the peptide binding site underlying target specificity, as well as the dimer interface of Grb 7 SH2. Dimer formation of Grb7 was determined to be in the μM range using analytical ultracentrifugation for both full-length Grb7 and the SH2 domain alone, suggesting the SH2 domain forms the basis of a physiological dimer. ITC measurements of the interaction of the G7-18NATE peptide with the Grb7 SH2 domain revealed that it binds with a binding affinity of K_d = ~35.7 μM and NMR spectroscopy titration experiments revealed that peptide binding causes perturbations to both the ligand binding surface of the Grb7 SH2 domain as well as to the dimer interface, suggesting that dimerisation of Grb7 is impacted on by peptide binding.

Conclusion

Together the data allow us to propose a model of the Grb7 SH2 domain/G7-18NATE interaction and to rationalize the basis for the observed binding specificity and affinity. We propose that the current study will assist with the development of second generation Grb7 SH2 domain inhibitors, potentially leading to novel inhibitors of cancer cell migration and invasion.