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.     

Origins of the high stability of an in vitro-selected cold-shock protein

In previous work, we had identified stabilized forms of the cold-shock protein Bs-CspB from Bacillus subtilis in a combinatorial library by an in vitro selection procedure. In this library, the sequence positions 2, 3, 46, 64, 66, and 67 had been randomized, because Bs-CspB differs from the naturally thermostable homolog Bc-Csp from Bacillus caldolyticus, among others, at these six positions. For the most stable selected variant, the midpoint of thermal unfolding (tM) increased by 28.2 deg. C and the Gibbs free energy of unfolding (deltaG(D)) by 19 kJ/mol. Here, we analyzed by site-directed mutagenesis how the selected residues contribute individually to this strong stabilization. Val3 and Val66, which replace Glu3 and Glu66 of wild-type Bs-CspB, each contribute about 7 kJ/mol to stability, the Thr64Arg substitution contributes 4.5 kJ/mol, and 3.2 kJ/mol originate from the Ala46Leu replacement. Gly67 at the carboxy terminus is unimportant for stability, the Arg selected at position 2 is overall slightly destabilizing but improves the coulombic interactions. The best variant differs from Bc-Csp at all six positions; nevertheless, natural and in vitro selection followed similar principles. In both cases, negatively charged residues at the adjacent positions 3 and 66 are avoided, and a positively charged residue is introduced into this area of the protein surface. Its exact location is unimportant. It can be at position 3, as in the thermophilic Bc-Csp, or at positions 2 or 64, as in the most stable selected variant. These positively charged residues contribute to stability not by engaging in pairwise coulombic interactions with a specific carboxyl group, but by generally improving the charge distribution in this particular region of the protein surface. These coulombic effects contribute significantly to the thermostability of the cold-shock proteins. They are only weakly interdependent and best explained by the presence of a flexible ion network at the protein surface. Our results emphasize that surface positions are very good candidates for optimizing protein stability.     

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.     

Role of the chain termini for the folding transition state of the cold shock protein.

Residues Arg3 and Leu66 are crucially important for the enhanced stability of the cold shock protein Bc-Csp from the thermophile Bacillus caldolyticus relative to its homologue Bs-CspB from the mesophile Bacillus subtilis. Arg3, which replaces Glu3 of Bs-CspB, accounts for two-thirds of the stability difference and for the entire difference in Coulombic interactions between the two proteins. Leu66, which replaces Glu66 of Bs-CspB, contributes additional hydrophobic interactions. To elucidate the role of these two residues near the chain termini for the rapid folding of the cold shock proteins, we performed an extensive mutational analysis of the folding kinetics to characterize interactions between residues 3, 46, and 66 in the transition state of folding. We employed a pressure-jump apparatus which allows folding to be followed over a broad range of temperatures and urea concentrations in the time range of microseconds to minutes. The N-terminal region folds early, and the interactions that originate from residue 3 are present to a large extent in the transition state already. They include a hydrophobic contribution, a general electrostatic stabilization by the positive charge of Arg3 in Bc-Csp, and a pairwise Coulombic repulsion with Glu46 in the Arg3Glu variant. The C-terminus appears to be largely unfolded in the transition state. The interactions of Leu66, including those with the already structured N-terminal region, are established only after passage through the transition state. The N- and C-termini of the cold shock proteins thus contribute differently to the folding kinetics, although they are very close in space in the folded protein.     

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.     

Electrostatic interactions in leucine zippers: thermodynamic analysis of the contributions of Glu and His residues and the effect of mutating salt bridges

