Construction of stabilized proteins by combinatorial consensus mutagenesis

We constructed stabilized variants of beta-lactamase (BLA) from Enterobacter cloacae by combinatorial recruitment of consensus mutations. By aligning the sequences of 38 BLA homologs, we identified 29 positions where the E.cloacae gene differs from the consensus sequence of lactamases and constructed combinatorial libraries using mixtures of mutagenic oligonucleotides encompassing all 29 positions. Screening of 90 random isolates from these libraries identified 15 variants with significantly increased thermostability. The stability of these isolates suggest that all tested mutations make additive contributions to protein stability. A statistical analysis of sequence and stability data identified 11 mutations that made stabilizing contributions and eight mutations that destabilized the protein. A second-generation library recombining these 11 stabilizing mutations led to the identification of BLA variants that showed further stabilization. The most stable variant had a mid-point of thermal denaturation (Tm) that was 9.1 degrees C higher than the starting molecule and contained eight consensus mutations. Incubation of three stabilized BLA variants with several proteases showed that all tested isolates have significantly increased resistance to proteolysis. Our data demonstrate that combinatorial consensus mutagenesis (CCM) allows the rapid generation of protein variants with improved thermal and proteolytic stability.     

Limitations of yeast surface display in engineering proteins of high thermostability

Engineering proteins that can fold to unique structures remains a challenge. Protein stability has previously been engineered via the observed correlation between thermal stability and eukaryotic secretion level. To explore the limits of an expression-based approach, variants of the highly thermostable three-helix bundle protein alpha3D were studied using yeast surface display. A library of alpha3D mutants was created to explore the possible correlation of protein stability and fold with expression level. Five efficiently expressed mutants were then purified and further studied biochemically. Despite their differences in stability, most mutants expressed at levels comparable with that of wild-type alpha3D. Two other related sequences (alpha3A and alpha3B) that form collapsed, stable molten globules but lack a uniquely folded structure were similarly expressed at high levels by yeast display. Together these observations suggest that the quality control system in yeast is unable to discriminate between well-folded proteins of high stability and molten globules. The present study, therefore, suggests that an optimization of the surface display efficiency on yeast may yield proteins that are thermally and chemically stable yet are poorly folded.     

Enhancing the thermal stability of mitochondrial cytochrome b5 by introducing a structural motif characteristic of the less stable microsomal isoform

Outer mitochondrial membrane cytochrome b5 (OM b5) is the most thermostable cytochrome b5 isoform presently known. Herein, we show that OM b5 thermal stability is substantially enhanced by swapping an apparently invariant motif in its heme-independent folding core with the corresponding motif characteristic of its less stable evolutionary relative, microsomal cytochrome b5 (Mc b5). The motif swap involved replacing two residues, Arg15 with His and Glu20 with Ser, thereby introducing a Glu11-His15-Ser20 H-bonding triad on the protein surface along with a His15/Trp22 pi-stacking interaction. The ferric and ferrous forms of the OM b5 R15H/E20S double mutant have thermal denaturation midpoints (Tm values) of approximately 93 degrees C and approximately 104 degrees C, respectively. A 15 degrees C increase in apoprotein Tm plays a key role in the holoprotein thermal stability enhancement, and is achieved by one of the most common natural mechanisms for stabilization of thermophilic versus mesophilic proteins: raising the unfolding free energy along the entire stability curve.     

Generation and analysis of proline mutants in protein G

The pyrrolidine ring of the amino acid proline reduces the conformational freedom of the protein backbone in its unfolded form and thus enhances protein stability. The strategy of inserting proline into regions of the protein where it does not perturb the structure has been utilized to stabilize many different proteins including enzymes. However, most of these efforts have been based on trial and error, rather than rational design. Here, we try to understand proline's effect on protein stability by introducing proline mutations into various regions of the B1 domain of Streptococcal protein G. We also applied the Optimization of Rotamers By Iterative Techniques computational protein design program, using two different solvation models, to determine the extent to which it could predict the stabilizing and destabilizing effects of prolines. Use of a surface area dependent solvation model resulted in a modest correlation between the experimental free energy of folding and computed energies; on the other hand, use of a Gaussian solvent exclusion model led to significant positive correlation. Including a backbone conformational entropy term to the computational energies increases the statistical significance of the correlation between the experimental stabilities and both solvation models.     

