Effects of changing the interaction between subdomains on the thermostability of Bacillus neutral proteases

NMR Structure and Biophysical Characterization of Thermophilic Single-Stranded DNA Binding Protein from Sulfolobus Solfataricus

Proteins from Sulfolobus solfataricus (S. solfataricus), an extremophile, are active even at high temperatures. The single-stranded DNA (ssDNA) binding protein of S. solfataricus (SsoSSB) is overexpressed to protect ssDNA during DNA metabolism. Although SsoSSB has the potential to be applied in various areas, its structural and ssDNA binding properties at high temperatures have not been studied. We present the solution structure, backbone dynamics, and ssDNA binding properties of SsoSSB at 50 °C. The overall structure is consistent with the structures previously studied at room temperature. However, the loop between the first two β sheets, which is flexible and is expected to undergo conformational change upon ssDNA binding, shows a difference from the ssDNA bound structure. The ssDNA binding ability was maintained at high temperature, but different interactions were observed depending on the temperature. Backbone dynamics at high temperature showed that the rigidity of the structured region was well maintained. The investigation of an N-terminal deletion mutant revealed that it is important for maintaining thermostability, structure, and ssDNA binding ability. The structural and dynamic properties of SsoSSB observed at high temperature can provide information on the behavior of proteins in thermophiles at the molecular level and guide the development of new experimental techniques.

Application of Rigidity Theory to the Thermostabilization of Lipase A from Bacillus subtilis

Protein thermostability is a crucial factor for biotechnological enzyme applications. Protein engineering studies aimed at improving thermostability have successfully applied both directed evolution and rational design. However, for rational approaches, the major challenge remains the prediction of mutation sites and optimal amino acid substitutions. Recently, we showed that such mutation sites can be identified as structural weak spots by rigidity theory-based thermal unfolding simulations of proteins. Here, we describe and validate a unique, ensemble-based, yet highly efficient strategy to predict optimal amino acid substitutions at structural weak spots for improving a protein’s thermostability. For this, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state has been found as the cause for thermostability. When applied prospectively to lipase A from Bacillus subtilis, we achieved both a high success rate (25% over all experimentally tested mutations, which raises to 60% if small-to-large residue mutations and mutations in the active site are excluded) in predicting significantly thermostabilized lipase variants and a remarkably large increase in those variants’ thermostability (up to 6.6°C) based on single amino acid mutations. When considering negative controls in addition and evaluating the performance of our approach as a binary classifier, the accuracy is 63% and increases to 83% if small-to-large residue mutations and mutations in the active site are excluded. The gain in precision (predictive value for increased thermostability) over random classification is 1.6-fold (2.4-fold). Furthermore, an increase in thermostability predicted by our approach significantly points to increased experimental thermostability (p < 0.05). These results suggest that our strategy is a valuable complement to existing methods for rational protein design aimed at improving thermostability.

Protein thermostability is a crucial factor for biotechnological enzyme applications. However, performance studies of computational approaches for predicting effects of mutations on protein (thermo)stability have suggested that there is still room for improvement. We describe and validate a novel and unique strategy to predict optimal amino acid substitutions at structural weak spots. At variance with other rational approaches, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state is the underlying cause for thermostability. When applied prospectively on lipase LipA from Bacillus subtilis, a high success rate in predicting thermostabilized lipase variants and a remarkably large increase in their thermostability is achieved. This demonstrates the value of the novel strategy, which extends the existing portfolio of methods for rational protein design.

Evaluation and directed evolution for thermostability improvement of a GH 13 thermostable α-glucosidase from Thermus thermophilus TC11

Background

Thermal stable α-glucosidases with transglycosylation activity could be applied to the industrial production of oligosaccharides as well as conjugation of sugars to biologically useful materials. Therefore, α-glucosidases isolated from thermophiles have gained attention over the past decade. In this study, the characterization of a highly thermostable α-glucosidase and its thermostability improved mutant from newly isolated strain Thermus thermophilus TC11 were investigated.

Results

