Protein flexibility in psychrophilic and mesophilic trypsins. Evidence of evolutionary conservation of protein dynamics in trypsin-like serine-proteases

To shed light on the molecular features related to cold-adaptation in serine-proteases, we have carried out molecular dynamics simulations of homologous mesophilic and psychrophilic trypsins, with particular attention to evaluation of intramolecular interactions and flexibility. Psychrophilic trypsins present fewer interdomain interactions and enhanced localized flexibility in regions close to the catalytic site. Notably, these regions fit well with the pattern of protein flexibility previously reported for psychrophilic elastases. Our results indicate that specific sites within the serine-protease fold can be considered hot spots of cold-adaptation and that psychrophilic trypsins and elastases have independently discovered similar molecular strategies to optimize flexibility at low temperatures.     

Exploring Local Flexibility/Rigidity in Psychrophilic and Mesophilic Carbonic Anhydrases

Molecular flexibility and rigidity are required to determine the function and specificity of protein molecules. Some psychrophilic enzymes demonstrate a higher catalytic efficiency at low temperatures, compared to the efficiency demonstrated by their meso/thermophilic homologous. The emerging picture suggests that such enzymes have an improved flexibility of the structural catalytic components, whereas other protein regions far from functional sites may be even more rigid than those of their mesophilic counterparts. To gain a deeper insight in the analysis of the activity-flexibility/rigidity relationship in protein structure, psychrophilic carbonic anhydrase of the Antarctic teleost Chionodraco hamatus has been compared with carbonic anhydrase II of Bos taurus through fluorescence studies, three-dimensional modeling, and activity analyses. Data demonstrated that the cold-adapted enzyme exhibits an increased catalytic efficiency at low and moderate temperatures and, more interestingly, a local flexibility in the region that controls the correct folding of the catalytic architecture, as well as a rigidity in the hydrophobic core. The opposite result was observed in the mesophilic counterpart. These results suggest a clear relationship between the activity and the presence of flexible and rigid protein substructures that may be useful in rational molecular and drug design of a class of enzymes playing a key role in pathologic processes.     

Structural and functional adaptations to extreme temperatures in psychrophilic,mesophilic, and thermophilic DNA ligases

Psychrophiles, host of permanently cold habitats, display metabolic fluxes comparable to those exhibited by mesophilic organisms at moderate temperatures. These organisms have evolved by producing, among other peculiarities, cold-active enzymes that have the properties to cope with the reduction of chemical reaction rates induced by low temperatures. The emerging picture suggests that these enzymes display a high catalytic efficiency at low temperatures through an improved flexibility of the structural components involved in the catalytic cycle, whereas other protein regions, if not implicated in catalysis, may be even more rigid than their mesophilic counterparts. In return, the increased flexibility leads to a decreased stability of psychrophilic enzymes. In order to gain further advances in the analysis of the activity/flexibility/stability concept, psychrophilic, mesophilic, and thermophilic DNA ligases have been compared by three-dimensional-modeling studies, as well as regards their activity, surface hydrophobicity, structural permeability, conformational stabilities, and irreversible thermal unfolding. These data show that the cold-adapted DNA ligase is characterized by an increased activity at low and moderate temperatures, an overall destabilization of the molecular edifice, especially at the active site, and a high conformational flexibility. The opposite trend is observed in the mesophilic and thermophilic counterparts, the latter being characterized by a reduced low temperature activity, high stability and reduced flexibility. These results strongly suggest a complex relationship between activity, flexibility and stability. In addition, they also indicate that in cold-adapted enzymes, the driving force for denaturation is a large entropy change.     

Dynamic properties of extremophilic subtilisin-like serine-proteases

