A systematic comparative analysis of 21 psychrophilic enzymes belonging to different structural families from prokaryotic and eukaryotic organisms is reported. The sequences of these enzymes were multiply aligned to 427 homologous proteins from mesophiles and thermophiles. The net flux of amino acid exchanges from meso/thermophilic to psychrophilic enzymes was measured. To assign the observed preferred exchanges to different structural environments, such as secondary structure, solvent accessibility and subunit interfaces, homology modeling was utilized to predict the secondary structure and accessibility of amino acid residues for the psychrophilic enzymes for which no experimental three-dimensional structure is available. Our results show a clear tendency for the charged residues Arg and Glu to be replaced at exposed sites on alpha-helices by Lys and Ala, respectively, in the direction from 'hot' to 'cold' enzymes. Val is replaced by Ala at buried regions in alpha-helices. Compositional analysis of psychrophilic enzymes shows a significant increase in Ala and Asn and a decrease in Arg at exposed sites. Buried sites in beta-strands tend to be depleted of VAL: Possible implications of the observed structural variations for protein stability and engineering are discussed.