Protein engineering of the high-alkaline serine protease PB92 from Bacillus alcalophilus: functional and structural consequences of mutation at the S4 substrate binding pocket

Serine endoproteases such as trypsins and subtilisins are knownto have an extended substrate binding region that interactswith residues P6 to P3' of a substrate. In order to investigatethe structural and functional effects of replacing residuesat the S4 substrate binding pocket, the serine protease fromthe alkalophilic Bacillus strain PB92, which shows homologywith the subtilisins, was mutated at positions 102 and 126–128.Substitution of Val102 by Trp results in a 12–fold increasein activity towards succinyl-L-Ala-L-Ala-L-Pro-L-Phe-p-nitroanilide(sAAPFpNA). An X-ray structure analysis of the V102W mutantshows that the Trp side chain occupies a hydrophobic pocketat the surface of the molecule leaving a narrow crevice forthe P4 residue of a substrate. Better binding of sAAPFpNA bythe mutant compared with the wild type protein as indicatedby the kinetic data might be due to the hydrophobic interactionof Ala P4 of the substrate with the introduced Trp102 side chain.The observed difference in binding of sAAPFpNA by protease PB92and thermitase, both of which possess a Trp at position 102,is probably related to the amino acid substitutions at positions105 and 126 (in the protease PB92 numbering).Kinetic data forthe variants obtained by random mutation of residues Serl26,Prol27 and Serl28 reveal that the activity towards sAAPFpNAincreases when a hydrophobic residue is introduced at position126. An X-ray diffraction analysis was carried out for the threeprotease PB92 mutants which have residues Serl26-Prol27-Serl28replaced by Met-Ala-Gly(‘MAG’ mutant), Phe-Gln-Ser(‘FQS’ mutant) and Asn-Ser-Ala (‘NSA’mutant). Met 126 and Phel26 in the crystal structures of thecorresponding mutants are fixed in the same hydrophobic environmentas Trp102 in the V102W mutant.In contrast, Asnl26 in the ‘NSA’mutant is completely disordered in both crystal forms for whichthe structure has been determined. According to our kineticmeasurements none of the mutants with Met, Phe, Leu or Val atposition 126 binds sAAPFpNA better than the wild type enzyme.Resultsof the site-directed mutagenesis at position 127 imply thatpossible interaction of this residue with a substrate has almostno effect on activity towards sAAPFpNA and casein.     

The refined crystal structure of subtilisin Carlsberg at 2.5 ? resolution

We report here the X-ray crystal structure of native subtilisinCarlsberg, solved at 2.5 ? resolution by molecular replacementand refined by restrained least squares to a crystal-lographicresidual of 0.206. we compare this structure to the crystal structure of subtilisin BPN'. We find that, despite82 amino acid substitutions and one deletion in subtilisin Carlsbergrelative to subtilisin BPN', the structures of these enzymesare remarkably similar. We calculate an r.m.s. difference betweenequivalent a-carbon positions in subtilisin Carlsberg and subtilisinBPN' of only 0.55 ?. This confirms previous reports of extensivestructural bomology between these two subtilisins based on X-raycrystal structures of the complex of eglin-c with subtilisinCarlsberg [McPhalen,C.A., Schnebli.H.P. and James,M.N.G. (1985)FEBS Lett., 188, 55; Bode,W., Papamokos,E. and Musil,D. (1987)Eur. J. Biochem., 166, 673-692]. In addition, we find that thenative active sites of subtilisins Carlsberg and BPN' are virtuallyidentical. While conservative substitutions at residues 217and 156 may have subtle effects on the environments of substrate-bindingsites SI' and SI respectively, we find no obvious structuralcorrelate for reports that subtilisins Carlsberg and BPN' differin their recognition of model substrates. In particular, wefind no evidence that the hydrophobic binding pocket SI in subtilisinCarlsberg is ‘deeper’, ‘narrower’ or'less polar' than the corresponding binding site hi subtilisinBPN' [Karasaki and Ohno (1978) J. Biochem., Tokyo, 84, 531–538].     

Molecular mechanics analysis of inhibitor binding to HIV-1 protease

