Crystal structure of a subtilisin-like serine proteinase from a psychrotrophic Vibrio species reveals structural aspects of cold adaptation

Exploring the molecular mechanism of cold-adaption of an alkaline protease mutant by molecular dynamics simulations and residue interaction network

Psychrophilic proteases have attracted enormous attention in past decades, due to their high catalytic activity at low temperatures in a wide range of industrial processes, especially in the detergent and leather industries. Among them, H5 is an alkaline protease mutant, which featuring psychrophilic-like behavior, but the reasons that H5 with higher activity at low temperatures are still poorly understood. Herein, the molecular dynamics (MD) simulations combined with residue interaction network (RIN) were utilized to investigate the mechanisms of the cold-adaption of mutant H5. The results demonstrated that two loops involved in the substrate binding G100-S104 and S125-S129 in H5 had higher mobility, and the distance enlargement between the two loops modulated the substrate's accessibility compared with wild type counterpart. Besides, H5 enhanced conformational flexibility by weakening salt bridges and increasing interaction with the solvent. In particular, the absence of Lys251-Asp197-Arg247 salt bridge network may contribute to the structural mobility. Based on the free energy landscape and molecular mechanics Poisson-Boltzmann surface area of the wild type and H5, it was elucidated that H5 possessed a large population of interconvertible conformations, resulting in the weaker substrate binding free energy. The calculated RIN topology parameters such as the average degree, average cluster coefficient, and average path length further verified that the mutant H5 attenuated residue-to-residue interactions. The investigation of the mechanisms by which how the residue mutation affects the stability and activity of enzymes provides a theoretical basis for the development of cold-adapted protease.

Cold-Adapted Proteases: An Efficient and Energy-Saving Biocatalyst

The modern biotechnology industry has a demand for macromolecules that can function in extreme environments. One example is cold-adapted proteases, possessing advantages such as maintaining high catalytic efficiency at low temperature and low energy input during production and inactivation. Meanwhile, cold-adapted proteases are characterised by sustainability, environmental protection, and energy conservation; therefore, they hold significant economic and ecological value regarding resource utilisation and the global biogeochemical cycle. Recently, the development and application of cold-adapted proteases have gained gaining increasing attention; however, their applications potential has not yet been fully developed, which has seriously restricted the promotion and application of cold-adapted proteases in the industry. This article introduces the source, related enzymology characteristics, cold resistance mechanism, and the structure-function relationship of cold-adapted proteases in detail. This is in addition to discussing related biotechnologies to improve stability, emphasise application potential in clinical medical research, and the constraints of the further developing of cold-adapted proteases. This article provides a reference for future research and the development of cold-adapted proteases.

Discrimination of psychrophilic enzymes using machine learning algorithms with amino acid composition descriptor

Introduction

Psychrophilic enzymes are a class of macromolecules with high catalytic activity at low temperatures. Cold-active enzymes possessing eco-friendly and cost-effective properties, are of huge potential application in detergent, textiles, environmental remediation, pharmaceutical as well as food industry. Compared with the time-consuming and labor-intensive experiments, computational modeling especially the machine learning (ML) algorithm is a high-throughput screening tool to identify psychrophilic enzymes efficiently.

Methods

In this study, the influence of 4 ML methods (support vector machines, K-nearest neighbor, random forest, and naïve Bayes), and three descriptors, i.e., amino acid composition (AAC), dipeptide combinations (DPC), and AAC + DPC on the model performance were systematically analyzed.

Results and discussion

Among the 4 ML methods, the support vector machine model based on the AAC descriptor using 5-fold cross-validation achieved the best prediction accuracy with 80.6%. The AAC outperformed than the DPC and AAC + DPC descriptors regardless of the ML methods used. In addition, amino acid frequencies between psychrophilic and non-psychrophilic proteins revealed that higher frequencies of Ala, Gly, Ser, and Thr, and lower frequencies of Glu, Lys, Arg, Ile,Val, and Leu could be related to the protein psychrophilicity. Further, ternary models were also developed that could classify psychrophilic, mesophilic, and thermophilic proteins effectively. The predictive accuracy of the ternary classification model using AAC descriptor via the support vector machine algorithm was 75.8%. These findings would enhance our insight into the cold-adaption mechanisms of psychrophilic proteins and aid in the design of engineered cold-active enzymes. Moreover, the proposed model could be used as a screening tool to identify novel cold-adapted proteins.

Identification of Exoenzymes Secreted by Entomopathogenic Fungus Beauveria pseudobassiana RGM 2184 and Their Effect on the Degradation of Cocoons and Pupae of Quarantine Pest Lobesia botrana

Beauveria pseudobassiana RGM 2184 has shown 80% maximum efficacy against the pest Lobesia botrana in the autumn and winter seasons. This suggests that the strain possesses an interesting battery of enzymes that are cold-adapted to penetrate the thick and hydrophobic cocoon of L. botrana. In this study, screening of the proteolytic, lipolytic, and chitinolytic activity of enzyme extracts secreted by the RGM 2184 strain was carried out in various culture media. The enzyme extracts with the highest activity were subjected to zymography and mass spectrometry. These analyses allowed the identification of two proteases, two lipases, and three chitinases. Comparative analysis indicated that the degree of similarity between these enzymes was substantially reduced when the highest degree of taxonomic relatedness between RGM 2184 and the entomopathogenic fungus strain was at the family level. These results suggest that there is a wide variety of exoenzymes in entomopathogenic fungi species belonging to the order Hypocreales. On the other hand, exoenzyme extract exposure of cocoons and pupae of L. botrana provoked damage at 10 °C. Additionally, an analysis of the amino acid composition of the RGM 2184 exoenzyme grouped them close to the cold-adapted protein cluster. These results support the use of this strain to control pests in autumn and winter. Additionally, these antecedents can form a scaffold for the future characterization of these exoenzymes along with the optimization of the strain’s biocontrol ability by overexpressing them.

