Skip to main content
  • Research article
  • Open access
  • Published:

Comparative proteome analysis of psychrophilic versus mesophilic bacterial species: Insights into the molecular basis of cold adaptation of proteins



Cold adapted or psychrophilic organisms grow at low temperatures, where most of other organisms cannot grow. This adaptation requires a vast array of sequence, structural and physiological adjustments. To understand the molecular basis of cold adaptation of proteins, we analyzed proteomes of psychrophilic and mesophilic bacterial species and compared the differences in amino acid composition and substitution patterns to investigate their likely association with growth temperatures.


In psychrophilic bacteria, serine, aspartic acid, threonine and alanine are overrepresented in the coil regions of secondary structures, whilst glutamic acid and leucine are underrepresented in the helical regions. Compared to mesophiles, psychrophiles comprise a significantly higher proportion of amino acids that contribute to higher protein flexibility in the coil regions of proteins, such as those with tiny/small or neutral side chains. Amino acids with aliphatic, basic, aromatic and hydrophilic side chains are underrepresented in the helical regions of proteins of psychrophiles. The patterns of amino acid substitutions between the orthologous proteins of psychrophiles versus mesophiles are significantly different for several amino acids when compared to their substitutions in orthologous proteins of within the mesophiles or psychrophiles.


Current results provide quantitative substitution preferences (or avoidance) of amino acids that lead to the adaptation of proteins to cold temperatures. These finding would help future efforts in selecting mutations for rational design of proteins with enhanced psychrophilic properties.


Microorganisms that live under forbidding conditions are called extremophiles, whose discovery points out the unique adaptability of primitive life-forms. These microorganisms are grouped according to their optimal growth conditions in which they exist such as acidophiles (exhibiting optimum growth in acidic pH conditions), alkaliphiles (thriving in alkaline pH conditions), barophiles (surviving under great pressures), endoliths (living in deep inside rocks), halophiles (thriving in high salt concentrations), psychrophiles (optimal temperature below 20°C), and the thermophiles (optimal temperature between 45–80°C), hyperthermophiles (optimal temperature above 80°C) [1]. The largest coverage of known extremophile conditions of the earth's biosphere is below 10°C. For example, three fourths of earth is covered by oceans, which maintain an average temperature of one to three degrees centigrade. Furthermore, the vast land areas of the Arctic and Antarctic are permanently frozen throughout the year [1]. Other few examples of cryo habitats include cold deserts, high alpine soils, sea ice, cold caves, marine sediments, permafrost soils, glacier, snow etc.

The majority of known psychrophiles belong to varieties of archaea and bacteria, and a few species of yeast, fungi and algae [2]. The ability to thrive at life-endangering effects of low temperatures, close to freezing point of water, requires a vast array of adaptations from all their cellular components, including their membranes, energy-generating systems, protein synthesis machinery, biodegradative enzymes and the components responsible for nutrient uptake etc., to maintain metabolism, sustain growth and reproduction compatible with life in these low temperature conditions [3, 4]. Having evolved with special mechanisms, the psychrophiles successfully colonized these niches [2, 5]. Psychrophilic proteins display sequences and structures comparable with those of their meso and (hyper) thermophilic homolog's, especially enzymes with their ability to work efficiently as catalysts at low temperatures [6]. The thermolability of these proteins at moderate temperatures warrant tremendous industrial applications in biotechnology, bioremediation, food, textiles, detergents bio-catalysis under low-water conditions and detergents etc [59].

Due to above facts, historically starting from mid-1970's, much attention was paid mainly to sequence and structural attributes contributing to adaptation of proteins (mainly enzymes) to high temperature conditions. Many investigators have compared sequence and structure-based parameters among thermophilic and mesophilic proteins [10]. With the advent of pioneering efforts in late 1990's in solving three dimensional structures of cryophilic enzymes such as alpha-amylase [11]; alkaline protease [12]; triose phosphate isomerase [13]; malate dehydrogenase [14] from Antarctic microorganisms, and due to handful of available structures in the protein data bank (PDB), groups have focused to address the structural basis of proteins in cold adaptation [4, 1520].

The steady increase in sequencing of proteomes of extremophiles has opened many new avenues in understanding adaptations to extreme conditions [16, 2125]. A comprehensive comparison of global amino acid preferences and substitution patterns as deduced from proteomes of different organisms is now possible [2628]. Using homologous sequences, clustering along with various statistical methods; we conducted an extensive analysis of proteomes of psychrophilic, mesophilic, thermophilic and hyperthermophilic microorganisms to examine a possible correlation of amino acid substitution patterns with adaptation to their respective optimal growth conditions. In this manuscript we discuss the results from comparative analysis of fully sequenced proteomes of six members from each of psychrophilic and mesophilic organisms.


On average we analyzed 2,816 proteins with 875,219 amino acids per proteome of mesophiles and 3665 proteins with 1,169,678 amino acids per proteome of psychrophiles. The amino acid (AA) frequencies given in Table 1 show that some of the AA differed significantly in psychrophile proteomes when compared to mesophile proteomes. When compared to psychrophiles, the mesophile proteomes show larger standard deviation for residues indicating that the six proteomes of mesophiles we used are considerably more divergent than the proteomes of psychrophiles. The frequencies of individual amino acids as well as property groups were further analyzed with student t-test.

