Sequence and structural parameters enhancing adaptation of proteins to low temperatures

A systematic analysis compared sequence and structural parameters distributions between 13 pairs of psychrophilic and mesophilic proteins for elucidating the cold adaptation parameters. The results of statistical test (t-test) revealed that helical content, tight turn content, disulfide bonds and hydrogen bonds do not show significant difference between psychrophilic and mesophilic proteins. However, it was demonstrated in this study that a larger proportion of open beta-turn in psychrophilic proteins is an effective parameter in specific activity at low temperature. In addition, substitution of amino acids of charged and aliphatic groups with amino acids of tiny and small groups in protein chains, tight turns and alpha-helices in the direction from mesophilic to psychrophilic proteins is one of the mechanisms of low temperature adaptation. Such sequence and structural parameter differences would help to develop a strategy for designing cold-adapted proteins.     

Molecular characterization of cold adaptation based on ortholog protein sequences from <Emphasis Type="Italic">Vibrionaceae</Emphasis> species

A set of 298 protein families from psychrophilic Vibrio salmonicida was compiled to identify genotypic characteristics that discern it from orthologous sequences from the mesophilic Vibrio/Photobacterium branch of the gamma-Proteobacteria (Vibrionaceae family). In our comparative exploration we employed alignment based bioinformatical and statistical methods. Interesting information was found in the substitution matrices, and the pattern of asymmetries in the amino acid substitution process. Together with the compositional difference, they identified the amino acids Ile, Asn, Ala and Gln as those having the most psycrophilic involvement. Ile and Asn are enhanced whereas Gln and Ala are suppressed. The inflexible Pro residue is also suppressed in loop regions, as expected in a flexible structure. The dataset were also classified and analysed according to the predicted subcellular location, and we made an additional study of 183 intracellular and 65 membrane proteins. Our results revealed that the psychrophilic proteins have similar hydrophobic and charge contributions in the core of the protein as mesophilic proteins, while the solvent-exposed surface area is significantly more hydrophobic. In addition, the psychrophilic intracellular (but not the membrane) proteins are significantly more negatively charged at the surface. Our analysis supports the hypothesis of preference for more flexible amino acids at the molecular surface. Life in cold climate seems to be obtained through many minor structural modifications rather than certain amino acids substitutions.     

Protein psychrophilicity: role of residual structural properties in adaptation of proteins to low temperatures

In order to investigate the structural distribution responsible for protein psychrophilicity, a systematic comparative analysis of 13 pairs of psychrophilic and mesophilic proteins is reported. Three kinds of residue structural states such as exposed, intermediate and buried were considered for analyzing the structural patterns of single amino acids and amino acids in different groups. The statistical test revealed that higher frequency in exposed state of Ala, higher frequency in intermediate state of His, lower frequency in buried state of Lys, lower frequency in exposed state of Gln, higher frequency in exposed state and in intermediate state of Thr, higher frequency in exposed and intermediate state of tiny and small amino acids groups could be critical factors related with protein psychrophilicity. Such structure-based differences of residual properties would help to develop a strategy for designing psychrophilic proteins.     

Adaptation of model proteins from cold to hot environments involves continuous and small adjustments of average parameters related to amino acid composition