Thermostabilization of VPR, a kinetically stable cold adapted subtilase, via multiple proline substitutions into surface loops

Protein stability is a widely studied topic, there are still aspects however that need addressing. In this paper we examined the effects of multiple proline substitutions into loop regions of the kinetically stable proteinase K-like serine protease VPR, using the thermostable structural homologue AQUI as a template. Four locations for proline substitutions were chosen to imitate the structure of AQUI. Variants were produced and characterized using differential scanning calorimetry (DSC), circular dichroism (CD), steady state fluorescence, acrylamide fluorescence quenching and thermal inactivation experiments. The final product VPR_ΔC_N3P/I5P/N238P/T265P was greatly stabilized which was achieved without any noticeable detrimental effects to the catalytic efficiency of the enzyme. This stabilization seems to be derived from the conformation restrictive properties of the proline residue in its ability to act as an anchor point and strengthen pre-existing interactions within the protein and allowing for these interactions to prevail when thermal energy is applied to the system. In addition, the results underline the importance of the synergy between distant local protein motions needed to result in stabilizing effects and thus giving an insight into the nature of the stability of VPR, its unfolding landscape and how proline residues can infer kinetic stability onto protein structures.

Improved expression, purification and characterization of VPR, a cold active subtilisin-like serine proteinase and the effects of calcium on expression and stability

Cloning into a pET 11a vector, followed by high-level expression of the cold adapted subtilase, VPR, utilizing the rhamnose titratable T7 system of Lemo21, resulted in a dramatic increase of soluble protein compared to the older system used. Expression optimization clearly shows the importance of calcium in the medium after induction, both for stability of the proteinase and cell health. Characterization of the purified enzyme obtained in a redesigned purification protocol which removed apparent RNA contaminants, resulted in a significantly higher value for k_cat than previously reported. The new recombinant protein exhibited slightly lower stability against thermal denaturation and thermal inactivation. Our results also indicate that two of the calcium binding sites have apparent binding constants in the mM range. Binding of calcium to the weaker of those two sites only affects resistance of the enzyme against irreversible thermal inactivation. Differential scanning calorimetry revealed a non-two-state denaturation process, with indication of presence of intermediates caused by unfolding of calcium binding motifs.

Comparative evolutionary histories of fungal proteases reveal gene gains in the mycoparasitic and nematode-parasitic fungus Clonostachys rosea

Background

The ascomycete fungus Clonostachys rosea (order Hypocreales) can control several important plant diseases caused by plant pathogenic fungi and nematodes. Subtilisin-like serine proteases are considered to play an important role in pathogenesis in entomopathogenic, mycoparasitic, and nematophagous fungi used for biological control. In this study, we analysed the evolutionary histories of protease gene families, and investigated sequence divergence and regulation of serine protease genes in C. rosea.

Results

Proteases of selected hypocrealean fungal species were classified into families based on the MEROPS peptidase database. The highest number of protease genes (590) was found in Fusarium solani, followed by C. rosea with 576 genes. Analysis of gene family evolution identified non-random changes in gene copy numbers in the five serine protease gene families S1A, S8A, S9X, S12 and S33. Four families, S1A, S8A, S9X, and S33, displayed gene gains in C. rosea. A gene-tree / species-tree reconciliation analysis of the S8A family revealed that the gene copy number increase in C. rosea was primarily associated with the S08.054 (proteinase K) subgroup. In addition, regulatory and predicted structural differences, including twelve sites evolving under positive selection, among eighteen C. rosea S8A serine protease paralog genes were also observed. The C. rosea S8A serine protease gene prs6 was induced during interaction with the plant pathogenic species F. graminearum.

Conclusions

Non-random increases in S8A, S9X and S33 serine protease gene numbers in the mycoparasitic species C. rosea, Trichoderma atroviride and T. virens suggests an involvement in fungal-fungal interactions. Regulatory and predicted structural differences between C. rosea S8A paralogs indicate that functional diversification is driving the observed increase in gene copy numbers. The induction of prs6 expression in C. rosea during confrontation with F. graminearum suggests an involvement of the corresponding protease in fungal-fungal interactions. The results pinpoint the importance of serine proteases for ecological niche adaptation in C. rosea, including a potential role in the mycoparasitic attack on fungal prey.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1291-1) contains supplementary material, which is available to authorized users.

Local structural motifs in proteins: Detection and characterization of fragments inserted in helices

The global/local fold of protein structures is stabilized by a variety of specific interactions. A primary role in this context is played by hydrogen bonds. In order to identify novel motifs in proteins, we searched Protein Data Bank structures looking for backbone H-bonds formed by NH groups of two (or more) consecutive residues with consecutive CO groups of distant residues in the sequence. The present analysis unravels the occurrence of recurrent structural motifs that, to the best of our knowledge, had not been characterized in literature. Indeed, these H-bonding patterns are found (i) in a specific parallel β-sheet capping, (ii) in linking of β-hairpins to α-helices, and (iii) in α-helix insertions. Interestingly, structural analyses of these motifs indicate that Gly residues frequently occupy prominent positions. The formation of these motifs is likely favored by the limited propensity of Gly to be embodied in helices/sheets. Of particular interest is the motif corresponding to insertions in helices that was detected in 1% of analyzed structures. Inserted fragments may assume different structures and aminoacid compositions and usually display diversified evolutionary conservation. Since inserted regions are physically separated from the rest of the protein structure, they represent hot spots for ad-hoc protein functionalization.

Molecular attributes and apoptosis-inducing activities of a putative serine protease isolated from Tiger Milk mushroom (Lignosus rhinocerus) sclerotium against breast cancer cells in vitro

