Proteases of enhanced stability: characterization of a thermostable variant of subtilisin

A procedure has been developed for the isolation and identification of mutants in the bacterial serine protease subtilisin that exhibit enhanced thermal stability. The cloned subtilisin BPN' gene from Bacillus amyloliquefaciens was treated with bisulfite, a chemical mutagen that deaminates cytosine to uracil in single-stranded DNA. Strains containing the cloned, mutagenized subtilisin gene which produced subtilisin with enhanced thermal stability were selected by a simple plate assay procedure which screens for esterase activity on nitrocellulose filters after preincubation at elevated temperatures. One thermostable subtilisin variant, designated 7150, has been fully characterized and found to differ from wild-type subtilisin by a single substitution of Ser for Asn at position 218. The 7150 enzyme was found to undergo thermal inactivation at one-fourth the rate of the wild-type enzyme when incubated at elevated temperatures. Moreover, the mid-point in the thermally induced transition from the folded to unfolded state was found to be 2.4-3.9 degrees C higher for 7150 as determined by differential scanning calorimetry under a variety of conditions. The refined, 1.8-A crystal structures of the wild-type and 7150 subtilisin have been compared in detail, leading to the conclusion that slight improvements in hydrogen bond parameters in the vicinity of position 218 result in the enhanced thermal stability of 7150.     

Isolation of subtilisin pro-sequence mutations that affect formation of active protease by localized random polymerase chain reaction mutagenesis

In order to analyze the role of the pro-sequence in folding of the alkaline serine protease subtilisin, localized random mutagenesis using the polymerase chain reaction with Taq DNA polymerase was employed to obtain mutations in the pro-sequence which prevent production of active protease. The unique aspect of this procedure is that random mutations can be easily generated in vitro over large but defined regions of a specific gene. The method was applied to a 458-base pair fragment encompassing the coding region of the pro-sequence of subtilisin, a region of the protein which has been shown to be required for proper folding. Protease-deficient mutants containing a variety of amino acid substitutions were isolated with a frequency of 4.3%. From analysis of these mutants, four independent amino acid substitution mutations in the pro-sequence were identified. The present results demonstrate that polymerase chain reaction is an efficient and simple method for obtaining random mutations within a localized region of a given gene.     

A natural variant of Bacillus licheniformis alpha-amylase isolated from flour mill wastewaters sheds light on the origin of high thermostability

AIMS: Understanding the origin of high thermostability exhibited by the alpha-amylase produced by a natural strain of Bacillus licheniformis. METHODS AND RESULTS: The MSH320 alpha-amylase gene has been cloned from a native strain of B. licheniformis isolated from flour mill wastewaters in Kashan, central Iran, and its nucleotide sequence was determined (GenBank Accession Number AF438149). Whereas previously cloned B. licheniformisalpha-amylase (BLA) genes are nearly identical, the MSH320 gene coding sequence presents only 93% identity with the reference 'wild-type' BLA gene, most of the nucleotide changes leading to silent mutations. Amino acid substitutions occurred at 19 of the 483 residues of the matured protein, distributed all along the protein sequence. Nevertheless, the natural BLA variant presents thermoinactivation kinetics similar to that of the reference BLA. Protein modelling and structural predictions at the substitution sites suggest that half of the mutations may have a significant stabilizing or destabilizing effect on the protein structure. Compensatory mutations thus occurred in the natural variant in order to maintain thermostability to the level of the reference enzyme. CONCLUSIONS: The exceptional high thermostability of BLA, although produced by a nonthermophilic organism, is not fortuitous but subject to a selective pressure still at work in natural environments. SIGNIFICANCE AND IMPACT OF THE STUDY: BLA thermal performances are not naturally maximized and can be substantially improved by protein engineering.     

Escherichia coli SecB stimulates export without maintaining export competence of ribose-binding protein signal sequence mutants.

Ribose-binding protein (RBP) is exported to the periplasm of Escherichia coli via the general export pathway. An rbsB-lacZ gene fusion was constructed and used to select mutants defective in RBP export. The spontaneous Lac+ mutants isolated in this selection contained either single-amino-acid substitutions or a deletion of the RBP signal sequence. Intact rbsB genes containing eight different point mutations in the signal sequence were reconstructed, and the effects of the mutations on RBP export were examined. Most of the mutations caused severe defects in RBP export. In addition, different suppressor mutations in SecY/PrlA protein were analyzed for their effects on the export of RBP signal sequence mutants in the presence or absence of SecB. Several RBP signal sequence mutants were efficiently suppressed, but others were not suppressed. Export of an RBP signal sequence mutant in prlA mutant strains was partially dependent on SecB, which is in contrast to the SecB independence of wild-type RBP export. However, the kinetics of export of an RBP signal sequence mutant point to a rapid loss of pre-RBP export competence, which occurs in strains containing or lacking SecB. These results suggest that SecB does not stabilize the export-competent conformation of RBP and may affect translocation by stabilizing the binding of pre-RBP at the translocation site.     

