Folding pathway of the b1 domain of protein G explored by multiscale modeling

The understanding of the folding mechanisms of single-domain proteins is an essential step in the understanding of protein folding in general. Recently, we developed a mesoscopic CA-CB side-chain protein model, which was successfully applied in protein structure prediction, studies of protein thermodynamics, and modeling of protein complexes. In this research, this model is employed in a detailed characterization of the folding process of a simple globular protein, the B1 domain of IgG-binding protein G (GB1). There is a vast body of experimental facts and theoretical findings for this protein. Performing unbiased, ab initio simulations, we demonstrated that the GB1 folding proceeds via the formation of an extended folding nucleus, followed by slow structure fine-tuning. Remarkably, a subset of native interactions drives the folding from the very beginning. The emerging comprehensive picture of GB1 folding perfectly matches and extends the previous experimental and theoretical studies.     

Fuzzy-oil-drop hydrophobic force field--a model to represent late-stage folding (in silico) of lysozyme

A model of hydrophobic collapse (in silico), which is generally considered to be the driving force for protein folding, is presented in this work. The model introduces the external field in the form of a fuzzy-oil-drop assumed to represent the environment. The drop is expressed in the form of a three-dimensional Gauss function. The usual probability value is assumed to represent the hydrophobicity distribution in the three-dimensional space of the virtual environment. The differences between this idealized hydrophobicity distribution and the one represented by the folded polypeptide chain is the parameter to be minimized in the structure optimization procedure. The size of fuzzy-oil-drop is critical for the folding process. A strong correlation between protein length and the dimension of the native and early-stage molecular form was found on the basis of single-domain proteins analysis. A previously presented early-stage folding (in silico) model was used to create the starting structure for the procedure of late-stage folding of lysozyme. The results of simulation were found to be promising, although additional improvements for the formation of beta-structure and disulfide bonds as well as the participation of natural ligand in folding process seem to be necessary.     

Simulation of the thermal denaturation of hen egg white lysozyme: trapping the molten globule state.

In the study of protein folding, much attention has focused on the characterization of folding intermediates. We report here molecular dynamics simulations in which the initial stages of the thermal denaturation of hen egg white lysozyme in aqueous solution are examined in detail. It is found that lysozyme unfolds in a two-stage process with the initial formation a quasi-stable state in which significant rearrangement of the secondary structure takes place. No evidence for distinct folding domains was found. The simulations suggest that the formation of well-defined secondary structure occurs after the initial collapse of the peptide chain and thus tend against the framework model of protein folding.     

Membrane protein folding: beyond the two stage model

The folding of alpha-helical membrane proteins has previously been described using the two stage model, in which the membrane insertion of independently stable alpha-helices is followed by their mutual interactions within the membrane to give higher order folding and oligomerization. Given recent advances in our understanding of membrane protein structure it has become apparent that in some cases the model may not fully represent the folding process. Here we present a three stage model which gives considerations to ligand binding, folding of extramembranous loops, insertion of peripheral domains and the formation of quaternary structure.     

Is the unfolded state the Rosetta Stone of the protein folding problem?

Solving the protein folding problem is one of the most challenging tasks in the post genomic era. Identification of folding-initiation sites is very important in order to understand the protein folding mechanism. Detection of residual structure in unfolded proteins can yield important clues to the initiation sites in protein folding. A substantial number of studied proteins possess residual structure in hydrophobic regions clustered together in the protein core. These stable structures can work as seeds in the folding process. In addition, local preferences for secondary structure in the form of turns for beta-sheet initiation and helical turns for alpha-helix formation can guide the folding reaction. In this respect the unfolded states, studied at increasing structural resolution, can be the Rosetta Stone of the protein folding problem.     

Fuzzy-Oil-Drop Hydrophobic Force Field—A Model to Represent Late-stage Folding (In Silico) of Lysozyme

Abstract

