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.     

Stabilization of the cold shock protein CspB from Bacillus subtilis by evolutionary optimization of Coulombic interactions

The bacterial cold shock proteins (Csp) are used by both experimentalists and theoreticians as model systems for analyzing the Coulombic contributions to protein stability. We employ Proside, a method of directed evolution, to identify stabilized variants of Bs-CspB from Bacillus subtilis. Proside links the increased protease resistance of stabilized protein variants to the infectivity of a filamentous phage. Here, three cspB libraries were used for in vitro selections to explore the stabilizing potential of charged amino acids in Bs-CspB. In the first library codons for nine selected surface residues were partially randomized, in the second one random mutations were introduced non-specifically by error-prone PCR, and in the third one the spontaneous mutation rate of the phage in Escherichia coli was used. Stabilizing mutations were found at the surface positions 1, 3, 46, 48, 65, and 66. The contributions of these mutations to stability were characterized by analyzing them individually and in combination. The best combination (M1R, E3K, K65I, and E66L) increased the midpoint of thermal unfolding of Bs-CspB from 53.8 to 85.0 degrees C. The effects of most mutations are strongly context dependent. A good example is provided by the E3R mutation. It is strongly stabilizing (DeltaDeltaGD=11.1kJ mol(-1)) in the wild-type protein, but destabilizing (DeltaDeltaGD=-4.0kJ mol(-1)) in the A46K/S48R/E66L variant. The stabilizations by charge mutations did not correlate well with the corresponding changes in the protein net charge, and they could not be ascribed to the formation of ion pairs. Previous theoretical analyses did not identify the stabilization caused by the mutations at positions 1, 46, and 48. Also, electrostatics calculations based on protein net charge or charge asymmetry did not predict well the stability changes that occur when charged residues in Bs-CspB are mutated. It remains a challenge to model the Coulombic interactions of charged residues in a protein and to determine their contributions to the Gibbs free energy of protein folding.     

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.     

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.     

In-vitro selection of highly stabilized protein variants with optimized surface

Thermostable proteins are of prime importance in protein science, but it has remained difficult to develop general strategies for stabilizing a protein. Site-directed mutagenesis based on comparisons with thermophilic homologs is rarely successful because the sequence differences are too numerous and dominated by neutral mutations. Here we used a method of directed evolution to increase the stability of a mesophilic protein, the cold shock protein Bs-CspB from Bacillus subtilis. It differs from its thermophilic counterpart Bc-Csp from Bacillus caldolyticus at 12 surface-exposed positions. To elucidate the stabilizing potential of exposed amino acid residues, six of these variant positions were randomized by saturation mutagenesis, the corresponding library of sequences was inserted into the gene-3-protein of the filamentous phage fd, and stabilized variants were selected by the Proside technique. Proside links the increased protease resistance of stabilized protein variants with the infectivity of the phage. Many strongly stabilized variants of Bs-CspB were identified in two selections, one in the presence of a denaturant and the other at elevated temperature. Several of them are significantly more stable than the naturally thermostable homolog Bc-Csp, and the best variant reaches Tm-Csp (the homolog from the hyperthermophile Thermotoga maritima) in stability. Remarkably, this variant differs from Tm-Csp at five and from Bc-Csp at all six randomized positions. This indicates that proteins can be strongly stabilized by many different sets of surface mutations, and Proside selects them efficiently from large libraries. The course of the selection could be directed by the conditions. In an ionic denaturant non-polar surface interactions were optimized, whereas at elevated temperature variants with improved electrostatics were selected, pointing to two different strategies for stabilization at protein surfaces.     

Structure of a protein G helix variant suggests the importance of helix propensity and helix dipole interactions in protein design

Six helix surface positions of protein G (Gbeta1) were redesigned using a computational protein design algorithm, resulting in the five fold mutant Gbeta1m2. Gbeta1m2 is well folded with a circular dichroism spectrum nearly identical to that of Gbeta1, and a melting temperature of 91 degrees C, approximately 6 degrees C higher than that of Gbeta1. The crystal structure of Gbeta1m2 was solved to 2.0 A resolution by molecular replacement. The absence of hydrogen bond or salt bridge interactions between the designed residues in Gbeta1m2 suggests that the increased stability of Gbeta1m2 is due to increased helix propensity and more favorable helix dipole interactions.     

