Selection of stably folded proteins by phage-display with proteolysis

Protein stabilization with retained function of monellin using a split GFP system

Sweet proteins are an unexploited resource in the search for non-carbohydrate sweeteners mainly due to their low stability towards heating. Variants of the sweet protein monellin, with increased stability, were derived by an in vivo screening method based on the thermodynamic linkage between fragment complementation and protein stability. This approach depends on the correlation between mutational effects on the affinity between protein fragments and the stability of the intact protein. By linking the two fragments of monellin to the split GFP (green fluorescent protein) system, reconstitution of GFP was promoted and moderately fluorescent colonies were obtained. Two separate random libraries were produced for the monellin chains and the mutant clones were ranked based on fluorescence intensity. Mutants with increased affinity between the fragments, and subsequently increased stability, caused increased fluorescence intensity of split GFP. Single chain monellin variants of the top-ranked mutants for each chain, S76Y in the A-chain and W3C + R39G in the B-chain and all combinations thereof, were expressed and the increase in stability was verified by temperature denaturation studies using circular dichroism spectroscopy. Functionality studies showed that mutant S76Y has retained sweetness and has potential use within the food industry.

Computational design and selections for an engineered, thermostable terpene synthase

Terpenoids include structurally diverse antibiotics, flavorings, and fragrances. Engineering terpene synthases for control over the synthesis of such compounds represents a long sought goal. We report computational design, selections, and assays of a thermostable mutant of tobacco 5-epi-aristolochene synthase (TEAS) for the catalysis of carbocation cyclization reactions at elevated temperatures. Selection for thermostability included proteolytic digestion followed by capture of intact proteins. Unlike the wild-type enzyme, the mutant TEAS retains enzymatic activity at 65°C. The thermostable terpene synthase variant denatures above 80°C, approximately twice the temperature of the wild-type enzyme.

High-throughput analysis of the protein sequence-stability landscape using a quantitative yeast surface two-hybrid system and fragment reconstitution

Stability evaluation of many mutants can lead to a better understanding of the sequence determinants of a structural motif and of factors governing protein stability and protein evolution. The traditional biophysical analysis of protein stability is low throughput, limiting our ability to widely explore sequence space in a quantitative manner. In this study, we have developed a high-throughput library screening method for quantifying stability changes, which is based on protein fragment reconstitution and yeast surface display. Our method exploits the thermodynamic linkage between protein stability and fragment reconstitution and the ability of the yeast surface display technique to quantitatively evaluate protein-protein interactions. The method was applied to a fibronectin type III (FN3) domain. Characterization of fragment reconstitution was facilitated by the co-expression of two FN3 fragments, thus establishing a yeast surface two-hybrid method. Importantly, our method does not rely on competition between clones and thus eliminates a common limitation of high-throughput selection methods in which the most stable variants are recovered predominantly. Thus, it allows for the isolation of sequences that exhibit a desired level of stability. We identified more than 100 unique sequences for a beta-bulge motif, which was significantly more informative than natural sequences of the FN3 family in revealing the sequence determinants for the beta-bulge. Our method provides a powerful means for the rapid assessment of the stability of many variants, for the systematic assessment of the contribution of different factors to protein stability, and for enhancement of the protein stability.

Molecular diversity in engineered protein libraries

Engineered protein libraries, defined here as a collection of different mutant variants of a single specific protein, are intentionally designed to be rich in molecular diversity and can span ranges from as little as 400 different variants to greater than 10(12) members per library. The goal of engineering libraries is to generate new protein variants, identified upon screening, that possess desired novel properties. Exploitation of the natural organization of the genetic code has led to 'focused' libraries that are lower in overall complexity yet biased towards variants with preferred biophysical properties. An emerging trend, in which computational algorithms are blended with in vivo screens, is also leading towards greater and more rapid success in the field of protein design.

When the Surface Tells What Lies Beneath: Combinatorial Phage-display Mutagenesis Reveals Complex Networks of Surface–Core Interactions in the Pacifastin Protease Inhibitor Family

Pacifastin protease inhibitors are small cysteine-rich motifs of approximately 35 residues that were discovered in arthropods. The family is divided into two related groups on the basis of the composition of their minimalist inner core. In group I, the core is governed by a Lys10-Trp26 interaction, while in group II it is organized around Phe10. Group I inhibitors exhibit intriguing taxon specificity: potent arthropod-trypsin inhibitors from this group are almost inactive against vertebrate enzymes. The group I member SGPI-1 and the group II member SGPI-2 are extensively studied inhibitors. SGPI-1 is taxon-selective, while SGPI-2 is not. Individual mutations failed to explain the causes underlying this difference. We deciphered this phenomenon using comprehensive combinatorial mutagenesis and phage display. We produced a complete chimeric SGPI-1 / SGPI-2 inhibitor-phage library, in which the two sequences were shuffled at the highest possible resolution of individual residues. The library was selected for binding to bovine trypsin and crayfish trypsin. Sequence analysis of the selectants revealed that taxon specificity is due to an intra-molecular functional coupling between a surface loop and the Lys10-Trp26 core. Five SGPI-2 surface residues transplanted into SGPI-1 resulted in a variant that retained the "taxon-specific" core, but potently inhibited both vertebrate and arthropod enzymes. An additional rational point mutation resulted in a picomolar inhibitor of both trypsins. Our results challenge the generally accepted view that surface residues are the exclusive source of selectivity for canonical inhibitors. Moreover, we provide important insights into general principles underlying the structure-function properties of small disulfide-rich polypeptides, molecules that exist at the borderline between peptides and proteins.

