Molecular simulations of cotranslational protein folding: fragment stabilities, folding cooperativity, and trapping in the ribosome

Although molecular simulation methods have yielded valuable insights into mechanistic aspects of protein refolding in vitro, they have up to now not been used to model the folding of proteins as they are actually synthesized by the ribosome. To address this issue, we report here simulation studies of three model proteins: chymotrypsin inhibitor 2 (CI2), barnase, and Semliki forest virus protein (SFVP), and directly compare their folding during ribosome-mediated synthesis with their refolding from random, denatured conformations. To calibrate the methodology, simulations are first compared with in vitro data on the folding stabilities of N-terminal fragments of CI2 and barnase; the simulations reproduce the fact that both the stability and thermal folding cooperativity increase as fragments increase in length. Coupled simulations of synthesis and folding for the same two proteins are then described, showing that both fold essentially post-translationally, with mechanisms effectively identical to those for refolding. In both cases, confinement of the nascent polypeptide chain within the ribosome tunnel does not appear to promote significant formation of native structure during synthesis; there are however clear indications that the formation of structure within the nascent chain is sensitive to location within the ribosome tunnel, being subject to both gain and loss as the chain lengthens. Interestingly, simulations in which CI2 is artificially stabilized show a pronounced tendency to become trapped within the tunnel in partially folded conformations: non-cooperative folding, therefore, appears in the simulations to exert a detrimental effect on the rate at which fully folded conformations are formed. Finally, simulations of the two-domain protease module of SFVP, which experimentally folds cotranslationally, indicate that for multi-domain proteins, ribosome-mediated folding may follow different pathways from those taken during refolding. Taken together, these studies provide a first step toward developing more realistic methods for simulating protein folding as it occurs in vivo.     

Fluorescence energy transfer indicates similar transient and equilibrium intermediates in staphylococcal nuclease folding

Fluorescence resonance energy transfer (FRET) is one of the few methods available to measure the rate at which a folding protein collapses. Using staphylococcal nuclease in which a cysteine residue was engineered in place of Lys64, permitted FRET measurements of the distance between the donor tryptophan 140 and 5-[[2-[(iodoacetyl)-amino]ethyl]amino]naphthalene-1-sulfonic acid-labeled Cys64. These measurements were undertaken on both equilibrium partially folded intermediates at low pH (A states), as well as transient intermediates during stopped-flow refolding. The results indicate that there is an initial collapse of the protein in the deadtime of the stopped-flow instrument, corresponding to a regain of approximately 60% of the native signal, followed by three slower transients. This is in contrast to circular dichroism measurements which show only 20-25% regain of the native secondary structure in the burst phase. Thus hydrophobic collapse precedes the formation of substantial secondary structure. The first two detected transient intermediate species have FRET properties essentially identical with those of the previously characterized equilibrium A state intermediates, suggesting similar structures between the equilibrium and transient intermediates.The effects of anions on the folding of acid-unfolded staphylococcal nuclease, and urea on the unfolding of the resulting A states, indicates that in folding the protein becomes compact prior to formation of major secondary structure, whereas in unfolding the protein expands prior to major loss of secondary structure. Comparison of the kinetics of refolding of staphylococcal nuclease, monitored by FRET, and for a proline-free variant, indicate that folding occurs via two partially folded intermediates leading to a native-like species with one (or more) proline residues in a non-native conformation. For the A states an excellent correlation between compactness measured by FRET, and compactness determined from small-angle X-ray scattering, was observed. Further, a linear relationship between compactness and free energy of unfolding was noted. Formation of soluble aggregates of the A states led to dramatic enhancement of the FRET, consistent with intermolecular fluorescence energy transfer.     

Aggregation events occur prior to stable intermediate formation during refolding of interleukin 1beta