Electrostatic interactions play a complex role in stabilizing proteins. Here, we present a rigorous thermodynamic analysis of the contribution of individual Glu and His residues to the relative pH-dependent stability of the designed disulfide-linked leucine zipper AB(SS). The contribution of an ionized side-chain to the pH-dependent stability is related to the shift of the pK(a) induced by folding of the coiled coil structure. pK(a)(F) values of ten Glu and two His side-chains in folded AB(SS) and the corresponding pK(a)(U) values in unfolded peptides with partial sequences of AB(SS) were determined by 1H NMR spectroscopy: of four Glu residues not involved in ion pairing, two are destabilizing (-5.6 kJ mol(-1)) and two are interacting with the positive alpha-helix dipoles and are thus stabilizing (+3.8 kJ mol(-1)) in charged form. The two His residues positioned in the C-terminal moiety of AB(SS) interact with the negative alpha-helix dipoles resulting in net stabilization of the coiled coil conformation carrying charged His (-2.6 kJ mol(-1)). Of the six Glu residues involved in inter-helical salt bridges, three are destabilizing and three are stabilizing in charged form, the net contribution of salt-bridged Glu side-chains being destabilizing (-1.1 kJ mol(-1)). The sum of the individual contributions of protonated Glu and His to the higher stability of AB(SS) at acidic pH (-5.4 kJ mol(-1)) agrees with the difference in stability determined by thermal unfolding at pH 8 and pH 2 (-5.3 kJ mol(-1)). To confirm salt bridge formation, the positive charge of the basic partner residue of one stabilizing and one destabilizing Glu was removed by isosteric mutations (Lys-->norleucine, Arg-->norvaline). Both mutations destabilize the coiled coil conformation at neutral pH and increase the pK(a) of the formerly ion-paired Glu side-chain, verifying the formation of a salt bridge even in the case where a charged side-chain is destabilizing. Because removing charges by a double mutation cycle mainly discloses the immediate charge-charge effect, mutational analysis tends to overestimate the overall energetic contribution of salt bridges to protein stability.     

Electrostatic stabilization of a thermophilic cold shock protein.

The cold shock protein Bc-Csp from the thermophile Bacillus caldolyticus differs from its mesophilic homolog Bs-CspB from Bacillus subtilis by 15.8 kJ mol(-1) in the Gibbs free energy of denaturation (DeltaG(D)). The two proteins vary in sequence at 12 positions but only two of them, Arg3 and Leu66 of Bc-Csp, which replace Glu3 and Glu66 of Bs-CspB, are responsible for the additional stability of Bc-Csp. These two positions are near the ends of the protein chain, but close to each other in the three-dimensional structure. The Glu3Arg exchange alone changed the stability by more than 11 kJ mol(-1). Here, we elucidated the molecular origins of the stability difference between the two proteins by a mutational analysis. Electrostatic contributions to stability were characterized by measuring the thermodynamic stabilities of many variants as a function of salt concentration. Double and triple mutant analyses indicate that the stabilization by the Glu3Arg exchange originates from three sources. Improved hydrophobic interactions of the aliphatic moiety of Arg3 contribute about 4 kJ mol(-1). Another 4 kJ mol(-1) is gained from the relief of a pairwise electrostatic repulsion between Glu3 and Glu66, as in the mesophilic protein, and 3 kJ mol(-1) originate from a general electrostatic stabilization by the positive charge of Arg3, which is not caused by a pairwise interaction. Mutations of all potential partners for an ion pair within a radius of 10 A around Arg3 had only marginal effects on stability. The Glu3-->Arg3 charge reversal thus optimizes ionic interactions at the protein surface by both local and global effects. However, it cannot convert the coulombic repulsion with another Glu residue into a corresponding attraction. Avoidance of unfavorable coulombic repulsions is probably a much simpler route to thermostability than the creation of stabilizing surface ion pairs, which can form only at the expense of conformational entropy.     

Changing the determinants of protein stability from covalent to non-covalent interactions by in vitro evolution: a structural and energetic analysis

The three disulfide bonds of the gene-3-protein of the phage fd are essential for the conformational stability of this protein, and it unfolds when they are removed by reduction or mutation. Previously, we used an iterative in vitro selection strategy to generate a stable and functional form of the gene-3-protein without these disulfides. It yielded optimal replacements for the disulfide bonds as well as several stabilizing second-site mutations. The best selected variant showed a higher thermal stability compared with the disulfide-bonded wild-type protein. Here, we investigated the molecular basis of this strong stabilization by solving the crystal structure of this variant and by analyzing the contributions to the conformational stability of the selected mutations individually. They could mostly be explained by improved side-chain packing. The R29W substitution alone increased the midpoint of the thermal unfolding transition by 14 deg and the conformational stability by about 25 kJ mol^− 1. This key mutation (i) removed a charged side chain that forms a buried salt bridge in the disulfide-containing wild-type protein, (ii) optimized the local packing with the residues that replace the C46-C53 disulfide and (iii) improved the domain interactions. Apparently, certain residues in proteins indeed play key roles for stability.     

Role of electrostatic interactions for the stability and folding behavior of cold shock protein

