Coordinate regulation of energy transduction modules in Halobacterium sp. analyzed by a global systems approach

Dimethyl multiplexed labeling combined with microcolumn separation and MS analysis for time course study in proteomics

Stable-isotope labeling coupled with liquid-phase separation and MS analysis is a powerful technique for comparative proteomics. We developed a dimethyl labeling strategy (Anal. Chem. 2003, 75, 6843-6852 and J. Proteome Res. 2005, 4, 101-108) to label peptide N-terminus and epsilon-amino groups of Lys with water-soluble formaldehyde via reductive methylation, and an isotopic pair of formaldehyde is used for binary labeling on two sets of samples. In this study, this approach is extended to a four sample labeling by combining the binary isotopic reagents of formaldehyde (d0, d2) and the binary isotopic reducing reagents, sodium cyanoborohydride (d0, d3). To ensure sufficient mass difference, this multiplexed labeling is coupled with endoproteinase Lys-C instead of trypsin for digestion, resulting in at least two labeling sites with a mass difference of 4 Da for each pair of peptide digest. Moreover, multiplex dimethyl labeling was proved to have no significant isotopic effect during RP LC elution. This method was further applied for monitoring Lys-C digestion using hemoglobin as a model. Dimethyl labeled digests derived from seven time points (1-30 h) were grouped into two sets of sample mixtures, separated by nano-LC to reduce the complexity, and then analyzed by ESI-MS/MS. The temporal study reveals that Lys-C digestion was completed in 10-15 h for all detected peptides. The multiplex dimethyl method has not only provided a simultaneous detection mean for four sample sets but has also conserved all the advantages associated with the original binary method.

Quantitative Profiling of the Membrane Proteome in a Halophilic Archaeon

We present a large scale quantitation study of the membrane proteome from Halobacterium salinarum. To overcome problems generally encountered with membrane proteins, we established a membrane preparation protocol that allows the application of most proteomic techniques originally developed for soluble proteins. Proteins were quantified using two complementary approaches. For gel-based quantitation, DIGE labeling was combined with two-dimensional gel electrophoresis on an improved 16-benzyldimethyl-n-hexadecylammonium chloride/SDS system. MS-based quantitation was carried out by combining gel-free separation with the recently developed isotope-coded protein labeling technique. Good correlations between these two independent quantitation strategies were obtained. From computational analysis we conclude that labeling of free amino groups by isotope-coded protein labeling (Lys and free N termini) is better suited for membrane proteins than Cys-based labeling strategies but that quantitation of integral membrane proteins remains cumbersome compared with soluble proteins. Nevertheless we could quantify 155 membrane proteins; 101 of these had transmembrane domains. We compared two growth states that strongly affect the energy supply of the cells: aerobic versus anaerobic/phototrophic conditions. The photosynthetic protein bacteriorhodopsin is the most highly regulated protein. As expected, several other membrane proteins involved in aerobic or anaerobic energy metabolism were found to be regulated, but in total, however, the number of regulated proteins is rather small.

Salt stress in Desulfovibrio vulgaris Hildenborough: an integrated genomics approach

The ability of Desulfovibrio vulgaris Hildenborough to reduce, and therefore contain, toxic and radioactive metal waste has made all factors that affect the physiology of this organism of great interest. Increased salinity is an important and frequent fluctuation faced by D. vulgaris in its natural habitat. In liquid culture, exposure to excess salt resulted in striking elongation of D. vulgaris cells. Using data from transcriptomics, proteomics, metabolite assays, phospholipid fatty acid profiling, and electron microscopy, we used a systems approach to explore the effects of excess NaCl on D. vulgaris. In this study we demonstrated that import of osmoprotectants, such as glycine betaine and ectoine, is the primary mechanism used by D. vulgaris to counter hyperionic stress. Several efflux systems were also highly up-regulated, as was the ATP synthesis pathway. Increases in the levels of both RNA and DNA helicases suggested that salt stress affected the stability of nucleic acid base pairing. An overall increase in the level of branched fatty acids indicated that there were changes in cell wall fluidity. The immediate response to salt stress included up-regulation of chemotaxis genes, although flagellar biosynthesis was down-regulated. Other down-regulated systems included lactate uptake permeases and ABC transport systems. The results of an extensive NaCl stress analysis were compared with microarray data from a KCl stress analysis, and unlike many other bacteria, D. vulgaris responded similarly to the two stresses. Integration of data from multiple methods allowed us to develop a conceptual model for the salt stress response in D. vulgaris that can be compared to those in other microorganisms.

Mining the acute respiratory distress syndrome proteome: identification of the insulin-like growth factor (IGF)/IGF-binding protein-3 pathway in acute lung injury

To obtain a more complete protein profile of the airspace milieu in acute respiratory distress syndrome (ARDS) and to identify new mediators, we analyzed bronchoalveolar lavage fluid (BALF) by shotgun proteomics. Using BALF from three patients, we identified a total of 870 different proteins, a nearly 10-fold increase from previous reports. Among the proteins identified were known markers of lung injury, such as surfactant, proteases, and serum proteins. We also identified several biologically interesting proteins not previously identified in patients with ARDS, including insulin-like growth factor-binding protein-3 (IGFBP-3). Because of the known role of IGFBP-3 in regulating cell survival, we measured IGFBP-3 levels by enzyme-linked immunosorbent assay in ARDS BALF. Normal controls had low levels of IGFBP-3, whereas patients with early ARDS had a significant increase in IGFBP-3. The IGF pathway, acting through the type 1 IGF-receptor, repressed apoptosis of lung fibroblasts but not lung epithelial cells. Furthermore, depletion of IGF in ARDS BALF led to enhanced fibroblast apoptosis. Our data suggest that the IGFBP-3/IGF pathway is involved in the pathogenesis of lung injury, illustrating the power of shotgun proteomics to catalog proteins present in complex biological fluids, such as BALF, from which hypotheses can be developed and tested.

Quantitative analysis of 2,3,7,8-tetrachlorodibenzo-p-dioxin-induced proteome alterations in 5L rat hepatoma cells using isotope-coded protein labels

