Conservation analysis and structure prediction of the SH2 family of phosphotyrosine binding domains

Cytoplasmic tail of transmembrane Dscam controls antibacterial responses by regulating cell proliferation-related genes in hemocytes of Chinese mitten crab (Eriocheir sinensis)

In arthropods, the involvement of Dscam (Down syndrome cell adhesion molecule) in innate immunity has been extensively demonstrated. Its cytoplasmic tail contains multiple conserved functional sites, which indicates its involvement in different intracellular signaling pathways. In this study, we focused on the role of the cytoplasmic tail of Dscam in the Chinese mitten crab (Eriocheir sinensis) immune defense. In the group with cytoplasmic tail knockdown (the site was located on constant exons 37 and 38), 3885 differentially expressed genes (DEGs) were identified. The DEGs were enriched in small molecule binding, protein-containing complex binding, and immunity-related pathways. The expression of selected genes were validated using quantitative real-time reverse transcription PCR. We identified key Cell cycle, Janus kinase (JAK)-signal transducer, activator of transcription (STAT) and mitogen-activated protein kinase (MAPK) signaling pathway genes, the results indicated that the cytoplasmic tail of Dscam controls antibacterial responses by regulating cell proliferation-related genes in hemocytes.

Targeted protein posttranslational modifications by chemically induced proximity for cancer therapy

Post-translational modifications (PTMs) regulate all aspects of protein function. Therefore, upstream regulators of PTMs, such as kinases, acetyltransferases, or methyltransferases, are potential therapeutic targets for human diseases, including cancer. To date, multiple inhibitors and/or agonists of these PTM upstream regulators are in clinical use, while others are still in development. However, these upstream regulators control not only the PTMs of disease-related target proteins but also other disease-irrelevant substrate proteins. Thus, nontargeted perturbing activities may introduce unwanted off-target toxicity issues that limit the use of these drugs in successful clinical applications. Therefore, alternative drugs that solely regulate a specific PTM of the disease-relevant protein target may provide a more precise effect in treating disease with relatively low side effects. To this end, chemically induced proximity has recently emerged as a powerful research tool, and several chemical inducers of proximity (CIPs) have been used to target and regulate protein ubiquitination, phosphorylation, acetylation, and glycosylation. These CIPs have a high potential to be translated into clinical drugs and several examples such as PROTACs and MGDs are now in clinical trials. Hence, more CIPs need to be developed to cover all types of PTMs, such as methylation and palmitoylation, thus providing a full spectrum of tools to regulate protein PTM in basic research and also in clinical application for effective cancer treatment.

Recombinant pebulin protein, a type 2 ribosome-inactivating protein isolated from dwarf elder (Sambucus ebulus L.) shows anticancer and antifungal activities in vitro

In this study, encoding sequence of a new type 2 RIP (pebulin) was isolated and cloned from dwarf elder (Sambucus ebulus L.) native to the northern regions of Iran. The nucleotide sequence of pebulin was ligated to the pET-28a(+) expression plasmid and cloned into the E. coli strain BL21 (DE3) in order to express heterologously of recombinant protein. The recombinant pebulin protein was mainly produced in the form of insoluble inclusion bodies probably because to absence of N-glycosylation process in E. coli. Therefore, in order to increase the expression of recombinant protein in soluble form, co-expression of the target protein with the pG-Tf2 chaperone plasmid and incubation of bacterial culture under low temperature were used to enhance solubility and accumulation of recombinant protein. After purification of the recombinant protein using affinity chromatography method, the bioactivity of pebulin was analyzed by hemagglutination, anticancer, and antifungal assays. The results of the hemagglutination assay showed that purified pebulin agglutinated erythrocytes in all human blood groups. In addition, pebulin considerably inhibited the proliferation of cancer cell lines MCF-7 and HT-29 in a time- and dose-dependent manner and indicated remarkably growth-inhibiting effect against the plant pathogenic fungi such as Alternaria solani and Fusarium oxysporum.

Cullin-RING Ligase 5: Functional characterization and its role in human cancers

Cullin-RING ligase 5 (CRL5) is a multi-protein complex and consists of a scaffold protien cullin 5, a RING protein RBX2 (also known as ROC2 or SAG), adaptor proteins Elongin B/C, and a substrate receptor protein SOCS. Through targeting a variety of substrates for proteasomal degradation or modulating various protein-protein interactions, CRL5 is involved in regulation of many biological processes, such as cytokine signal transduction, inflammation, viral infection, and oncogenesis. As many substrates of CRL5 are well-known oncoproteins or tumor suppressors, abnormal regulation of CRL5 is commonly found in human cancers. In this review, we first briefly introduce each of CRL5 components, and then discuss the biological processes regulated by four members of SOCS-box-containing substrate receptor family through substrate degradation. We next describe how CRL5 is hijacked by a variety of viral proteins to degrade host anti-viral proteins, which facilitates virus infection. We further discuss the regulation of CUL5 and its various roles in human cancers, acting as either a tumor suppressor or an oncoprotein in a context-dependent manner. Finally, we propose novel insights for future perspectives on the validation of cullin5 and other CRL5 components as potential targets, and possible targeting strategies to discover CRL5 inhibitors for anti-cancer and anti-virus therapies.

