An in silico proteomics screen to predict and prioritize protein-protein interactions dependent on post-translationally modified motifs.

Motivation: The development of proteomic methods for the characterization of domain/motif interactions has greatly expanded our understanding of signal transduction. However, proteomics-based binding screens have limitations including that the queried tissue or cell type may not harbor all potential interacting partners or post-translational modifications (PTMs) required for the interaction. Therefore, we sought a generalizable, complementary in silico approach to identify potentially novel motif and PTM-dependent binding partners of high priority. Results: We used as an initial example the interaction between the Src homology 2 (SH2) domains of the adaptor proteins CT10 regulator of kinase (CRK) and CRK-like (CRKL) and phosphorylated-YXXP motifs. Employing well-curated, publicly-available resources, we scored and prioritized potential CRK/CRKL-SH2 interactors possessing signature characteristics of known interacting partners. Our approach gave high priority scores to 102 of the >9000 YXXP motif-containing proteins. Within this 102 were 21 of the 25 curated CRK/CRKL-SH2-binding partners showing a more than 80-fold enrichment. Several predicted interactors were validated biochemically. To demonstrate generalized applicability, we used our workflow to predict protein-protein interactions dependent upon motif-specific arginine methylation. Our data demonstrate the applicability of our approach to, conceivably, any modular binding domain that recognizes a specific post-translationally modified motif. Supplementary information: Supplementary data are available at Bioinformatics online.

Comparative phosphoproteomic analysis of neonatal and adult murine brain.

Developmental processes are governed by diverse regulatory mechanisms including a suite of signaling pathways employing reversible phosphorylation. With the advent of large-scale phosphoproteomics, it is now possible to identify thousands of phosphorylation sites from tissues at distinct developmental stages. We describe here the identification of over 6000 nonredundant phosphorylation sites from neonatal murine brain. When compared to nearly three times the number of phosphorylation sites identified from 3-week-old murine brain, remarkably one-third of the neonatal sites were unique. This fraction only dropped to one-quarter when allowing the site to stray plus or minus 15 residues. This provides evidence for considerable change in the profiles of developmentally regulated phosphoproteomes. Using quantitative MS we characterized a novel phosphorylation site (Ser265) identified uniquely in the neonatal brain on doublecortin (Dcx), a protein essential for proper mammalian brain development. While the relative levels of Dcx and phospho-Ser265 Dcx between embryonic and neonatal brain were similar, their levels fell precipitously by postnatal day 21, as did phospho-Ser297, a site required for proper neuronal migration. Both sites lie near the microtubule-binding domain and may provide functionally similar regulation via different kinases.

SE: an algorithm for deriving sequence alignment from a pair of superimposed structures.

BACKGROUND: Generating sequence alignments from superimposed structures is an important part of many structure comparison programs. The accuracy of the alignment affects structure recognition, classification and possibly function prediction. Many programs use a dynamic programming algorithm to generate the sequence alignment from superimposed structures. However, this procedure requires using a gap penalty and, depending on the value of the penalty used, can introduce spurious gaps and misalignments. Here we present a new algorithm, Seed Extension (SE), for generating the sequence alignment from a pair of superimposed structures. The SE algorithm first finds "seeds", which are the pairs of residues, one from each structure, that meet certain stringent criteria for being structurally equivalent. Three consecutive seeds form a seed segment, which is extended along the diagonal of the alignment matrix in both directions. Distance and the amino acid type similarity between the residues are used to resolve conflicts that arise during extension of more than one diagonal. The manually curated alignments in the Conserved Domain Database were used as the standard to assess the quality of the sequence alignments. RESULTS: SE gave an average accuracy of 95.9% over 582 pairs of superimposed proteins tested, while CHIMERA, LSQMAN, and DP extracted from SHEBA, which all use a dynamic programming algorithm, yielded 89.9%, 90.2% and 91.0%, respectively. For pairs of proteins with low sequence or structural similarity, SE produced alignments up to 18% more accurate on average than the next best scoring program. Improvement was most pronounced when the two superimposed structures contained equivalent helices or beta-strands that crossed at an angle. When the SE algorithm was implemented in SHEBA to replace the dynamic programming routine, the alignment accuracy improved by 10% on average for structure pairs with RMSD between 2 and 4 A. SE also used considerably less CPU time than DP. CONCLUSION: The Seed Extension algorithm is fast and, without using a gap penalty, produces more accurate sequence alignments from superimposed structures than three other programs tested that use dynamic programming algorithm.