Background

The highly valued medicinal tiger milk mushroom (also known as Lignosus rhinocerus) has the ability to cure numerous ailments. Its anticancer activities are well explored, and recently a partially purified cytotoxic protein fraction termed F5 from the mushroom’s sclerotial cold water extract consisting mainly of fungal serine proteases was found to exhibit potent selective cytotoxicity against a human breast adenocarcinoma cell line (MCF7) with IC₅₀ value of 3.00 μg/ml. However, characterization of its cell death-inducing activity has yet to be established.

Methods

The mechanism involved in the cytotoxic activities of F5 against MCF7 cells was elucidated by flow cytometry-based apoptosis detection, caspases activity measurement, and expression profiling of apoptosis markers by western blotting. Molecular attributes of F5 were further mined from L. rhinocerus’s published genome and transcriptome for future exploration.

Results and Discussion

Apoptosis induction in MCF7 cells by F5 may involve a cross-talk between the extrinsic and intrinsic apoptotic pathways with upregulation of caspase-8 and -9 activities and a marked decrease of Bcl-2. On the other hand, the levels of pro-apoptotic Bax, BID, and cleaved BID were increased accompanied by observable actin cleavage. At gene level, F5 composed of three predicted non-synonymous single nucleotide polymorphisms (T > C) and an alternative 5′ splice site.

Conclusions

Findings from this study provide an advanced framework for further investigations on cancer therapeutics development from L. rhinocerus.

Insights into the role of electrostatics in temperature adaptation: a comparative study of psychrophilic, mesophilic, and thermophilic subtilisin-like serine proteases

To investigate the role of electrostatics in different temperature adaptations, we performed a comparative study on subtilisin-like serine proteases from psychrophilic Vibrio sp. PA-44 (VPR), mesophilic Engyodontium album (Tritirachium album) (PRK), and thermophilic Thermus aquaticus (AQN) using multiple-replica molecular dynamics (MD) simulations combined with continuum electrostatics calculations. The results reveal that although salt bridges are not a crucial factor in determining the overall thermostability of these three proteases, they on average provide the greatest, moderate, and least electrostatic stabilization to AQN, PRK, and VPR, respectively, at the respective organism growth temperatures. Most salt bridges in AQN are effectively stabilizing and thus contribute to maintaining the overall structural stability at 343 K, while nearly half of the salt bridges in VPR interconvert between being stabilizing and being destabilizing, likely aiding in enhancing the local conformational flexibility at 283 K. The individual salt bridges, salt-bridge networks, and calcium ions contribute differentially to local stability and flexibility of these three enzyme structures, depending on their spatial distributions and electrostatic strengths. The shared negatively charged surface potential at the active center of the three enzymes may provide the active-center flexibility necessary for nucleophilic attack and proton transfer. The differences in distributions of the electro-negative, electro-positive, and electro-neutral potentials, particularly over the back surfaces of the three proteases, may modulate/affect not only protein solubility and thermostability but also structural stability and flexibility/rigidity. These results demonstrate that electrostatics contributes to both heat and cold adaptation of subtilisin-like serine proteases through fine-tuning, either globally or locally, the structural stability and conformational flexibility/rigidity, thus providing a foundation for further engineering and mutagenesis studies.

Differently charged surface patches contribute to temperature adaptation of subtilisin-like serine proteases through affecting/modulating the protein solubility and thermostability and the structural flexibility/rigidity/stability.

The discovery of novel heat-stable keratinases from Meiothermus taiwanensis WR-220 and other extremophiles

Billions of tons of keratin bio-wastes are generated by poultry industry annually but discarded that result in serious environmental pollution. Keratinase is a broad spectrum protease with the unique ability to degrade keratin, providing an eco-friendly way to convert keratin wastes to valuable amino acids. In this report, a feather-degrading thermophilic bacterium, Meiothermus taiwanensis WR-220, was investigated due to its ability to apparently complete feather decay at 65 °C in two days. By genomics, proteomics, and biochemical approaches, the extracellular heat-stable keratinase (MtaKer) from M. taiwanensis WR-220 was identified. The recombinant MtaKer (rMtaKer) possesses keratinolytic activities at temperatures ranging from 25 to 75 °C and pH from 4 to 11, with a maximum keratinolytic activity at 65 °C and pH 10. The phylogenetic and structural analysis revealed that MtaKer shares low sequence identity but high structural similarity with known keratinases. Accordingly, our findings have enabled the discovery of more keratinases from other extremophiles, Thermus and Deinococcus. Proteins encoded in the extremophiles shall be evolved to be functional in the extreme conditions. Hence, our study expands the current boundary of hunting keratinases that can tolerate extreme conditions for keratin wastes biorecycle and other industrial applications.

Molecular motions and free-energy landscape of serine proteinase K in relation to its cold-adaptation: a comparative molecular dynamics simulation study and the underlying mechanisms

The physicochemical bases for enzyme cold-adaptation remain elusive.

Discovery, Molecular Mechanisms, and Industrial Applications of Cold-Active Enzymes

Cold-active enzymes constitute an attractive resource for biotechnological applications. Their high catalytic activity at temperatures below 25°C makes them excellent biocatalysts that eliminate the need of heating processes hampering the quality, sustainability, and cost-effectiveness of industrial production. Here we provide a review of the isolation and characterization of novel cold-active enzymes from microorganisms inhabiting different environments, including a revision of the latest techniques that have been used for accomplishing these paramount tasks. We address the progress made in the overexpression and purification of cold-adapted enzymes, the evolutionary and molecular basis of their high activity at low temperatures and the experimental and computational techniques used for their identification, along with protein engineering endeavors based on these observations to improve some of the properties of cold-adapted enzymes to better suit specific applications. We finally focus on examples of the evaluation of their potential use as biocatalysts under conditions that reproduce the challenges imposed by the use of solvents and additives in industrial processes and of the successful use of cold-adapted enzymes in biotechnological and industrial applications.