Table 1 The composition of individual amino acids and property groups in protein sequences of psychrophilic and mesophilic proteomes.

Amino acid composition preferences

The t-test results demonstrate significant preferences in frequencies of amino acid occurrences and property groups in psychrophilic proteomes as compared to mesophilic proteomes or vice versa (Table 1). The compositional trend of AA is somewhat similar in both types of genomes. However, as indicated by t-values from Table 1, there are a few AA residues such as A, D, S and T, significantly preferred in psychrophiles as compared to mesophiles. On the other hand, AA residues E and L are significantly less favored in psychrophile proteomes. When comparing frequencies of occurrences of property groups of AAs, we observe that tiny/small and neutral amino acid groups are significantly preferred in psychrophiles where as charged, basic, aromatic and hydrophilic groups are significantly less favored as shown by their corresponding t-values in Table 1. When we compared the AA compositions of the sequences in alignments of respective orthologous proteins alone (the data not given), we observed similar trends.

Secondary Structural Elements

The composition of AA of psychrophilic and mesophilic proteomes in three major secondary structural elements, α-helices, β-sheets and coils, is given in Table 2. Collectively taken, the psychrophilic proteomes contain significantly less number of residues (~2%) in the α-helices and significantly more number of residues (~2%) in the coil regions. The majority of amino acids exhibit similar compositions in either of the two genome sequences. However, amino acids E, F L, N and Y show significantly low frequencies in α-helices of psychrophilic proteomes and amino acids A, D, G, S, T, and V are significantly high in the coil region of psychrophilic proteomes. The amino acid, E is significantly low in the coil region of psychrophilic proteomes. Except in an increase in Alanine residues, β-sheets of psychrophile proteomes did not show any significant changes as compared to mesophiles.

Table 2 Distribution of amino acids and property group parameters in the predicted secondary structural elements of psychrophilic and mesophilic proteomes.

When evaluated for the frequencies of occurrence of property groups of amino acids, the majority of them show significantly low frequencies in helices of psychrophilic proteomes. Except in tiny group of AAs, the β-sheet regions of psychrophile proteomes did not show any significant changes (Table 2). The tiny, small, hydrophobic, neutral, acidic, aliphatic, and non-polar amino acid groups showed significantly high frequencies in the coil/loop regions of psychrophilic proteomes.

Comparative Proteome Analysis

Towards identification of residue substitutions, likely to have undergone in psychrophilic proteins as the species adapted to cold temperatures, a comparative proteome analysis was performed on the basis of amino acid substitutions occurred between the orthologous protein sequences of psychrophile and mesophile proteomes. The orthologous sequence pairs for a protein sequence with significant length coverage only were considered for the analysis. Coverage of hits with respect to each proteome is shown in Table 3. On average 16.3% psychrophile proteins have orthologous proteins in mesophiles we used. This was cross-checked using mesophile proteomes as query sequences and the psychrophile proteomes as subject sequences, where we found 13.9% orthologs (Table 3). We used sequence alignments of these sequence pairs to compute substitutions of amino acids between mesophilic and psychrophilic proteomes and vice versa and the obtained values were averaged. On average, ten to twenty five percent of sequences from individual proteomes exhibited best hit homologues from members of other thermal groups (Table 3). This percent depended on the size of the proteome under consideration. The higher the number of proteins in a query proteome the higher percentage of hits from the subject proteomes searched. This may be because some of the paralogous sequences selecting the same protein as its ortholog (this redundancy was removed in our final data). We also considered homologous proteins among the psychrophiles and mesophiles, and calculated the substitutions within them to use as background substitution frequencies [see Eqn. (ii) and (iii)]. We observed, on average 16.7% and 17.1% orthologous proteins within the psychrophilic and mesophilic proteomes, respectively (Table 3).

Table 3 Percentage of orthologs (best hits) in Psychrophiles and Mesophiles with respect to total number of sequences (size) in each of organisms considered.

Amino Acid Substitution Patterns

Log odd scores (LOS) of AA substitutions were calculated using the frequency of occurrence of substitutions among orthologous proteins of psychrophiles and mesophiles by normalizing with frequency of occurrence of substitutions within the proteomes of the same temperature sensitive group. In LOS calculation the substitutions influenced by factors other than temperature are nullified and values represent true substitution due to cold adaptation of species. In Table 4 we show LOSMeso values computed using equation (ii) as described in methods. These values clearly show that psychrophilic proteins avoid containing the amino acids E, F, K, N and Y. On the other hand they prefer containing the residues A, D, G, S, and T as compared to mesophile proteins. The individual values in Table 4 show that to what extant certain substitutions are favored or avoided depending on corresponding LOSMeso score +ve or -ve, respectively. For example the W in mesophiles mutating to S in the psychrophiles is highly favored with LOS scores of 10.9. On the other hand G in mesophiles mutating to K in psychrophiles is avoided with LOS score of -12.4. The LOSPsychro [Eqn. (iii)] scores calculated using substitution frequencies among psychrophiles as normalizing factor have shown similar results (data not shown).