The growth temperature adaptation of six model proteins has been studied in 42 microorganisms belonging to eubacterial and archaeal kingdoms, covering optimum growth temperatures from 7 to 103 degrees C. The selected proteins include three elongation factors involved in translation, the enzymes glyceraldehyde-3-phosphate dehydrogenase and superoxide dismutase, the cell division protein FtsZ. The common strategy of protein adaptation from cold to hot environments implies the occurrence of small changes in the amino acid composition, without altering the overall structure of the macromolecule. These continuous adjustments were investigated through parameters related to the amino acid composition of each protein. The average value per residue of mass, volume and accessible surface area allowed an evaluation of the usage of bulky residues, whereas the average hydrophobicity reflected that of hydrophobic residues. The specific proportion of bulky and hydrophobic residues in each protein almost linearly increased with the temperature of the host microorganism. This finding agrees with the structural and functional properties exhibited by proteins in differently adapted sources, thus explaining the great compactness or the high flexibility exhibited by (hyper)thermophilic or psychrophilic proteins, respectively. Indeed, heat-adapted proteins incline toward the usage of heavier-size and more hydrophobic residues with respect to mesophiles, whereas the cold-adapted macromolecules show the opposite behavior with a certain preference for smaller-size and less hydrophobic residues. An investigation on the different increase of bulky residues along with the growth temperature observed in the six model proteins suggests the relevance of the possible different role and/or structure organization played by protein domains. The significance of the linear correlations between growth temperature and parameters related to the amino acid composition improved when the analysis was collectively carried out on all model proteins.     

Structural adaptation to low temperatures--analysis of the subunit interface of oligomeric psychrophilic enzymes

Enzymes from psychrophiles show higher catalytic efficiency in the 0-20 degrees C temperature range and often lower thermostability in comparison with meso/thermophilic homologs. Physical and chemical characterization of these enzymes is currently underway in order to understand the molecular basis of cold adaptation. Psychrophilic enzymes are often characterized by higher flexibility, which allows for better interaction with substrates, and by a lower activation energy requirement in comparison with meso/thermophilic counterparts. In their tertiary structure, psychrophilic enzymes present fewer stabilizing interactions, longer and more hydrophilic loops, higher glycine content, and lower proline and arginine content. In this study, a comparative analysis of the structural characteristics of the interfaces between oligomeric psychrophilic enzyme subunits was carried out. Crystallographic structures of oligomeric psychrophilic enzymes, and their meso/thermophilic homologs belonging to five different protein families, were retrieved from the Protein Data Bank. The following structural parameters were calculated: overall and core interface area, characterization of polar/apolar contributions to the interface, hydrophobic contact area, quantity of ion pairs and hydrogen bonds between monomers, internal area and total volume of non-solvent-exposed cavities at the interface, and average packing of interface residues. These properties were compared to those of meso/thermophilic enzymes. The results were analyzed using Student's t-test. The most significant differences between psychrophilic and mesophilic proteins were found in the number of ion pairs and hydrogen bonds, and in the apolarity of their subunit interface. Interestingly, the number of ion pairs at the interface shows an opposite adaptation to those occurring at the monomer core and surface.     

Factors determining the instability of the conformation of human mutant hemoglobins. A method for evaluating the instability

A method that allows to estimate the influence of single amino acid substitutions on the stability of proteins with known three-dimensional structure was developed. The role of 20 different physico-chemical properties of amino acids and some characteristics of the three-dimensional structure in the destabilization of protein structure caused by mutations was analyzed. Analysis of 54 proteins led us to reveal the characteristics of the three-dimensional structure. The parameters of the method were calculated on the basis of experimental data for 330 mutant human hemoglobins. The prediction accuracy of the method is on average 81%. The method may be used for the analysis of the total mutational spectra of the proteins and also for the estimation of effects of the single amino acid substitution on the protein structure in gene-engineering experiments.     

Amino acid similarity coefficients for protein modeling and sequence alignment derived from main-chain folding angles

A set of "similarity-parameters" was calculated that reflects the influence of the proteinogenic amino acids on the structure of the protein backbone. The parameters were derived from a detailed analysis of the amino acid specific main-chain torsion angle distributions as they are found in proteins (highly resolved protein structures from the Brookhaven Protein Data Bank). The purpose of these parameters is threefold: (1) they should help in estimating the structural effect of an amino acid substitution during the design of new mutants in protein-engineering; (2) in modeling by homology they should mark places in the protein where changes in the folding are expected; and (3) they should form a scoring matrix in protein sequence alignment superior to identity scoring. The usability of the "structure derived correlation matrix (SCM)" for these purposes is assessed and demonstrated for some examples in the paper.     

木聚糖酶与XIP型木聚糖酶抑制蛋白相互作用分子机制的研究进展