Crystallographic structures of HIV protease with three differentpeptide-mimetic inhibitors were subjected to energy minimizationusing molecular mechanics, the minimized structures analyzedand the inhibitor binding energies calculated. Partial chargeassignment for the hydrogen bonded catalytic aspartk acids,Asp25 and -25', was in good agreement with charge calculationsusing semi-empirical molecular orbital methods. Root mean squaredeviations on minimization were small and similar for both subunitsin the protease dimer. The surface loops, which had the largestB factors, changed most on minimization; the hydrophobic coreand the inhibitor binding site showed little change. The distance-dependentdielectric of D(r) = 4r was found to be preferable to D(r) =r. Distance restraints were applied for the intermolecular hydrogenbonds to maintain the conformation of the inhibitor bindingsite. Using the dielectric of D(r) = 4r, the calculated interactionenergy of the three inhibitors with the protease ranged from–53 to –56 kcal/mol. The groups of the inhibitorswere changed to add or remove a ‘transition state analogue’hydroxyl group, and the loss in energy on the removal of thisgroup was calculated to be 0.9–1.7 kcal/mol. This wouldrepresent 19–36% of the total measured difference in bindingenergy between the inhibitors JG365 and MVT-101.     

A hybrid of bovine pancreatic ribonuclease and human angiogenin: an external loop as a module controlling substrate specificity?

A comparison of the sequences of three homologous ribonucleases(RNase A, angiogenin and bovine seminal RNase) identifies threesurface loops that are highly variable between the three proteins.Two hypotheses were contrasted: (i) that this variation mightbe responsible for the different catalytic activities of thethree proteins; and (ii) that this variation is simply an exampleof surface loops undergoing rapid neutral divergence in sequence.Three hybrids of angiogenin and bovine pancreatic ribonuclease(RNase) A were prepared where regions in these loops taken fromangiogenin were inserted into RNase A. Two of the three hybridshad unremarkable catalytic properties. However, the RNase Amutant containing residues 63–74 of angiogenin had greatlydiminished catalytic activity against uridylyl-(3' – 5')-adenosine(UpA), and slightly increased catalytic activity as an inhibitorof translation in vitro. Both catalytic behaviors are characteristicof angiogenin. This is one of the first examples of an engineeredexternal loop in a protein. Further, these results are complementaryto those recently obtained from the complementary experiment,where residues 59–70 of RNase were inserted into angiogenin[Harper and Vallee (1989) Biochemistry, 28, 1875–1884].Thus, the external loop in residues 63–74 of RNase A appearsto behave, at least in part, as an interchangeable ‘module’that influences substrate specificity in an enzyme in a waythat is isolated from the influences of other regions in theprotein.     

A structure-function relationship for the calcium affinities of regulatory proteins containing 'EF-hand' pairs

Using a series of homologous calcium-binding proteins, a quantitativestructure–activity relationship (QSAR), log(1/K_d) = –18.986– 1.6278(X₁) + 0.7981(X₂) + 0.2312(X₃), has been established,which relates the calcium-binding affinities (1/K_d) of the regulatoryproteins with (i) the net ligand charge (X₁) of the two calciumbinding loops, (ii) the hydrophobicity (X₂) of the ß-sheetsegment of the loops and (iii) the hydrophobicity (X₃) of thefour ‘EF-hand’ helices. It is found that the bindingaffinities are influenced by the ‘EF-hand’ pairrather than the individual ‘EF-hands’. The QSAR,in addition to explaining satisfactorily the large variationin the observed calcium affinities, can predict the affinitiesof the ‘EF-hand’ pairs in other proteins from theamino acid sequence and can also account for the changes inthe affinities caused by substitution in the hydrophobic and/ormetal-coordinating residues. Thus, this relationship can beemployed in protein design and engineering. The method is potentiallyuseful in the development of similar relationships for the bindingof other proteins to substrates, inhibitors, drugs and co-factors.     

Molecular modeling of coiled-coil {alpha}-tropomyosin: analysis of staggered and in register helix--helix interactions

In register and staggered models of tropomyosin coiled-coilwere built from X-ray C coordinates and refined via moleculardynamics. The two models show similar structural features withthe X-ray structure of GCN4 leucine zipper. Empirical energeticmethods used to compare the in register and staggered modelsindicate that both are equally probable. The two models havesimilar profiles of solvation free energy of folding for residuesat positions a and d of the repeating heptad, indicating thatresidues at these positions are as well buried in an in registerstructure as in a staggered one. Neither the in register northe 14 residues staggered structure can be ruled out based onhydrophobic or e–g' (g–e') electrostatic interactionswhich are not able to distinguish between the two models andare therefore not selective. However, the eg–b'c' electrostaticinteractions, although smaller in magnitude, are in favor ofthe in register model. Furthermore, analysis of hydrophobicand electrostatic interactions along the tropomyosin sequenceshows that bulky residues in positions a and d prevent the formationof inter-chain salt bridges.     