A point mutation, lysine 97 --> isoleucine (K97I), in a surface loop in the beta-sheet protein interleukin 1beta (IL-1beta), exhibits increased levels of inclusion body (IB) formation relative to the wild-type protein (WT) when expressed in Escherichia coli. Despite the common observation that less stable proteins are often found in IBs, K97I is more stable than WT. We examined the folding pathway of the mutant and wild-type proteins at pH 6.5 and 25 degrees C with manual-mixing and stopped-flow optical spectroscopy to determine whether changes in the properties of transiently populated species in vitro correlate with the observation of increased aggregation in vivo. The refolding reactions of the WT and K97I proteins are both described by three exponential processes. Two exponential processes characterize fast events (0.1-1.0 s) in folding while the third exponential process correlates with a slow (70 s) single pathway to and from the native state. The K97I replacement affects the earlier steps in the refolding pathway. Aggregation, absent in the WT refolding reaction, occurs in K97I above a critical protein concentration of 18 microM. This observation is consistent with an initial nucleation step mediating protein aggregation. Stopped-flow kinetic studies of the K97I aggregation process demonstrate that K97I aggregates most rapidly during the earliest refolding times, when unfolded protein conformers remain highly populated and the concentration of folding intermediates is low. Folding and aggregation studies together support a model in which the formation of stable folding intermediates afford protection against further K97I aggregation.     

Folding of intracellular retinol and retinoic acid binding proteins

The folding mechanisms of cellular retinol binding protein II (CRBP II), cellular retinoic acid binding protein I (CRABP I), and cellular retinoic acid binding protein II (CRABP II) were examined. These beta-sheet proteins have very similar structures and higher sequence homologies than most proteins in this diverse family. They have similar stabilities and show completely reversible folding at equilibrium with urea as a denaturant. The unfolding kinetics of these proteins were monitored during folding and unfolding by circular dichroism (CD) and fluorescence. During unfolding, CRABP II showed no intermediates, CRABP I had an intermediate with nativelike secondary structure, and CRBP II had an intermediate that lacked secondary structure. The refolding kinetics of these proteins were more similar. Each protein showed a burst-phase change in intensity by both CD and fluorescence, followed by a single observed phase by both CD and fluorescence and one or two additional refolding phases by fluorescence. The fluorescence spectral properties of the intermediate states were similar and suggested a gradual increase in the amount of native tertiary structure present for each step in a sequential path. However, the rates of folding differed by as much as 3 orders of magnitude and were slower than those expected from the contact order and topology of these proteins. As such, proteins with the same final structure may not follow the same route to the native state.     

蛋白质的氧化重折叠

经过近几十年来广泛而深入的研究,蛋白质氧化重折叠的机制已得到相当详细的阐明。1在已研究过的蛋白质中,大多数蛋白质都是沿着多途径而非单一、特定的途径进行氧化重折叠,这与折叠能量景观学说是一致的。2正是氨基酸残基间的天然相互作用而不是非天然的相互作用控制蛋白质的折叠过程。这一结论与含非天然二硫键的折叠中间体在牛胰蛋白酶抑制剂（BPTI）折叠中所起的重要作用并非相互排斥,因为后者仅仅是进行链内二硫键重排的化学反应所必需,与控制肽链折叠无直接关系。3根据对BPTI的研究,二硫键曾被认为仅仅具有稳定蛋白质天然结构的作用,既不决定折叠途径也不决定其三维构象。这一观点不适用于其它蛋白质。对凝乳酶原的研究表明,天然二硫键的形成是恢复天然构象的前提。天然二硫键的形成与肽键的正确折叠相辅相成,更具有普遍意义。4在氧化重折叠的早期,二硫键的形成基本上是一个随机过程,随着肽链的折叠二硫键的形成越来越受折叠中间体构象的限制。提高重组蛋白质的复性产率是生物技术领域中的一个巨大的挑战。除了分子聚集外,在折叠过程中所形成的二硫键错配分子是导致低复性率的另一个主要原因。氧化重折叠机制的阐明为解决此问题提供了有益的启示。如上所述,在折叠的后期,二硫键的形成决定于折叠中间体的构象,类天然、有柔性的结构有利于天然二硫键形成和正确折叠,具有这类结构的分子为有效的折叠中间体,最终都能转变为天然产物;而无效折叠中间体往往具有稳定的结构,使巯基、二硫键内埋妨碍二硫键重排,并因能垒的障碍不利于进一步折叠。因此,降低无效折叠中间体的稳定性使之转变为有效折叠中间体是提高含二硫键蛋白质复性率的一条基本原则,实验证明,碱性pH、低温、降低蛋白质稳定性的试剂、蛋白质二硫键异构酶、改变蛋白质一级结构是实现这一原则的有效手段。此外,这里还就氧化重折叠的基础和应用研究的前景进行了讨论。     