Xylanase inhibitor protein (XIP)型木聚糖酶抑制蛋白对大部分GH10、GH11家族的真菌木聚糖酶具有抑制作用,但是却不能抑制细菌来源和植物自身所产生的木聚糖酶。XIP型木聚糖酶抑制蛋白对木聚糖酶的抑制作用主要是通过模拟底物接触酶的活性位点,迅速阻塞底物进入活性位点区域的通道。然而,在对XIP型木聚糖酶抑制蛋白具有抗性的GH10和GH11木聚糖酶晶体结构中,连接二级结构的Loop构象严重阻碍了XIP型木聚糖酶抑制蛋白的抑制功能。与对XIP型木聚糖酶抑制蛋白敏感的木聚糖酶相比,氨基酸残基的插入突变导致抗性木聚糖酶的Loop具有明显的凸出构象;而在GH11家族抗性木聚糖酶中,"拇指"结构中部分氨基酸的替换致使XIP型木聚糖酶抑制蛋白与"拇指"结构无法形成稳固的氢键和疏水建,从而削弱XIP的抑制作用。     

Physicochemical evolution and molecular adaptation of the cetacean and artiodactyl cytochrome b proteins

Cetaceans have most likely experienced metabolic shifts since evolutionarily diverging from their terrestrial ancestors, shifts that may be reflected in the proteins such as cytochrome b that are responsible for metabolic efficiency. However, accepted statistical methods for detecting molecular adaptation are largely biased against even moderately conservative proteins because the primary criterion involves a comparison of nonsynonymous and synonymous substitution rates (dN/dS); they do not allow for the possibility that adaptation may come in the form of very few amino acid changes. We apply the MM01 model to the possible molecular adaptation of cytochrome b among cetaceans because it does not rely on a dN/dS ratio, instead evaluating positive selection in terms of the amino acid properties that comprise protein phenotypes that selection at the molecular level may act upon. We also apply the codon-degeneracy model (CDM), which focuses on evaluating overall patterns of nucleotide substitution in terms of base exchange, codon position, and synonymy to estimate the overall effect of selection. Using these relatively new models, we characterize the molecular adaptation that has occurred in the cetacean cytochrome b protein by comparing revealed amino acid replacement patterns to those found among artiodactyls, the modern terrestrial mammals found to be most closely related to cetaceans. Our findings suggest that several regions of the cetacean cytochrome b protein have experienced molecular adaptation. Also, these adaptations are spatially associated with domain structure, protein function, and the structure and function of the cytochrome bc(1) complex and its constituents. We also have found a general correlation between the results of the analytical software programs TreeSAAP (which implements the MM01 model) and CDM (which implements the codon-degeneracy model).     

Two types of amino acid substitutions in protein evolution

Summary The frequency of amino acid substitutions, relative to the frequency expected by chance, decreases linearly with the increase in physico-chemical differences between amino acid pairs involved in a substitution. This correlation does not apply to abnormal human hemoglobins. Since abnormal hemoglobins mostly reflect the process of mutation rather than selection, the correlation manifest during protein evolution between substitution frequency and physico-chemical difference in amino acids can be attributed to natural selection. Outside of abnormal proteins, the correlation also does not apply to certain regions of proteins characterized by rapid rates of substitution. In these cases again, except for the largest physico-chemical differences between amino acid pairs, the substitution frequencies seem to be independent of the physico-chemical parameters. The limination of the substituents involving the largest physicochemical differences can once more be attributed to natural selection. For smaller physico-chemical differences, natural selection, if it is operating in the polypeptide regions, must be based on parameters other than those examined.     

Crystal structure of a cold-adapted class C beta-lactamase

