In vitro evolution of a hyperstable Gbeta1 variant

An in-vitro selection strategy was used to obtain strongly stabilized variants of the beta1 domain of protein G (Gbeta1). In a two-step approach, first candidate positions with a high potential for stabilization were identified in Gbeta1 libraries that were created by error-prone PCR, and then, after randomization of these positions by saturation mutagenesis, strongly stabilized variants were selected. For both steps the in-vitro selection method Proside was employed. Proside links the stability of a protein with the infectivity of a filamentous phage. Ultimately, residues from the two best selected variants were combined in a single Gbeta1 molecule. This variant with the four mutations E15V, T16L, T18I, and N37L showed an increase of 35.1 degrees C in the transition midpoint and of 28.5 kJ mol(-1) (at 70 degrees C) in the Gibbs free energy of stabilization. It was considerably more stable than the best variant from a previous Proside selection, in which positions were randomized that had originally been identified by computational design. Only a single substitution (T18I) was found in both selections. The best variants from the present selection showed a higher cooperativity of thermal unfolding, as indicated by an increase in the enthalpy of unfolding by about 60 kJ mol(-1). This increase is apparently correlated with the presence of Leu residues that were selected at the positions 16 and 37.     

Optimization of the gbeta1 domain by computational design and by in vitro evolution: structural and energetic basis of stabilization

Computational design and in vitro evolution are major strategies for stabilizing proteins. For the four critical positions 16, 18, 25, and 29 of the B domain of the streptococcal protein G (Gbeta1), they identified the same optimal residues at positions 16 and 25, but not at 18 and 29. Here we analyzed the energetic contributions of the residues from these two approaches by single and double mutant analyses and determined crystal structures for a variant from the calculation (I16/L18/E25/K29) and from the selection (I16/I18/E25/F29). The structural analysis explains the observed differences in stabilization. Residues 16, 18, and 29 line an invagination, which results from a packing defect between the helix and the beta-sheet of Gbeta1. In all stabilized variants, residues with larger side-chains occur at these positions and packing is improved. In the selected variant, packing is better optimized than in the computed variant. Such differences in side-chain packing strongly affect stability but are difficult to evaluate by computation.     

Optimized variants of the cold shock protein from in vitro selection: structural basis of their high thermostability

The bacterial cold shock proteins (Csp) are widely used as models for the experimental and computational analysis of protein stability. In a previous study, in vitro evolution was employed to identify strongly stabilizing mutations in Bs-CspB from Bacillus subtilis. The best variant found by this approach contained the mutations M1R, E3K and K65I, which raised the midpoint of thermal unfolding of Bs-CspB from 53.8 degrees C to 83.7 degrees C, and increased the Gibbs free energy of stabilization by 20.9 kJ mol(-1). Another selected variant with the two mutations A46K and S48R was stabilized by 11.1 kJ mol(-1). To elucidate the molecular basis of these stabilizations, we determined the crystal structures of these two Bs-CspB variants. The mutated residues are generally well ordered and provide additional stabilizing interactions, such as charge interactions, additional hydrogen bonds and improved side-chain packing. Several mutations improve the electrostatic interactions, either by the removal of unfavorable charges (E3K) or by compensating their destabilizing interactions (A46K, S48R). The stabilizing mutations are clustered at a contiguous surface area of Bs-CspB, which apparently is critically important for the stability of the beta-barrel structure but not well optimized in the wild-type protein.     

Retarded protein folding of deficient human alpha 1-antitrypsin D256V and L41P variants

alpha(1)-Antitrypsin is the most abundant protease inhibitor in plasma and is the archetype of the serine protease inhibitor superfamily. Genetic variants of human alpha(1)-antitrypsin are associated with early-onset emphysema and liver cirrhosis. However, the detailed molecular mechanism for the pathogenicity of most variant alpha(1)-antitrypsin molecules is not known. Here we examined the structural basis of a dozen deficient alpha(1)-antitrypsin variants. Unlike most alpha(1)-antitrypsin variants, which were unstable, D256V and L41P variants exhibited extremely retarded protein folding as compared with the wild-type molecule. Once folded, however, the stability and inhibitory activity of these variant proteins were comparable to those of the wild-type molecule. Retarded protein folding may promote protein aggregation by allowing the accumulation of aggregation-prone folding intermediates. Repeated observations of retarded protein folding indicate that it is an important mechanism causing alpha(1)-antitrypsin deficiency by variant molecules, which have to fold into the metastable native form to be functional.     