Hydrogen exchange in BPTI variants that do not share a common disulfide bond.

Bovine pancreatic trypsin inhibitor (BPTI) is stabilized by 3 disulfide bonds, between cysteines 30-51, 5-55, and 14-38. To better understand the influence of disulfide bonds on local protein structure and dynamics, we have measured amide proton exchange rates in 2 folded variants of BPTI, [5-55]Ala and [30-51; 14-38]V5A55, which share no common disulfide bonds. These proteins resemble disulfide-bonded intermediates that accumulate in the BPTI folding pathway. Essentially the same amide hydrogens are protected from exchange in both of the BPTI variants studied here as in native BPTI, demonstrating that the variants adopt fully folded, native-like structures in solution. However, the most highly protected amide protons in each variant differ, and are contained within the sequences of previously studied peptide models of related BPTI folding intermediates containing either the 5-55 or the 30-51 disulfide bond.     

An early intermediate in the folding of ribonuclease A is protected against cleavage by pepsin

Folding of bovine pancreatic ribonuclease A (RNase A) is a sequential process which involves the formation of well-populated structural intermediates under suitable conditions. Two intermediates have been detected on the major slow-refolding pathway of RNase A: a late intermediate (IN) which already resembles the native protein in a number of properties and a rapidly formed early intermediate (I1) which shows extensive hydrogen-bonded secondary structure. Here competition experiments between refolding and proteolytic cleavage of the peptide chain are described which yield information about the decrease in accessibility of particular proteolytic cleavage sites during the folding process. Results obtained with pepsin as a proteolytic probe of folding indicate that the primary cleavage site for pepsin, Phe-120-Asp-121, becomes inaccessible early in the course of refolding, if folding is carried out under conditions which effectively stabilize the native state. Under marginally stable conditions, folding is very slow, and protection against peptic cleavage is not detectable prior to the final formation of native protein. The comparison with amide proton exchange experiments suggests that the protection against peptic cleavage occurs during the formation and/or stabilization of hydrogen-bonded secondary structure in the early intermediate (I1). We conclude that the carboxy-terminal region of the RNase peptide chain, which is known to be important for the stability of the folded protein, may also be relevant for early steps of refolding.     

Peptide models of four possible insulin folding intermediates with two disulfides

The single-chain insulin (PIP) can spontaneously fold into native structure through preferred kinetic intermediates. During refolding, pairing of the first disulfide A20-B19 is highly specific, whereas pairing of the second disulfide is likely random because two two-disulfide intermediates have been trapped. To get more details of pairing property of the second disulfide, four model peptides of possible folding intermediates with two disulfides were prepared by protein engineering, and their properties were analyzed. The four model peptides were named [A20-B19, A7-B7]PIP, [A20-B19, A6-B7]PIP, [A20-B19, A6-A11]PIP, and [A20-B19, A7-A11]PIP according to their remaining disulfides. The four model peptides all adopt partially folded structure with moderate conformational differences. In redox buffer, the disulfides of the model peptides are more easily reduced than those of the wild-type PIP. During in vitro refolding, the reduced model peptides share similar relative folding rates but different folding yields: The refolding efficiency of the reduced [A20-B19, A7-A11]PIP is about threefold lower than that of the other three peptides. The present results indicate that the folding intermediates corresponding to the present model peptides all adopt partially folded conformation, and can be formed during PIP refolding, but the chance of forming the intermediate with disulfide [A20-B19, A7-A11] is much lower than that of forming the other three intermediates.     

