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.     

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.     

A proline switch controls folding and domain interactions in the gene-3-protein of the filamentous phage fd

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.     

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.     

Energetic Communication between Functional Sites of the Gene-3-Protein during Infection by Phage fd

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.     

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.     

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

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

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.     

Unlocking of the filamentous bacteriophage virion during infection is mediated by the C domain of pIII

Protein III (pIII) of filamentous phage is required for both the beginning and the end of the phage life cycle. The infection starts by binding of the N-terminal N2 and N1 domains to the primary and secondary host receptors, F pilus and TolA protein, respectively, whereas the life cycle terminates by the C-terminal domain-mediated release of the membrane-anchored virion from the cell. It has been assumed that the role of the C-terminal domain of pIII in the infection is that of a tether for the receptor-binding domains N1N2 to the main body of the virion. In a poorly understood process that follows receptor binding, the virion disassembles as its protein(s) become integrated into the host inner membrane, resulting in the phage genome entry into the bacterial cytoplasm. To begin revealing the mechanism of this process, we showed that tethering the functional N1N2 receptor-binding domain to the virion via termination-incompetent C domain abolishes infection. This infection defect cannot be complemented by in trans supply of the functional C domain. Therefore, the C domain of pIII acts in concert with the receptor-binding domains to mediate the post receptor binding events in the infection. Based on these findings, we propose a model in which binding of the N1 domain to the periplasmic portion of TolA, the secondary receptor, triggers in cis a conformational change in the C domain, and that this change opens or unlocks the pIII end of the virion, allowing the entry phase of infection to proceed. To our knowledge, this is the first virus that uses the same protein domain both for the insertion into and release from the host membrane.     

Conformational stability and activity of ribonuclease T1 with zero, one, and two intact disulfide bonds

Ribonuclease T1 has two disulfide bonds linking cysteine residues 2-10 and 6-103. We have prepared a derivative of ribonuclease T1 in which one disulfide bond is broken and the cysteine residues carboxymethylated, (2-10)-RCM-T1, and three derivatives in which both disulfides are broken and the cysteine residues reduced, R-T1, carboxamidomethylated, RCAM-T1, or carboxymethylated, RCM-T1. The RNA hydrolyzing activity of these proteins has been measured, and urea and thermal denaturation studies have been used to determine conformational stability. The activity, melting temperature, and conformational stability of the proteins are: ribonuclease T1 (100%, 59.3 degrees C, 10.2 kcal/mol), (2-10)-RCM-T1 (86%, 53.3 degrees C, 6.8 kcal/mol), R-T1 (53%, 27.2 degrees C, 3.0 kcal/mol), RCAM-T1 (43%, 21.2 degrees C, 1.5 kcal/mol), and RCM-T1 (35%, 16.6 degrees C, 0.9 kcal/mol). Thus, the conformational stability is decreased by 3.4 kcal/mol when one disulfide bond is broken and by 7.2-9.3 kcal/mol when both disulfide bonds are broken. It is quite remarkable that RNase T1 can fold and function with both disulfide bonds broken and the cysteine residues carboxymethylated. The large decrease in the stability is due mainly to an increase in the conformational entropy of the unfolded protein which results when the constraints of the disulfide bonds on the flexibility are removed. We propose a new equation for predicting the effect of a cross-link on the conformational entropy of a protein: delta Sconf = -2.1 - (3/2)R 1n n, where n is the number of residues between the side chains which are cross-linked. This equation gives much better agreement with experimental results than other forms of this equation which have been used previously.     

Early events in the disulfide-coupled folding of BPTI.

Recent studies of the refolding of reduced bovine pancreatic trypsin inhibitor (BPTI) have shown that a previously unidentified intermediate with a single disulfide is formed much more rapidly than any other one-disulfide species. This intermediate contains a disulfide that is present in the native protein (between Cys14 and 38), but it is thermodynamically less stable than the other two intermediates with single native disulfides. To characterize the role of the [14-38] intermediate and the factors that favor its formation, detailed kinetic and mutational analyses of the early disulfide-formation steps were carried out. The results of these studies indicate that the formation of [14-38] from the fully reduced protein is favored by both local electrostatic effects, which enhance the reactivities of the Cys14 and 38 thiols, and conformational tendencies that are diminished by the addition of urea and are enhanced at lower temperatures. At 25 degrees C and pH 7.3, approximately 35% of the reduced molecules were found to initially form the 14-38 disulfide, but the majority of these molecules then undergo intramolecular rearrangements to generate non-native disulfides, and subsequently the more stable intermediates with native disulfides. Amino acid replacements, other than those involving Cys residues, were generally found to have only small effects on either the rate of forming [14-38] or its thermodynamic stability, even though many of the same substitutions greatly destabilized the native protein and other disulfide-bonded intermediates. In addition, those replacements that did decrease the steady-state concentration of [14-38] did not adversely affect further folding and disulfide formation. These results suggest that the weak and transient interactions that are often detected in unfolded proteins and early folding intermediates may, in some cases, not persist or promote subsequent folding steps.     

A 3-disulfide mutant of mouse prion protein expression, oxidative folding, reductive unfolding, conformational stability, aggregation and isomerization

The structure of wild-type mouse prion protein mPrP(23-231) consists of two distinctive segments with approximately equal size, a disordered and flexible N-terminal domain encompassing residues 23-124 and a largely structured C-terminal domain containing about 40% of helical structure and stabilized by one disulfide bond (Cys(178)-Cys(213)). We have expressed a mPrP mutant with 4 Ala/Ser-->Cys replacements, two each at the N-(Cys(36), Cys(112)) and C-(Cys(134), Cys(169)) domains. Our specific aims are to study the interaction between N- and C-domains of mPrP during the oxidative folding and to produce stabilized isomers of mPrP for further analysis. Oxidative folding of fully reduced mutant, mPrP(6C), generates one predominant 3-disulfide isomer, designated as N-mPrP(3SS), which comprises the native disulfide (Cys(178)-Cys(213)) and two non-native disulfide bonds (Cys(36)-Cys(134) and Cys(112)-Cys(169)) that covalently connect the N- and C-domains. In comparison to wild-type mPrP(23-231), N-mPrP(3SS) exhibits an indistinguishable CD spectra, a similar conformational stability in the absence of thiol and a reduced ability to aggregate. In the presence of thiol catalyst and denaturant, N-mPrP(3SS) unfolds and generates diverse isomers that are amenable to further isolation, structural and functional analysis.     

NMR structure of a variant human prion protein with two disulfide bridges

The nuclear magnetic resonance structure of the globular domain with residues 121-230 of a variant human prion protein with two disulfide bonds, hPrP(M166C/E221C), shows the same global fold as wild-type hPrP(121-230). It contains three alpha-helices of residues 144-154, 173-194 and 200-228, an anti-parallel beta-sheet of residues 128-131 and 161-164, and the disulfides Cys166-Cys221 and Cys179-Cys214. The engineered extra disulfide bond in the presumed "protein X"-binding site is accommodated with slight, strictly localized conformational changes. High compatibility of hPrP with insertion of a second disulfide bridge in the protein X epitope was further substantiated by model calculations with additional variant structures. The ease with which the hPrP structure can accommodate a variety of locations for a second disulfide bond within the presumed protein X-binding epitope suggests a functional role for the extensive perturbation by a natural second disulfide bond of the corresponding region in the human doppel protein.     

The folding mechanism of a two-domain protein: folding kinetics and domain docking of the gene-3 protein of phage fd

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.     

A conformational unfolding reaction activates phage fd for the infection of Escherichia coli

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.     

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.     

Enthalpic and entropic contributions mediate the role of disulfide bonds on the conformational stability of interleukin-4

The role of disulfide bridges in the structure, stability, and folding pathways of proteins has been the subject of wide interest in the fields of protein design and engineering. However, the relative importance of entropic and enthalpic contributions for the stabilization of proteins provided by disulfides is not always clear. Here, we perform a detailed analysis of the role of disulfides in the conformational stability of human Interleukin-4 (IL4), a four-helix bundle protein. In order to evaluate the contribution of two out of the three disulfides to the structure and stability of IL4, two IL4 mutants, C3T-IL4 and C24T-IL4, were used. NMR and ANS binding experiments were compatible with altered dynamics and an increase of the nonpolar solvent-accessible surface area of the folded state of the mutant proteins. Chemical and thermal unfolding experiments followed by fluorescence and circular dichroism revealed that both mutant proteins have lower conformational stability than the wild-type protein. Transition temperatures of unfolding decreased 14 degrees C for C3T-IL4 and 10 degrees C for C24T-IL4, when compared to WT-IL4, and the conformational stability, at 25 degrees C, decreased 4.9 kcal/mol for C3T-IL4 and 3.2 kcal/mol for C24T-IL4. Interestingly, both the enthalpy and the entropy of unfolding, at the transition temperature, decreased in the mutant proteins. Moreover, a smaller change in heat capacity of unfolding was also observed for the mutants. Thus, disulfide bridges in IL4 play a critical role in maintaining the thermodynamic stability and core packing of the helix bundle.     

Elimination of a cis-Proline-Containing Loop and Turn Optimization Stabilizes a Protein and Accelerates Its Folding

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.     

A phage display system for studying the sequence determinants of protein folding.

We have developed a phage display system that provides a means to select variants of the IgG binding domain of peptostreptococcal protein L that fold from large combinatorial libraries. The premise underlying the selection scheme is that binding of protein L to IgG requires that the protein be properly folded. Using a combination of molecular biological and biophysical methods, we show that this assumption is valid. First, the phage selection procedure strongly selects against a point mutation in protein L that disrupts folding but is not in the IgG binding interface. Second, variants recovered from a library in which the first third of protein L was randomized are properly folded. The degree of sequence variation in the selected population is striking: the variants have as many as nine substitutions in the 14 residues that were mutagenized. The approach provides a selection for "foldedness" that is potentially applicable to any small binding protein.     

Disulfide bonds of phospholipase A2 from bee venom yield discrete contributions to its conformational stability

Disulfide bonds are known to be crucial for protein stability. To probe the contribution of each of the five disulfide bonds (C9-C31, C30-C70, C37-C63, C61-C95, and C105-C113) in bee venom phospholipase A₂ to stability, variants with deleted disulfide bonds were produced by substituting two serine residues for each pair of cysteine residues. The mutations started from the pseudo-wild-type variant (pWT) with the mutation I1A (Markert et al., Biotechnol. Bioeng. 98 (2007) 48-59). All variants were expressed in Escherichia coli, refolded from inclusion bodies and purified as pWT. The activity of the variants ranged from 12 to 82% of pWT. From the transition curves of guanidine hydrochloride-induced unfolding, the contributions of the individual disulfide bonds to conformational stability were estimated. They increased in the sequence C9-C31 < C105-C113 < C30-C70 ≈ C37-C63 < C61-C95. For two disulfide bonds (C9-C31, C105-C113) the effects were confirmed on additionally produced variants with the substitution of cysteine by alanine. Despite distinct differences in stability, all variants showed similar cooperativity in unfolding. Selected variants were also probed for proteolytic stability toward thermolysin. The removal of disulfide bonds increased the proteolytic susceptibility of the native proteins in the same way as the stability decreased. From the comparison of the results with literature data on phospholipase A₂ from bovine pancreas possessing seven disulfide bonds, it was concluded that conserved disulfide bonds in homologous proteins fulfill related functions in conformational stability.