Non-native interactions in protein folding intermediates: molecular dynamics simulations of hen lysozyme

Molecular dynamics simulations of protein denaturation can complement and extend experimental studies of protein folding by providing atomic-level structural information about conformational transitions and any conformational states along the unfolding pathway. Previous unfolding simulations of hen egg-white lysozyme have resulted in intermediate structures with an unfolded alpha-domain and a structured beta-domain, which is inconsistent with experiment. In contrast, the beta-domain unfolded first in the two simulations presented here leaving a structured alpha-domain. Following this, intermediate states were identified that differ with respect to the packing of the helices and the elements of non-native structure adopted. The non-native structure is critical for explaining many of the experimental observations. Overall, the pooled ensemble of these intermediates is in agreement with the experimental data for the major kinetic intermediate, suggesting that the kinetic intermediate may be made up of distinct, but rapidly interconverting, partially folded conformations distinguished primarily by differences in helix packing.     

Interaction of GroEL with a highly structured folding intermediate: iterative binding cycles do not involve unfolding

The GroE proteins are molecular chaperones involved in protein folding. The general mechanism by which they facilitate folding is still enigmatic. One of the central open questions is the conformation of the GroEL-bound nonnative protein. Several suggestions have been made concerning the folding stage at which a protein can interact with GroEL. Furthermore, the possibility exists that binding of the nonnative protein to GroEL results in its unfolding. We have addressed these issues that are basic for understanding the GroE-mediated folding cycle by using folding intermediates of an Fab antibody fragment as molecular probes to define the binding properties of GroEL. We show that, in addition to binding to an early folding intermediate, GroEL is able to recognize and interact with a late quaternary-structured folding intermediate (Dc) without measurably unfolding it. Thus, the prerequisite for binding is not a certain folding stage of a nonnative protein. In contrast, general surface properties of nonnative proteins seem to be crucial for binding. Furthermore, unfolding of a highly structured intermediate does not necessarily occur upon binding to GroEL. Folding of Dc in the presence of GroEL and ATP involves cycles of binding and release. Because in this system no off-pathway reactions or kinetic traps are involved, a quantitative analysis of the reactivation kinetics observed is possible. Our results indicate that the association reaction of Dc and GroEL in the presence of ATP is rather slow, whereas in the absence of ATP association is several orders of magnitude more efficient. Therefore, it seems that ATP functions by inhibiting reassociation rather than promoting release of the bound substrate.     

The burst-phase intermediate in the refolding of beta-lactoglobulin studied by stopped-flow circular dichroism and absorption spectroscopy

The kinetics of the guanidine hydrochloride-induced unfolding and refolding of bovine beta-lactoglobulin, a predominantly beta-sheet protein in the native state, have been studied by stopped-flow circular dichroism and absorption measurements at pH 3.2 and 4.5 degrees C. The refolding reaction was a complex process composed of different kinetic phases, while the unfolding was a single-phase reaction. Most notably, a burst-phase intermediate of refolding, which was formed during the dead time of stopped-flow measurements (approximately 18 ms), showed more intense ellipticity signals in the peptide region below 240 nm than the native state, yielding overshoot behavior in the refolding curves. We have investigated the spectral properties and structural stability of the burst-phase intermediate and also the structural properties in the unfolded state in 4.0 M guanidine hydrochloride of the protein and its disulfide-cleaved derivative. The main conclusions are: (1) the more intense ellipticity of the intermediate in the peptide region arises from formation of non-native alpha-helical structure in the intermediate, apparently suggesting that the folding of beta-lactoglobulin is not represented by a simple sequential mechanism. (2) The burst-phase intermediate has, however, a number of properties in common with the folding intermediates or with the molten globule states of other globular proteins whose folding reactions are known to be represented by the sequential model. These properties include: the presence of the secondary structure without the specific tertiary structure; formation of a hydrophobic core; broad unfolding transition of the intermediate; and rapidity of formation of the intermediate. The burst-phase intermediate of beta-lactoglobulin is thus classified as the same species as the molten globule state. (3) The circular dichroism spectra of beta-lactoglobulin and its disulfide-cleaved derivative in 4.0 M guanidine hydrochloride suggests the presence of the residual beta-structure in the unfolded state and the stabilization of the beta-structure by disulfide bonds. Thus; if this residual beta-structure is part of the native beta-structure and forms a folding initiation site, the folding reaction of beta-lactoglobulin may not necessarily be inconsistent with the sequential model. The non-native alpha-helices in the burst-phase intermediate may be formed in an immature part of the protein molecule because of the local alpha-helical propensity in this part.     