The impacts of three charged‐residue‐involved mutations, E46A, R3E, and R3E/L66E, on the thermostability and folding behavior of the cold shock protein from the themophile Bacillus caldolyticus (Bc‐Csp) were investigated by using a modified Gō‐like model, in which the nonspecific electrostatic interactions of charged residues were taken into account. Our simulation results show that the wild‐type Bc‐Csp and its three mutants are all two‐sate folders, which is consistent with the experimental observations. It is found that these three mutations all lead to a decrease of protein thermodynamical stability, and the effect of R3E mutation is the strongest. The lower stability of these three mutants is due to the increase of the enthalpy of the folded state and the entropy of the unfolded state. Using this model, we also studied the folding kinetics and the folding/unfolding pathway of the wild‐type Bc‐Csp as well as its three mutants and then discussed the effects of electrostatic interactions on the folding kinetics. The results indicate that the substitutions at positions 3 and 46 largely decrease the folding kinetics, whereas the mutation of residue 66 only slightly decreases the folding rate. This result agrees well with the experimental observations. It is also found that these mutations have little effects on the folding transition state and the folding pathway, in which the N‐terminal β sheet folds earlier than the C‐terminal region. We also investigated the detailed unfolding pathway and found that it is really the reverse of the folding pathway, providing the validity of our simulation results. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

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.     

Contribution of surface salt bridges to protein stability: guidelines for protein engineering

The small globular protein, ubiquitin, contains a pair of oppositely charged residues, K11 and E34, that according to the three-dimensional structure are located on the surface of this protein with a spatial orientation characteristic of a salt bridge. We investigated the strength of this salt bridge and its contribution to the global stability of the ubiquitin molecule. Using the "double mutant cycle" analysis, the strength of the pairwise interactions between K11 and E34 was estimated to be favorable by 3.6kJ/mol. Further, the salt bridge of the reverse orientation, i.e. E11/K34, can be formed and is found to have a strength (3.8kJ/mol) similar to that of the K11/E34 pair. However, the global stability of the K11/E34 variant of ubiquitin is 2.2kJ/mol higher than that of the E11/K34 variant. The difference in the contribution of the opposing salt bridge orientations to the overall stability of the ubiquitin molecule is attributed to the difference in the charge-charge interactions between residues forming the salt bridge and the rest of the ionizable groups in this protein. On the basis of these results, we concluded that surface salt bridges are stabilizing, but their contribution to the overall protein stability is strongly context-dependent, with charge-charge interactions being the largest determinant. Analysis of 16 salt bridges from six different proteins, for which detailed experimental data on energetics have been reported, support the conclusions made from the analysis of the salt bridge in ubiquitin. Implications of these findings for engineering proteins with enhanced thermostability are discussed.     

The role of loop 6/7 in folding and functional performance of Na,K-ATPase

Alanine substitutions were made for 15 amino acids in the cytoplasmic loop between transmembrane helices 6 and 7 (L6/7) of the human alpha(1)-subunit of Na,K-ATPase. Most mutations reduced Na,K-ATPase activity by less than 50%; however, the mutations R834A, R837A, and R848A reduced Na,K-ATPase activity by 75, 89, and 66%, respectively. Steady-state phosphoenzyme formation from ATP was reduced in mutants R834A, R837A, and R848A, and R837A also had a faster E(2)P --> E(2) dephosphorylation rate compared with the wild-type enzyme. Effects of L6/7 mutations on the phosphorylation domain of the protein were also demonstrated by (18)O exchange, which showed that intrinsic rate constants for P(i) binding and/or reaction with the protein were altered. Although most L6/7 mutations had no effect on the interaction of Na(+) or K(+) with Na,K-ATPase, the E825A, E828A, R834A, and R837A mutations reduced the apparent affinity of the enzyme for both Na(+) and K(+) by 1.5-3-fold. 1-Bromo-2,4,6-tris(methylisothiouronium)benzene (Br-TITU(3+)), a competitive antagonist of Rb(+) and Na(+) occlusion, was used to test whether charged residues in L6/7 are involved in binding monovalent cations and cation antagonists. Br-TITU(3+) inhibited ouabain binding to wild type Na,K-ATPase with an IC(50) of 30 microM. Ouabain binding to the E825A, E828A, R834A, or R837A mutants was still inhibited by Br-TITU(3+), indicating that Br-TITU(3+) does not bind to charged residues in L6/7. This observation makes it unlikely that L6/7 functions as a cytoplasmic cation binding site in Na,K-ATPase, and together with the effects of L6/7 mutations on phosphate interactions with the enzyme suggests that L6/7 is important in stabilizing the phosphorylation domain and its relationship to the ion binding sites of the protein.     