Interactions between immunoglobulin-like and catalytic modules in Clostridium thermocellum cellulosomal cellobiohydrolase CbhA

Cellobiohydrolase CbhA from Clostridium thermocellum cellulosome is a multi-modular protein composed starting from the N-terminus of a carbohydrate-binding module (CBM) of family 4, an immunoglobulin(Ig)-like module, a catalytic module of family 9 glycoside hydrolases (GH9), X1(1) and X1(2) modules, a CBM of family 3 and a dockerin module. Deletion of the Ig-like module from the Ig-GH9 construct results in complete inactivation of the GH9 module. The crystal structure of the Ig-GH9 module pair reveals the existence of an extensive module interface composed of over 40 amino acid residues of both modules and maintained through a large number of hydrophilic and hydrophobic interactions. To investigate the importance of these interactions between the two modules, we compared the secondary and tertiary structures and thermostabilities of the individual Ig-like and GH9 modules and the Ig-GH9 module pair using both circular dichroism (CD) spectroscopy and differential scanning calorimetry (DSC). Thr230, Asp262 and Asp264 of the Ig-like module are located in the module interface of the Ig-GH9 module pair and are suggested to be important in 'communication' between the modules. These residues were mutated to alanyl residues. The structure, stability and catalytic properties of the native Ig-GH9 and its D264A and T230A/D262A mutants were compared. The results indicate that despite being able to fold relatively independently, the Ig-like and GH9 modules interact and these interactions affect the final fold and stability of each module. Mutations of one or two amino acid residues lead to destabilization and change of the mechanism of thermal unfolding of the polypeptides. The enzymatic properties of native Ig-GH9, D264A and T230A/D262A mutants are similar. The results indicate that inactivation of the GH9 module occurs as a result of multiple structural disturbances finally affecting the topology of the catalytic center.     

Revisiting the correlation between proteins' thermoresistance and organisms' thermophilicity

The possibility to rationally design protein mutants that remain structured and active at high temperatures strongly depends on a better understanding of the mechanisms of protein thermostability. Studies devoted to this issue often rely on the living temperature (T(env)) of the host organism rather than on the melting temperature (T(m)) of the analyzed protein. To investigate the scale of this approximation, we probed the relationship between T(m) and T(env) on a dataset of 127 proteins, and found a much weaker correlation than previously expected: the correlation coefficient is equal to 0.59 and the regression line is T(m) approximately 42.9 degrees C + 0.62T(env). To illustrate the effect of using T(env) rather than T(m) to analyze protein thermoresistance, we derive statistical distance potentials, describing Glu-Arg and Asp-Arg salt bridges, from protein structure sets with high or low T(m) or T(env). The results show that the more favorable nature of salt bridges, relative to other interactions, at high temperatures is more clear-cut when defining thermoresistance in terms of T(m). The T(env)-based sets nevertheless remain informative.     

Engineering a chimeric pyrroloquinoline quinone glucose dehydrogenase: improvement of EDTA tolerance, thermal stability and substrate specificity

An engineered Escherichia coli PQQ glucose dehydrogenase (PQQGDH)with improved enzymatic characteristics was constructed by substitutingand combining the gene-encoding protein regions responsiblefor EDTA tolerance, thermal stability and substrate specificity.The protein region responsible for complete EDTA tolerance inAcinetobacter calcoaceticus, which is recognized as the indicatorof high stability in co-factor binding, was elucidated. Theregion is located between 32 and 59% from the N-terminus ofA.calcoaceticus PQQGDH(A27 region) and also corressponds tothe same position from 32 to 59% from the N-terminus in E.coliPQQGDH, though E.coli PQQGDH is EDTA sensitive. We previouslyreported that the C-terminal 3% region of A.calcoaceticus (A3region) played an important role in the increase of thermalstability, and that His775Asn substitution in E.coli PQQGDHresulted in an increase in the substrate specificity of E.coliPQQGDH towards glucose. Based on these findings, chimeric and/ormutated PQQGDHs, E97A3 H775N, E32A27E41 H782N, E32A27E38A3 andE32A27E38A3 H782N were constructed to investigate the compatibilityof two protein regions and one amino acid substitution. His775substitution to Asn corresponded to His782 substitution to Asn(H782N) in chimeric enzymes harbouring the A27 region. Sinceall the chimeric PQQGDHs harbouring the A27 region were EDTAtolerant, the A27 region was found to be compatible with theother region and substituted amino acid responsible for theimprovement of enzymatic properties. The contribution of theA3 region to thermal stability complemented the decrease inthe thermal stability due to the His775 or His782 substitutionto Asn. E32A27E38A3 H782N, which harbours all the above mentionedthree regions, showed improved EDTA tolerance, thermal stabilityand substrate specificity. These results suggested a strategyfor the construction of a semi-artificial enzyme by substitutingand combining the gene-encoding protein regions responsiblefor the improvement of enzyme characteristics. The characteristicsof constructed chimeric PQQGDH are discussed based on the predictedmodel, ß-propeller structure.     