The kinetics of bovine growth hormone folding are consistent with a framework model

The framework model of protein folding requires the hydrogen-bonded secondary structure to be formed early in folding (i.e. the formation of secondary structure precedes the tertiary structure) (Kim, P. S., and Baldwin, R. L. (1982) Annu. Rev. Biochem. 51, 459-489). To test the framework model directly the kinetics of bovine growth hormone (bGH) folding were compared utilizing two methods of detection, one that measures the secondary structure (far ultraviolet circular dichroism) and another that measures the tertiary structure (near ultraviolet absorbance). The results demonstrate that, under identical folding conditions, the kinetics observed by far ultraviolet circular dichroism are faster than those observed by ultraviolet absorption. The faster kinetics observed by circular dichroism indicate the existence of a helix-containing intermediate which is consistent with the framework model. The effect of protein concentration and denaturant concentration on the kinetics of refolding were studied. The rate of refolding measured by absorbance and circular dichroism was dependent on protein concentration. The protein concentration dependence on refolding is due to the transient formation of an associated intermediate. The concentration dependence of folding is taken as evidence that folding is a sequential process with partially folded monomers responsible for the observed association effect. At dilute protein concentrations the refolding can be studied independent of the association phenomena. The growth hormones utilized in this study were derived from Escherichia coli through recombinant DNA technology and from bovine pituitaries. The pituitary-derived bGH has been shown to be heterogeneous at the NH2 terminus (Lorenson, M. F., and Ellis, S. (1975) Endocrinology 96, 833-838), whereas the recombinant DNA-derived bGH contains a single NH2 terminus. No differences in the folding kinetics between the recombinant DNA and pituitary derived-bGH were observed. It is concluded that the heterogeneity of the NH2 terminus of growth hormone obtained from bovine pituitaries does not affect the observed in vitro folding kinetics.     

The Folding Trajectory of RNase H Is Dominated by Its Topology and Not Local Stability: A Protein Engineering Study of Variants that Fold via Two-State and Three-State Mechanisms

Proteins can sample a variety of partially folded conformations during the transition between the unfolded and native states. Some proteins never significantly populate these high-energy states and fold by an apparently two-state process. However, many proteins populate detectable, partially folded forms during the folding process. The role of such intermediates is a matter of considerable debate. A single amino acid change can convert Escherichia coli ribonuclease H from a three-state folder that populates a kinetic intermediate to one that folds in an apparent two-state fashion. We have compared the folding trajectories of the three-state RNase H and the two-state RNase H, proteins with the same native-state topology but altered regional stability, using a protein engineering approach. Our data suggest that both versions of RNase H fold through a similar trajectory with similar high-energy conformations. Mutations in the core and the periphery of the protein affect similar aspects of folding for both variants, suggesting a common trajectory with folding of the core region followed by the folding of the periphery. Our results suggest that formation of specific partially folded conformations may be a general feature of protein folding that can promote, rather than hinder, efficient folding.     

Mechanistic studies of the folding of human lysozyme and the origin of amyloidogenic behavior in its disease-related variants.

