NMR monitoring of accumulation and folding of 15N-labeled protein overexpressed in Pichia pastoris.

Postgenomic studies have led to an increasing demand for isotope-labeled proteins. We present a method for producing large quantities of truly native (15)N-labeled protein. Based on the secretion capabilities of the yeast Pichia pastoris, the recombinant protein is easily purified in a single step as it is secreted. Control of all nitrogen sources permits very high labeling yields. As a result, accumulation and folding of the recombinant protein can be monitored by heteronuclear NMR without purification. Comparison of sample spectra with the spectrum of the purified recombinant protein allows detection of the secreted protein in the culture and monitoring of its folding, from the start of the induction phase. The detection limit for a (15)N-labeled protein is estimated as 20 microM and corresponds, for a 10-kDa protein, to a load of 40 mg/liter in the fermentor. This concentration is reached by most reported preparations in P. pastoris. Further concentration by ultrafiltration would compensate for lower production. This procedure may be useful in many structural genomics and combinatorial chemistry screening projects where most protein productions meet the requirements for this method.     

Unfolding Studies of <Emphasis Type="Italic">Escherichia coli</Emphasis> Maltodextrin Glucosidase Monitored by Fluorescence Spectroscopy

Equilibrium unfolding of a 69-kDa monomeric Escherichia coli maltodextrin glucosidase (MalZ) was studied using intrinsic and extrinsic fluorescence spectroscopy. The unfolding transition of MalZ followed a three-state process, involving the formation of a stable intermediate state having more exposed hydrophobic surface. It was found that the protein structure can be easily perturbed by low concentration of guanidium hydrochloride (GdnHCl) and, at a GdnHCl concentration of 2 M, MalZ was denatured completely. The active site of the protein also has been proved to be sensitive to a low concentration of GdnHCl since MalZ deactivated at 0.5 M GdnHCl completely. The surface hydrophobicity and ANS-binding site of the protein have been determined to be 150.7 and 0.24, respectively. Perhaps the formation of the stable unfolding intermediate, having higher surface hydrophobicity, may be one of the reasons for aggregation of MalZ and its recognition by chaperonin GroEL during the assisted folding pathway.     

Energy barriers, cooperativity, and hidden intermediates in the folding of small proteins

Current theoretical views of the folding process of small proteins (< approximately 100 amino acids) postulate that the landscape of potential mean force (PMF) for the formation of the native state has a funnel shape and that the free energy barrier to folding arises from the chain configurational entropy only. However, recent theoretical studies on the formation of hydrophobic clusters with explicit water suggest that a barrier should exist on the PMF of folding, consistent with the fact that protein folding generally involves a large positive activation enthalpy at room temperature. In addition, high-resolution structural studies of the hidden partially unfolded intermediates have revealed the existence of non-native interactions, suggesting that the correction of the non-native interactions during folding should also lead to barriers on PMF. To explore the effect of a PMF barrier on the folding behavior of proteins, we modified Zwanzig's model for protein folding with an uphill landscape of PMF for the formation of transition states. We found that the modified model for short peptide segments can satisfy the thermodynamic and kinetic criteria for an apparently two-state folding. Since the Levinthal paradox can be solved by a stepwise folding of short peptide segments, a landscape of PMF with a locally uphill search for the transition state and cooperative stabilization of folding intermediates/native state is able to explain the available experimental results for small proteins. We speculate that the existence of cooperative hidden folding intermediates in small proteins could be the consequence of the highly specific structures of the native state, which are selected by evolution to perform specific functions and fold in a biologically meaningful time scale.     

Distinct roles of the N-terminal and C-terminal precursor domains in the biogenesis of the Bordetella pertussis filamentous hemagglutinin.