The recombinant α-glucosidase (TtAG) from Thermus thermophilus TC11 was expressed in Escherichia coli BL21 (DE3) and purified. The purified enzyme had a molecular mass of 184 kDa and consisted of 59-kDa subunits; it showed hydrolytic activity for pNP-α-d-glucopyranoside (pNPG), sucrose, trehalose, panose, and isomaltooligosaccharides and very low activity for maltose. The highest specific activity of 288.96 U/mg was observed for pNPG at 90 °C and pH 5.0; Pb²⁺ provided a 20 % activity increase. TtAG was stable at 70 °C for more than 7 h and had a half-life of 195 min at 80 °C and 130 min at 90 °C. Transglycosylation activity was also observed with sucrose and trehalose as substrates. TtAG showed differences on substrate specificity, transglycosylation, multimerization, effects of metal ions and optimal pH from other reported Thermus α-glucosidases. One single-substitution TtAG mutant Q10Y with improved thermostability was also obtained from random mutagenesis library. The site-saturation mutagenesis and structural modelling analysis indicated that Q10Y substitution stabilized TtAG structure via additional hydrogen bonding and hydrophobic interactions.

Conclusion

Our findings indicate that TtAG is a highly thermostable and more acidic α-glucosidase distinct from other reported Thermus α-glucosidases. And this work also provides new insights into the catalytic and thermal tolerance mechanisms of α-glucosidases, which may guide molecular engineering of α-glucosidase and other thermostable enzymes for industrial application.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-015-0197-x) contains supplementary material, which is available to authorized users.

Stabilization of cyclohexanone monooxygenase by a computationally designed disulfide bond spanning only one residue

•

Cyclohexanone monooxygenase was stabilized by an in silico designed disulfide bond.

•

Stabilizing disulfide bonds were successfully designed based on a model structure.

•

The half-life at 30 °C was increased 12-fold for the mutant enzyme.

•

The apparent melting point was increased by 6 °C for the mutant enzyme.

•

The most stabilizing disulfide bond spans only one residue.

Enzyme stability is an important parameter in biocatalytic applications, and there is a strong need for efficient methods to generate robust enzymes. We investigated whether stabilizing disulfide bonds can be computationally designed based on a model structure. In our approach, unlike in previous disulfide engineering studies, short bonds spanning only a few residues were included. We used cyclohexanone monooxygenase (CHMO), a Baeyer–Villiger monooxygenase (BVMO) from Acinetobacter sp. NCIMB9871 as the target enzyme. This enzyme has been the prototype BVMO for many biocatalytic studies even though it is notoriously labile. After creating a small library of mutant enzymes with introduced cysteine pairs and subsequent screening for improved thermostability, three stabilizing disulfide bonds were identified. The introduced disulfide bonds are all within 12 Å of each other, suggesting this particular region is critical for unfolding. This study shows that stabilizing disulfide bonds do not have to span many residues, as the most stabilizing disulfide bond, L323C–A325C, spans only one residue while it stabilizes the enzyme, as shown by a 6 °C increase in its apparent melting temperature.

Enhancing the thermostability of alpha-glucosidase from Thermoanaerobacter tengcongensis MB4 by single proline substitution

Thermostability can be increased by introducing prolines at suitable sites in target proteins. In this study, we compared five thermostable alpha-glucosidases and the moderate thermostable alpha-glucosidase (TtGluA) from Thermoanaerobacter tengcongensis MB4. Based on the amino acid sequence alignment, four sites (Leu152, Asn208, Lys285, and Thr430) of TtGluA were chosen for proline substitution to improve its thermostability. Thermostability of mutants L152P, K285P, and T430P increased evidently, but no thermostability improvement was observed for N208P. Compared to the wild-type enzyme, T(50)(15) of T430P had a rise of 2 degrees C without distinct loss of activity. However, T(50)(15) values of L152P and K285P increased 2 degrees C and 10.5 degrees C, respectively, while retaining activity of only 26.6% and 24.9% of wild-type enzyme. The K(m) of L152P, K285P, T430P and wild-type enzyme was 1.61, 0.32, 1.64, and 1.08 mM, respectively. These indicate that the selected sites are not only important for the thermostability but also related to the substrate binding and catalytic activity of TtGluA. The CD spectra analysis of the improved mutants and wild-type enzyme showed no distinct changes in their secondary structures. Combining analysis of secondary structure prediction and 3D structure modeling, the proline substitution at the three sites stabilized TtGluA possibly by reducing the flexibility of loop and random coil or (and) increasing the hydrophobic effect at these strategic regions with no evident structure change.

Rational engineering of enzyme stability

