A novel C1q-domain-containing protein from Zhikong scallop Chlamys farreri with lipopolysaccharide binding activity

The C1q-domain-containing (C1qDC) proteins are a family of proteins characterized by a globular C1q (gC1q) domain in their C-terminus. They are involved in various processes of vertebrates and supposed to be an important pattern recognition receptor in innate immunity of invertebrates. In this study, a novel member of C1q-domain-containing protein family was identified from Zhikong scallop Chlamys farreri (designated as CfC1qDC) by expressed sequence tag (EST) and rapid amplification of cDNA ends (RACE) approaches. The full-length cDNA of CfC1qDC was of 777 bp, consisting of a T-terminal untranslated region (UTR) of 62 bp and a 3' UTR of 178 bp with a polyadenylation signal sequence AATAAA and a poly (A) tail. The CfC1qDC cDNA encoded a polypeptide of 178 amino acids, including a signal peptide and a C1q-domain of 158 amino acids with the theoretical isoelectric point of 5.19 and the predicted molecular weight of 17.2 kDa. The C1q-domain in CfC1qDC exhibited homology with those in sialic acid binding lectin from mollusks and C1qDC proteins from higher vertebrates. The typical 10 beta-strand jelly-roll folding topology structure of C1q-domain and the residues essential for effective packing of the hydrophobic core were well conserved in CfC1qDC. By fluorescent quantitative real-time PCR, mRNA transcripts of CfC1qDC were mainly detected in kidney, mantle, adductor muscle and gill, and also marginally detectable in hemocytes. In the bacterial challenge experiment, after the scallops were challenged by Listonella anguillarum, there was a significant up-regulation in the relative expression level of CfC1qDC and at 6 h post-injection, the mRNA expression reached the maximum level and was 4.55-fold higher than that of control scallops. Similarly, the expression of CfC1qDC mRNA in mixed primary cultures of hemocytes stimulated by lipopolysaccharides (LPS) was up-regulated and reached the maximum level at 6 h post-stimulation, and then dropped back to the original level gradually. In order to investigate its function, the cDNA fragment encoding the mature peptide of CfC1qDC was recombined and expressed in Escherichia coli BL21 (DE3). The recombinant CfC1qDC protein displayed a significantly strong activity to bind LIDS from E. coli, although no obvious antibacterial or agglutinating activity toward Gram-negative bacteria E. coli JM109, L. anguillarum and Gram-positive bacteria Micrococcus luteus was observed. These results suggested that CfC1qDC was absolutely a novel member of the C1qDC protein family and was involved in the recognition of invading microorganisms probably as a pattern recognition molecule in mollusk. (c) 2008 Elsevier Ltd. All rights reserved.

Molecular cloning of an invertebrate goose-type lysozyme gene from Chlamys farreri, and lytic activity of the recombinant protein

Lysozyme is a widely distributed hydrolase possessing lytic activity against bacterial peptidoglycan, which enables it to protect the host against pathogenic infection. In the present study, the cDNA of an invertebrate goose-type lysozyme (designated CFLysG) was cloned from Zhikong scallop Chlamys farreri by expressed sequence tag (EST) and rapid amplification of cDNA ends (RACE) techniques. The full-length cDNA of CFLysG consisted of 829 nucleotides with a canonical polyadenylation signal sequence AATAAA and a poly(A) tail, and an open reading frame (ORF) of 603 bp encoding a polypeptide of 200 amino acid residues with a predicted molecular weight of 21.92 kDa and theoretical isoelectric point of 7.76. The high similarity of CFLysG with goose-type (g-type) lysozymes in vertebrate indicated that CFLysG should be an invertebrate counterpart of g-type lysozyme family, which suggested that the origin of g-type lysozyme preceded the emergence of urochordates and even preceded the emergence of deuterostomes. Similar to most g-type lysozymes, CFLysG possessed all conserved features critical for the fundamental structure and function of g-type lysozymes, such as three catalytic residues (Glu 82, Asp 97, Asp 108). By Northern blot analysis, mRNA transcript of CFLysG was found to be most abundantly expressed in the tissues of gills, hepatopancreas and gonad, weakly expressed in the tissues of haemocytes and mantle, while undetectable in the adductor muscle. These results suggested that CFLysG could possess combined features of both the immune and digestive adaptive lysozymes. To gain insight into the in vitro lytic activities of CFLysG, the mature peptide coding region was cloned into Pichia pastoris for heterogeneous expression. Recombinant CFLysG showed inhibitive effect on the growth of both Gram-positive and Gram-negative bacteria with more potent activities against Gram-positive bacteria, which indicated the involvement of CFLysG in the innate immunity of C. farreri. (c) 2006 Elsevier Ltd. All rights reserved.