The 220-kDa Bordetella pertussis filamentous hemagglutinin (FHA) is the major exported protein found in culture supernatants. The structural gene of FHA has a coding potential for a 367-kDa protein, and the mature form constitutes the N-terminal 60% of the 367-kDa precursor. The C-terminal domain of the precursor was found to be important for the high-level secretion of full-length FHA but not of truncated analogs (80 kDa or less). The secretion of full-length and truncated FHA polypeptides requires the presence of the approximately 100-amino-acid N-terminal domain and the outer membrane protein FhaC, homologous to the N-terminal domains of the Serratia marcescens and Proteus mirabilis hemolysins and their accessory proteins, respectively. By analogy to these hemolysins, it is likely that the N-terminal domain of the FHA precursor interacts, directly or indirectly, with the accessory protein during FHA biogenesis. However, immunogenicity and antigenicity studies suggest that the N-terminal domain of FHA is masked by its C-terminal domain and therefore should not be available for its interactions with FhaC. These observations suggest a model in which the C-terminal domain of the FHA precursor may play a role as an intramolecular chaperone to prevent premature folding of the protein. Both heparin binding and hemagglutination are expressed by the N-terminal half of FHA, indicating that this domain contains important functional regions of the molecule.     

Nucleation and the transition state of the SH3 domain

We present a verified computational model of the SH3 domain transition state (TS) ensemble. This model was built for three separate SH3 domains using experimental phi-values as structural constraints in all-atom protein folding simulations. While averaging over all conformations incorrectly considers non-TS conformations as transition states, quantifying structures as pre-TS, TS, and post-TS by measurement of their transmission coefficient ("probability to fold", or p(fold)) allows for rigorous conclusions regarding the structure of the folding nucleus and a full mechanistic analysis of the folding process. Through analysis of the TS, we observe a highly polarized nucleus in which many residues are solvent-exposed. Mechanistic analysis suggests the hydrophobic core forms largely after an early nucleation step. SH3 presents an ideal system for studying the nucleation-condensation mechanism and highlights the synergistic relationship between experiment and simulation in the study of protein folding.     

Isolation of point mutations that affect the folding of the H chain of human ferritin in E.coli.

We have approached the problem of folding and assembly of the heavy (H) chain of human ferritin by isolating point mutations that affect this process. Apoferritin is an ideal model system to approach the problem of protein folding and assembly into multimeric structures. We have developed a recombinant hybrid molecule that allows us to select for ferritin mutants in which the folding-assembly process is altered or completely impaired. The selection procedure is based on a recombinant protein which consists of a fusion between the H chain of human ferritin and the alpha-peptide of beta-galactosidase. In the wild type situation, the alpha-peptide domain is segregated inside the apoferritin shell upon assembly and is unable to interact with the substrate and perform its enzymic function. We show that by selecting for mutations that restore beta-galactosidase activity we are able to identify ferritin mutations that affect the folding-assembly process. The selective procedure was applied to the analysis of the amino acid side chains that are important for the attainment of the correct conformation of the carboxy-terminal E helix in the 4-fold axis.     

Protein folding stabilizing time measurement: a direct folding process and three-dimensional random walk simulation

Protein particles undergo Brownian motion and collisions in solution. The diffusive collisions may lead to aggregation. For proteins to fold successfully the process has to occur quickly and before significant collision takes place. The speed of protein folding was deduced by studying the correlation time of a lysozyme refolding process from autocorrelation function analysis of the mean collision time and aggregation/soluble ratio of protein. It is a measure of time before which an aggregate can be formed and also is the time measure for a protein to fold into a stable state. We report on the protein folding stabilizing time of a lysozyme system to be 25.5-27.5 micros (<+/-4%) between 295 and 279K via direct folding experimental studies, supported by a three-dimensional random walk simulation of diffusion-limited aggregation model. Aggregation is suppressed when the protein is folded to a stable form. Spontaneous folding and diffusion-limited aggregation are antagonistic in nature. Meanwhile, the resultant aggresome, suggested by Raman and mass spectroscopy, may be formed by cross-linkages of disulfide bonds and hydrophobic interactions.     

Sub-microsecond protein folding