A single mutation Gln142Lys doubles the catalytic activity of VPR, a cold adapted subtilisin-like serine proteinase

Structural comparisons of the cold adapted subtilase VPR and its thermophilic homologue, aqualysin I (AQUI) indicated the presence of additional salt bridges in the latter. Few of those appear to contribute significantly to thermal stability of AQUI. This includes a putative salt bridge between residues Lys142 and Glu172 as its deletion did not have any significant effect on its stability or activity (Jónsdóttir et al. (2014)). Insertion of this putative salt bridge into the structure of VPR, in a double mutant (VPRΔC_Q142K/S172E), however was detrimental to the stability of the enzyme. Incorporation of either the Q142K or S172E mutations into VPR, were found to significantly affect the catalytic properties of the enzyme. The single mutation Q142K was highly effective, as it increased the kcat and kcat/Km more than twofold. When the Q142K mutation was inserted into a thermostabilized, but a low activity mutant of VPR (VPRΔC_N3P/I5P), the activity increased about tenfold in terms of kcat and kcat/Km, while retaining the stability of the mutant. Molecular dynamics simulations of the single mutants were carried out to provide structural rationale for these experimental observations. Based on root mean square fluctuation (RMSF) profiles, the two mutants were more flexible in certain regions of the structure and the Q142K mutant had the highest overall flexibility of the three enzymes. The results suggest that weakening of specific H-bonds resulting from the mutations may be propagated over some distance giving rise to higher flexibility in the active site regions of the enzyme, causing higher catalytic activity in the mutants.

Shotgun proteomic analysis of tiger milk mushroom (Lignosus rhinocerotis) and the isolation of a cytotoxic fungal serine protease from its sclerotium

ETHNOPHARMACOLOGICAL RELEVANCE

The sclerotium of Lignosus rhinocerotis (Cooke) Ryvarden (tiger milk mushroom) has been traditionally used as a complementary and alternative medicine for cancer treatment by the local communities of Southeast Asia. Despite the continuous research interest in its antiproliferative activity, the identity of the bioactive compound(s) responsible has yet to be determined. This study aims to bridge the gap in existing research literature by using proteomics approach for investigation of the nature of the anticancer substance of L. rhinocerotis.

AIM OF THE STUDY

To elucidate the proteome of L. rhinocerotis TM02 sclerotium by protein mass spectrometry and to further isolate and identify the cytotoxic component(s) bearing anticancer potential.

MATERIALS AND METHODS

The proteome of L. rhinocerotis sclerotium was analyzed by label-free quantitative shotgun proteomics, using 1D-SDS-PAGE coupled with nano-ESI-LC-MS/MS based on the availability of its genome-sequence database. The cytotoxicity of L. rhinocerotis sclerotial extracts against human breast adenocarcinoma cells (MCF7) were assessed by MTT cytotoxicity assay prior to successive purification steps by a combination of gel filtration chromatography, ammonium sulfate precipitation, and anion exchange chromatography. Bioactive compound(s) in the extracts was identified by shotgun proteomics and N-terminal protein sequencing.

RESULTS

Several proteins with interesting biological activities including lectins, fungal immunomodulatory proteins, and several antioxidant proteins were identified from the proteome of L. rhinocerotis. A cytotoxic protein fraction (termed F5) which was partially purified from its sclerotial cold water extract F5 shows two distinct bands of 31 and 36 kDa in reducing SDS-PAGE and exhibited potent selective cytotoxicity against MCF7 cells with IC50 value of 3.00 ± 1.01 μg/ml. Both bands were identified to be serine protease by LC-MS/MS analysis. Phenylmethylsulfonyl fluoride, a specific serine protease inhibitor, inhibited both the proteolytic activity and cytotoxicity of F5, suggesting that the cytotoxicity of F5 is related to its protease activity.

CONCLUSIONS

This study provides the first comprehensive and semi-quantitative profiling of the proteome of L. rhinocerotis sclerotium. Further investigation into its selective cytotoxicity shows that a serine protease-like protein, termed F5, may be targeted for new anticancer agent development.

Advances in protease engineering for laundry detergents

Proteases are essential ingredients in modern laundry detergents. Over the past 30 years, subtilisin proteases employed in the laundry detergent industry have been engineered by directed evolution and rational design to tailor their properties towards industrial demands. This comprehensive review discusses recent success stories in subtilisin protease engineering. Advances in protease engineering for laundry detergents comprise simultaneous improvement of thermal resistance and activity at low temperatures, a rational strategy to modulate pH profiles, and a general hypothesis for how to increase promiscuous activity towards the production of peroxycarboxylic acids as mild bleaching agents. The three protease engineering campaigns presented provide in-depth analysis of protease properties and have identified principles that can be applied to improve or generate enzyme variants for industrial applications beyond laundry detergents.

Highly conserved salt bridge stabilizes a proteinase K subfamily enzyme, Aqualysin I, from Thermus aquaticus YT-1