Copy Number Variation of the SHE Gene in Sheep and Its Association with Economic Traits

Simple Summary

Src Homology 2 Domain Containing E (SHE) is a protein coding gene, and member of the SH2 domain-containing family. Sequencing revealed a 2000 bp copy number variation in the SHE gene. There is overlap between SHE copy number variation (CNV) and quantitative trait loci related to milk fat percentage and bone density. A total of 750 sheep, including Chaka sheep (CKS), Hu sheep (HS), Small Tail Han sheep (STHS) and Large Tail Han sheep (LTHS) were available to detect the CNV of the SHE gene and correlate these gene variations with economic traits. The results showed that there were more individuals with SHE copy number loss in CKS and HS than in STHS and LTHS. Association analysis showed that gain and normal copy number types performed better in body length (p < 0.05), circumference of cannon bone (p < 0.05), heart girth (p < 0.05), chest width (p < 0.05) and high at the cross (p < 0.05) in CKS, HS and STHS. Chi-square analyses found significant variation in the CNV of the SHE gene, so it varies greatly between varieties. These findings clarified the relationship between the CNV of the SHE gene and the economic traits in these four kinds of sheep, and provide a reference for sheep breeding.

Abstract

Copy number variation (CNV) caused by gene rearrangement is an important part of genomic structural variation. We found that the copy number variation region of the Src Homology 2 Domain Containing E (SHE) gene correlates with a quantitative trait locus of sheep related to milk fat percentage and bone density. The aim of our study was to detect the copy number variation of the SHE gene in four sheep breeds and to conduct a correlation analysis with economic traits, hoping to provide some reference for sheep breeding. In this study, we examined 750 sheep from four Chinese breeds: Chaka sheep (CKS), Hu sheep (HS), Large Tail Han sheep (LTHS) and Small Tail Han sheep (STHS). We used qPCR to evaluate the copy number of the SHE gene, and then used general linear models to analyze the associations between CNV and economic traits. The results showed that there were more individuals with SHE copy number loss in CKS and HS than in STHS and LTHS individuals. Association analyses showed that gain and normal copy number types were correlated to body length, circumference of cannon bone, heart girth, chest width and high at the cross in CKS, HS and STHS (p < 0.05), but this association was not observed for LTHS. Chi-square values (χ²) found prominent differences in CNV distribution among the studied breeds. Overall, the CNV of the SHE gene may be an important consideration for sheep molecular breeding.

Novel mutations and their genotype-phenotype correlations in patients with Noonan syndrome, using next-generation sequencing

PURPOSE

Noonan Syndrome (NS) is an autosomal dominant disorder with many variable and heterogeneous conditions. The genetic basis for 20-30% of cases is still unknown. This study evaluates Iranian Noonan patients both clinically and genetically for the first time.

MATERIALS/METHODS

Mutational analysis of PTPN11 gene was performed in 15 Iranian patients, using PCR and Sanger sequencing at phase one. Then, as phase two, Next Generation Sequencing (NGS) in the form of targeted resequencing was utilized for analysis of exons from other related genes. Homology modelling for the novel founded mutations was performed as well. The genotype, phenotype correlation was done according to the molecular findings and clinical features.

RESULTS

Previously reported mutation (p.N308D) in some patients and a novel mutation (p.D155N) in one of the patients were identified in phase one. After applying NGS methods, known and new variants were found in four patients in other genes, including: CBL (p. V904I), KRAS (p. L53W), SOS1 (p. I1302V), and SOS1 (p. R552G). Structural studies of two deduced novel mutations in related genes revealed deficiencies in the mutated proteins. Following genotype, phenotype correlation, a new pattern of the presence of intellectual disability in two patients was registered.

CONCLUSIONS

NS shows strong variable expressivity along the high genetic heterogeneity especially in distinct populations and ethnic groups. Also possibly unknown other causative genes may be exist. Obviously, more comprehensive and new technologies like NGS methods are the best choice for detection of molecular defects in patients for genotype, phenotype correlation and disease management.

Global transformation of erythrocyte properties via engagement of an SH2-like sequence in band 3

Src homology 2 (SH2) domains are composed of weakly conserved sequences of ∼100 aa that bind phosphotyrosines in signaling proteins and thereby mediate intra- and intermolecular protein-protein interactions. In exploring the mechanism whereby tyrosine phosphorylation of the erythrocyte anion transporter, band 3, triggers membrane destabilization, vesiculation, and fragmentation, we discovered a SH2 signature motif positioned between membrane-spanning helices 4 and 5. Evidence that this exposed cytoplasmic sequence contributes to a functional SH2-like domain is provided by observations that: (i) it contains the most conserved sequence of SH2 domains, GSFLVR; (ii) it binds the tyrosine phosphorylated cytoplasmic domain of band 3 (cdb3-PO₄) with K_d = 14 nM; (iii) binding of cdb3-PO₄ to erythrocyte membranes is inhibited both by antibodies against the SH2 signature sequence and dephosphorylation of cdb3-PO₄; (iv) label transfer experiments demonstrate the covalent transfer of photoactivatable biotin from isolated cdb3-PO₄ (but not cdb3) to band 3 in erythrocyte membranes; and (v) phosphorylation-induced binding of cdb3-PO₄ to the membrane-spanning domain of band 3 in intact cells causes global changes in membrane properties, including (i) displacement of a glycolytic enzyme complex from the membrane, (ii) inhibition of anion transport, and (iii) rupture of the band 3-ankyrin bridge connecting the spectrin-based cytoskeleton to the membrane. Because SH2-like motifs are not retrieved by normal homology searches for SH2 domains, but can be found in many tyrosine kinase-regulated transport proteins using modified search programs, we suggest that related cases of membrane transport proteins containing similar motifs are widespread in nature where they participate in regulation of cell properties.