Selection strategies for improved biocatalysts

Enzymes have become an attractive alternative to conventional catalysts in numerous industrial processes. However, their properties do not always meet the criteria of the application of interest. Directed evolution is a powerful tool for adopting the characteristics of an enzyme. However, selection of the evolved variants is a critical step, and therefore new strategies to enable selection of the desired enzymatic activity have been developed. This review focuses on these novel strategies for selecting enzymes from large libraries, in particular those that are used in the synthesis of pharmaceutical intermediates and pharmaceuticals.

Exploiting Elements of Transcriptional Machinery to Enhance Protein Stability

The correlation between protein structure and function is well established, yet the role stability/flexibility plays in protein function is being explored. Here, we describe an in vivo screen in which the thermal stability of a test protein is correlated directly to the transcriptional regulation of a reporter gene. The screen readout is independent of the function of the test protein, proteolytic resistance, solubility or propensity to aggregate indiscriminately, and is thus dependent solely on the overall stability of the test protein. The system entails the use of an engineered chimeric construct that consists of three covalently linked domains; a constant N-terminal DNA-binding domain, a variable central test protein, and a constant C-terminal transcriptional activation domain. The test proteins are mutant variants of the beta1 domain of streptococcal protein G that span fairly evenly a thermal stability range from as low as 38 degrees C to above 100 degrees C. When the chimeric construct contains a test variant of low thermal stability, the reporter gene is up-regulated to a greater extent relative to that of more stable/less flexible variants. A panel of nine Gbeta1 mutant variants was used to benchmark the screen, and spectroscopic methods were employed to characterize the thermal and structural properties of each variant accurately. The screen was combined with in silico methods to interrogate a library of randomized variants for selection of mutants of greater structural integrity.

In Vitro Evolution of a Hyperstable Gβ1 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.

Genetic algorithms as a tool for helix design – computational and experimental studies on prion protein helix 1

Evolutionary computing is a general optimization mechanism successfully implemented for a variety of numeric problems in a variety of fields, including structural biology. We here present an evolutionary approach to optimize helix stability in peptides and proteins employing the AGADIR energy function for helix stability as scoring function. With the ability to apply masks determining positions, which are to remain constant or fixed to a certain class of amino acids, our algorithm is capable of developing stable helical scaffolds containing a wide variety of structural and functional amino acid patterns. The algorithm showed good convergence behaviour in all tested cases and can be parameterized in a wide variety of ways. We have applied our algorithm for the optimization of the stability of prion protein helix 1, a structural element of the prion protein which is thought to play a crucial role in the conformational transition from the cellular to the pathogenic form of the prion protein, and which therefore poses an interesting target for pharmacological as well as genetic engineering approaches to counter the as of yet uncurable prion diseases. NMR spectroscopic investigations of selected stabilizing and destabilizing mutations found by our algorithm could demonstrate its ability to create stabilized variants of secondary structure elements.

Energy barriers, cooperativity, and hidden intermediates in the folding of small proteins

Current theoretical views of the folding process of small proteins (< approximately 100 amino acids) postulate that the landscape of potential mean force (PMF) for the formation of the native state has a funnel shape and that the free energy barrier to folding arises from the chain configurational entropy only. However, recent theoretical studies on the formation of hydrophobic clusters with explicit water suggest that a barrier should exist on the PMF of folding, consistent with the fact that protein folding generally involves a large positive activation enthalpy at room temperature. In addition, high-resolution structural studies of the hidden partially unfolded intermediates have revealed the existence of non-native interactions, suggesting that the correction of the non-native interactions during folding should also lead to barriers on PMF. To explore the effect of a PMF barrier on the folding behavior of proteins, we modified Zwanzig's model for protein folding with an uphill landscape of PMF for the formation of transition states. We found that the modified model for short peptide segments can satisfy the thermodynamic and kinetic criteria for an apparently two-state folding. Since the Levinthal paradox can be solved by a stepwise folding of short peptide segments, a landscape of PMF with a locally uphill search for the transition state and cooperative stabilization of folding intermediates/native state is able to explain the available experimental results for small proteins. We speculate that the existence of cooperative hidden folding intermediates in small proteins could be the consequence of the highly specific structures of the native state, which are selected by evolution to perform specific functions and fold in a biologically meaningful time scale.

Isolation of monomeric human V(H)s by a phage selection

