A proteome-integrated, carbon source dependent genetic regulatory network in Saccharomyces cerevisiae.

Integrated regulatory networks can be powerful tools to examine and test properties of cellular systems, such as modelling environmental effects on the molecular bioeconomy, where protein levels are altered in response to changes in growth conditions. Although extensive regulatory pathways and protein interaction data sets exist which represent such networks, few have formally considered quantitative proteomics data to validate and extend them. We generate and consider such data here using a label-free proteomics strategy to quantify alterations in protein abundance for S. cerevisiae when grown on minimal media using glucose, galactose, maltose and trehalose as sole carbon sources. Using a high quality-controlled subset of proteins observed to be differentially abundant, we constructed a proteome-informed network, comprising 1850 transcription factor interactions and 37 chaperone interactions, which defines the major changes in the cellular proteome when growing under different carbon sources. Analysis of the differentially abundant proteins involved in the regulatory network pointed to their significant roles in specific metabolic pathways and function, including glucose homeostasis, amino acid biosynthesis, and carbohydrate metabolic process. We noted strong statistical enrichment in the differentially abundant proteome of targets of known transcription factors associated with stress responses and altered carbon metabolism. This shows how such integrated analysis can lend further experimental support to annotated regulatory interactions, since the proteomic changes capture both magnitude and direction of gene expression change at the level of the affected proteins. Overall this study highlights the power of quantitative proteomics to help define regulatory systems pertinent to environmental conditions.

Cross species proteomics.

Proteomics has advanced in leaps and bounds over the past couple of decades. However, the continuing dependency of mass spectrometry-based protein identification on the searching of spectra against protein sequence databases limits many proteomics experiments. If there is no sequenced genome for a given species, then cross species proteomics is required, attempting to identify proteins across the species boundary, typically using the sequenced genome of a closely related species. Unlike sequence searching for homologues, the proteomics equivalent is confounded by small differences in amino acid sequences, leading to large differences in peptide masses; this renders mass matching of peptides and their product ions difficult. Therefore, the phylogenetic distance between the two species and the attendant level of conservation between the homologous proteins play a huge part in determining the extent of protein identification that is possible across the species boundary. In this chapter, we review the cross species challenge itself, as well as various approaches taken to deal with it and the success met with in past studies. This is followed by recommendations of best practice and suggestions to researchers facing this challenge as well as a final section predicting developments, which may help improve cross species proteomics in the future.

Characterising alternate splicing and tissue specific expression in the chicken from ESTs.

Alternate splicing is believed to produce the greatest diversity in transcriptional complexity and function in eukaryotic species. In this study, we present an analysis of alternative splicing events that occur in the chicken, using the recently sequenced genomic sequence and over 580,000 EST sequences mapped back to the genome. A carefully controlled EST-to-genome mapping pipeline is presented, based around the EXONERATE program using the est2genome model, which also considers several quality control steps to filter out erroneous matches. The data is then used to estimate the level of alternate splicing events with respect to Ensembl predicted transcripts. The EST-genome mappings are characterised at the exon level, in order to classify individual splicing events and provide estimates of novel transcripts not currently annotated by the Ensembl genome database. This is the first large scale analysis of this kind in an avian species, and suggests that chicken displays a similar level of alternate splicing as that found in other higher vertebrates such as human and mouse, both in terms of the number of genes that undergo alternate splicing events, and the average number of transcripts produced per gene. The EST data suggests alternate splicing may occur in some 50-60% of the chicken gene set and with an average of around 2.3 transcripts per gene which undergo this process. The EST data is also used to look at gene and transcript usage in the tissues sequenced in embryonic and adult libraries. Genes which display notable biases were analysed in more detail, including twinfilin-2 and embryonic heavy chain myosin. This also highlights several as yet functionally un-annotated genes which appear to be important in embryonic tissues and also undergo alternate splicing events. The analysis also demonstrates some of the difficulties involved in using EST-based data to annotate transcriptional activity in eukaryotic genes, where a broad spectrum of tissues and a large number of sequenced transcripts are required in order to fully characterise alternate splicing and differential expression.

New insights into the transmission biology of urinary schistosomiasis in Zanzibar.

A better understanding of the transmission biology of urinary schistosomiasis in Zanzibar, Tanzania was only possible after the development of molecular DNA markers for identification of Bulinus africanus group snails, the potential intermediate hosts of Schistosoma haematobium. Hitherto, identification of natural populations of B. globosus and B. nasutus was problematic and the intermediate host status and distribution of either species remained speculative. By recourse to molecular markers, snail distribution maps could be drawn, revealing an allopatric distribution and, more importantly, leading to the discovery that B. nasutus played no role in transmission. Indeed, in Unguja the area of active transmission of S. haematobium to humans is confined within the distribution of B. globosus. This strong relationship may prove useful for predicting the distribution of urinary schistosomiasis within Zanzibar and, if snail schistosome compatibilities persist, in other areas nearby, e.g. coastal Tanzania and Kenya. The transmission biology of urinary schistosomiasis in Zanzibar is reviewed, the paper reports on ongoing malacological studies in Zanzibar and Kenya and finally closes by posing the question whether medical malacology forms an essential component associated with mass-scale chemotherapy control programmes.

Identification of snails within the Bulinus africanus group from East Africa by multiplex SNaPshot trade mark analysis of single nucleotide polymorphisms within the cytochrome oxidase subunit I.

Identification of populations of Bulinus nasutus and B. globosus from East Africa is unreliable using characters of the shell. In this paper, a molecular method of identification is presented for each species based on DNA sequence variation within the mitochondrial cytochrome oxidase subunit I (COI) as detected by a novel multiplexed SNaPshotTM assay. In total, snails from 7 localities from coastal Kenya were typed using this assay and variation within shell morphology was compared to reference material from Zanzibar. Four locations were found to contain B. nasutus and 2 locations were found to contain B. globosus. A mixed population containing both B. nasutus and B. globosus was found at Kinango. Morphometric variation between samples was considerable and UPGMA cluster analysis failed to differentiate species. The multiplex SNaPshotTM assay is an important development for more precise methods of identification of B. africanus group snails. The assay could be further broadened for identification of other snail intermediate host species.

Identification of snails within the Bulinus africanus group from East Africa by multiplex SNaPshotTManalysis of single nucleotide polymorphisms within the cytochrome oxidase subunit I

Identification of populations of Bulinus nasutus and B. globosus from East Africa is unreliable using characters of the shell. In this paper, a molecular method of identification is presented for each species based on DNA sequence variation within the mitochondrial cytochrome oxidase subunit I (COI) as detected by a novel multiplexed SNaPshotTM assay. In total, snails from 7 localities from coastal Kenya were typed using this assay and variation within shell morphology was compared to reference material from Zanzibar. Four locations were found to contain B. nasutus and 2 locations were found to contain B. globosus. A mixed population containing both B. nasutus and B. globosus was found at Kinango. Morphometric variation between samples was considerable and UPGMA cluster analysis failed to differentiate species. The multiplex SNaPshotTM assay is an important development for more precise methods of identification of B. africanus group snails. The assay could be further broadened for identification of other snail intermediate host species

A combination of chemical derivatisation and improved bioinformatic tools optimises protein identification for proteomics.

The identification of individual protein species within an organism's proteome has been optimised by increasing the information produced from mass spectral analysis through the chemical derivatisation of tryptic peptides and the development of new software tools. Peptide fragments are subjected to two forms of derivatisation. First, lysine residues are converted to homoarginine moieties by guanidination. This procedure has two advantages, first, it usually identifies the C-terminal amino acid of the tryptic peptide and also greatly increases the total information content of the mass spectrum by improving the signal response of C-terminal lysine fragments. Second, an Edman-type phenylthiocarbamoyl (PTC) modification is carried out on the N-terminal amino acid. The renders the first peptide bond highly susceptible to cleavage during mass spectrometry (MS) analysis and consequently allows the ready identification of the N-terminal residue. The utility of the procedure has been demonstrated by developing novel bioinformatic tools to exploit the additional mass spectral data in the identification of proteome proteins from the yeast Saccharomyces cerevisiae. With this combination of novel chemistry and bioinformatics, it should be possible to identify unambiguously any yeast protein spot or band from either two-dimensional or one-dimensional electropheretograms.