In an effort to contribute to a better understanding of the hepatic toxicity of the ubiquitous environmental pollutant and hepatocarcinogen 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), a comprehensive quantitative proteome analysis was performed on 5L rat hepatoma cells exposed to 1 nM TCDD for 8 h. Changes in the abundances of individual protein species in TCDD-treated cells as compared to untreated cells were analysed using the nongel-based isotope-coded protein label (ICPL) method [Schmidt, A., Kellermann, J., Lottspeich, F., Proteomics 2005, 5, 4-15]. 89 proteins were identified as up- or down-regulated by TCDD. For the majority of the altered proteins, an impact of TCDD on their abundance had not been known before. Due to the physicochemical properties or the translational regulation of a large number of the affected proteins, their alteration would have escaped detection by gel-based methods for proteome analysis and by standard mRNA expression profiling, respectively. The identified proteins with TCDD-altered abundance include several proteins implicated in cell cycle regulation, growth factor signalling and the control of apoptosis. The results thus provide new starting-points for the investigation of specific aspects of the toxicity and carcinogenicity of dioxin in liver.

Integrated biclustering of heterogeneous genome-wide datasets for the inference of global regulatory networks

Background

The learning of global genetic regulatory networks from expression data is a severely under-constrained problem that is aided by reducing the dimensionality of the search space by means of clustering genes into putatively co-regulated groups, as opposed to those that are simply co-expressed. Be cause genes may be co-regulated only across a subset of all observed experimental conditions, biclustering (clustering of genes and conditions) is more appropriate than standard clustering. Co-regulated genes are also often functionally (physically, spatially, genetically, and/or evolutionarily) associated, and such a priori known or pre-computed associations can provide support for appropriately grouping genes. One important association is the presence of one or more common cis-regulatory motifs. In organisms where these motifs are not known, their de novo detection, integrated into the clustering algorithm, can help to guide the process towards more biologically parsimonious solutions.

Results

We have developed an algorithm, cMonkey, that detects putative co-regulated gene groupings by integrating the biclustering of gene expression data and various functional associations with the de novo detection of sequence motifs.

Conclusion

We have applied this procedure to the archaeon Halobacterium NRC-1, as part of our efforts to decipher its regulatory network. In addition, we used cMonkey on public data for three organisms in the other two domains of life: Helicobacter pylori, Saccharomyces cerevisiae, and Escherichia coli. The biclusters detected by cMonkey both recapitulated known biology and enabled novel predictions (some for Halobacterium were subsequently confirmed in the laboratory). For example, it identified the bacteriorhodopsin regulon, assigned additional genes to this regulon with apparently unrelated function, and detected its known promoter motif. We have performed a thorough comparison of cMonkey results against other clustering methods, and find that cMonkey biclusters are more parsimonious with all available evidence for co-regulation.

A systems view of haloarchaeal strategies to withstand stress from transition metals

Given that transition metals are essential cofactors in central biological processes, misallocation of the wrong metal ion to a metalloprotein can have resounding and often detrimental effects on diverse aspects of cellular physiology. Therefore, in an attempt to characterize unique and shared responses to chemically similar metals, we have reconstructed physiological behaviors of Halobacterium NRC-1, an archaeal halophile, in sublethal levels of Mn(II), Fe(II), Co(II), Ni(II), Cu(II), and Zn(II). Over 20% of all genes responded transiently within minutes of exposure to Fe(II), perhaps reflecting immediate large-scale physiological adjustments to maintain homeostasis. At steady state, each transition metal induced growth arrest, attempts to minimize oxidative stress, toxic ion scavenging, increased protein turnover and DNA repair, and modulation of active ion transport. While several of these constitute generalized stress responses, up-regulation of active efflux of Co(II), Ni(II), Cu(II), and Zn(II), down-regulation of Mn(II) uptake and up-regulation of Fe(II) chelation, confer resistance to the respective metals. We have synthesized all of these discoveries into a unified systems-level model to provide an integrated perspective of responses to six transition metals with emphasis on experimentally verified regulatory mechanisms. Finally, through comparisons across global transcriptional responses to different metals, we provide insights into putative in vivo metal selectivity of metalloregulatory proteins and demonstrate that a systems approach can help rapidly unravel novel metabolic potential and regulatory programs of poorly studied organisms.

The Inferelator: an algorithm for learning parsimonious regulatory networks from systems-biology data sets de novo

We present a method (the Inferelator) for deriving genome-wide transcriptional regulatory interactions, and apply the method to predict a large portion of the regulatory network of the archaeon Halobacterium NRC-1. The Inferelator uses regression and variable selection to identify transcriptional influences on genes based on the integration of genome annotation and expression data. The learned network successfully predicted Halobacterium's global expression under novel perturbations with predictive power similar to that seen over training data. Several specific regulatory predictions were experimentally tested and verified.

Post-genomics of the model haloarchaeon Halobacterium sp. NRC-1

Halobacteriumsp. NRC-1 is an extremely halophilic archaeon that is easily cultured and genetically tractable. Since its genome sequence was completed in 2000, a combination of genetic, transcriptomic, proteomic, and bioinformatic approaches have provided insights into both its extremophilic lifestyle as well as fundamental cellular processes common to all life forms. Here, we review post-genomic research on this archaeon, including investigations of DNA replication and repair systems, phototrophic, anaerobic, and other physiological capabilities, acidity of the proteome for function at high salinity, and role of lateral gene transfer in its evolution.

From genomes to function: haloarchaea as model organisms

Haloarchaea are adapted to high-salt environments and accumulate equally high salt concentrations in the cytoplasm. The genomes of representatives of six haloarchaeal genera have been fully or partially sequenced, allowing the analysis of haloarchaeal propertiesin silico. Transcriptome and proteome analyses have been established forHalobacterium salinarumandHaloferax volcanii. Genetic systems are available including methods that allow the fast in-frame deletion or modification of chromosomal genes. The high-efficiency transformation system ofHf. volcaniiallows the isolation of genes essential for a biological process by complementation of loss-of-function mutants. For the analysis of haloarchaeal biology many molecular genetic, biochemical, structural and cell biological methods have been adapted to application at high salt concentrations. Recently it has become clear that several different mechanisms allow the adaptation of proteins to the high salt concentration of the cytoplasm. Taken together, the wealth of techniques available make haloarchaea excellent archaeal model species.

From replication to cultivation: hot news from Haloarchaea