Stabilizing the subtilisin BPN' pro-domain by phage display selection: How restrictive is the amino acid code for maximum protein stability?

We have devised a procedure using monovalent phage display to select for stable mutants in the pro-domain of the serine protease, subtilisin BPN'. In complex with subtilisin, the pro-domain assumes a compact structure with a four-stranded antiparallel beta-sheet and two three-turn alpha-helices. When isolated, however, the pro-domain is 97% unfolded. These experiments use combinatorial mutagenesis to select for stabilizing amino acid combinations at a particular structural locus and determine how many combinations are close to the maximum protein stability. The selection for stability is based on the fact that the independent stability of the pro-domain is very low and that binding to subtilisin is thermodynamically linked to folding. Two libraries of mutant pro-domains were constructed and analyzed to determine how many combinations of amino acids at a particular structural locus result in the maximum stability. A library comprises all combinations of four amino acids at a structural locus. Previous studies using combinatorial genetics have shown that many different combinations of amino acids can be accommodated in a selected locus without destroying function. The present results indicate that the number of sequence combinations at a structural locus, which are close to the maximum stability, is small. The most striking example is a selection at an interior locus of the pro-domain. After two rounds of phagemid selection, one amino acid combination is found in 40% of sequenced mutants. The most frequently selected mutant has a deltaG(unfolding) = 4 kcal/mol at 25 degrees C, an increase of 6 kcal/mol relative to the naturally occurring sequence. Some implications of these results on the amount of sequence information needed to specify a unique tertiary fold are discussed. Apart from possible implications on the folding code, the phage display selection described here should be useful in optimizing the stability of other proteins, which can be displayed on the phage surface.     

Protein engineering of subtilisin

The serine protease subtilisin is an important industrial enzyme as well as a model for understanding the enormous rate enhancements affected by enzymes. For these reasons along with the timely cloning of the gene, ease of expression and purification and availability of atomic resolution structures, subtilisin became a model system for protein engineering studies in the 1980s. Fifteen years later, mutations in well over 50% of the 275 amino acids of subtilisin have been reported in the scientific literature. Most subtilisin engineering has involved catalytic amino acids, substrate binding regions and stabilizing mutations. Stability has been the property of subtilisin which has been most amenable to enhancement, yet perhaps least understood. This review will give a brief overview of the subtilisin engineering field, critically review what has been learned about subtilisin stability from protein engineering experiments and conclude with some speculation about the prospects for future subtilisin engineering.     

Replacement of the Bacillus subtilis subtilisin structural gene with an In vitro-derived deletion mutation.

The entire subtilisin structural gene from Bacillus subtilis I168 has been cloned, and its nucleotide sequence has been determined. When expressed on a high-copy-number shuttle vector, a fivefold increase in serine protease activity was observed. The DNA sequence of the gene is 80% homologous to the Bacillus amyloliquefaciens subtilisin structural gene, and the translated mature coding sequence is 85% homologous to the published protein sequence of subtilisin BPN'. The chloramphenicol resistance determinant of a plasmid integrated at the subtilisin locus was mapped by PBS1 transduction and was found to be linked to glyB (83%) and argC (60%), but not with metC or purB . The chromosomal locus containing the wild-type subtilisin allele was replaced with an in vitro-derived allele of the gene (delta apr-684) that contained a 684-base-pair deletion. The technique used for introducing the deletion is a variation of the gene replacement methods used in Saccharomyces cerevisiae and Escherichia coli. When used in B. subtilis, deletion mutants could be directly screened among the transformants. Physiological characterization of the delta apr-684 mutation revealed no discernable effect on the formation of heat-resistant endospores, but strains carrying the mutation produced only 10% of wild-type serine protease activity. A model is presented that outlines the pathway for plasmid integration and deletion formation in B. subtilis.     

Protease-deficient Bacillus subtilis host strains for production of Staphylococcal protein A.

We constructed strains of Bacillus subtilis which produced very low levels of extracellular proteases. These strains carried insertion or deletion mutations in the subtilisin structural gene (apr) which were constructed in vitro by using the cloned gene. The methods used to construct the mutations involved the use of a plasmid vector which allowed the selection of chromosomal integrates and their subsequent excision by homologous recombination to effect replacement of the chromosomal apr gene by a derivative carrying an inactivating insert with a selectable marker (a cat gene conferring chloramphenicol resistance). The strains produced no subtilisin, no detectable extracellular metalloprotease activity, and residual extracellular serine protease levels as low as 0.5% of that of the standard strain from which they were derived. The strains proved to be superior host strains for the production of staphylococcal protein A, accumulating higher levels of intact protein than do previously available B. subtilis strains.     