Effect of polymorphisms on ligand binding by mouse major urinary proteins.

Mouse urine contains an abundance of major urinary proteins, lipocalins, whose roles include slow release of semiochemicals. These proteins are highly polymorphic, with small sequence differences between individual members. In this study, we purified to homogeneity four of these proteins from two strains of inbred mice and characterized them by mass spectrometry. This analysis has led to the discovery of another variant in this group of proteins. Three of the polymorphic variants that map to the surface have no effect on the binding of a fluorescent probe in the binding cavity, but the fourth, characterized by a Phe to Val substitution in the cavity, shows a substantially lower affinity and fluorescence yield for the probe. These results are interpreted in light of the known crystal structure of the protein and molecular modeling calculations, which rationalize the experimental findings. This work raises the possibility that the calyx-binding site can show specificity for different ligands, the implications of which on pheromone binding and chemical communication are discussed.

Marked differences between metalloproteases meprin A and B in substrate and peptide bond specificity.

Meprin A and B are highly regulated, secreted, and cell-surface metalloendopeptidases that are abundantly expressed in the kidney and intestine. Meprin oligomers consist of evolutionarily related alpha and/or beta subunits. The work herein was carried out to identify bioactive peptides and proteins that are susceptible to hydrolysis by mouse meprins and kinetically characterize the hydrolysis. Gastrin-releasing peptide fragment 14-27 and gastrin 17, regulatory molecules of the gastrointestinal tract, were found to be the best peptide substrates for meprin A and B, respectively. Peptide libraries and a variety of naturally occurring peptides revealed that the meprin beta subunit has a clear preference for acidic amino acids in the P1 and P1' sites of substrates. The meprin alpha subunit selected for small (e.g. serine, alanine) or hydrophobic (e.g. phenylalanine) residues in the P1 and P1' sites, and proline was the most preferred amino acid at the P2' position. Thus, although the meprin alpha and beta subunits share 55% amino acid identity within the protease domain and are normally localized at the same tissue cell surfaces, they have very different substrate and peptide bond specificities indicating different functions. Homology models of the mouse meprin alpha and beta protease domains, based on the astacin crystal structure, revealed active site differences that can account for the marked differences in substrate specificity of the two subunits.

Bioinformatic assessment of mass spectrometric chemical derivatisation techniques for proteome database searching.

Identification of proteins from the mass spectra of peptide fragments generated by proteolytic cleavage using database searching has become one of the most powerful techniques in proteome science, capable of rapid and efficient protein identification. Using computer simulation, we have studied how the application of chemical derivatisation techniques may improve the efficiency of protein identification from mass spectrometric data. These approaches enhance ion yield and lead to the promotion of specific ions and fragments, yielding additional database search information. The impact of three alternative techniques has been assessed by searching representative proteome databases for both single proteins and simple protein mixtures. For example, by reliably promoting fragmentation of singly-charged peptide ions at aspartic acid residues after homoarginine derivatisation, 82% of yeast proteins can be unambiguously identified from a single typical peptide-mass datum, with a measured mass accuracy of 50 ppm, by using the associated secondary ion data. The extra search information also provides a means to confidently identify proteins in protein mixtures where only limited data are available. Furthermore, the inclusion of limited sequence information for the peptides can compensate and exceed the search efficiency available via high accuracy searches of around 5 ppm, suggesting that this is a potentially useful approach for simple protein mixtures routinely obtained from two-dimensional gels.

Sequence search algorithms for single pass sequence identification: does one size fit all?