Haloarchaea have developed into model organisms that are utilized to study many biological processes. Examples are the mechanisms of chromosome maintenance, gene expression and its regulation, protein export and degradation, and motility and sensing. In addition to the analysis of model species like Halobacterium salinarum and Haloferax volcanii, natural communities have been characterized. Halophilic Archaea were found in low-salt environments and are thus more widespread than previously thought.

Heat shock response of Archaeoglobus fulgidus

The heat shock response of the hyperthermophilic archaeon Archaeoglobus fulgidus strain VC-16 was studied using whole-genome microarrays. On the basis of the resulting expression profiles, approximately 350 of the 2,410 open reading frames (ORFs) (ca. 14%) exhibited increased or decreased transcript abundance. These span a range of cell functions, including energy production, amino acid metabolism, and signal transduction, where the majority are uncharacterized. One ORF called AF1298 was identified that contains a putative helix-turn-helix DNA binding motif. The gene product, HSR1, was expressed and purified from Escherichia coli and was used to characterize specific DNA recognition regions upstream of two A. fulgidus genes, AF1298 and AF1971. The results indicate that AF1298 is autoregulated and is part of an operon with two downstream genes that encode a small heat shock protein, Hsp20, and cdc48, an AAA+ ATPase. The DNase I footprints using HSR1 suggest the presence of a cis-binding motif upstream of AF1298 consisting of CTAAC-N5-GTTAG. Since AF1298 is negatively regulated in response to heat shock and encodes a protein only distantly related to the N-terminal DNA binding domain of Phr of Pyrococcus furiosus, these results suggest that HSR1 and Phr may belong to an evolutionarily diverse protein family involved in heat shock regulation in hyperthermophilic and mesophilic Archaea organisms.

Integrative genomic and proteomic analysis of prostate cancer reveals signatures of metastatic progression

Molecular profiling of cancer at the transcript level has become routine. Large-scale analysis of proteomic alterations during cancer progression has been a more daunting task. Here, we employed high-throughput immunoblotting in order to interrogate tissue extracts derived from prostate cancer. We identified 64 proteins that were altered in prostate cancer relative to benign prostate and 156 additional proteins that were altered in metastatic disease. An integrative analysis of this compendium of proteomic alterations and transcriptomic data was performed, revealing only 48%-64% concordance between protein and transcript levels. Importantly, differential proteomic alterations between metastatic and clinically localized prostate cancer that mapped concordantly to gene transcripts served as predictors of clinical outcome in prostate cancer as well as other solid tumors.

Global whole-cell FTICR mass spectrometric proteomics analysis of the heat shock response in the radioresistant bacterium Deinococcus radiodurans

The results of previous studies indicated that D. radiodurans mounts a regulated protective response to heat shock, and that expression of more than 130 genes, including classical chaperones such as the groESL and dnaKJ operons and proteases such as clpB are induced in response to elevated temperature. In addition, previous qualitative whole-cell mass spectrometric studies conducted under heat shock conditions indicated global changes in the D. radiodurans proteome. To enable the discovery of novel heat shock inducible proteins as well as gain greater biological insight into the classical heat shock response at the protein level, we undertook the global whole-cell FTICR mass spectrometric proteomics study reported here. We have greatly increased the power of this approach by conducting a large number of replicate experiments in addition to taking a semiquantitative approach to data analysis, finding good reproducibility between replicates. Through this analysis, we have identified with high confidence a core set of classical heat shock proteins whose expression increases dramatically and reproducibly in response to elevated temperature. In addition, we have found that the heat shock proteome includes a large number of induced proteins that have not been identified previously as heat responsive, and have therefore been designated as candidate responders. Finally, our results are consistent with the hypothesis that elevated temperature stress could lead to cross-protection against other related stresses.

The biological limitations of transcriptomics in elucidating stress and stress responses

Global analysis of mRNA abundance via genomic arrays (i.e. transcriptomics or transcriptional profiling) is one approach to finding the genes that matter to organisms undergoing environmental stress. In evolutionary analyses of stress, mRNA abundance is often invoked as a proxy for the protein activity that may underlie variation in fitness. To provoke discussion of the utility and sensible application of this valuable approach, this manuscript examines the adequacy of mRNA abundance as a proxy for protein activity, fitness and stress. Published work to date suggests that mRNA abundance typically provides little information on protein activity and fitness and cannot substitute for detailed functional and ecological analyses of candidate genes. While the transcriptional profile can be an exquisitely sensitive indicator of stress, simpler indicators will often suffice. In view of this outcome, transcriptomics should undergo careful cost-benefit analysis before investigators deploy it in studies of stress responses and their evolution.

Deconvolution of a complex target using DNA aptamers

In vitro selection of single-stranded nucleic acid aptamers from large random sequence libraries is now a straightforward process particularly when screening with a single target molecule. These libraries contain considerable shape diversity as evident by the successful isolation of aptamers that bind with high affinity and specificity to chemically diverse targets. We propose that aptamer libraries contain sufficient shape diversity to allow deconvolution of a complex mixture of targets. Using unfractionated human plasma as our experimental model, we aim to develop methods to obtain aptamers against as many proteins as possible. To begin, it is critical that we understand how aptamer populations change with increasing rounds of in vitro selection when using complex mixtures. Our results show that sequence representation in the selected population changes dramatically with increasing rounds of selection. Certain aptamer families were apparent after only three selection rounds. Two additional cycles saw a decline in the relative abundance of these families and the emergence of yet another family that accounted for more than 60% of sequences in the pool. To overcome this population convergence, an aptamer-based target depletion method was developed, and the library screen was repeated. The previous dominant family effectively disappeared from the selected populations but was replaced by other aptamer families. Insights gained from these initial experiments are now being applied in the creation of second generation plasma protein screens and also to the analysis of other complex biological targets.

A uniform proteomics MS/MS analysis platform utilizing open XML file formats

The analysis of tandem mass (MS/MS) data to identify and quantify proteins is hampered by the heterogeneity of file formats at the raw spectral data, peptide identification, and protein identification levels. Different mass spectrometers output their raw spectral data in a variety of proprietary formats, and alternative methods that assign peptides to MS/MS spectra and infer protein identifications from those peptide assignments each write their results in different formats. Here we describe an MS/MS analysis platform, the Trans-Proteomic Pipeline, which makes use of open XML file formats for storage of data at the raw spectral data, peptide, and protein levels. This platform enables uniform analysis and exchange of MS/MS data generated from a variety of different instruments, and assigned peptides using a variety of different database search programs. We demonstrate this by applying the pipeline to data sets generated by ThermoFinnigan LCQ, ABI 4700 MALDI-TOF/TOF, and Waters Q-TOF instruments, and searched in turn using SEQUEST, Mascot, and COMET.

