The files in the subdirectories of
extdata support the examples in the package documentation and vignettes.
abundance contain protein abundance data:
stressis a data frame listing proteins identified in selected proteomic stress response experiments. The names of proteins begin at row 3, and columns are all the same length (padded as necessary at the bottom by
NAs). Names correspond to ordered locus names (for Sce) or gene names (for Eco). The column names identify the experiments, the first row contains the name of the organism (Sce or Eco) and the third row has the reference key for the source of the data (listed in
AA03.csvhas reference abundances for 71 proteins taken from Fig. 3 of Anderson and Anderson, 2002 (as corrected in Anderson and Anderson, 2003). The columns with data taken from these sources are type (hemoglobin, plasma, tissue, or interleukin), description (name used in the original figure), log10(pg/ml) (upper limit of abundance interval shown in Anderson and Anderson, 2003, log10 of concentration in pg/ml). The additional columns are data derived from a search of the SWISS-PROT/UniProtKB database based on the descriptions of the proteins: name (nominal UniProtKB name for this protein), name2 (other UniProtKB names(s) that could apply to the protein), and note (notes based on searching for a protein of this description). The amino acid compositions of all proteins whose names are not NA are included in
abbrvcolumn for the proteins contains the description given by Anderson and Anderson, 2003, followed by (in parentheses) the UniProtKB accession number. Annotated initiator methionines (e.g. for ferritin, myoglobin, ENOG), signal peptides or propeptides were removed from the proteins (except where they are not annotated in UniProtKB: IGHG1, IGHA1, IGHD, MBP). In cases were multiple isoforms are present in UniProtKB (e.g. Albumin) only the first isoform was taken. In the case of C4 Complement (CO4A) and C5 Complement (CO5), the amino acid composition of only the alpha chains are listed. In the case of the protein described as iC3b, the amino acid sequence is taken to be that of Complement C3c alpha' chain fragment 1 from CO3, and is given the name CO3.C3c. The non-membrane (soluble) chains of TNF-binding protein (TNR1A) and TNF-alpha (TNFA) were used. Rantes, MIP-1 beta and MIP-1 alpha were taken from C-C motif chemokines (CCL5, CCL4, CCL3 respectively). C-peptide was taken from the corresponding annotation for insulin and here is named INS.C. See
read.exprfor examples that use this file.
ISR+08.csvhas columns excerpted from Additional File 2 of Ishihama et al. (2008) for protein abundances in E. coli cytosol. The columns in this file are ID (Swiss-Prot ID), accession (Swiss-Prot accession), emPAI (exponentially modified protein abundance index), copynumber (emPAI-derived copy number/cell), GRAVY (Kyte-Doolittel), FunCat (FunCat class description), PSORT (PSORT localisation), ribosomal (yes/no). See
read.exprfor examples that use this file.
yeastgfp.csv.xzHas 28 columns; the names of the first five are
GFP visualized?, and
abundance. The remaining columns correspond to the 23 subcellular localizations considered in the YeastGFP project (Huh et al., 2003 and Ghaemmaghami et al., 2003) and hold values of either
Ffor each protein. yeastgfp.csv was downloaded on 2007-02-01 from http://yeastgfp.ucsf.edu using the Advanced Search, setting options to download the entire dataset and to include localization table and abundance, sorted by orf number. See
yeastgfpfor examples that use this file.
bison contain BLAST results and taxonomic information for a metagenome:
bisonP...are partial tabular BLAST results for proteins in the Bison Pool Environmental Genome. Protein sequences predicted in the metagenome were downloaded from the Joint Genome Institute's IMG/M system on 2009-05-13. The target database for the searches was constructed from microbial protein sequences in National Center for Biotechnology Information (NCBI) RefSeq database version 57, representing 7415 microbial genomes. The ‘blastall’ command was used with the default setting for E value cuttoff (10.0) and options to make a tabular output file consisting of the top 20 hits for each query sequence. The function
read.blastwas used to extract only those hits with E values less than or equal to 1e-5 and with sequence similarity (percent identity) at least 30 percent, and to keep only the first hit for each query sequence. The function
write.blastwas used to save partial BLAST files (only selected columns). The files provided with CHNOSZ contain the first 5,000 hits for each sampling site at Bison Pool, representing between about 7 to 15 percent of the first BLAST hits after similarity and E value filtering.
gi.taxid.txt.xzis a table that lists the sequence identifiers (gi numbers) that appear in the example BLAST files (see above), together with the corresponding taxon ids used in the NCBI databases. This file is not a subset of the complete ‘gi_taxid_prot.dmp.gz’ available at ftp://ftp.ncbi.nih.gov/pub/taxonomy/ but instead is a subset of ‘gi.taxid.txt’ generated from the RefSeq release catalog using ‘gencat.sh’ in the
id.blastfor an example that uses this file and the BLAST files described above.
cpetc contain heat capacity data and other thermodynamic properties:
PM90.csvHeat capacities of four unfolded aqueous proteins taken from Privalov and Makhatadze, 1990. Names of proteins are in the first column, temperature in °C in the second, and heat capacities in J mol^-1 K^-1 in the third. See
ionize.aafor an example that uses this file.
RH95.csvHeat capacity data for iron taken from Robie and Hemingway, 1995. Temperature in Kelvin is in the first column, heat capacity in J K^-1 mol^-1 in the second. See
subcrtfor an example that uses this file.
RT71.csvpH titration measurements for unfolded lysozyme (LYSC_CHICK) taken from Roxby and Tanford, 1971. pH is in the first column, net charge in the second. See
ionize.aafor an example that uses this file.
SOJSH.csvExperimental equilibrium constants for the reaction NaCl(aq) = Na+ + Cl- as a function of temperature and pressure taken from Fig. 1 of Shock et al., 1992. Data were extracted from the figure using g3data (http://www.frantz.fi/software/g3data.php). See
waterfor an example that uses this file.
V.CH4.HWM96.csvApparent molar heat capacities and volumes of CH4 in dilute aqueous solutions reported by Hnedkovsky and Wood, 1997 and Hnedkovsky et al., 1996. See
EOSregressfor examples that use these files.
BKM60_Fig7.datEh-pH values for normal, wet and waterlogged soils from Fig. 7 of Baas Becking et al., 1960. See the ‘anintro’ vignette for an example that uses this file.
SC10_Rainbow.csvValues of temperature (°C), pH and logarithms of activity of CO2, H2, NH4+, H2S and CH4 for mixing of seawater and hydrothermal fluid at Rainbow field (Mid-Atlantic Ridge), taken from Shock and Canovas, 2010.
fasta contain protein sequences:
HTCC1062.faa.xzis a FASTA file of 1354 protein sequences in the organism Pelagibacter ubique HTCC1062 downloaded from the NCBI RefSeq collection on 2009-04-12. The search term was Protein: txid335992[Organism:noexp] AND "refseq"[Filter]. See
revisitfor examples that use this file.
EF-Tu.alnconsists of aligned sequences (394 amino acids) of elongation factor Tu (EF-Tu). The sequences correspond to those taken from UniProtKB for ECOLI (Escherichia coli), THETH (Thermus thermophilus) and THEMA (Thermotoga maritima), and reconstructed ancestral sequences taken from Gaucher et al., 2003 (maximum likelihood bacterial stem and mesophilic bacterial stem, and alternative bacterial stem). See the ‘formation’ vignette for an example that uses this file.
protein contain protein composition data for model organisms.
read.expr for examples that use these files.
Sce.csv.xzData frame of amino acid composition of 6716 proteins from the Saccharomyces Genome Database (SGD). Values in the first three columns are the
ORFnames of proteins,
GENEnames. The remaining twenty columns (
VAL) contain the numbers of the respective amino acids in each protein. The sources of data for Sce.csv are the files protein_properties.tab and SGD_features.tab (for the gene names), downloaded from http://www.yeastgenome.org on 2013-08-24.
Eco.csv.xzAmino acid compositions of 4407 proteins in Escherichia coli strain K12. Format is the one used
thermo$protein, with columns protein holding the gene name, organism set to ECOLI, and
abbrvholding the UniProt ID. The source of data is the file ECOLI.fas downloaded from the HAMAP (High-quality Automated and Manual Annotation of microbial Proteomes system) FTP site (Gattiker et al., 2003) on 2010-09-25 (old URL:
refseq contain code and results of processing NCBI Reference Sequences (RefSeq) for microbial proteins, using RefSeq release 61 of 2013-09-09:
README.txtInstructions for producing the data files.
gencat.shBash script to extract microbial protein records from the RefSeq catalog.
gi.taxid.txtOutput from above. The complete file is too large to distribute with CHNOSZ, but a portion is included in
extdata/bisonto support processing example BLAST files for the Bison Pool metagenome (based on RefSeq 57, 2013-01-08).
mkfaa.shCombine the contents of .faa.gz files into a single FASTA file (to use e.g. for making a BLAST database).
protein.refseq.RCalculate average amino acid composition of all proteins for each organism identified by a taxonomic ID.
trim_refseq.RKeep only selected organism names (reduces number of taxa from 6758 to 779, helps to control package size).
protein_refseq.csv.xzOutput from above. See example in
taxid.names.RGenerate a table of scientific names for the provided taxids. Requires the complete
nodes.dmpfrom NCBI taxonomy files.
taxid_names.csv.xzOutput from above. NOTE: For backward compatibility with the example BLAST files for the Bison Pool metagenome, the packaged file merges records for taxids found in either RefSeq 57 or 61. Certain taxids in release 57 were not located in the current RefSeq catalog, probably related to the transition to the “WP” multispecies accessions (ftp://ftp.ncbi.nlm.nih.gov/refseq/release/announcements/WP-proteins-06.10.2013.pdf). See example for
taxonomy contain example taxonomic data files:
nodes.dmpare excerpts of the taxonomy files available on the NCBI ftp site (ftp://ftp.ncbi.nih.gov/pub/taxonomy/taxdump.tar.gz, accessed 2010-02-15). These example files contain only the entries for Escherichia coli K-12, Saccharomyces cerevisiae, Homo sapiens, Pyrococcus furisosus and Methanocaldococcus jannaschii (taxids 83333, 4932, 9606, 186497, 243232) and the higher-ranking nodes (genus, family, etc.) in the respective lineages. See
taxonomyfor examples that use this file.
thermo contain additional thermodynamic data and group additivity definitions:
OBIGT-2.csvcontains supplementary thermodynamic data in the same format as the primary database in
data/OBIGT.csv. Data for some entries in the primary database are taken from different literature sources in this file. The default action of
add.obigtis to add the contents of this file to CHNOSZ's working database in
diagramand the code of
anim.TCAfor examples that use this file.
obigt_check.csvcontains the results of running
check.obigtto check the internal consistency of entries in the primary and supplementary databases.
groups_big.csvGroup contribution matrix: five structural groups on the columns ([-CH3],[-CH2-],[-CH2OH],[-CO-],[-COOH]) and 24 compounds on the rows (alkanes, alcohols, ketones, acids, multiply substituted compounds).
groups_small.csvGroup contribution matrix: twelve bond-specific groups on the columns, and 25 compounds on the rows (as above, plus isocitrate). Group identity and naming conventions adapted from Benson and Buss (1958) and Domalski and Hearing (1993). See the ‘xadditivity’ vignette for examples that use this file and
RH98_Table15.csvGroup stoichiometries for high molecular weight crystalline and liquid organic compounds taken from Table 15 of Richard and Helgeson, 1998. The first three columns have the
state(cr or liq). The remaining columns have the numbers of each group in the compound; the names of the groups (columns) correspond to species in
thermo$obigt. The compound named 5a(H),14a(H)-cholestane in the paper has been changed to 5a(H),14b(H)-cholestane here to match the group stoichiometry given in the table. See
RH2obigtfor a function that uses this file.
DLEN67.csvStandard Gibbs energies of formation, in kcal/mol, from Dayhoff et al., 1967, for nitrogen (N2) plus 17 compounds shown in Fig. 2 of Dayhoff et al., 1964, at 300, 500, 700 and 1000 K.
Anderson, N. L. and Anderson, N. G. (2002) The human plasma proteome: History, character and diagnostic prospects. Molecular and Cellular Proteomics 1, 845–867. http://dx.doi.org/10.1074/mcp.R200007-MCP200
Anderson, N. L. and Anderson, N. G. (2003) The human plasma proteome: History, character and diagnostic prospects (Vol. 1 (2002) 845-867). Molecular and Cellular Proteomics 2, 50. http://dx.doi.org/10.1074/mcp.A300001-MCP200
Baas Becking, L. G. M., Kaplan, I. R. and Moore, D. (1960) Limits of the natural environment in terms of pH and oxidation-reduction potentials. Journal of Geology 68(3), 243–284. http://www.jstor.org/stable/30059218
Benson, S. W. and Buss, J. H. (1958) Additivity rules for the estimation of molecular properties. Thermodynamic properties. J. Chem. Phys. 29, 546–572. http://dx.doi.org/10.1063/1.1744539
Dayhoff, M. O. and Lippincott, E. R. and Eck, R. V. (1964) Thermodynamic Equilibria In Prebiological Atmospheres. Science 146, 1461–1464. http://dx.doi.org/10.1126/science.146.3650.1461
Dayhoff, M. O. and Lippincott, E. R., Eck, R. V. and Nagarajan (1967) Thermodynamic Equilibrium In Prebiological Atmospheres of C, H, O, N, P, S, and Cl. Report SP-3040, National Aeronautics and Space Administration. http://ntrs.nasa.gov/search.jsp?R=19670017966
Domalski, E. S. and Hearing, E. D. (1993) Estimation of the thermodynamic properties of C-H-N-O-S-Halogen compounds at 298.15 K J. Phys. Chem. Ref. Data 22, 805–1159. http://dx.doi.org/10.1063/1.555927
Gattiker, A., Michoud, K., Rivoire, C., Auchincloss, A. H., Coudert, E., Lima, T., Kersey, P., Pagni, M., Sigrist, C. J. A., Lachaize, C., Veuthey, A.-L., Gasteiger, E. and Bairoch, A. (2003) Automatic annotation of microbial proteomes in Swiss-Prot. Comput. Biol. Chem. 27, 49–58. http://dx.doi.org/10.1016/S1476-9271(02)00094-4
Gaucher, E. A., Thomson, J. M., Burgan, M. F. and Benner, S. A (2003) Inferring the palaeoenvironment of ancient bacteria on the basis of resurrected proteins. Nature 425(6955), 285–288.
Ghaemmaghami, S., Huh, W., Bower, K., Howson, R. W., Belle, A., Dephoure, N., O'Shea, E. K. and Weissman, J. S. (2003) Global analysis of protein expression in yeast. Nature 425(6959), 737–741.
Huh, W. K., Falvo, J. V., Gerke, L. C., Carroll, A. S., Howson, R. W., Weissman, J. S. and O'Shea, E. K. (2003) Global analysis of protein localization in budding yeast. Nature 425(6959), 686–691.
HAMAP system. HAMAP FTP directory, ftp://ftp.expasy.org/databases/hamap/
Hnedkovsky, L., Wood, R. H. and Majer, V. (1996) Volumes of aqueous solutions of CH4, CO2, H2S, and NH3 at temperatures from 298.15 K to 705 K and pressures to 35 MPa. J. Chem. Thermodyn. 28, 125–142. http://dx.doi.org/10.1006/jcht.1996.0011
Hnedkovsky, L. and Wood, R. H. (1997) Apparent molar heat capacities of aqueous solutions of CH4, CO2, H2S, and NH3 at temperatures from 304 K to 704 K at a pressure of 28 MPa. J. Chem. Thermodyn. 29, 731–747. http://dx.doi.org/10.1006/jcht.1997.0192
Ishihama, Y., Schmidt, T., Rappsilber, J., Mann, M., Hartl, F. U., Kerner, M. J. and Frishman, D. (2008) Protein abundance profiling of the Escherichia coli cytosol. BMC Genomics 9:102. http://dx.doi.org/10.1186/1471-2164-9-102
Joint Genome Institute (2007) Bison Pool Environmental Genome. Protein sequence files downloaded from IMG/M (http://img.jgi.doe.gov/cgi-bin/m/main.cgi?section=FindGenomes&page=findGenomes)
Privalov, P. L. and Makhatadze, G. I. (1990) Heat capacity of proteins. II. Partial molar heat capacity of the unfolded polypeptide chain of proteins: Protein unfolding effects. J. Mol. Biol. 213, 385–391. http://dx.doi.org/10.1016/S0022-2836(05)80198-6
Richard, L. and Helgeson, H. C. (1998) Calculation of the thermodynamic properties at elevated temperatures and pressures of saturated and aromatic high molecular weight solid and liquid hydrocarbons in kerogen, bitumen, petroleum, and other organic matter of biogeochemical interest. Geochim. Cosmochim. Acta 62, 3591–3636. http://dx.doi.org/10.1016/S0016-7037(97)00345-1
Robie, R. A. and Hemingway, B. S. (1995) Thermodynamic Properties of Minerals and Related Substances at 298.15 K and 1 Bar (10^5 Pascals) Pressure and at Higher Temperatures. U. S. Geol. Surv., Bull. 2131, 461 p. http://www.worldcat.org/oclc/32590140
Roxby, R. and Tanford, C. (1971) Hydrogen ion titration curve of lysozyme in 6 M guanidine hydrochloride. Biochemistry 10, 3348–3352. http://dx.doi.org/10.1021/bi00794a005
SGD project. Saccharomyces Genome Database, http://www.yeastgenome.org
Shock, E. L., Oelkers, E. H., Johnson, J. W., Sverjensky, D. A. and Helgeson, H. C. (1992) Calculation of the thermodynamic properties of aqueous species at high pressures and temperatures: Effective electrostatic radii, dissociation constants and standard partial molal properties to 1000 °C and 5 kbar. J. Chem. Soc. Faraday Trans. 88, 803–826. http://dx.doi.org/10.1039/FT9928800803
Shock, E. and Canovas, P. (2010) The potential for abiotic organic synthesis and biosynthesis at seafloor hydrothermal systems. Geofluids 10, 161–192. http://dx.doi.org/10.1111/j.1468-8123.2010.00277.x
YeastGFP project. Yeast GFP Fusion Localization Database, http://yeastgfp.ucsf.edu; Current location: http://yeastgfp.yeastgenome.org