Evolutionary stabilization of the gene-3-protein of phage fd reveals the principles that govern the thermodynamic stability of two-domain proteins

The gene-3-protein (G3P) of filamentous phage is essential for their propagation. It consists of three domains. The CT domain anchors G3P in the phage coat, the N2 domain binds to the F pilus of Escherichia coli and thus initiates infection, and the N1 domain continues by interacting with the TolA receptor. Phage are thus only infective when the three domains of G3P are tightly linked, and this requirement is exploited by Proside, an in vitro selection method for proteins with increased stability. In Proside, a repertoire of variants of the protein to be stabilized is inserted between the N2 and the CT domains of G3P. Stabilized variants can be selected because they resist cleavage by a protease and thus maintain the essential linkage between the domains. The method is limited by the proteolytic stability of G3P itself. We improved the stability of G3P by subjecting the phage without a guest protein to rounds of random in vivo mutagenesis and proteolytic Proside selections. Variants of G3P with one to four mutations were selected, and the temperature at which the corresponding phage became accessible for a protease increased in a stepwise manner from 40 degrees C to almost 60 degrees C. The N1-N2 fragments of wild-type gene-3-protein and of the four selected variants were purified and their stabilities towards thermal and denaturant-induced unfolding were determined. In the biphasic transitions of these proteins domain dissociation and unfolding of N2 occur in a concerted reaction in the first step, followed by the independent unfolding of domain N1 in the second step. N2 is thus less stable than N1, and it unfolds when the interactions with N1 are broken. The strongest stabilizations were caused by mutations in domain N2, in particular in its hinge subdomain, which provides many stabilizing interactions between the N1 and N2 domains. These results reveal how the individual domains and their assembly contribute to the overall stability of two-domain proteins and how mutations are optimally placed to improve the stability of such proteins.     

Structure‐based design of combinatorial mutagenesis libraries

The development of protein variants with improved properties (thermostability, binding affinity, catalytic activity, etc.) has greatly benefited from the application of high‐throughput screens evaluating large, diverse combinatorial libraries. At the same time, since only a very limited portion of sequence space can be experimentally constructed and tested, an attractive possibility is to use computational protein design to focus libraries on a productive portion of the space. We present a general‐purpose method, called “Structure‐based Optimization of Combinatorial Mutagenesis ” (SOCoM ), which can optimize arbitrarily large combinatorial mutagenesis libraries directly based on structural energies of their constituents. SOCoM chooses both positions and substitutions, employing a combinatorial optimization framework based on library‐averaged energy potentials in order to avoid explicitly modeling every variant in every possible library. In case study applications to green fluorescent protein, β‐lactamase, and lipase A, SOCoM optimizes relatively small, focused libraries whose variants achieve energies comparable to or better than previous library design efforts, as well as larger libraries (previously not designable by structure‐based methods) whose variants cover greater diversity while still maintaining substantially better energies than would be achieved by representative random library approaches. By allowing the creation of large‐scale combinatorial libraries based on structural calculations, SOCoM promises to increase the scope of applicability of computational protein design and improve the hit rate of discovering beneficial variants. While designs presented here focus on variant stability (predicted by total energy), SOCoM can readily incorporate other structure‐based assessments, such as the energy gap between alternative conformational or bound states.     

High-affinity binders selected from designed ankyrin repeat protein libraries

We report here the evolution of ankyrin repeat (AR) proteins in vitro for specific, high-affinity target binding. Using a consensus design strategy, we generated combinatorial libraries of AR proteins of varying repeat numbers with diversified binding surfaces. Libraries of two and three repeats, flanked by 'capping repeats,' were used in ribosome-display selections against maltose binding protein (MBP) and two eukaryotic kinases. We rapidly enriched target-specific binders with affinities in the low nanomolar range and determined the crystal structure of one of the selected AR proteins in complex with MBP at 2.3 A resolution. The interaction relies on the randomized positions of the designed AR protein and is comparable to natural, heterodimeric protein-protein interactions. Thus, our AR protein libraries are valuable sources for binding molecules and, because of the very favorable biophysical properties of the designed AR proteins, an attractive alternative to antibody libraries.     