In this study, the crystal structure of a class C beta-lactamase from a psychrophilic organism, Pseudomonas fluorescens, has been refined to 2.2 A resolution. It is one of the few solved crystal structures of psychrophilic proteins. The structure was compared with those of homologous mesophilic enzymes and of another, modeled, psychrophilic protein. The elucidation of the 3D structure of this enzyme provides additional insights into the features involved in cold adaptation. Structure comparison of the psychrophilic and mesophilic beta-lactamases shows that electrostatics seems to play a major role in low-temperature adaptation, with a lower total number of ionic interactions for cold enzymes. The psychrophilic enzymes are also characterized by a decreased number of hydrogen bonds, a lower content of prolines, and a lower percentage of arginines in comparison with lysines. All these features make the structure more flexible so that the enzyme can behave as an efficient catalyst at low temperatures.     

Approaches to predicting effects of single amino acid substitutions on the function of a protein

The relative activities of 313 mutants of the gene V protein of bacteriophage f1, assayed in vivo, have been used to evaluate two approaches to predicting the effects of single amino acid substitutions on the function of a protein. First, we tested methods that only depend on the properties of the wild-type and substituting amino acids. None of the properties or measures of the functional equivalence of amino acids we tested, including the frequency of exchange of amino acids among homologous proteins as well as changes in side-chain size, hydrophobicity, and charge, were found to be more than weakly correlated with the activities of mutants. The principal reason for this poor correlation was found to be that the effect of a particular substitution varies considerably from site to site. We then tested an approach using the activities of several mutants with substitutions at a site to predict the activity of another mutant, and we find that this is a relatively good indicator of whether the other mutant at that site will be functional. A predictive scheme was developed that combines the weak information from the models depending on the properties of the wild-type and substituting amino acids with the stronger information from the tolerance of a site to substitution. Although this scheme requires no knowledge of the structure of a mutant protein, it is useful in predicting the activities of mutants.     

Substitution of aspartic acid with glutamic acid increases the unfolding transition temperature of a protein

Proteins from thermophiles are more stable than those from mesophiles. Several factors have been suggested as causes for this greater stability, but no general rule has been found. The amino acid composition of thermophile proteins indicates that the content of polar amino acids such as Asn, Gln, Ser, and Thr is lower, and that of charged amino acids such as Arg, Glu, and Lys is higher than in mesophile proteins. Among charged amino acids, however, the content of Asp is even lower in thermophile proteins than in mesophile proteins. To investigate the reasons for the lower occurrence of Asp compared to Glu in thermophile proteins, Glu was substituted with Asp in a hyperthermophile protein, MjTRX, and Asp was substituted with Glu in a mesophile protein, ETRX. Each substitution of Glu with Asp decreased the Tm of MjTRX by about 2 degrees C, while each substitution of Asp with Glu increased the Tm of ETRX by about 1.5 degrees C. The change of Tm destabilizes the MjTRX by 0.55 kcal/mol and stabilizes the ETRX by 0.45 kcal/mol in free energy.     

Nucleus-targeting domain of the matrix protein (M1) of influenza virus.

The matrix protein M1 of influenza virus A/WSN/33 was shown by immunofluorescent staining to be transported into the nuclei of transfected cells without requiring other viral proteins. We postulated the existence of a potential signal sequence at amino acids 101 to 105 (RKLKR) that is required for nuclear localization of the M1 protein. When CV1 cells were transfected with recombinant vectors expressing the entire M1 protein (amino acids 1 to 252) or just the first 112 N-terminal amino acids, both the complete M1 protein and the truncated M1 protein were transported to the nucleus. In contrast, expression in CV1 cells of vectors coding for M1 proteins with deletions from amino acids 77 to 202 or amino acids 1 to 134 resulted only in cytoplasmic immunofluorescent staining of these truncated M1 proteins without protein being transported to the nucleus. Moreover, no nuclear membrane translocation occurred when CV1 cells were transfected with recombinant vectors expressing M1 proteins with deletions of amino acids 101 to 105 or with substitution at amino acids 101 to 105 of SNLNS for RKLKR. Furthermore, a synthetic oligopeptide corresponding to M1 protein amino acids 90 to 108 was also transported into isolated nuclei derived from CV1 cells, whereas oligopeptides corresponding to amino acid sequences 25 to 40, 67 to 81, and 135 to 164 were not transported into the isolated cell nuclei. These data suggest that the amino acid sequence 101RKLKR105 is the nuclear localization signal of the M1 protein.     