A model of hydrophobic collapse (in silico), which is generally considered to be the driving force for protein folding, is presented in this work. The model introduces the external field in the form of a fuzzy-oil-drop assumed to represent the environment. The drop is expressed in the form of a three-dimensional Gauss function. The usual probability value is assumed to represent the hydrophobicity distribution in the three-dimensional space of the virtual environment. The differences between this idealized hydrophobicity distribution and the one represented by the folded polypeptide chain is the parameter to be minimized in the structure optimization procedure. The size of fuzzy-oil-drop is critical for the folding process. A strong correlation between protein length and the dimension of the native and early-stage molecular form was found on the basis of single-domain proteins analysis. A previously presented early-stage folding (in silico) model was used to create the starting structure for the procedure of late-stage folding of lysozyme. The results of simulation were found to be promising, although additional improvements for the formation of β-structure and disulfide bonds as well as the participation of natural ligand in folding process seem to be necessary.     

Is protein folding hierarchic? II. Folding intermediates and transition states

The folding reactions of some small proteins show clear evidence of a hierarchic process, whereas others, lacking detectable intermediates, do not. Evidence from folding intermediates and transition states suggests that folding begins locally, and that the formation of native secondary structure precedes the formation of tertiary interactions, not the reverse. Some notable examples in the literature have been interpreted to the contrary. For these examples, we have simulated the local structures that form when folding begins by using the LINUS program with nonlocal interactions turned off. Our results support a hierarchic model of protein folding.     

Oxidative folding of Amaranthus alpha-amylase inhibitor: disulfide bond formation and conformational folding

Oxidative folding is the fusion of native disulfide bond formation with conformational folding. This complex process is guided by two types of interactions: first, covalent interactions between cysteine residues, which transform into native disulfide bridges, and second, non-covalent interactions giving rise to secondary and tertiary protein structure. The aim of this work is to understand both types of interactions in the oxidative folding of Amaranthus alpha-amylase inhibitor (AAI) by providing information both at the level of individual disulfide species and at the level of amino acid residue conformation. The cystine-knot disulfides of AAI protein are stabilized in an interdependent manner, and the oxidative folding is characterized by a high heterogeneity of one-, two-, and three-disulfide intermediates. The formation of the most abundant species, the main folding intermediate, is favored over other species even in the absence of non-covalent sequential preferences. Time-resolved NMR and photochemically induced dynamic nuclear polarization spectroscopies were used to follow the oxidative folding at the level of amino acid residue conformation. Because this is the first time that a complete oxidative folding process has been monitored with these two techniques, their results were compared with those obtained at the level of an individual disulfide species. The techniques proved to be valuable for the study of conformational developments and aromatic accessibility changes along oxidative folding pathways. A detailed picture of the oxidative folding of AAI provides a model study that combines different biochemical and biophysical techniques for a fuller understanding of a complex process.     

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.     

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.     

Comparing the folding and misfolding energy landscapes of phosphoglycerate kinase

Partitioning of polypeptides between protein folding and amyloid formation is of outstanding pathophysiological importance. Using yeast phosphoglycerate kinase as model, here we identify the features of the energy landscape that decide the fate of the protein: folding or amyloidogenesis. Structure formation was initiated from the acid-unfolded state, and monitored by fluorescence from 10 ms to 20 days. Solvent conditions were gradually shifted between folding and amyloidogenesis, and the properties of the energy landscape governing structure formation were reconstructed. A gradual transition of the energy landscape between folding and amyloid formation was observed. In the early steps of both folding and misfolding, the protein searches through a hierarchically structured energy landscape to form a molten globule in a few seconds. Depending on the conditions, this intermediate either folds to the native state in a few minutes, or forms amyloid fibers in several days. As conditions are changed from folding to misfolding, the barrier separating the molten globule and native states increases, although the barrier to the amyloid does not change. In the meantime, the native state also becomes more unstable and the amyloid more stable. We conclude that the lower region of the energy landscape determines the final protein structure.     

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.     

Stepwise formation of alpha-helices during cytochrome c folding

Two models have been proposed to describe the folding pathways of proteins. The framework model assumes the initial formation of the secondary structures whereas the hydrophobic collapse model supposes their formation after the collapse of backbone structures. To differentiate between these models for real proteins, we have developed a novel CD spectrometer that enables us to observe the submillisecond time frame of protein folding and have characterized the timing of secondary structure formation in the folding process of cytochrome c (cyt c). We found that approximately 20% of the native helical content was organized in the first phase of folding, which is completed within milliseconds. Furthermore, we suggest the presence of a second intermediate, which has alpha-helical content resembling that of the molten globule state. Our results indicate that many of the alpha-helices are organized after collapse in the folding mechanism of cyt c.     