During the past 15 years there has been a continuous flow of reports describing proteins stabilized by the introduction of mutations. These reports span a period from pioneering rational design work on small enzymes such as T4 lysozyme and barnase to protein design, and directed evolution. Concomitantly, the purification and characterization of naturally occurring hyperstable proteins has added to our understanding of protein stability. Along the way, many strategies for rational protein stabilization have been proposed, some of which (e.g. entropic stabilization by introduction of prolines or disulfide bridges) have reasonable success rates. On the other hand, comparative studies and efforts in directed evolution have revealed that there are many mutational strategies that lead to high stability, some of which are not easy to define and rationalize. Recent developments in the field include increasing awareness of the importance of the protein surface for stability, as well as the notion that normally a very limited number of mutations can yield a large increase in stability. Another development concerns the notion that there is a fundamental difference between the "laboratory stability" of small pure proteins that unfold reversibly and completely at high temperatures and "industrial stability", which is usually governed by partial unfolding processes followed by some kind of irreversible inactivation process (e.g. aggregation). Provided that one has sufficient knowledge of the mechanism of thermal inactivation, successful and efficient rational stabilization of enzymes can be achieved.

The effect of changing the hydrophobic S1' subsite of thermolysin-like proteases on substrate specificity

The hydrophobic S1' subsite is one of the major determinants of the substrate specificity of thermolysin and related M4 family proteases. In the thermolysin-like protease (TLP) produced by Bacillus stearothermophilus (TLP-ste), the hydrophobic S1' subsite is mainly formed by Phe130, Phe133, Val139 and Leu202. In the present study, we have examined the effects of replacing Leu202 by smaller (Gly, Ala, Val) and larger (Phe, Tyr) hydrophobic residues. The mutational effects showed that the wild-type S1' pocket is optimal for binding leucine side chains. Reduction of the size of residue 202 resulted in a higher efficiency towards substrates with Phe in the P1' position. Rather unexpectedly, the Leu202-->Phe and Leu202-->Tyr mutations, which were expected to decrease the size of the S1' subsite, resulted in a large increase in activity towards dipeptide substrates with Phe in the P1' position. This is probably due to the fact that 202Phe and 202Tyr adopt a second possible rotamer that opens up the subsite compared to Leu202, and also favours interactions with the substrate. To validate these results, we constructed variants of thermolysin with changes in the S1' subsite. Thermolysin and TLP-ste variants with identical S1' subsites were highly similar in terms of their preference for Phe vs. Leu in the P1' position.

Substrate specificity in the highly heterogeneous M4 peptidase family is determined by a small subset of amino acids

The members of the M4 peptidase family are involved in processes as diverse as pathogenicity and industrial applications. For the first time a number of M4 family members, also known as thermolysin-like proteases, has been characterized with an identical substrate set and a uniform set of assay conditions. Characterization with peptide substrates as well as high performance liquid chromatography analysis of beta-casein digests shows that the M4 family is a homogeneous family in terms of catalysis, even though there is a significant degree of amino acid sequence variation. The results of this study show that differences in substrate specificity within the M4 family do not correlate with overall sequence differences but depend on a small number of identifiable amino acids. Indeed, molecular modeling followed by site-directed mutagenesis of one of the substrate binding pocket residues of the thermolysin-like proteases of Bacillus stearothermophilus converted the catalytic characteristics of this variant into that of thermolysin.

Directed evolution converts subtilisin E into a functional equivalent of thermitase

We used directed evolution to convert Bacillus subtilis subtilisin E into an enzyme functionally equivalent to its thermophilic homolog thermitase from Thermoactinomyces vulgaris. Five generations of random mutagenesis, recombination and screening created subtilisin E 5-3H5, whose half-life at 83 degrees C (3.5 min) and temperature optimum for activity (Topt, 76 degrees C) are identical with those of thermitase. The Topt of the evolved enzyme is 17 degrees C higher and its half-life at 65 degrees C is >200 times that of wild-type subtilisin E. In addition, 5-3H5 is more active towards the hydrolysis of succinyl-Ala-Ala-Pro-Phe-p-nitroanilide than wild-type at all temperatures from 10 to 90 degrees C. Thermitase differs from subtilisin E at 157 amino acid positions. However, only eight amino acid substitutions were sufficient to convert subtilisin E into an enzyme equally thermostable. The eight substitutions, which include known stabilizing mutations (N218S, N76D) and also several not previously reported, are distributed over the surface of the enzyme. Only two (N218S, N181D) are found in thermitase. Directed evolution provides a powerful tool to unveil mechanisms of thermal adaptation and is an effective and efficient approach to increasing thermostability without compromising enzyme activity.