TCI temperamental scores in bulimia nervosa patients and normal women with and without diet experiences

The tailspike protein (TSP) of bacteriophage P22 is a homotrimeric multifunctional protein responsible for cell attachment and hydrolysis of the Salmonella typhimurium host cell receptor. Despite the folding of TSP involves the formation of thermolabile intermediates, the mature protein is extremely resistant to heat and detergent denaturation. We have analyzed the thermal resistance and unfolding pathway of two mutant, functional TSPs carrying end-terminal peptide fusions. Whereas the C-terminal fusion has minor effects on the TSP stability, the presence of a 23-mer foreign peptide at the N terminus (protein ATSP) results in a significant enhancement of the thermal resistance by retarding the first transition step of the unfolding process. At 65 degrees C and in 2% SDS, the unfolding rate constant for the transition from the native to the unfolding intermediate is 9.3 x 10(-4) s(-1) for ATSP versus 1.7 x 10(-3) s(-1) for wild-type TSP. On the other hand, the electrophoretic mobility of ATSP intermediates is greatly affected, proving structural modifications induced by the fused peptide. These results suggest a critical participation of the N-terminal domain in the unfolding kinetic barriers generated during the TSP denaturation pathway.     

Protein folding coupled to disulphide bond formation

Protein folding that is coupled to disulphide bond formation has many experimental advantages. In particular, the kinetic roles and importance of all the disulphide intermediates can be determined, usually unambiguously. This contrasts with other types of protein folding, where the roles of any intermediates detected are usually not established. Nevertheless, there is considerable confusion in the literature about even the best-characterized disulphide folding pathways. This article attempts to set the record straight.     

Folding mechanism of the alpha-subunit of tryptophan synthase, an alpha/beta barrel protein: global analysis highlights the interconversion of multiple native, intermediate, and unfolded forms through parallel channels

A variety of techniques have been used to investigate the urea-induced kinetic folding mechanism of the alpha-subunit of tryptophan synthase from Escherichia coli. A distinctive property of this 29 kDa alpha/beta barrel protein is the presence of two stable equilibrium intermediates, populated at approximately 3 and 5 M urea. The refolding process displays multiple kinetic phases whose lifetimes span the submillisecond to greater than 100 s time scale; unfolding studies yield two relaxation times on the order of 10-100 s. In an effort to understand the populations and structural properties of both the stable and transient intermediates, stopped-flow, manual-mixing, and equilibrium circular dichroism data were globally fit to various kinetic models. Refolding and unfolding experiments from various initial urea concentrations as well as forward and reverse double-jump experiments were critical for model discrimination. The simplest kinetic model that is consistent with all of the available data involves four slowly interconverting unfolded forms that collapse within 5 ms to a marginally stable intermediate with significant secondary structure. This early intermediate is an off-pathway species that must unfold to populate a set of four on-pathway intermediates that correspond to the 3 M urea equilibrium intermediate. Reequilibrations among these conformers act as rate-limiting steps in folding for a majority of the population. A fraction of the native conformation appears in less than 1 s at 25 degrees C, demonstrating that even large proteins can rapidly traverse a complex energy surface.     

Mutational analysis of the BPTI folding pathway: II. Effects of aromatic-->leucine substitutions on folding kinetics and thermodynamics