Membrane Binding and Modulation of the PDZ Domain of PICK1

Scaffolding proteins serve to assemble protein complexes in dynamic processes by means of specific protein-protein and protein-lipid binding domains. Many of these domains bind either proteins or lipids exclusively; however, it has become increasingly evident that certain domains are capable of binding both. Especially, many PDZ domains, which are highly abundant protein-protein binding domains, bind lipids and membranes. Here we provide an overview of recent large-scale studies trying to generalize and rationalize the binding patterns as well as specificity of PDZ domains towards membrane lipids. Moreover, we review how these PDZ-membrane interactions are regulated in the case of the synaptic scaffolding protein PICK1 and how this might affect cellular localization and function.

Molecular characterization of an h-type thioredoxin gene from grape

Thioredoxins (Trxs), as small ubiquitous proteins, participate in dithiol-disulfide exchange reactions. In contrast to other organisms, plants have a complex family of Trxs, which contains seven different Trx types: f, h, m, o, x, y, and z. The h-type Trx consists of multiple forms that are involved in different processes. A full-length cDNA coding for a Trx h, designated VvTrx h2, was isolated and cloned from grape (Vitis vinifera L. cv. White Seedless) berry tissue by RT-PCR technique. Nucleotide sequence analysis revealed 561 nucleotides in length encoded for a protein of 114 amino acid residues. The deduced polypeptide sequence harbors a typical catalytic site, WCGPC and its calculated molecular mass and its predicted isoelectric point are 12.79 and 5.06 kDa, respectively. The threedimensional modeling and docking studies allow for the proposal that VvTrx h2 could be reduced by a NADP-thioredoxin reductase rather than glutaredoxin, as shown for its ortholog from Arabidopsis. The deduced amino acid sequence showed a high degree of similarity to Trx h isoforms from other sources. Phylogenetic studies indicated that VvTrx h2 gene is related to h-type Trx subgroup I. Semi-quantitative RT-PCR analysis revealed that the VvTrx h2 gene was expressed in all plant tissues at different developmental stages.

Isolation, identification and sequence analysis of a thioredoxin h gene, a member of subgroup III of h-type Trxs from grape (Vitis vinifera L. cv. Askari)

Thioredoxins (Trxs) are small ubiquitous proteins which play a regulatory role in a variety of cellular processes. In contrast to other organisms, plants have a great number of Trx types, consisting of six well-defined groups: f, m, x, and y in chloroplasts, o in mitochondria, and h mainly in cytosol. A full-length cDNA, designated VvCxxS2, encoding Trx h polypeptide was isolated and cloned from grape (Vitis vinifera L. cv. Askari) berries organ by reverse transcription polymerase chain reaction (RT-PCR). The cDNA was 381 bp nucleotides in length with a deduced amino acid of 126 residues, possessing a WCIPS active site, which belongs to the subgroup III of h-type Trxs based on phylogenetic analysis. The calculated molecular mass and the predicted isoelectric point of the deduced polypeptide are 14.25 kDa and 4.68, respectively. Nucleotide sequence analysis of genomic DNA fragment of VvCxxS2 gene revealed that this gene possesses two introns at positions identical to the previously sequenced Trx h genes. A modeling analysis indicated that VvCxxS2 shares a common structure with other Trxs, and is preferably reduced by Grx rather than NADPH-dependent thioredoxin reductase (NTR). The deduced protein sequence showed a high similarity to Trx h from other plants, in particular from castor bean (Ricinus communis), Betula pendula and sweet orange (Citrus sinensis). Semiquantitative RT-PCR experiments indicated that the transcripts of VvCxxS2 gene are present in all plant organs and different developmental stages. In addition, the higher expression of the VvCxxS2 gene was observed in berry organ as compared to the other organs.

CAERUS: predicting CAncER oUtcomeS using relationship between protein structural information, protein networks, gene expression data, and mutation data

