BACKGROUND
life exists over a wide temperature range, from as low as -15°c to as high as 122°c  <cit> .
on the upper end of the temperature spectrum, thermophiles and hyperthermophiles have been studied extensively by the scientific community, particularly the molecular mechanisms that support protein structure and function at high temperatures.
for example, compact and strong hydrophobic packing is typically found in most cores of thermophilic proteins, which increases the energy needed to unfold the protein, making it possible for thermophilic proteins to retain native structure at high temperatures  <cit> .
indeed, a strong correlation exists between high core packing density and thermostability  <cit> .
we are also intrigued by the obverse - organisms that survive optimally at cold temperatures; harsh environments with restricted molecular mobility and reduced reaction kinetics that hinder myriad cellular and biomolecular processes  <cit> .
psychrophiles, "cold-loving" microorganisms, have adapted to life at low temperatures by using a variety of mechanisms.
these include the production of anti-freeze and cold-shock proteins, alterations in membrane composition, and overexpression of proteins that destabilize dna structures, among other mechanisms  <cit> .
some of the more interesting adaptations are found in psychrophilic enzymes.
in particular, psychrophilic enzymes typically have a higher occurrence of nonpolar residues on their surface, which destabilizes water structure around the enzyme such that water properties are considerably distinct from bulk water, and a lower abundance of arginine and proline residues, which increases backbone flexibility.
moreover, fewer weak interactions, such as hydrogen bonds, reduce the degree of packing within the enzyme core and make psychrophilic enzymes more conformationally flexible .
increased flexibility has been proposed as a major molecular mechanism for the evolution of cold adapted enzymes, a hypothesis supported by spectroscopic analysis  simulations  <cit> ).
psychrophilic enzymes have been shown to lose more entropy upon transition state activation than their mesophilic counterparts, suggesting that the psychrophilic proteins exist in a more disordered ground state  <cit> .
this intuitively makes sense, since in order for psychrophilic enzymes to be more active at low temperatures they must have more kinetic energy  <cit> .
nonetheless, the structural basis of this flexibility remains unclear.
for example, an md simulations study of one-on-one comparisons of trypsin from atlantic salmon and its mesophilic counterpart  <cit>  found that the psychrophilic trypsin possesses higher flexibility close to the active site parts.
interestingly, the active site is where there is the most significant difference in amino acid composition between mesophilic and psychrophilic trypsins.
however, other psychrophilic enzymes exhibit regions of flexibility in places distant from the active site.
for example, psychrophilic uracil-dna glycosylases  <cit>  have fewer strong ion-pair interactions near the c-terminus than their mesophilic counterparts.
this local difference results in increased flexibility in the psychrophilic enzymes.
psychrophilic alpha-amylases  <cit>  show greater overall flexibility than the mesophilic/thermophilic counterparts.
this is consistent with thermal unfolding experiments, where the psychrophilic enzymes unfold at low temperatures cooperatively without intermediates due to fewer stabilizing interactions, while thermophilic amylases show intermediates during unfolding, indicating that there are regions with greater rigidity than others.
the different strategies that psychrophilic enzymes use to adapt to colder environments have resulted in a number of divergent viewpoints on the influence of local and global protein flexibility on cold adaptation.
an interesting structural property of proteins, particularly as it relates to thermostability, is cavity number and size.
protein thermostability appears to diminish when cavities are created in the protein  <cit> .
such cavities represent packing defects in the protein core.
in the case of psychrophilic proteins, only a few studies have addressed these cavities, and then only with a small set of proteins  <cit> .
no clear relationship was observed between psychrophily and total void-volume in the enzyme core, as might be expected from the strong hydrophobic cores present in thermophilic proteins.
however, an increased number of three-dimensional structures of psychrophilic proteins have become available and offer the opportunity to revisit such a structural analysis.
here we test the hypothesis that cavity properties in psychrophilic enzymes endow these proteins with the increased conformationally flexibility necessary to function optimally at low temperatures.
in the current work we investigated differences in the average cavity volume and crystallographic waters between psychrophilic and mesophilic homologs.
a non-redundant set of  <dig> psychrophilic enzymes was examined, with each paired to a homologous mesophilic enzyme , and all with high resolution crystal structures.
in addition to counting cavities and calculating void volumes, we evaluated amino acid frequencies in residue positions surrounding cavities.
we present evidence that the average cavity size of psychrophilic enzymes is larger, and contains more surrounding hydrophilic groups, than their mesophilic counterparts.
these findings support a hypothesis that psychrophilic enzymes may have a predisposition to having more water molecules within their cavities, which consequently increases enzyme conformational dynamics, leading to greater activity under the rigidifying cold environment.
the results provide a potential strategy by which optimal temperature could be altered.
RESULTS
database construction
a database of psychrophilic and mesophilic homologous proteins was constructed from several literature references and an in-house pipeline script .
all  <dig> psychrophilic enzymes from siddiqui et al.
<cit>  were included.
additionally, eight psychrophilic enzymes were included from the ncbi database through query using the ncbi entrez tool.
a total of  <dig> mesophilic proteins that were homologous to the  <dig> psychrophilic proteins were then obtained with dalilite.
the mesophilic proteins were then rejected if they had resolution >  <dig>  Å, did not belong to the same enzyme family as their homologous psychrophilic enzyme, or had less than 30% sequence identity.
dataset of  <dig> psychrophilic/mesophilic protein pairs, pdb code, species name, and structure resolution are provided.
b-values
psychrophilic enzymes sacrifice conformational stability to become more flexible, and with greater intrinsic flexibility they remain catalytically active at lower temperatures  <cit> .
nonetheless, the structural basis for this flexibility remains unclear.
to investigate psychrophilic protein flexibility on a molecular level, we carried out a statistical comparison of crystallographic b-values between psychrophilic and homologous mesophilic proteins.
b-values are measurements of atomic fluctuations in protein crystals.
atomic fluctuations can be divided into dynamic motions and crystal lattice defects, but only the former is relevant to flexibility  <cit> .
to isolate the atomic motion component from the lattice component, we first eliminated outliers from the raw b-values  and then normalized the resulting b-values based on the overall mean and standard deviation.
this approach, which has been used previously  <cit>  factors out overall differences by assuming the same overall distribution of atomic flexibility across proteins, allowing us to examine the regional distribution of b-value differences.
using the normalized b-values and a structure-based alignment, we calculated Δb'i = b'i-b'i for all aligned positions i, as observed in additional file  <dig>  figure s <dig> for the pair 1elt and 1eai.
relatively high values of Δb'i indicate regions of higher flexibility in the psychrophilic protein relative to the mesophilic homolog.
distributions of Δb'i across secondary structure as defined by the dictionary of protein secondary structure  revealed interesting correlations .
a plot of the calculated Δb' values for each secondary structure can be seen in figure  <dig> of the distribution of all  <dig> homolog pairs.
the strand  and the 5-turn  secondary structures were significantly more flexible in psychrophilic proteins as compared to their mesophilic counterparts.
following a similar methodology, siglioccolo et al.
<cit>  also analyzed b'-values of psychrophilic and mesophilic proteins grouped by common secondary structures .
they observed that β-sheets and turns tend to be more flexible in psychrophilic proteins, relative to helices.
the results hold true for this new, larger set of homologous pairs.
because overall flexibility differences were factored out by normalization, we cannot distinguish an increased flexibility in β-sheets from a decreased flexibility in α-helices using this data; however, it is worth noting that β-sheets tend to sit in the core of a protein where lattice defects have less of an effect, leaving intrinsic flexibility as the likely cause of higher b-values.
cavity volume and morphology
b-values suggest, but do not inform, structural differences between a psychrophilic protein and its mesophilic counterpart, which prompted us to investigate the intrinsic characteristics of the protein core.
thermophilic globular proteins are known to have a highly compact hydrophobic core  <cit> .
hydrophobic packing in psychrophilic proteins has not been widely studied, but the increased presence of amino acids with smaller side-chains inside the protein, points to a weak hydrophobic core.
the internal loose packing in the protein core relates to intrinsic flexibility, which has been previously noted  <cit> .
paired proteins were analyzed by comparing total void volume , total pocket volume , total cavity volume, total number of voids , total number of pockets, and total number of cavities.
the pockets and cavities were analyzed at probe radii ranging from  <dig>  Å to  <dig>  Å. the data was not normalized by molecular weight because for the most part each protein pair has similar molecular weights.
although it was expected that psychrophilic enzymes would have larger total cavity volumes when compared to their mesophilic counterparts, 40% of the psychrophilic proteins examined had smaller total cavity volumes throughout the entire range of probe sizes.
also, analysis of total number of cavities showed no significant differences between the two populations.
no significant differences were found for either total void volume  or total pocket volume and number.
however, a deeper analysis revealed a correlation wherein the average cavity volume is significantly larger  in psychrophilic proteins at probe sizes equivalent to the size of a water molecule .
this was observed in two distinct ways: 1) when a psychrophilic enzyme had a smaller total cavity void volume than its mesophilic homolog, then the number of cavities in the former was less than in the latter, resulting in a larger average cavity volume in the psychrophilic enzyme; and 2) when a psychrophilic enzyme had more cavities than the homologous mesophilic enzyme, the total cavity volume was larger in the psychrophilic enzyme.
psychrophilic enzymes had an average cavity volume of  <dig> Å <dig> compared to an average cavity volume of  <dig> Å <dig> found in mesophilic enzymes at the  <dig>  Å probe size.
figure  <dig> shows the distribution of the difference in cavity average size between each paired psychrophilic and mesophilic enzymes for all  <dig> pairs at different probe sizes.
our observation that cavities of the size of a water molecule are present in greater volume in psychrophilic enzymes, suggests that water may be bound to psychrophilic enzymes in greater numbers, enhancing the internal solvation.
amino acid frequency surrounding psychrophilic cavities
to better understand cavity composition, inward-facing side chains of amino acids that surround the cavities at probe size  <dig>  Å were compared between psychrophilic and mesophilic proteins .
outward facing side chains were not considered because the backbone chemistry is the same in the two populations.
amino acids were classified into four groups: nonpolar, hydrophobic, basic, and acidic.
it was observed that the average number of acidic amino acid residues lining psychrophilic protein cavities exceed the average number observed in homologous mesophilic protein cavities .
moreover, there are fewer hydrophobic side chains surrounding psychrophilic cavities than mesophilic cavities .
finally, cavities in psychrophilic proteins are more predominant in regions containing turns and coils, in contrast to mesophilic protein cavities .
this study reveals that cavities in psychrophilic enzymes found at probe size  <dig>  to  <dig>  Å:  are statistically significantly different than in homologous mesophilic enzymes; and  their cavity surfaces contain a higher proportion of acidic amino acids.
several studies have found that cavities with these specific characteristics allow water molecules to exist stably within them  <cit> .
two key characteristics, in particular, are high polarity and excess space for movement.
for example, park and saven  <cit>  analyzed  <dig>  buried water molecules contained in  <dig> different high resolution protein structures, revealing that these water molecules formed hydrogen bonds with polar atoms, predominantly near residues that compose turns or coils.
chen and stites  <cit>  obtained a similar result both experimentally and computationally with staphylococcal nuclease, wherein a water molecule stabilized glu- <dig>  these authors concluded that water molecules in polar cavities make more stable hydrogen bonds with the cavity walls and have longer residence times than in hydrophobic cavities.
however, the presence of water and its function in protein cavities remains unclear.
it is well established that overall hydration increases the flexibility of the protein  <cit> .
for example, armstrong et al.
<cit> , using electron paramagnetic resonance  spectroscopy, showed the importance of water hydration in the core of apomyoglobin and its role in protein transition between several structural conformational states, presumably by acting as a lubricant.
the dynamics of water inside the protein core affects protein thermostability  <cit> .
another methodology used to understand the role of water in enzymes, specifically the role of water in catalysis, involves nonaqueous systems.
specifically, adding water to a final concentration of 1%  in tetrahydrofuran resulted in a significant increase in the proteolytic activity of subtilisin carlsberg, concomitant with an increase in active site mobility as determined by epr  <cit> .
using molecular dynamics simulations, tarek and tobias  <cit>  demonstrated that higher levels of hydration contribute to increased in protein motion.
similarly, rupley et al.
<cit> , showed that a specific level of hydration is required by proteins to be active .
these findings point to an important correlation between hydration, protein flexibility, and enzymatic activity.
many studies have proposed the importance of cavity hydration and its relation to higher protein flexibility.
fischer et al.
<cit>  theoretically calculated the vibration entropy of bovine pancreatic trypsin inhibitor with bound and unbound water- <dig> .
bound water molecules increase the vibrational entropy of the protein, which could also be thought of as an increase in protein flexibility  <cit> .
garcia and hummer  <cit>  used md simulations to observe that water molecules inside proteins slowly exchange with the solvent, and when the molecules escape or penetrate the protein, they cause dynamic fluctuations.
counting crystallographic waters
the observation that cavities in psychrophiles favor water-sized objects should imply that more waters are bound to the cold-adapted enzymes.
to verify this, we isolated the buried crystallographic waters in all  <dig> enzymes  by removing the surface-exposed waters and any waters that were buried by surface exposed waters, and then counting them.
a correlation was observed  in the number of buried waters between psychrophilic and mesophilic homologs, showing a consistency across crystals and crystallographers in modeling ordered solvent .
there is a visible trend towards more water molecules in psychrophilic enzymes; however, the significance is marginal .
interestingly, counter examples in this set appear to possess large numbers of buried waters.
in addition to 1a <dig>  if we remove the five pairs of psychrophilic-mesophilic homologs with >  <dig> buried waters , then  <dig> out of the remaining  <dig> pairs fit the trend of higher water content in psychrophilic enzymes than in mesophilic enzymes.
nonetheless, ignoring such data pruning, the difference is significant using either the binomial test  or the paired t-test .
moreover, we may not be counting all cavity water molecules because of possible experimental errors or water in cavities being highly mobile   <cit> .
therefore, a more sensitive approach, e.g., md simulations, should be considered to test the hypothesis that psychrophilic enzymes bind more water than their mesophilic counterparts.
CONCLUSIONS
in this study we explored a variety of structural bioinformatic metrics to seek a structural explanation for cold adaptation in enzymes.
the most significant structural differences are an increase in the size of the water-sized cavities and a trend in amino acid composition towards carboxylic acids in these cavities.
through an additional consensus of measures, including a significant increase in crystallographic b-values in β-sheets and turns, and a marginally significant increase in the number of buried crystallographic waters, we can conclude that psychrophilic enzymes tend to be more solvated in the core as compared to mesophilic enzymes.
in particular, the evidence that psychrophilic cavities are well characterized by a water-sized probe suggests that mutations that reshaped internal cavities to fit water may have led to more bound water, which in turn led to an increased flexibility in the core, consistent with water-protein literature.
statistical metrics suffer from small data sets, crystallographic variability, and the fact that multiple mechanisms for cold adaptation exist.
nonetheless, our results point to a common hypothesis that can now be tested by protein design experiments .
cavities are not necessarily the only element of a protein that endows psychrophilic proteins with cold adaptation, but structure-based differences in cavities may reveal themselves to be critical to cold adaptation and might help to design enzymes capable of functioning at low temperatures.
similar results may be obtained in other rigidifying environments, such as organic solvents, polymeric materials  <cit> , or protein-nanomaterial conjugate materials  <cit> .

 life exists over a wide temperature range, from as low as -15°c to as high as 122°c. on the upper end of the temperature spectrum, thermophiles and hyperthermophiles have been studied extensively by the scientific community, particularly the molecular mechanisms that support protein structure and function at high temperatures. for example, compact and strong hydrophobic packing is typically found in most cores of thermophilic proteins, which increases the energy needed to unfold the protein, making it possible for thermophilic proteins to retain native structure at high temperatures. indeed, a strong correlation exists between high core packing density and thermostability. we are also intrigued by the obverse - organisms that survive optimally at cold temperatures; harsh environments with restricted molecular mobility and reduced reaction kinetics that hinder myriad cellular and biomolecular processes. psychrophiles, "cold-loving" microorganisms, have adapted to life at low temperatures by using a variety of mechanisms. these include the production of anti-freeze and cold-shock proteins, alterations in membrane composition, and overexpression of proteins that destabilize dna structures, among other mechanisms. some of the more interesting adaptations are found in psychrophilic enzymes. in particular, psychrophilic enzymes typically have a higher occurrence of nonpolar residues on their surface, which destabilizes water structure around the enzyme such that water properties are considerably distinct from bulk water, and a lower abundance of arginine and proline residues, which increases backbone flexibility. moreover, fewer weak interactions, such as hydrogen bonds, reduce the degree of packing within the enzyme core and make psychrophilic enzymes more conformationally flexible. increased flexibility has been proposed as a major molecular mechanism for the evolution of cold adapted enzymes, a hypothesis supported by spectroscopic analysis simulations ). psychrophilic enzymes have been shown to lose more entropy upon transition state activation than their mesophilic counterparts, suggesting that the psychrophilic proteins exist in a more disordered ground state. this intuitively makes sense, since in order for psychrophilic enzymes to be more active at low temperatures they must have more kinetic energy. nonetheless, the structural basis of this flexibility remains unclear. for example, an md simulations study of one-on-one comparisons of trypsin from atlantic salmon and its mesophilic counterpart found that the psychrophilic trypsin possesses higher flexibility close to the active site parts. interestingly, the active site is where there is the most significant difference in amino acid composition between mesophilic and psychrophilic trypsins. however, other psychrophilic enzymes exhibit regions of flexibility in places distant from the active site. for example, psychrophilic uracil-dna glycosylases have fewer strong ion-pair interactions near the c-terminus than their mesophilic counterparts. this local difference results in increased flexibility in the psychrophilic enzymes. psychrophilic alpha-amylases show greater overall flexibility than the mesophilic/thermophilic counterparts. this is consistent with thermal unfolding experiments, where the psychrophilic enzymes unfold at low temperatures cooperatively without intermediates due to fewer stabilizing interactions, while thermophilic amylases show intermediates during unfolding, indicating that there are regions with greater rigidity than others. the different strategies that psychrophilic enzymes use to adapt to colder environments have resulted in a number of divergent viewpoints on the influence of local and global protein flexibility on cold adaptation. an interesting structural property of proteins, particularly as it relates to thermostability, is cavity number and size. protein thermostability appears to diminish when cavities are created in the protein. such cavities represent packing defects in the protein core. in the case of psychrophilic proteins, only a few studies have addressed these cavities, and then only with a small set of proteins. no clear relationship was observed between psychrophily and total void-volume in the enzyme core, as might be expected from the strong hydrophobic cores present in thermophilic proteins. however, an increased number of three-dimensional structures of psychrophilic proteins have become available and offer the opportunity to revisit such a structural analysis. here we test the hypothesis that cavity properties in psychrophilic enzymes endow these proteins with the increased conformationally flexibility necessary to function optimally at low temperatures. in the current work we investigated differences in the average cavity volume and crystallographic waters between psychrophilic and mesophilic homologs. a non-redundant set of psychrophilic enzymes was examined, with each paired to a homologous mesophilic enzyme, and all with high resolution crystal structures. in addition to counting cavities and calculating void volumes, we evaluated amino acid frequencies in residue positions surrounding cavities. we present evidence that the average cavity size of psychrophilic enzymes is larger, and contains more surrounding hydrophilic groups, than their mesophilic counterparts. these findings support a hypothesis that psychrophilic enzymes may have a predisposition to having more water molecules within their cavities, which consequently increases enzyme conformational dynamics, leading to greater activity under the rigidifying cold environment. the results provide a potential strategy by which optimal temperature could be altered.  database construction a database of psychrophilic and mesophilic homologous proteins was constructed from several literature references and an in-house pipeline script. all psychrophilic enzymes from siddiqui et al. were included. additionally, eight psychrophilic enzymes were included from the ncbi database through query using the ncbi entrez tool. a total of mesophilic proteins that were homologous to the psychrophilic proteins were then obtained with dalilite. the mesophilic proteins were then rejected if they had resolution > Å, did not belong to the same enzyme family as their homologous psychrophilic enzyme, or had less than 30% sequence identity. dataset of psychrophilic/mesophilic protein pairs, pdb code, species name, and structure resolution are provided. b-values psychrophilic enzymes sacrifice conformational stability to become more flexible, and with greater intrinsic flexibility they remain catalytically active at lower temperatures. nonetheless, the structural basis for this flexibility remains unclear. to investigate psychrophilic protein flexibility on a molecular level, we carried out a statistical comparison of crystallographic b-values between psychrophilic and homologous mesophilic proteins. b-values are measurements of atomic fluctuations in protein crystals. atomic fluctuations can be divided into dynamic motions and crystal lattice defects, but only the former is relevant to flexibility. to isolate the atomic motion component from the lattice component, we first eliminated outliers from the raw b-values and then normalized the resulting b-values based on the overall mean and standard deviation. this approach, which has been used previously factors out overall differences by assuming the same overall distribution of atomic flexibility across proteins, allowing us to examine the regional distribution of b-value differences. using the normalized b-values and a structure-based alignment, we calculated Δb'i = b'i-b'i for all aligned positions i, as observed in additional file figure s for the pair 1elt and 1eai. relatively high values of Δb'i indicate regions of higher flexibility in the psychrophilic protein relative to the mesophilic homolog. distributions of Δb'i across secondary structure as defined by the dictionary of protein secondary structure revealed interesting correlations. a plot of the calculated Δb' values for each secondary structure can be seen in figure of the distribution of all homolog pairs. the strand and the 5-turn secondary structures were significantly more flexible in psychrophilic proteins as compared to their mesophilic counterparts. following a similar methodology, siglioccolo et al. also analyzed b'-values of psychrophilic and mesophilic proteins grouped by common secondary structures. they observed that β-sheets and turns tend to be more flexible in psychrophilic proteins, relative to helices. the results hold true for this new, larger set of homologous pairs. because overall flexibility differences were factored out by normalization, we cannot distinguish an increased flexibility in β-sheets from a decreased flexibility in α-helices using this data; however, it is worth noting that β-sheets tend to sit in the core of a protein where lattice defects have less of an effect, leaving intrinsic flexibility as the likely cause of higher b-values. cavity volume and morphology b-values suggest, but do not inform, structural differences between a psychrophilic protein and its mesophilic counterpart, which prompted us to investigate the intrinsic characteristics of the protein core. thermophilic globular proteins are known to have a highly compact hydrophobic core. hydrophobic packing in psychrophilic proteins has not been widely studied, but the increased presence of amino acids with smaller side-chains inside the protein, points to a weak hydrophobic core. the internal loose packing in the protein core relates to intrinsic flexibility, which has been previously noted. paired proteins were analyzed by comparing total void volume, total pocket volume, total cavity volume, total number of voids, total number of pockets, and total number of cavities. the pockets and cavities were analyzed at probe radii ranging from Å to Å. the data was not normalized by molecular weight because for the most part each protein pair has similar molecular weights. although it was expected that psychrophilic enzymes would have larger total cavity volumes when compared to their mesophilic counterparts, 40% of the psychrophilic proteins examined had smaller total cavity volumes throughout the entire range of probe sizes. also, analysis of total number of cavities showed no significant differences between the two populations. no significant differences were found for either total void volume or total pocket volume and number. however, a deeper analysis revealed a correlation wherein the average cavity volume is significantly larger in psychrophilic proteins at probe sizes equivalent to the size of a water molecule. this was observed in two distinct ways: 1) when a psychrophilic enzyme had a smaller total cavity void volume than its mesophilic homolog, then the number of cavities in the former was less than in the latter, resulting in a larger average cavity volume in the psychrophilic enzyme; and 2) when a psychrophilic enzyme had more cavities than the homologous mesophilic enzyme, the total cavity volume was larger in the psychrophilic enzyme. psychrophilic enzymes had an average cavity volume of Å compared to an average cavity volume of Å found in mesophilic enzymes at the Å probe size. figure shows the distribution of the difference in cavity average size between each paired psychrophilic and mesophilic enzymes for all pairs at different probe sizes. our observation that cavities of the size of a water molecule are present in greater volume in psychrophilic enzymes, suggests that water may be bound to psychrophilic enzymes in greater numbers, enhancing the internal solvation. amino acid frequency surrounding psychrophilic cavities to better understand cavity composition, inward-facing side chains of amino acids that surround the cavities at probe size Å were compared between psychrophilic and mesophilic proteins. outward facing side chains were not considered because the backbone chemistry is the same in the two populations. amino acids were classified into four groups: nonpolar, hydrophobic, basic, and acidic. it was observed that the average number of acidic amino acid residues lining psychrophilic protein cavities exceed the average number observed in homologous mesophilic protein cavities. moreover, there are fewer hydrophobic side chains surrounding psychrophilic cavities than mesophilic cavities. finally, cavities in psychrophilic proteins are more predominant in regions containing turns and coils, in contrast to mesophilic protein cavities. this study reveals that cavities in psychrophilic enzymes found at probe size to Å: are statistically significantly different than in homologous mesophilic enzymes; and their cavity surfaces contain a higher proportion of acidic amino acids. several studies have found that cavities with these specific characteristics allow water molecules to exist stably within them. two key characteristics, in particular, are high polarity and excess space for movement. for example, park and saven analyzed buried water molecules contained in different high resolution protein structures, revealing that these water molecules formed hydrogen bonds with polar atoms, predominantly near residues that compose turns or coils. chen and stites obtained a similar result both experimentally and computationally with staphylococcal nuclease, wherein a water molecule stabilized glu- these authors concluded that water molecules in polar cavities make more stable hydrogen bonds with the cavity walls and have longer residence times than in hydrophobic cavities. however, the presence of water and its function in protein cavities remains unclear. it is well established that overall hydration increases the flexibility of the protein. for example, armstrong et al., using electron paramagnetic resonance spectroscopy, showed the importance of water hydration in the core of apomyoglobin and its role in protein transition between several structural conformational states, presumably by acting as a lubricant. the dynamics of water inside the protein core affects protein thermostability. another methodology used to understand the role of water in enzymes, specifically the role of water in catalysis, involves nonaqueous systems. specifically, adding water to a final concentration of 1% in tetrahydrofuran resulted in a significant increase in the proteolytic activity of subtilisin carlsberg, concomitant with an increase in active site mobility as determined by epr. using molecular dynamics simulations, tarek and tobias demonstrated that higher levels of hydration contribute to increased in protein motion. similarly, rupley et al., showed that a specific level of hydration is required by proteins to be active. these findings point to an important correlation between hydration, protein flexibility, and enzymatic activity. many studies have proposed the importance of cavity hydration and its relation to higher protein flexibility. fischer et al. theoretically calculated the vibration entropy of bovine pancreatic trypsin inhibitor with bound and unbound water-. bound water molecules increase the vibrational entropy of the protein, which could also be thought of as an increase in protein flexibility. garcia and hummer used md simulations to observe that water molecules inside proteins slowly exchange with the solvent, and when the molecules escape or penetrate the protein, they cause dynamic fluctuations. counting crystallographic waters the observation that cavities in psychrophiles favor water-sized objects should imply that more waters are bound to the cold-adapted enzymes. to verify this, we isolated the buried crystallographic waters in all enzymes by removing the surface-exposed waters and any waters that were buried by surface exposed waters, and then counting them. a correlation was observed in the number of buried waters between psychrophilic and mesophilic homologs, showing a consistency across crystals and crystallographers in modeling ordered solvent. there is a visible trend towards more water molecules in psychrophilic enzymes; however, the significance is marginal. interestingly, counter examples in this set appear to possess large numbers of buried waters. in addition to 1a if we remove the five pairs of psychrophilic-mesophilic homologs with > buried waters, then out of the remaining pairs fit the trend of higher water content in psychrophilic enzymes than in mesophilic enzymes. nonetheless, ignoring such data pruning, the difference is significant using either the binomial test or the paired t-test. moreover, we may not be counting all cavity water molecules because of possible experimental errors or water in cavities being highly mobile. therefore, a more sensitive approach, e.g., md simulations, should be considered to test the hypothesis that psychrophilic enzymes bind more water than their mesophilic counterparts.  in this study we explored a variety of structural bioinformatic metrics to seek a structural explanation for cold adaptation in enzymes. the most significant structural differences are an increase in the size of the water-sized cavities and a trend in amino acid composition towards carboxylic acids in these cavities. through an additional consensus of measures, including a significant increase in crystallographic b-values in β-sheets and turns, and a marginally significant increase in the number of buried crystallographic waters, we can conclude that psychrophilic enzymes tend to be more solvated in the core as compared to mesophilic enzymes. in particular, the evidence that psychrophilic cavities are well characterized by a water-sized probe suggests that mutations that reshaped internal cavities to fit water may have led to more bound water, which in turn led to an increased flexibility in the core, consistent with water-protein literature. statistical metrics suffer from small data sets, crystallographic variability, and the fact that multiple mechanisms for cold adaptation exist. nonetheless, our results point to a common hypothesis that can now be tested by protein design experiments. cavities are not necessarily the only element of a protein that endows psychrophilic proteins with cold adaptation, but structure-based differences in cavities may reveal themselves to be critical to cold adaptation and might help to design enzymes capable of functioning at low temperatures. similar results may be obtained in other rigidifying environments, such as organic solvents, polymeric materials, or protein-nanomaterial conjugate materials. 
