Predicting protein interaction interfaces from protein sequences: Case studies of subtilisin and phycocyanin

Identification of protein interaction interfaces is very important for understanding the molecular mechanisms underlying biological phenomena. Here, we present a novel method for predicting protein interaction interfaces from sequences by using PAM matrix (PIFPAM). Sequence alignments for interacting proteins were constructed and parsed into segments using sliding windows. By calculating distance matrix for each segment, the correlation coefficients between segments were estimated. The interaction interfaces were predicted by extracting highly correlated segment pairs from the correlation map. The predictions achieved an accuracy 0.41-0.71 for eight intraprotein interaction examples, and 0.07-0.60 for four interprotein interaction examples. Compared with three previously published methods, PIFPAM predicted more contacting site pairs for 11 out of the 12 example proteins, and predicted at least 34% more contacting site pairs for eight proteins of them. The factors affecting the predictions were also analyzed. Since PIFPAM uses only the alignments of the two interacting proteins as input, it is especially useful when no three-dimensional protein structure data are available.

The responsive expression of heat shock protein 22 gene in zhikong scallop Chlamys farreri against a bacterial challenge

HSP22 is a member of a small HSP subfamily contributing to the growth, transformation and apoptosis of the cell as well as acting as a molecular chaperone. In the present study, CfHSP22 cDNA was cloned from Chlamys farreri by the rapid amplification of cDNA ends technique. The full-length cDNA of CfHSP22 was of 1279 bp, consisting of a 5'-terminal untranslated region (5'UTR) of 122 bp, a 3'UTR of 581 bp with a canonical polyadenylation signal sequence AATAAA and a poly( A) tail, and an open reading frame of 576 bp encoding a polypeptide with a molecular mass of 22.21 kDa and a predicted isoelectric point of 9.69. There was an alpha-crystallin domain, a hallmark of the sHSP subfamily, in the C-terminus, and the deduced amino acid sequence of CfHSP22 showed high similarity to previously identified HSP22s. CfHSP22 was constitutively expressed in the haemocyte, muscle, kidney, gonad, gill, heart and hepatopancreas, and the expression level in the hepatopancreas was higher than that in the other tissues. CfHSP22 transcription was up-regulated and reached a maximal level at 12 h after the bacterial challenge, and then declined progressively to the original level at 48 h. These results suggested that CfHSP22 perhaps play a critical role in response to the bacterial challenge in haemocytes of scallop C. farreri.

A close phylogenetic relationship between Sipuncula and Annelida evidenced from the complete mitochondrial genome sequence of Phascolosoma esculenta

Background: There are many advantages to the application of complete mitochondrial (mt) genomes in the accurate reconstruction of phylogenetic relationships in Metazoa. Although over one thousand metazoan genomes have been sequenced, the taxonomic sampling is highly biased, left with many phyla without a single representative of complete mitochondrial genome. Sipuncula (peanut worms or star worms) is a small taxon of worm-like marine organisms with an uncertain phylogenetic position. In this report, we present the mitochondrial genome sequence of Phascolosoma esculenta, the first complete mitochondrial genome of the phylum. Results: The mitochondrial genome of P. esculenta is 15,494 bp in length. The coding strand consists of 32.1% A, 21.5% C, 13.0% G, and 33.4% T bases (AT = 65.5%; AT skew = -0.019; GC skew = -0.248). It contains thirteen protein-coding genes (PCGs) with 3,709 codons in total, twenty-two transfer RNA genes, two ribosomal RNA genes and a non-coding AT-rich region (AT = 74.2%). All of the 37 identified genes are transcribed from the same DNA strand. Compared with the typical set of metazoan mt genomes, sipunculid lacks trnR but has an additional trnM. Maximum Likelihood and Bayesian analyses of the protein sequences show that Myzostomida, Sipuncula and Annelida (including echiurans and pogonophorans) form a monophyletic group, which supports a closer relationship between Sipuncula and Annelida than with Mollusca, Brachiopoda, and some other lophotrochozoan groups. Conclusion: This is the first report of a complete mitochondrial genome as a representative within the phylum Sipuncula. It shares many more similar features with the four known annelid and one echiuran mtDNAs. Firstly, sipunculans and annelids share quite similar gene order in the mitochondrial genome, with all 37 genes located on the same strand; secondly, phylogenetic analyses based on the concatenated protein sequences also strongly support the sipunculan + annelid clade (including echiurans and pogonophorans). Hence annelid "key-characters" including segmentation may be more labile than previously assumed.