Carcinogenesis is a complex process with multiple genetic and environmental factors contributing to the development of one or more tumors. Understanding the underlying mechanism of this process and identifying related markers to assess the outcome of this process would lead to more directed treatment and thus significantly reduce the mortality rate of cancers. Recently, molecular diagnostics and prognostics based on the identification of patterns within gene expression profiles in the context of protein interaction networks were reported. However, the predictive performances of these approaches were limited. In this study we propose a novel integrated approach, named CAERUS, for the identification of gene signatures to predict cancer outcomes based on the domain interaction network in human proteome. We first developed a model to score each protein by quantifying the domain connections to its interacting partners and the somatic mutations present in the domain. We then defined proteins as gene signatures if their scores were above a preset threshold. Next, for each gene signature, we quantified the correlation of the expression levels between this gene signature and its neighboring proteins. The results of the quantification in each patient were then used to predict cancer outcome by a modified naïve Bayes classifier. In this study we achieved a favorable accuracy of 88.3%, sensitivity of 87.2%, and specificity of 88.9% on a set of well-documented gene expression profiles of 253 consecutive breast cancer patients with different outcomes. We also compiled a list of cancer-associated gene signatures and domains, which provided testable hypotheses for further experimental investigation. Our approach proved successful on different independent breast cancer data sets as well as an ovarian cancer data set. This study constitutes the first predictive method to classify cancer outcomes based on the relationship between the domain organization and protein network.

It is widely known that cancer is a complex process in which a large number of genes appear to be involved. Through experimental approaches, some oncogenes and tumor suppressors have been identified as playing important roles in the signaling and the regulatory pathways. However, we have not fully understood the complete mechanism of how cancer develops and how it leads to different disease outcomes (aggressive/dangerous or non-aggressive/less-dangerous). In order to identify a list of gene signatures and better predict cancer outcome, we developed an integrated and systematical approach by investigating gene expression profiling alternation caused by disruptions between protein-protein interactions and domain-domain interactions in the human interactome. Our approach achieves the favorable predictive performance if tested on a set of well-documented breast cancer patients, which suggests that the disrupted interactome is important to determine patient prognosis. Our approach is robust if tested on other independent data sets. This work provides a promising prognostic tool to classify different cancer outcomes.

Ligand-binding PAS domains in a genomic, cellular, and structural context

Per-Arnt-Sim (PAS) domains occur in proteins from all kingdoms of life. In the bacterial kingdom, PAS domains are commonly positioned at the amino terminus of signaling proteins such as sensor histidine kinases, cyclic-di-GMP synthases/hydrolases, and methyl-accepting chemotaxis proteins. Although these domains are highly divergent at the primary sequence level, the structures of dozens of PAS domains across a broad section of sequence space have been solved, revealing a conserved three-dimensional architecture. An all-versus-all alignment of 63 PAS structures demonstrates that the PAS domain family forms structural clades on the basis of two principal variables: (a) topological location inside or outside the plasma membrane and (b) the class of small molecule that they bind. The binding of a chemically diverse range of small-molecule metabolites is a hallmark of the PAS domain family. PAS ligand binding either functions as a primary cue to initiate a cellular signaling response or provides the domain with the capacity to respond to secondary physical or chemical signals such as gas molecules, redox potential, or photons. This review synthesizes the current state of knowledge of the structural foundations and evolution of ligand recognition and binding by PAS domains.

CLE14/CLE20 peptides may interact with CLAVATA2/CORYNE receptor-like kinases to irreversibly inhibit cell division in the root meristem of Arabidopsis

Towards an understanding of the interacting nature of the CLAVATA (CLV) complex, we predicted the 3D structures of CLV3/ESR-related (CLE) peptides and the ectodomain of their potential receptor proteins/kinases, and docking models of these molecules. The results show that the ectodomain of CLV1 can form homodimers and that the 12-/13-amino-acid CLV3 peptide fits into the binding clefts of the CLV1 dimers. Our results also demonstrate that the receptor domain of CORYNE (CRN), a recently identified receptor-like kinase, binds tightly to the ectodomain of CLV2, and this likely leads to an increased possibility for docking with CLV1. Furthermore, our docking models reveal that two CRN-CLV2 ectodomain heterodimers are able to form a tetramer receptor complex. Peptides of CLV3, CLE14, CLE19, and CLE20 are also able to bind a potential CLV2-CRN heterodimer or heterotetramer complex. Using a cell-division reporter line, we found that synthetic 12-amino-acid CLE14 and CLE20 peptides inhibit, irreversibly, root growth by reducing cell division rates in the root apical meristem, resulting in a short-root phenotype. Intriguingly, we observed that exogenous application of cytokinin can partially rescue the short-root phenotype induced by over-expression of either CLE14 or CLE20 in planta. However, cytokinin treatment does not rescue the short-root phenotype caused by exogenous application of the synthetic CLE14/CLE20 peptides, suggesting a requirement for a condition provided only in living plants. These results therefore imply that the CLE14/CLE20 peptides may act through the CLV2-CRN receptor kinase, and that their availabilities and/or abundances may be affected by cytokinin activity in planta.

Protein microarray assay for the screening of SH3 domain interactions

Analysis of cellular signal transduction processes increasingly focuses on the systematic characterization of complete protein interaction networks. Understanding the interplay of signaling components enables insight into the molecular basis of diverse diseases such as cancer. This paves the way for the rational design of specific therapeutics. Protein interactions are often mediated by conserved modular domains, e.g., SH3-domains, which recognize proline-rich sequences in their cognate ligands. In the course of this study, different microarray formats (reactive silane monolayers and nitrocellulose on glass slides) and assay work flows were evaluated to develop a microarray based screening assay that permits the reliable identification of interactions between certain target proteins with a set of SH3 domains. Nine representative SH3 domains which were produced and purified as GST-fusion proteins were spotted on the microarray substrates and probed with two well-characterized ligands, the Nef protein from HIV-1 and the human protein Sam68. The best results from these low-density model arrays were obtained with nitrocellulose slides. We show that a straightforward and highly robust detection of ligand binding is achieved by staining with a fluorescently labeled antibody directed against the N-terminal His-tag attached to these proteins. The optimized assay protocol reported here allows for the identification of SH3-interactions with high reproducibility and adequate signal-to-background and signal-to-noise ratios, as well as the quantitative determination of relative binding affinities.