The unfolding and refolding properties of human lysozyme and two amyloidogenic variants (Ile56Thr and Asp67His) have been studied by stopped-flow fluorescence and hydrogen exchange pulse labeling coupled with mass spectrometry. The unfolding of each protein in 5.4 M guanidine hydrochloride (GuHCl) is well described as a two-state process, but the rates of unfolding of the Ile56Thr variant and the Asp67His variant in 5.4 M GuHCl are ca. 30 and 160 times greater, respectively, than that of the wild type. The refolding of all three proteins in 0.54 M GuHCl at pH 5.0 proceeds through persistent intermediates, revealed by multistep kinetics in fluorescence experiments and by the detection of well-defined populations in quenched-flow hydrogen exchange experiments. These findings are consistent with a predominant mechanism for refolding of human lysozyme in which one of the structural domains (the alpha-domain) is formed in two distinct steps and is followed by the folding of the other domain (the beta-domain) prior to the assembly of the two domains to form the native structure. The refolding kinetics of the Asp67His variant are closely similar to those of the wild-type protein, consistent with the location of this mutation in an outer loop of the beta-domain which gains native structure only toward the end of the refolding process. By contrast, the Ile56Thr mutation is located at the base of the beta-domain and is involved in the domain interface. The refolding of the alpha-domain is unaffected by this substitution, but the latter has the effect of dramatically slowing the folding of the beta-domain and the final assembly of the native structure. These studies suggest that the amyloidogenic nature of the lysozyme variants arises from a decrease in the stability of the native fold relative to partially folded intermediates. The origin of this instability is different in the two variants, being caused in one case primarily by a reduction in the folding rate and in the other by an increase in the unfolding rate. In both cases this results in a low population of soluble partially folded species that can aggregate in a slow and controlled manner to form amyloid fibrils.     

A salt-induced kinetic intermediate is on a new parallel pathway of lysozyme folding.

Lysozyme folds through two competing pathways. A fast pathway leads directly from a collapsed state to the native protein, whereas folding on a slow pathway proceeds through a partially folded intermediate (I(1)). At NaCl concentrations above 100 mM, a second transient intermediate (I(2)) is induced as judged by the appearance of an additional apparent rate constant in the refolding kinetics. Monitoring the time course of native molecules and of both intermediates shows that the NaCl-induced state (I(2)) is located on neither of the two folding pathways observed at low-salt concentrations. These results suggest that I(2) is a metastable high-energy intermediate at low-ionic strength and is located on a third folding pathway. The folding landscape of lysozyme seems to be complex with several high-energy intermediates located on parallel folding routes. However, the experiments show no evidence for partially folded states on the fast direct pathway.     

Estimation of parameters in a multi-affinity-state model for haemoglobin from oxygen binding data in whole blood and in concentrated haemoglobin solutions.

The reduced fragment of pancreatic trypsin inhibitor lacking the six C-terminal residues, which is produced by cyanogen bromide cleavage, formed a seemingly random mixture of disulphide bonds under refolding conditions where normal pancreatic trypsin inhibitor refolds correctly and quantitatively. This illustrates the importance of the C-terminal residues in folding of the normal protein, the uniqueness of the normal folded conformation, and the apparently central role in protein folding of long-range interactions between residues distant in the primary structure.The intact polypeptide chain of reduced pancreatic trypsin inhibitor in which the methionine residue normally at position 52 had been converted to homoserine refolded slightly less readily than the normal reduced compound. This was observed to be due to an altered spectrum of single-disulphide intermediates: the normally predominant intermediate with the 30–51 disulphide bond was less stable by about 0.8 kcal/mol relative to the other normal single-disulphide intermediates. The other steps in refolding appeared to be normal, although the refolded protein was observed to be susceptible to an unexplained reaction with iodoacetate.     

Slow formation of aggregation‐resistant β‐sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding‐related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large β‐helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all‐β‐sheet protein allows detailed analysis of the formation of β‐sheet structure in larger proteins. Using a combination of fluorescence and far‐UV circular dichroism spectroscopy, we show that the pertactin β‐helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and β‐sheet‐rich topology, pertactin refolding is reversible and not complicated by off‐pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate‐limiting step. Furthermore, site‐specific labeling experiments indicate that the β‐helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, β‐sheet‐rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, β‐sheet‐rich refolding intermediates. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Energetic coupling between native-state prolyl isomerization and conformational protein folding

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

Formation of a native-like subdomain in a partially folded intermediate of bovine pancreatic trypsin inhibitor.