Application of metabolome data in functional genomics: A conceptual strategy

A gene with yet unknown physiological function can be studied by changing its expression level followed by analysis of the resulting phenotype. This type of functional genomics study can be complicated by the occurrence of 'silent mutations', the phenotypes of which are not easily observable in terms of metabolic fluxes (e.g., the growth rate). Nevertheless, genetic alteration may give rise to significant yet complicated changes in the metabolome. We propose here a conceptual functional genomics strategy based on microbial metabolome data, which identifies changes in in vivo enzyme activities in the mutants. These predicted changes are used to formulate hypotheses to infer unknown gene functions. The required metabolome data can be obtained solely from high-throughput mass spectrometry analysis, which provides the following in vivo information: (1) the metabolite concentrations in the reference and the mutant strain; (2) the metabolic fluxes in both strains and (3) the enzyme kinetic parameters of the reference strain. We demonstrate in silico that changes in enzyme activities can be accurately predicted by this approach, even in 'silent mutants'.

Genomic analysis of anaerobic respiration in the archaeon Halobacterium sp. strain NRC-1: dimethyl sulfoxide and trimethylamine N-oxide as terminal electron acceptors

We have investigated anaerobic respiration of the archaeal model organism Halobacterium sp. strain NRC-1 by using phenotypic and genetic analysis, bioinformatics, and transcriptome analysis. NRC-1 was found to grow on either dimethyl sulfoxide (DMSO) or trimethylamine N-oxide (TMAO) as the sole terminal electron acceptor, with a doubling time of 1 day. An operon, dmsREABCD, encoding a putative regulatory protein, DmsR, a molybdopterin oxidoreductase of the DMSO reductase family (DmsEABC), and a molecular chaperone (DmsD) was identified by bioinformatics and confirmed as a transcriptional unit by reverse transcriptase PCR analysis. dmsR, dmsA, and dmsD in-frame deletion mutants were individually constructed. Phenotypic analysis demonstrated that dmsR, dmsA, and dmsD are required for anaerobic respiration on DMSO and TMAO. The requirement for dmsR, whose predicted product contains a DNA-binding domain similar to that of the Bat family of activators (COG3413), indicated that it functions as an activator. A cysteine-rich domain was found in the dmsR gene, which may be involved in oxygen sensing. Microarray analysis using a whole-genome 60-mer oligonucleotide array showed that the dms operon is induced during anaerobic respiration. Comparison of dmsR+ and DeltadmsR strains by use of microarrays showed that the induction of the dmsEABCD operon is dependent on a functional dmsR gene, consistent with its action as a transcriptional activator. Our results clearly establish the genes required for anaerobic respiration using DMSO and TMAO in an archaeon for the first time.

An efficient proteomics based strategy for the functional characterization of a novel halophilic enzyme fromHalobacterium salinarum

The extremely halophilic archaeon, Halobacterium salinarum grows in environments containing over 25% NaCl. The enzymes of this organism have thus been adapted to be active and stable in hypersaline conditions, which makes them strong candidates as robust industrial enzymes. In this study, the proteomics approach was applied to screen novel halophilic enzymes. We focused initially on proteins that are differentially expressed under different salt concentrations in culture media. After two-dimensional gel electrophoresis over a pH 3.5-4.5 range, 29 differentially expressed protein spots were identified by tandem mass spectrometry and six of these had no similarity to preexisting genes of known function. To predict the function of them, we used various bioinformatic methods. Among other proteins, we selected Vng0487h, which showed a high similarity to acetyltransferases. As a step toward assaying the enzymatic activity of this protein, we cloned the Vng0487h gene of H. salinarum and expressed and purified the recombinant protein with a glutathione-S-transferase (GST) tag in Escherichia coli. Using a GST-pulldown assay, a protein fragment derived from E. coli could interact with recombinant Vng0487h, and was identified to be the ribosomal protein L3. This protein showed high sequence homology with ribosomal protein L7/12 from E. coli and ribosomal protein L13p from H. salinarum. This suggests that Vng0487h acetylates a subunit of ribosomal protein, possibly L13p, in H. salinarum. During the present study, an efficient procedure was established to screen novel halophilic enzymes, and to predict and assess their functions.

Analysis of proteasome-dependent proteolysis in Haloferax volcanii cells, using short-lived green fluorescent proteins

Proteasomes are energy-dependent proteases that are central to the quality control and regulated turnover of proteins in eukaryotic cells. Dissection of this proteolytic pathway in archaea, however, has been hampered by the lack of substrates that are easily detected in whole cells. In the present study, we developed a convenient reporter system by functional expression of a green fluorescent protein variant with C-terminal fusions in the haloarchaeon Haloferax volcanii. The levels of this reporter protein correlated with whole-cell fluorescence that was readily detected in culture. Accumulation of the reporter protein was dependent on the sequence of the C-terminal amino acid fusion, as well as the presence of an irreversible, proteasome-specific inhibitor (clasto-lactacystin beta-lactone). This inhibitor was highly specific for H. volcanii 20S proteasomes, with a Ki of approximately 40 nM. In contrast, phenylmethanesulfonyl fluoride did not influence the levels of fluorescent reporter protein or inhibit 20S proteasomes. Together, these findings provide a powerful tool for the elucidation of protein substrate recognition motifs and the identification of new genes which may be involved in the proteasome pathway of archaea.

Increased quantitative proteome coverage with13C/12C-based, acid-cleavable isotope-coded affinity tag reagent and modified data acquisition scheme