V3 loop of the human immunodeficiency virus type 1 Env protein: interpreting sequence variability.

Two different states of human immunodeficiency virus type 1 are apparent in the asymptomatic and late stages of infection. Important determinants associated with these two states have been found within the V3 loop of the viral Env protein. In this study, two large data sets of published V3 sequences were analyzed to identify patterns of sequence variability that would correspond to these two states of the virus. We were especially interested in the pattern of basic amino acid substitutions, since the presence of basic amino acids in V3 has been shown to change virus tropism in cell culture. Four features of the sequence heterogeneity in V3 were observed: (i) approximately 70% of all nonconservative basic substitutions occur at four positions in V3, and V3 sequences with a basic substitution in at least one of these four positions contain approximately 95% of all nonconservative basic substitutions; (ii) substitution patterns within V3 are influenced by the identity of the amino acid at position 25; (iii) sequence polymorphisms account for a significant fraction of uncharged amino acid substitutions at several positions in V3, and sequence heterogeneity other than these polymorphisms is most significant at two positions near the tip of V3; and (iv) sequence heterogeneity in V3 (in addition to the basic amino acid substitutions) is approximately twofold greater in V3 sequences that contain basic amino acid substitutions. By using this sequence analysis, we were able to identify distinct groups of V3 sequences in infected patients that appear to correspond to these two virus states. The identification of these discrete sequence patterns in vivo demonstrates how the V3 sequence can be used as a genetic marker for studying the two states of human immunodeficiency virus type 1.     

Predicting isomorphic residue replacements for protein design

The replaceability of amino acids as reflected in their neighbourhood selectivity is analyzed in this paper. Neighbourhood selectivity was found to be more characteristic of the individual amino acids than the other commonly used parameters. The residue replacement rank obtained was compared with other proposed ranks and was tested on naturally accepted point mutations in homologous proteins.     

Mutation of amino acids in the alpha 1,3-fucosyltransferase motif affects enzyme activity and Km for donor and acceptor substrates

Alpha 1,3-fucosyltransferases (FucT) share a conserved amino acid sequence designated the alpha 1,3 FucT motif that has been proposed to be important for nucleotide sugar binding. To evaluate the importance of the amino acids in this motif, each of the alpha 1,3 FucT motif amino acids was replaced with alanine (alanine scanning mutagenesis) in human FucT VI, and the resulting mutant proteins were analyzed for enzyme activity and kinetically characterized in those cases in which the mutant protein had sufficient activity. Two of the mutant proteins were inactive, six had less than 1% of wild-type activity, and four had approximately 10-50% of wild-type enzyme activity. Three of the mutant proteins with significant enzyme activity had substantially larger Km (5 to 15 times) for GDP-fucose than FucT VI wild-type enzyme. The fourth mutant protein with significant enzyme activity (S249A) had a Km at least 10 times larger than wild-type FucT VI for the acceptor substrate, with only a slightly larger (2-3 times) Km for GDP-fucose. Thus mutation of any of the amino acids within the alpha 1,3 FucT motif to Ala affects alpha 1,3-FucT activity, and substitution of Ala for some of the alpha 1,3 FucT motif amino acids results in proteins with altered kinetic constants for both the acceptor and donor substrates. Secondary structure prediction suggests a helix-loop-helix fold for the alpha 1,3 FucT motif, which can be used to rationalize the effects of mutations in terms of 3D structure.     

A topological model for the haemolysin translocator protein HlyD