The rates of the individual steps in the disulfide-coupled folding and unfolding of eight BPTI variants, each containing a single aromatic to leucine amino acid replacement, were measured. From this analysis, the contributions of the four phenylalanine and four tyrosine residues to the stabilities of the native protein and the disulfide-bonded folding intermediates were determined. While the substitutions were found to destabilize the native protein by 2 to 7 kcal/mol, they had significantly smaller effects on the intermediates that represent the earlier stages of folding, even when the site of the substitution was located within the ordered regions of the intermediates. These results suggest that stabilizing interactions contribute less to conformational stability in the context of a partially folded intermediate than in a fully folded native protein, perhaps because of decreased cooperativity among the individual interactions. The kinetic analysis also provides new information about the transition states associated with the slowest steps in folding and unfolding, supporting previous suggestions that these transition states are extensively unfolded. Although the substitutions caused large changes in the distribution of folding intermediates and in the rates of some steps in the folding pathway, the kinetically-preferred pathway for all of the variants involved intramolecular disulfide rearrangements, as observed previously for the wild-type protein. These results suggest that the predominance of the rearrangement mechanism reflects conformational constraints present relatively early in the folding pathway.     

Refolding rate of stability-enhanced cytochrome c is independent of thermodynamic driving force

N52I iso-2 cytochrome c is a variant of yeast iso-2 cytochrome c in which asparagine substitutes for isoleucine 52 in an alpha helical segment composed of residues 49-56. The N52I substitution results in a significant increase in both stability and cooperativity of equilibrium unfolding, and acts as a "global suppressor" of destabilizing mutations. The equilibrium m-value for denaturant-induced unfolding of N52I iso-2 increases by 30%, a surprisingly large amount for a single residue substitution. The folding/unfolding kinetics for N52I iso-2 have been measured by stopped-flow mixing and by manual mixing, and are compared to the kinetics of folding/unfolding of wild-type protein, iso-2 cytochrome c. The results show that the observable folding rate and the guanidine hydrochloride dependence of the folding rate are the same for iso-2 and N52I iso-2, despite the greater thermodynamic stability of N52I iso-2. Thus, there is no linear free-energy relationship between mutation-induced changes in stability and observable refolding rates. However, for N52I iso-2 the unfolding rate is slower and the guanidine hydrochloride dependence of the unfolding rate is smaller than for iso-2. The differences in the denaturant dependence of the unfolding rates suggest that the N52I substitution decreases the change in the solvent accessible hydrophobic surface between the native state and the transition state. Two aspects of the results are inconsistent with a two-state folding/unfolding mechanism and imply the presence of folding intermediates: (1) observable refolding rate constants calculated from the two-state mechanism by combining equilibrium data and unfolding rate measurements deviate from the observed refolding rate constants; (2) kinetically unresolved signal changes ("burst phase") are observed for both N52I iso-2 and iso-2 refolding. The "burst phase" amplitude is larger for N52I iso-2 than for iso-2, suggesting that the intermediates formed during the "burst phase" are stabilized by the N52I substitution.     

SCL 90-R interpretation and brain tumour: a correction factor?

Two models are being considered for the mechanism of chaperonin-assisted protein folding in E. coli: (i) GroEL/GroES act primarily by enclosing substrate polypeptide in a folding cage in which aggregation is prevented during folding. (ii) GroEL mediates the repetitive unfolding of misfolded polypeptides, returning them onto a productive folding track. Both models are not mutually exclusive, but studies with the polypeptide-binding domain of GroEL have suggested that unfolding is the primary mechanism, enclosure being unnecessary. Here we investigate the capacity of the isolated apical polypeptide-binding domain to functionally replace the complete GroEL/GroES system. We show that the apical domain binds aggregation-sensitive polypeptides but cannot significantly assist their refolding in vitro and fails to replace the groEL gene or to complement defects of groEL mutants in vivo. A single-ring version of GroEL cannot substitute for GroEL. These results strongly support the view that sequestration of aggregation-prone intermediates in a folding cage is an important element of the chaperonin mechanism.     

Synergy between simulation and experiment in describing the energy landscape of protein folding