Increased folding stability of TEM-1 beta-lactamase by in vitro selection

In vitro selections of stabilized proteins lead to more robust enzymes and, at the same time, yield novel insights into the principles of protein stability. We employed Proside, a method of in vitro selection, to find stabilized variants of TEM-1 β-lactamase from Escherichia coli. Proside links the increased protease resistance of stabilized proteins to the infectivity of a filamentous phage. Several libraries of TEM-1 β-lactamase variants were generated by error-prone PCR, and variants with increased protease resistance were obtained by raising temperature or guanidinium chloride concentration during proteolytic selections. Despite the small size of phage libraries, several strongly stabilizing mutations could be obtained, and a manual combination of the best shifted the profiles for thermal unfolding and temperature-dependent inactivation of β-lactamase by almost 20 °C to a higher temperature. The wild-type protein unfolds in two stages: from the native state via an intermediate of the molten-globule type to the unfolded form. In the course of the selections, the native protein was stabilized by 27 kJ mol^− 1 relative to the intermediate and the cooperativity of unfolding was strongly increased. Three of our stabilizing replacements (M182T, A224V, and R275L) had been identified independently in naturally occurring β-lactamase variants with extended substrate spectrum. In these variants, they acted as global suppressors of destabilizations caused by the mutations in the active site. The comparison between the crystal structure of our best variant and the crystal structure of the wild-type protein indicates that most of the selected mutations optimize helices and their packing. The stabilization by the E147G substitution is remarkable. It removes steric strain that originates from an overly tight packing of two helices in the wild-type protein. Such unfavorable van der Waals repulsions are not easily identified in crystal structures or by computational approaches, but they strongly reduce the conformational stability of a protein.     

Folding and stability of the leucine-rich repeat domain of internalin B from Listeri monocytogenes

Internalin B (InlB), a surface protein of the human pathogen Listeria monocytogenes, promotes invasion into various host cell types by inducing phagocytosis of the entire bacterium. The N-terminal half of InlB (residues 36-321, InlB321), which is sufficient for this process, contains a central leucine-rich repeat (LRR) domain that is flanked by a small alpha-helical cap and an immunoglobulin (Ig)-like domain. Here we investigated the spectroscopic properties, stability and folding of InlB321 and of a shorter variant lacking the Ig-like domain (InlB248). The circular dichroism spectra of both protein variants in the far ultraviolet region are very similar, with a characteristic minimum found at approximately 200 nm, possibly resulting from the high 3(10)-helical content in the LRR domain. Upon addition of chemical denaturants, both variants unfold in single transitions with unusually high cooperativity that are fully reversible and best described by two-state equilibria. The free energies of GdmCl-induced unfolding determined from transitions at 20 degrees C are 9.9(+/-0.8)kcal/mol for InlB321 and 5.4(+/-0.4)kcal/mol for InlB248. InlB321 is also more stable against thermal denaturation, as observed by scanning calorimetry. This suggests, that the Ig-like domain, which presumably does not directly interact with the host cell receptor during bacterial invasion, plays a critical role for the in vivo stability of InlB.     

The N-terminal helix is a post-assembly clamp in the bacterial outer membrane protein PagP