Chemical synthesis of the RGD-protein decorsin: Pro-->Ala replacement reduces protein thermostability

Decorsin is a 39-residue polypeptide chain, crosslinked by three disulfide bridges, that strongly inhibits platelet aggregation. We report the chemical synthesis and characterization of analogs of decorsin with the aim of investigating the role of proline residues in protein structure, stability and biological activity. Decorsin analogs have been synthesized in which one (P23A and P24A decorsin) or two (P23,24A decorsin) proline residues have been substituted by alanine. The crude synthetic polypeptides were purified by reversed-phase HPLC in their reduced form and allowed to refold oxidatively to their disulfide-crosslinked species. The homogeneity of the synthetic mini-proteins, and also the correct pairing of the three disulfide bridges, were established by a number of analytical criteria, including fingerprinting analysis of the refolded synthetic analogs by using thermolysin and proteinase K as proteolytic enzymes. Replacement of proline by alanine results in a significant and cumulative decrease of the high thermal stability (Tm 74 degrees C) of native decorsin. The mono-substituted analogs display a Tm of 66-67 degrees C, while the double-substituted analog a Tm of 50 degrees C. On the other hand, the overall secondary and tertiary structures were not affected by the Pro-->Ala exchanges, as judged from circular dichroism measurements. Platelet aggregation assays established that the proline substitutions do not impair significantly the biological activity of decorsin. The results of this study clearly indicate that proline residues contribute significantly to the protein thermal stability. Our results are in line with the 'proline rule', previously advanced for explaining the unusual thermal stability of thermophilic enzymes, which usually show an enhanced content of proline residues with respect to their mesophilic counterparts.     

Random exchanges of non-conserved amino acid residues among four parental termite cellulases by family shuffling improved thermostability

There have been two major problems preventing applications of termite cellulases; one was difficulty for their hetelologous overexpression, and another is their low thermostability. We previously achieved adaptation of termite cellulase genes to an overexpression system of Escherichia coli by family shuffling of four orthologous cDNAs (Biosci. Biotechnol. Biochem., 2005; 69: 1711-1720). Using the adapted mutant cDNAs as parental genes combined with native-form cDNAs, we performed further family shuffling and obtained mutant cDNAs, which gave enzymes with improved thermostability. The best-evolved clone (PA68) was improved by 10 degrees C in maximum stability (retaining 90% original activity for 30 min incubation) from the parental enzymes, and kept 54% of its original activity for 150 min at 50 degrees C, whereas the most thermostable enzyme amongst the parents (A18) retained 30% of its original activity. PA68 showed 889 (micromoles of reducing sugars/min/mg of protein) in V(max) and 560 (micromoles of reducing sugars/min/mg of protein) in the specific activity against carboxymethylcellulose, which corresponds to 9.8 and 13.1 times of those of one of the ancestral enzymes rRsEG. In summary, we improved thermostability of the termite cellulase and increased the V(max) value and specific activity by combining only cDNAs encoding enzymes adapted for normal temperatures.     

Selection of a thermostable variant of chloramphenicol acetyltransferase (Cat-86)