SmC2P1, a C2 domain protein from Scophthalmus maximus that binds bacterial pathogen-infected lymphocytes and reduces bacterial survival

C2 domains are protein structural modules found in many eukaryotic proteins involved in signal transduction, membrane trafficking, and immune defense. Most of the studied C2 domain-containing proteins are multi-domained in structure, in which the C2 domain is an independently folded motif and plays an essential role in calcium-dependent membrane-targeting. Although C2 domains isolated from intact proteins have been studied for biological functions, no study on natural proteins containing C2 domain only has been documented. In this study, we identified a Scophthalmus maximus protein SmC2P1 that is comprised of a single C2 domain and lacks any other apparent domain structures. The deduced amino acid sequence of SmC2P1 contains 129 residues and shares 36-38% identities with the C2 domains of the perforins of several fish species. Like typical C2 domains, SmC2P1 is predicted to organize into eight beta-strands with a Ca2+-binding site located in inter-strand loops. SmC2P1 expression was detected, in deceasing order, in liver, spleen, blood, brain, muscle, kidney, gill, and heart. Experimental challenge of turbot with a bacterial pathogen significantly upregulated SmC2P1 expression in kidney in a time-dependent manner. Recombinant SmC2P1 purified from yeast exhibits no hemolytic activity but binds to pathogen-infected kidney lymphocytes in the presence of calcium. Furthermore, interaction of recombinant SmC2P1 with bacterium-infected lymphocytes reduced bacterial survival. These results indicate that SmC2P1 is a functional protein that is involved in host immune defense against bacterial infection. (C) 2010 Elsevier Ltd. All rights reserved.

Synthesizing Regularity Exposing Attributes in Large Protein Databases

This thesis describes a system that synthesizes regularity exposing attributes from large protein databases. After processing primary and secondary structure data, this system discovers an amino acid representation that captures what are thought to be the three most important amino acid characteristics (size, charge, and hydrophobicity) for tertiary structure prediction. A neural network trained using this 16 bit representation achieves a performance accuracy on the secondary structure prediction problem that is comparable to the one achieved by a neural network trained using the standard 24 bit amino acid representation. In addition, the thesis describes bounds on secondary structure prediction accuracy, derived using an optimal learning algorithm and the probably approximately correct (PAC) model.

Cloning and structure analysis of hydrogenase gene from Chlamydomonas reinhardtii SE

The cDNA of Chlamydomonas reinhardtii SE encoding hydrogenase (HydA2) was obtained from the total RNA of C reinhardtii SE by RT-PCR. The DNA of hydrogenase was amplified by PCR from the genomic DNA of C reinhardtii SE. The cDNA and DNA of hydrogenase were sequenced, respectively. The structure of hydrogenase gene was analyzed by biology software. The open reading frame predicts that the hydrogenase is composed of 3584 bp encoding 505 amino acids in length with a predicted M.W. of 53.69 kDa. Ten exons (including 1518 bp) and nine introns (including 2066 bp) have been found in the hydrogenase, and there were two potential N-glycosylate sites, eight protein kinase C phosphorylation site, eight casein kinase H phosphorylation site and one sulphorylation in the sequence. The theory pI was 6.15. Total number of negatively charged residues (Asp + Glu) and positively charged residues (Arg + Lys) were 55 and 61, respectively. (c) 2005 Elsevier Ltd. All rights reserved.

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

Application of Inductive Logic Programming to Structure-Based Drug Design