A membrane-associated thioredoxin required for plant growth moves from cell to cell, suggestive of a role in intercellular communication

Thioredoxins (Trxs) are small ubiquitous regulatory disulfide proteins. Plants have an unusually complex complement of Trxs composed of six well-defined types (Trxs f, m, x, y, h, and o) that reside in different cell compartments and function in an array of processes. The extraplastidic h type consists of multiple members that in general have resisted isolation of a specific phenotype. In analyzing mutant lines in Arabidopsis thaliana, we identified a phenotype of dwarf plants with short roots and small yellowish leaves for AtTrx h9 (henceforth, Trx h9), a member of the Arabidopsis Trx h family. Trx h9 was found to be associated with the plasma membrane and to move from cell to cell. Controls conducted in conjunction with the localization of Trx h9 uncovered another h-type Trx in mitochondria (Trx h2) and a Trx in plastids earlier described as a cytosolic form in tomato. Analysis of Trx h9 revealed a 17-amino acid N-terminal extension in which the second Gly (Gly(2)) and fourth cysteine (Cys(4)) were highly conserved. Mutagenesis experiments demonstrated that Gly(2) was required for membrane binding, possibly via myristoylation. Both Gly(2) and Cys(4) were needed for movement, the latter seemingly for protein structure and palmitoylation. A three-dimensional model was consistent with these predictions as well as with earlier evidence showing that a poplar ortholog is reduced by a glutaredoxin rather than NADP-thioredoxin reductase. In demonstrating the membrane location and intercellular mobility of Trx h9, the present results extend the known boundaries of Trx and suggest a role in cell-to-cell communication.

PepCyber:P~PEP: a database of human protein protein interactions mediated by phosphoprotein-binding domains

Phosphoprotein-binding domains (PPBDs) mediate many important cellular and molecular processes. Ten PPBDs have been known to exist in the human proteome, namely, 14-3-3, BRCT, C2, FHA, MH2, PBD, PTB, SH2, WD-40 and WW. PepCyber:P∼PEP is a newly constructed database specialized in documenting human PPBD-containing proteins and PPBD-mediated interactions. Our motivation is to provide the research community with a rich information source emphasizing the reported, experimentally validated data for specific PPBD–PPEP interactions. This information is not only useful for designing, comparing and validating the relevant experiments, but it also serves as a knowledge-base for computationally constructing systems signaling pathways and networks. PepCyber:P∼PEP is accessible through the URL, http://www.pepcyber.org/PPEP/. The current release of the database contains 7044 PPBD-mediated interactions involving 337 PPBD-containing proteins and 1123 substrate proteins.

Early-stage folding in proteins (in silico) sequence-to-structure relation

A sequence-to-structure library has been created based on the complete PDB database. The tetrapeptide was selected as a unit representing a well-defined structural motif. Seven structural forms were introduced for structure classification. The early-stage folding conformations were used as the objects for structure analysis and classification. The degree of determinability was estimated for the sequence-to-structure and structure-to-sequence relations. Probability calculus and informational entropy were applied for quantitative estimation of the mutual relation between them. The structural motifs representing different forms of loops and bends were found to favor particular sequences in structure-to-sequence analysis.

Evaluation and improvement of multiple sequence methods for protein secondary structure prediction

A new dataset of 396 protein domains is developed and used to evaluate the performance of the protein secondary structure prediction algorithms DSC, PHD, NNSSP, and PREDATOR. The maximum theoretical Q3 accuracy for combination of these methods is shown to be 78%. A simple consensus prediction on the 396 domains, with automatically generated multiple sequence alignments gives an average Q3 prediction accuracy of 72.9%. This is a 1% improvement over PHD, which was the best single method evaluated. Segment Overlap Accuracy (SOV) is 75.4% for the consensus method on the 396-protein set. The secondary structure definition method DSSP defines 8 states, but these are reduced by most authors to 3 for prediction. Application of the different published 8- to 3-state reduction methods shows variation of over 3% on apparent prediction accuracy. This suggests that care should be taken to compare methods by the same reduction method. Two new sequence datasets (CB513 and CB251) are derived which are suitable for cross-validation of secondary structure prediction methods without artifacts due to internal homology. A fully automatic World Wide Web service that predicts protein secondary structure by a combination of methods is available via http://barton.ebi.ac.uk/.

Structure and function of phosphoinositide 3-kinases

Phosphoinositide kinases (PI3Ks) play an important role in mitogenic signaling and cell survival, cytoskeletal remodeling, metabolic control and vesicular trafficking. Here we summarize the structure-function relationships delineating the activation process of class I PI3Ks involving various domains of adapter subunits, Ras, and interacting proteins. The resulting product, PtdIns(3,4,5)P3, targets Akt/protein kinase B (PKB), Bruton's tyrosine kinase (Btk), phosphoinositide-dependent kinases (PDK), integrin-linked kinase (ILK), atypical protein kinases C (PKC), phospholipase Cgamma and more. Surface receptor-activated PI3Ks function in mammals, insects, nematodes and slime mold, but not yeast. While many members of the class II family have been identified and characterized biochemically, it is presently unknown how these C2-domain containing PI3Ks are activated, and which PI substrate they phosphorylate in vivo. PtdIns 3-P is produced by Vps34p/class III PI3Ks and operates via the PtdIns 3-P-binding proteins early endosomal antigen (EEA1), yeast Vac1p, Vps27p, Pip1p in lysosomal protein targeting. Besides the production of D3 phosphorylated lipids, PI3Ks have an intrinsic protein kinase activity. For trimeric GTP-binding protein-activated PI3Kgamma, protein kinase activity seems to be sufficient to trigger mitogen-activated protein kinase (MAPK). Recent disruption of PI3K genes in slime mold, Caenorhabditis elegans, Drosophila melanogaster and mice further underlines the importance of PI3K signaling systems and elucidates the role of PI3K signaling in multicellular organisms.

Identification of phosphopeptide ligands for the Src-homology 2 (SH2) domain of Grb2 by phage display

We report here on the identification of phosphopetide ligands which interact with the Src-homology 2 (SH2) domain of the adapter protein Grb2 by screening a random peptide library established on phage. Phage were phosphorylated in vitro at an invariant tyrosine residue by a mixture of phosphotyrosine kinases c-Src, Blk and Syk. Selection of binding motifs was carried out by interaction of the library with the recombinant SH2 domain of Grb2 expressed as a glutathione S-transferase (GST) fusion protein. Several subsequent cycles of selection led to the enrichment of phage which bound to the GST-Grb2 SH2 domain only when previously phosphorylated. Sequence analysis revealed that all of the selected phage displayed peptides with the consensus motif Y*M/ENW (Y* denotes phosphotyrosine). One of these peptides, bearing the Y*ENW motif, bound the Grb2 SH2 domain with a threefold higher affinity than the peptide motif Y*VNV derived from the natural ligand Shc. Thus, phage display can be employed to rapidly identify high affinity ligands to SH2 domains.

Predicted structure of the extracellular region of ligand-gated ion-channel receptors shows SH2-like and SH3-like domains forming the ligand-binding site

Fast synaptic neurotransmission is mediated by ligand-gated ion-channel (LGIC) receptors, which include receptors for acetylcholine, serotonin, GABA, glycine, and glutamate. LGICs are pentamers with extracellular ligand-binding domains and form integral membrane ion channels that are selective for cations (acetylcholine and serotonin 5HT3 receptors) or anions (GABAA and glycine receptors and the invertebrate glutamate-binding chloride channel). They form a protein superfamily with no sequence similarity to any protein of known structure. Using a 1D-3D structure mapping approach, we have modeled the extracellular ligand-binding domain based on a significant match with the SH2 and SH3 domains of the biotin repressor structure. Refinement of the model based on knowledge of the large family of SH2 and SH3 structures, sequence alignments, and use of structure templates for loop building, allows the prediction of both monomer and pentamer models. These are consistent with medium-resolution electron microscopy structures and with experimental structure/function data from ligand-binding, antibody-binding, mutagenesis, protein-labeling and subunit-linking studies, and glycosylation sites. Also, the predicted polarity of the channel pore calculated from electrostatic potential maps of pentamer models of superfamily members is consistent with known ion selectivities. Using the glycine receptor alpha 1 subunit, which forms homopentamers, the monomeric and pentameric models define the agonist and antagonist (strychnine) binding sites to a deep crevice formed by an extended loop, which includes the invariant disulfide bridge, between the SH2 and SH3 domains. A detailed binding site for strychnine is reported that is in strong agreement with known structure/function data. A site for interaction of the extracellular ligand-binding domain with the activation of the M2 transmembrane helix is also suggested.

Protein secondary structure prediction using local alignments

The accuracy of secondary structure prediction methods has been improved significantly by the use of aligned protein sequences. The PHD method and the NNSSP method reach 71 to 72% of sustained overall three-state accuracy when multiple sequence alignments are with neural networks and nearest-neighbor algorithms, respectively. We introduce a variant of the nearest-neighbor approach that can achieve similar accuracy using a single sequence as the query input. We compute the 50 best non-intersecting local alignments of the query sequence with each sequence from a set of proteins with known 3D structures. Each position of the query sequence is aligned with the database amino acids in alpha-helical, beta-strand or coil states. The prediction type of secondary structure is selected as the type of aligned position with the maximal total score. On the dataset of 124 non-membrane non-homologous proteins, used earlier as a benchmark for secondary structure predictions, our method reaches an overall three-state accuracy of 71.2%. The performance accuracy is verified by an additional test on 461 non-homologous proteins giving an accuracy of 71.0%. The main strength of the method is the high level of prediction accuracy for proteins without any known homolog. Using multiple sequence alignments as input the method has a prediction accuracy of 73.5%. Prediction of secondary structure by the SSPAL method is available via Baylor College of Medicine World Wide Web server.

Interaction of phosphorylated FcepsilonRIgamma immunoglobulin receptor tyrosine activation motif-based peptides with dual and single SH2 domains of p72syk. Assessment of binding parameters and real time binding kinetics

To examine the characteristics of the interaction of the FcepsilonRIgamma ITAM with the SH2 domains of p72(syk), the binding of an 125I-labeled dual phosphorylated FcepsilonRIgamma ITAM-based peptide to the p72(syk) SH2 domains was monitored utilizing a novel scintillation proximity based assay. The Kd for this interaction, determined from the saturation binding isotherm, was 1.4 nM. This high affinity binding was reflected in the rapid rate of association for the peptide binding to the SH2 domains. Competition studies utilizing a soluble C-terminal SH2 domain knockout and N-terminal SH2 domain knockouts revealed that both domains contribute cooperatively to the high affinity binding. Unlabeled dual phosphorylated peptide competed with the 125I-labeled peptide for binding to the dual p72(syk) SH2 domains with an IC50 value of 4.8 nM. Monophosphorylated 24-mer FcepsilonRIgamma ITAM peptides, and phosphotyrosine also competed for binding, but with substantially higher IC50 values. This, and other data discussed, suggest that high affinity binding requires both tyrosine residues to be phosphorylated and that the preferred binding orientation of the ITAM is such that the N-terminal phosphotyrosine occupies the C-terminal SH2 domain and the C-terminal phosphotyrosine occupies the N-terminal SH2 domain.

Modulation of the SH2 binding specificity and kinase activity of Src by tyrosine phosphorylation within its SH2 domain

The Src family of kinases are held in an inactive state by interaction of their SH2 domain with a C-terminal phosphotyrosine. Dephosphorylation of this site can reactivate Src; however, recent evidence suggests that activation can also occur without dephosphorylation. In this study, platelet-derived growth factor receptor phosphorylation of Src on Tyr-213 specifically blocked binding of its SH2 domain to a phosphopeptide corresponding to the C-terminal regulatory sequence, while binding to other sequences, such as the platelet-derived growth factor receptor or a peptide from the epidermal growth factor receptor, was unaffected. Consequently, Src was activated over 50-fold. This is the first demonstration of regulation of a SH2 domain specificity by post-translational modification and is likely to be a general mechanism for regulation of all Src-like kinases.

Helical fold prediction for the cyclin box

The smooth progression of the eukaryotic cell cycle relies on the periodic activation of members of a family of cell cycle kinases by regulatory proteins called cyclins. Outside of the cell cycle, cyclin homologs play important roles in regulating the assembly of transcription complexes; distant structural relatives of the conserved cyclin core or "box" can also function as general transcription factors (like TFIIB) or survive embedded in the chain of the tumor suppressor, retinoblastoma protein. The present work attempts the prediction of the canonical secondary, supersecondary, and tertiary fold of the minimal cyclin box domain using a combination of techniques that make use of the evolutionary information captured in a multiple alignment of homolog sequences. A tandem set of closely packed, helical modules are predicted to form the cyclin box domain.

A simple and fast approach to prediction of protein secondary structure from multiply aligned sequences with accuracy above 70%

To improve secondary structure predictions in protein sequences, the information residing in multiple sequence alignments of substituted but structurally related proteins is exploited. A database comprised of 70 protein families and a total of 2,500 sequences, some of which were aligned by tertiary structural superpositions, was used to calculate residue exchange weight matrices within alpha-helical, beta-strand, and coil substructures, respectively. Secondary structure predictions were made based on the observed residue substitutions in local regions of the multiple alignments and the largest possible associated exchange weights in each of the three matrix types. Comparison of the observed and predicted secondary structure on a per-residue basis yielded a mean accuracy of 72.2%. Individual alpha-helix, beta-strand, and coil states were respectively predicted at 66.7, and 75.8% correctness, representing a well-balanced three-state prediction. The accuracy level, verified by cross-validation through jack-knife tests on all protein families, dropped, on average, to only 70.9%, indicating the rigor of the prediction procedure. On the basis of robustness, conceptual clarity, accuracy, and executable efficiency, the method has considerable advantage, especially with its sole reliance on amino acid substitutions within structurally related proteins.

Evaluation of current techniques for ab initio protein structure prediction

The results of a protein structure prediction contest are reviewed. Twelve different groups entered predictions on 14 proteins of known sequence whose structures had been determined but not yet disseminated to the scientific community. Thus, these represent true tests of the current state of structure prediction methodologies. From this work, it is clear that accurate tertiary structure prediction is not yet possible. However, protein fold and motif prediction are possible when the motif is recognizable similar to another known structure. Internal symmetry and the information inherent in an aligned family of homologous sequences facilitate predictive efforts. Novel folds remain a major challenge for prediction efforts.

A predicted consensus structure for the protein kinase C2 homology (C2H) domain, the repeating unit of synaptotagmin

A secondary structure has been predicted for the protein kinase C2 regulatory domain found in homologous form in synaptotagmin, some phospholipases, and some GTP activated proteins. The proposed structure is built from seven consecutive beta strands followed by a terminal alpha helix. Considerations of overall surface exposure of individual secondary structural elements suggest that these are packed into a 2-sheet beta sandwich structure, with one of only three of the many possible folds being preferred.