Quantitative protein profiling using the isotope-coded affinity tag (ICAT) method and tandem mass spectrometry (MS) enables the pair-wise comparison of protein expression levels in biological samples. A new version of the ICAT reagent with an acid-cleavable bond, which allows removal of the biotin moiety prior to MS and which utilizes (13)C substitution for (12)C in the heavy-ICAT reagent rather than (2)H (for (1)H) as in the original reagent, was investigated. We developed and validated an MS data acquisition strategy using this new reagent that results in an increased number of protein identifications per experiment, without losing the accuracy of protein quantification. This was achieved by following a single survey (precursor) ion scan and serial collision induced dissociations (CIDs) of four different precursor ions observed in the prior survey scan. This strategy is common to many high-performance liquid chromatography-electrospray ionization (HPLC-ESI)-MS shotgun proteomic strategies, but heretofore not to ICAT experiments. This advance is possible because the new ICAT reagent uses (13)C as the "heavy" element rather than (2)H, thus, eliminating the slight delay in retention time of ICAT-labeled "light" peptides on a C18-based HPLC separation that occurs with (2)H and (1)H. Analyses using this new scheme of an ICAT-labeled trypsin-digested six protein mixture as well as a tryptic digest of a total yeast lysate, indicated that about two times more proteins were identified in a single analysis, and that there was no loss in accuracy of quantification.

Genome sequence of Haloarcula marismortui: a halophilic archaeon from the Dead Sea

We report the complete sequence of the 4,274,642-bp genome of Haloarcula marismortui, a halophilic archaeal isolate from the Dead Sea. The genome is organized into nine circular replicons of varying G+C compositions ranging from 54% to 62%. Comparison of the genome architectures of Halobacterium sp. NRC-1 and H. marismortui suggests a common ancestor for the two organisms and a genome of significantly reduced size in the former. Both of these halophilic archaea use the same strategy of high surface negative charge of folded proteins as means to circumvent the salting-out phenomenon in a hypersaline cytoplasm. A multitiered annotation approach, including primary sequence similarities, protein family signatures, structure prediction, and a protein function association network, has assigned putative functions for at least 58% of the 4242 predicted proteins, a far larger number than is usually achieved in most newly sequenced microorganisms. Among these assigned functions were genes encoding six opsins, 19 MCP and/or HAMP domain signal transducers, and an unusually large number of environmental response regulators-nearly five times as many as those encoded in Halobacterium sp. NRC-1--suggesting H. marismortui is significantly more physiologically capable of exploiting diverse environments. In comparing the physiologies of the two halophilic archaea, in addition to the expected extensive similarity, we discovered several differences in their metabolic strategies and physiological responses such as distinct pathways for arginine breakdown in each halophile. Finally, as expected from the larger genome, H. marismortui encodes many more functions and seems to have fewer nutritional requirements for survival than does Halobacterium sp. NRC-1.

Transcriptional regulation in Archaea

During the past few decades, it has become clear that microorganisms can thrive under the most diverse conditions, including extremes of temperature, pressure, salinity and pH. Most of these extremophilic organisms belong to the third domain of life, that of the Archaea. The organisms of this domain are of particular interest because most informational systems that are associated with archaeal genomes and their expression are reminiscent of those seen in Eucarya, whereas, most of their metabolic aspects are similar to those of Bacteria. A better understanding of the regulatory mechanisms of gene expression in Archaea will, therefore, help to integrate the body of knowledge regarding the regulatory mechanisms that underlie gene expression in all three domains of life.

Metabolite profiling: from diagnostics to systems biology

The concept of metabolite profiling has been around for several decades, but only recent technical innovations have allowed metabolite profiling to be carried out on a large scale - with respect to both the number of metabolites measured and the number of experiments carried out. As a result, the power of metabolite profiling as a technology platform for diagnostics, and the research areas of gene-function analysis and systems biology, is now beginning to be fully realized.

Systems biology, proteomics, and the future of health care: toward predictive, preventative, and personalized medicine

The emergence of systems biology is bringing forth a new set of challenges for advancing science and technology. Defining ways of studying biological systems on a global level, integrating large and disparate data types, and dealing with the infrastructural changes necessary to carry out systems biology, are just a few of the extraordinary tasks of this growing discipline. Despite these challenges, the impact of systems biology will be far-reaching, and significant progress has already been made. Moving forward, the issue of how to use systems biology to improve the health of individuals must be a priority. It is becoming increasingly apparent that the field of systems biology and one of its important disciplines, proteomics, will have a major role in creating a predictive, preventative, and personalized approach to medicine. In this review, we define systems biology, discuss the current capabilities of proteomics and highlight some of the necessary milestones for moving systems biology and proteomics into mainstream health care.

Insights into the relation between mrna and protein expression patterns: II. Experimental observations in Escherichia coli

There is a need for improved appreciation of the importance of genome-wide mRNA and protein expression measurements and their role in understanding translation and in relation to genome-wide mathematical frameworks for gene expression regulation. We investigated the use of a high-density microarray technique for mRNA expression analysis and a two-dimensional protein electrophoresis-tandem mass spectrometry method for protein analysis to monitor changes in gene expression. We applied these analytical tools in the context of an environmental perturbation of Escherichia coli cells-the addition of varying amounts of IPTG. We also tested the application of these tools to the study of a genetic perturbation of Escherichia coli cells-the ability of certain strains to hypersecrete the hemolysin protein. We observed a lack of correspondence between mRNA and protein expression profiles. Although our data do not include measurements on all expressed genes (because the ability to measure protein expression profiles is limiting), we observed that the qualitative and quantitative behavior of the measurements of a subset of expressed genes is similar to the behavior of the entire system. The change in observed average mRNA and protein amplification factors for 77 and 52 genes coincided with the observed change in mRNA amplification factor for the entire system. Furthermore, we found that the use of relative changes in expression could be used to elucidate mechanisms of gene expression regulation for the system studied, even when measurements were made on a small subset of the system.

Comprehensive de novo structure prediction in a systems-biology context for the archaea Halobacterium sp. NRC-1

BACKGROUND

Large fractions of all fully sequenced genomes code for proteins of unknown function. Annotating these proteins of unknown function remains a critical bottleneck for systems biology and is crucial to understanding the biological relevance of genome-wide changes in mRNA and protein expression, protein-protein and protein-DNA interactions. The work reported here demonstrates that de novo structure prediction is now a viable option for providing general function information for many proteins of unknown function.

RESULTS

We have used Rosetta de novo structure prediction to predict three-dimensional structures for 1,185 proteins and protein domains (<150 residues in length) found in Halobacterium NRC-1, a widely studied halophilic archaeon. Predicted structures were searched against the Protein Data Bank to identify fold similarities and extrapolate putative functions. They were analyzed in the context of a predicted association network composed of several sources of functional associations such as: predicted protein interactions, predicted operons, phylogenetic profile similarity and domain fusion. To illustrate this approach, we highlight three cases where our combined procedure has provided novel insights into our understanding of chemotaxis, possible prophage remnants in Halobacterium NRC-1 and archaeal transcriptional regulators.