Enot, D. and King, R. D. (2003) Application of Inductive Logic Programming to Structure-Based Drug Design. 7th European Conference on Principles and Practice of Knowledge Discovery in Databases (PKDD '03). Springer LNAI 2838 p156-167

Rubisco small subunit, chlorophyll a/b-binding protein and sucrose : fructan-6-fructosyl transferase gene expression and sugar status in single barley leaf cells in situ. Cell type specificity and induction by light

Chungui Lu, Olga A. Koroleva, John F. Farrar, Joe Gallagher, Chris J. Pollock, and A. Deri Tomos (2002). Rubisco small subunit, chlorophyll a/b-binding protein and sucrose : fructan-6-fructosyl transferase gene expression and sugar status in single barley leaf cells in situ. Cell type specificity and induction by light. Plant Physiology, 130 (3) pp.1335-1348 Sponsorship: BBSRC RAE2008

Integration of relational and hierarchical network information for protein function prediction

BACKGROUND:In the current climate of high-throughput computational biology, the inference of a protein's function from related measurements, such as protein-protein interaction relations, has become a canonical task. Most existing technologies pursue this task as a classification problem, on a term-by-term basis, for each term in a database, such as the Gene Ontology (GO) database, a popular rigorous vocabulary for biological functions. However, ontology structures are essentially hierarchies, with certain top to bottom annotation rules which protein function predictions should in principle follow. Currently, the most common approach to imposing these hierarchical constraints on network-based classifiers is through the use of transitive closure to predictions.RESULTS:We propose a probabilistic framework to integrate information in relational data, in the form of a protein-protein interaction network, and a hierarchically structured database of terms, in the form of the GO database, for the purpose of protein function prediction. At the heart of our framework is a factorization of local neighborhood information in the protein-protein interaction network across successive ancestral terms in the GO hierarchy. We introduce a classifier within this framework, with computationally efficient implementation, that produces GO-term predictions that naturally obey a hierarchical 'true-path' consistency from root to leaves, without the need for further post-processing.CONCLUSION:A cross-validation study, using data from the yeast Saccharomyces cerevisiae, shows our method offers substantial improvements over both standard 'guilt-by-association' (i.e., Nearest-Neighbor) and more refined Markov random field methods, whether in their original form or when post-processed to artificially impose 'true-path' consistency. Further analysis of the results indicates that these improvements are associated with increased predictive capabilities (i.e., increased positive predictive value), and that this increase is consistent uniformly with GO-term depth. Additional in silico validation on a collection of new annotations recently added to GO confirms the advantages suggested by the cross-validation study. Taken as a whole, our results show that a hierarchical approach to network-based protein function prediction, that exploits the ontological structure of protein annotation databases in a principled manner, can offer substantial advantages over the successive application of 'flat' network-based methods.

Tyrosine phosphorylation of G protein alpha subunits by pp60c-src.

A number of lines of evidence suggest that cross-talk exists between the cellular signal transduction pathways involving tyrosine phosphorylation catalyzed by members of the pp60c-src kinase family and those mediated by guanine nucleotide regulatory proteins (G proteins). In this study, we explore the possibility that direct interactions between pp60c-src and G proteins may occur with functional consequences. Preparations of pp60c-src isolated by immunoprecipitation phosphorylate on tyrosine residues the purified G-protein alpha subunits (G alpha) of several heterotrimeric G proteins. Phosphorylation is highly dependent on G-protein conformation, and G alpha(GDP) uncomplexed by beta gamma subunits appears to be the preferred substrate. In functional studies, phosphorylation of stimulatory G alpha (G alpha s) modestly increases the rate of binding of guanosine 5'-[gamma-[35S]thio]triphosphate to Gs as well as the receptor-stimulated steady-state rate of GTP hydrolysis by Gs. Heterotrimeric G proteins may represent a previously unappreciated class of potential substrates for pp60c-src.

The receptor kinase family: primary structure of rhodopsin kinase reveals similarities to the beta-adrenergic receptor kinase.

Light-dependent deactivation of rhodopsin as well as homologous desensitization of beta-adrenergic receptors involves receptor phosphorylation that is mediated by the highly specific protein kinases rhodopsin kinase (RK) and beta-adrenergic receptor kinase (beta ARK), respectively. We report here the cloning of a complementary DNA for RK. The deduced amino acid sequence shows a high degree of homology to beta ARK. In a phylogenetic tree constructed by comparing the catalytic domains of several protein kinases, RK and beta ARK are located on a branch close to, but separate from the cyclic nucleotide-dependent protein kinase and protein kinase C subfamilies. From the common structural features we conclude that both RK and beta ARK are members of a newly delineated gene family of guanine nucleotide-binding protein (G protein)-coupled receptor kinases that may function in diverse pathways to regulate the function of such receptors.

Novel Computational Protein Design Algorithms with Applications to Cystic Fibrosis and HIV

Proteins are essential components of cells and are crucial for catalyzing reactions, signaling, recognition, motility, recycling, and structural stability. This diversity of function suggests that nature is only scratching the surface of protein functional space. Protein function is determined by structure, which in turn is determined predominantly by amino acid sequence. Protein design aims to explore protein sequence and conformational space to design novel proteins with new or improved function. The vast number of possible protein sequences makes exploring the space a challenging problem.

Computational structure-based protein design (CSPD) allows for the rational design of proteins. Because of the large search space, CSPD methods must balance search accuracy and modeling simplifications. We have developed algorithms that allow for the accurate and efficient search of protein conformational space. Specifically, we focus on algorithms that maintain provability, account for protein flexibility, and use ensemble-based rankings. We present several novel algorithms for incorporating improved flexibility into CSPD with continuous rotamers. We applied these algorithms to two biomedically important design problems. We designed peptide inhibitors of the cystic fibrosis agonist CAL that were able to restore function of the vital cystic fibrosis protein CFTR. We also designed improved HIV antibodies and nanobodies to combat HIV infections.

Statistical analysis of crystallization database links protein physico-chemical features with crystallization mechanisms.

X-ray crystallography is the predominant method for obtaining atomic-scale information about biological macromolecules. Despite the success of the technique, obtaining well diffracting crystals still critically limits going from protein to structure. In practice, the crystallization process proceeds through knowledge-informed empiricism. Better physico-chemical understanding remains elusive because of the large number of variables involved, hence little guidance is available to systematically identify solution conditions that promote crystallization. To help determine relationships between macromolecular properties and their crystallization propensity, we have trained statistical models on samples for 182 proteins supplied by the Northeast Structural Genomics consortium. Gaussian processes, which capture trends beyond the reach of linear statistical models, distinguish between two main physico-chemical mechanisms driving crystallization. One is characterized by low levels of side chain entropy and has been extensively reported in the literature. The other identifies specific electrostatic interactions not previously described in the crystallization context. Because evidence for two distinct mechanisms can be gleaned both from crystal contacts and from solution conditions leading to successful crystallization, the model offers future avenues for optimizing crystallization screens based on partial structural information. The availability of crystallization data coupled with structural outcomes analyzed through state-of-the-art statistical models may thus guide macromolecular crystallization toward a more rational basis.

Conformational kinetics reveals affinities of protein conformational states.

Most biological reactions rely on interplay between binding and changes in both macromolecular structure and dynamics. Practical understanding of this interplay requires detection of critical intermediates and determination of their binding and conformational characteristics. However, many of these species are only transiently present and they have often been overlooked in mechanistic studies of reactions that couple binding to conformational change. We monitored the kinetics of ligand-induced conformational changes in a small protein using six different ligands. We analyzed the kinetic data to simultaneously determine both binding affinities for the conformational states and the rate constants of conformational change. The approach we used is sufficiently robust to determine the affinities of three conformational states and detect even modest differences in the protein's affinities for relatively similar ligands. Ligand binding favors higher-affinity conformational states by increasing forward conformational rate constants and/or decreasing reverse conformational rate constants. The amounts by which forward rate constants increase and reverse rate constants decrease are proportional to the ratio of affinities of the conformational states. We also show that both the affinity ratio and another parameter, which quantifies the changes in conformational rate constants upon ligand binding, are strong determinants of the mechanism (conformational selection and/or induced fit) of molecular recognition. Our results highlight the utility of analyzing the kinetics of conformational changes to determine affinities that cannot be determined from equilibrium experiments. Most importantly, they demonstrate an inextricable link between conformational dynamics and the binding affinities of conformational states.

Modulation of bacterial outer membrane vesicle production by envelope structure and content.

BACKGROUND: Vesiculation is a ubiquitous secretion process of Gram-negative bacteria, where outer membrane vesicles (OMVs) are small spherical particles on the order of 50 to 250 nm composed of outer membrane (OM) and lumenal periplasmic content. Vesicle functions have been elucidated in some detail, showing their importance in virulence factor secretion, bacterial survival, and biofilm formation in pathogenesis. Furthermore, OMVs serve as an envelope stress response, protecting the secreting bacteria from internal protein misfolding stress, as well as external envelope stressors. Despite their important functional roles very little is known about the regulation and mechanism of vesicle production. Based on the envelope architecture and prior characterization of the hypervesiculation phenotypes for mutants lacking the lipoprotein, Lpp, which is involved in the covalent OM-peptidoglycan (PG) crosslinks, it is expected that an inverse relationship exists between OMV production and PG-crosslinked Lpp. RESULTS: In this study, we found that subtle modifications of PG remodeling and crosslinking modulate OMV production, inversely correlating with bound Lpp levels. However, this inverse relationship was not found in strains in which OMV production is driven by an increase in "periplasmic pressure" resulting from the accumulation of protein, PG fragments, or lipopolysaccharide. In addition, the characterization of an nlpA deletion in backgrounds lacking either Lpp- or OmpA-mediated envelope crosslinks demonstrated a novel role for NlpA in envelope architecture. CONCLUSIONS: From this work, we conclude that OMV production can be driven by distinct Lpp concentration-dependent and Lpp concentration-independent pathways.