we examine twenty homologous enzyme pairs from psychrophiles and mesophiles to investigate flexibility as a key characteristic for cold adaptation. thus, our results point to the important role enzyme flexibility plays in adaptation to cold environments. their enzymes are active at cold temperatures by being structurally more flexible than mesophilic enzymes. psychrophiles, cold-adapted organisms, have adapted to live at low temperatures by using a variety of mechanisms. enzyme cavities, identified using castp at various probe sizes, indicate that psychrophilic enzymes have larger average cavity sizes at probe radii of - Å, sufficient for water molecules. these findings suggest that embedded water molecules may play a significant role in cavity flexibility, and therefore, overall protein flexibility. even though, there are some indications of the possible structural mechanisms by which psychrophilic enzymes are catalytic active at cold temperatures, there is not a generalized structural property common to all psychrophilic enzymes. comparing psychrophilic to mesophilic protein b-factors reveals that psychrophilic enzymes are more flexible in 5-turn and strand secondary structures. furthermore, amino acid side chains lining these cavities show an increased frequency of acidic groups in psychrophilic enzymes.
 psychrophiles, cold-adapted organisms, have adapted to live at low temperatures by using a variety of mechanisms. their enzymes are active at cold temperatures by being structurally more flexible than mesophilic enzymes. even though, there are some indications of the possible structural mechanisms by which psychrophilic enzymes are catalytic active at cold temperatures, there is not a generalized structural property common to all psychrophilic enzymes.  we examine twenty homologous enzyme pairs from psychrophiles and mesophiles to investigate flexibility as a key characteristic for cold adaptation. b-factors in protein x-ray structures are one way to measure flexibility. comparing psychrophilic to mesophilic protein b-factors reveals that psychrophilic enzymes are more flexible in 5-turn and strand secondary structures. enzyme cavities, identified using castp at various probe sizes, indicate that psychrophilic enzymes have larger average cavity sizes at probe radii of - Å, sufficient for water molecules. furthermore, amino acid side chains lining these cavities show an increased frequency of acidic groups in psychrophilic enzymes.  these findings suggest that embedded water molecules may play a significant role in cavity flexibility, and therefore, overall protein flexibility. thus, our results point to the important role enzyme flexibility plays in adaptation to cold environments. 