CONCLUSIONS

Simultaneous analysis of the association network, coordinated mRNA level changes in microarray experiments and genome-wide structure prediction has allowed us to glean significant biological insights into the roles of several Halobacterium NRC-1 proteins of previously unknown function, and significantly reduce the number of proteins encoded in the genome of this haloarchaeon for which no annotation is available.

A statistical model for identifying proteins by tandem mass spectrometry

A statistical model is presented for computing probabilities that proteins are present in a sample on the basis of peptides assigned to tandem mass (MS/MS) spectra acquired from a proteolytic digest of the sample. Peptides that correspond to more than a single protein in the sequence database are apportioned among all corresponding proteins, and a minimal protein list sufficient to account for the observed peptide assignments is derived using the expectation-maximization algorithm. Using peptide assignments to spectra generated from a sample of 18 purified proteins, as well as complex H. influenzae and Halobacterium samples, the model is shown to produce probabilities that are accurate and have high power to discriminate correct from incorrect protein identifications. This method allows filtering of large-scale proteomics data sets with predictable sensitivity and false positive identification error rates. Fast, consistent, and transparent, it provides a standard for publishing large-scale protein identification data sets in the literature and for comparing the results obtained from different experiments.

Integrated genomic and proteomic analyses of gene expression in Mammalian cells

Using DNA microarrays together with quantitative proteomic techniques (ICAT reagents, two-dimensional DIGE, and MS), we evaluated the correlation of mRNA and protein levels in two hematopoietic cell lines representing distinct stages of myeloid differentiation, as well as in the livers of mice treated for different periods of time with three different peroxisome proliferative activated receptor agonists. We observe that the differential expression of mRNA (up or down) can capture at most 40% of the variation of protein expression. Although the overall pattern of protein expression is similar to that of mRNA expression, the incongruent expression between mRNAs and proteins emphasize the importance of posttranscriptional regulatory mechanisms in cellular development or perturbation that can be unveiled only through integrated analyses of both proteins and mRNAs.

Absolute Quantification of Specific Proteins in Complex Mixtures Using Visible Isotope-Coded Affinity Tags

The identification of proteins in complex mixtures is most useful when quantitative information is also obtained. We describe a new type of protein tagging reagent called the visible isotope-coded affinity tag (VICAT) which allows the absolute amount of a target protein or proteins to be quantified in a complex biological sample such as a eukaryotic cell lysate. VICAT reagents tag thiol groups of cysteines or thioacetylated amino groups and introduce into the tryptic peptide a biotin affinity handle, a visible moiety for tracking the chromatographic location of the target peptide by a method other than mass spectrometry, a photocleavable linker for removing a portion of the tag, and an isotope tag for distinguishing sample and internal standard peptides. We demonstrate the use of VICAT reagents together with isoelectric focusing of peptides on an immobilized gel strip followed by combined micro-liquid chromatography/electrospray ionization mass spectrometry operating in selected reaction monitoring mode to determine the absolute abundance of a specific protein, human group V phospholipase A(2), in eukaryotic cell lysates. It is found that human lung macrophages contain 66 fmol of this protein per 100 microg of cell protein. Western blot analysis of human group V phospholipase A(2) in macrophages gave inconclusive data. VICAT reagents should be useful for numerous applications including the analysis of candidate disease markers in complex mixtures such as serum.

Systems level insights into the stress response to UV radiation in the halophilic archaeon Halobacterium NRC-1

We report a remarkably high UV-radiation resistance in the extremely halophilic archaeon Halobacterium NRC-1 withstanding up to 110 J/m2 with no loss of viability. Gene knockout analysis in two putative photolyase-like genes (phr1 and phr2) implicated only phr2 in photoreactivation. The UV-response was further characterized by analyzing simultaneously, along with gene function and protein interactions inferred through comparative genomics approaches, mRNA changes for all 2400 genes during light and dark repair. In addition to photoreactivation, three other putative repair mechanisms were identified including d(CTAG) methylation-directed mismatch repair, four oxidative damage repair enzymes, and two proteases for eliminating damaged proteins. Moreover, a UV-induced down-regulation of many important metabolic functions was observed during light repair and seems to be a phenomenon shared by all three domains of life. The systems analysis has facilitated the assignment of putative functions to 26 of 33 key proteins in the UV response through sequence-based methods and/or similarities of their predicted three-dimensional structures to known structures in the PDB. Finally, the systems analysis has raised, through the integration of experimentally determined and computationally inferred data, many experimentally testable hypotheses that describe the metabolic and regulatory networks of Halobacterium NRC-1.

Chemical probes and tandem mass spectrometry: a strategy for the quantitative analysis of proteomes and subproteomes

Quantitative proteome profiling using mass spectrometry and stable isotope dilution is being widely applied for the functional analysis of biological systems and for the detection of clinical, diagnostic or prognostic marker proteins. Because of the enormous complexity of proteomes, their comprehensive analysis is unlikely to be routinely achieved in the near future. However, in recent years, significant progress has been achieved focusing quantitative proteomic analyses on specific protein classes or subproteomes that are rich in biologically or clinically important information. Such projects typically combine the use of chemical probes that are specific for a targeted group of proteins and may contain stable isotope signatures for accurate quantification with automated tandem mass spectrometry and bioinformatics tools for data analysis. In this review, we summarize technical and conceptual advances in quantitative subproteome profiling based on tandem mass spectrometry and chemical probes.

Diversity of bacteriorhodopsins in different hypersaline waters from a single Spanish saltern

Haloarchaeal rhodopsins are a diverse group of transmembrane proteins that use light energy to drive several different cellular processes. Two rhodopsins, bacteriorhodopsin and halorhodopsins, are H+ and Cl- ion pumps, respectively, and two rhodopsins, sensory rhodopsin I and II, regulate phototaxis. Bacteriorhodopsin is of special interest as it is a non-chlorophyll-based type of phototrophy (i.e. generation of chemical energy from light energy). However, very little is known about the diversity and distribution of rhodopsin genes in hypersaline environments. Here, we have used environmental PCR and cloning techniques to directly retrieve rhodopsin genes from three different salinity ponds located in a sea salt manufacturing facility near Alicante, Spain. Our survey resulted in the discovery of previously concealed variation including what is hypothesized to be bacteriorhodopsin genes from the uncultivated square morphotype that dominates these environments. In some instances, identical genes were discovered in seemingly different habitats suggesting that some haloarchaea are present over widely varying concentrations of salt.

Low-pass sequencing for microbial comparative genomics

BACKGROUND

We studied four extremely halophilic archaea by low-pass shotgun sequencing: (1) the metabolically versatile Haloarcula marismortui; (2) the non-pigmented Natrialba asiatica; (3) the psychrophile Halorubrum lacusprofundi and (4) the Dead Sea isolate Halobaculum gomorrense. Approximately one thousand single pass genomic sequences per genome were obtained. The data were analyzed by comparative genomic analyses using the completed Halobacterium sp. NRC-1 genome as a reference. Low-pass shotgun sequencing is a simple, inexpensive, and rapid approach that can readily be performed on any cultured microbe.

RESULTS

As expected, the four archaeal halophiles analyzed exhibit both bacterial and eukaryotic characteristics as well as uniquely archaeal traits. All five halophiles exhibit greater than sixty percent GC content and low isoelectric points (pI) for their predicted proteins. Multiple insertion sequence (IS) elements, often involved in genome rearrangements, were identified in H. lacusprofundi and H. marismortui. The core biological functions that govern cellular and genetic mechanisms of H. sp. NRC-1 appear to be conserved in these four other halophiles. Multiple TATA box binding protein (TBP) and transcription factor IIB (TFB) homologs were identified from most of the four shotgunned halophiles. The reconstructed molecular tree of all five halophiles shows a large divergence between these species, but with the closest relationship being between H. sp. NRC-1 and H. lacusprofundi.

CONCLUSION

Despite the diverse habitats of these species, all five halophiles share (1) high GC content and (2) low protein isoelectric points, which are characteristics associated with environmental exposure to UV radiation and hypersalinity, respectively. Identification of multiple IS elements in the genome of H. lacusprofundi and H. marismortui suggest that genome structure and dynamic genome reorganization might be similar to that previously observed in the IS-element rich genome of H. sp. NRC-1. Identification of multiple TBP and TFB homologs in these four halophiles are consistent with the hypothesis that different types of complex transcriptional regulation may occur through multiple TBP-TFB combinations in response to rapidly changing environmental conditions. Low-pass shotgun sequence analyses of genomes permit extensive and diverse analyses, and should be generally useful for comparative microbial genomics.

Cytoscape: a software environment for integrated models of biomolecular interaction networks

Cytoscape is an open source software project for integrating biomolecular interaction networks with high-throughput expression data and other molecular states into a unified conceptual framework. Although applicable to any system of molecular components and interactions, Cytoscape is most powerful when used in conjunction with large databases of protein-protein, protein-DNA, and genetic interactions that are increasingly available for humans and model organisms. Cytoscape's software Core provides basic functionality to layout and query the network; to visually integrate the network with expression profiles, phenotypes, and other molecular states; and to link the network to databases of functional annotations. The Core is extensible through a straightforward plug-in architecture, allowing rapid development of additional computational analyses and features. Several case studies of Cytoscape plug-ins are surveyed, including a search for interaction pathways correlating with changes in gene expression, a study of protein complexes involved in cellular recovery to DNA damage, inference of a combined physical/functional interaction network for Halobacterium, and an interface to detailed stochastic/kinetic gene regulatory models.

Complementary Analysis of the Mycobacterium tuberculosis Proteome by Two-dimensional Electrophoresis and Isotope-coded Affinity Tag Technology

Classical proteomics combined two-dimensional gel electrophoresis (2-DE) for the separation and quantification of proteins in a complex mixture with mass spectrometric identification of selected proteins. More recently, the combination of liquid chromatography (LC), stable isotope tagging, and tandem mass spectrometry (MS/MS) has emerged as an alternative quantitative proteomics technology. We have analyzed the proteome of Mycobacterium tuberculosis, a major human pathogen comprising about 4,000 genes, by (i) 2-DE and mass spectrometry (MS) and by (ii) the isotope-coded affinity tag (ICAT) reagent method and MS/MS. The data obtained by either technology were compared with respect to their selectivity for certain protein types and classes and with respect to the accuracy of quantification. Initial datasets of 60,000 peptide MS/MS spectra and 1,800 spots for the ICAT-LC/MS and 2-DE/MS methods, respectively, were reduced to 280 and 108 conclusively identified and quantified proteins, respectively. ICAT-LC/MS showed a clear bias for high M(r) proteins and was complemented by the 2-DE/MS method, which showed a preference for low M(r) proteins and also identified cysteine-free proteins that were transparent to the ICAT-LC/MS method. Relative quantification between two strains of the M. tuberculosis complex also revealed that the two technologies provide complementary quantitative information; whereas the ICAT-LC/MS method quantifies the sum of the protein species of one gene product, the 2-DE/MS method quantifies at the level of resolved protein species, including post-translationally modified and processed polypeptides. Our data indicate that different proteomic technologies applied to the same sample provide complementary types of information that contribute to a more complete understanding of the biological system studied.

Cytoscape: a software environment for integrated models of biomolecular interaction networks

Cytoscape is an open source software project for integrating biomolecular interaction networks with high-throughput expression data and other molecular states into a unified conceptual framework. Although applicable to any system of molecular components and interactions, Cytoscape is most powerful when used in conjunction with large databases of protein-protein, protein-DNA, and genetic interactions that are increasingly available for humans and model organisms. Cytoscape's software Core provides basic functionality to layout and query the network; to visually integrate the network with expression profiles, phenotypes, and other molecular states; and to link the network to databases of functional annotations. The Core is extensible through a straightforward plug-in architecture, allowing rapid development of additional computational analyses and features. Several case studies of Cytoscape plug-ins are surveyed, including a search for interaction pathways correlating with changes in gene expression, a study of protein complexes involved in cellular recovery to DNA damage, inference of a combined physical/functional interaction network for Halobacterium, and an interface to detailed stochastic/kinetic gene regulatory models.

Cytoscape: a software environment for integrated models of biomolecular interaction networks