Rapid protein identification using monolithic enzymatic microreactor and LC-ESI-MS/MS

A monolithic enzymatic microreactor was prepared in a fused-silica capillary by in situ polymerization of acrylamide, glycidyl methacrylate (GMA) and ethylene dimethacrylate (EDMA) in the presence of a binary porogenic mixture of dodecanol and cyclohexanol, followed by ammonia solution treatment, glutaraldehyde activation and trypsin modification. The choice of acrylamide as co-monomer was found useful to improve the efficiency of trypsin modification, thus, to increase the enzyme activity. The optimized microreactor offered very low back pressure, enabling the fast digestion of proteins flowing through the reactor. The performance of the monolithic microreactor was demonstrated with the digestion of cytochrome c at high flow rate. The digests were then characterized by CE and HPLC-MS/MS with the sequence coverage of 57.7%. The digestion efficiency was found over 230 times as high as that of the conventional method. in addition, for the first time, protein digestion carried out in a mixture of water and ACN was compared with the conventional aqueous reaction using MS/MS detection, and the former solution was found more compatible and more efficient for protein digestion.

Confirmation of Data Mining Based Predictions of Protein Function

King, R. D. and Wise, P. H. and Clare, A. (2004) Confirmation of Data Mining Based Predictions of Protein Function. Bioinformatics 20(7), 1110-1118

Homology Induction: the use of machine learning to improve sequence similarity searches