Protein folding includes oligomerization - examples from the endoplasmic reticulum and cytosol

A correct three-dimensional structure is a prerequisite for protein functionality, and therefore for life. Thus, it is not surprising that our cells are packed with proteins that assist protein folding, the process in which the native three-dimensional structure is formed. In general, plasma membrane and secreted proteins, as well as those residing in compartments along the endocytic and exocytic pathways, fold and oligomerize in the endoplasmic reticulum. The proteins residing in the endoplasmic reticulum are specialized in the folding of this subset of proteins, which renders this compartment a protein-folding factory. This review focuses on protein folding in the endoplasmic reticulum, and discusses the challenge of oligomer formation in the endoplasmic reticulum as well as the cytosol.     

A computer model to dynamically simulate protein folding: studies with crambin

The current work describes a simplified representation of protein structure with uses in the simulation of protein folding. The model assumes that a protein can be represented by a freely rotating rigid chain with a single atom approximating the effect of each side chain. Potentials describing the attraction or repulsion between different types of amino acids are determined directly from the distribution of amino acids in the database of known protein structures. The optimization technique of simulated annealing has been used to dynamically sample the conformations available to this simple model, allowing the protein to evolve from an extended, random coil into a compact globular structure. Many characteristics expected of true proteins, such as the sequence-dependent formation of secondary structure, the partitioning of hydrophobic residues, and specific disulfide pairing, are reproduced by the simulation, suggesting the model may accurately simulate the folding process.     

An analysis of protein folding pathways

We have developed a model of the protein folding process based on three primary assumptions: that burying of hydrophobic area is the dominant contribution to the relative free energy of a conformation, that a record of the folding process is largely preserved in the final structure, and that the denatured state is a random coil. Detailed folding pathways are identified for 19 protein structures. The picture of the folding process that emerges from this analysis is one of nucleation by regions of 8-16 residues. Nucleation sites then lead to larger structures by two mechanisms: propagation and diffusion/collision. A Monte Carlo simulation is used to follow the folding pathway when propagation is the dominant mechanism. Because detailed pathways are derived for each protein, the models are susceptible to experimental verification.     

Multi-State Proteins: Approach Allowing Experimental Determination of the Formation Order of Structure Elements in the Green Fluorescent Protein

The most complex problem in studying multi-state protein folding is the determination of the sequence of formation of protein intermediate states. A far more complex issue is to determine at what stages of protein folding its various parts (secondary structure elements) develop. The structure and properties of different intermediate states depend in particular on these parts. An experimental approach, named μ-analysis, which allows understanding the order of formation of structural elements upon folding of a multi-state protein was used in this study. In this approach the same elements of the protein secondary structure are “tested” by substitutions of single hydrophobic amino acids and by incorporation of cysteine bridges. Single substitutions of hydrophobic amino acids contribute to yielding information on the late stages of protein folding while incorporation of ss-bridges allows obtaining data on the initial stages of folding. As a result of such an μ-analysis, we have determined the order of formation of beta-hairpins upon folding of the green fluorescent protein.     

What determines protein folding type? An investigation of intrinsic structural properties and its implications for understanding folding mechanisms

Protein folding experiments demonstrate that the folding behaviors of many proteins can be roughly classified into two types: two-state kinetics and multi-state kinetics. Although the two types of protein folding kinetics have been observed for a long time, what determines the folding type of a protein is still largely unclear. The present work performed a comparative study based on a dataset of 43 two-state and 42 multi-state folders at different levels of proteins' intrinsic properties from the simplest sequence length to native structure topology. The results show that protein's amino acids composition and the long-range interaction-based topological complexity rather than secondary structure contents are the major determinants of protein folding type. Furthermore, a sequence-based folding type prediction achieved an accuracy of more than 80%. These findings implicate that there is no clear boundary between secondary and tertiary structure formation during the protein folding process and support the existence of a continuum of folding mechanism between the two ends of hierarchic and nucleation folding scenarios.