Electrostatics significantly affect the stability of designed homeodomain variants

The role of electrostatic interactions in determining the stability of designed proteins was studied by constructing and analyzing a set of designed variants of the Drosophila engrailed homeodomain. Computational redesign of 29 surface positions results in a 25-fold mutant with moderate stability, similar to the wild-type protein. Incorporating helix dipole and N-capping considerations into the design algorithm by restricting amino acid composition at the helix termini and N-capping positions yields a ninefold mutant of the initial design (a 23-fold mutant of wild-type) that is over 3 kcal mol(-1) more stable than the protein resulting from the unbiased design. Four additional proteins were constructed and analyzed to isolate the effects of helix dipole and N-capping interactions in each helix. Based on the results of urea-denaturation experiments and calculations using the finite difference Poisson-Boltzmann method, both classes of interaction are found to increase the stability of the designed proteins significantly. The simple electrostatic model used in the optimization of rotamers by iterative techniques (ORBIT) force-field, which is similar to the electrostatic models used in other protein design force-fields, is unable to predict the experimentally determined stabilities of the designed variants. The helix dipole and N-capping restrictions provide a simple but effective method to incorporate two types of electrostatic interactions that impact protein stability significantly.     

A stable disulfide-free gene-3-protein of phage fd generated by in vitro evolution

Disulfide bonds provide major contributions to the conformational stability of proteins, and their cleavage often leads to unfolding. The gene-3-protein of the filamentous phage fd contains two disulfides in its N1 domain and one in its N2 domain, and these three disulfide bonds are essential for the stability of this protein. Here, we employed in vitro evolution to generate a disulfide-free variant of the N1-N2 protein with a high conformational stability. The gene-3-protein is essential for the phage infectivity, and we exploited this requirement for a proteolytic selection of stabilized protein variants from phage libraries. First, optimal replacements for individual disulfide bonds were identified in libraries, in which the corresponding cysteine codons were randomized. Then stabilizing amino acid replacements at non-cysteine positions were selected from libraries that were created by error-prone PCR. This stepwise procedure led to variants of N1-N2 that are devoid of all three disulfide bonds but stable and functional. The best variant without disulfide bonds showed a much higher conformational stability than the disulfide-containing wild-type form of the gene-3-protein. Despite the loss of all three disulfide bonds, the midpoints of the thermal transitions were increased from 48.5 degrees C to 67.0 degrees C for the N2 domain and from 60.0 degrees C to 78.7 degrees C for the N1 domain. The major loss in conformational stability caused by the removal of the disulfides was thus over-compensated by strongly improved non-covalent interactions. The stabilized variants were less infectious than the wild-type protein, probably because the domain mobility was reduced. Only a small fraction of the sequence space could be accessed by using libraries created by error-prone PCR, but still many strongly stabilized variants could be identified. This is encouraging and indicates that proteins can be stabilized by mutations in many different ways.     

Sequence plasticity in the antigen-binding site of a therapeutic anti-HER2 antibody

We have examined the plasticity of the antigen-combining site of a high-affinity antibody. In phage-displayed Fab libraries, selected CDR positions and one FR position of the humanized anti-Her2 antibody hu4D5 were substituted with all 20 amino acids. Antigen-binding selections were used to enrich for high-affinity variants, and a large number of sequences were obtained prior to convergence of the selected pool to a small set of clones. As expected, sequence variability of the antigen-binding site is overall diminished compared to known IgG sequences; however, certain positions retain much higher variability than others. The sequence variability map of the hu4D5 binding site is compared with a map derived from previous alanine-scanning of the antibody. Affinities of soluble Fab fragments for antigen confirm that multiple variants were selected with high affinity for antigen, including one variant with a single point mutation that was about threefold improved in affinity compared to the parental hu4D5. Interestingly, this mutation is one of the most radical in terms of changing side-chain chemistry (Trp for Asp) and occurs at the most plastic site as calculated by the Wu-Kabat variability coefficient. Thus variability mapping yields information about the antibody-antigen interaction that is useful and complementary to that obtained by alanine scanning mutagenesis.     