The Escherichia coli outer membrane beta-barrel enzyme PagP and its homologues are unique in that the eight-stranded barrel is tilted by about 25 degrees with respect to the membrane normal and is preceded by a 19-residue amphipathic alpha-helix. To investigate the role of this helix in the folding and stability of PagP, mutants were generated in which the helix was deleted (Delta(1-19)), or in which residues predicted to be involved in helix-barrel interactions were altered (W17A or R59L). The ability of the variants to insert into detergent micelles or liposomes was studied in vitro using circular dichroism, fluorescence, Fourier transform infrared spectroscopy, electrophoretic mobility and gain of enzyme activity. The data show that PagP, initially unfolded in 5% (w/v) perfluoro-octanoic acid or 6 M guanidinium chloride, inserts spontaneously and folds quantitatively to an active conformation into detergent micelles of cyclofos-7 or into large vesicles of diC(12:0)-phosphatidylcholine (diC(12:0)PC), respectively, the latter in the presence of 7 M urea. Successful refolding of all variants into both micelles and liposomes ruled out an essential role for the helix or helix-barrel interactions in folding and membrane insertion. Measurements of thermal stability indicated that the variants R59L, W17A/R59L and Delta(1-19) were destabilised substantially compared with wild-type PagP. However, in contrast to the other variants, destabilisation of the W17A variant relative to wild-type PagP was much greater in liposomes than in micelles. Analysis of the kinetics of folding and unfolding of all variants in diC(12:0)PC liposomes suggested that this destabilisation arises predominantly from an increased dissociation of the refolded variant proteins from the lipid-inserted state. The data support the view that the helix of PagP is not required for folding and assembly, but instead acts as a clamp, stabilising membrane-inserted PagP after folding and docking with the membrane are complete.     

Evolution of Stability in a Cold-Active Enzyme Elicits Specificity Relaxation and Highlights Substrate-Related Effects on Temperature Adaptation

Molecular aspects of thermal adaptation of proteins were studied by following the co-evolution of temperature dependence, conformational stability, and substrate specificity in a cold-active lipase modified via directed evolution. We found that the evolution of kinetic stability was accompanied by a relaxation in substrate specificity. Moreover, temperature dependence and selectivity turned out to be mutually dependent. While the wild-type protein was strictly specific for short-chain triglycerides (C4) in the temperature range 10-50 °C and displayed highest activity in the cold, its stabilized variant was able to accept C8 and C12 molecules and its selectivity was temperature dependent. We could not detect any improvement in the overall structural robustness of the mutant when the structure was challenged by temperature or chemical denaturants. There is, however, strong evidence for local stabilization effects in the active-site region provided by two independent approaches. Differential scanning fluorimetry revealed that the exposure of hydrophobic patches (as the active site is) precedes denaturation, and molecular dynamics simulations confirmed that stability was obtained by restriction of the mobility of the lid, a flexible structure that regulates the access to the enzyme active site and influences its stability. This reduction of lid movements is suggested to be accompanied by a concomitant increase in the mobility of other protein regions, thus accounting for the observed broadening of substrate specificity.     

Enhancement of transglycosylation activity by construction of chimeras between mesophilic and thermophilic beta-glucosidase

The family 3 beta-glucosidase from Thermotoga maritima is a highly thermostable enzyme (85 degrees C) that displays transglycosylation activity. In contrast, the beta-glucosidase from Cellvibrio gilvus is mesophilic (35 degrees C) and displays no such transglycosylation activity. Both enzymes consist of two domains, an N-terminal and a C-terminal domain, and the amino acid identities between the two enzymes in these domains are 32.4 and 36.4%, respectively. In an attempt to identify the molecular basis underpinning the display of transglycosylation activity and the requirements for thermal stability, eight chimeric genes were constructed by shuffling the two parental beta-glucosidase genes at four selected borders, two in the N-terminal domain and two in the C-terminal domain. Of the eight chimeric genes constructed, only two chimeric enzymes (Tm578/606Cg and Tm638/666Cg) gave catalytically active forms and these were the ones shuffled in the C-terminal domain. For these active chimeric enzymes, 80% (Tm578/606Cg) and 88% (Tm638/666Cg) of their amino acid sequences originated from T. maritima. With regard to their thermal profiles, the two active chimeric enzymes, Tm578/606Cg and Tm638/666Cg, displayed profiles intermediate to those of the two parental enzymes as they were optimally active at 65 and 70 degrees C, respectively. These two chimeric enzymes were optimally active at pH 4.1 and 3.9, which is closer to that observed for the T. maritima enzyme (pH 3.2-3.5) than that for the C. gilvus enzyme (pH 6.2-6.5). Kinetic parameters for the chimeric enzymes were investigated with five different substrates including pNP-beta-D-glucopyranoside. The kinetic parameters obtained for the chimeric enzymes were closer to those of the T. maritima enzyme than to those of the C. gilvus enzyme. Transglycosylation activity was observed for both chimeric enzymes and the activity of the Tm578/606Cg chimera was at a level twice that observed with the T. maritima enzyme. This study is an effective demonstration of the usefulness of chimeric enzymes in altering the characteristics of an enzyme.     