Early steps in the unfolding of thermolysin-like proteases

Several series of site-directed mutations in thermolysin-like proteases are presented that show remarkable nonadditivity in their effect on thermal stability. A simple model is proposed that relates this nonadditivity to the occurrence of independent partial unfolding processes that occur in parallel at elevated temperatures. To prove this model, a thermolysin-like protease was designed in which two mutations located approximately 35 A apart in the structure individually exert small stabilizing effects of 2.3 and 4. 1 degreesC, respectively, but when combined stabilize the protease by 14.6 degreesC. This overadditivity, which follows directly from the model, confirms that unfolding of this engineered protease starts in parallel at two different regions of the protein.

Mutational analysis of a surface area that is critical for the thermal stability of thermolysin-like proteases

Site-directed mutagenesis was used to assess the contribution of individual residues and a bound calcium in the 55-69 region of the thermolysin-like protease of Bacillus stearothermophilus (TLP-ste) to thermal stability. The importance of the 55-69 region was reflected by finding that almost all mutations had drastic effects on stability. These effects (both stabilizing and destabilizing) were obtained by mutations affecting main chain flexibility, as well as by mutations affecting the interaction between the 55-69 region and the rest of the protease molecule. The calcium-dependency of stability could be largely abolished by mutating one of its ligands (Asp57 or Asp59). In the case of the Asp57-->Ser mutation, the accompanying loss in stability was modest compared with the effects of other destabilizing mutations or the effects of (combinations of) stabilizing mutations. The detailed knowledge of the stability-determining region of TLP-ste permits effective rational design of stabilizing mutations, which, presumably, are also useful for related TLP such as thermolysin. This is demonstrated by the successful design of a stabilizing salt bridge involving residues 65 and 11.

Extreme stabilization of a thermolysin-like protease by an engineered disulfide bond

The thermal inactivation of broad specificity proteases such as thermolysin and subtilisin is initiated by partial unfolding processes that render the enzyme susceptible to autolysis. Previous studies have revealed that a surface-located region in the N-terminal domain of the thermolysin-like protease produced by Bacillus stearothermophilus is crucial for thermal stability. In this region a disulfide bridge between residues 8 and 60 was designed by molecular modelling, and the corresponding single and double cysteine mutants were constructed. The disulfide bridge was spontaneously formed in vivo and resulted in a drastic stabilization of the enzyme. This stabilization presents one of the very few examples of successful stabilization of a broad specificity protease by a designed disulfide bond. We propose that the success of the present stabilization strategy is the result of the localization and mutation of an area of the molecule involved in the partial unfolding processes that determine thermal stability.

TRAJAN: a tool for analyzing trajectories from molecular simulations

Molecular dynamics simulations of biological systems are notoriously difficult to analyze because of the complexity of the information that contain. We describe a new method for analyzing trajectories from simulations in order to extract important features of the motion. The trajectory of each atom is partitioned into conformation wells, in each of which its motion is assumed to be predominantly harmonic oscillation around an average position. Thus each atom moves anharmonically through a sequence of wells during the trajectory. The movement of atoms between wells, their ellipsoids of motion within each well, and correlations in the motion of atoms are quantified and can be visualized with molecular graphics. The TRAJAN analysis procedure is applicable to trajectories from equilibrium and nonequilibrium simulations and is not restricted to molecular dynamics simulations. Its application is demonstrated for a range of model systems.

Structural determinants of the stability of thermolysin-like proteinases

Thermolysin is a member of a family of homologous proteinases which differ in their resistance to thermally induced unfolding and subsequent autolytic degradation. Site-directed mutagenesis studies of the thermolysin-like proteinase (TLP) from Bacillus stearothermophilus (TLP-ste) show that its reduced resistance to thermally induced autolysis, as compared to thermolysin, is due to only some of the 44 naturally occurring amino-acid differences between them. In fact TLP-ste becomes more resistant than thermolysin by mutation of just a few of these amino-acids. The crucial differences are all localized to a solvent-exposed region in the N-terminal domain of TLP-ste.

Modeling studies of the change in conformation required for cleavage of limited proteolytic sites