Chagas disease vector blood meal sources identified by protein mass spectrometry.

Chagas disease is a complex vector borne parasitic disease involving blood feeding Triatominae (Hemiptera: Reduviidae) insects, also known as kissing bugs, and the vertebrates they feed on. This disease has tremendous impacts on millions of people and is a global health problem. The etiological agent of Chagas disease, Trypanosoma cruzi (Kinetoplastea: Trypanosomatida: Trypanosomatidae), is deposited on the mammalian host in the insect's feces during a blood meal, and enters the host's blood stream through mucous membranes or a break in the skin. Identifying the blood meal sources of triatomine vectors is critical in understanding Chagas disease transmission dynamics, can lead to identification of other vertebrates important in the transmission cycle, and aids management decisions. The latter is particularly important as there is little in the way of effective therapeutics for Chagas disease. Several techniques, mostly DNA-based, are available for blood meal identification. However, further methods are needed, particularly when sample conditions lead to low-quality DNA or to assess the risk of human cross-contamination. We demonstrate a proteomics-based approach, using liquid chromatography tandem mass spectrometry (LC-MS/MS) to identify host-specific hemoglobin peptides for blood meal identification in mouse blood control samples and apply LC-MS/MS for the first time to Triatoma dimidiata insect vectors, tracing blood sources to species. In contrast to most proteins, hemoglobin, stabilized by iron, is incredibly stable even being preserved through geologic time. We compared blood stored with and without an anticoagulant and examined field-collected insect specimens stored in suboptimal conditions such as at room temperature for long periods of time. To our knowledge, this is the first study using LC-MS/MS on field-collected arthropod disease vectors to identify blood meal composition, and where blood meal identification was confirmed with more traditional DNA-based methods. We also demonstrate the potential of synthetic peptide standards to estimate relative amounts of hemoglobin acquired when insects feed on multiple blood sources. These LC-MS/MS methods can contribute to developing Ecohealth control strategies for Chagas disease transmission and can be applied to other arthropod disease vectors.

Metagenomic investigation of the microbial diversity in a chrysotile asbestos mine pit pond, Lowell, Vermont, USA.

Here we report on a metagenomics investigation of the microbial diversity in a serpentine-hosted aquatic habitat created by chrysotile asbestos mining activity at the Vermont Asbestos Group (VAG) Mine in northern Vermont, USA. The now-abandoned VAG Mine on Belvidere Mountain in the towns of Eden and Lowell includes three open-pit quarries, a flooded pit, mill buildings, roads, and > 26 million metric tons of eroding mine waste that contribute alkaline mine drainage to the surrounding watershed. Metagenomes and water chemistry originated from aquatic samples taken at three depths (0.5 m, 3.5 m, and 25 m) along the water column at three distinct, offshore sites within the mine's flooded pit (near 44°46'00.7673″, - 72°31'36.2699″; UTM NAD 83 Zone 18 T 0695720 E, 4960030 N). Whole metagenome shotgun Illumina paired-end sequences were quality trimmed and analyzed based on a translated nucleotide search of NCBI-NR protein database and lowest common ancestor taxonomic assignments. Our results show strata within the pit pond water column can be distinguished by taxonomic composition and distribution, pH, temperature, conductivity, light intensity, and concentrations of dissolved oxygen. At the phylum level, metagenomes from 0.5 m and 3.5 m contained a similar distribution of taxa and were dominated by Actinobacteria (46% and 53% of reads, respectively), Proteobacteria (45% and 38%, respectively), and Bacteroidetes (7% in both). The metagenomes from 25 m showed a greater diversity of phyla and a different distribution of reads than the two upper strata: Proteobacteria (60%), Actinobacteria (18%), Planctomycetes, (10%), Bacteroidetes (5%) and Cyanobacteria (2.5%), Armatimonadetes (< 1%), Verrucomicrobia (< 1%), Firmicutes (< 1%), and Nitrospirae (< 1%). Raw metagenome sequence data from each sample reside in NCBI's Short Read Archive (SRA ID: SRP056095) and are accessible through NCBI BioProject PRJNA277916.

Assessment of CASP6 predictions for new and nearly new fold targets.