Experimental data from protein engineering studies and NMR spectroscopy have been used by theoreticians to develop algorithms for helix propensity and to benchmark computer simulations of folding pathways and energy landscapes. Molecular dynamic simulations of the unfolding of chymotrypsin inhibitor 2 (CI2) have provided detailed structural models of the transition state ensemble for unfolding/folding of the protein. We now have used the simulated transition state structures to design faster folding mutants of CI2. The models pinpoint a number of unfavorable local interactions at the carboxyl terminus of the single alpha-helix and in the protease-binding loop region of CI2. By removing these interactions or replacing them with stabilizing ones, we have increased the rate of folding of the protein up to 40-fold (tau = 0.4 ms). This correspondence, and other examples of agreement between experiment and theory in general, Phi-values and molecular dynamics simulations, in particular, suggest that significant progress has been made toward describing complete folding pathways at atomic resolution by combining experiment and simulation.     

Folding intermediates of wild-type and mutants of barnase. I. Use of phi-value analysis and m-values to probe the cooperative nature of the folding pre-equilibrium

It is difficult to determine whether transient folding intermediates have a cooperative (or first-order) folding transition without measuring their rates of formation directly. An intermediate I could be formed by a second-order transition from a denatured state D that is progressively changed into I as conditions are changed. We have not been able to monitor the rate of formation of the folding intermediate of barnase directly, but have analysed its reactivity and the equilibrium constant for its formation over a combination of wide ranges of temperature, concentration of denaturant and structural variation. Phase diagrams have been constructed for wild-type and 16 mutant proteins to map out the nature of the energy landscape of the denatured state. The free energy of unfolding of I, delta GD-I, changes with [urea] according to a highly cooperative transition. Further, mD-I (= delta delta GD-I/delta [urea]) for wild-type and several mutants is relatively insensitive to temperature, as would be expected for an intermediate that is formed cooperatively, rather than one that melts out according to a second-order transition. The phi-values for the formation of I change abruptly through the folding transitions rather than have the smooth changes expected for a second-order transition. There is a subset of mutants for which both mD-I and phi-value analysis indicate that a second intermediate becomes populated close to the melting temperatures of the native proteins. The folding intermediate of barnase is, thus, a relatively discrete and compact entity which is formed cooperatively.     

Structural characterisation of apoflavodoxin shows that the location of the stable nucleus differs among proteins with a flavodoxin-like topology

The structural characteristics of Azotobacter vinelandii apoflavodoxin II have been determined using multidimensional NMR spectroscopy. Apoflavodoxin has a stable, well-ordered core but its flavin binding region is flexible. The local stability of apoflavodoxin was probed using hydrogen/deuterium exchange measurements. The existence of an apoflavodoxin equilibrium folding intermediate is inferred from the non-coincidence of CD and fluorescence unfolding curves obtained for the guanidinium hydrochloride induced unfolding of apoflavodoxin. We suggest that the structured part of the putative intermediate is composed of the elements of secondary structure which have the slowest exchanging amide protons in the native protein. These elements are strands beta1, beta3, beta4 and beta5a and helices alpha4 and alpha5. We propose that it is a general feature of flavodoxins that the stable nucleus resides in the C-terminal part of these proteins. The results on flavodoxin are compared with those on two sequentially unrelated proteins sharing the flavodoxin-like fold: Che Y and cutinase. It is shown that the stable nucleus is found in different parts of the flavodoxin-like topology.     

Chimeras of human lysozyme and alpha-lactalbumin: an interesting tool for studying partially folded states during protein folding

Protein folding is an extremely active field of research where biology, chemistry, computer science and physics meet. Although the study of protein-folding intermediates in general and equilibrium intermediates in particular has grown considerably in recent years, many questions regarding the conformational state and the structural features of the various partially folded intermediate states remain unanswered. Performing kinetic measurements on proteins that have had their structures modified by site-directed mutagenesis, the so-called protein-engineering method, is an obvious way to gain fine structural information. In the present review, this method has been applied to a variety of proteins belonging to the lysozyme/alpha-lactalbumin family. Besides recombinants obtained by point mutations of individual critical residues, chimeric proteins in which whole structural elements (10-25 residues) from alpha-lactalbumin were inserted into a human lysozyme matrix are examined. The conformational properties of the equilibrium intermediate states are discussed together with the structural characterization of the partially unfolded states encountered in the kinetic folding pathway.     

Assembly characteristics of flagellar cap protein HAP2 of Salmonella: decamer and pentamer in the pH-sensitive equilibrium