We have investigated the structure, equilibria, and folding kinetics of an engineered 35-residue subdomain of the chicken villin headpiece, an ultrafast-folding protein. Substitution of two buried lysine residues by norleucine residues stabilizes the protein by 1 kcal/mol and increases the folding rate sixfold, as measured by nanosecond laser T-jump. The folding rate at 300 K is (0.7 micros)(-1) with little or no temperature dependence, making this protein the first sub-microsecond folder, with a rate only twofold slower than the theoretically predicted speed limit. Using the 70 ns process to obtain the effective diffusion coefficient, the free energy barrier height is estimated from Kramers theory to be less than approximately 1 kcal/mol. X-ray crystallographic determination at 1A resolution shows no significant change in structure compared to the single-norleucine-substituted molecule and suggests that the increased stability is electrostatic in origin. The ultrafast folding rate, very accurate X-ray structure, and small size make this engineered villin subdomain an ideal system for simulation by atomistic molecular dynamics with explicit solvent.     

A novel method for study of protein folding kinetics by monitoring diffusion coefficient in time domain

Molecular diffusion process after the photo-induced electron injection to ferric cytochrome c (Fe(III) cyt c) in guanidine hydrochloride (GdnHCl) 3.5 M buffer solution is studied by the time-resolved transient grating technique. Circular dichroism studies have revealed that Fe(III) cyt c is unfolded under this condition but the reduced form, Fe(II) cyt c, is folded. Hence, this pulsed laser-induced reduction should initiate the folding process of cyt c. The observed transient grating signal shows prominent features, which have never been observed before. Based on several characteristic points, we concluded that the apparent diffusion coefficient (D) of Fe(II) cyt c after the reduction is time dependent, which must be associated with the protein folding dynamics. This time-dependent apparent D should reflect either the continuous time development of the hydrodynamic radius or population change of the unfolded and folded states during the folding dynamics. This is the first observation of the time-dependent apparent D during any chemical reaction, and this time-dependent measurement of D should be a unique and powerful way to study the protein folding kinetics from a viewpoint of the protein's shape or the protein-water intermolecular interaction.     

Recombinant 70-kDa protein from the amino-terminal region of rat fibronectin inhibits binding of fibronectin to cells and bacteria

Binding of fibronectin to substrate-attached cells and to Staphylococcus aureus is mediated by the amino-terminal 70-kDa portion of fibronectin. The 70-kDa amino-terminus is composed of nine type I and two type II internal homology units, each containing two intrachain disulfide bonds. The exact structural features of the 70-kDa amino-terminus that are necessary for binding to cells and bacteria are not known. We characterized a recombinant 70-kDa protein from the amino-terminus of rat fibronectin using a baculovirus expression system. Recombinant 70-kDa (r70kDa) protein was easily purified in high amounts from the conditioned medium by affinity chromatography on gelatin-agarose. Secretion was much less when N-linked glycosylation was blocked by tunicamycin. Like the native fragment, the r70kDa protein contains intrachain disulfide bonds. In addition, the r70kDa protein was indistinguishable from the nonrecombinant 70-kDa fragment in its ability to compete for binding sites on fibroblasts and S. aureus. Thus, the r70kDa protein retains the important functional characteristics of the native fragment. This expression system is well adapted to studying the structural features important for the interaction of 70-kDa protein with cells.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

The nature of protein folding pathways: The classical versus the new view

Summary Pulsed hydrogen exchange and other studies of the kinetic refolding pathways of several small proteins have established that folding intermediates with native-like secondary structures are well populated, but these studies have also shown that the folding kinetics are not well synchronized. Older studies of the kinetics of formation of the native protein, monitored by optical probes, indicate that the folding kinetics should be synchronized. The model commonly used in these studies is the simple sequential model, which postulates a unique folding pathway with defined and sequential intermediates. Theories of the folding process and Monte Carlo simulations of folding suggest that neither the folding pathway nor the set of folding intermediates is unique, and that folding intermediates accumulate because of kinetic traps caused by partial misfolding. Recent experiments with cytochrome c lend support to this new view of folding pathways. These different views of the folding process are discussed. Misfolding and consequent slowing down of the folding process as a result of cis-trans isomerization about prolyl peptide bonds in the unfolded protein are well known; isomerization occurs before refolding is initiated. The occurrence of equilibrium intermediates on the kinetic folding pathways of some proteins, such as -lactalbumin and apomyoglobin, argues that these intermediates are not caused by kinetic traps but rather are stable intermediates under certain conditions, and this conclusion is consistent with a sequential model of folding. Folding reactions with successive kinetic intermediates, in which late intermediates are more highly folded than early intermediates, indicate that folding is hierarchical. New experiments that test the predictions of the classical and the new views are needed.     