Making tissue-type plasminogen activator more fibrin specific

The fibrin specificity of tissue-type plasminogen activatorcan be increased by mutagenesis within at least four sites inthe protease domain. These sites include residue I276, the newN-terminus formed by conversion to a two-chain structure, residueson either side of the active site cleft, KHRR 296–299or DDD 364–366, a charged surface involved in fibrin interactions,which includes residues H432, R434, D460, R462 and a loop structure,PQANL 466–470, near the fibrin-binding patch. Variantswith mutations at any of these sites have low fibrinogen-stimulatedactivity, whereas fibrin-stimulated activity is at least normal.Kinetic analysis reveals that mutations at these positions reducethe k_cat in the presence of fibrinogen, but leave the moleculeswith normal kinetic constants in the presence of fibrin. A significantexception is found at positions 296–299, where the presenceof fibrin manifests significant increases in both k_cat and K_m.Combinations of mutations at these sites appear to be additivewith respect to fibrin specificity.     

Cumulative stabilizing effects of glycine to alanine substitutions in Bacillus subtilis neutral protease

Oligonucleotide-directed mutagenesis has been used to replaceglycine residues by alanine in neutral protease from Bacillussubtilis. One Gly to Ala substitution (G147A) was located ina helical region of the protein, while the other (G189A) wasin a loop. The effects of mutational substitutions on the functional,conformational and stability properties of the enzyme have beeninvestigated using enzymatic assays and spectroscopic measurements.Single substitutions of both G1y147 and Gly189 with Ala residuesaffect the enzyme kinetic properties using synthetic peptidesas substrates. When Gly replacements were concurrently introducedat both positions, the kinetic characteristics of the doublemutant were roughly intermediate between those of the two singlemutants, and similar to those of the wild-type protease. Bothmutants G147A and G189A were found to be more stable towardsirreversible thermal inactivation/unfolding than the wild-typespecies. Moreover, the stabilizing effect of the Gly to Alasubstitution was roughly additive in the double mutant G147A/G189A,which shows a 3.2°C increase in T_m with respect to the wild-typeprotein. These findings indicate that the Gly to Ala substitutioncan be used as a strategy to stabilize globular proteins. Thepossible mechanisms of protein stabilization are also discussed.     

Conformational modeling of substrate binding to endocellulase E2 from Thermomonospora fusca

Molecular mechanics calculations have been used to place a cellotetraosesubstrate into the active site of the crystallographlcally determinedstructure of endocellulase E2 from Thermomonospora fusca. Inthe lowest energy model structure, the second residue of thesubstrate oligosaccharide is tilted away from the planar ribbongeometry of cellulose as it is in the X–ray structureof the E2_cd–cellobiose cocrystal. This tilt is the resultof the topology of the binding site, and results in severalstrong carbohydrate–protein hydrogen bonds. The tiltingproduces a twisting of the glycosidic linkage of the cleavagesite between residues two and three. In the predicted enzyme–substratecomplex both of the Asp residues believed to function in generalacid and base roles in the previously proposed model for themechanism are distant from the bond being cleaved. Moleculardynamics simulations of the complex were conducted, and whilethe putative catalytic Asp residues remained distant from thecleavage site, the proton of Tyr73 briefly came within van derWaals contact of the linkage oxygen.     

Knowledge-based modeling of the serine protease triad into non-proteases