Karwath, A. King, R. Homology induction: the use of machine learning to improve sequence similarity searches. BMC Bioinformatics. 23rd April 2002. 3:11 Additional File Describes the title organims species declaration in one string [http://www.biomedcentral.com/content/supplementary/1471- 2105-3-11-S1.doc] Sponsorship: Andreas Karwath and Ross D. King were supported by the EPSRC grant GR/L62849.

A Phylogenetic Mixture Model for Detecting Pattern-Heterogeneity in Gene Sequence or Character-State Data

Mark Pagel, Andrew Meade (2004). A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data. Systematic Biology, 53(4), 571-581. RAE2008

A binding site for the cyclic adenosine 3',5'-monophosphate-response element-binding protein as a regulatory element in the grp78 promoter.

The 78-kDa glucose-regulated protein (GRP78) is ubiquitously expressed in many cell types. Its promoter contains multiple protein-binding sites and functional elements. In this study we examined a high affinity protein-binding site spanning bp -198 to -180 of the rat grp78 promoter, using nuclear extracts from both B-lymphoid and HeLa cells. This region contains a sequence TGACGTGA which, with the exception of one base, is identical to the cAMP-response element (CRE). Site-directed mutagenesis reveals that this sequence functions as a major basal level regulatory element in hamster fibroblast cells and is also necessary to maintain high promoter activity under stress-induced conditions. By gel mobility shift analysis, we detect two specific protein complexes. The major specific complex I, while immunologically distinct from the 42-kDa CRE-binding protein (CREB), binds most strongly to the grp site, but also exhibits affinity for the CRE consensus sequence. As such, complex I may consist of other members of the CREB/activating transcription factor protein family. The minor specific complex II consists of CREB or a protein antigenically related to it. A nonspecific complex III consists of the Ku autoantigen, an abundant 70- to 80-kDa protein complex in HeLa nuclear extracts. By cotransfection experiments, we demonstrate that in F9 teratocarcinoma cells, the grp78 promoter can be transactivated by the phosphorylated CREB or when the CREB-transfected cells are treated with the calcium ionophore A23187. The differential regulation of the grp78 gene by cAMP in specific cell types and tissues is discussed.

Mass spectrometry-based thermal shift assay for protein-ligand binding analysis.

Described here is a mass spectrometry-based screening assay for the detection of protein-ligand binding interactions in multicomponent protein mixtures. The assay utilizes an oxidation labeling protocol that involves using hydrogen peroxide to selectively oxidize methionine residues in proteins in order to probe the solvent accessibility of these residues as a function of temperature. The extent to which methionine residues in a protein are oxidized after specified reaction times at a range of temperatures is determined in a MALDI analysis of the intact proteins and/or an LC-MS analysis of tryptic peptide fragments generated after the oxidation reaction is quenched. Ultimately, the mass spectral data is used to construct thermal denaturation curves for the detected proteins. In this proof-of-principle work, the protocol is applied to a four-protein model mixture comprised of ubiquitin, ribonuclease A (RNaseA), cyclophilin A (CypA), and bovine carbonic anhydrase II (BCAII). The new protocol's ability to detect protein-ligand binding interactions by comparing thermal denaturation data obtained in the absence and in the presence of ligand is demonstrated using cyclosporin A (CsA) as a test ligand. The known binding interaction between CsA and CypA was detected using both the MALDI- and LC-MS-based readouts described here.

Interaction of Cryptococcus neoformans Rim101 and protein kinase A regulates capsule.

Cryptococcus neoformans is a prevalent human fungal pathogen that must survive within various tissues in order to establish a human infection. We have identified the C. neoformans Rim101 transcription factor, a highly conserved pH-response regulator in many fungal species. The rim101 multiply sign in circle mutant strain displays growth defects similar to other fungal species in the presence of alkaline pH, increased salt concentrations, and iron limitation. However, the rim101 multiply sign in circle strain is also characterized by a striking defect in capsule, an important virulence-associated phenotype. This capsular defect is likely due to alterations in polysaccharide attachment to the cell surface, not in polysaccharide biosynthesis. In contrast to many other C. neoformans capsule-defective strains, the rim101 multiply sign in circle mutant is hypervirulent in animal models of cryptococcosis. Whereas Rim101 activation in other fungal species occurs through the conserved Rim pathway, we demonstrate that C. neoformans Rim101 is also activated by the cAMP/PKA pathway. We report here that C. neoformans uses PKA and the Rim pathway to regulate the localization, activation, and processing of the Rim101 transcription factor. We also demonstrate specific host-relevant activating conditions for Rim101 cleavage, showing that C. neoformans has co-opted conserved signaling pathways to respond to the specific niche within the infected host. These results establish a novel mechanism for Rim101 activation and the integration of two conserved signaling cascades in response to host environmental conditions.

The G-protein-coupled receptor kinases beta ARK1 and beta ARK2 are widely distributed at synapses in rat brain.

The beta-adrenergic receptor kinase (beta ARK) phosphorylates the agonist-occupied beta-adrenergic receptor to promote rapid receptor uncoupling from Gs, thereby attenuating adenylyl cyclase activity. Beta ARK-mediated receptor desensitization may reflect a general molecular mechanism operative on many G-protein-coupled receptor systems and, particularly, synaptic neurotransmitter receptors. Two distinct cDNAs encoding beta ARK isozymes were isolated from rat brain and sequenced. The regional and cellular distributions of these two gene products, termed beta ARK1 and beta ARK2, were determined in brain by in situ hybridization and by immunohistochemistry at the light and electron microscopic levels. The beta ARK isozymes were found to be expressed primarily in neurons distributed throughout the CNS. Ultrastructurally, beta ARK1 and beta ARK2 immunoreactivities were present both in association with postsynaptic densities and, presynaptically, with axon terminals. The beta ARK isozymes have a regional and subcellular distribution consistent with a general role in the desensitization of synaptic receptors.

Independent beta-arrestin 2 and G protein-mediated pathways for angiotensin II activation of extracellular signal-regulated kinases 1 and 2.

Stimulation of a mutant angiotensin type 1A receptor (DRY/AAY) with angiotensin II (Ang II) or of a wild-type receptor with an Ang II analog ([sarcosine1,Ile4,Ile8]Ang II) fails to activate classical heterotrimeric G protein signaling but does lead to recruitment of beta-arrestin 2-GFP and activation of extracellular signal-regulated kinases 1 and 2 (ERK1/2) (maximum stimulation approximately 50% of wild type). This G protein-independent activation of mitogen-activated protein kinase is abolished by depletion of cellular beta-arrestin 2 but is unaffected by the PKC inhibitor Ro-31-8425. In parallel, stimulation of the wild-type angiotensin type 1A receptor with Ang II robustly stimulates ERK1/2 activation with approximately 60% of the response blocked by the PKC inhibitor (G protein dependent) and the rest of the response blocked by depletion of cellular beta-arrestin 2 by small interfering RNA (beta-arrestin dependent). These findings imply the existence of independent G protein- and beta-arrestin 2-mediated pathways leading to ERK1/2 activation and the existence of distinct "active" conformations of a seven-membrane-spanning receptor coupled to each.

G protein signaling and vein graft intimal hyperplasia: reduction of intimal hyperplasia in vein grafts by a Gbetagamma inhibitor suggests a major role of G protein signaling in lesion development.

Vein grafting results in the development of intimal hyperplasia with accompanying changes in guanine nucleotide-binding (G) protein expression and function. Several serum mitogens that act through G protein-coupled receptors, such as lysophosphatidic acid, stimulate proliferative pathways that are dependent on the G protein betagamma subunit (Gbetagamma)-mediated activation of p21ras. This study examines the role of Gbetagamma signaling in intimal hyperplasia by targeting a gene encoding a specific Gbetagamma inhibitor in an experimental rabbit vein graft model. This inhibitor, the carboxyl terminus of the beta-adrenergic receptor kinase (betaARK(CT)), contains a Gbetagamma-binding domain. Vein graft intimal hyperplasia was significantly reduced by 37% (P<0.01), and physiological studies demonstrated that the normal alterations in G protein coupling phenotypically seen in this model were blocked by betaARK(CT) treatment. Thus, it appears that Gbetagamma-mediated pathways play a major role in intimal hyperplasia and that targeting inhibitors of Gbetagamma signaling offers novel intraoperative therapeutic modalities to inhibit the development of vein graft intimal hyperplasia and subsequent vein graft failure.

beta2-Adrenergic receptor regulation by GIT1, a G protein-coupled receptor kinase-associated ADP ribosylation factor GTPase-activating protein.

G protein-coupled receptor activation leads to the membrane recruitment and activation of G protein-coupled receptor kinases, which phosphorylate receptors and lead to their inactivation. We have identified a novel G protein-coupled receptor kinase-interacting protein, GIT1, that is a GTPase-activating protein (GAP) for the ADP ribosylation factor (ARF) family of small GTP-binding proteins. Overexpression of GIT1 leads to reduced beta2-adrenergic receptor signaling and increased receptor phosphorylation, which result from reduced receptor internalization and resensitization. These cellular effects of GIT1 require its intact ARF GAP activity and do not reflect regulation of GRK kinase activity. These results suggest an essential role for ARF proteins in regulating beta2-adrenergic receptor endocytosis. Moreover, they provide a mechanism for integration of receptor activation and endocytosis through regulation of ARF protein activation by GRK-mediated recruitment of the GIT1 ARF GAP to the plasma membrane.

G-protein-coupled receptor genes as protooncogenes: constitutively activating mutation of the alpha 1B-adrenergic receptor enhances mitogenesis and tumorigenicity.

The alpha 1B-adrenergic receptor (alpha 1B-ADR) is a member of the G-protein-coupled family of transmembrane receptors. When transfected into Rat-1 and NIH 3T3 fibroblasts, this receptor induces focus formation in an agonist-dependent manner. Focus-derived, transformed fibroblasts exhibit high levels of functional alpha 1B-ADR expression, demonstrate a catecholamine-induced enhancement in the rate of cellular proliferation, and are tumorigenic when injected into nude mice. Induction of neoplastic transformation by the alpha 1B-ADR, therefore, identifies this normal cellular gene as a protooncogene. Mutational alteration of this receptor can lead to activation of this protooncogene, resulting in an enhanced ability of agonist to induce focus formation with a decreased latency and quantitative increase in transformed foci. In contrast to cells expressing the wild-type alpha 1B-ADR, focus formation in "oncomutant"-expressing cell lines appears constitutively activated with the generation of foci in unstimulated cells. Further, these cell lines exhibit near-maximal rates of proliferation even in the absence of catecholamine supplementation. They also demonstrate an enhanced ability for tumor generation in nude mice with a decreased period of latency compared with cells expressing the wild-type receptor. Thus, the alpha 1B-ADR gene can, when overexpressed and activated, function as an oncogene inducing neoplastic transformation. Mutational alteration of this receptor gene can result in the activation of this protooncogene, enhancing its oncogenic potential. These findings suggest that analogous spontaneously occurring mutations in this class of receptor proteins could play a key role in the induction or progression of neoplastic transformation and atherosclerosis.