Table 4 Log odd scores (LOS) of amino acid substitutions calculated using the equation (ii).

Student t-test was further applied to evaluate level of significance of LOS substitution scores and these data were shown in Additional files (Additional file 1 &2). Table Additional file 1(a) shows t-values for the LOSMeso scores and Additional file 1(b) shows t-values for the LOSPsychro scores. The LOSMeso scores (Table 4) with their corresponding t-values (Additional file 1) greater than 1.37 or less than -1.37 could be considered as significantly preferred or avoided, respectively, at 90% confidence level (shaded in color, Table 4). It can be seen from the table that there are about 45% substitutions that are shown to be significant. Further, the LOSMeso for AA property group substitutions are given in Table 5 and their corresponding t-values are given in Additional file 2. It is clearly seen from Table 5 that there is high preference for tiny, small and neutral AAs whereas charged (including both basic and acidic) aromatic and hydrophilic AAs are avoided significantly in the psychrophiles as compared to mesophiles.

Table 5 Log odd scores (LOS) of amino acid substitutions calculated using the Eqn (ii).


Our objective in this study was to analyze systematically the compositional variation and substitution preferences of amino acids in proteomes of psychrophiles compared to the proteomes of mesophiles to investigate general proteome wide characteristics for cold adaptation. We considered total compositional differences in proteomes as well as compositional differences in their orthologous proteins alone. We performed analysis at different levels, through simple amino acid compositions, student t-test and finally by substitution patterns in their orthologous proteins. Some of the methods we used were previously applied [26, 2931] but not to the complete proteome analysis.

In psychrophiles individual residue compositions show that there is a significant preference for A, D, S and T content and significant avoidance of E and L content and moderate preference for G and avoidance for F and K content (Table 1). All these residue preferences and avoidance directly show a strong correlation with respect to avoidance for helical content in psychrophiles, as S, D and G are helix breakers [32] and T is a helix indifferent. Likewise, the presence of E tends to favor formation of helical structures and L tends to stabilize helical structures [32] that are highly avoided in psychrophiles. Amino acid D is observed to be unstable at high temperatures and therefore its frequency observed to decrease as optimal growth temperature of organisms increase [33]. Reverse trends are observed for E to counter the trend in favor of making ion pair interactions to form salt-bridges at higher temperatures [34, 35]. Helix destabilizing beta-branched residues (I, T and V) are preferred in beta sheets and loop regions of psychrophilic proteins [36, 37]. The substitution pattern in the orthologous proteins of two temperature groups show several interesting features that are not readily seen in the simple AA compositions. On other hand, they strongly support observed differences in compositions apart from giving additional insights, such as what specific substitutions were more favored or avoided as shown by LOS values.

Nonpolar Amino Acids

Our results in this study confirm that overall composition of nonpolar AA group did not show much difference in proteins belonging to psychrophiles and mesophiles, but there is a decrease in nonpolar AA group frequency in helices and a significant increase in loop regions among proteins of psychrophiles. This is in accordance with earlier findings that there are more nonpolar amino acids on the exposed surface area of the majority of psychrophilic proteins [35, 38] as more loops are observed on surface regions. Among nonpolar residues, I, L and V belong to aliphatic group of residues, which are significantly reduced to favor protein flexibility. It has been widely accepted that the aliphatic amino acids would contribute to the hydrophobic interaction for maintaining conformational stability and rigidity in core region of the proteins [36, 39, 40]. Observed low average hydropathy and low aliphatic residues in psychrophiles are mainly contributed by significantly low in L composition [35].

Tiny and Small Amino acids

Tiny and small amino acids are those with short side chains and are unable to participate in long range interactions among secondary structural elements and are usually confined to form local interactions. Overall their compositions are significantly increased in beta sheets and loops of psychrophilic proteins over mesophilic counter parts. This is also clear from substitutions observed in orthologous proteins. The amino acid G is devoid of side chain, is more flexible with greater rotational freedom, is capable of making cavities in the core parts of the proteins structures [39] and was shown to be in less frequency in thermophiles [41]. Whereas P with pyrolidine ring structure has restricted conformations and was shown to occur in higher frequency in thermophiles [40]. Our LOS values confirm that the amino acid G is preferred and the P is avoided in psychrophiles as compared to mesophiles.

Charged amino acids

Charged residues are polar and hydrophilic. They contribute to ion pair electrostatic interactions that are important binding force for maintaining conformational stability in surface of the proteins [37, 39, 42]. The more charged residues were found in thermophilic proteins than in mesophilic proteins [36, 40]. The charged AA (especially basic and hydrophilic) residues are significantly avoided in psychrophiles (Table 5). Present analysis also supports the notion that these residues are significantly avoided in psychrophiles as observed from compositions as well as LOS scores. The most striking feature of psychrophilic proteins is an increase in amino acid D and decrease in amino acid E over mesophiles. The charged residues in mesophiles are mainly replaced with small and tiny residues in psychrophiles.

Aromatic Amino acids