This is a report of the assessment of the predictions made for the CASP6 protein structure prediction experiment conducted in 2004 in the New Fold (NF) category. There were nine protein domains that were judged to have new folds (NF) and 16 for which a similar structure was known but the sequence similarity was judged to be too low for them to be easily recognized (FR/A). We selected all NF targets and eight of the 16 FR/A targets judged to be at the borderline between NF and FR/A for evaluation in the NF category. A total of 165 prediction groups submitted over 7400 structural models for these targets. The quality of these models was evaluated using the GDT_TS scores of the structural similarity detection program LGA and by visual inspection of the top-scoring models. The best models submitted bore an overall similarity to the target structure for three or four of the nine NF targets and for all but one of the FR/A targets. High-scoring models for the NF targets were submitted by several different groups. When both the NF and FR/A targets were considered, Baker group dominated by submitting best models for seven of the 17 targets, but 14 other groups also managed to submit best models for one or more targets.

Evaluation of domain prediction in CASP6.

We present an analysis of the domain boundary prediction, a new category, in the sixth community-wide experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP6). There were 1011 predictions submitted for 63 targets. Each prediction was compared to the set of domains defined manually by visual inspection of the experimental structure. The comparison was scored using a new domain prediction scoring scheme. As the definition of a domain is subjective, many targets were assigned alternate definitions. For such targets, each prediction was compared with all different definitions and the best score was chosen. The predictors found it difficult to accurately predict domain boundaries when the target protein contained many domains or domains made of multiple sequence segments. The CBRC-DR (P0536) and Sternberg (P0237) groups were the most successful among human experts, while Baker-Rossettadom (P0353) and Baker-Robetta-Ginzu (P0421) did well among servers.

Quantitative comparison of CrkL-SH3 binding proteins from embryonic murine brain and liver: Implications for developmental signaling and the quantification of protein species variants in bottom-up proteomics.

A major aim of proteomics is to comprehensively identify and quantify all protein species variants from a given biological source. However, in spite of its tremendous utility, bottom-up proteomic strategies can do little to provide true quantification of distinct whole protein species variants given its reliance on proteolysis. This is particularly true when molecular size information is lost as in gel-free proteomics. Crk and CrkL comprise a family of adaptor proteins that couple upstream phosphotyrosine signals to downstream effectors by virtue of their SH2 and SH3 domains respectively. Here we compare the identification and quantification of CrkL-SH3 binding partners between embryonic murine brain and liver. We also uncover and quantify tissue-specific variants in CrkL-SH3 binding proteins.

Characterization of the B cell epitopes associated with a truncated form of Pseudomonas exotoxin (PE38) used to make immunotoxins for the treatment of cancer patients.

Recombinant immunotoxins composed of an Ab Fv fragment joined to a truncated portion of Pseudomonas exotoxin A (termed PE38) have been evaluated in clinical trials for the treatment of various human cancers. Immunotoxin therapy is very effective in hairy cell leukemia and also has activity in other hemological malignancies; however, a neutralizing Ab response to PE38 in patients with solid tumors prevents repeated treatments to maximize the benefit. In this study, we analyze the murine Ab response as a model to study the B cell epitopes associated with PE38. Sixty distinct mAbs to PE38 were characterized. Mutual competitive binding of the mAbs indicated the presence of 7 major epitope groups and 13 subgroups. The competition pattern indicated that the epitopes are discrete and could not be reproduced using a computer simulation program that created epitopes out of random surface residues on PE38. Using sera from immunotoxin-treated patients, the formation of human Abs to each of the topographical epitopes was demonstrated. One epitope subgroup, E1a, was identified as the principal neutralizing epitope. The location of each epitope on PE38 was determined by preparing 41 mutants of PE38 in which bulky surface residues were mutated to either alanine or glycine. All 7 major epitope groups and 9 of 13 epitope subgroups were identified by 14 different mutants and these retained high cytotoxic activity. Our results indicate that a relatively small number of discrete immunogenic sites are associated with PE38, most of which can be eliminated by point mutations.

High expression of a cytokeratin-associated protein in many cancers.

We have described previously a cDNA library made from membrane-bound polysomal mRNA prepared from breast and prostate cancer cell lines. The library is highly enriched for cDNAs encoding membrane proteins, secreted proteins, and cytokeratins. To characterize this library, 25,277 cDNA clones were sequenced and aligned with various databases; 1,439 clones did not align with known genes. From this set of clones we identified a previously uncharacterized gene encoding a 334-aa protein. Although protein structural motif prediction programs indicate that the gene encodes a membrane protein comprising a signal sequence, a series of leucine-rich repeats, and a single transmembrane domain with a cytoplasmic tail, confocal microscopy of MCF7 breast cancer cells demonstrates that the protein is not directly associated with the plasma membrane or intracellular membranes but instead colocalizes with intermediate filaments and cytokeratins within the cell. Immunofluorescence studies also show that protein expression is increased greatly in mitotic MCF7 cells, and immunohistochemistry demonstrates its expression in human breast cancer cells. Analysis of mRNA levels in 25 different normal tissues by RT-PCR shows that this gene is expressed highly in normal prostate and salivary gland, very weakly in colon, pancreas, and intestine, and not at all in other tissues. RT-PCR studies on human cancer samples show that the RNA is expressed highly in many cancer cell lines and cancer specimens, including 26 of 33 human breast cancers, 3 of 3 prostate cancers, 3 of 3 colon cancers, and 3 of 3 pancreatic cancers. We name the protein CAPC, cytokeratin-associated protein in cancer.

Mutants of immunotoxin anti-Tac(dsFv)-PE38 with variable number of lysine residues as candidates for site-specific chemical modification. 1. Properties of mutant molecules.

Chemical modification of proteins with substances such as poly(ethylene glycol) can add useful properties to proteins. Currently PEGylation is done in a random manner utilizing amino residues dispersed throughout a protein. For proteins such as immunotoxins, which have several different functional domains, random modification leads to inactivation. To determine if we could produce an immunotoxin with a diminished number of lysine residues so that chemical modification could be restricted to certain regions of the protein, we chose the recombinant immunotoxin anti-Tac(dsFv)-PE38 that has 13 lysine residues in the Fv portion and 3 in the toxin. We prepared a series of mutants with 0-12 lysines in the Fv and 0 or 3 in the toxin. Almost all of these molecules retain full biological activity. Our data indicate that replacement of lysine residues can be achieve without loss of biological potency. These molecules are a useful starting point to carry out site-specific PEGylation experiments.

Discovery of the breast cancer gene BASE using a molecular approach to enrich for genes encoding membrane and secreted proteins.

To identify unknown membrane proteins that could be used as targets for breast and prostate cancer immunotherapies and secreted proteins to be used as diagnostic markers, a cDNA library was generated from membrane-associated polyribosomal RNA derived from four breast cancer cell lines, one normal breast cell line, and a prostate cancer cell line. The membrane-associated polyribosomal cDNA library was subtracted with RNA from normal brain, liver, lung, kidney, and muscle. Of the 15,581 clones sequenced from the subtracted cDNA library, sequences from 10,506 clones map to known genes, but 5,075 sequences, representing 3,181 unique transcripts, are not associated with known genes. As one example, we experimentally investigated expression of a previously uncharacterized breast cancer gene that encodes a secreted protein designated BASE (breast cancer and salivary gland expression). BASE is expressed in many breast cancers but not in essential normal tissues including the five organs used for subtraction. Further analysis of this library should yield additional gene products of use in the diagnosis or treatment of breast or prostate cancer.

PATE, a gene expressed in prostate cancer, normal prostate, and testis, identified by a functional genomic approach.

To identify target antigens for prostate cancer therapy, we have combined computer-based screening of the human expressed sequence tag database and experimental expression analysis to identify genes that are expressed in normal prostate and prostate cancer but not in essential human tissues. Using this approach, we identified a gene that is expressed specifically in prostate cancer, normal prostate, and testis. The gene has a 1.5-kb transcript that encodes a protein of 14 kDa. We named this gene PATE (expressed in prostate and testis). In situ hybridization shows that PATE mRNA is expressed in the epithelial cells of prostate cancers and in normal prostate. Transfection of the PATE cDNA with a Myc epitope tag into NIH 3T3 cells and subsequent cell fractionation analysis shows that the PATE protein is localized in the membrane fraction of the cell. Analysis of the amino acid sequence of PATE shows that it has structural similarities to a group of proteins known as three-finger toxins, which includes the extracellular domain of the type beta transforming growth factor receptor. Restricted expression of PATE makes it a potential candidate for the immunotherapy of prostate cancer.