The moderate thermophile Bacillus stearothermophilus was usedas a host in which to detect more thermostable variants of theB.pumilus chloramphenicol acetyltransferase (Cat-86) protein.Seventeen mutants were isolated and detected by their abilityto grow in the presence of chloramphenicol at a previously restrictivetemperature (58°C). The genes encoding these proteins weresequenced; all 17 mutants carried the same C to T transitionthat conferred an amino acid substitution of alanine by valineat position 203 of the protein sequence. The wild-type and onemutant Cat-86 protein were purified to homogeneity using affinitychromatography, and kinetic and thermal stability studies wereundertaken. Both enzymes had similar sp. act. in the regionof 215 U/mg, with K_m values for chloramphenicol in the range13.8–15.4 µM and for acetyl CoA in the range 13.6–15.5µM. The A203V mutant shows greater stability than thewild-type Cat-86 protein at temperatures above 50°C andappears to pass through a transition state between 48 and 50°C.     

A simple dual selection for functionally active mutants of Plasmodium falciparum dihydrofolate reductase with improved solubility

Sufficient solubility of the active protein in aqueous solution is a prerequisite for crystallization and other structural studies of proteins. In this study, we have developed a simple and effective in vivo screening system to select for functionally active proteins with increased solubility by using Plasmodium falciparum dihydrofolate reductase (pfDHFR), a well-known malarial drug target, as a model. Prior to the dual selection process, pfDHFR was fused to green fluorescent protein (GFP), which served as a reporter for solubility. The fusion gene was used as a template for construction of mutated DNA libraries of pfDHFR. Two amino acids with large hydrophobic side chains (Y35 and F37) located on the surface of pfDHFR were selected for site-specific mutagenesis. Additionally, the entire pfDHFR gene was randomly mutated using error-prone PCR. During the first step of the dual selection, mutants with functionally active pfDHFR were selected from two libraries by using bacterial complementation assay. Fluorescence signals of active mutants were subsequently measured and five mutants with increased GFP signal, namely Y35Q + F37R, Y35L + F37T, Y35G + F37L and Y35L + F37R from the site-specific mutant library and K27E from the random mutant library, were recovered. The mutants were expressed, purified and characterized as monofunctional pfDHFR following excision of GFP. Our studies indicated that all mutant pfDHFRs exhibited kinetic properties similar to that of the wild-type protein. For comparison of protein solubility, the maximum concentrations of mutant enzymes prior to aggregation were determined. All mutants selected in this study exhibited 3- to 6-fold increases in protein solubility compared with the wild-type protein, which readily aggregated at 2 mg/ml. The dual selection system we have developed should be useful for engineering functionally active protein mutants with sufficient solubility for functional/structural studies and other applications.     

A33scFv-green fluorescent protein, a recombinant single-chain fusion protein for tumor targeting

Chemical conjugates of monoclonal antibodies with fluorophores or enzymes have long been used for diagnostic purposes and experimental therapeutic approaches. Recombinant technology allows for the design and expression of tailored genuine fusion proteins, providing defined molecules as to size, molar ratios of the functional components and stability. The production of functional protein, however, is often limited or impossible due to refolding and solubility problems. Here, we report on the production of a soluble recombinant fusion construct, A33scFv-green fluorescent protein (A33scFv::GFP) in Pichia pastoris. A33scFv is a single-chain antibody recognizing the A33 antigen, which is expressed by approximately 95% of colorectal carcinomas and has become a focus of pre-clinical and clinical investigation. The fusion partner GFP was selected both as an experimental tool for functional studies of the A33 antigen and as a potential diagnostic for colon cancer detection and therapy planning. Pichia pastoris yeast strains were transformed with A33scFv::GFP cDNA under the methanol-inducible AOX1 promotor. The construct was properly expressed and secreted into culture supernatants as a soluble protein, which was bifunctional without additional renaturation or solubilization steps. The crude protein solution was purified by affinity chromatography. Surface plasmon resonance, flow cytometry and fluorescence microscopy on sections of normal and cancerous colon tissue revealed specific binding and the applicability of this fusion protein for diagnostic purposes. In addition, the biodistribution of A33scFv::GFP was analyzed in mice bearing A33-positive tumor xenografts, confirming specific tumor targeting.     

Thermostable variants of bovine {beta}-lactoglobulin