SH2 domain structure and function

An emerging theme in both the biology of signal transduction and the biochemistry of proteins has been the modular function of small protein domains. In some cases these can directly regulate catalytic activity. In others, they serve to interconnect important regulatory proteins. SH2 (src homology 2) domains represent some of the best studied models. Originally identified on the basis of homology in src and fps [1], SH2s are elements that ordinarily respond to tyrosine phosphorylation by binding the phosphorylated sequence. As such, they are key elements in tyrosine kinase regulation of cellular processes. Because SH2 interactions result from phosphorylation, such elements provide a regulatable circuitry along which signals can be transmitted in a timely manner. Because the regulation is based on a common mechanism, signal generators can target several different proteins coordinately. The PDGF receptor (PDGFr), for example, may interact with as many as ten different elements [2,3]. There are a number of excellent reviews on SH2 domains available [4-11]. This discussion will try to show how genetic, biochemical and biophysical results can be integrated in a satisfying way.

A consensus prediction of the secondary structure for the 6-phospho-beta-D-galactosidase superfamily

Two separate unrefined models for the secondary structure of two subfamilies of the 6-phospho-beta-D-galactosidase superfamily were independently constructed by examining patterns of variation and conservation within homologous protein sequences, assigning surface, interior, parsing, and active site residues to positions in the alignment, and identifying periodicities in these. A consensus model for the secondary structure of the entire superfamily was then built. The prediction tests the limits of an unrefined prediction made using this approach in a large protein with substantial functional and sequence divergence within the family. The protein belongs to the (alpha-beta class), with the core beta strands aligned parallel. The supersecondary structural elements that are readily identified in this model is a parallel beta sheet built by strands C, D, and E, with helices 2 and 3 connecting strands (C+D) and (D+E), respectively, and an analogous beta-alpha unit (strand G and helix 7) toward the end of the sequence. The resemblance of the supersecondary model to the tertiary structure formed by 8-fold alpha-beta barrel proteins is almost certainly not coincidental.

Protein-protein interactions: methods for detection and analysis

The function and activity of a protein are often modulated by other proteins with which it interacts. This review is intended as a practical guide to the analysis of such protein-protein interactions. We discuss biochemical methods such as protein affinity chromatography, affinity blotting, coimmunoprecipitation, and cross-linking; molecular biological methods such as protein probing, the two-hybrid system, and phage display: and genetic methods such as the isolation of extragenic suppressors, synthetic mutants, and unlinked noncomplementing mutants. We next describe how binding affinities can be evaluated by techniques including protein affinity chromatography, sedimentation, gel filtration, fluorescence methods, solid-phase sampling of equilibrium solutions, and surface plasmon resonance. Finally, three examples of well-characterized domains involved in multiple protein-protein interactions are examined. The emphasis of the discussion is on variations in the approaches, concerns in evaluating the results, and advantages and disadvantages of the techniques.

Predicted secondary and supersecondary structure for the serine-threonine-specific protein phosphatase family

A bona fide consensus prediction for the secondary and supersecondary structure of the serine-threonine specific protein phosphatases is presented. The prediction includes assignments of active site segments, an internal helix, and a region of possible 3(10) helical structure. An experimental structure for a member of this family of proteins should appear shortly, allowing this prediction to be evaluated.

Predicting the conformation of proteins from sequences. Progress and future progress

Recent progress in structure prediction has allowed bona fide predictions, those made and published before an experimental structure is determined, to be remarkably accurate. The most successful methods rely on an analysis of patterns of conservation and variation within homologous protein sequences, extract tertiary structural information before secondary structure is predicted, and avoid 'three state per residue scores' as a tool for evaluating a prediction, focusing instead on efforts to understand why a prediction is successful when it is successful, and why it fails when it fails.

A secondary structure model of the integrin alpha subunit N-terminal domain based on analysis of multiple alignments

The integrins are alpha/beta heterodimeric proteins which mediate cell-matrix and cell-cell interactions. Current data indicate that the N-terminal moiety of the alpha subunit is involved in ligand binding. This region of the receptor is made up of a seven-fold repeated sequence of unknown structure which contains EF-hand-like putative divalent cation-binding sites. Recent studies have shown that multiple sequence alignments can be analysed to yield secondary structure predictions. Therefore, to obtain a model structure for the integrin alpha subunit N-terminal domain repeat, a large alignment of the seven repeats from sixteen integrin sequences was generated. Two methods of analysis were used: First, Chou and Fasman and Garnier, Osguthorpe and Robson predictions were carried out for individual sequences and the consensus predictions derived. Consensus hydrophobicity and chain flexibility data were also used to provide additional data. Second, sites of conservation and variation were analysed by a computer program STAMA (STructure After Multiple Alignment) to yield a secondary structure prediction. The two analyses gave essentially the same predicted structure: undefined region, loop, alpha-helix, beta-strand, divalent cation-binding loop, beta-strand, putative turn, loop, beta-strand. This is the first model structure to be presented for an integrin domain. Its implications for integrin function are discussed.