The Asp–His–Ser triad of serine proteases has beenregarded, in the present study, as an independent catalyticmotif, because in nature it has been incorporated at the activesites of enzymes as diverse as the serine proteases and thelipases. Incorporating this motif into non-protease scaffolds,by rational design and mutagenesis, might lead to the generationof novel catalysts. As an aid to such experiments, a knowledge-basedcomputer modeling procedure has been developed to model theprotease Asp–His–Ser triad into non-proteases. Catalytictriads from a set of trypsin family proteases have been analyzedand criteria that characterize the geometry of the triads havebeen obtained. Using these criteria, the modeling procedurefirst identifies sites in non-proteases that are suitable formodeling the protease triad. H-bonded Asp–His–Sertriads, that mimic the protease catalytic triad in geometry,are then modeled in at these sites, provided it is stereochemicallypossible to do so. Thus non-protease sites at which H-bondedAsp–His–Ser triads are successfully modeled in maybe considered for mutagenesis experiments that aim at introducingthe protease triad into non-proteases. The triad modeling procedurehas been used to identify sites for introducing the proteasetriad in three binding proteins and an immunoglobulin. A scoringfunction, depending on inter-residue distances, solvent accessibilityand the substitution potential of amino acid residues at themodeling sites in the host proteins, has been used to assessthe quality of the model triads.     

The effect of cavity-filling mutations on the thermostability of Bacillus stearothermophilus neutral protease

Cavities in the hydrophobic core of the neutral protease ofBacillus stearothermophilus were analyzed using a threedimensionalmodel that was inferred from the crystal structure of thermolysin,the highly homologous neutral protease of B.thermoproteolyticus(85% sequence identity). Site–directed mutagenesis wasused to fill some of these cavities, thereby improving hydrophobicpacking in the protein interior. The mutations had small effectson the thermostability, even after drastic changes, such asLeu284Trp and Met168Trp. The effects on T50, the temperatureat which 50% of the enzyme is irreversibly inactivated in 30min, ranged from 0.0 to +0.4°C. These results can be explainedby assuming that the mutations have positive and negative structuraleffects of approximately the same magnitude. Alternatively,it could be envisaged that the local unfolding steps, whichrender the enzyme susceptible towards autolysis and which arerate limiting in the process of thermal inactivation, are onlyslightly affected by alterations in the hydrophobic core.     

Prediction of the antigenic sites of the cystic fibrosis transmembrane conductance regulator protein by molecular modelling

Antibodies are powerful tools for studying the in situ localizationand physiology of proteins. The prediction of epitopes by molecularmodelling has been used successfully for the papilloma virus,and valuable antibodies have been raised [Muller et aL (1990)J. Gen. Virol, 71, 2709–2717]. We have improved the modellingapproach to allow us to predict epitopes from the primary sequencesof the cystic fibrosis transmembrane conductance regulator.The procedure involves searching for fragments of primary sequenceslikely to make amphipathic secondary structures, which are hydrophilicenough to be at the surface of the folded protein and thus accessibleto antibodies. Amphipathic helices were predicted using themethods of Berzofsky, Eisenberg and Jahnig. Their hydrophobichydrophilicinterface was calculated and drawn, and used to predict theorientation of the helices at the surface of the native protein.Amino acids involved in turns were selected using the algorithmof Eisenberg. Tertiary structures were calculated using ‘FOLDING’,a software developed by R.Brasseur for the prediction of smallprotein structures [Brasseur (1995) J. MoL Graphics, in press].We selected sequences that folded as turns with at least fiveprotruding polar residues. One important property of antibodiesis selectivity. To optimize the selectivity of the raised antibodies,each sequence was screened for similarity (FASTA) to the proteinsequences from several databanks. Ubiquitous sequences werediscarded. This approach led to the identification of 13 potentialepitopes in the cystic fibrosis transmembrane conductance regulator:seven helices and six loops.     

The cystine knot of a squash-type protease inhibitor as a structural scaffold for Escherichia coli cell surface display of conformationally constrained peptides