Analysis of the context dependent sequence requirements of active site residues in the metallo-beta-lactamase IMP-1

The metallo-beta-lactamase IMP-1 catalyzes the hydrolysis of a broad range of beta-lactam antibiotics to provide bacterial resistance to these compounds. In this study, 29 amino acid residue positions in and near the active-site pocket of the IMP-1 enzyme were randomized individually by site-directed mutagenesis of the corresponding codons in the bla(IMP-1) gene. The 29 random libraries were used to identify positions that are critical for the catalytic and substrate-specific properties of the IMP-1 enzyme. Mutants from each of the random libraries were selected for the ability to confer to Escherichia coli resistance to ampicillin, cefotaxime, imipenem or cephaloridine. The DNA sequence of several functional mutants was determined for each of the substrates. Comparison of the sequences of mutants obtained from the different antibiotic selections indicates the sequence requirements for each position in the context of each substrate. The zinc-chelating residues in the active site were found to be essential for hydrolysis of all antibiotics tested. Several positions, however, displayed context-dependent sequence requirements, in that they were essential for one substrate(s) but not others. The most striking examples included Lys69, Asp84, Lys224, Pro225, Gly232, Asn233, Asp236 and Ser262. In addition, comparison of the results for all 29 positions indicates that hydrolysis of imipenem, cephaloridine and ampicillin has stringent sequence requirements, while the requirements for hydrolysis of cefotaxime are more relaxed. This suggests that more information is required to specify active-site pockets that carry out imipenem, cephaloridine or ampicillin hydrolysis than one that catalyzes cefotaxime hydrolysis.     

G protein β1γ2 subunits purification and their interaction with adenylyl cyclase

A preliminary study on the interaction of G protein (guanine triphosphate binding protein) b1g2 subunits and their coupled components in cell signal transduction was conducted in vitro. The insect cell lines, Sf9 (Spodoptera frugiperda) and H5 (Trichoplusia ni) were used to express the recombinant protein Gβ₁γ₂. The cell membrane containing Gβ₁γ₂ was isolated through affinity chromatography column with Ni-NTA agarose by FPLC method, and the highly purified protein was obtained. The adenylyl cyclase 2 (AC2) activity assay showed that the purified Gβ₁γ₂ could significantly stimulate AC₂ activity. The interaction of β₁γ₂ subunits of G protein with the cytoplasmic tail of various mammalian adenylyl cyclases was monitored by BIAcore technology using NTA sensor chip, which relies on the phenomenon of surface plasmon resonance (SPR). The experiments showed the direct binding of Gβ₁γ₂ to the cytoplasmic tail C2 domain of AC2. The specific binding domain of AC2 with Gβ₁γ₂ was the same as AC2 activity domain which was stimulated by Gβ₁γ₂.     

Designing specific protein-protein interactions using computation, experimental library screening, or integrated methods

Given the importance of protein-protein interactions for nearly all biological processes, the design of protein affinity reagents for use in research, diagnosis or therapy is an important endeavor. Engineered proteins would ideally have high specificities for their intended targets, but achieving interaction specificity by design can be challenging. There are two major approaches to protein design or redesign. Most commonly, proteins and peptides are engineered using experimental library screening and/or in vitro evolution. An alternative approach involves using protein structure and computational modeling to rationally choose sequences predicted to have desirable properties. Computational design has successfully produced novel proteins with enhanced stability, desired interactions and enzymatic function. Here we review the strengths and limitations of experimental library screening and computational structure-based design, giving examples where these methods have been applied to designing protein interaction specificity. We highlight recent studies that demonstrate strategies for combining computational modeling with library screening. The computational methods provide focused libraries predicted to be enriched in sequences with the properties of interest. Such integrated approaches represent a promising way to increase the efficiency of protein design and to engineer complex functionality such as interaction specificity.     