Human V(H) domains are promising molecules in applications involving antibodies, in particular, immunotherapy because of their human origin. However, they are, in general, prone to aggregation. Therefore, various strategies have been employed to acquire monomeric human V(H)s. We had previously discovered that filamentous phages displaying engineered monomeric V(H) domains gave rise to significantly larger plaques on bacterial lawns than phages displaying wild type V(H)s with aggregation tendencies. Using plaque size as the selection criterion and a phage-displayed naïve human V(H) library we identified 15 V(H)s that were monomeric. Additionally, the V(H)s demonstrated good expression yields, good refolding properties following thermal denaturation, resistance to aggregation during long incubation at 37 degrees C, and to trypsin at 37 degrees C. These 15 V(H)s should serve as good scaffolds for developing immunotherapeutics, and the selection method employed here should have general utility for isolating proteins with desirable biophysical properties.

High-affinity fragment complementation of a fibronectin type III domain and its application to stability enhancement

The tenth fibronectin type III (FN3) domain of human fibronectin (FNfn10), a prototype of the ubiquitous FN3 domain, is a small, monomeric beta-sandwich protein. In this study, we have bisected FNfn10 in each loop to generate a total of six fragment pairs. We found that fragment pairs bisected at multiple loops of FNfn10 show complementation in vivo as tested with a yeast two-hybrid system. The dissociation constant of these fragment pairs determined in vitro were as low as 3 nM, resulting in one of the tightest fragment complementation systems reported so far. Furthermore, we show that the affinity of fragment complementation is correlated with the stability of the uncut parent protein. Exploring this correlation, we screened a yeast two-hybrid library of one fragment and identified mutations that suppress the effect of a destabilizing mutation in the other fragment. One of the identified mutations significantly increased the stability of the uncut wild-type protein, proving that fragment complementation can be used as a novel strategy for the selection of proteins with enhanced stability.

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.

Detection and structure determination of an equilibrium unfolding intermediate of Rd-apocytochrome b562: native fold with non-native hydrophobic interactions

The absence of detectable kinetic and equilibrium folding intermediates by optical probes is commonly taken to indicate that protein folding is a two-state process. However, for some small proteins with apparent two-state behavior, unfolding intermediates have been identified in native-state hydrogen exchange or kinetic unfolding experiments monitored by nuclear magnetic resonance. Rd-apocytochrome b(562), a four-helix bundle, is one such protein. Here, we found another unfolding intermediate for Rd-apocytochrome b(562). It is based on a cooperative transition of (15)N chemical shifts of amide protons as a function of urea concentrations before the global unfolding. We have solved the high-resolution structure of the protein at 2.8 M urea, which is after this cooperative transition but before the global unfolding. All four helices remained intact, but a number of hydrophobic core residues repacked. This intermediate provides a possible structural interpretation for the kinetic unfolding intermediates observed using nuclear magnetic resonance methods for several proteins and has important implications for theoretical studies of protein folding.

Phage-display as a tool for quantifying protein stability determinants

To address questions of protein stability, researchers have increasingly turned to combinatorial approaches that permit the rapid analysis of libraries of protein variants. Phage-display has proved to be a powerful tool for analyzing protein stability due to the large library size and the robustness of the phage particle to a variety of denaturing conditions. With the B1 domain of protein G (GB1) and a camelid heavy chain antibody as model systems, we are using phage-display libraries to experimentally address questions that have generally been addressed in silico, either through computational studies or statistical analysis of known protein structures. One effort has focused on identifying novel solutions to repacking the hydrophobic core of GB1, while maintaining stability comparable to the wild type protein. In a second study, a small set of substitutions in complimentarity-determining region 3 was found to stabilize the framework of the camelid antibody. Another major focus has been to obtain quantitative data on beta-sheet stability determinants. We have successfully adapted a phage-display method for quantitating affinities of protein variants (shotgun alanine scanning) to analysis of GB1 stability. Using this method, we have analyzed the energetic contributions of cross-strand side chain-side chain interactions. Finally, we discuss parameters to consider in using phage-display to discriminate subtle stability differences among fully folded variants. Overall, this method provides a fast approach for quantitatively addressing biophysical questions.

Combinatorial approaches to protein stability and structure

Why do proteins adopt the conformations that they do, and what determines their stabilities? While we have come to some understanding of the forces that underlie protein architecture, a precise, predictive, physicochemical explanation is still elusive. Two obstacles to addressing these questions are the unfathomable vastness of protein sequence space, and the difficulty in making direct physical measurements on large numbers of protein variants. Here, we review combinatorial methods that have been applied to problems in protein biophysics over the last 15 years. The effects of hydrophobic core composition, the most important determinant of structure and stability, are still poorly understood. Particular attention is given to core composition as addressed by library methods. Increasingly useful screens and selections, in combination with modern high-throughput approaches borrowed from genomics and proteomics efforts, are making the empirical, statistical correlation between sequence and structure a tractable problem for the coming years.