Search

To initiate infection of Escherichia coli, phage fd uses its gene-3-protein (G3P) to bind first to an F pilus and then to the TolA protein at the cell surface. G3P is normally auto-inhibited because a tight interaction between the two N-terminal domains N1 and N2 buries the TolA binding site. Binding of N2 to the pilus activates G3P by initiating long-range conformational changes that are relayed to the domain interface and to a proline timer. We discovered that the 23-28 loop of the N1 domain is critical for propagating these conformational signals. The analysis of the stability and the folding dynamics of G3P variants with a shortened loop combined with TolA interaction studies and phage infection experiments reveal how the contact between the N2 domain and the 23-28 loop of N1 is energetically linked with the interdomain region and the proline timer and how it affects phage infectivity. Our results illustrate how conformational transitions and prolyl cis/trans isomerization can be coupled energetically and how conformational signals to and from prolines can be propagated over long distances in proteins.

Parvulins are small prolyl isomerases and serve as catalytic domains of folding enzymes. SurA (survival protein A) from the periplasm of Escherichia coli consists of an inactive (Par1) and an active (Par2) parvulin domain as well as a chaperone domain. In the absence of the chaperone domain, the folding activity of Par2 is virtually abolished. We created a chimeric protein by inserting the chaperone domain of SlyD, an unrelated folding enzyme from the FKBP family, into a loop of the isolated Par2 domain of SurA. This increased its folding activity 450-fold to a value higher than the activity of SurA, in which Par2 is linked with its natural chaperone domain. In the presence of both the natural and the foreign chaperone domain, the folding activity of Par2 was 1500-fold increased. Related and unrelated chaperone domains thus are similarly efficient in enhancing the folding activity of the prolyl isomerase Par2. A sequence analysis of various chaperone domains suggests that clusters of exposed methionine residues in mobile chain regions might be important for a generic interaction with unfolded protein chains. This binding is highly dynamic to allow frequent transfer of folding protein chains between chaperone and catalytic domains.

Folding enzymes often use distinct domains for the binding of substrate proteins (“chaperone domains”) and for the catalysis of slow folding reactions such as disulfide formation or prolyl isomerization. The human prolyl isomerase FKBP12 is a small single-domain protein without a chaperone domain. Its very low folding activity could previously be increased by inserting the chaperone domain from the homolog SlyD (sensitive-to-lysis protein D) of Escherichia coli. We now inserted three unrelated chaperone domains into human FKBP12: the apical domain of the chaperonin GroEL from E. coli, the chaperone domain of protein disulfide isomerase from yeast, or the chaperone domain of SurA from the periplasm of E. coli. All three conveyed FKBP12 with a high affinity for unfolded proteins and increased its folding activity. Substrate binding and release of the chimeric folding enzymes were found to be very fast. This allows rapid substrate transfer from the chaperone domain to the catalytic domain and ensures efficient rebinding of protein chains that were unable to complete folding. The advantage of having separate sites, first for generic protein binding and then for specific catalysis, explains why our construction of the artificial folding enzymes with foreign chaperone domains was successful.

The filamentous phage fd uses its gene 3 protein (G3P) to target Escherichia coli cells in a two-step process. First, the N2 domain of G3P attaches to an F pilus, and then the N1 domain binds to TolA-C. N1 and N2 are tightly associated, rendering the phage robust but noninfectious because the binding site for TolA-C is buried at the domain interface. Binding of N2 to the F pilus initiates partial unfolding, domain disassembly, and prolyl cis-to-trans isomerization in the hinge between N1 and N2. This activates the phage, and trans-Pro213 maintains this state long enough for N1 to reach TolA-C. Phage IF1 targets I pili, and its G3P contains also an N1 domain and an N2 domain. The pilus-binding N2 domains of the phages IF1 and fd are unrelated, and the N1 domains share a 31% sequence identity. We show that N2 of phage IF1 mediates binding to the I pilus, and that N1 targets TolA. Crystallographic and NMR analyses of the complex between N1 and TolA-C indicate that phage IF1 interacts with the same site on TolA-C as phage fd. In IF1-G3P, N1 and N2 are independently folding units, however, and the TolA binding site on N1 is permanently accessible. Activation by unfolding and prolyl isomerization, as in the case of phage fd, is not observed. In IF1-G3P, the absence of stabilizing domain interactions is compensated for by a strong increase in the stabilities of the individual domains. Apparently, these closely related filamentous phages evolved different mechanisms to reconcile robustness with high infectivity.

In the N2 domain of the gene-3-protein of phage fd, two consecutive β-strands are connected by a mobile loop of seven residues (157–163). The stability of this loop is low, and the Asp160-Pro161 bond at its tip shows conformational heterogeneity with 90% being in the cis and 10% in the trans form. The refolding kinetics of N2 are complex because the molecules with cis or trans isomers at Pro161 both fold to native-like conformations, albeit with different rates. We employed consensus design to shorten the seven-residue irregular loop around Pro161 to a four-residue type I′ turn without a proline. This increased the conformational stability of N2 by almost 10 kJ mol− 1 and abolished the complexity of the folding kinetics. Turn sequences obtained from in vitro selections for increased stability strongly resembled those derived from the consensus design. Two other type I′ turns of N2 could also be stabilized by consensus design. For all three turns, the gain in stability originates from an increase in the rate of refolding. The turns form native-like structures early during refolding and thus stabilize the folding transition state. The crystal structure of the variant with all three stabilized turns confirms that the 157–163 loop was in fact shortened to a type I′ turn and that the other turns maintained their type I′ conformation after sequence optimization.

Prolyl cis/trans isomerizations determine the rates of many protein-folding reactions, and they can serve as molecular switches and timers. The energy required to shift the prolyl cis/trans equilibrium during these processes originates from conformational reactions that are linked structurally and energetically with prolyl isomerization. We used the N2 domain of the gene-3-protein of phage fd to elucidate how such an energetic linkage develops in the course of folding. The Asp160-Pro161 bond at the tip of a beta hairpin of N2 is cis in the crystal structure, but in fact, it exists as a mixture of conformers in folded N2. During refolding, about 10 kJ mol(-1) of conformational energy becomes available for a 75-fold shift of the cis/trans equilibrium constant at Pro161, from 7/93 in the unfolded to 90/10 in the folded form. We combined single- and double-mixing kinetic experiments with a mutational analysis to identify the structural origin of this proline shift energy and to elucidate the molecular path for the transfer of this energy to Pro161. It originates largely, if not entirely, from the two-stranded beta sheet at the base of the Pro161 hairpin. The two strands improve their stabilizing interactions when Pro161 is cis, and this stabilization is propagated to Pro161, because the connector peptides between the beta strands and Pro161 are native-like folded when Pro161 is cis. In the presence of a trans-Pro161, the connector peptides are locally unfolded, and thus, Pro161 is structurally and energetically uncoupled from the beta sheet. Such interrelations between local folding and prolyl isomerization and the potential modulation by prolyl isomerases might also be used to break and reestablish slow communication pathways in proteins.

Summary SlyD (sensitive to lysis D) is a putative folding helper from the bacterial cytosol and harbors prolyl isomerase and chaperone activities. We determined the solution NMR structure of a truncated version of SlyD (1–165) from Escherichia coli (SlyD⁎) that lacks the presumably unstructured C-terminal tail. SlyD⁎ consists of two well-separated domains: the FKBP domain, which harbors the prolyl isomerase activity, and the insert-in-flap (IF) domain, which harbors the chaperone activity. The IF domain is inserted into a loop of the FKBP domain near the prolyl isomerase active site. The NMR structure of SlyD⁎ showed no distinct orientation of the two domains relative to each other. In the FKBP domain, Tyr68 points into the active site, which might explain the lowered intrinsic prolyl isomerase activity and the much lower FK506 binding affinity of the protein compared with archetype human FKBP12 (human FK506 binding protein with 12 kDa). The thermodynamics and kinetics of substrate binding by SlyD⁎ were quantified by fluorescence resonance energy transfer. NMR titration experiments revealed that the IF domain recognizes and binds unfolded or partially folded proteins and peptides. Insulin aggregation is markedly slowed by SlyD⁎ as evidenced by two-dimensional NMR spectroscopy in real time, probably due to SlyD⁎ binding to denatured insulin. The capacity of the IF domain to establish an initial encounter–collision complex, together with the flexible orientation of the two interacting domains, makes SlyD⁎ a very powerful catalyst of protein folding.

SlyD, the sensitive-to-lysis protein from Escherichia coli, consists of two domains. They are not arranged successively along the protein chain, but one domain, the "insert-in-flap" (IF) domain, is inserted internally as a guest into a surface loop of the host domain, which is a prolyl isomerase of the FK506 binding protein (FKBP) type. We used SlyD as a model to elucidate how such a domain insertion affects the stability and folding mechanism of the host and the guest domain. For these studies, the two-domain protein was compared with a single-domain variant SlyDDeltaIF, SlyD* without the chaperone domain (residues 1-69 and 130-165) in which the IF domain was removed and replaced by a short loop, as present in human FKBP12. Equilibrium unfolding and folding kinetics followed an apparent two-state mechanism in the absence and in the presence of the IF domain. The inserted domain decreased, however, the stability of the host domain in the transition region and decelerated its refolding reaction by about 10-fold. This originates from the interruption of the chain connectivity by the IF domain and its inherent instability. To monitor folding processes in this domain selectively, a Trp residue was introduced as fluorescent probe. Kinetic double-mixing experiments revealed that, in intact SlyD, the IF domain folds and unfolds about 1000-fold more rapidly than the FKBP domain, and that it is strongly stabilized when linked with the folded FKBP domain. The unfolding limbs of the kinetic chevrons of SlyD show a strong downward curvature. This deviation from linearity is not caused by a transition-state movement, as often assumed, but by the accumulation of a silent unfolding intermediate at high denaturant concentrations. In this kinetic intermediate, the FKBP domain is still folded, whereas the IF domain is already unfolded.

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.

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.

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.

Unfolding usually leads to the loss of the biological function of a protein. Here, we show that an unfolding reaction activates the gene-3-protein of the filamentous phage fd for its function during the infection of Escherichia coli. Before infection, the gene-3-protein is in a fully folded locked form, in which the binding site for the phage receptor TolA is buried at the domain interface. To expose this binding site, the gene-3-protein must be activated, and previously we identified the cis-to-trans isomerization at Pro213 in the hinge region between the two domains as a key step of activation. We now report that Pro213 isomerization destabilizes the protein and leads to a loss of folded structure, presumably in the hinge region. The partially unfolded form of the gene-3-protein is metastable, and trans-Pro213 arrests the protein in this activated form for an extended time, long enough to find the receptor TolA. The partial unfolding and its timing by prolyl isomerization are essential for the biological function.

The p19INK4d protein consists of five ankyrin repeats (ANK) and controls the human cell cycle by inhibiting the cyclin D-dependent kinases (CDK) 4 and 6. We investigated the folding of p19INK4d by urea-induced unfolding transitions, kinetic analyses of unfolding and refolding, including double-mixing experiments and a special assay for folding intermediates. Folding is a sequential two-step reaction via a hyperfluorescent on-pathway intermediate. This intermediate is present under all conditions, during unfolding, refolding and at equilibrium. The folding mechanism was confirmed by a quantitative global fit of a consistent set of equilibrium and kinetic data revealing the thermodynamics and intrinsic folding rates of the different states. Surprisingly, the N ↔ I transition is much faster compared to the I ↔ U transition. The urea-dependence of the intrinsic folding rates causes population of the intermediate at equilibrium close to the transition midpoint. NMR detected hydrogen/deuterium exchange and the analysis of truncated variants showed that the C-terminal repeats ANK3-5 are already folded in the on-pathway intermediate, whereas the N-terminal repeats 1 and 2 are not folded. We suggest that during refolding, repeats ANK3–ANK5 first form the scaffold for the subsequent assembly of repeats ANK1 and ANK2. The binding function of p19INK4d resides in the latter repeats. We propose that the graded stability and the facile unfolding of repeats 1 and 2 is a prerequisite for the down-regulation of the inhibitory activity of p19INK4d during the cell-cycle.

The catalytic activity of human FKBP12 as a prolyl isomerase is high towards short peptides, but very low in proline-limited protein folding reactions. In contrast, the SlyD proteins, which are members of the FKBP family, are highly active as folding enzymes. They contain an extra "insert-in-flap" or IF domain near the prolyl isomerase active site. The excision of this domain did not affect the prolyl isomerase activity of SlyD from Escherichia coli towards short peptide substrates but abolished its catalytic activity in proline-limited protein folding reactions. The reciprocal insertion of the IF domain of SlyD into human FKBP12 increased its folding activity 200-fold and generated a folding catalyst that is more active than SlyD itself. The IF domain binds to refolding protein chains and thus functions as a chaperone module. A prolyl isomerase catalytic site and a separate chaperone site with an adapted affinity for refolding protein chains are the key elements for a productive coupling between the catalysis of prolyl isomerization and conformational folding in the enzymatic mechanisms of SlyD and other prolyl isomerases, such as trigger factor and FkpA.

The signal adapter protein c-CrkII from chicken but not from human uses isomerization at Pro238 in the SH3C domain to regulate the activity of the SH3N domain. The different behavior of human and chicken c-CrkII originates from only two differences in sequence, at positions 239 after Pro238 and 272 in the N-Src loop of SH3C. We analyzed the kinetics of substrate binding to SH3N and an assay for its coupling with Pro238 isomerization in SH3C to identify the molecular path from Pro238 to the substrate binding site of SH3N. The trans→cis isomerization at Pro238 and a relocation of Phe239 re-organize the energetics of a hydrophobic cluster in the N-Src loop of SH3C and re-shape this region to optimize its interactions with SH3N. Concomitantly, the backbone becomes strained at Met272. We suggest that, in human c-CrkII, movement at position 239 and strain at position 272 are not tolerated because the β-branched residues Ile239 and Val272 restrain the backbone mobility and thus destabilize the cis Pro238 form.

The cold shock protein CspB shows a five-stranded beta-sheet structure, and it folds rapidly via a native-like transition state. A previous Phi value analysis showed that most of the residues with Phi values close to one reside in strand beta1, and two of them, Lys5 and Lys7 are partially exposed charged residues. To elucidate how coulombic interactions of these two residues contribute to the energetic organisation of the folding transition state we performed comparative folding experiments in the presence of an ionic denaturant (guanidinium chloride) and a non-ionic denaturant (urea) and a double-mutant analysis. Lys5 contributes 6.6 kJ mol(-1) to the stability of the transition state, and half of it originates from screenable coulombic interactions. Lys7 contributes 5.3 kJ mol(-1), and 3.4 kJ mol(-1) of it are screened by salt. In the folded protein Lys7 interacts with Asp25, and the screenable coulombic interaction between these two residues is fully formed in the transition state. This suggests that long-range coulombic interactions such as those originating from Lys5 and Lys7 of CspB can be important for organizing and stabilizing native-like structure early in protein folding.

The two-state folding reaction of the cold shock protein from Bacillus caldolyticus (Bc-Csp) is preceded by a rapid chain collapse. A fast shortening of intra-protein distances was revealed by Förster resonance energy transfer (FRET) measurements with protein variants that carried individual pairs of donor and acceptor chromophores at various positions along the polypeptide chain. Here we investigated the specificity of this rapid compaction. Energy transfer experiments that probed the stretching of strand beta2 and the close approach between the strands beta1 and beta2 revealed that the beta1-beta2 hairpin is barely formed in the collapsed form, although it is native-like in the folding transition state of Bc-Csp. The time course of the collapse could not be resolved by pressure or temperature jump experiments, indicating that the collapsed and extended forms are not separated by an energy barrier. The co-solute (NH4)2SO4 stabilizes both native Bc-Csp and the collapsed form, which suggests that the large hydrated SO4(2-) ions are excluded from the surface of the collapsed form in a similar fashion as they are excluded from folded Bc-Csp. Ethylene glycol increases the stability of proteins because it is excluded preferentially from the backbone, which is accessible in the unfolded state. The collapsed form of Bc-Csp resembles the unfolded form in its interaction with ethylene glycol, suggesting that in the collapsed form the backbone is still accessible to water and small molecules. Our results thus rule out that the collapsed form is a folding intermediate with native-like chain topology. It is better described as a mixture of compact conformations that belong to the unfolded state ensemble. However, some of its structural elements are reminiscent of the native protein.

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.

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.

The gene-3 protein (G3P) of filamentous phages is essential for the infection of Escherichia coli. The carboxy-terminal domain anchors this protein in the phage coat, whereas the two amino-terminal domains N1 and N2 protrude from the phage surface. We analyzed the folding mechanism of the two-domain fragment N1-N2 of G3P (G3P(*)) and the interplay between folding and domain assembly. For this analysis, a variant of G3P(*) was used that contained four stabilizing mutations (IIHY-G3P(*)). The observed refolding kinetics extend from 10 ms to several hours. Domain N1 refolds very rapidly (with a time constant of 9.4 ms at 0.5 M guanidinium chloride, 25 degrees C) both as a part of IIHY-G3P(*) and as an isolated protein fragment. The refolding of domain N2 is slower and involves two reactions with time constants of seven seconds and 42 seconds. These folding reactions of the individual domains are followed by a very slow, spectroscopically silent docking process, which shows a time constant of 6200 seconds. This reaction was detected by a kinetic unfolding assay for native molecules. Before docking, N1 and N2 unfold fast and independently, after docking they unfold slowly in a correlated fashion. A high energy barrier is thus created by domain docking, which protects G3P kinetically against unfolding. The slow domain docking is possibly important for the infection of E.coli by the phage. Upon binding to the F pilus, the N2 domain separates from N1 and the binding site for TolA on domain N1 is exposed. Since domain reassembly is so slow, this binding site remains accessible until pilus retraction has brought N1 close to TolA on the bacterial surface.

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.

Prolyl isomerizations are intrinsically slow processes. They determine the rates of many protein folding reactions and control regulatory events in folded proteins. Prolyl isomerases are able to catalyze these isomerizations, and thus, they have the potential to assist protein folding and to modulate protein function. Here, we provide examples for how prolyl isomerizations limit protein folding and are accelerated by prolyl isomerases and how native-state prolyl isomerizations regulate protein functions. The roles of prolines in protein folding and protein function are closely interrelated because both of them depend on the coupling between cis/trans isomerization and conformational changes that can involve extended regions of a protein.

The cold-shock protein CspB folds rapidly in a N ⇄ U two-state reaction via a transition state that is about 90% native in its interactions with denaturants and water. This suggested that the energy barrier to unfolding is overcome by processes occurring in the protein itself, rather than in the solvent. Nevertheless, CspB unfolding depends on the solvent viscosity. We determined the activation volumes of unfolding and refolding by pressure-jump and high-pressure stopped-flow techniques in the presence of various denaturants. The results obtained by these methods agree well. The activation volume of unfolding is positive ( Δ V NU ‡ =16(±4) ml / mol ) and virtually independent of the nature and the concentration of the denaturant. We suggest that in the transition state the protein is expanded and water molecules start to invade the hydrophobic core. They have, however, not yet established favorable interactions to compensate for the loss of intra-protein interactions. The activation volume of refolding is positive as well ( Δ V UN ‡ =53(±6) ml / mol ) and, above 3 M urea, independent of the concentration of the denaturant. At low concentrations of urea or guanidinium thiocyanate, Δ V UN ‡ decreases significantly, suggesting that compact unfolded forms become populated under these conditions.

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 (ΔGD) 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.

Mia40 (a mitochondrial import and assembly protein) catalyzes disulfide bond formation in proteins in the mitochondrial intermembrane space. By using Cox17 (a mitochondrial copper-binding protein) as a natural substrate, we discovered that, in the presence of Mia40, the formation of native disulfides is strongly favored. The catalytic mechanism of Mia40 involves a functional interplay between the chaperone site and the catalytic disulfide. Mia40 forms a specific native disulfide in Cox17 much more rapidly than other disulfides, in particular, non-native ones, which originates from the recently described high affinity for hydrophobic regions near target cysteines and the long lifetime of the mixed disulfide. In addition to its thiol oxidase function, Mia40 is active also as a disulfide reductase and isomerase. We found that species with inadvertently formed incorrect disulfides are rebound by Mia40 and reshuffled, revealing a proofreading mechanism that is steered by the conformational folding of the substrate protein.

The signaling protein CrkII switches between forms with high or low binding affinity. Both phosphorylation and native-state prolyl isomerization were suggested to regulate the transition between these forms. Here we analyzed how phosphorylation at Tyr222 and Tyr252 and the Pro238Ala substitution affect signal transfer of human and chicken CrkII to a downstream target. Human CrkII is regulated by phosphorylation only, but chicken CrkII is regulated by Pro238 trans→cis isomerization and by Tyr222 phosphorylation. Surprisingly, they act in an independent fashion. Apparently, the allosteric transition to a low-activity form can be induced by phosphorylation or prolyl isomerization located at distant sites in CrkII.

The cold shock proteins Bc-Csp from the thermophile Bacillus caldolyticus and Bs-CspB from the mesophile Bacillus subtilis differ significantly in their conformational stability, although the two proteins differ by only 12 out of 67 amino acid residues. The three-dimensional structure of these small and compact β-barrel proteins without disulfide bonds, cis-proline residues or tightly bound cofactors is very similar. Previous work has shown that Bc-Csp displays a twofold increase in the free energy of stabilization relative to its homolog Bs-CspB, and indicated that electrostatic interactions are, in part, responsible for this effect. It was further described that the stability difference is almost exclusively due to surface-exposed charged residues at sequence positions 3 and 66 of Bc-Csp and Bs-CspB, whereas all other amino acid changes between both proteins have no net effect on stability. To investigate how two surface residues determine the stability of Bc-Csp, Arg3 and Leu66 were replaced by glutamic acid, corresponding to the Bs-CspB sequence. The crystal structures of the resultant protein variants, Bc-Csp R3E and Bc-Csp L66E, were determined at 1.4 A and 1.27 A resolution, and refined to R values of 13.9 % and 15.8 %, respectively. Both structures closely resemble Bc-Csp in their global fold and show different hydrogen bonding and salt-bridge patterns when two independent molecules in the asymmetric unit of the crystal are compared. To extend the study to neighbored residues that help determine the surface charge around Arg3 and Leu66, the mutant proteins Bc-Csp E46A, Bc-Csp R3E/E46A/L66E and Bc-Csp V64T/L66E/67A were crystallized. Their structures were determined at resolutions of 1.8 A, 1.32 A and 1.8 A and refined to R values of 18.5 %, 13.8 % and 19.3 %, respectively. A systematic comparison of the crystal structures of all forms of the B. caldolyticus cold shock protein shows varying patterns of hydrogen bonds and electrostatic interactions around residues 3 and 66. Thermal destabilization of the protein by mutation appears to correlate with the extent of an acidic surface patch near the C-terminal carboxylate group.

The bacterial cold shock proteins are small compact beta-barrel proteins without disulfide bonds, cis-proline residues or tightly bound cofactors. Bc-Csp, the cold shock protein from the thermophile Bacillus caldolyticus shows a twofold increase in the free energy of stabilization relative to its homolog Bs-CspB from the mesophile Bacillus subtilis, although the two proteins differ by only 12 out of 67 amino acid residues. This pair of cold shock proteins thus represents a good system to study the atomic determinants of protein thermostability. Bs-CspB and Bc-Csp both unfold reversibly in cooperative transitions with T(M) values of 49.0 degrees C and 77.3 degrees C, respectively, at pH 7.0. Addition of 0.5 M salt stabilizes Bs-CspB but destabilizes Bc-Csp. To understand these differences at the structural level, the crystal structure of Bc-Csp was determined at 1.17 A resolution and refined to R=12.5% (R(free)=17.9%). The molecular structures of Bc-Csp and Bs-CspB are virtually identical in the central beta-sheet and in the binding region for nucleic acids. Significant differences are found in the distribution of surface charges including a sodium ion binding site present in Bc-Csp, which was not observed in the crystal structure of the Bs-CspB. Electrostatic interactions are overall favorable for Bc-Csp, but unfavorable for Bs-CspB. They provide the major source for the increased thermostability of Bc-Csp. This can be explained based on the atomic-resolution crystal structure of Bc-Csp. It identifies a number of potentially stabilizing ionic interactions including a cation-binding site and reveals significant changes in the electrostatic surface potential.

The rate-limiting step during the refolding of S54G/P55N ribonuclease T1 is determined by the slow trans-->cis prolyl isomerisation of Pro39. We investigated the refolding of this variant by one-dimensional (1D) and two-dimensional (2D) real-time NMR spectroscopy, initiated by a tenfold dilution from 6 M guanidine hydrochloride at 10 degreesC. Two intermediates could be resolved with the 1D approach. The minor intermediate, which is only present early during refolding, is largely unfolded. The major intermediate, with an incorrect trans Pro39 peptide bond, is highly structured with 33 amide protons showing native chemical shifts and native NOE patterns. They could be assigned in a real-time 2D-NOESY (nuclear Overhauser enhancement spectroscopy) by using a new assignment strategy to generate positive and negative signal intensities for native and non-native NOE cross-peaks, respectively. Surprisingly, amide protons with non-native environments are located not only close to Tyr38-Pro39, but are spread throughout the entire protein, including the C-terminal part of the alpha-helix, beta-strands 3 and 4 and several loop regions. Native secondary and tertiary structure was found for the major intermediate in the N-terminal beta-strands 1 and 2 and the C terminus (connected by the disulfide bonds), the N-terminal part of the alpha-helix, and the loops between beta-strands 4/5 and 5/6. Implications of these native and non-native structure elements of the intermediate for the refolding of S54G/P55N ribonuclease T1 and for cis/trans isomerizations are discussed.

The trigger factor is associated with bacterial ribosomes and catalyzes proline-limited protein folding reactions. Its folding activity is very high and conserved in evolution, as shown for the homologous enzymes from Escherichia coli and Mycoplasma genitalium. The folding protein substrate (a variant of ribonuclease T1) binds with high affinity to the trigger factors, and permanently unfolded proteins are strong, competitive inhibitors. We used this inhibition to characterize the substrate binding sites of the trigger factors. Unfolded α-lactalbumin binds very tightly and inhibits the trigger factor from M. genitalium with a KI value of 50 nM. The binding of inhibitory proteins is independent of proline residues, as shown for unfolded tendamistat, which binds to the trigger factor with equal affinity in the presence and in the absence of its three proline residues. The good inhibition by a non-folding variant of ribonuclease T1 that lacks Pro39 showed that this proline, at which the catalysis of folding occurs, is dispensable for substrate binding. The trigger factors cannot catalyze prolyl isomerization when proteins are partially folded already. They preferentially recognize unstructured protein chains, which bind with high affinity to a site distinct from the catalytic prolyl isomerase center in the FKBP domain.

Recently a new family of prolyl isomerases was discovered, which is unrelated with the cyclophilins or the FK-506 binding proteins. Parvulin, the smallest member of this new family, is a protein with only 92 residues, but parvulin-like domains occur in several large proteins that are apparently involved in protein folding or activation processes. We show here that, in addition to its activity in assays with proline-containing tetrapeptides, parvulin catalyzes the proline-limited folding of a variant of ribonuclease T1 with a kcat/Km value of 30,000 M-1 s-1. This value is much smaller than the kcat/Km value of 1.1x10(7) M-1 s-1 determined for the parvulin-catalyzed prolyl isomerization in the tetrapeptide succinyl-Ala-Leu-Pro-Phe-4-nitroanilide. Parvulin itself unfolds and refolds reversibly in a simple two-state reaction with a Gibbs free energy of stabilization of 28 kJ/mol at 10 degrees C. Most of the unfolded parvulin molecules refold in a slow reaction that involves prolyl isomerization and is catalyzed by cyclophilin, another prolyl isomerase. Moreover, parvulin accelerates its own refolding in an autocatalytic fashion, and the rate of refolding increases tenfold when the concentration of parvulin is increased from 0.5 to 3.0 microM. Apparently, small single-domain prolyl isomerases catalyze prolyl isomerization much better in short peptides than in protein folding reactions, presumably because the prolyl bonds are less accessible in refolding protein chains. It is possible that the additional domains of the large prolyl isomerases improve the affinity for protein substrates.

Electrostatic interactions between charged residues and the helix dipole in a protein were investigated by protein engineering methods. In ribonuclease T1, two surface-exposed acidic residues (Glu28 and Asp29) are located near the carboxyl terminus of the alpha-helix between residues 13 and 29. They were replaced, individually and in concert, by the uncharged amides Gln28 and Asn29, and the stabilities of the wild-type protein and its variants were determined as a function of pH. The effects of the two mutations are additive. Either one leads to a marginal destabilization by 0.7 kJ/mol at pH 2 but to a strong stabilization by about 3.2 kJ/mol at pH 7. This suggests that the deprotonations of Glu28 and Asp29 reduce the free energy of stabilization of folded ribonuclease T1 by about 4 kJ/mol each. This destabilization is probably caused by unfavorable electrostatic interactions of Glu28 and Asp29 with the negative end of the helix dipole. The activation energies for the unfolding of the different variants of ribonuclease T1 change in parallel with the differences in the thermodynamic stability when the pH is varied. This indicates that the unfavorable electrostatic interactions of Glu28 and Asp29 are lost very early in unfolding, and are not present in the activated state of unfolding.

In wild-type ribonuclease T1 the peptide bond between Tyr38 and Pro39 is in the cis conformation. When Pro39 is replaced by an alanine this cis conformation is retained, and a non-prolyl cis Tyr38-Ala39 peptide bond is generated. We employed a stopped-flow double-mixing technique to investigate the kinetics of the cis-->trans isomerization of this peptide bond in the unfolding and the trans-->cis isomerization in the refolding of Pro39Ala-ribonuclease T1. In 6.0 M GdmCl (pH 1.6) and 25 degrees C the protein unfolds rapidly with a time constant of 20 ms, followed by Tyr38-Ala39 cis-->trans isomerization. This reaction shows a time constant of 730 ms and is about 60-fold faster than the isomerization of the Tyr38-Pro39 bond in the wild-type protein. Unfolded molecules with the Tyr38-Ala39 bond still in the native-like cis conformation accumulate transiently for a short time after unfolding is initiated, and they can refold very rapidly to the native state with a time constant of 290 ms (in 1.0 M GdmCl, pH 4.6, 25 degrees C). After more than three seconds of unfolding virtually all protein molecules contain an incorrect trans Tyr38-Ala39 bond and refolding is decelerated approximately 1000-fold, because Tyr38-Ala39 trans-->cis re-isomerization is very slow and, with its time constant of 480 s, determines the overall rate of refolding. Due to the coupling of the cis-trans equilibrium with protein folding it was possible to measure the kinetic parameters of the isomerization of a non-prolyl peptide bond in a protein. Previously this could not be accomplished, because the trans isomer is strongly preferred for unsubstituted peptide bonds in oligopeptides under virtually all conditions. Our data indicate that the kinetics of Tyr38-Pro39 and of Tyr38-Ala39 isomerization differ predominantly in the rate of the cis-->trans, rather than of the trans-->cis reaction. The rate of the trans-->cis reaction is, however, measured during refolding and may be influenced by the formation of ordered protein structure.

Disulfide bonds in a folding protein chain are equivalent to prematurely formed native-like tertiary interactions. We investigated whether the mechanism of protein folding is changed by the presence of disulfide bonds. As a model we used the S54G/P55N-variant of ribonuclease T1, a protein with two disulfide bonds and a single cis proline (Pro39), and we measured both the direct and the proline-limited folding reactions before and after breaking of the disulfide bonds. The folding kinetics were compared under refolding conditions, in the regions of the area-induced unfolding transitions of the two forms, and under folding conditions. The kinetics in the transition regions were analyzed on the basis of a three-species mechanism and all microscopic rate constants of folding and of prolyl isomerization could be determined as a function of the area concentration from the measured rates and amplitudes. These kinetic analyses indicated that the disulfide bonds can be rather unfavorable for the folding of S54G/P55N-ribonuclease T1. Under strongly native conditions they retard the rate-limiting trans → cis isomerization of Pro39 because they allow the rapid formation of partially ordered structure prior to the proline-limited refolding reaction. Under unfolding conditions the isomerization of Pro39 is not affected. The direct unfolding and refolding reactions in the transition region of polypeptide chains with correct prolyl isomers are also decelerated when the disulfide bonds are present. These changes in the folding kinetics are possibly related to the decrease in chain flexibility that is caused by the disulfide bonds. A high flexibility is probably important throughout folding, and in the case of ribonuclease T1 a premature locking of tertiary contacts by intact disulfide bonds can interfere unfavorably with both the direct and the proline-limited folding reactions.

The refolding of ribonuclease T1 involves two major slow processes that exhibit properties of prolyl isomerization reactions. A comparison of the wild-type protein and a designed variant where the cis Ser54-Pro55 bond was replaced by a Gly54-Asn55 bond indicated that the faster of these reactions is the isomerization of Pro55. Here we report the replacement of the other cis proline of ribonuclease T1 at position 39 by alanine. The Pro39Ala variant is similar to the wild-type protein in secondary and tertiary structure, and the enzymatic activity towards RNA and a dinucleotide substrate remains almost unchanged. The fluorescence emission of the single Trp59 is lowered by the Pro39Ala substitution, probably because Trp59 is in close contact to Pro39 in wild-type ribonuclease T1. Unlike the substitution of cis Pro55, the Pro39Ala mutation is strongly destabilizing and reduces the Gibbs free energy of the folded protein by about 20 kJ/mol. Pro39 is buried in native RNase T1 and located near the active site. The observed destabilization could originate from the presence of a cis alanyl bond in the Pro39Ala variant or from a local distortion caused by the incorporation of a trans alanyl peptide bond in the interior of the protein. In the refolding kinetics the replacement of Pro39 leads to a disappearance of the fast-refolding species. Refolding still involves two consecutive slow steps. The first and faster step could be the isomerization of the remaining cis Pro55. The second, very slow step is a novel reaction that appears to have no counterpart in the refolding of the wild-type protein. All mutant molecules must undergo this reaction before reaching the native state. These major changes in the folding kinetics strongly indicate that cis-Pro39 is indeed of major importance for the folding of the wild-type protein. They indicate, moreover, that some new feature of protein folding kinetics is observed in these studies of the Pro39Ala variant.

The SlyD (sensitive to lysis D) protein of Escherichia coli is a folding enzyme with a chaperone domain and a prolyl isomerase domain of the FK506 binding protein type. Here we investigated how the two domains and their interplay are optimized for function in protein folding. Unfolded protein molecules initially form a highly dynamic complex with the chaperone domain of SlyD, and they are then transferred to the prolyl isomerase domain. The turnover number of the prolyl isomerase site is very high and guarantees that, after transfer, prolyl peptide bonds in substrate proteins are isomerized very rapidly. The Michaelis constant of catalyzed folding reflects the substrate affinity of the chaperone domain, and the turnover number is presumably determined by the rate of productive substrate transfer from the chaperone to the prolyl isomerase site and by the intrinsic propensity of the refolding protein chain to leave the active site with the native prolyl isomer. The efficiency of substrate transfer is high because dissociation from the chaperone site is very fast and because the two sites are close to each other. Protein molecules that left the prolyl isomerase site with an incorrect prolyl isomer can rapidly be re-bound by the chaperone domain because the association rate is very high as well.

The folding and unfolding kinetics within the transition region were measured for RNase A and for RNase T1. The data were used to evaluate the theoretical models for the influence of prolyl isomerization on the observed folding kinetics. These two proteins were selected, since the folding reaction of RNase A is faster than prolyl isomerization, whereas in RNase T1, folding is slower than isomerization in the transition region. Folding of RNase T1 was investigated for three variants with different numbers of cis prolyl residues. The results indicate that in the transition region the folding rates are indeed strongly dependent on the number of prolyl residues. The variant of RNase T1 that contains only one cis prolyl residue folds about ten times faster than two variants that contain two cis prolyl residues. For both RNase A and RNase T1, the apparent rates of folding and unfolding as well as the corresponding amplitudes depend on the concentration of denaturant in a manner that was predicted by the model calculations. When refolding was started from the fast-folding species, additional kinetic phases could be observed in the transition region for both proteins. The obtained values could be used to calculate the microscopic rate constants of folding and isomerization on the basis of theoretical models.

Kinetic models were developed to describe the influence of prolyl peptide bond isomerization on the kinetics of reversible protein folding for cases in which structural intermediates do not occur. In the simulations, the number of prolyl residues and the relative rates of folding and isomerization were varied. The experimentally observed rate constants were found to be identical with the intrinsic rate constants of folding and isomerization only when folding remains much faster than prolyl isomerization throughout the transition region. When the rate of folding becomes similar to or lower than the rate of isomerization, the observed kinetic parameters are complex functions of all microscopic rate constants. In particular, the observed folding rates in the transition region decrease with the number of prolyl residues. Pseudo two-state kinetics with single folding and unfolding reactions are observed in several cases, although the apparent folding rates depend strongly on prolyl isomerization reactions in the unfolded chain. This virtual simplicity can easily lead to misinterpretation of kinetic data. Additional phases can be resolved when refolding is started from the fast-folding species (UF). The coupling between folding and prolyl peptide bond isomerization also modifies the dependence on denaturant concentration of the apparent rate constants of folding. We suggest several tests to detect and characterize the contributions of folding and isomerization steps to the observed folding kinetics.

In previous work, a strongly stabilized variant of the beta1 domain of streptococcal protein G (Gbeta1) was obtained by an in vitro selection method. This variant, termed Gbeta1-M2, contains the four substitutions E15V, T16L, T18I, and N37L. Here we elucidated the molecular basis of the observed strong stabilizations. The contributions of these four residues were analyzed individually and in various combinations, additional selections with focused Gbeta1 gene libraries were performed, and the crystal structure of Gbeta1-M2 was determined. All single substitutions (E15V, T16L, T18I, and N37L) stabilize wild-type Gbeta1 by contributions of between 1.6 and 6.0 kJ mol(-1) (at 70 degrees C). Hydrophobic residues at positions 16 and 37 provide the major contribution to stabilization by enlarging the hydrophobic core of Gbeta1. They also increase the tendency to form dimers, as shown by dependence on the concentration of apparent molecular mass in analytical ultracentrifugation, by concentration-dependent stability, and by a strongly increased van't Hoff enthalpy of unfolding. The 0.88-A crystal structure of Gbeta1-M2 and NMR measurements in solution provide the explanation for the observed dimer formation. It involves a head-to-head arrangement of two Gbeta1-M2 molecules via six intermolecular hydrogen bonds between the two beta strands 2 and 2' and an adjacent self-complementary hydrophobic surface area, which is created by the T16L and N37L substitutions and a large 120 degrees rotation of the Tyr33 side chain. This removal of hydrophilic groups and the malleability of the created hydrophobic surface provide the basis for the dimer formation of stabilized Gbeta1 variants.

Unfolded ribonuclease (RNase) from porcine pancreas consists of a mixture of fast and slow-refolding species. The equilibrium distribution of these species differs strongly from other homologous RNases, because an additional proline residue is present at position 115 of the porcine protein. The major slow-folding species of porcine RNase contains incorrect proline isomers at Pro93 and at Pro114-Pro115. Both positions are presumably part of β-turn structures in the native protein, as deduced from the structure of the homologous bovine RNase A. The folding kinetics of these molecules depend strongly on the conditions used. Under unfavorable conditions (near the unfolding transition), refolding is virtually blocked by the presence of the incorrect proline peptide bonds and partially folded intermediates with incorrect isomers could not be detected. As a consequence, folding is very slow under such conditions and the re-isomerization of Pro114-Pro115 is the first and ratelimiting step of folding. Under strongly native conditions (such as in the presence of ammonium sulfate), refolding is much faster. A largely folded intermediate accumulates with the turns around Pro93 and Pro114-Pro115 still in the non-native conformation. These results suggest that incorrect proline isomers strongly influence protein folding and that, under favorable conditions, the polypeptide chain can fold with two β-turns locked into a non-native conformation. We conclude, therefore, that early formation of correct turn structure is not necessarily required for protein folding. However, the presence of incorrect turns, locked-in by non-native proline isomers, strongly decreases the rate of refolding. Alternative pathways of folding exist. The choice of pathway depends on the number and distribution of incorrect proline isomers and on the folding conditions.

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 (Gβ1), 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 β-sheet of Gβ1. 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.

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 °C to 83.7 °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 β-barrel structure but not well optimized in the wild-type protein.

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.

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.

The envelope proteins of human immunodeficiency virus (HIV) and human T-cell lymphotrophic virus (HTLV) mediate cell attachment and membrane fusion. For HIV-1, the precursor protein gp160 is cleaved proteolytically into two fragments, the surface-associated receptor binding subunit gp120 and the membrane spanning subunit gp41, which is involved in membrane fusion during virus entry. Soluble and immunoreactive variants of gp41 are essential for the reliable diagnosis of HIV-1 infections. Hitherto, gp41 was solubilized by adding detergents, or in acidic or alkaline solvents. We find that covalent fusions with SlyD or FkpA, two homodimeric Escherichia coli chaperones with peptidyl-prolyl isomerase activity, solubilize retroviral envelope proteins without compromising their immunological reactivity. gp41 from HIV-1, gp36 from HIV-2 and gp21 from HTLV could be expressed in large amounts in the Escherichia coli cytosol when fused with one or two subunits of SlyD or FkpA. The fusion proteins could be easily isolated and refolded, and showed high solubility and immunoreactivity, thus providing sensitive and reliable tools for diagnostic applications. Covalent fusions with SlyD or FkpA might be valuable generic tools for the solubilization and activation of aggregation-prone proteins.

The cold shock protein Bc-Csp folds very rapidly in a reaction that is well described by a kinetic two-state mechanism without intermediates. We measured the shortening of six intra-protein distances during folding by Forster resonance energy transfer (FRET) in combination with stopped-flow experiments. Single tryptophan residues were engineered into the protein as the donors, and single 5-(((acetylamino)ethyl)amino)naphthalene-1-sulfonate (AEDANS) residues were placed as the acceptors at solvent-exposed sites of Bc-Csp. Their R0 value of about 22 A was well suited for following distance changes during the folding of this protein with a high sensitivity. The mutagenesis and the labeling did not alter the refolding kinetics. The changes in energy transfer during folding were monitored by both donor and acceptor emission and reciprocal effects were found. In two cases the donor-acceptor distances were similar in the unfolded and the folded state and, as a consequence, the kinetic changes in energy transfer upon folding were very small. For four donor/acceptor pairs we found that > or =50% of the increase in energy transfer upon folding occurred prior to the rate-limiting step of folding. This reveals that about half of the shortening of the intra-molecular distances upon folding has occurred already before the rate-limiting step and suggests that the fast two-state folding reaction of Bc-Csp is preceded by a very rapid collapse.

The amino-terminal domains N1 and N2 of the gene-3-protein of phage fd form a bilobal structural and functional entity that protrudes from the phage tip. Domain N2 initiates the infection of Escherichia coli by binding to the F pilus. This binding results in the dissociation of the two domains and allows N1 to interact with the TolA receptor at the cell surface. The refolding of the N1-N2 fragment begins with the folding of domain N1, which takes a few milliseconds, followed by the folding of domain N2, which is complete within five minutes. The subsequent domain assembly is unusually slow and shows a time-constant of 6200 s at 25 degrees C. We found that the rate of this reaction is controlled by the trans to cis isomerization of the Gln212-Pro213 bond in the hinge subdomain of N2, a region that provides many interactions between N1 and N2 in the gene-3-protein. The substitution of Pro213 by Gly accelerated domain association 30-fold and revealed that the folding of the two individual domains and their assembly are indeed sequential steps in the refolding of the gene-3-protein. In the course of infection, the domains must separate to expose the binding site for TolA on domain N1. The kinetic block of domain reassembly caused by Pro213 isomerization could ensure that after the initial binding of N2 to the F pilus the open state persists until N1 and TolA are close enough for their mutual interaction. Pro213 isomerization might thus serve as a slow conformational switch in the function of the gene-3-protein.

Continuum electrostatic models are used to examine in detail the mechanism of protein stabilization and destabilization due to salt near physiological concentrations. Three wild-type cold shock proteins taken from mesophilic, thermophilic, and hyperthermophilic bacteria are studied using these methods. The model is validated by comparison with experimental data collected for these proteins. In addition, a number of single point mutants and three designed sequences are examined. The results from this study demonstrate that the sensitivity of protein stability toward salt is correlated with thermostability in the cold shock protein family. The calculations indicate that the mesophile is stabilized by the presence of salt while the thermophile and hyperthermophile are destabilized. A decomposition of the salt influence at a residue level permits identification of regions of the protein sequences that contribute toward the observed salt-dependent stability. This model is used to rationalize the effect of various point mutations with regard to sensitivity toward salt. Finally, it is demonstrated that designed cold shock protein variants exhibit electrostatic properties similar to the natural thermophilic and hyperthermophilic proteins.

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.

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.