Bioinformatic tools have become essential to biologists in their quest to understand the vast quantities of sequence data, and now whole genomes, which are being produced at an ever increasing rate. Much of these sequence data are single-pass sequences, such as sample sequences from organisms closely related to other organisms of interest which have already been sequenced, or cDNAs or expressed sequence tags (ESTs). These single-pass sequences often contain errors, including frameshifts, which complicate the identification of homologues, especially at the protein level. Therefore, sequence searches with this type of data are often performed at the nucleotide level. The most commonly used sequence search algorithms for the identification of homologues are Washington University's and the National Center for Biotechnology Information's (NCBI) versions of the BLAST suites of tools, which are to be found on websites all over the world. The work reported here examines the use of these tools for comparing sample sequence datasets to a known genome. It shows that care must be taken when choosing the parameters to use with the BLAST algorithms. NCBI's version of gapped BLASTn gives much shorter, and sometimes different, top alignments to those found using Washington University's version of BLASTn (which also allows for gaps), when both are used with their default parameters. Most of the differences in performance were found to be due to the choices of default parameters rather than underlying differences between the two algorithms. Washington University's version, used with defaults, compares very favourably with the results obtained using the accurate but computationally intensive Smith-Waterman algorithm.

Conceptual modelling of genomic information.

MOTIVATION: Genome sequencing projects are making available complete records of the genetic make-up of organisms. These core data sets are themselves complex, and present challenges to those who seek to store, analyse and present the information. However, in addition to the sequence data, high throughput experiments are making available distinctive new data sets on protein interactions, the phenotypic consequences of gene deletions, and on the transcriptome, proteome, and metabolome. The effective description and management of such data is of considerable importance to bioinformatics in the post-genomic era. The provision of clear and intuitive models of complex information is surprisingly challenging, and this paper presents conceptual models for a range of important emerging information resources in bioinformatics. It is hoped that these can be of benefit to bioinformaticians as they attempt to integrate genetic and phenotypic data with that from genomic sequences, in order to both assign gene functions and elucidate the different pathways of gene action and interaction. RESULTS: This paper presents a collection of conceptual (i.e. implementation-independent) data models for genomic data. These conceptual models are amenable to (more or less direct) implementation on different computing platforms.

Proteolysis of native proteins. Trapping of a reaction intermediate.

When limited proteolysis of the mouse major urinary proteins by trypsin was stopped by rapid denaturation of the proteinase, a covalent adduct of the two proteins was observed. The formation of this complex required active trypsin, was favored at low pH, and could be reversed by the addition of covalent or non-covalent trypsin inhibitors. Electrospray mass spectrometry of the complex demonstrated that it was an acyl-enzyme complex, formed after an unusual exopeptidase attack on the C-terminal-Arg-Glu-OH sequence by trypsin. The complex could sequester over 50% of the trypsin in a digestion mixture, and as anticipated, the protein was an effective trypsin inhibitor.

Assessment of conformational parameters as predictors of limited proteolytic sites in native protein structures.

Despite the importance of limited proteolysis in biological systems it is often difficult to rationalize why a proteinase hydrolyses a particular bond, given a simple sequence specificity alone. Understanding of the structural properties limiting the proteolysis represents a first step on the pathway to control and manipulation of this phenomena. An expanded set of nick-sites in proteins of known tertiary structure, cut by both narrow and broad specificity proteinases, has been generated yielding a robust data set of strictly limited sites. A critical evaluation of an expanded set of conformational parameters revealed a strong correlation with limited proteolytic sites, although they are only modest predictors in isolation. The overall predictive power is significantly improved when the conformational parameters are combined in a weighted predictive scheme that permits their relative importance to be compared via a Metropolis search protocol. A subset of the parameters performs equally well demonstrating the key determinants of susceptibility. The derived predictive algorithm has been made available via the internet. Its utility for predicting other surface-correlated features is also discussed.

Use of fortuitous in vitro mutations in a synthetic Cu-Zn superoxide dismutase gene to illuminate protein structure-function relationships.