The proteinase K subfamily enzymes, thermophilic Aqualysin I (AQN) from Thermus aquaticus YT-1 and psychrophilic serine protease (VPR) from Vibrio sp. PA-44, have six and seven salt bridges, respectively. To understand the possible significance of salt bridges in the thermal stability of AQN, we prepared mutant proteins in which amino acid residues participating in salt bridges common to proteinase K subfamily members and intrinsic to AQN were replaced to disrupt the bridges one at a time. Disruption of a salt bridge common to proteinase K subfamily enzymes in the D183N mutant resulted in a significant reduction in thermal stability, and a massive change in the content of the secondary structure was observed, even at 70°C, in the circular dichroism (CD) analysis. These results indicate that the common salt bridge Asp183-Arg12 is important in maintaining the conformation of proteinase K subfamily enzymes and suggest the importance of proximity between the regions around Asp183 and the N-terminal region around Arg12. Of the three mutants that lack an AQN intrinsic salt bridge, D212N was more prone to unfolding at 80°C than the wild-type enzyme. Similarly, D17N and E237Q were less thermostable than the wild-type enzyme, although this may be partially due to increased autolysis. The AQN intrinsic salt bridges appear to confer additional thermal stability to this enzyme. These findings will further our understanding of the factors involved in stabilizing protein structure.

Comparative assessment of multiple criteria for the in silico prediction of cross-reactivity of proteins to known allergens

Genetically modified crops are becoming important components of a sustainable food supply and must be brought to market efficiently while also safeguarding the public from cross-reactivity of novel proteins to known allergens. Bioinformatic assessments can help to identify proteins warranting further experimental checks for cross-reactivity. This study is a large-scale in silico evaluation of assessment criteria, including searches for: alignments between a query and an allergen having ≥ 35% identity over a length ≥ 80; any sequence (of some minimum length) found in both a query and an allergen; any alignment between a query and an allergen with an E-value below some threshold. The criteria and an allergen database (AllergenOnline) are used to assess 27,243 Viridiplantae proteins for potential allergenicity. (A protein is classed as a "real allergen" if it exceeds a test-specific level of identity to an AllergenOnline entry; assessment of real allergens in the query set is against a reduced database from which the identifying allergen has been removed.) Each criterion's ability to minimize false positives without increasing false negative levels of current methods is determined. At best, the data show a reduction in false positives to ∼6% (from ∼10% under current methods) without any increase in false negatives.

Psychrophily and catalysis

Polar and other low temperature environments are characterized by a low content in energy and this factor has a strong incidence on living organisms which populate these rather common habitats. Indeed, low temperatures have a negative effect on ectothermic populations since they can affect their growth, reaction rates of biochemical reactions, membrane permeability, diffusion rates, action potentials, protein folding, nucleic acids dynamics and other temperature-dependent biochemical processes. Since the discovery that these ecosystems, contrary to what was initially expected, sustain a rather high density and broad diversity of living organisms, increasing efforts have been dedicated to the understanding of the molecular mechanisms involved in their successful adaptation to apparently unfavorable physical conditions. The first question that comes to mind is: How do these organisms compensate for the exponential decrease of reaction rate when temperature is lowered? As most of the chemical reactions that occur in living organisms are catalyzed by enzymes, the kinetic and thermodynamic properties of cold-adapted enzymes have been investigated. Presently, many crystallographic structures of these enzymes have been elucidated and allowed for a rather clear view of their adaptation to cold. They are characterized by a high specific activity at low and moderate temperatures and a rather low thermal stability, which induces a high flexibility that prevents the freezing effect of low temperatures on structure dynamics. These enzymes also display a low activation enthalpy that renders them less dependent on temperature fluctuations. This is accompanied by a larger negative value of the activation entropy, thus giving evidence of a more disordered ground state. Appropriate folding kinetics is apparently secured through a large expression of trigger factors and peptidyl–prolyl cis/trans-isomerases.

Psychrophilic enzymes: from folding to function and biotechnology

Psychrophiles thriving permanently at near-zero temperatures synthesize cold-active enzymes to sustain their cell cycle. Genome sequences, proteomic, and transcriptomic studies suggest various adaptive features to maintain adequate translation and proper protein folding under cold conditions. Most psychrophilic enzymes optimize a high activity at low temperature at the expense of substrate affinity, therefore reducing the free energy barrier of the transition state. Furthermore, a weak temperature dependence of activity ensures moderate reduction of the catalytic activity in the cold. In these naturally evolved enzymes, the optimization to low temperature activity is reached via destabilization of the structures bearing the active site or by destabilization of the whole molecule. This involves a reduction in the number and strength of all types of weak interactions or the disappearance of stability factors, resulting in improved dynamics of active site residues in the cold. These enzymes are already used in many biotechnological applications requiring high activity at mild temperatures or fast heat-inactivation rate. Several open questions in the field are also highlighted.

Gly or Ala substitutions for Pro(210)Thr(211)Asn(212) at the β8-β9 turn of subtilisin Carlsberg increase the catalytic rate and decrease thermostability