The investigation of the structural determinants of enzymatic temperature adaptation is a crucial pre-requisite both in terms of fundamental research and industrial applications to develop new biocatalysts active at different temperature ranges. In several cases, the differences related to cold- or warm-adaptation are related to subtle structural and aminoacidic differences at the molecular level, often hard to detect. In this context, we present a comparative study of psychrophilic, mesophilic and thermophilic subtilisin-like serine proteases by all-atom molecular dynamics (MD) simulations in explicit solvent using a multiple-replica approach. Our results strongly enforce the current view on localized flexibility in crucial functional regions for cold-adapted serine proteases and point out a different optimization and usage of salt-bridge interactions and networks in cold- and warm-adapted enzymes. The analyses allow to identify a subset of structural and dynamic features strictly associated to cold adaptation and which change from cold- to heat-active subtilisins. In particular, the thermophilic subtilisin presents a high affinity calcium binding site which is not structurally conserved in the mesophilic and psychrophilic counterparts, which, as it turns out from the MD analyses, at the same position show a stable salt bridge network and no stabilizing intra-molecular interactions, respectively. These aspects, along with differential flexibility in regions close to the active site or substrate binding pocket, can be an indication of evolution at this protein site toward a lower stability moving from high to low temperature conditions.     

"Cold spots" in protein cold adaptation: Insights from normalized atomic displacement parameters (B'-factors)

Many analyses published in the last decade suggest that enzymes isolated from cold-adapted organisms are characterized by a higher flexibility of their molecular structure. Recently, it has been argued that all cold-adapted enzymes with catalytic efficiency greater than that of their mesophilic counterparts display local flexibility or rigidity that are likely to cooperate, each acting on specific areas of the enzyme structure. Here we report an analysis of the normalized thermal B-factor distributions in psychrophilic proteins compared with those of their mesophilic and thermophilic counterparts with the aim to detect statistically significant local variations of relative backbone flexibility possibly linked to cold adaptation. We utilized a strategy based mainly on intra-family comparison of local distribution of normalized B-factors. After careful statistical treatment of data, the picture emerging from our results suggests that the distribution of the flexibility in psychrophilic enzymes is locally more heterogeneous than in their respective mesophilic homologues.     

Some like it cold: biocatalysis at low temperatures

In the last few years, increased attention has been focused on a class of organisms called psychrophiles. These organisms, hosts of permanently cold habitats, often display metabolic fluxes more or less comparable to those exhibited by mesophilic organisms at moderate temperatures. Psychrophiles have evolved by producing, among other peculiarities, "cold-adapted" enzymes which have the properties to cope with the reduction of chemical reaction rates induced by low temperatures. Thermal compensation in these enzymes is reached, in most cases, through a high catalytic efficiency associated, however, with a low thermal stability. Thanks to recent advances provided by X-ray crystallography, structure modelling, protein engineering and biophysical studies, the adaptation strategies are beginning to be understood. The emerging picture suggests that psychrophilic enzymes are characterized by an improved flexibility of the structural components involved in the catalytic cycle, whereas other protein regions, if not implicated in catalysis, may be even more rigid than their mesophilic counterparts. Due to their attractive properties, i.e., a high specific activity and a low thermal stability, these enzymes constitute a tremendous potential for fundamental research and biotechnological applications.     

Cold adaptation of enzyme reaction rates

A major issue for organisms living at extreme temperatures is to preserve both stability and activity of their enzymes. Cold-adapted enzymes generally have a reduced thermal stability, to counteract freezing, and show a lower enthalpy and a more negative entropy of activation compared to mesophilic and thermophilic homologues. Such a balance of thermodynamic activation parameters can make the reaction rate decrease more linearly, rather than exponentially, as the temperature is lowered, but the structural basis for rate optimization toward low working temperatures remains unclear. In order to computationally address this problem, it is clear that reaction simulations rather than standard molecular dynamics calculations are needed. We have thus carried out extensive computer simulations of the keto-enol(ate) isomerization steps in differently adapted citrate synthases to explore the structure-function relationships behind catalytic rate adaptation to different temperatures. The calculations reproduce the absolute rates of the psychrophilic and mesophilic enzymes at 300 K, as well as the lower enthalpy and more negative entropy of activation of the cold-adapted enzyme, where the latter simulation result is obtained from high-precision Arrhenius plots. The overall catalytic effect originates from electrostatic stabilization of the transition state and enolate and the reduction of reorganization free energy. The simulations, however, show psychrophilic, mesophilic, and hyperthermophilic citrate synthases to have increasingly stronger electrostatic stabilization of the transition state, while the energetic penalty in terms of internal protein interactions follows the reverse order with the cold-adapted enzyme having the most favorable energy term. The lower activation enthalpy and more negative activation entropy observed for cold-adapted enzymes are found to be associated with a decreased protein stiffness. The origin of this effect is, however, not localized to the active site but to other regions of the protein structure.     

Comparative thermal unfolding study of psychrophilic and mesophilic subtilisin-like serine proteases by molecular dynamics simulations

Molecular dynamics (MD) simulations of a subtilisin-like serine protease VPR from the psychrophilic marine bacterium Vibrio sp. PA-44 and its mesophilic homologue, proteinase K (PRK), have been performed for 20 ns at four different temperatures (300, 373, 473, and 573 K). The comparative analyses of MD trajectories reveal that at almost all temperatures, VPR exhibits greater structural fluctuations/deviations, more unstable regular secondary structural elements, and higher global flexibility than PRK. Although these two proteases follow similar unfolding pathways at high temperatures, VPR initiates unfolding at a lower temperature and unfolds faster at the same high temperatures than PRK. These observations collectively indicate that VPR is less stable and more heat-labile than PRK. Analyses of the structural/geometrical properties reveal that, when compared to PRK, VPR has larger radius of gyration (Rg), less intramolecular contacts and hydrogen bonds (HBs), more protein-solvent HBs, and smaller burial of nonpolar area and larger exposure of polar area. These suggest that the increased flexibility of VPR would be most likely caused by its reduced intramolecular interactions and more favourable protein-solvent interactions arising from the larger exposure of the polar area, whereas the enhanced stability of PRK could be ascribed to its increased intramolecular interactions arising from the better optimized hydrophobicity. The factors responsible for the significant differences in local flexibility between these two proteases were also analyzed and ascertained. This study provides insights into molecular basis of thermostability of homologous serine proteases adapted to different temperatures.     

Cold adaptation of enzymes: structural, kinetic and microcalorimetric characterizations of an aminopeptidase from the Arctic psychrophile Colwellia psychrerythraea and of human leukotriene A(4) hydrolase

The relationships between structure, activity, stability and flexibility of a cold-adapted aminopeptidase produced by a psychrophilic marine bacterium have been investigated in comparison with a mesophilic structural and functional human homolog. Differential scanning calorimetry, fluorescence monitoring of thermal- and guanidine hydrochloride-induced unfolding and fluorescence quenching were used to show that the cold-adapted enzyme is characterized by a high activity at low temperatures, a low structural stability versus thermal and chemical denaturants and a greater structural permeability to a quenching agent relative to the mesophilic homolog. These findings support the hypothesis that cold-adapted enzymes maintain their activity at low temperatures as a result of increased global or local structural flexibility, which results in low stability. Analysis of the thermodynamic parameters of irreversible thermal unfolding suggests that entropy-driven factors are responsible for the fast unfolding rate of the cold-adapted aminopeptidase. A reduced number of proline residues, a lower degree of hydrophobic residue burial and a decreased surface accessibility of charged residues may be responsible for this effect. On the other hand, the reduction in enthalpy-driven interactions is the primary determinant of the weak conformational stability.     

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.     

Stepwise adaptations to low temperature as revealed by multiple mutants of psychrophilic α-amylase from Antarctic Bacterium

The mutants Mut5 and Mut5CC from a psychrophilic α-amylase bear representative stabilizing interactions found in the heat-stable porcine pancreatic α-amylase but lacking in the cold-active enzyme from an Antarctic bacterium. From an evolutionary perspective, these mutants can be regarded as structural intermediates between the psychrophilic and the mesophilic enzymes. We found that these engineered interactions improve all the investigated parameters related to protein stability as follows: compactness; kinetically driven stability; thermodynamic stability; resistance toward chemical denaturation, and the kinetics of unfolding/refolding. Concomitantly to this improved stability, both mutants have lost the kinetic optimization to low temperature activity displayed by the parent psychrophilic enzyme. These results provide strong experimental support to the hypothesis assuming that the disappearance of stabilizing interactions in psychrophilic enzymes increases the amplitude of concerted motions required by catalysis and the dynamics of active site residues at low temperature, leading to a higher activity.     

Comparative void-volume analysis of psychrophilic and mesophilic enzymes: Structural bioinformatics of psychrophilic enzymes reveals sources of core flexibility

Background

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.

Results

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 1.4-1.5 Å, sufficient for water molecules. Furthermore, amino acid side chains lining these cavities show an increased frequency of acidic groups in psychrophilic enzymes.

Conclusions

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.

     

Study and design of stability in GH5 cellulases

Thermostable enzymes that hydrolyze lignocellulosic materials provide potential advantages in process configuration and enhancement of production efficiency over their mesophilic counterparts in the bioethanol industry. In this study, the dynamics of β-1,4-endoglucanases (EC: 3.2.1.4) from family 5 of glycoside hydrolases (GH5) were investigated computationally. The conformational flexibility of 12 GH5 cellulases, ranging from psychrophilic to hyperthermophilic, was investigated by molecular dynamics (MD) simulations at elevated temperatures. The results indicated that the protein flexibility and optimum activity temperatures are appreciably correlated. Intra-protein interactions, packing density and solvent accessible area were further examined in crystal structures to investigate factors that are possibly involved in higher rigidity of thermostable cellulases. The MD simulations and the rules learned from analyses of stabilizing factors were used in design of mutations toward the thermostabilization of cellulase C, one of the GH5 endoglucanases. This enzyme was successfully stabilized both chemically and thermally by introduction of a new disulfide cross-link to its highly mobile 56-amino acid subdomain.     

Structural, kinetic, and calorimetric characterization of the cold-active phosphoglycerate kinase from the antarctic Pseudomonas sp. TACII18

The gene encoding the phosphoglycerate kinase (PGK) from the Antarctic Pseudomonas sp. TACII18 has been cloned and found to be inserted between the genes encoding for glyceraldhyde-3-phosphate dehydrogenase and fructose aldolase. The His-tagged and the native recombinant PGK from the psychrophilic Pseudomonas were expressed in Escherichia coli. The wild-type and the native recombinant enzymes displayed identical properties, such as a decreased thermostability and a 2-fold higher catalytic efficiency at 25 degrees C when compared with the mesophilic PGK from yeast. These properties, which reflect typical features of cold-adapted enzymes, were strongly altered in the His-tagged recombinant PGK. The structural model of the psychrophilic PGK indicated that a key determinant of its low stability is the reduced number of salt bridges, surface charges, and aromatic interactions when compared with mesophilic and thermophilic PGK. Differential scanning calorimetry of the psychrophilic PGK revealed unusual variations in its conformational stability for the free and substrate-bound forms. In the free form, a heat-labile and a thermostable domain unfold independently. It is proposed that the heat-labile domain acts as a destabilizing domain, providing the required flexibility around the active site for catalysis at low temperatures.     

Comparative molecular dynamics of mesophilic and psychrophilic protein homologues studied by 1.2 ns simulations

It is well established that the dynamic motion of proteins plays an important functional role, and that the adaptation of a protein molecule to its environment requires optimization of internal non-covalent interactions and protein-solvent interactions. Serine proteinases in general, and trypsin in particular has been used as a model system in exploring possible structural features for cold adaptation. In this study, a 500 p.s. and a 1200 p.s. molecular dynamics (MD) simulation at 300 K of both anionic salmon trypsin and cationic bovine trypsin are analyzed in terms of molecular flexibility, internal non-covalent interactions and protein-solvent interactions. The present MD simulations do not indicate any increased flexibility of the cold adapted enzyme on an overall basis. However, the apparent higher flexibility and deformability of the active site of anionic salmon trypsin may lower the activation energy for ligand binding and for catalysis, and might be a reason for the increased binding affinity and catalytic efficiency compared to cationic bovine trypsin.     

Enzymes Produced by a Pseudomonas Species Which Inactivate Inhibitors of Certain Viral Hemagglutinins. I. Identification and Purification of a Proteinase and Phospholipase C

Filtrates from cultures of a psychrophilic Pseudomonas species, which inactivate serum inhibitors of certain viral hemagglutinins, were shown to contain both lecithinase (phospholipase C) and a proteolytic enzyme with elastase activity. The bacterium was cultivated under conditions favoring production of the respective enzymes, and the enzymes were purified by ammonium sulfate precipitation followed by column chromatography or by gel filtration. The elastase was obtained in crystalline form and was recrystallized. It has properties similar to those of a number of other bacterial elastases but is more heat-labile than most. Although a high degree of purification was achieved for the lecithinase, as evidenced by an increase in specific activity, it was not obtained in crystalline form. Partially purified preparations of the lecithinase had extremely high activity compared to that of commercial preparations of phospholipase C from Clostridium welchii.     

Structure of phenylalanine hydroxylase from Colwellia psychrerythraea 34H, a monomeric cold active enzyme with local flexibility around the active site and high overall stability

The characteristic of cold-adapted enzymes, high catalytic efficiency at low temperatures, is often associated with low thermostability and high flexibility. In this context, we analyzed the catalytic properties and solved the crystal structure of phenylalanine hydroxylase from the psychrophilic bacterium Colwellia psychrerythraea 34H (CpPAH). CpPAH displays highest activity with tetrahydrobiopterin (BH(4)) as cofactor and at 25 degrees C (15 degrees C above the optimal growth temperature). Although the enzyme is monomeric with a single L-Phe-binding site, the substrate binds cooperatively. In comparison with PAH from mesophilic bacteria and mammalian organisms, CpPAH shows elevated [S(0.5)](L-Phe) (= 1.1 +/- 0.1 mm) and K(m)(BH(4))(= 0.3 +/- 0.1 mm), as well as high catalytic efficiency at 10 degrees C. However, the half-inactivation and denaturation temperature is only slightly lowered (T(m) approximately 52 degrees C; where T(m) is half-denaturation temperature), in contrast to other cold-adapted enzymes. The crystal structure shows regions of local flexibility close to the highly solvent accessible binding sites for BH(4) (Gly(87)/Phe(88)/Gly(89)) and l-Phe (Tyr(114)-Pro(118)). Normal mode and COREX analysis also detect these and other areas with high flexibility. Greater mobility around the active site and disrupted hydrogen bonding abilities for the cofactor appear to represent cold-adaptive properties that do not markedly affect the thermostability of CpPAH.     

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.     

Specificity, stability, and potency of monocyclic beta-lactam inhibitors of human leucocyte elastase.

Stable, potent, highly specific, time-dependent monocyclic beta-lactam inhibitors of human leucocyte elastase (HLE) are described. The heavily substituted beta-lactams are stable under physiological conditions including in the presence of enzymes of the digestive tract. The beta-lactams were unstable in base. At pH 11.3 and 37 degrees C they were hydrolyzed with half-lives of 1.5-2 h. Hydrolysis produced characteristic products including the substituent originally at C-4 of the lactam ring, a substituted urea, and products resulting from decarboxylation of the acid after ring opening. The most potent beta-lactam displayed only 2-fold less activity versus HLE than alpha 1PI, the natural proteinaceous inhibitor. The compounds were more potent against the human and primate PMN elastases than versus either the dog or rat enzymes. Differences in the structure-activity relationships of the human versus the rat enzymes suggest significant differences between these two functionally similar enzymes. The specificity of these compounds toward HLE versus porcine pancreatic elastase (PPE) is consistent with the differences in substrate specificity reported for these enzymes [Zimmerman & Ashe (1977) Biochim. Biophys. Acta 480, 241-245]. These differences suggest that the alkyl substitutions at C-3 of the lactam ring bind in the S1 specificity pocket of these enzymes. The dependence of the stereochemistry at C-4 suggests additional differences between HLE and PPE. Most of the compounds do not inhibit other esterases or human proteases. Weak, time-dependent inhibition of human cathepsin G and alpha-chymotrypsin by one compound suggested a binding mode to these enzymes that places the N-1 substitution in the S1 pocket.     

Temperature adaptation analysis of a psychrophilic mannanase through structural,functional and molecular dynamics simulation

The present paper reports structure prediction and analysis of a psychrophilic β-mannanase from Glaciozyma antarctica PI12 yeast. A threading method was used for 3D structure prediction of the enzyme using the MODELLER 9v12 program regarding its low sequence identity (<30%). The constructed model has been used in a comparative study to analyse its cold adaptation mechanism using other mesophilic, thermophilic, and hyperthermophilic mannanases. The structural and molecular dynamics analysis suggests that flexibility of the enzyme is increased through different structural characteristics, and therefore, the possibility of efficient catalytic reactions is provided at cold environment. These characteristics are the presence of longer loops, broken or shorter strands and helices, a lower number of salt bridges and hydrogen bonds, a higher exposure of the hydrophobic side chains to the solvent and an increased total solvent accessible surface area. Furthermore, the high catalytic efficiency and structural flexibility of the psychrophilic mannanase was supported by the results of principal component analysis.