Summary A topological model for HlyD is proposed that is based on results obtained with gene fusions of lacZ and phoA to hlyD. Active H1yD-LacZ fusion proteins were only generated when lacZ was fused to hlyD. within the first 180 by (60 amino acids). H1yD-PhoA proteins exhibiting alkaline phosphatase (AP) activity were obtained when phoA was inserted into hlyD. between nucleotides 262 (behind amino acid position 87) and 1405 (behind amino acid position 468, only 10 amino acids away from the C-terminus of HlyD Active insertions of phoA into the middle region of hlyD. were not observed on in vivo transposition but such fusions exhibiting AP activity could be constructed by in vitro techniques. A fusion protein that carried the PhoA part close to the C-terminal end of HlyD proved to be the most stable HlyD-PhoA fusion protein. In contrast to the other, rather unstable, HlyD-PhoA⁺ fusions, no proteolytic degradation product of this HlyD-PhoA protein was observed and nearly all the alkaline phosphatase activity was membrane bound. Protease accessibility and cell fractionation experiments indicated that the alkaline phosphatase moiety of this fusion protein was located in the periplasm as for all other HlyD-PhoA⁺ proteins. These data and computer-assisted predictions suggest a topological model for HlyD with the N-terminal 60 amino acids located in the cytoplasm, a single transmembrane segment from amino acids 60 to 80 and a large periplasmic region extending from amino acid 80 to the C-terminus. Neither the HlyD fusion proteins obtained nor a mutant HlyD protein that had lost the last 10 amino acids from the C-terminus of HlyD exhibited translocator activity for HlyA or other reporter proteins carrying the HlyA signal sequence. The C-terminal 10 amino acids of HlyD showed significant similarity with the corresponding sequences of other HlyD-related proteins involved in protein secretion.     

Tertiary structural constraints on protein evolutionary diversity: templates, key residues and structure prediction

The pattern of residue substitution in divergently evolving families of globular proteins is highly variable. At each position in a fold there are constraints on the identities of amino acids from both the three-dimensional structure and the function of the protein. To characterize and quantify the structural constraints, we have made a comparative analysis of families of homologous globular proteins. Residues are classified according to amino acid type, secondary structure, accessibility of the sidechain, and existence of hydrogen bonds from sidechain to other sidechains or peptide carbonyl or amide functions. There are distinct patterns of substitution especially where residues are both solvent inaccessible and hydrogen bonded through their sidechains. The patterns of residue substitution can be used to construct templates or to identify 'key' residues if one or more structures are known. Conversely, analysis of conversation and substitution across a large family of aligned sequences in terms of substitution profiles can allow prediction of tertiary environment or indicate a functional role. Similar analyses can be used to test the validity of putative structures if several homologous sequences are available.     

RUBISCO ADAPTATION TO LOW TEMPERATURES: A COMPARATIVE STUDY IN PSYCHROPHILIC AND MESOPHILIC UNICELLULAR ALGAE

Some properties of the ribulose-1,5-bisphosphate carboxylase/oxygenase (RUBISCO) from two psychrophilic Chloromonas species have been investigated in relation to their adaptation to cold environments. Contrary to the situation usually encountered with psychrophilic enzymes, the carboxylase activity of both purified "cold" RUBISCO enzymes was lower at low temperatures than that found with the enzyme of the mesophilic alga Chlamydomonas reinhardtii Dangeard. Moreover, the apparent optimal temperature for RUBISCO carboxylase activity was similar for psychrophilic and mesophilic enzymes. Psychrophilic RUBISCOs, however, showed a greater thermosensitivity than the C. reinhardtii enzyme. Genes encoding small and large subunits of RUBISCO from one psychrophilic isolate were sequenced. Comparison of the deduced amino acid sequences to those of higher plants and green algae revealed the substitution of a very highly conserved residue (cysteine₂₄₇ → serine in the large subunit) that could be responsible, at least in part, for the increased thermosensitivity of the "cold" enzyme. Interestingly, the relative amount of RUBISCO subunits found in the psychrophilic isolates was about twice as high as the amount observed in C. reinhardtii and five other mesophilic algae. The high production of a key enzyme to counterbalance its poor catalytic efficiency at low temperature could constitute a novel type of adaptive mechanism to cold environments.