Protease-deficient Bacillus subtilis host strains for production of Staphylococcal protein A

We constructed strains of Bacillus subtilis which produced very low levels of extracellular proteases. These strains carried insertion or deletion mutations in the subtilisin structural gene (apr) which were constructed in vitro by using the cloned gene. The methods used to construct the mutations involved the use of a plasmid vector which allowed the selection of chromosomal integrates and their subsequent excision by homologous recombination to effect replacement of the chromosomal apr gene by a derivative carrying an inactivating insert with a selectable marker (a cat gene conferring chloramphenicol resistance). The strains produced no subtilisin, no detectable extracellular metalloprotease activity, and residual extracellular serine protease levels as low as 0.5% of that of the standard strain from which they were derived. The strains proved to be superior host strains for the production of staphylococcal protein A, accumulating higher levels of intact protein than do previously available B. subtilis strains.     

Structural basis of selection and thermostability of laboratory evolved Bacillus subtilis lipase

Variation in gene sequences generated by directed evolution approaches often does not assure a minimalist design for obtaining a desired property in proteins. While screening for enhanced thermostability, structural information was utilized in selecting mutations that are generated by error-prone PCR. By this approach we have increased the half-life of denaturation by 300-fold compared to the wild-type Bacillus subtilis lipase through three point mutations generated by only two cycles of error-prone PCR. At lower temperatures the activity parameters of the thermostable mutants are unaltered. High-resolution crystal structures of the mutants show subtle changes, which include stacking of tyrosine residues, peptide plane flipping and a better anchoring of the terminus, that challenge rational design and explain the structural basis for enhanced thermostability. The approach may offer an efficient and minimalist solution for the enhancement of a desired property of a protein.     

Cloning and structure of the hem A gene of Escherichia coli K-12

An Escherichia coli gene, which complements two independent hemA mutants of E. coli, has been cloned onto a multi-copy plasmid and both its strands have been sequenced. Both complemented mutants produce 5-aminolevulinic acid (ALA) and display fluorescence after 24h. The cloned sequence appears to encode a 46-kDa protein, which when produced in the maxicell procedure is processed to a 41-kDa protein as determined by sodium dodecyl sulfate-polyacrylamide-gel electrophoresis. The amino acid sequence of the cloned gene product shows no significant homologies with any cloned ALA synthase, nor with any protein, in two E. coli databanks. A second cloned gene fragment, which has its coding region 34 bp away from the coding region of the gene that complements hemA, has been identified as part of protein release factor 1(RF1), thus confirming the location of hemA at min 26.7 and mapping it precisely near RF1. We have shown that E. coli utilizes the intact five-carbon chain of glutamate for the synthesis of ALA [Li et al., J Bacteriol. 171 (1989b) 2547-2552].     

Structural basis of thermostability. Analysis of stabilizing mutations in subtilisin BPN'

The crystal structures of two thermally stabilized subtilisin BPN' variants, S63 and S88, are reported here at 1.8 and 1.9 A resolution, respectively. The micromolar affinity calcium binding site (site A) has been deleted (Delta75-83) in these variants, enabling the activity and thermostability measurements in chelating conditions. Each of the variants includes mutations known previously to increase the thermostability of calcium-independent subtilisin in addition to new stabilizing mutations. S63 has eight amino acid replacements: D41A, M50F, A73L, Q206W, Y217K, N218S, S221C, and Q271E. S63 has 75-fold greater stability than wild type subtilisin in chelating conditions (10 mm EDTA). The other variant, S88, has ten site-specific changes: Q2K, S3C, P5S, K43N, M50F, A73L, Q206C, Y217K, N218S, and Q271E. The two new cysteines form a disulfide bond, and S88 has 1000 times greater stability than wild type subtilisin in chelating conditions. Comparisons of the two new crystal structures (S63 in space group P2(1) with A cell constants 41.2, 78.1, 36.7, and beta = 114.6 degrees and S88 in space group P2(1)2(1)2(1) with cell constants 54.2, 60.4, and 82.7) with previous structures of subtilisin BPN' reveal that the principal changes are in the N-terminal region. The structural bases of the stabilization effects of the new mutations Q2K, S3C, P5S, D41A, Q206C, and Q206W are generally apparent. The effects are attributed to the new disulfide cross-link and to improved hydrophobic packing, new hydrogen bonds, and other rearrangements in the N-terminal region.     