In the folding of bovine pancreatic trypsin inhibitor (BPTI), the single-disulfide intermediate [30-51] plays a key role. We have investigated a recombinant analog of [30-51] using a 2-dimensional nuclear magnetic resonance (2D-NMR). This recombinant analog, named [30-51]Ala, contains a disulfide bond between Cys-30 and Cys-51, but contains alanine in place of the other cysteines in BPTI to prevent the formation of other intermediates. By 2D-NMR, [30-51]Ala consists of 2 regions-one folded and one predominantly unfolded. The folded region resembles a previously characterized peptide model of [30-51], named P alpha P beta, that contains a native-like subdomain with tertiary packing. The unfolded region includes the first 14 N-terminal residues of [30-51] and is as unfolded as an isolated peptide containing these residues. Using protein dissection, we demonstrate that the folded and unfolded regions of [30-51]Ala are structurally independent. The partially folded structure of [30-51]Ala explains many of the properties of authentic [30-51] in the folding pathway of BPTI. Moreover, direct structural characterization of [30-51]Ala has revealed that a crucial step in the folding pathway of BPTI coincides with the formation of a native-like subdomain, supporting models for protein folding that emphasize the formation of cooperatively folded subdomains.     

The group II chaperonin Mm-Cpn binds and refolds human γD crystallin

Chaperonins assist in the folding of nascent and misfolded proteins, though the mechanism of folding within the lumen of the chaperonin remains poorly understood. The archeal chaperonin from Methanococcus marapaludis, Mm-Cpn, shares the eightfold double barrel structure with other group II chaperonins, including the eukaryotic TRiC/CCT, required for actin and tubulin folding. However, Mm-Cpn is composed of a single species subunit, similar to group I chaperonin GroEL, rather than the eight subunit species needed for TRiC/CCT. Features of the β-sheet fold have been identified as sites of recognition by group II chaperonins. The crystallins, the major components of the vertebrate eye lens, are β-sheet proteins with two homologous Greek key domains. During refolding in vitro a partially folded intermediate is populated, and partitions between productive folding and off-pathway aggregation. We report here that in the presence of physiological concentrations of ATP, Mm-Cpn suppressed the aggregation of HγD-Crys by binding the partially folded intermediate. The complex was sufficiently stable to permit recovery by size exclusion chromatography. In the presence of ATP, Mm-Cpn promoted the refolding of the HγD-Crys intermediates to the native state. The ability of Mm-Cpn to bind and refold a human β-sheet protein suggests that Mm-Cpn may be useful as a simplified model for the substrate recognition mechanism of TRiC/CCT.     

Folding and unfolding of gammaTIM monomers and dimers

Kinetic simulations of the folding and unfolding of triosephosphate isomerase (TIM) from yeast were conducted using a single monomer gammaTIM polypeptide chain that folds as a monomer and two gammaTIM chains that fold to the native dimer structure. The basic protein model used was a minimalist Gō model using the native structure to determine attractive energies in the protein chain. For each simulation type--monomer unfolding, monomer refolding, dimer unfolding, and dimer refolding--thirty simulations were conducted, successfully capturing each reaction in full. Analysis of the simulations demonstrates four main conclusions. First, all four simulation types have a similar "folding order", i.e., they have similar structures in intermediate stages of folding between the unfolded and folded state. Second, despite this similarity, different intermediate stages are more or less populated in the four different simulations, with 1), no intermediates populated in monomer unfolding; 2), two intermediates populated with beta(2)-beta(4) and beta(1)-beta(5) regions folded in monomer refolding; 3), two intermediates populated with beta(2)-beta(3) and beta(2)-beta(4) regions folded in dimer unfolding; and 4), two intermediates populated with beta(1)-beta(5) and beta(1)-beta(5) + beta(6) + beta(7) + beta(8) regions folded in dimer refolding. Third, simulations demonstrate that dimer binding and unbinding can occur early in the folding process before complete monomer-chain folding. Fourth, excellent agreement is found between the simulations and MPAX (misincorporation proton alkyl exchange) experiments. In total, this agreement demonstrates that the computational Gō model is accurate for gammaTIM and that the energy landscape of gammaTIM appears funneled to the native state.