Optimization of combinatorial mutagenesis

Protein engineering by combinatorial site-directed mutagenesis evaluates a portion of the sequence space near a target protein, seeking variants with improved properties (e.g., stability, activity, immunogenicity). In order to improve the hit-rate of beneficial variants in such mutagenesis libraries, we develop methods to select optimal positions and corresponding sets of the mutations that will be used, in all combinations, in constructing a library for experimental evaluation. Our approach, OCoM (Optimization of Combinatorial Mutagenesis), encompasses both degenerate oligonucleotides and specified point mutations, and can be directed accordingly by requirements of experimental cost and library size. It evaluates the quality of the resulting library by one- and two-body sequence potentials, averaged over the variants. To ensure that it is not simply recapitulating extant sequences, it balances the quality of a library with an explicit evaluation of the novelty of its members. We show that, despite dealing with a combinatorial set of variants, in our approach the resulting library optimization problem is actually isomorphic to single-variant optimization. By the same token, this means that the two-body sequence potential results in an NP-hard optimization problem. We present an efficient dynamic programming algorithm for the one-body case and a practically-efficient integer programming approach for the general two-body case. We demonstrate the effectiveness of our approach in designing libraries for three different case study proteins targeted by previous combinatorial libraries--a green fluorescent protein, a cytochrome P450, and a beta lactamase. We found that OCoM worked quite efficiently in practice, requiring only 1 hour even for the massive design problem of selecting 18 mutations to generate 10? variants of a 443-residue P450. We demonstrate the general ability of OCoM in enabling the protein engineer to explore and evaluate trade-offs between quality and novelty as well as library construction technique, and identify optimal libraries for experimental evaluation.     

白介素18结合蛋白研究进展

白介素18结合蛋白（IL-18BP）是新近发现的一种糖蛋白,属免疫球蛋白超家族成员,目前已发现有6种IL-18BP的同工蛋白,IL-18BP可以在体内外有效抑制IL-18的作用。因而被认为是IL-18的天然拮抗剂,另外发现几种痘病毒编码的蛋白质与IL-18BP有高度同源性,其病毒产物可减弱IL-18诱导的Th1反应,利用IL-18BP拮抗IL-18的作用进行基因治疗将为某些自身免疫性疾病的治疗开辟一条新的途径。     

Accurate computer-based design of a new backbone conformation in the second turn of protein L.

The rational design of loops and turns is a key step towards creating proteins with new functions. We used a computational design procedure to create new backbone conformations in the second turn of protein L. The Protein Data Bank was searched for alternative turn conformations, and sequences optimal for these turns in the context of protein L were identified using a Monte Carlo search procedure and an energy function that favors close packing. Two variants containing 12 and 14 mutations were found to be as stable as wild-type protein L. The crystal structure of one of the variants has been solved at a resolution of 1.9 A, and the backbone conformation in the second turn is remarkably close to that of the in silico model (1.1 A RMSD) while it differs significantly from that of wild-type protein L (the turn residues are displaced by an average of 7.2 A). The folding rates of the redesigned proteins are greater than that of the wild-type protein and in contrast to wild-type protein L the second beta-turn appears to be formed at the rate limiting step in folding.     

Coordinate expression of the alpha and beta subunits of heterotrimeric G proteins involves regulation of protein degradation in CHO cells

Individual cell types express a characteristic balance between heterotrimeric G protein alpha and betagamma subunits, but little is known about the regulatory mechanism. We systemically examined the regulatory mechanism in CHO cells. We found that expression of Galphas, Galphai2, and Galphaq proteins increased in direct proportion to the increase of Gbeta1gamma2 overexpressed transiently. Expression of Gbeta protein also increased following overexpression of Galphas, Galphai2, and Galphaq. The Gbetagamma overexpression stimulated degradation of Gbeta in contrast to reduction of Galphas degradation. We conclude that coordinate expression of the G protein subunits involves regulation of protein degradation via proteasome in CHO cells.