The Ecballium elaterium trypsin inhibitor II (EETI-II), a memberof the squash family of protease inhibitors, is composed of28 amino acid residues and is a potent inhibitor of trypsin.Its compact structure is defined by a triple-stranded antiparallelß-sheet, which is held together by three intramoleculardisulfide bonds forming a cystine knot. In order to explorethe potential of the EETI-II peptide to serve as a structuralscaffold for the presentation of randomized oligopeptides, weconstructed two EETI-II derivatives, where the six-residue inhibitorloop was replaced by a 13-residue epitope of Sendai virus L-proteinand by a 17-residue epitope from human bone Gla-protein. EETI-IIand derived variants were produced via fusion to maltose bindingprotein MalE. By secretion of the fusion into the periplasmicspace, fully oxidized and correctly folded EETI-II was obtainedin high yield. EETI-II and derived variants could be presentedon the Escherichia coli outer membrane by fusion to truncatedLpp'–OmpA', which comprises the first nine residues ofmature lipoprotein plus the membrane spanning ß-strandfrom residues 46–66 of OmpA protein. Gene expression wasunder control of the strong and tightly regulated tetA promoter/operator.Cell viability was found to be drastically reduced by high levelexpression of Lpp'–OmpA'–EETI-II fusion protein.To restore cell viability, net accumulation of fusion proteinin the outer membrane was reduced to a tolerable level by introductionof an amber codon at position 9 of the lpp' sequence and utilizingan amber suppressor strain as expression host. Cells expressingEETI-II variants containing an epitope were shown to be surfacelabeled with the respective monoclonal antibody by indirectimmunofluorescence corroborating the cell surface exposure ofthe epitope sequences embedded in the EETI-II cystine knot scaffold.Cells displaying a particular epitope sequence could be enriched10⁷-fold by combining magnetic cell sorting with fluorescence-activatedcell sorting. These results demonstrate that E.coli cell surfacedisplay of conformationally constrained peptides tethered tothe EETI-II cystine knot scaffold has the potential to becomean effective technique for the rapid isolation of small peptidemolecules from combinatorial libraries that bind with high affinityto acceptor molecules.     

Determination of the pKa values of titratable groups of an antigen-antibody complex, HyHEL-5-hen egg lysozyme

The titration behavior of the ionizable residues of the HyHEL-5–henegg lysozyme complex and its individual components has beenstudied using continuum electrostatic calculations. Severalresidues of HyHEL-5 had pK_a values shifted away from model valuesfor isolated residues by more than three pH units. Shifts awayfrom the model values were smaller for the residues of hen egglysozyme. A moderate variation in the pK_a values of the titratablegroups was observed upon increase of the ionic strength from0 to 100 mM, amounting to 1–2 pH units in most cases.Under physiological conditions, the net charge of HyHEL-5 wasopposite that for hen egg lysozyme. Several residues, includingthose involved in the Arg–Glu salt bridges that have beenproposed to be important in antibody-antigen binding, had pK_avalues that were changed significantly upon binding. The maintitration event upon antibody-antigen binding appears to beloss of a proton from residue GluH50 of the Fv molecule. Thelimitations of our calculation methods and the role they mightplay in the design of antibodies for use in assays, sensorsand separations are discussed     

Studies of the pore-forming domain of a voltage-gated potassium channel protein

Recent mutagenesis studies nave identified a stretch of aminoacid residues which form the ion-selective pore of the voltage-gatedpotassium channel. It has been suggested that this sequenceof amino acids forms a ß-barrel structure making upthe structure of the ion-selective pore [Hartman,H.A., Kirsch,G.E.,DreweJ.A., Taglialatela.M., Joho.R.H. and Brown,A.M. (1991)Science, 251, 942–944; YeUen.G., Jurman,M.E., Abramson,T.and MacKinnon,R. (1991) Science, 251, 939–942; Yool,AJ.and Schwarz.T.L. (1991) Nature, 349, 700–704]. We havesynthesized a polypeptide corresponding to this amino add sequence(residues 431–449 of the ShA potassium channel from Drosophila).A tetrameric version of this sequence was also synthesized byUnking together four of these peptldes onto a branching lysinecore. Fourier transform infrared (FT-LR) and circular dichroism(CD) spectroscopy have been used to investigate the structureof these peptides after their reconstitution into lyso phos-phatidylcholinemicelles and lipid bilayers composed of dimyristoyl phosphatidyfcholineand dimyristoyl phosphatidyl-glycerol. The spectroscopic studiesshow that these peptides are predominantly a-helical in theselipid environments. When Incorporated into planar lipid bilayersboth peptides induce ion channel activity. Molecular modellingstudies based upon the propensity of these peptides to forman -helical secondary structure in a hydrophobfc environmentare described. These results are discussed in the light of recentmutagenesis and binding studies of the Drosophila Shaker potassiumion channel protein     

A study of structural determinants in the interleukin-1 fold