Pressure-exploration of the 33-kDa protein from the spinach photosystem II particle.

The 33-kDa protein isolated from the spinach photosystem II particle is an ideal model to explore high-pressure protein-unfolding. The protein has a very low free energy as previously reported by chemical unfolding studies, suggesting that it must be easy to modulate its unfolding transition by rather mild pressure. Moreover, the protein molecule consists of only one tryptophan residue (Trp241) and eight tyrosine residues, which can be conveniently used to probe the protein conformation and structural changes under pressure using either fluorescence spectroscopy or fourth derivative UV absorbance spectroscopy. The different experimental methods used in the present study indicate that at 20 degrees C and pH 6, the 33-kDa protein shows a reversible two-state unfolding transition from atmospheric pressure to about 180 MPa. This value is much lower than those found for the unfolding of most proteins studied so far. The unfolding transition induces a large red shift of the maximum fluorescence emission of 34 nm (from 316 nm to 350 nm). The change in standard free energy (DeltaGo) and in volume (DeltaV) for the transition at pH 6.0 and 20 degrees C are -14.6 kJ.mol-1 and -120 mL.mol-1, respectively, in which the DeltaGo value is consistent with that obtained by chemical denaturation. We found that pressure-induced protein unfolding is promoted by elevated temperatures, which seem largely attributed to the decrease in the absolute value of DeltaGo (only a minor variation was observed for the DeltaV value). However, the promotion of the unfolding by alkaline pH seems mainly related to the increase in DeltaV without any significant changes in DeltaGo. It was also found that NaCl significantly protects the protein from pressure-induced unfolding. In the presence of 1 M NaCl, the pressure needed to induce the half-unfold of the protein is shifted to a higher value (shift of 75 MPa) in comparison with that observed without NaCl. Interestingly, in the presence of NaCl, the value of DeltaV is significantly reduced whilst that of DeltaGo remains as before. The unfolding-refolding kinetics of the protein has also been studied by pressure-jump, in which it was revealed that both reactions are a two-state transition process with a relatively slow relaxation time of about 102 s.     

Are residues in a protein folding nucleus evolutionarily conserved?

Protein is the working molecule of the cell, and evolution is the hallmark of life. It is important to understand how protein folding and evolution influence each other. Several studies correlating experimental measurement of residue participation in folding nucleus and sequence conservation have reached different conclusions. These studies are based on assessment of sequence conservation at folding nucleus sites using entropy or relative entropy measurement derived from multiple sequence alignment. Here we report analysis of conservation of folding nucleus using an evolutionary model alternative to entropy-based approaches. We employ a continuous time Markov model of codon substitution to distinguish mutation fixed by evolution and mutation fixed by chance. This model takes into account bias in codon frequency, bias-favoring transition over transversion, as well as explicit phylogenetic information. We measure selection pressure using the ratio omega of synonymous versus non-synonymous substitution at individual residue site. The omega-values are estimated using the PAML method, a maximum-likelihood estimator. Our results show that there is little correlation between the extent of kinetic participation in protein folding nucleus as measured by experimental phi-value and selection pressure as measured by omega-value. In addition, two randomization tests failed to show that folding nucleus residues are significantly more conserved than the whole protein, or the median omega value of all residues in the protein. These results suggest that at the level of codon substitution, there is no indication that folding nucleus residues are significantly more conserved than other residues. We further reconstruct candidate ancestral residues of the folding nucleus and suggest possible test tube mutation studies for testing folding behavior of ancient folding nucleus.     

Identification of the hepatocyte Na+-dependent bile acid transport protein using monoclonal antibodies