The thermal stability of bovine ß-lactoglobulin (BLG)has been enhanced by the introduction of an additional disulfidebond. Wild-type BLG has two disulfide bonds, C106–C119and C66–C160, with a free cysteine at position 121. Wehave designed, with the aid of molecular modeling calculations,two mutants of a recombinant BLG (rBLG), L104C and A132C. Moleculardynamics simulations were performed at 300K to study the effectof these alterations on the conformation of the protein. Thesemutants were then created by site-directed mutagenesis and purifiedfrom Escherichia coli carrying a tac expression vector usinga two-step renaturation method. Formation of disulfide linkagesin the correct arrangement, as designed, was confirmed by peptidemapping. In contrast to wild-type rBLG, which polymerizes attemperatures >65°C, neither of the mutant proteins polymerized.The conformational stability of the L104C and A132C mutant proteinsagainst thermal denaturation has been substantially increased(8- 10°C) as compared with wild-type rBLG. Furthermore,the A132C rBLG exhibits an enhanced stability against denaturationby guanidine hydrocnloride as compared with the wild-type orL104C rBLG     

Emulsifying performance of modular beta-sandwich proteins: the hydrophobic moment and conformational stability

Our understanding of protein emulsifying properties is largely based on analysis of emulsifiers found in milk and seed. The 9th-10th type III fibronectin domain pair retains full biological activity following emulsification-encapsulation into polyester microspheres, for controlled delivery, but the conformational criteria determining emulsification efficiency (EE) are unknown. Here, we have generated a series of mutants of this beta-sandwich protein, changing the hydrophobic moment and conformational stability, to investigate the structure-emulsification relationship. Predictive modelling of the hydrophobic moment of beta-strands and mutations known to increase conformational stability were used to generate the series. The proteins were tested for their emulsion stability and EE for oil-in-water mixtures. We show that the stabilization of emulsions by beta-sandwich proteins is best predicted by conformational stability during equilibrium denaturation in ionic surfactant. In contrast, the EE of these proteins is inversely related to an increase in their surface hydrophobicity following unfolding in surfactant. We also describe a novel beta-sandwich emulsifier with strong EE. The requirement for interdomain flexibility to achieve maximum emulsion stability and EE is also shown. This work increases our understanding of the mechanisms involved in protein emulsification and will be of use to the microencapsulation of proteins into polyester microspheres via emulsion-extraction protocols.     

Design of fully active FGF-1 variants with increased stability

Fibroblast growth factor 1 is a powerful mitogen playing an important role in morphogenesis, angiogenesis and wound healing and is therefore of potential medical interest. Using homologous sequence and structure comparisons, we designed and constructed 16 mutants of FGF-1 with increased thermodynamic stability, as determined by chemical and heat denaturation. For multiple mutants, additive effects on stability were observed, providing mutants up to 7.8 degrees C more stable than the wild-type. None of the introduced mutations affected any FGF-1 biological activities, such as stimulation of DNA synthesis, MAP kinase activation and binding to the FGF receptor on the cell surface. Our study provides a good starting point to improve the stability of FGF-1 in the context of its wide potential therapeutic applications. We showed that a homology approach is an effective method to change the thermodynamic properties of the protein without altering its function.     

Characterization of active-site aromatic residues in xylanase A from Streptomyces lividans

The role of four aromatic residues (W85, Y172, W266 and W274)in the structure–function relationship in xylanase A fromStreptomyces lividans (XlnA) was investigated by site-directedmutagenesis where each residue was subjected to three substitutions(W85A/H/F; W266A/H/F; W274A/H/F and Y172A/F/S). These four aminoacids are highly conserved among family 10 xylanases and structuraldata have implicated them in substrate binding at the activesite. Far-UV circular dichroism spectroscopy was used to showthat the overall structure of XlnA was not affected by any ofthese mutations. High-performance liquid chromatographic analysisof the hydrolysis products of birchwood xylan and xylopentaoseshowed that mutation of these aromatic residues did not alterthe enzyme's mode of action. As expected, though, it did reducethe affinity of XlnA for birchwood xylan. A comparison of thekinetic parameters of different mutants at the same positiondemonstrated the importance of the aromatic nature of W85, Y172and W274 in substrate binding. Replacement of these residuesby a phenylalanine resulted in mutant proteins with a K_M closerto that of the wild-type protein in comparison with the othermutations analyzed. The kinetic analysis of the mutant proteinsat position W266 indicated that this amino acid is importantfor both substrate binding and efficient catalysis by XlnA.These studies also demonstrated the crucial role of these activesite aromatic residues for the thermal stability of XlnA.     