Cytoscape is an open source software project for integrating biomolecular interaction networks with high-throughput expression data and other molecular states into a unified conceptual framework. Although applicable to any system of molecular components and interactions, Cytoscape is most powerful when used in conjunction with large databases of protein-protein, protein-DNA, and genetic interactions that are increasingly available for humans and model organisms. Cytoscape's software Core provides basic functionality to layout and query the network; to visually integrate the network with expression profiles, phenotypes, and other molecular states; and to link the network to databases of functional annotations. The Core is extensible through a straightforward plug-in architecture, allowing rapid development of additional computational analyses and features. Several case studies of Cytoscape plug-ins are surveyed, including a search for interaction pathways correlating with changes in gene expression, a study of protein complexes involved in cellular recovery to DNA damage, inference of a combined physical/functional interaction network for Halobacterium, and an interface to detailed stochastic/kinetic gene regulatory models.

Proteomic analysis of Pseudomonas aeruginosa grown under magnesium limitation

In this study, large-scale qualitative and quantitative proteomic technology was applied to the analysis of the opportunistic bacterial pathogen Pseudomonas aeruginosa grown under magnesium limitation, an environmental condition previously shown to induce expression of various virulence factors. For quantitative analysis, whole cell and membrane proteins were differentially labeled with isotope-coded affinity tag (ICAT) reagents and ICAT reagent-labeled peptides were separated by two-dimensional chromatography prior to analysis by electrospray ionization-tandem mass spectrometry (ESI-MS/MS) in an ion trap mass spectrometer (ITMS). To increase the number of protein identifications, gas-phase fractionation (GPF) in the m/z dimension was employed for analysis of ICAT peptides derived from whole cell extracts. The experiments confirmed expression of 1331 P. aeruginosa proteins of which 145 were differentially expressed upon limitation of magnesium. A number of conserved Gram-negative magnesium stress-response proteins involved in bacterial virulence were among the most abundant proteins induced in low magnesium. Comparative ICAT analysis of membrane versus whole cell protein indicated that growth of P. aeruginosa in low magnesium resulted in altered subcellular compartmentalization of large enzyme complexes such as ribosomes. This result was confirmed by 2-D PAGE analysis of P. aeruginosa outer membrane proteins. This study shows that large-scale quantitative proteomic technology can be successfully applied to the analysis of whole bacteria and to the discovery of functionally relevant biologic phenotypes.

Quantitative proteomic analysis of chromatin-associated factors

A method to identify and quantify chromatin-associated proteins has been developed and applied to the analysis of changes in chromatin-associated proteins induced by Myc oncoprotein expression in human B lymphocytes. Chromatin-enriched fractions were isolated by differential detergent/salt extraction and analyzed by ICAT reagent labeling, multi-dimensional chromatography and tandem mass spectrometry. Many known chromatin-associated regulatory factors were identified and quantified. The method will be widely applicable to various biological systems and reveal changes in chromatin-associated regulatory factors that underlie biological phenomena.

Whole-genome DNA microarray analysis of a hyperthermophile and an archaeon: Pyrococcus furiosus grown on carbohydrates or peptides

The first complete-genome DNA microarray was constructed for a hyperthermophile or a nonhalophilic archaeon by using the 2,065 open reading frames (ORFs) that have been annotated in the genome of Pyrococcus furiosus (optimal growth temperature, 100 degrees C). This was used to determine relative transcript levels in cells grown at 95 degrees C with either peptides or a carbohydrate (maltose) used as the primary carbon source. Approximately 20% (398 of 2065) of the ORFs did not appear to be significantly expressed under either growth condition. Of the remaining 1,667 ORFs, the expression of 125 of them (8%) differed by more than fivefold between the two cultures, and 82 of the 125 (65%) appear to be part of operons, indicating extensive coordinate regulation. Of the 27 operons that are regulated, 5 of them encode (conserved) hypothetical proteins. A total of 18 operons are up-regulated (greater than fivefold) in maltose-grown cells, including those responsible for maltose transport and for the biosynthesis of 12 amino acids, of ornithine, and of citric acid cycle intermediate products. A total of nine operons are up-regulated (greater than fivefold) in peptide-grown cells, including those encoding enzymes involved in the production of acyl and aryl acids and 2-ketoacids, which are used for energy conservation. Analyses of the spent growth media confirmed the production of branched-chain and aromatic acids during growth on peptides. In addition, six nonlinked enzymes in the pathways of sugar metabolism were regulated more than fivefold--three in maltose-grown cells that are unique to the unusual glycolytic pathway and three in peptide-grown cells that are unique to gluconeogenesis. The catalytic activities of 16 metabolic enzymes whose expression appeared to be highly regulated in the two cell types correlated very well with the microarray data. The degree of coordinate regulation revealed by the microarray data was unanticipated and shows that P. furiosus can readily adapt to a change in its primary carbon source.

Mass spectrometry-based proteomics

Recent successes illustrate the role of mass spectrometry-based proteomics as an indispensable tool for molecular and cellular biology and for the emerging field of systems biology. These include the study of protein-protein interactions via affinity-based isolations on a small and proteome-wide scale, the mapping of numerous organelles, the concurrent description of the malaria parasite genome and proteome, and the generation of quantitative protein profiles from diverse species. The ability of mass spectrometry to identify and, increasingly, to precisely quantify thousands of proteins from complex samples can be expected to impact broadly on biology and medicine.

Advances in quantitative proteomics via stable isotope tagging and mass spectrometry

The high-throughput identification and accurate quantification of proteins are essential components of proteomic strategies for studying cellular functions and processes. Techniques that are largely based on stable isotope protein or peptide labeling and automated tandem mass spectrometry are increasingly being applied in quantitative proteomic studies. Over the past year, significant progress has been made toward improving and diversifying these technologies with respect to the methods for stable isotope labeling, process automation and data processing and analysis. Advances in stable isotope protein labeling and recent biological studies that used stable isotope based quantitative proteomics techniques are reviewed.

The digital code of DNA

The discovery of the structure of DNA transformed biology profoundly, catalysing the sequencing of the human genome and engendering a new view of biology as an information science. Two features of DNA structure account for much of its remarkable impact on science: its digital nature and its complementarity, whereby one strand of the helix binds perfectly with its partner. DNA has two types of digital information--the genes that encode proteins, which are the molecular machines of life, and the gene regulatory networks that specify the behaviour of the genes.