Mutational analysis of the bacteriophage phi X174 replication origin

Bacteriophage phi X174 mutants within the 30 base-pair replication origin were constructed using oligodeoxynucleotide-directed mutagenesis. A total of 18 viable base substitution mutants at 13 different positions within the origin region were obtained. The majority of these ori mutants have a plaque morphology and burst size comparable to that of wild-type phi X174. Two phi X174 ori mutants with a reduced growth ability spontaneously acquired additional mutations that enhanced the growth rate. The additional mutation was located at the same site as the original mutation or was located in the N-terminal part of the gene A protein. This latter secondary mutation is responsible for a better binding and/or recognition of the gene A protein to the mutated origin. In a Darwinian experiment wild-type phi X174 outgrows all phi X174 ori mutants, indicating the superiority of the wild-type ori sequence for the reproduction of bacteriophage phi 174. Insertions and deletions were constructed at different positions within the phi X174 replication origin cloned in a plasmid. Small insertions and deletions in the A + T-rich spacer region do not inhibit phi X174 gene A protein cleavage in vitro, but severely impair packaging of single-stranded plasmid DNA in viral coats.     

Cover Image,Volume 116, Number 8, August 2019

Nattokinase (NK) is a serine protease of the subtilisin family; as a potent fibrinolytic enzyme, it is potentially useful for thrombosis therapy. For NK to be applied as an oral medicine for the treatment of cardiovascular diseases, it must overcome the extremely acidic environments of the gastrointestinal tract despite its limited acidic stability. In this study, three strategies were adopted to improve the acid resistance of NK: (a) Surface charge engineering, (b) sequence alignment, and (c) mutation based on the literature. Eleven variants were constructed and four single-point mutations were screened out for their distinctive catalytic properties: Q59E increased the specific activity; S78T improved the acid stability; Y217K enhanced the acid and thermal stabilities; and N218D improved the thermostability. Based on these observations, multipoint variants were constructed and characterized, and one variant with better acid stability, catalytic efficiency, and thermostability was obtained. Molecular dynamics simulation was carried out to clarify the molecular mechanism of the increased stability of S78T and Y217K mutants under acidic conditions. This study explored effective strategies to engineer acid resistance of NK; moreover, the NK variants with better catalytic properties found in this study have potential applications for the medical industry.     

Identification, physical map location and sequence of the denV gene from bacteriophage T4

The denV gene from bacteriophage T4, which codes for endonuclease V, a small DNA repair enzyme, has been cloned and identified by an approach combining DNA sequencing and genetics, independent of the phenotypic effect of the cloned gene. Appropriate DenV+ and DenV- deletion mutants were mapped physically to define precisely a region encompassing the denV gene. This region was sequenced in order to identify a protein-coding sequence of the correct size for the denV gene (400-500 bp). Finally, identification was confirmed by sequencing the corresponding fragments cloned from four genetically and phenotypically well-characterized denV mutants. The denV gene is located at 64 kb on the T4 genome, adjacent to the ipII gene, and codes for a basic protein of 138 amino acids with a deduced molecular weight of 16,078.     

Thermostable Bacillus subtilis lipases: in vitro evolution and structural insight

In vitro evolution methods are now being routinely used to identify protein variants with novel and enhanced properties that are difficult to achieve using rational design. However, one of the limitations is in screening for beneficial mutants through several generations due to the occurrence of neutral/negative mutations occurring in the background of positive ones. While evolving a lipase in vitro from mesophilic Bacillus subtilis to generate thermostable variants, we have designed protocols that combine stringent three-tier testing, sequencing and stability assessments on the protein at the end of each generation. This strategy resulted in a total of six stabilizing mutations in just two generations with three mutations per generation. Each of the six mutants when evaluated individually contributed additively to thermostability. A combination of all of them resulted in the best variant that shows a remarkable 15 °C shift in melting temperature and a millionfold decrease in the thermal inactivation rate with only a marginal increase of 3 kcal mol⁻¹ in free energy of stabilization. Notably, in addition to the dramatic shift in optimum temperature by 20 °C, the activity has increased two- to fivefold in the temperature range 25-65 °C. High-resolution crystal structures of three of the mutants, each with 5° increments in melting temperature, reveal the structural basis of these mutations in attaining higher thermostability. The structures highlight the importance of water-mediated ionic networks on the protein surface in imparting thermostability. Saturation mutagenesis at each of the six positions did not result in enhanced thermostability in almost all the cases, confirming the crucial role played by each mutation as revealed through the structural study. Overall, our study presents an efficient strategy that can be employed in directed evolution approaches employed for obtaining improved properties of proteins.