Psychrophilic proteins and their secondary structural elements show significant decrease in aromatic amino acids. Especially F and Y that are capable of binding to cationic amino acid side chains of K and R in forming cation-π interactions and play important role in stabilizing three dimensional structure of proteins [4345]. Our studies show that, a significant decrease in aromatic residues occurs by substitution with tiny/small and neutral amino acids in psychrophiles (Table 5). Finally, our LOS values, combining with PAM or BLOSUM mutual substitution scores, assist in selecting suitable mutations in designing mesophilic proteins to optimally function in cold temperatures or vice versa. For example substitution of E in mesophilic proteins with D may not change chemical characteristics significantly but may result in optimizing protein function in cold temperatures. In similar lines we observe that amino acid K is highly avoided and R is preferred in psychrophiles although both being basic amino acids.


We analyzed compositions of individual amino acid residues, amino acid groups and their distribution pattern in secondary structures and then computed and quantified their substitution patterns and directional preferences between mesophiles and psychrophiles. Significant differences in composition of amino acid residues were observed between the mesophilic and psychrophilic proteins as summarized below: (i) we observed an increase in frequency of individual amino acids like A, D, S and T that avoid helices and a decrease in E and L amino acids; (ii) There is an increase in small/tiny and neutral group residues which contribute to protein flexibility and a decrease in charged amino acids, particularly basic as well as hydrophilic residues that contribute to ionic interactions; (iii) there is a decrease in aromatic amino acids residues that contribute to the cation-π interactions; (iv) there is a decrease in aliphatic residues which provide good covering and masking to produce hydrophobic pockets that are involved in stabilizing protein structure; (v) there is a reduction in amino acid preferences for helices and an increase of coil forming residues; (vi) we also observed a significant level of substitutions of aliphatic and charged amino acids by tiny/small or neutral amino acids; (vii) the results from this analysis, especially significant t-values of LOS substitution pairs, can be used as a knowledge base in rational design of mutations for engineering of mesophilic proteins to function optimally in cold temperatures or vice versa.


Proteome sequences from six members each of available completely sequenced species of psychrophiles and mesophiles (listed in the foot note of Table 1) are collected. They were selected randomly from independent genuses of mesophiles and 6 of 9 available completely sequenced genomes of psychrophiles to control plausible variations from phylogenetic non-independence (PNI) as related species may share the similar traits due to shared ancestry. These proteome sequences were downloaded from the NCBI proteome project server in the fasta format. The growth temperatures of these species were obtained from NCBI [46] and/or PGTdb (Prokaryotic Growth Temperature database) [47]. We computed the frequencies of amino acid residues in the protein sequences of psychrophilic and mesophilic proteomes. We also grouped the amino acids into 12 property groups [48] as follows: Acidic amino acids group include D and E; aliphatic: I, L and V; aromatic: H, F, W and Y; basic: R, H, and K; charged: R, D, E, H and K; hydrophilic: D, E, K, N, Q and R; hydrophobic: A, C, F, I, L, M, V, W and Y neutral: G, Q, H, S and T; non-polar: A, C, G, I, L, M, F, P, V, W and Y; polar: R, N, D, E, Q, H, K, S and T; small: A, C, D, G, N, P, S, T, V and tiny: A, C, G, S and T. The sum of frequencies of amino acids that fall in each property group are calculated for psychrophiles and mesophiles and are also compared (Table 1). Some of the amino acids are included in more than one property group.

Student t-test

To compare the means of two groups of data, t-test is essentially a good tool for the signal-to-noise metaphor used in research analysis. IIn present analysis, we compared the mean frequencies of single amino acids, 12 different property groups of amino acids and three secondary structural elements from all protein sequences of psychrophilic and mesophilic proteomes considered in this analysis. The t-values are calculated as follows:

t = F P s y c h r o F M e s o ( V a r P s y c h r o / n P s y c h r o ) + ( V a r M e s o / n M e s o ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiDaqNaeyypa0tcfa4aaSaaaeaacqWGgbGrdaWgaaqaaiabdcfaqjabdohaZjabdMha5jabdogaJjabdIgaOjabdkhaYjabd+gaVbqabaGaeyOeI0IaemOray0aaSbaaeaacqWGnbqtcqWGLbqzcqWGZbWCcqWGVbWBaeqaaaqaamaakaaabaWaaSGbaeaacqGGOaakcqWGwbGvcqWGHbqycqWGYbGCdaWgaaqaaiabdcfaqjabdohaZjabdMha5jabdogaJjabdIgaOjabdkhaYjabd+gaVbqabaaabaGaemOBa42aaSbaaeaacqWGqbaucqWGZbWCcqWG5bqEcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBaeqaaaaacqGGPaqkcqGHRaWkdaWcgaqaaiabcIcaOiabdAfawjabdggaHjabdkhaYnaaBaaabaGaemyta0KaemyzauMaem4CamNaem4Ba8gabeaaaeaacqWGUbGBdaWgaaqaaiabd2eanjabdwgaLjabdohaZjabd+gaVbqabaaaaiabcMcaPaqabaaaaaaa@6ECC@

Where Var Psychro and Var Meso are the variance of residues or property groups; (F Psychro ) and (F Meso ) are mean frequencies of psychrophilic and mesophilic proteomes respectively. The n Psychro and n Meso are the total number of psychrophilic and mesophilic proteomes investigated in this study, respectively. Based on student's t-distribution table of significance, critical values for such t-test at various probabilities are as follows (see Table 6):

Table 6 t-test critical values

If t-value is positive and greater than critical value at 10% probability (1.372) then the mean frequency (F Psychro ) of psychrophilic proteomes is significantly greater than that of the mesophilic proteomes (F Meso ) at 90% or higher confidence level. If the frequency of residue or property group t-value is negative and less than -1.372 then the mean frequency of psychrophilic proteomes (F Psychro ) is significantly less than that of mesophilic proteomes (F Meso ) at 90% or higher confidence level [18, 40, 49].

Secondary Structure Prediction

We predicted secondary structural elements in protein sequences using GTOP [50] and/or PSIPRED [51]. We used these predictions to compute frequencies of different amino acids and property groups of residues in three major secondary structural regions, helix (H), strand (E) and coil (C). PSIPRED is a highly reliable secondary structure prediction method with ~83% reported prediction accuracy. We have also tested its prediction accuracy on some of the known psychrophilic and mesophilic proteins to see if there is any significant difference in its prediction for psychrophilic proteins. A total of about 25 proteins, approximately 5000 residues, from each we observed that PSIPRED prediction accuracy was 78.42% and 80.89% for psychrophilic and mesophilic proteins, respectively.

Comparative Analysis of Amino Acid Substitutions

All protein sequences from each mesophilic species in dataset were searched against each proteome of psychrophylic species and vice versa, using BLASTP [52] with 10-3 expectation value cutoff and considerable length coverage. We picked up pair-wise alignments obtained from BLAST results of each protein sequence in a query proteome that showed best hit homolog (ortholog) in the subject proteome. The pair-wise alignments (exclusion of gapped regions) were parsed to calculate amino acid substitution counts between the two proteins from respective proteomes. The substitution counts were normalized to total amino acids present in their respective proteomes pairs individually and finally to all the pairs. The resultant frequency of substitutions was further used to calculate two types of likelihood log odd scores (LOS):

L O S M e s o = L o g F ( X M e s o Y P s y c h r o ) F ( Y M e s o X M e s o ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemitaWKaem4ta8Kaem4uam1aaSbaaSqaaiabd2eanjabdwgaLjabdohaZjabd+gaVbqabaGccqGH9aqpcqWGmbatcqWGVbWBcqWGNbWzjuaGdaWcaaqaaiabdAeagjabcIcaOiabdIfaynaaBaaabaGaemyta0KaemyzauMaem4CamNaem4Ba8gabeaacqGHsgIRcqWGzbqwdaWgaaqaaiabdcfaqjabdohaZjabdMha5jabdogaJjabdIgaOjabdkhaYjabd+gaVbqabaGaeiykaKcabaGaemOrayKaeiikaGIaemywaK1aaSbaaeaacqWGnbqtcqWGLbqzcqWGZbWCcqWGVbWBaeqaaiabgkziUkabdIfaynaaBaaabaGaemyta0KaemyzauMaem4CamNaem4Ba8gabeaacqGGPaqkaaaaaa@62E6@
L O S P s y c h r o = L o g F ( X P s y c h r o Y M e s o ) F ( X P s y c h r o Y P s y c h r o ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemitaWKaem4ta8Kaem4uam1aaSbaaSqaaiabdcfaqjabdohaZjabdMha5jabdogaJjabdIgaOjabdkhaYjabd+gaVbqabaGccqGH9aqpcqWGmbatcqWGVbWBcqWGNbWzjuaGdaWcaaqaaiabdAeagjabcIcaOiabdIfaynaaBaaabaGaemiuaaLaem4CamNaemyEaKNaem4yamMaemiAaGMaemOCaiNaem4Ba8gabeaacqGHsgIRcqWGzbqwdaWgaaqaaiabd2eanjabdwgaLjabdohaZjabd+gaVbqabaGaeiykaKcabaGaemOrayKaeiikaGIaemiwaG1aaSbaaeaacqWGqbaucqWGZbWCcqWG5bqEcqWGJbWycqWGObaAcqWGYbGCcqWGVbWBaeqaaiabgkziUkabdMfaznaaBaaabaGaemiuaaLaem4CamNaemyEaKNaem4yamMaemiAaGMaemOCaiNaem4Ba8gabeaacqGGPaqkaaaaaa@6FAF@

Where F(X Meso Y Psychro ) represent normalized frequency of amino acid X in mesophile substituted by an amino acid Y in psychrophile. The LOS values are calculated by using background substitution frequencies among the mesophilic and/or psychrophilic proteomes in the denominator. The LOS scores, therefore, indicated the pattern of substitutions that are predominantly due to their thermal adaptation and therefore minimize the effect of substitutions due to any speciation events in the evolution process.


  1. Madigan MT, Marrs BL: Extremophiles. Sci Am. 1997, 276 (4): 82-87.

    Article  CAS  PubMed  Google Scholar 

  2. Margesin R, Schinner F, Marx JC, Gerday C: Psychrophiles: From Biodiversity to Biotechnology. 2008, London: Springer, XVI:

    Chapter  Google Scholar 

  3. D'Amico S, Collins T, Marx JC, Feller G, Gerday C: Psychrophilic microorganisms: challenges for life. EMBO Rep. 2006, 7 (4): 385-389. 10.1038/sj.embor.7400662.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Feller G, Gerday C: Psychrophilic enzymes: hot topics in cold adaptation. Nat Rev Microbiol. 2003, 1 (3): 200-208. 10.1038/nrmicro773.

    Article  CAS  PubMed  Google Scholar 

  5. Atomi H: Recent progress towards the application of hyperthermophiles and their enzymes. Curr Opin Chem Biol. 2005, 9 (2): 166-173. 10.1016/j.cbpa.2005.02.013.

    Article  CAS  PubMed  Google Scholar 

  6. Tutino ML, Duilio A, Parrilli R, Remaut E, Sannia G, Marino G: A novel replication element from an Antarctic plasmid as a tool for the expression of proteins at low temperature. Extremophiles. 2001, 5 (4): 257-264. 10.1007/s007920100203.

    Article  CAS  PubMed  Google Scholar 

  7. Cavicchioli R, Siddiqui KS, Andrews D, Sowers KR: Low-temperature extremophiles and their applications. Curr Opin Biotechnol. 2002, 13 (3): 253-261. 10.1016/S0958-1669(02)00317-8.

    Article  CAS  PubMed  Google Scholar 

  8. Marshall CJ: Cold-adapted enzymes. Trends Biotechnol. 1997, 15 (9): 359-364. 10.1016/S0167-7799(97)01086-X.

    Article  CAS  PubMed  Google Scholar 

  9. Russell NJ: Toward a molecular understanding of cold activity of enzymes from psychrophiles. Extremophiles. 2000, 4 (2): 83-90. 10.1007/s007920050141.

    Article  CAS  PubMed  Google Scholar 

  10. Zhou XX, Wang YB, Pan YJ, Li WF: Differences in amino acids composition and coupling patterns between mesophilic and thermophilic proteins. Amino Acids. 2007, 34 (1): 25-33. 10.1007/s00726-007-0589-x.

    Article  PubMed  Google Scholar 

  11. Aghajari N, Feller G, Gerday C, Haser R: Crystallization and preliminary X-ray diffraction studies of alpha-amylase from the antarctic psychrophile Alteromonas haloplanctis A23. Protein Sci. 1996, 5 (10): 2128-2129. 10.1002/pro.5560051021.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Villeret V, Chessa JP, Gerday C, Van Beeumen J: Preliminary crystal structure determination of the alkaline protease from the Antarctic psychrophile Pseudomonas aeruginosa. Protein Sci. 1997, 6 (11): 2462-2464. 10.1002/pro.5560061121.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Alvarez M, Zeelen JP, Mainfroid V, Rentier-Delrue F, Martial JA, Wyns L, Wierenga RK, Maes D: Triose-phosphate isomerase (TIM) of the psychrophilic bacterium Vibrio marinus. Kinetic and structural properties. J Biol Chem. 1998, 273 (4): 2199-2206. 10.1074/jbc.273.4.2199.

    Article  CAS  PubMed  Google Scholar 

  14. Kim SY, Hwang KY, Kim SH, Sung HC, Han YS, Cho Y: Structural basis for cold adaptation. Sequence, biochemical properties, and crystal structure of malate dehydrogenase from a psychrophile Aquaspirillium arcticum. J Biol Chem. 1999, 274 (17): 11761-11767. 10.1074/jbc.274.17.11761.

    Article  CAS  PubMed  Google Scholar 

  15. Gianese G, Bossa F, Pascarella S: Comparative structural analysis of psychrophilic and meso- and thermophilic enzymes. Proteins. 2002, 47 (2): 236-249. 10.1002/prot.10084.

    Article  CAS  PubMed  Google Scholar 

  16. Tronelli D, Maugini E, Bossa F, Pascarella S: Structural adaptation to low temperatures – analysis of the subunit interface of oligomeric psychrophilic enzymes. Febs J. 2007, 274 (17): 4595-4608. 10.1111/j.1742-4658.2007.05988.x.

    Article  CAS  PubMed  Google Scholar 

  17. Jenney FE, Adams MW: The impact of extremophiles on structural genomics (and vice versa). Extremophiles. 2008, 12 (1): 39-50. 10.1007/s00792-007-0087-9.

    Article  CAS  PubMed  Google Scholar 

  18. Jahandideh S, Abdolmaleki P, Jahandideh M, Barzegari Asadabadi E: Sequence and structural parameters enhancing adaptation of proteins to low temperatures. J Theor Biol. 2007, 246 (1): 159-166. 10.1016/j.jtbi.2006.12.008.

    Article  CAS  PubMed  Google Scholar 

  19. Jahandideh M, Barkooie SM, Jahandideh S, Abdolmaleki P, Movahedi MM, Hoseini S, Asadabadi EB, Jouni FJ, Karami Z, Firoozabadi NH: Elucidating the protein cold-adaptation: Investigation of the parameters enhancing protein psychrophilicity. J Theor Biol. 2008, 255 (1): 113-118. 10.1016/j.jtbi.2008.07.034.

    Article  CAS  PubMed  Google Scholar 

  20. Goldstein RA: Amino-acid interactions in psychrophiles, mesophiles, thermophiles, and hyperthermophiles: insights from the quasi-chemical approximation. Protein Sci. 2007, 16 (9): 1887-1895. 10.1110/ps.072947007.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Medigue C, Krin E, Pascal G, Barbe V, Bernsel A, Bertin PN, Cheung F, Cruveiller S, D'Amico S, Duilio A: Coping with cold: the genome of the versatile marine Antarctica bacterium Pseudoalteromonas haloplanktis TAC125. Genome Res. 2005, 15 (10): 1325-1335. 10.1101/gr.4126905.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Methe BA, Nelson KE, Deming JW, Momen B, Melamud E, Zhang X, Moult J, Madupu R, Nelson WC, Dodson RJ: The psychrophilic lifestyle as revealed by the genome sequence of Colwellia psychrerythraea 34H through genomic and proteomic analyses. Proc Natl Acad Sci USA. 2005, 102 (31): 10913-10918. 10.1073/pnas.0504766102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Bult CJ, White O, Olsen GJ, Zhou L, Fleischmann RD, Sutton GG, Blake JA, FitzGerald LM, Clayton RA, Gocayne JD: Complete genome sequence of the methanogenic archaeon, Methanococcus jannaschii. Science. 1996, 273 (5278): 1058-1073. 10.1126/science.273.5278.1058.

    Article  CAS  PubMed  Google Scholar 

  24. DeLong EF: Extreme genomes. Genome Biol. 2000, 1 (6): REVIEWS1029-10.1186/gb-2000-1-6-reviews1029.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. De Vendittis E, Castellano I, Cotugno R, Ruocco MR, Raimo G, Masullo M: Adaptation of model proteins from cold to hot environments involves continuous and small adjustments of average parameters related to amino acid composition. J Theor Biol. 2008, 250 (1): 156-171. 10.1016/j.jtbi.2007.09.006.

    Article  CAS  PubMed  Google Scholar 

  26. Saunders NF, Thomas T, Curmi PM, Mattick JS, Kuczek E, Slade R, Davis J, Franzmann PD, Boone D, Rusterholtz K: Mechanisms of thermal adaptation revealed from the genomes of the Antarctic Archaea Methanogenium frigidum and Methanococcoides burtonii. Genome Res. 2003, 13 (7): 1580-1588. 10.1101/gr.1180903.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Pasamontes A, Garcia-Vallve S: Use of a multi-way method to analyze the amino acid composition of a conserved group of orthologous proteins in prokaryotes. BMC Bioinformatics. 2006, 7: 257-10.1186/1471-2105-7-257.

    Article  PubMed Central  PubMed  Google Scholar 

  28. Nakashima H, Fukuchi S, Nishikawa K: Compositional Changes in RNA, DNA and Proteins for Bacterial Adaptation to Higher and Lower Temperatures. J Biochem. 2003, 3 (4): 507-513. 10.1093/jb/mvg067.

    Article  Google Scholar 

  29. Chan SK, Hsing M, Hormozdiari F, Cherkasov A: Relationship between insertion/deletion (indel) frequency of proteins and essentiality. BMC Bioinformatics. 2007, 8: 227-10.1186/1471-2105-8-227.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Kreil DP, Ouzounis CA: Identification of thermophilic species by the amino acid compositions deduced from their genomes. Nucleic Acids Res. 2001, 29 (7): 1608-1615. 10.1093/nar/29.7.1608.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Tekaia F, Yeramian E, Dujon B: Amino acid composition of genomes, lifestyles of organisms, and evolutionary trends: a global picture with correspondence analysis. Gene. 2002, 297 (1–2): 51-60. 10.1016/S0378-1119(02)00871-5.

    Article  CAS  PubMed  Google Scholar 

  32. Lewis PN, Go N, Go M, Kotelchuck D, Scheraga HA: Helix probability profiles of denatured proteins and their correlation with native structures. Proc Natl Acad Sci USA. 1970, 65 (4): 810-815. 10.1073/pnas.65.4.810.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Szilagyi A, Zavodszky P: Structural differences between mesophilic, moderately thermophilic and extremely thermophilic protein subunits: results of a comprehensive survey. Structure. 2000, 8 (5): 493-504. 10.1016/S0969-2126(00)00133-7.

    Article  CAS  PubMed  Google Scholar 

  34. Gromiha MM, Oobatake M, Sarai A: Important amino acid properties for enhanced thermostability from mesophilic to thermophilic proteins. Biophys Chem. 1999, 82 (1): 51-67. 10.1016/S0301-4622(99)00103-9.

    Article  CAS  PubMed  Google Scholar 

  35. Zhou XX, Wang YB, Pan YJ, Li WF: Differences in amino acids composition and coupling patterns between mesophilic and thermophilic proteins. Amino Acids. 2008, 34 (1): 25-33. 10.1007/s00726-007-0589-x.

    Article  CAS  PubMed  Google Scholar 

  36. Chakravarty S, Varadarajan R: Elucidation of determinants of protein stability through genome sequence analysis. FEBS Lett. 2000, 470 (1): 65-69. 10.1016/S0014-5793(00)01267-9.

    Article  CAS  PubMed  Google Scholar 

  37. Kumar S, Tsai CJ, Nussinov R: Factors enhancing protein thermostability. Protein Eng. 2000, 13 (3): 179-191. 10.1093/protein/13.3.179.

    Article  CAS  PubMed  Google Scholar 

  38. Siddiqui KS, Cavicchioli R: Cold-Adapted Enzymes. Annu Rev Biochem. 2006

    Google Scholar 

  39. Creighton TE: Proteins: structures and molecular properties. 1997, New York: Freeman & Company, 2

    Google Scholar 

  40. Pack SP, Yoo YJ: Protein thermostability: structure-based difference of amino acid between thermophilic and mesophilic proteins. J Biotechnol. 2004, 111 (3): 269-277. 10.1016/j.jbiotec.2004.01.018.

    Article  CAS  PubMed  Google Scholar 

  41. Panasik N, Brenchley JE, Farber GK: Distributions of structural features contributing to thermostability in mesophilic and thermophilic alpha/beta barrel glycosyl hydrolases. Biochim Biophys Acta. 2000, 1543 (1): 189-201.

    Article  CAS  PubMed  Google Scholar 

  42. Vogt G, Woell S, Argos P: Protein thermal stability, hydrogen bonds, and ion pairs. J Mol Biol. 1997, 269 (4): 631-643. 10.1006/jmbi.1997.1042.

    Article  CAS  PubMed  Google Scholar 

  43. Gromiha MM, Suwa M: Structural analysis of residues involving cation-pi interactions in different folding types of membrane proteins. Int J Biol Macromol. 2005, 35 (1–2): 55-62. 10.1016/j.ijbiomac.2004.12.001.

    Article  CAS  PubMed  Google Scholar 

  44. Gallivan JP, Dougherty DA: Cation-pi interactions in structural biology. Proc Natl Acad Sci USA. 1999, 96 (17): 9459-9464. 10.1073/pnas.96.17.9459.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Gromiha MM, Santhosh C, Ahmad S: Structural analysis of cation-pi interactions in DNA binding proteins. Int J Biol Macromol. 2004, 34 (3): 203-211. 10.1016/j.ijbiomac.2004.04.003.

    Article  CAS  PubMed  Google Scholar 

  46. Jenuth JP: The NCBI. Publicly available tools and resources on the Web. Methods Mol Biol. 2000, 132: 301-312.

    CAS  PubMed  Google Scholar 

  47. Huang SL, Wu LC, Liang HK, Pan KT, Horng JT, Ko MT: PGTdb: a database providing growth temperatures of prokaryotes. Bioinformatics. 2004, 20 (2): 276-278. 10.1093/bioinformatics/btg403.

    Article  CAS  PubMed  Google Scholar 

  48. Taylor WR: The classification of amino acid conservation. J Theor Biol. 1986, 119 (2): 205-218. 10.1016/S0022-5193(86)80075-3.

    Article  CAS  PubMed  Google Scholar 

  49. Vogt G, Etzold T, Argos P: An assessment of amino acid exchange matrices in aligning protein sequences: the twilight zone revisited. J Mol Biol. 1995, 249 (4): 816-831. 10.1006/jmbi.1995.0340.

    Article  CAS  PubMed  Google Scholar 

  50. Kawabata T, Fukuchi S, Homma K, Ota M, Araki J, Ito T, Ichiyoshi N, Nishikawa K: GTOP: a database of protein structures predicted from genome sequences. Nucleic Acids Res. 2002, 30 (1): 294-298. 10.1093/nar/30.1.294.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Jones DT: Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol. 1999, 292 (2): 195-202. 10.1006/jmbi.1999.3091.

    Article  CAS  PubMed  Google Scholar 

  52. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


This work was supported by a grant from Howard Hughes Medical Institute to Queens College, CUNY, and Queens College Research Enhancement Grant. We also acknowledge Dormitory Authority, State of New York.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Boojala Vijay B Reddy.

Additional information

Authors' contributions

MRPR and BVBR conceived of the study and participated in its design and coordination. MRPR carried out the analysis and authored the first draft of this manuscript. BVBR was involved in useful discussions and provided comments and revisions to the final version of this text. Both authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: T-test values for amino acid substitution preferences between proteins of Mesophilic and Psychrophilic proteomes. The t-values calculated using LOS substitution scores of mutation frequencies normalized with substitution frequencies within the (a) mesophiles (m2m) and (b) psychrophiles (p2p). (DOC 156 KB)


Additional file 2: T-test values for amino acid property groups' substitution preferences between proteins of Mesophilic and Psychrophilic proteomes. The t-values calculated using LOS substitution scores of mutation frequencies normalized with substitution frequencies within the mesophiles (DOC 60 KB)

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Metpally, R.P.R., Reddy, B.V.B. Comparative proteome analysis of psychrophilic versus mesophilic bacterial species: Insights into the molecular basis of cold adaptation of proteins. BMC Genomics 10, 11 (2009).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: