Reversible Unfolding of a Thermophilic Membrane Protein in Phospholipid/Detergent Mixed Micelles

Folding mechanisms and stability of membrane proteins are poorly understood because of the known difficulties in finding experimental conditions under which reversible denaturation could be possible. In this work, we describe the equilibrium unfolding of Archaeoglobus fulgidus CopA, an 804-residue α-helical membrane protein that is involved in transporting Cu⁺ throughout biological membranes. The incubation of CopA reconstituted in phospholipid/detergent mixed micelles with high concentrations of guanidinium hydrochloride induced a reversible decrease in fluorescence quantum yield, far-UV ellipticity, and loss of ATPase and phosphatase activities. Refolding of CopA from this unfolded state led to recovery of full biological activity and all the structural features of the native enzyme. CopA unfolding showed typical characteristics of a two-state process, with ΔG_w° = 12.9 kJ mol^− 1, m = 4.1 kJ mol^− 1 M^− 1, C_m = 3 M, and ΔCp_w° = 0.93 kJ mol^− 1 K^− 1. These results point out to a fine-tuning mechanism for improving protein stability. Circular dichroism spectroscopic analysis of the unfolded state shows that most of the secondary and tertiary structures were disrupted. The fraction of Trp fluorescence accessible to soluble quenchers shifted from 0.52 in the native state to 0.96 in the unfolded state, with a significant spectral redshift. Also, hydrophobic patches in CopA, mainly located in the transmembrane region, were disrupted as indicated by 1-anilino-naphtalene-8-sulfonate fluorescence. Nevertheless, the unfolded state had a small but detectable amount of residual structure, which might play a key role in both CopA folding and adaptation for working at high temperatures.     

The Contribution of a Covalently Bound Cofactor to the Folding and Thermodynamic Stability of an Integral Membrane Protein

The factors controlling the stability, folding, and dynamics of integral membrane proteins are not fully understood. The high stability of the membrane protein bacteriorhodopsin (bR), an archetypal member of the rhodopsin photoreceptor family, has been ascribed to its covalently bound retinal cofactor. We investigate here the role of this cofactor in the thermodynamic stability and folding kinetics of bR. Multiple spectroscopic probes were used to determine the kinetics and energetics of protein folding in mixed lipid/detergent micelles in the presence and absence of retinal. The presence of retinal increases extrapolated values for the overall unfolding free energy from 6.3 ± 0.4 kcal mol^− 1 to 23.4 ± 1.5 kcal mol^− 1 at zero denaturant, suggesting that the cofactor contributes 17.1 kcal mol^− 1 towards the overall stability of bR. In addition, the cooperativity of equilibrium unfolding curves is markedly reduced in the absence of retinal with overall m-values decreasing from 31.0 ± 2.0 kcal mol^− 1 to 10.9 ± 1.0 kcal mol^− 1, indicating that the folded state of the apoprotein is less compact than the equivalent for the holoprotein. This change in the denaturant response means that the difference in the unfolding free energy at a denaturant concentration midway between the two unfolding curves is only ca 3-6 kcal mol^− 1. Kinetic data show that the decrease in stability upon removal of retinal is associated with an increase in the apparent intrinsic rate constant of unfolding, k_u^H2O, from ~1 × 10^− 16 s^− 1 to ~1 × 10^− 4 s^− 1 at 25 °C. This correlates with a decrease in the unfolding activation energy by 16.3 kcal mol^− 1 in the apoprotein, extrapolated to zero SDS. These results suggest that changes in bR stability induced by retinal binding are mediated solely by changes in the activation barrier for unfolding. The results are consistent with a model in which bR is kinetically stabilized via a very slow rate of unfolding arising from protein-retinal interactions that increase the rigidity and compactness of the polypeptide chain.     

Desolvation and Development of Specific Hydrophobic Core Packing during Im7 Folding

Development of a tightly packed hydrophobic core drives the folding of water-soluble globular proteins and is a key determinant of protein stability. Despite this, there remains much to be learnt about how and when the hydrophobic core becomes desolvated and tightly packed during protein folding. We have used the bacterial immunity protein Im7 to examine the specificity of hydrophobic core packing during folding. This small, four-helix protein has previously been shown to fold via a compact three-helical intermediate state. Here, overpacking substitutions, in which residue side-chain size is increased, were used to examine the specificity and malleability of core packing in the folding intermediate and rate-limiting transition state. In parallel, polar groups were introduced into the Im7 hydrophobic core via Val→Thr or Phe→Tyr substitutions and used to determine the solvation status of core residues at different stages of folding. Over 30 Im7 variants were created allowing both series of substitutions to cover all regions of the protein structure. Φ-value analysis demonstrated that the major changes in Im7 core solvation occur prior to the population of the folding intermediate, with key regions involved in docking of the short helix III remaining solvent-exposed until after the rate-limiting transition state has been traversed. In contrast, overpacking core residues revealed that some regions of the native Im7 core are remarkably malleable to increases in side-chain volume. Overpacking residues in other regions of the Im7 core result in substantial (> 2.5 kJ mol^− 1) destabilisation of the native structure or even prevents efficient folding to the native state. This study provides new insights into Im7 folding; demonstrating that whilst desolvation occurs early during folding, adoption of a specifically packed core is achieved only at the very last step in the folding mechanism.     

Relative Domain Folding and Stability of a Membrane Transport Protein

There is a limited understanding of the folding of multidomain membrane proteins. Lactose permease (LacY) of Escherichia coli is an archetypal member of the major facilitator superfamily of membrane transport proteins, which contain two domains of six transmembrane helices each. We exploit chemical denaturation to determine the unfolding free energy of LacY and employ Trp residues as site-specific thermodynamic probes. Single Trp LacY mutants are created with the individual Trps situated at mirror image positions on the two LacY domains. The changes in Trp fluorescence induced by urea denaturation are used to construct denaturation curves from which unfolding free energies can be determined. The majority of the single Trp tracers report the same stability and an unfolding free energy of approximately + 2 kcal mol^− 1. There is one exception; the fluorescence of W33 at the cytoplasmic end of helix I on the N domain is unaffected by urea. In contrast, the equivalent position on the first helix, VII, of the C-terminal domain exhibits wild-type stability, with the single Trp tracer at position 243 on helix VII reporting an unfolding free energy of + 2 kcal mol^− 1. This indicates that the region of the N domain of LacY at position 33 on helix I has enhanced stability to urea, when compared the corresponding location at the start of the C domain. We also find evidence for a potential network of stabilising interactions across the domain interface, which reduces accessibility to the hydrophilic substrate binding pocket between the two domains.     

SEA domain autoproteolysis accelerated by conformational strain: energetic aspects

A subclass of proteins with the SEA (sea urchin sperm protein, enterokinase, and agrin) domain fold exists as heterodimers generated by autoproteolytic cleavage within a characteristic G^− 1S^+ 1VVV sequence. Autoproteolysis occurs by a nucleophilic attack of the serine hydroxyl on the vicinal glycine carbonyl followed by an N → O acyl shift and hydrolysis of the resulting ester. The reaction has been suggested to be accelerated by the straining of the scissile peptide bond upon protein folding. In an accompanying article, we report the mechanism; in this article, we provide further key evidence and account for the energetics of coupled protein folding and autoproteolysis. Cleavage of the GPR116 domain and that of the MUC1 SEA domain occur with half-life (t_½) values of 12 and 18 min, respectively, with lowering of the free energy of the activation barrier by ∼ 10 kcal mol^− 1 compared with uncatalyzed hydrolysis. The free energies of unfolding of the GPR116 and MUC1 SEA domains were measured to ∼ 11 and ∼ 15 kcal mol^− 1, respectively, but ∼ 7 kcal mol^− 1 of conformational energy is partitioned as strain over the scissile peptide bond in the precursor to catalyze autoproteolysis by substrate destabilization. A straining energy of ∼ 7 kcal mol^− 1 was measured by using both a pre-equilibrium model to analyze stability and cleavage kinetics data obtained with the GPR116 SEA domain destabilized by core mutations or urea addition, as well as the difference in thermodynamic stabilities of the MUC1 SEA precursor mutant S1098A (with a G^− 1A^+ 1VVV motif) and the wild-type protein. The results imply that cleavage by N → O acyl shift alone would proceed with a t_½ of ∼ 2.3 years, which is too slow to be biochemically effective. A subsequent review of structural data on other self-cleaving proteins suggests that conformational strain of the scissile peptide bond may be a common mechanism of autoproteolysis.     

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.     

Conformational stability of α-amylase from malted sorghum (Sorghum bicolor): Reversible unfolding by denaturants

α-Amylase from Sorghum bicolor, is reversibly unfolded by chemical denaturants at pH 7.0 in 50 mM Hepes containing 13.6 mM calcium and 15 mM DTT. The isothermal equilibrium unfolding at 27 °C is characterized by two state transition with ΔG (H₂O) of 16.5 kJ mol⁻¹ and 22 kJ mol⁻¹, respectively, at pH 4.8 and pH 7.0 for GuHCl and ΔG (H₂O) of 25.2 kJ mol⁻¹ at pH 4.8 for urea. The conformational stability indicators such as the change in excess heat capacity (ΔC_p), the unfolding enthalpy (H_g) and the temperature at ΔG = 0 (T_g) are 17.9 ± 0.7 kJ mol⁻¹ K⁻¹, 501.2 ± 18.2 kJ mol⁻¹ and 337.3 ± 6.9 K at pH 4.8 and 14.3 ± 0.5 kJ mol⁻¹ K⁻¹, 509.3 ± 21.7 kJ mol⁻¹ and 345.4 ± 4.8 K at pH 7.0, respectively. The reactivity of the conserved cysteine residues, during unfolding, indicates that unfolding starts from the ‘B’ domain of the enzyme. The oxidation of cysteine residues, during unfolding, can be prevented by the addition of DTT. The conserved cysteine residues are essential for enzyme activity but not for the secondary and tertiary fold acquired during refolding of the denatured enzyme. The pH dependent stability described by ΔG (H₂O) and the effect of salt on urea induced unfolding confirm the role of electrostatic interactions in enzyme stability.     

Osmolyte effect on the stability and folding of a hyperthermophilic protein

Proteins are known to be stabilized by naturally occurring osmolytes such as amino acids, sugars, and methylamines. Here, we examine the effect of trimethylamine-N-oxide (TMAO) on the conformational stability of ribonuclease HII from a hyperthermophile, Thermococcus kodakaraensis (Tk-RNase HII), which inherently possesses high conformational stability. Heat- and guanidine hydrochloride-induced unfolding experiments demonstrated that the conformational stability of Tk-RNase HII in the presence of 0.5M TMAO was higher than that in the absence of TMAO at all examined temperatures. TMAO affected the unfolding and refolding kinetics of Tk-RNase HII to a similar extent. These results indicate that proteins are universally stabilized by osmolytes, regardless of their robustness, and suggest a stabilization mechanism by osmolytes, caused by the unfavorable interaction of osmolytes with protein backbones in the denatured state. Our results also imply that the basic protein folding principle is not dependent on protein stability and evolution.     

Kinetically robust monomeric protein from a hyperthermophile

Equilibrium and kinetic studies were carried out under denaturation conditions to clarify the energetic features of the high stability of a monomeric protein, ribonuclease HII, from a hyperthermophile, Thermococcus kodakaraensis (Tk-RNase HII). Guanidine hydrochloride (GdnHCl)-induced unfolding and refolding were measured with circular dichroism at 220 nm, and heat-induced denaturation was studied with differential scanning calorimetry. Both GdnHCl- and heat-induced denaturation are very reversible. It was difficult to obtain the equilibrated unfolding curve of Tk-RNase HII below 40 degrees C, because of the remarkably slow unfolding. The two-state unfolding and refolding reactions attained equilibrium at 50 degrees C after 2 weeks. The Gibbs energy change of GdnHCl-induced unfolding (DeltaG(H(2)O)) at 50 degrees C was 43.6 kJ mol(-1). The denaturation temperature in the DSC measurement shifted as a function of the scan rate; the denaturation temperature at a scan rate of 90 degrees C h(-1) was higher than at a scan rate of 5 degrees C h(-1). The unfolding and refolding kinetics of Tk-RNase HII were approximated as a first-order reaction. The ln k(u) and ln k(r) values depended linearly on the denaturant concentration between 10 and 50 degrees C. The DeltaG(H(2)O) value obtained from the rate constant in water using the two-state model at 50 degrees C, 44.5 kJ mol(-1), was coincident with that from the equilibrium study, 43.6 kJ mol(-1), suggesting the two-state folding of Tk-RNase HII. The values for the rate constant in water of the unfolding for Tk-RNase HII were much smaller than those of E. coli RNase HI and Thermus thermophilus RNase HI, which has a denaturation temperature similar to that of Tk-RNase HII. In contrast, little difference was observed in the refolding rates among these proteins. These results indicate that the stabilization mechanism of monomeric protein from a hyperthermophile, Tk-RNase HII, with reversible two-state folding is characterized by remarkably slow unfolding.     

The folding process of acylphosphatase from Escherichia coli is remarkably accelerated by the presence of a disulfide bond

The acylphosphatase from Escherichia coli (EcoAcP) is the first AcP so far studied with a disulfide bond. A mutational variant of the enzyme lacking the disulfide bond has been produced by substituting the two cysteine residues with alanine (EcoAcP mutational variant C5A/C49A, mutEcoAcP). The native states of the two protein variants are similar, as shown by far-UV and near-UV circular dichroism and dynamic light-scattering measurements. From unfolding experiments at equilibrium using intrinsic fluorescence and far-UV circular dichroism as probes, EcoAcP shows an increased conformational stability as compared with mutEcoAcP. The wild-type protein folds according to a two-state model with a very fast rate constant (k_F^H2O = 72,600 s^− 1), while mutEcoAcP folds ca 1500-fold slower, via the accumulation of a partially folded species. The correlation between the hydrophobicity of the polypeptide chain and the folding rate, found previously in the AcP-like structural family, is maintained only when considering the mutant but not the wild-type protein, which folds much faster than expected from this correlation. Similarly, the correlation between the relative contact order and the folding rate holds only for mutEcoAcP. The correlation also holds for EcoAcP, provided the relative contact order value is recalculated by considering the disulfide bridge as an alternate path for the backbone to determine the shortest sequence separation between contacting residues. These results indicate that the presence of a disulfide bond in a protein is an important determinant of the folding rate and allows its contribution to be determined in quantitative terms.     

Proline effect on the thermostability and slow unfolding of a hyperthermophilic protein

Ribonuclease HII from hyperthermophile Thermococcus kodakaraensis (Tk-RNase HII) is a robust monomeric protein under kinetic control, which possesses some proline residues at the N-terminal of alpha-helices. Proline residue at the N-terminal of an alpha-helix is thought to stabilize a protein. In this work, the thermostability and folding kinetics of Tk-RNase HII were measured for mutant proteins in which a proline residue is introduced (Xaa to Pro) or removed (Pro to Ala) at the N-terminal of alpha-helices. In the folding experiments, the mutant proteins examined exhibit little influence on the remarkably slow unfolding of Tk-RNase HII. In contrast, E111P and K199P exhibit some thermostabilization, whereas P46A, P70A and P174A have some thermodestabilization. E111P/K199P and P46A/P70A double mutations cause cumulative changes in stability. We conclude that the proline effect on protein thermostability is observed in a hyperthermophilic protein, but each proline residue at the N-terminal of an alpha-helix slightly contributes to the thermostability. The present results also mean that even a natural hyperthermophilic protein can acquire improved thermostability.     

Molecular rearrangement in POR macrodomains as a reason for the blue shift of chlorophyllide fluorescence observed after phototransformation

The activation energy and activation volume of the spectral blue shift subsequent to protochlorophyllide phototransformation (called Shibata shift in intact leaves) were studied in prolamellar body (PLB) and prothylakoid-(PT)-enriched membrane fractions prepared from dark-grown wheat (Triticum aestivum, L.) leaves. The measurements were done at 20, 30 and 40 °C and at various pressure values. The activation energy values were 181 ± 8 kJ mol^− 1 and 188 ± 6 kJ mol^− 1 for the PLBs and the PTs, respectively. The pressure stabilized the structure of the NADPH:protochlorophyllide oxidoreductase (POR) macrodomains; it prevented or slowed down the blue shift. There were no significant differences between the activation volumes of PLBs and PTs at 30 or 40 °C giving values around 100-125 ml mol^− 1 which correspond to changes in the tertiary structure of proteins but also resemble the volume changes occurring during the disaggregation of protein dimers or oligomers, or during dissociation of peripheral membrane proteins from membranes. The small differences in the activation parameters of PLBs and PTs indicate that molecular rearrangements inside the POR macrodomains are the primary reasons of the fluorescence blue shift; however, their lipid microenvironment must be also important in the initialization of the shift.     

NMR unfolding studies on a liver bile acid binding protein reveal a global two-state unfolding and localized singular behaviors

The folding properties of a bile acid binding protein, belonging to a subfamily of the fatty acid binding proteins, have been here investigated both by hydrogen exchange measurements, using the SOFAST NMR approach, and urea denaturation experiments. The urea unfolding profiles of individual residues, acting as single probes, were simultaneously analyzed through a global fit, according to a two-state unfolding model. The resulting conformational stability ΔG_U(H₂O) = 7.2 ± 0.25 kcal mol⁻¹ is in good agreement with hydrogen exchange stability ΔG_op. While the majority of protein residues satisfy this model, few amino-acids display a singular behavior, not directly amenable to the presence of a folding intermediate, as reported for other fatty acid binding proteins. These residues are part of a protein patch characterized by enhanced plasticity. To explain this singular behavior a tentative model has been proposed which takes into account the interplay between the dynamic features and the formation of transient aggregates. A functional role for this plasticity, related to translocation across the nuclear membrane, is discussed.     

Structure, folding and stability of FimA, the main structural subunit of type 1 pili from uropathogenic Escherichia coli strains

Filamentous type 1 pili are responsible for attachment of uropathogenic Escherichia coli strains to host cells. They consist of a linear tip fibrillum and a helical rod formed by up to 3000 copies of the main structural pilus subunit FimA. The subunits in the pilus interact via donor strand complementation, where the incomplete, immunoglobulin-like fold of each subunit is complemented by an N-terminal donor strand of the subsequent subunit. Here, we show that folding of FimA occurs at an extremely slow rate (half-life: 1.6 h) and is catalyzed more than 400-fold by the pilus chaperone FimC. Moreover, FimA is capable of intramolecular self-complementation via its own donor strand, as evidenced by the loss of folding competence upon donor strand deletion. Folded FimA is an assembly-incompetent monomer of low thermodynamic stability (− 10.1 kJ mol^− 1) that can be rescued for pilus assembly at 37 °C because FimC selectively pulls the fraction of unfolded FimA molecules from the FimA folding equilibrium and allows FimA refolding on its surface. Elongation of FimA at the C-terminus by its own donor strand generated a self-complemented variant (FimAa) with alternative folding possibilities that spontaneously adopts the more stable conformation (− 85.0 kJ mol^− 1) in which the C-terminal donor strand is inserted in the opposite orientation relative to that in FimA. The solved NMR structure of FimAa revealed extensive β-sheet hydrogen bonding between the FimA pilin domain and the C-terminal donor strand and provides the basis for reconstruction of an atomic model of the pilus rod.     

Dimer Formation of a Stabilized Gβ1 Variant: A Structural and Energetic Analysis

In previous work, a strongly stabilized variant of the β1 domain of streptococcal protein G (Gβ1) was obtained by an in vitro selection method. This variant, termed Gβ1-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 Gβ1 gene libraries were performed, and the crystal structure of Gβ1-M2 was determined. All single substitutions (E15V, T16L, T18I, and N37L) stabilize wild-type Gβ1 by contributions of between 1.6 and 6.0 kJ mol^− 1 (at 70 °C). Hydrophobic residues at positions 16 and 37 provide the major contribution to stabilization by enlarging the hydrophobic core of Gβ1. 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-Å crystal structure of Gβ1-M2 and NMR measurements in solution provide the explanation for the observed dimer formation. It involves a head-to-head arrangement of two Gβ1-M2 molecules via six intermolecular hydrogen bonds between the two β 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° 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 Gβ1 variants.     

Identification of a Mechanical Rheostat in the Hydrophobic Core of Protein L

The ability of proteins and their complexes to withstand or respond to mechanical stimuli is vital for cells to maintain their structural organisation, to relay external signals and to facilitate unfolding and remodelling. Force spectroscopy using the atomic force microscope allows the behaviour of single protein molecules under an applied extension to be investigated and their mechanical strength to be quantified. protein L, a simple model protein, displays moderate mechanical strength and is thought to unfold by the shearing of two mechanical sub-domains. Here, we investigate the importance of side-chain packing for the mechanical strength of protein L by measuring the mechanical strength of a series of protein L variants containing single conservative hydrophobic volume deletion mutants. Of the five thermodynamically destabilised variants characterised, only one residue (I60V) close to the interface between two mechanical sub-domains was found to differ in mechanical properties to wild type (ΔF_I60V-WT = − 36 pN at 447 nm s^− 1, Δx_uI60V-WT = 0.2 nm). Φ-value analysis of the unfolding data revealed a highly native transition state. To test whether the number of hydrophobic contacts across the mechanical interface does affect the mechanical strength of protein L, we measured the mechanical properties of two further variants. protein L L10F, which increases core packing but does not enhance interfacial contacts, increased mechanical strength by 13 ± 11 pN at 447 nm s^− 1. By contrast, protein L I60F, which increases both core and cross-interface contacts, increased mechanical strength by 72 ± 13 pN at 447 nm s^− 1. These data suggest a method by which nature can evolve a varied mechanical response from a limited number of topologies and demonstrate a generic but facile method by which the mechanical strength of proteins can be rationally modified.     

Thermodynamic and kinetic characterization of apoHmpH, a fast-folding bacterial globin

Despite the widespread presence of the globin fold in most living organisms, only eukaryotic globins have been employed as model proteins in folding/stability studies so far. This work introduces the first thermodynamic and kinetic characterization of a prokaryotic globin, that is, the apo form of the heme-binding domain of flavohemoglobin (apoHmpH) from Escherichia coli. This bacterial globin has a widely different sequence but nearly identical structure to its eukaryotic analogues. We show that apoHmpH is a well-folded monomeric protein with moderate stability at room temperature [apparent ΔG°_UN(w) = − 3.1 ± 0.3 kcal mol^− 1; m_UN = − 1.7 kcal mol^− 1 M^− 1] and predominant α-helical structure. Remarkably, apoHmpH is the fastest-folding globin known to date, as it refolds about 4- to 16-fold more rapidly than its eukaryotic analogues (e.g., sperm whale apomyoglobin and soybean apoleghemoglobin), populating a compact kinetic intermediate (β_I = 0.9 ± 0.2) with significant helical content. Additionally, the single Trp120 (located in the native H helix) becomes locked into a fully native-like environment within 6 ms, suggesting that this residue and its closest spatial neighbors complete their folding at ultrafast (submillisecond) speed. In summary, apoHmpH is a bacterial globin that shares the general folding scheme (i.e., a rapid burst phase followed by slower rate-determining phases) of its eukaryotic analogues but displays an overall faster folding and a kinetic intermediate with some fully native-like traits. This study supports the view that the general folding features of bacterial and eukaryotic globins are preserved through evolution while kinetic details differ.     