Effect of Lys->Arg mutation on the thermal stability of Cu,Zn superoxide dismutase: influence on the monomer-dimer equilibrium

The thermal stability of two single (K3R, K67R) and one double(K3R-K67R) mutants of Xenopus laevis B Cu,Zn superoxide dismutasehas been studied to test LysArg substitution as an ‘electrostaticallyconservative’ strategy to increase protein stability.The K3R mutant displays an increased thermostability with respectto the wild-type enzyme, whilst a decreased stability was observedin the case of the K67R and K3R-K67R mutants. Concentrationdependence of the apparent inactivation constant (k_app) of thelatter mutants, as compared to that of the wild type enzymeand K3R mutant, indicates that their higher sensitivity to heatinactivation is due to a perturbation of the dimer association.These results are confirmed also by fluorescence anisotropymeasurements of the internal probe Tyr149. The possible roleof Arg67 in perturbing the dimer dissociation equilibrium towardthe monomeric form is discussed.     

Conservation of mechanism, variation of rate: folding kinetics of three homologous four-helix bundle proteins

The amino acid sequence of a protein determines both its final folded structure and the folding mechanism by which this structure is attained. The differences in folding behaviour between homologous proteins provide direct insights into the factors that influence both thermodynamic and kinetic properties. Here, we present a comprehensive thermodynamic and kinetic analysis of three homologous homodimeric four-helix bundle proteins. Previous studies with one member of this family, Rop, revealed that both its folding and unfolding behaviour were interesting and unusual: Rop folds (k(0)(f) = 29 s(-1)) and unfolds (k(0)(u) = 6 x 10(-7) s(-1)) extremely slowly for a protein of its size that contains neither prolines nor disulphides in its folded structure. The homologues we discuss have significantly different stabilities and rates of folding and unfolding. However, the rate of protein folding directly correlates with stability for these homologous proteins: proteins with higher stability fold faster. Moreover, in spite of possessing differing thermodynamic and kinetic properties, the proteins all share a similar folding and unfolding mechanism. We discuss the properties of these naturally occurring Rop homologues in relation to previously characterized designed variants of Rop.     

Ribosome display of mammalian receptor domains

Many mammalian receptor domains, among them a large number of potential therapeutic target proteins, are highly aggregation-prone upon heterologous expression in bacteria. This severely limits functional studies of such receptor domains and also their engineering towards improved properties. One of these proteins is the Nogoreceptor, which plays a central role in mediating the inhibition of axon growth and functional recovery after injury of the adult mammalian central nervous system. We show here that the ligand binding domain of the Nogoreceptor folds to an active conformation in ternary ribosomal complexes, as formed in ribosome display. In these complexes the receptor is still connected, via a C-terminal tether, to the peptidyl tRNA in the ribosome and the mRNA also stays connected. The ribosome prevents aggregation of the protein, which aggregates as soon as the release from the ribosome is triggered. In contrast, no active receptor was observed in phage display, where aggregation appears to prevent incorporation of the protein into the phage coat. This strategy sets the stage for rapidly studying defined mutations of such aggregation-prone receptors in vitro and to improve their properties by in vitro evolution using the ribosome display technology.     

Improved peptide function from random mutagenesis over short 'windows'

We have applied random mutagenesis over short contiguous residuetracts (‘windows’) within an active peptide (the-peptide of ß-galactosidase) such that all windowresidues are replaced simultaneously. A novel technique usingmixed synthetic oligonucleotides and selection against an EcoKrestrictionsite has allowed the construction of libraries of mutants fortwo separate windows, sites A and B. Mutant phenotypes can beeasily assessed in vivoby a complementation test, and panelsof mutants have been quantitatively tested in vivoThis allowedthe rapid probing of structural requirements for each site.The two windows yielded markedly disparate results. Site B wasmuch less stringent in its sequence requirements for significantfunction than Site A, and mutants with improved function wereisolated at Site B alone. In addition, one Site B mutant withwild-type levels of activity showed enhanced stability to heator a protein denaturant. We propose that short tracts with thecharacteristics of Site B constitute ‘secondary’interaction sites which are more tolerant of sequence diversity.Random manipulation of such secondary sites is thus more likelyto yield up-mutations for standard or altered environments.Window mutagenesis can in principle be applied to any protein-proteinor protein-Ugand interaction.