Previous analyses of limited proteolytic sites within native, folded protein structures have shown that a significant conformational change is required in order to facilitate binding into the active site of the attacking proteinase. For the serine proteinases, the optimum conformation to match the proteinase binding-site geometry has been well characterized crystallographically by the conserved main-chain geometry of the reactive site loops of their protein inhibitors. A good substrate must adopt a conformation very similar to this "target" main-chain conformation prior to cleavage. Using a "loop-closure" modeling approach, we have tested the ability of a set of tryptic-limited proteolytic sites to achieve this target conformation and further tested their suitability for cleavage. The results show that in most cases, significant changes in the conformation of at least 12 residues are required. All the putative tryptic cleavage sites in 1 protein, elastase, were also modeled and tested to compare the results to the actual nicksite in that protein. These results strongly suggest that large local motions proximate to the scissile bond are required for proteolysis, and it is this ability to unfold locally without perturbing the overall protein conformation that is the prime determinant for limited proteolysis.

Protein stabilization by hydrophobic interactions at the surface

The contribution of the solvent-exposed residue 63 to thermal stability of the thermolysin-like neutral protease of Bacillus stearothermophilus was studied by analyzing the effect of twelve different amino acid substitutions at this position. The thermal stability of the enzyme was increased considerably by introducing Arg, Lys or bulky hydrophobic amino acids. In general, the effects of the mutations showed that hydrophobic contacts in this surface-located region of the protein are a major determinant of thermal stability. This observation contrasts with general concepts concerning the contribution of surface-located residues and surface hydrophobicity to protein stability and indicates new ways for protein stabilization by site-directed mutagenesis.

Prediction and analysis of structure, stability and unfolding of thermolysin-like proteases

Bacillus neutral proteases (NPs) form a group of well-characterized homologous enzymes, that exhibit large differences in thermostability. The three-dimensional (3D) structures of several of these enzymes have been modelled on the basis of the crystal structures of the NPs of B. thermoproteolyticus (thermolysin) and B. cereus. Several new techniques have been developed to improve the model-building procedures. Also a 'model-building by mutagenesis' strategy was used, in which mutants were designed just to shed light on parts of the structures that were particularly hard to model. The NP models have been used for the prediction of site-directed mutations aimed at improving the thermostability of the enzymes. Predictions were made using several novel computational techniques, such as position-specific rotamer searching, packing quality analysis and property-profile database searches. Many stabilizing mutations were predicted and produced: improvement of hydrogen bonding, exclusion of buried water molecules, capping helices, improvement of hydrophobic interactions and entropic stabilization have been applied successfully. At elevated temperatures NPs are irreversibly inactivated as a result of autolysis. It has been shown that this denaturation process is independent of the protease activity and concentration and that the inactivation follows first-order kinetics. From this it has been conjectured that local unfolding of (surface) loops, which renders the protein susceptible to autolysis, is the rate-limiting step. Despite the particular nature of the thermal denaturation process, normal rules for protein stability can be applied to NPs. However, rather than stabilizing the whole protein against global unfolding, only a small region has to be protected against local unfolding. In contrast to proteins in general, mutational effects in proteases are not additive and their magnitude is strongly dependent on the location of the mutation. Mutations that alter the stability of the NP by a large amount are located in a relatively weak region (or more precisely, they affect a local unfolding pathway with a relatively low free energy of activation). One weak region, that is supposedly important in the early steps of NP unfolding, has been determined in the NP of B. stearothermophilus. After eliminating this weakest link a drastic increase in thermostability was observed and the search for the second-weakest link, or the second-lowest energy local unfolding pathway is now in progress. Hopefully, this approach can be used to unravel the entire early phase of unfolding.

Stabilization of Bacillus stearothermophilus neutral protease by introduction of prolines

The thermostability of neutral proteases has been shown to depend on autolysis which presumably occurs in flexible regions of the protein. In an attempt to rigidify such a region in the neutral protease of Bacillus stearothermophilus, residues in the solvent-exposed 63-69 loop were replaced by proline. The mutations caused large positive (Ser-65-->Pro, Ala-69-->Pro) or negative (Thr-63-->Pro, Tyr-66-->Pro) changes in thermostability, which were explained on the basis of molecular modelling of the mutant proteins. The data show that the introduction of prolines at carefully selected positions in the protein can be a powerful method for stabilization.