Electrostatic interactions between transmembrane segments mediate folding of Shaker K+ channel subunits.

In voltage-dependent Shaker K+ channels, charged residues E293 in transmembrane segment S2 and R365, R368, and R371 in S4 contribute significantly to the gating charge movement that accompanies activation. Using an intragenic suppression strategy, we have now probed for structural interaction between transmembrane segments S2, S3, and S4 in Shaker channels. Charge reversal mutations of E283 in S2 and K374 in S4 disrupt maturation of the protein. Maturation was specifically and efficiently rescued by second-site charge reversal mutations, indicating that electrostatic interactions exist between E283 in S2 and R368 and R371 in S4, and between K374 in S4 and E293 in S2 and D316 in S3. Rescued subunits were incorporated into functional channels, demonstrating that a native structure was restored. Our data indicate that K374 interacts with E293 and D316 within the same subunit. These electrostatic interactions mediate the proper folding of the protein and are likely to persist in the native structure. Our results raise the possibility that the S4 segment is tilted relative to S2 and S3 in the voltage-sensing domain of Shaker channels. Such an arrangement might provide solvent access to voltage-sensing residues, which we find to be highly tolerant of mutations.     

The systematic substitutions around the conserved charged residues of the cytoplasmic loop of Na+-driven flagellar motor component PomA

PomA, a homolog of MotA in the H+-driven flagellar motor, is an essential component for torque generation in the Na+-driven flagellar motor. Previous studies suggested that two charged residues, R90 and E98, which are in the single cytoplasmic loop of MotA, are directly involved in this process. These residues are conserved in PomA of Vibrio alginolyticus as R88 and E96, respectively. To explore the role of these charged residues in the Na+-driven motor, we replaced them with other amino acids. However, unlike in the H+-driven motor, both of the single and the double PomA mutants were functional. Several other positively and negatively charged residues near R88 and E96, namely K89, E97 and E99, were neutralized. Motility was retained in a strain producing the R88A/K89A/E96Q/E97Q/E99Q (AAQQQ) PomA protein. The swimming speed of the AAQQQ strain was as fast as that of the wild-type PomA strain, but the direction of motor rotation was abnormally counterclockwise-biased. We could, however, isolate non-motile or poorly motile mutants when certain charged residues in PomA were reversed or neutralized. The charged residues at positions 88-99 of PomA may not be essential for torque generation in the Na+-driven motor and might play a role in motor function different from that of the equivalent residues of the H+-driven motor.     

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.     

Electrostatic contributions to the stability of a thermophilic cold shock protein

The thermophilic Bacillus caldolyticus cold shock protein (Bc-Csp) differs from the mesophilic Bacillus subtilis cold shock protein B (Bs-CspB) in 11 of the 66 residues. Stability measurements of Schmid and co-workers have implicated contributions of electrostatic interactions to the thermostability. To further elucidate the physical basis of the difference in stability, previously developed theoretical methods that treat electrostatic effects in both the folded and the unfolded states were used in this paper to study the effects of mutations, ionic strength, and temperature. For 27 mutations that narrow the difference in sequence between Bc-Csp and Bs-CspB, calculated changes in unfolding free energy (Delta G) and experimental results have a correlation coefficient of 0.98. Bc-Csp appears to use destabilization of the unfolded state by unfavorable charge-charge interactions as a mechanism for increasing stability. Accounting for the effects of ionic strength and temperature on the electrostatic free energies in both the folded and the unfolded states, explanations for two important experimental observations are presented. The disparate ionic strength dependences of Delta G for Bc-Csp and Bs-CspB were attributed to the difference in the total charges (-2e and -6e, respectively). A main contribution to the much higher unfolding entropy of Bs-CspB was found to come from the less favorable electrostatic interactions in the folded state. These results should provide insight for understanding the thermostability of other thermophilic proteins.     

Three charged amino acids in extracellular loop 1 are involved in maintaining the outer pore architecture of CFTR