A comparison of the primary structures among psychrophilic, mesophilic, and thermophilic subtilases revealed that the turn between the β8 and β9 strands (β8-β9 turn, BPN' numbering) of psychrophilic subtilases are more flexible than those of their mesophilic and thermophilic counterparts. To investigate the relationship between structure of this turn and enzyme activity as well as thermostability of mesophilic subtilisin Carlsberg (sC), we analyzed 6 mutants of sC with a single, double, or triple Gly or Ala substitutions for Pro(210)Thr(211)Asn(212) at the β8-β9 turn. Among the single Gly substitutions, the P210G substitution most significantly (1.5-fold) increased the specific activity on N-succinyl-Ala-Ala-Pro-Phe-p-nitroanilide (AAPF) substrate and 12-fold decreased the thermostability. All mutants tested showed the increased k(cat) for the AAPF substrate and reduced thermostability compared with the wild-type sC. The k(cat) values of the P210G, P210G/T211G, and P210G/T211G/N212G mutants were 1.5-, 1.7-, and 1.8-fold higher than that of the wild-type sC. There were significant positive correlations between k(cat) and thermal inactivation rates as well as k(cat) and K(m) of the wild-type and mutants. These results demonstrate that the structure of β8-β9 turn, despite its distance from the active site, has significant effects on the catalytic rate and thermostability of sC through a global network of intramolecular interactions and suggest that the lack of flexibility of this turn stabilizes the wild-type sC against thermal inactivation in compensation for some loss of catalytic activity.

Proteases from psychrotrophs: an overview

Proteases are hydrolytic enzymes which catalyze the total hydrolysis of proteins in to amino acids. Although proteolytic enzymes can be obtained from animals and plants but microorganisms are the preferred source for industrial applications in view of scientific and economical advantage. Among various groups of microbes, psychrotrophs are ideal candidates for enzymes production keeping in mind that enzymes active at low temperature and stable under alkaline condition, in presence of oxidants and detergents are in large demand as laundry additive. The proteases from psychrotrophs also find application in environmental bioremediation, food and molecular biology. During the previous two decades, proteases from psychrotrophs have received increased attention because of their wide range of applications, but the full potential of psychrotrophic proteases has not been exploited. This review focuses attention on the present status of knowledge on the production, optimization, molecular characteristics, applications, substrate specificity, and crystal structure of psychrotrophic proteases. The review will help in making strategies for exploitation of psychrotrophic protease resources and improvement of enzymes to obtain more robust proteases of industrial and biotechnological significance.

Molecular adaptations to psychrophily: the impact of 'omic' technologies

The ability of cold-adapted microorganisms (generally referred to as psychrophiles) to survive is the result of molecular evolution and adaptations which, together, counteract the potentially deleterious effects of low kinetic energy environments and the freezing of water. These physiological adaptations are seen at many levels. Against a background of detailed comparative protein structural analyses, the recent surge of psychrophile proteome, genome, metagenome and transcriptome sequence data has triggered a series of sophisticated analyses of changes in global protein composition. These studies have revealed consistent and statistically robust changes in amino acid composition, interpreted as evolutionary mechanisms designed to destabilise protein structures, as well as identifying the presence of novel genes involved in cold adaptation.

The crystal structures of the psychrophilic subtilisin S41 and the mesophilic subtilisin Sph reveal the same calcium-loaded state

We determine and compare the crystal structure of two proteases belonging to the subtilisin superfamily: S41, a cold-adapted serine protease produced by Antarctic bacilli, at 1.4 A resolution and Sph, a mesophilic serine protease produced by Bacillus sphaericus, at 0.8 A resolution. The purpose of this comparison was to find out whether multiple calcium ion binding is a molecular factor responsible for the adaptation of S41 to extreme low temperatures. We find that these two subtilisins have the same subtilisin fold with a root mean square between the two structures of 0.54 A. The final models for S41 and Sph include a calcium-loaded state of five ions bound to each of these two subtilisin molecules. None of these calcium-binding sites correlate with the high affinity known binding site (site A) found for other subtilisins. Structural analysis of the five calcium-binding sites found in these two crystal structures indicate that three of the binding sites have two side chains of an acidic residue coordinating the calcium ion, whereas the other two binding sites have either a main-chain carbonyl, or only one acidic residue side chain coordinating the calcium ion. Thus, we conclude that three of the sites are of high affinity toward calcium ions, whereas the other two are of low affinity. Because Sph is a mesophilic subtilisin and S41 is a psychrophilic subtilisin, but both crystal structures were found to bind five calcium ions, we suggest that multiple calcium ion binding is not responsible for the adaptation of S41 to low temperatures.

Effect of proline substitutions on stability and kinetic properties of a cold adapted subtilase

A cold adapted subtilisin-like serine proteinase from a Vibrio species is two amino acids shorter at the N-terminus than related enzymes adapted to higher temperatures and has a 15 residues' C-terminal extension relative to the highly homologous thermophilic enzyme aqualysin I from Thermus aquaticus. These enzymes are produced as pro-enzymes with an N-terminal chaperone sequence for correct folding and a C-terminal signal peptide for secretion, which are subsequently cleaved off by autocatalysis to give the mature enzyme. A truncated form of the Vibrio proteinase where the C-terminal extension was removed and two residues near the N-terminus were substituted with proline, to resemble the N- and C-terminal regions in aqualysin I, resulted in increased thermostability and diminished catalytic efficiency. The proline substitutions shift the site of autocatalytic cleavage at the N-terminus by two amino acids, apparently by rigidifying the terminal residues and support the formation of a beta-sheet that fixes the N-terminus to the main body of the protein.

Characteristics of mutants designed to incorporate a new ion pair into the structure of a cold adapted subtilisin-like serine proteinase

Structural comparisons of VPR, a subtilisin-like serine proteinase from a psychrotrophic Vibrio species and a thermophilic homologue, aqualysin I, have led us to hypothesize about the roles of different residues in the temperature adaptation of the enzymes. Some of these hypotheses are now being examined by analysis of mutants of the enzymes. The selected substitutions are believed to increase the stability of the cold adapted enzyme based on structural analysis of the thermostable structure. We report here on mutants, which were designed to incorporate an ion pair into the structure of VPR. The residues Asp17 and Arg259 are assumed to form an ion pair in aqualysin I. The cold adapted VPR contains Asn (Asn15) and Lys (Lys257) at corresponding sites in its structure. In VPR, Asn 15 is located on a surface loop with its side group pointing towards the side chain of Lys257. By substituting Asn15 by Asp (N15D) it was considered feasible that a salt bridge would form between the oppositely charged groups. To mimic further the putative salt bridge from the thermophile enzyme the corresponding double mutant (N15D/K257R) was also produced. The N15D mutation increased the thermal stability of VPR by approximately 3 degrees C, both in T(50%) and T(m). Addition of the K257R mutation did not however, increase the stability of the double mutant any further. Despite this stabilization of the VPR mutants the catalytic activity (k(cat)) against the substrate Suc-AAPF-NH-Np was increased in the mutants. Molecular dynamics simulations on wild type and the two mutant proteins suggested that indeed a salt bridge was formed in both cases. Furthermore, a truncated form of the N15D mutant (N15DDeltaC) was produced, lacking a 15 residue long C-terminal extended sequence not present in the thermophilic enzyme. In wild type VPR this supposedly moveable, negatively charged arm on the protein molecule might interfere with the new salt bridge introduced as a result of the N15D mutation. Removal of the C-terminal arm improved the thermal stability (T(m) approximately +1.5 degrees C) of the truncated enzyme (VPRDeltaC) as compared to the wild type VPR. Introduction of the N15D substitution into VPRDeltaC improved the thermal stability further by about 3 degrees C, or to about the same extent as in the wild type. However, contrary to what was observed for the wild type, the introduction of the putative salt bridge did not affect the catalytic properties (k(cat)) of the C-terminal truncated enzyme.

Structural insights into cold inactivation of tryptophanase and cold adaptation of subtilisin S41

A wide variety of enzymes can undergo a reversible loss of activity at low temperature, a process that is termed cold inactivation. This phenomenon is found in oligomeric enzymes such as tryptophanase (Trpase) and other pyridoxal phosphate dependent enzymes. On the other hand, cold-adapted, or psychrophilic enzymes, isolated from organisms able to thrive in permanently cold environments, have optimal activity at low temperature, which is associated with low thermal stability. Since cold inactivation may be considered "contradictory" to cold adaptation, we have looked into the amino acid sequences and the crystal structures of two families of enzymes, subtilisin and tryptophanase. Two cold adapted subtilisins, S41 and subtilisin-like protease from Vibrio, were compared to a mesophilic and a thermophilic subtilisins, as well as to four PLP-dependent enzymes in order to understand the specific surface residues, specific interactions, or any other molecular features that may be responsible for the differences in their tolerance to cold temperatures. The comparison between the psychrophilic and the mesophilic subtilisins revealed that the cold adapted subtilisins have a high content of acidic residues mainly found on their surface, making it charged. The analysis of the Trpases showed that they have a high content of hydrophobic residues on their surface. Thus, we suggest that the negatively charged residues on the surface of the subtilisins may be responsible for their cold adaptation, whereas the hydrophobic residues on the surface of monomeric Trpase molecules are responsible for the tetrameric assembly, and may account for their cold inactivation and dissociation.

Crystal structure of the cold-active aminopeptidase from Colwellia psychrerythraea, a close structural homologue of the human bifunctional leukotriene A4 hydrolase

The crystal structure of a cold-active aminopeptidase (ColAP) from Colwellia psychrerythraea strain 34H has been determined, extending the number of crystal structures of the M1 metallopeptidase family to four among the 436 members currently identified. In agreement with their sequence similarity, the overall structure of ColAP displayed a high correspondence with leukotriene A4 hydrolase (LTA4H), a human bifunctional enzyme that converts leukotriene A4 (LTA4) in the potent chemoattractant leukotriene B4. Indeed, both enzymes are composed of three domains, an N-terminal saddle-like domain, a catalytic thermolysin-like domain, and a less conserved C-terminal alpha-helical flat spiral domain. Together, these domains form a deep cavity harboring the zinc binding site formed by residues included in the conserved HEXXHX(18)H motif. A detailed structural comparison of these enzymes revealed several plausible determinants of ColAP cold adaptation. The main differences involve specific amino acid substitutions, loop content and solvent exposure, complexity and distribution of ion pairs, and differential domain flexibilities. Such elements may act synergistically to allow conformational flexibility needed for an efficient catalysis in cold environments. Furthermore, the region of ColAP corresponding to the aminopeptidase active site of LTA4H is much more conserved than the suggested LTA4 substrate binding region. This observation supports the hypothesis that this region of the LTA4H active site has evolved in order to fit the lipidic substrate.

The expanding diversity of serine hydrolases

Serine hydrolases use a hydroxyl of a serine, assisted by one or more other residues, to cleave peptide bonds. They belong to several different families whose general mechanism is well known. However, the subtle structural differences that have recently been observed across a variety of families shed light on their functional diversity, including variations in mechanism of action, differences in the modes of substrate binding, and substrate-assisted orientation of catalytic residues. Of particular interest are the Rhomboid family serine proteinases that are active within the plasma membrane, for which several new structures have been reported. Because these enzymes are involved in biological and pathological processes, many are becoming important targets of drug design.

Cold-adaptation in Sea-water-borne Signal Proteins: Sequence and NMR Structure of the Pheromone En-6 from the Antarctic Ciliate Euplotes nobilii

Ciliates of Euplotes species constitutively secrete pleiotropic protein pheromones, which are capable to function as prototypic autocrine growth factors as well as paracrine inducers of mating processes. This paper reports the amino acid sequence and the NMR structure of the pheromone En-6 isolated from the antarctic species Euplotes nobilii. The 63-residue En-6 polypeptide chain forms three alpha-helices in positions 18-25, 36-40 and 46-56, which are arranged in an up-down-up three-helix bundle forming the edges of a distorted trigonal pyramid. The base of the pyramid is covered by the N-terminal heptadecapeptide segment, which includes a 3(10)-turn of residues 3-6. This topology is covalently anchored by four long-range disulfide bonds. Comparison with the smaller pheromones of E. raikovi, a closely related species living in temperate waters, shows that the two-pheromone families have the same three-helix bundle architecture. It then appears that cold-adaptation of the En proteins is primarily related to increased lengths of the chain-terminal peptide segments and the surface-exposed loops connecting the regular secondary structures, and to the presence of solvent-exposed clusters of negatively charged side-chains.

Structural and Functional Properties of Isocitrate Dehydrogenase from the Psychrophilic Bacterium Desulfotalea psychrophila Reveal a Cold-active Enzyme with an Unusual High Thermal Stability

Isocitrate dehydrogenase (IDH) has been studied extensively due to its central role in the Krebs cycle, catalyzing the oxidative NAD(P)(+)-dependent decarboxylation of isocitrate to alpha-ketoglutarate and CO(2). Here, we present the first crystal structure of IDH from a psychrophilic bacterium, Desulfotalea psychrophila (DpIDH). The structural information is combined with a detailed biochemical characterization and a comparative study with IDHs from the mesophilic bacterium Desulfitobacterium hafniense (DhIDH), porcine (PcIDH), human cytosolic (HcIDH) and the hyperthermophilic Thermotoga maritima (TmIDH). DpIDH was found to have a higher melting temperature (T(m)=66.9 degrees C) than its mesophilic homologues and a suboptimal catalytic efficiency at low temperatures. The thermodynamic activation parameters indicated a disordered active site, as seen also for the drastic increase in K(m) for isocitrate at elevated temperatures. A methionine cluster situated at the dimeric interface between the two active sites and a cluster of destabilizing charged amino acids in a region close to the active site might explain the poor isocitrate affinity. On the other hand, DpIDH was optimized for interacting with NADP(+) and the crystal structure revealed unique interactions with the cofactor. The highly acidic surface, destabilizing charged residues, fewer ion pairs and reduced size of ionic networks in DpIDH suggest a flexible global structure. However, strategic placement of ionic interactions stabilizing the N and C termini, and additional ionic interactions in the clasp domain as well as two enlarged aromatic clusters might counteract the destabilizing interactions and promote the increased thermal stability. The structure analysis of DpIDH illustrates how psychrophilic enzymes can adjust their flexibility in dynamic regions during their catalytic cycle without compromising the global stability of the protein.

Effect of selected Ser/Ala and Xaa/Pro mutations on the stability and catalytic properties of a cold adapted subtilisin-like serine proteinase

A subtilisin-like serine proteinase from a psychrotrophic Vibrio species (VPR) shows distinct cold adapted traits regarding stability and catalytic properties, while sharing high sequence homology with enzymes adapted to higher temperatures. Based on comparisons of sequences and examination of 3D structural models of VPR and related enzymes of higher temperature origin, five sites were chosen to be subject to site directed mutagenesis. Three serine residues were substituted with alanine and two residues in loops were substituted with proline. The single mutations were combined to make double and triple mutants. The single Ser/Ala mutations had a moderately stabilizing effect and concomitantly decreased catalytic efficiency. Introducing a second Ser/Ala mutation did not have additive effect on stability; on the contrary a double Ser/Ala mutant had reduced stability with regard to both wild type and single mutants. The Xaa/Pro mutations stabilized the enzyme and did also tend to decrease the catalytic efficiency more than the Ser/Ala mutations.

Structural, thermodynamic, and mutational analyses of a psychrotrophic RNase HI

Ribonuclease (RNase) HI from the psychrotrophic bacterium Shewanella oneidensis MR-1 was overproduced in Escherichia coli, purified, and structurally and biochemically characterized. The amino acid sequence of MR-1 RNase HI is 67% identical to that of E. coli RNase HI. The crystal structure of MR-1 RNase HI determined at 2.0 A resolution was highly similar to that of E. coli RNase HI, except that the number of intramolecular ion pairs and the fraction of polar surface area of MR-1 RNase HI were reduced compared to those of E. coli RNase HI. The enzymatic properties of MR-1 RNase HI were similar to those of E. coli RNase HI. However, MR-1 RNase HI was much less stable than E. coli RNase HI. The stability of MR-1 RNase HI against heat inactivation was lower than that of E. coli RNase HI by 19 degrees C. The conformational stability of MR-1 RNase HI was thermodynamically analyzed by monitoring the CD values at 220 nm. MR-1 RNase HI was less stable than E. coli RNase HI by 22.4 degrees C in Tm and 12.5 kJ/mol in DeltaG(H2O). The thermodynamic stability curve of MR-1 RNase HI was characterized by a downward shift and increased curvature, which results in an increased DeltaCp value, compared to that of E. coli RNase HI. Site-directed mutagenesis studies suggest that the difference in the number of intramolecular ion pairs partly accounts for the difference in stability between MR-1 and E. coli RNases HI.

Characterizing structural features of cuticle-degrading proteases from fungi by molecular modeling

Background

Serine proteases secreted by nematode and insect pathogenic fungi are bio-control agents which have commercial potential for developing into effective bio-pesticides. A thorough understanding of the structural and functional features of these proteases would significantly assist with targeting the design of efficient bio-control agents.

Results

Structural models of serine proteases PR1 from entomophagous fungus, Ver112 and VCP1 from nematophagous fungi, have been modeled using the homology modeling technique based on the crystal coordinate of the proteinase K. In combination with multiple sequence alignment, these models suggest one similar calcium-binding site and two common disulfide bridges in the three cuticle-degrading enzymes. In addition, the predicted models of the three cuticle-degrading enzymes present an essentially identical backbone topology and similar geometric properties with the exception of a limited number of sites exhibiting relatively large local conformational differences only in some surface loops and the N-, C termini. However, they differ from each other in the electrostatic surface potential, in hydrophobicity and size of the S4 substrate-binding pocket, and in the number and distribution of hydrogen bonds and salt bridges within regions that are part of or in close proximity to the S2-loop.

Conclusion

These differences likely lead to variations in substrate specificity and catalytic efficiency among the three enzymes. Amino acid polymorphisms in cuticle-degrading enzymes were discussed with respect to functional effects and host preference. It is hoped that these structural models would provide a further basis for exploitation of these serine proteases from pathogenic fungi as effective bio-control agents.