Evolutionary protein stabilization in comparison with computational design

Two major strategies are currently used for stabilizing proteins: in vitro evolution and computational design. Here, we used gene libraries of the beta1 domain of the streptococcal protein G (Gbeta1) and Proside, an in vitro selection method, to identify stabilized variants of this protein. In the Gbeta1 libraries, the codons for the four boundary positions 16, 18, 25, and 29 were randomized. Many Gbeta1 variants with strongly increased thermal stabilities were found in 11 selections performed with five independent libraries. Previously, Mayo and co-workers used computational design to stabilize Gbeta1 by sequence optimization at the same positions. Their best variant ranked third within the panel of the selected variants. None of the ten computed sequences was found in the Proside selections, because several computed residues for positions 18 and 29 were not optimal for stability.     

High tolerance to mutations in a Chlamydia trachomatis peptide deformylase loop

AIM: To determine if and how a loop region in the peptide deformylase (PDF) of Chlamydia trachomatis regulates enzyme function.METHODS: Molecular dynamics simulation was used to study a structural model of the chlamydial PDF (cPDF) and predict the temperature factor per residue for the protein backbone atoms. Site-directed mutagenesis was performed to construct cPDF variants. Catalytic properties of the resulting variants were determined by an enzyme assay using formyl-Met-Ala-Ser as a substrate.RESULTS: In silico analysis predicted a significant increase in atomic motion in the DGELV sequence (residues 68-72) of a loop region in a cPDF mutant, which is resistant to PDF inhibitors due to two amino acid substitutions near the active site, as compared to wild-type cPDF. The D68R and D68R/E70R cPDF variants demonstrated significantly increased catalytic efficiency. The E70R mutant showed only slightly decreased efficiency. Although deletion of residues 68-72 resulted in a nearly threefold loss in substrate binding, this deficiency was compensated for by increased catalytic efficiency.CONCLUSION: Movement of the DGELV loop region is involved in a rate-limiting conformational change of the enzyme during catalysis. However, there is no stringent sequence requirement for this region for cPDF enzyme activity.     

Relationship between stability and flexibility in the most flexible region of Photinus pyralis luciferase

Firefly luciferase is a protein with a large N-terminal and a small C-terminal domain. B-factor analysis shows that its C-terminal is much more flexible than its N-terminal. Studies on hyperthermophile proteins have been shown that the increased thermal stability of hyperthermophile proteins is due to their enhanced conformational rigidity and the relationship between flexibility, stability and function in most of proteins is on debate. Two mutations (D474K and D476N) in the most flexible region of firefly luciferase were designed. Thermostability analysis shows that D476N mutation doesn't have any significant effect but D474K mutation destabilized protein. On the other hand, flexibility analysis using dynamic quenching and limited proteolysis demonstrates that D474K mutation became much more flexible than wild type although D476N doesn't have any significant difference. Intrinsic and ANS fluorescence studies demonstrate that D476N mutation is brought about by structural changes without significant effect on thermostability and flexibility. Molecular modeling reveals that disruption of a salt bridge between D⁴⁷⁴ and K⁴⁴⁵ accompanying with some H-bond deletion may be involved in destabilization of D474K mutant.     

Stabilization of yeast hexokinase A by polyol osmolytes: correlation with the physicochemical properties of aqueous solutions