The cystic fibrosis (CF) transmembrane conductance regulator (CFTR) bears six extracellular loops (ECL1–6); ECL1 is the site of several mutations associated with CF. Mutation R117H has been reported to reduce current amplitude, whereas D110H, E116K, and R117C/L/P may impair channel stability. We hypothesized that these amino acids might not be directly involved in ion conduction and permeation but may contribute to stabilizing the outer vestibule architecture in CFTR. We used cRNA injected oocytes combined with electrophysiological techniques to test this hypothesis. Mutants bearing cysteine at these sites were not functionally modified by extracellular MTS reagents and were blocked by GlyH-101 similarly to WT-CFTR. These results suggest that these three residues do not contribute directly to permeation in CFTR. In contrast, mutants D110R-, E116R-, and R117A-CFTR exhibited instability of the open state and significantly shortened burst duration compared with WT-CFTR and failed to be locked into the open state by AMP-PNP (adenosine 5′-(β,γ-imido) triphosphate); charge-retaining mutants showed mainly the full open state with comparably longer open burst duration. These interactions suggest that these ECL1 residues might be involved in maintaining the outer pore architecture of CFTR. A CFTR homology model suggested that E116 interacts with R104 in both the closed and open states, D110 interacts with K892 in the fully closed state, and R117 interacts with E1126 in the open state. These interactions were confirmed experimentally. The results suggest that D110, E116, and R117 may contribute to stabilizing the architecture of the outer pore of CFTR by interactions with other charged residues.     

Electrostatic forces involved in orienting Anabaena ferredoxin during binding to Anabaena ferredoxin:NADP+ reductase: site-specific mutagenesis, transient kinetic measurements, and electrostatic surface potentials.

Transient absorbance measurements following laser flash photolysis have been used to measure the rate constants for electron transfer (et) from reduced Anabaena ferredoxin (Fd) to wild-type and seven site-specific charge-reversal mutants of Anabaena ferredoxin:NADP+ reductase (FNR). These mutations have been designed to probe the importance of specific positively charged amino acid residues on the surface of the FNR molecule near the exposed edge of the FAD cofactor in the protein-protein interaction during et with Fd. The mutant proteins fall into two groups: overall, the K75E, R16E, and K72E mutants are most severely impaired in et, and the K138E, R264E, K290E, and K294E mutants are impaired to a lesser extent, although the degree of impairment varies with ionic strength. Binding constants for complex formation between the oxidized proteins and for the transient et complexes show that the severity of the alterations in et kinetics for the mutants correlate with decreased stabilities of the protein-protein complexes. Those mutated residues, which show the largest effects, are located in a region of the protein in which positive charge predominates, and charge reversals have large effects on the calculated local surface electrostatic potential. In contrast, K138, R264, K290, and K294 are located within or close to regions of intense negative potential, and therefore the introduction of additional negative charges have considerably smaller effects on the calculated surface potential. We attribute the relative changes in et kinetics and complex binding constants for these mutants to these characteristics of the surface charge distribution in FNR and conclude that the positively charged region of the FNR surface located in the vicinity of K75, R16, and K72 is especially important in the binding and orientation of Fd during electron transfer.     

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.     

残基突变对冷休克蛋白热稳定性影响的计算分析研究(英文)

残基突变是提高蛋白质热稳定性最直接有效的方式。在本文中,我们选取一对冷休克蛋白质作为研究对象,其中一个来自嗜温的Bacillus subtilis(Bs-CspB),另一个来自嗜热的Bacillus caldolyticus(Bc-Csp),这两个蛋白质在序列和结构上具有高度的相似性,但两者的耐热能力却相差很大。我们利用全原子模型计算残基突变前后蛋白质的自由能和氨基酸之间相互作用能的变化,分析残基突变对冷休克蛋白热稳定性的影响。通过对比两个蛋白质对应位置上残基的能量,我们成功鉴别出对Bc-Csp的高热稳定性有突出贡献的残基。我们计算了这些残基突变前后,该残基的静电相互作用和范德华相互作用的变化,以分析该残基对Bc-Csp高热稳定性的主要贡献。同时,我们分析了离子键对蛋白质热稳定性的贡献。我们的计算结果和实验结果吻合得很好,关键在于利用该方法可以详细地说明残基突变影响蛋白质热稳定性的根本原因。本文为研究残基突变对蛋白质热稳定性的影响提供了一种计算思路和方法,并有助于设计具有高耐热能力的蛋白质。