The folding of a 98 residue protein, muscle acylphosphatase (AcP), has been studied using a variety of techniques including circular dichroism, fluorescence and NMR spectroscopy following transfer of chemically denatured protein into refolding conditions. A low-amplitude phase, detected in concurrence with the main kinetic phase, corresponds to the folding of a minor population (13%) of molecules with one or both proline residues in a cis conformation, as shown from the sensitivity of its rate to peptidyl prolyl isomerase. The major phase of folding has the same kinetic characteristics regardless of the technique employed to monitor it. The plots of the natural logarithms of folding and unfolding rate constants versus urea concentration are linear over a broad range of urea concentrations. Moreover, the initial state formed rapidly after the initiation of refolding is highly unstructured, having a similar circular dichroism, intrinsic fluorescence and NMR spectrum as the protein denatured at high concentrations of urea. All these results indicate that AcP folds in a two-state manner without the accumulation of intermediates. Despite this, the folding of the protein is extremely slow. The rate constant of the major phase of folding in water, kfH2O, is 0.23 s-1 at 28 degreesC and, at urea concentrations above 1 M, the folding process is slower than the cis-trans proline isomerisation step. The slow refolding of this protein is therefore not the consequence of populated intermediates that can act as kinetic traps, but arises from a large intrinsic barrier in the folding reaction.     

Acquisition of native beta-strand topology during the rapid collapse phase of protein folding

The 98 residue C-terminal domain of the cell-surface receptor protein CD2 (CD2.D1) has a beta-sandwich fold belonging to the immunoglobulin superfamily but lacking the usual disulfide bridges. Kinetic studies on the folding/unfolding of CD2.D1 reveal that folding proceeds through a rapidly formed intermediate state [Parker, M. J., & Clarke, A. R. (1997) Biochemistry 36, 5786-5794]. To characterize the structural properties of this intermediate we have performed a series of amide hydrogen exchange studies using the pH competition method, in which folding and exchange are initiated simultaneously. The complex beta-sheet topology of this molecule makes it an ideal object for examining the acquisition of backbone hydrogen bonds made between sequence-local and sequence-distant segments of the chain during folding. The pattern of protected amides in the intermediate reveal that the essential features of the beta-sheet topology of CD2.D1 are defined early in the folding pathway, before the development of intimate side chain interactions characteristic of the native state. The results are discussed in light of current issues concerning the mechanistic relevance of kinetic protein folding intermediates.     

The folding kinetics and thermodynamics of the Fyn-SH3 domain

The equilibrium unfolding and the kinetic folding and unfolding of the 67 residue Fyn-SH3 domain have been investigated. Equilibrium unfolding experiments indicate that, despite the lack of both disulfide bonds and prosthetic groups, Fyn-SH3 is relatively stable with a free energy of folding of -6.0 +/- 0.6 kcal mol-1 at 20 degrees C. Kinetic experiments indicate that the domain refolds in a rapid two-state manner without significant population of intermediates (k = 94.3 s-1 in H2O at 20 degrees C). Despite the presence of two proline residues, the refolding of the domain is monophasic, and no significant proline isomerization-like refolding phase is observed. This can be attributed to an extremely low level of the incorrect (cis) isomer of the structurally important Pro134 residue in the protein denatured in 8 M guanidine hydrochloride. Analysis of the temperature and guanidine hydrochloride dependence of the folding rate suggests that the folding transition state of this protein is relatively well organized. A comparison with the refolding kinetics and thermodynamics of other homologous SH3 domains indicates that these exhibit an equivalent degree of transition state organization. This potentially arises from conservation of key features of the transition state conformation despite sometimes relatively low overall sequence identity. Such a comparison further suggests that relative thermodynamic stability is an important factor in determining the relative folding rates of natural proteins with a common fold, but that specific details of the amino acid sequence can also play a significant role in individual cases.     

Interactions in nonnative and truncated forms of staphylococcal nuclease as indicated by mutational free energy changes

Several mixed disulfide variants of staphylococcal nuclease have been produced by disulfide bond formation between nuclease V23C and methane, ethane, 1-propane, 1-n-butane, and 1-n-pentane thiols. Although CD spectroscopy shows that the native state is largely unperturbed, the stability toward urea-induced unfolding is highly dependent on the nature of the group at this position, with the methyl disulfide protein being the most stable. The variant produced by modification with iodoacetic acid, however, gives a CD spectrum indicative of an unfolded polypeptide. Thiol-disulfide exchange equilibrium constants between nuclease V23C and 2-hydroxyethyl disulfide have been measured as a function of urea concentration. Because thiol-disulfide exchange and unfolding are thermodynamically linked, the effects of a mutation (disulfide exchange) can be partitioned between various conformational states. In the case of unmodified V23C and the 2-hydroxyethyl protein mixed disulfide, significant effects in the nonnative states of nuclease are observed. Truncated forms of staphylococcal nuclease are thought to be partially folded and may be good models for early folding intermediates. We have characterized a truncated form of nuclease comprised of residues 1-135 with a V23C mutation after chemical modification of the cysteine residue. High-resolution size-exclusion chromatography indicates that modification brings about significant changes in the Stokes radius of the protein, and CD spectroscopy indicates considerable differences in the amount of secondary structure present. Measurement of the disulfide exchange equilibrium constant between this truncated protein and 2-hydroxyethyl disulfide indicate significant interactions between position 23 and the rest of the protein when the urea concentration is lower than 1.5 M.(ABSTRACT TRUNCATED AT 250 WORDS)     

Fast and one-step folding of closely and distantly related homologous proteins of a four-helix bundle family

Bovine acyl-coenzyme A binding protein is a four-helix bundle protein belonging to a group of homologous eukaryote proteins that binds medium and long-chain acyl-coenzyme A esters with a very high affinity. The three-dimensional structure of both the free and the ligated protein together with the folding kinetics have been described in detail for the bovine protein and with four new sequences reported here, a total of 16 closely related sequences ranging from yeasts and plants to human are known. The kinetics of folding and unfolding in different concentrations of guanidine hydrochloride together with equilibrium unfolding have been measured for bovine, rat and yeast acyl-coenzyme A binding protein. The bovine and rat sequences are closely related whereas the yeast is more distantly related to these. In addition to the three natural variants, kinetics of a bovine mutant protein, Tyr31 --> Asn, have been studied. Both the folding and unfolding rates in water of the yeast protein are 15 times faster than those of bovine. The folding rates in water of the two mammalian forms, rat and bovine, are similar, though still significantly different. A faster unfolding rate both for rat and the bovine mutant protein results from a lower stability of the native states of these. These hydrophobic regions, mini cores, have been identified in the three-dimensional structure of the bovine protein and found to be formed primarily by residues that have been conserved throughout the entire eukaryote evolution from yeasts to both plants and mammals as seen in the sample of 16 sequences. The conserved residues are found to stabilize helix-helix interactions and serve specific functional purposes for ligand binding. The fast one-step folding mechanism of ACBP has been shown to be a feature that seems to be maintained throughout evolution despite numerous differences in sequence and even dramatic differences in folding kinetics and protein stability. The protein study raises the question to what extent does the conserved hydrophobic residues provide a scaffold for an efficient one-step folding mechanism.     

Probing the partly folded states of proteins by limited proteolysis

The folding of a polypeptide chain of a relatively large globular protein into its unique three-dimensional and functionally active structure occurs via folding intermediates. These partly folded states of proteins are difficult to characterize, because they are usually short lived or exist as a distribution of possible conformers. A variety of experimental techniques and approaches have been utilized in recent years in numerous laboratories for characterizing folding intermediates that occur at equilibrium, including spectroscopic techniques, solution X-ray scattering, calorimetry and gel filtration chromatography, as well as genetic methods and theoretical calculations. In this review, we focus on the use of proteolytic enzymes as probes of the structure and dynamics of folding intermediates and we show that this simple biochemical technique can provide useful information, complementing that obtained by other commonly used techniques and approaches. The key result of the proteolysis experiments is that partly folded states (molten globules) of proteins can be sufficiently rigid to prevent extensive proteolysis and appear to maintain significant native-like structure.