The foldon substructure of staphylococcal nuclease

To search for submolecular foldon units, the spontaneous reversible unfolding and refolding of staphylococcal nuclease under native conditions was studied by a kinetic native-state hydrogen exchange (HX) method. As for other proteins, it appears that staphylococcal nuclease is designed as an assembly of well-integrated foldon units that may define steps in its folding pathway and may regulate some other functional properties. The HX results identify 34 amide hydrogens that exchange with solvent hydrogens under native conditions by way of large transient unfolding reactions. The HX data for each hydrogen measure the equilibrium stability (ΔG_HX) and the kinetic unfolding and refolding rates (k_op and k_cl) of the unfolding reaction that exposes it to exchange. These parameters separate the 34 identified residues into three distinct HX groupings. Two correspond to clearly defined structural units in the native protein, termed the blue and red foldons. The remaining HX grouping contains residues, not well separated by their HX parameters alone, that represent two other distinct structural units in the native protein, termed the green and yellow foldons. Among these four sets, a last unfolding foldon (blue) unfolds with a rate constant of 6 × 10^− 6 s^− 1 and free energy equal to the protein's global stability (10.0 kcal/mol). It represents part of the β-barrel, including mutually H-bonding residues in the β4 and β5 strands, a part of the β3 strand that H-bonds to β5, and residues at the N-terminus of the α2 helix that is capped by β5. A second foldon (green), which unfolds and refolds more rapidly and at slightly lower free energy, includes residues that define the rest of the native α2 helix and its C-terminal cap. A third foldon (yellow) defines the mutually H-bonded β1-β2-β3 meander, completing the native β-barrel, plus an adjacent part of the α1 helix. A final foldon (red) includes residues on remaining segments that are distant in sequence but nearly adjacent in the native protein. Although the structure of the partially unfolded forms closely mimics the native organization, four residues indicate the presence of some nonnative misfolding interactions. Because the unfolding parameters of many other residues are not determined, it seems likely that the concerted foldon units are more extensive than is shown by the 34 residues actually observed.     

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.     

Downhill versus Barrier-Limited Folding of BBL 1: Energetic and Structural Perturbation Effects upon Protonation of a Histidine of Unusually Low pKa

A dispersion of melting temperatures at pH 5.3 for individual residues of the BBL protein domain has been adduced as evidence for barrier-free downhill folding. Other members of the peripheral subunit domain family fold cooperatively at pH 7. To search for possible causes of anomalies in BBL's denaturation behavior, we measured the pH titration of individual residues by heteronuclear NMR. At 298 K, the pK_a of His142 was close to that of free histidine at 6.47 ± 0.04, while that of the more buried His166 was highly perturbed at 5.39 ± 0.02. Protonation of His166 is thus energetically unfavorable and destabilizes the protein by ∼ 1.5 kcal/mol. Changes in C^α secondary shifts at pH 5.3 showed a decrease in helicity of the C-terminus of helix 2, where His166 is located, which was accompanied by a measured decrease of 1.1 ± 0.2 kcal/mol in stability from pH 7 to 5.3. Protonation of His166 perturbs, therefore, the structure of BBL. Only ∼ 1% of the structurally perturbed state will be present at the biologically relevant pH 7.6. Experiments at pH 5.3 report on a near-equal mixture of the two different native states. Further, at this pH, small changes of pH and pK_a induced by changes in temperature will have near-maximal effects on pH-dependent conformational equilibria and on propagation of experimental error. Accordingly, conventional barrier-limited folding predicts some dispersion of measured thermal unfolding curves of individual residues at pH 5.3.