Monoclonal antibodies have been utilized to characterize the hepatocyte Na+-dependent bile acid transport system. Sinusoidal plasma membrane proteins in the 49-54-kDa range, which are thought to be components of this transport system, based on photo-affinity labeling and reconstitution studies, have been partially purified by affinity chromatography and utilized as an immunogen for the production of a panel of monoclonal antibodies (mAb). One of these mAbs, 25A-3, recognized both a 49- and a 54-kDa protein as assessed by immunoprecipitation. In addition, it was shown to protect the bile acid transport system from inhibition by 4,4'-diisothiocyanostilbene-2,2'-disulfonic acid (DIDS) in a dose-dependent manner. DIDS covalently labeled membrane proteins of 49 and 54 kDa, and this process could be significantly inhibited when performed in the presence of mAb 25A-3. Furthermore, the DIDS-labeled membrane proteins were immunoprecipitated by 25A-3. These results establish that one of these membrane components is the bile acid carrier protein. Another mAb (25D-1) which immunoprecipitated only a 49-kDa protein was shown to block the protective effect of 25A-3 on DIDS inhibition of bile acid transport. In addition both antibodies effected each other's binding capacity to hepatocytes and reacted with the same 49-kDa protein as established by sequential immunoprecipitation. Binding studies indicated that there are approximately 3.3 X 10(6) 49-kDa transport molecules/hepatocyte. These results firmly establish that the 49-kDa protein is the Na+-dependent hepatocyte bile acid transporter.     

Toward minimalist models of larger proteins: a ubiquitin-like protein

Our recently developed off-lattice bead model capable of simulating protein structures with mixed alpha/beta content has been extended to model the folding of a ubiquitin-like protein and provides a means for examining the more complex kinetics involved in the folding of larger proteins. Using trajectories generated from constant-temperature Langevin dynamics simulations and sampling with the multiple multi-histogram method over five-order parameters, we are able to characterize the free energy landscape for folding and find evidence for folding through compact intermediates. Our model reproduces the observation that the C-terminus loop structure in ubiquitin is the last to fold in the folding process and most likely plays a spectator role in the folding kinetics. The possibility of a productive metastable intermediate along the folding pathway consisting of collapsed states with no secondary structure, and of intermediates or transition structures involving secondary structural elements occurring early in the sequence, is also supported by our model. The kinetics of folding remain multi-exponential below the folding temperature, with glass-like kinetics appearing at T/T(f) approximately 0.86. This new physicochemical model, designed to be predictive, helps validate the value of modeling protein folding at this level of detail for genomic-scale studies, and motivates further studies of other protein topologies and the impact of more complex energy functions, such as the addition of solvation forces.     

Influence of the GroE molecular chaperone machine on the in vitro refolding of Escherichia coli beta-galactosidase.

We have studied the effect of the components of the GroE molecular chaperone machine on the refolding of the Escherichia coli enzyme beta-galactosidase, a tetrameric protein whose 116-kDa promoters should not completely fit within the central cavity of the GroEL toroid. In the absence of other additives, GroEL formed a weak complex with chemically denatured beta-galactosidase, reduced its propensity to aggregate, and increased the recovery yields of active enzyme twofold without altering its folding pathway. When present together with the chaperonin, ATP--and to a lesser extent AMP-PNP--reduced the recovery yields and led to the resumption of aggregation. The use of the complete chaperonin system (GroEL, GroES, and ATP) eliminated the GroEL-mediated increase in recovery and folding proceeded less efficiently than in buffer alone. This unusual behavior can be explained in terms of a chaperonin "buffering" effect and the different affinities of GroE complexes for denatured beta-galactosidase.     

Transient conformational remodeling of folding proteins by GroES—individually and in concert with GroEL

The commonly accepted dogma of the bacterial GroE chaperonin system entails protein folding mediated by cycles of several ATP-dependent sequential steps where GroEL interacts with the folding client protein. In contrast, we herein report GroES-mediated dynamic remodeling (expansion and compression) of two different protein substrates during folding: the endogenous substrate MreB and carbonic anhydrase (HCAII), a well-characterized protein folding model. GroES was also found to influence GroEL binding induced unfolding and compression of the client protein underlining the synergistic activity of both chaperonins, even in the absence of ATP. This previously unidentified activity by GroES should have important implications for understanding the chaperonin mechanism and cellular stress response. Our findings necessitate a revision of the GroEL/ES mechanism.