A completely synthetic bovine copper-zinc superoxide dismutase gene (Cu-ZnSOD), designed using the most favoured codons for expression in yeast, was constructed. Fortuitous mutations introduced while cloning the synthetic gene permitted the additional construction of four altered-polypeptide products representing two single (Pro121-->Leu and Gly128-->Asp), one double (Pro100-->Leu, Arg113-->Lys) and one triple (Pro100-->Leu, Arg113-->Lys, Pro121-->Leu) mutant. All five versions of the gene were expressed in a SOD-deficient Escherichia coli strain. The 'wild-type' version of the gene and the two single-mutants were expressed to equal extents (approximately 8% of total soluble protein). However, compared with the 'wild-type' enzyme, one single-mutant (Gly128-->Asp) showed almost twice as much dismutase activity whilst the other (Pro121-->Leu) exhibited only 70% of the 'wild-type' level. In contrast, the double and triple mutants showed diminished expression of the gene (approximately 1 and 3% of total soluble protein, respectively) and almost no detectable SOD activity. Polyclonal antibovine SOD antibody bound all the recombinant proteins, although some of the products showed decreased size and probably altered conformations. The 'wild-type' superoxide dismutase recombinant was correctly dimerized and possessed dismutase activity, as did the Gly128-->Asp mutant despite the change in charge. Mutations in the other three versions affected enzyme folding and activity. The effect of the different mutations appeared to be additive, with the Pro121-->Leu substitution leading to the apparent proteolytic degradation of the enzyme in vivo.

A functional role for protein cavities in domain: domain motions.

Motions between individual domains are known to play an important role in protein function. Protein cavities at domain interfaces have been suggested to facilitate such movements. Consequently, the cavity morphology in a set of multi-domain proteins has been critically examined. The conformational changes were well characterised by atomic resolution tertiary structures prior to and after domain motions. The results showed that interdomain cavities play a number of specific functional roles by either facilitating, or being otherwise involved with, domain: domain motions. Correspondingly, a higher fraction of cavity surface is observed at domain interfaces as compared to that buried within individual domains. Furthermore, interdomain cavity-forming residues were found to be highly conserved in terms of amino acid residue sequence and volume within their aligned protein families, more so than residues exclusive to the domain interface and intradomain cavities. These results provide substantial evidence of cavities fulfilling a specific functional role in multi-domain proteins.

Comparison of atomic solvation parametric sets: applicability and limitations in protein folding and binding.

Atomic solvation parameters (ASP) are widely used to estimate the solvation contribution to the thermodynamic stability of proteins as well as the free energy of association for protein-ligand complexes. They are also included in several molecular mechanics computer programs. In this work, a total of eight atomic solvation parametric sets has been employed to calculate the solvation contribution to the free energy of folding delta Gs for 17 proteins. A linear correlation between delta Gs and the number of residues in each protein was found for each ASP set. The calculations also revealed a great variety in the absolute value and in the sign of delta Gs values such that certain ASP sets predicted the unfolded state to be more stable than the folded, whereas others yield precisely the opposite. Further, the solvation contribution to the free energy of association of helix pairs and to the disassociation of loops (connection between secondary structural elements in proteins) from the protein tertiary structures were computed for each of the eight ASP sets and discrepancies were evident among them.

Detection of internal cavities in globular proteins.

We have undertaken a study of internal cavities in five protein structure groups, each containing different crystallographic structure determinations of the same protein, to understand better the nature of packing defects in protein tertiary architectures. Our results show that cavity detection and consistency of detection are highly dependent on probe and cavity size, cavity position within the globular protein and the local "quality' (r.m.s. deviation) of structural consistency within the group. The consistency of solvent placement within cavities has also been examined. We provide guidelines for estimating the likelihood of a given cavity to be an actual packing defect or to be a result of experimental error.

Evidence on close packing and cavities in proteins.

The packing of a protein's constituent atoms and the attendant constraints placed upon them form the basis of many attempts to understand and predict protein structure, stability, folding and even function. Although the significance of packing is yet to be fully comprehended, recent experimental and theoretical investigations have increased our understanding through the description of mutational effects on structure and stability, determination of the limits of packing constraints for both protein folding and structure prediction, and delineation of packing guidelines on the basis of observed cavities in the native protein folds. These advances and allowing protein modellers, engineers and designers to tackle their problems from a more rational perspective.