The structures of interleukin-1ß, basic fibroblastgrowth factor and Erythrina trypsin inhibitor have been analysedin order to determine whether the hydrophobic core remains conserved,even when the structures have extremely low sequence similarities.We find that there are significant differences in the way eachprotein achieves a satisfactory arrangement of core residuesand that positions which contribute to the core of one structureare not guaranteed to contribute to the integrity of another.Furthermore, the side-chain packing arrangements of these coreresidues vary significantly between the three structures. Duringthis analysis the side-chain rotamers for three independentlydetermined interleukin-1ß structures were also compared.It was found that although buried residues are generally inagreement the remaining residues frequently occupy differentrotamers in the three structures. This suggests that althoughmeaningful studies are possible for buried side-chains the resultsobtained from equivalent analyses of accessible residues shouldbe treated with caution. These results are discussed with specificreference to the optimization of side-chain packing in proteinsof known structure.     

Molecular dynamics model structures for the molten globule state of {alpha}-lactalbumin: aromatic residue clusters I and II

To model the molten globule structure of -lactalbumin, moleculardynamics (MD) simulations were carried out for the protein inexplicit water at high temperature. In these simulations, long-rangeCoulomb interactions were evaluated explicitly with an originalmethod (particle–particle and particle–cell: PPPC)to avoid artifacts caused by the cut-off. The MD simulationswere started from two initial conditions to verify that similarresults would be obtained. From the last 150 ps trajectoriesof the two MD simulations, two partially unfolded average structureswere obtained. These structures had the following common structuralfeatures which are characteristic of the molten globule state.The radii of gyration for these conformations were 7.4 and 9.6%larger than that of the native state. These values were almostthe same as the experimental value (9.6%) observed recentlyby small-angle X-ray scattering (Kataoka,M., Kuwajima,K., Tokunaga,F.and Goto,Y., 1997, Protein Sci., 6, 422–430). Furthermore,aromatic residues of clusters I and II in these structures werefar apart from each other except for Try103–Trp104. Thisresult is in good agreement with NMR experimental results forthe acid-denatured molten globule state (Alexandrescu et al.,1992, 1993); that is, NOE signals between the aromatic residueswere not observed, except for that of Try103–Trp104 inthe molten globule state. Other structural features of thesemodels for the molten globule state are discussed with referenceto native state structures.     

The {alpha}/{beta} hydrolase fold

We have identified a new protein fold—the /ßhydrolase fold—that is common to several hydrolytic enzymesof widely differing phylogenetic origin and catalytic function.The core of each enzyme is similar: an /ß sheet, notbarrel, of eight ß-sheets connected by -helices. Theseenzymes have diverged from a common ancestor so as to preservethe arrangement of the catalytic residues, not the binding site.They all have a catalytic triad, the elements of which are borneon loops which are the best-conserved structural features inthe fold. Only the histidine in the nucleophile-histidine-acidcatalytic triad is completely conserved, with the nucleophileand acid loops accommodating more than one type of amino acid.The unique topological and sequence arrangement of the triadresidues produces a catalytic triad which is, in a sense, amirror-image of the serine protease catalytic triad. There arenow four groups of enzymes which contain catalytic triads andwhich are related by convergent evolution towards a stable,useful active site: the eukaryotic serine proteases, the cysteineproteases, subtilisins and the /ß hydrolase fold enzymes.     

Increasing the thermostability of a neutral protease by replacing positively charged amino acids in the N-terminal turn of {alpha}-helices

The 247–260 and 289–299 -helices of Bacillus subtilisneutral protease have a lysine in their N-terminal turn. Theselysines were replaced by Ser or Asp in order to improve electrostaticinteractions with the -helix dipole. After replacing Lys bySer at positions 249 or 290, the thermostability of the enzymewas increased by 0.3 and 1.0°C, respectively. The Asp249and Asp290 mutants exhibited a stabilization of 0.6 and 1.2°C,respectively. The results show the feasibility of stabilizingenzymes by introducing favourable residues at the end of -helices.     

Introduction of a stabilizing 10 residue {beta}-hairpin in Bacillus subtilis neutral protease

A 10 residue ß-hairpin, which is characteristic ofthermostable Bacillus neutral proteases, was engineered intothe thermolabile neutral protease of Bacillus subtilis. Therecipient enzyme remained fully active after introduction ofthe loop. However, the mutant protein exhibited autocatalyticnicking and a 0.4°C decrease in thermostability. Two additionalpoint mutations designed to improve the interactions betweenthe enzyme surface and the introduced ß-hairpin resultedin reduced nicking and increased thermostability. After theintroduction of both additional mutations in the loopcontainingmutant, nicking was largely prevented and an increase in thermostabilityof 1.1°C was achieved.