Osmolytes of the polyol series are known to accumulate in biological systems under stress and stabilize the structures of a wide variety of proteins. While increased surface tension of aqueous solutions has been considered an important factor in protein stabilization effect, glycerol is an exception, lowering the surface tension of water. To clarify this anomalous effect, the effect of a series of polyols on the thermal stability of a highly thermolabile two domain protein yeast hexokinase A has been investigated by differential scanning calorimetry and by monitoring loss in the biological activity of the enzyme as a function of time. A larger increase in the T(m) of domain 1 compared with that of domain 2, varying linearly with the number of hydroxyl groups in polyols, has been observed, sorbitol being the best stabilizer against both thermal as well as urea denaturation. Polyols help retain the activity of the enzyme considerably and a good correlation of the increase in T(m) (DeltaT(m)) and the retention of activity with the increase in the surface tension of polyol solutions, except glycerol, which breaks this trend, has been observed. However, the DeltaT(m) values show a linear correlation with apparent molal heat capacity and volume of aqueous polyol solutions including glycerol. These results suggest that while bulk solution properties contribute significantly to protein stabilization, interfacial properties are not always a good indicator of the stabilizing effect. A subtle balance of various weak binding and exclusion effects of the osmolytes mediated by water further regulates the stabilizing effect. Understanding these aspects is critical in the rational design of stable protein formulations.     

Asp-196-->Ala mutant of Leuconostoc mesenteroides sucrose phosphorylase exhibits altered stereochemical course and kinetic mechanism of glucosyl transfer to and from phosphate

Mutagenesis of Asp-196 into Ala yielded an inactive variant of Leuconostoc mesenteroides sucrose phosphorylase (D196A). External azide partly complemented the catalytic defect in D196A with a second-order rate constant of 0.031 M-1 s-1 (pH 5, 30 degrees C) while formate, acetate and halides could not restore activity. The mutant utilized azide to convert alpha-D-glucose 1-phosphate into beta-D-glucose 1-azide, reflecting a change in stereochemical course of glucosyl transfer from alpha-retaining in wild-type to inverting in D196A. Phosphorolysis of beta-D-glucose 1-azide by D196A occurred through a ternary complex kinetic mechanism, in marked contrast to the wild-type whose reactions feature a common glucosyl enzyme intermediate and Ping-Pong kinetics. Therefore, Asp-196 is identified unambiguously as the catalytic nucleophile of sucrose phosphorylase, and its substitution by Ala forces the reaction to proceed via single nucleophilic displacement. D196A is not detectably active as alpha-glucosynthase.     

A polymorphic position in electron transfer flavoprotein modulates kinetic stability as evidenced by thermal stress

The electron transfer flavoprotein (ETF) is a hub interacting with at least 11 mitochondrial flavoenzymes and linking them to the respiratory chain. Here we report the effect of the ETFα-T/I171 polymorphism on protein conformation and kinetic stability under thermal stress. Although variants have comparable thermodynamic stabilities, kinetically their behavior is rather distinct as ETFα-T171 displays increased susceptibility to cofactor flavin adenine dinucleotide (FAD) loss and enhanced kinetics of inactivation during thermal stress. Mimicking a fever episode yields substantial activity loss. However, the presence of substoichiometric concentrations of GroEL is sufficient to act as an effective buffer against long-term thermal denaturation. Our investigations are compatible with the notion that the ETFα-T171 variant displays an altered conformational landscape that results in reduced protein function under thermal stress.     

The crystal structure of the cis-proline to glycine variant (P114G) of ribonuclease A

Replacement of a cis-proline by glycine at position 114 in ribonuclease A leads to a large decrease in thermal stability and simplifies the refolding kinetics. A crystallographic approach was used to determine whether the decrease in thermal stability results from the presence of a cis glycine peptide bond, or from a localized structural rearrangement caused by the isomerization of the mutated cis 114 peptide bond. The structure was solved at 2.0 A resolution and refined to an R-factor of 19.5% and an R(free) of 21.9%. The overall conformation of the protein was similar to that of wild-type ribonuclease A; however, there was a large localized rearrangement of the mutated loop (residues 110-117-a 9.3 A shift of the Calpha atom of residue 114). The peptide bond before Gly114 is in the trans configuration. Interestingly, a large anomalous difference density was found near residue 114, and was attributed to a bound cesium ion present in the crystallization experiment. The trans isomeric configuration of the peptide bond in the folded state of this mutant is consistent with the refolding kinetics previously reported, and the associated protein conformational change provides an explanation for the decreased thermal stability.