Role of arginine in the stabilization of proteins against aggregation

The amino acid arginine is frequently used as a solution additive to stabilize proteins against aggregation, especially in the process of protein refolding. Despite arginine's prevalence, the mechanism by which it stabilizes proteins is not presently understood. We propose that arginine deters aggregation by slowing protein-protein association reactions, with only a small concomitant effect on protein folding. The associated rate effect was observed experimentally in association of globular proteins (insulin and a monoclonal anti-insulin) and in refolding of carbonic anhydrase. We suggest that this effect arises because arginine is preferentially excluded from protein-protein encounter complexes but not from dissociated protein molecules. Such an effect is predicted by our gap effect theory [Baynes and Trout (2004) Biophys. J. 87, 1631] for "neutral crowder" additives such as arginine which are significantly larger than water but have only a small effect on the free energies of isolated protein molecules. The effect of arginine on refolding of carbonic anhydrase was also shown to be consistent with this hypothesis.     

Slow,Reversible, Coupled Folding and Binding of the Spectrin Tetramerization Domain

Many intrinsically disordered proteins (IDPs) are significantly unstructured under physiological conditions. A number of these IDPs have been shown to undergo coupled folding and binding reactions whereby they can gain structure upon association with an appropriate partner protein. In general, these systems display weaker binding affinities than do systems with association between completely structured domains, with micromolar K_d values appearing typical. One such system is the association between α- and β-spectrin, where two partially structured, incomplete domains associate to form a fully structured, three-helix bundle, the spectrin tetramerization domain. Here, we use this model system to demonstrate a method for fitting association and dissociation kinetic traces where, using typical biophysical concentrations, the association reactions are expected to be highly reversible. We elucidate the unusually slow, two-state kinetics of spectrin assembly in solution. The advantages of studying kinetics in this regime include the potential for gaining equilibrium constants as well as rate constants, and for performing experiments with low protein concentrations. We suggest that this approach would be particularly appropriate for high-throughput mutational analysis of two-state reversible binding processes.     

Does alpha-helix folding necessarily provide an energy source for the protein-lipid binding?

Lipid-induced alpha-helix folding, which occurs in many lipid surface-binding proteins and peptides such as apolipoproteins and synucleins, has been proposed to provide an energy source for protein-lipid interactions. We propose that in a system comprised of a phospholipid surface and a small polypeptide that is unfolded in solution and binds reversibly to lipid surface, helical folding involves expenditure of free energy as compared to a similar polypeptide that is alpha-helical in solution. This is a consequence of the entropic cost of helix folding that is illustrated in a simple thermodynamic model and exemplifies the general "key-into-lock" paradigm of protein-ligand binding. Even though this simple model does not explicitly address the protein-induced lipid re-arrangement and may not directly apply to large proteins that undergo significant tertiary structural changes upon lipid binding, it suggests that the notion of helix folding as an energy source for lipid binding should be treated with caution.     

Probing nature's knots: the folding pathway of a knotted homodimeric protein

The homodimeric protein YibK from Haemophilus influenzae belongs to a recently discovered superfamily of knotted proteins that has brought about a new protein-folding conundrum. Members of the alpha/beta-knot clan form deep trefoil knots in their native backbone structure, a topological feature that is currently unexplained in the protein-folding field. To help solve the puzzle of how a polypeptide chain can efficiently knot itself, the folding kinetics of YibK have been studied extensively and the results are reported here. Folding was monitored using probes for changes in both secondary and tertiary structure, and the monomer-dimer equilibrium was perturbed with a variety of solution conditions to allow characterisation of otherwise inaccessible states. Multiphasic kinetics were observed in the unfolding and refolding reactions of YibK, and under conditions where the dimer is favoured, dissociation and association were rate-limiting, respectively. A folding model consistent with all kinetic data is proposed: YibK appears to fold via two parallel pathways, partitioned by proline isomerisation events, to two distinct monomeric intermediates. These form a common third intermediate that is able to fold to native dimer. Kinetic simulations suggest that all intermediates are on-pathway. These results provide the valuable groundwork required to further understand how Nature codes for knot formation.     

Role of diffusion in the folding of the alpha subunit of tryptophan synthase from Escherichia coli

The rate-limiting step in the folding of the alpha subunit of tryptophan synthase has been proposed to be the association of two folding units. To probe the role of diffusion in this rate-limiting step, the urea-induced unfolding and refolding of the protein was examined in the presence of a number of viscosity-enhancing agents. The analysis was simplified by studying the effect of these agents on folding unit dissociation, the rate-limiting unfolding reaction, and the reverse of the rate-limiting step in refolding. In the presence of ethylene glycol, the relaxation times for unfolding to the same final conditions increased with increasing concentration of the cosolvent. When the effects of the cosolvent on protein stability were taken into account, the rates were found to show a unitary linear dependence on the viscosity of the solution. Similar results were obtained with glycerol and low concentrations of glucose, demonstrating that the effect is general and not specific to any viscogenic agent. These results clearly demonstrate that the rate-limiting folding unit association/dissociation reaction in the alpha subunit of tryptophan synthase involves a diffusional process.     

Effect of denaturant and protein concentrations upon protein refolding and aggregation: a simple lattice model.

We present a study of the competition between protein refolding and aggregation for simple lattice model proteins. The effect of solvent conditions (i.e., the denaturant concentration and the protein concentration) on the folding and aggregation behavior of a system of simple, two-dimensional lattice protein molecules has been investigated via (dynamic Monte Carlo simulations. The population profiles and aggregation propensities of the nine most populated intermediate configurations exhibit a complex dependence on the solution conditions that can be understood by considering the competition between intra- and interchain interactions. Some of these configurations are not even seen in isolated chain simulations; they are observed to be highly aggregation prone and are stabilized primarily by the aggregation reaction in multiple-chain systems. Aggregation arises from the association of partially folded intermediates rather than from the association of denatured random-coil states. The aggregation reaction dominates over the folding reaction at high protein concentration and low denaturant concentration, resulting in low refolding yields at those conditions. However, optimum folding conditions exist at which the refolding yield is a maximum, in agreement with some experimental observations.     

Protein self-association in the cell: a mechanism for fine tuning the level of macromolecular crowding?

A new role for protein self-association in the cell is discussed. An argument is advanced that when cellular protein is in its associated state the excluded volume component of the solution is minimized. Conversely, when cellular protein is in its dissociated state the excluded volume component of the solution is maximized. For proteins that make up a substantial fraction of the intracellular protein concentration, control of the self-association event thus presents itself as a means of regulating cellular processes that are influenced by different levels of volume exclusion. In this communication we examine how the control of protein association/dissociation might influence one such important process, namely the folding of a protein to a compact state.     

High-pressure NMR spectroscopy for characterizing folding intermediates and denatured states of proteins

Extensive structural studies using high-pressure NMR spectroscopy have recently been carried out on proteins, which potentially contribute to our understanding of the mechanisms of protein folding. Pressure shifts the conformational equilibrium from higher to lower volume conformers. If the pressure is varied, starting from the folded native structure, in many cases we observe intermediate conformers before the onset of total unfolding. This enables the investigation of details of the structure and thermodynamic characteristics of various intermediate conformers of proteins under equilibrium conditions. We can also examine pressure effects on the structure and stability of some typical denatured states such as helical denatured, molten globule, and unfolded states. The high-pressure NMR method can also be used to investigate association/dissociation equilibria of oligomeric or aggregated proteins. Beside direct observation of kinetic intermediates upon pressure jump, NMR structural investigations of equilibrium conformers under pressure provide information about the structures of kinetic intermediates during folding/unfolding reactions.     

The interaction of the GroEL chaperone with early kinetic intermediates of renaturing proteins inhibits the formation of their native structure

The main function of the chaperone GroEL is to prevent nonspecific association of nonnative protein chains and provide their correct folding. In the present work, the renaturation kinetics of three globular proteins (human alpha-lactalbumin, bovine carbonic anhydrase, and yeast phosphoglycerate kinase) in the presence of different molar excess of GroEL (up to 10-fold) was studied. It was shown that the formation of the native structure during the refolding of these proteins is retarded with an increase in GroEL molar excess due to the interaction of kinetic protein intermediates with the chaperone. Mg(2+)-ATP and Mg(2+)-ADP weaken this interaction and decrease the retarding effect of GroEL on the protein refolding kinetics. The theoretical modeling of protein folding in the presence of GroEL showed that the experimentally observed linear increase in the protein refolding half-time with increasing molar excess of GroEL must occur only when the protein adopts its native structure outside of GroEL (i.e. in the free state), while the refolding of the protein in the complex with GroEL is inhibited. The dissociation constants of GroEL complexed with the kinetic intermediates of the proteins studied were evaluated, and a simple mechanism of the functioning of GroEL as a molecular chaperone was proposed.     

The cooperativity of burst phase reactions explored.

The denaturant-dependence of the major, observable relaxation rates for folding (kobs) of ribonuclease HI from Escherichia coli (RNase H) and phage T4 lysozyme (T4L) reveal that, for both proteins, folding begins with the rapid and transient accumulation of intermediate species in a "burst phase" which precedes the rate-limiting formation of the native state; this is evidenced by a "rollover" in the folding limb of the rate profiles (kobs versus denaturant, or chevron plot). These rate profiles are most simply described by a three-state mechanism (unfolded-to-intermediate-to-native), which implies that the burst phase represents a transition between two distinct thermodynamic states. It is shown here that the equilibrium properties of these burst phase reactions can be equally well modeled by a mechanism involving a continuum of states where the free energy of each state is linearly related to its m-value (the parameter describing the linear relationship between free energy and denaturant). A numerical model is also developed to describe the time evolution of such a system, which exhibits nearly perfect exponential behavior. Both models emphasize how a continuum of states operating under a linear free energy relationship may behave like a two state system. Such a scheme finds experimental justification from an interpretation of recent native state hydrogen exchange data. The analytical model described for a continuum can account for the observed kinetic profiles of several other model proteins. The results, however, appear context specific, suggesting that burst phase reactions are not entirely random and non-specific. The results reported in this study have important implications for the concept of cooperativity in protein folding reactions.     

Monte Carlo study of substrate-induced folding and refolding of lattice proteins

Many proteins can switch from one conformation to another under the influence of an external driving force, such as the binding to a specific substrate. Using a simple lattice model we show that it is feasible to design protein-like lattice proteins that can have two different conformations, depending on whether or not they are bound to a substrate. We give three different examples of such substrate-induced refolding. In addition, we have explored substrate-induced folding of lattice proteins that do not fold when free in solution. We show that such proteins can bind with the same high specificity as prefolded protein, but have a considerably lower binding free energy. In this way proteins can bind to a substrate in a way that is highly specific, yet reversible.     

Designing a High Throughput Refolding Array Using a Combination of the GroEL Chaperonin and Osmolytes

Although GroE chaperonins and osmolytes had been used separately as protein folding aids, combining these two methods provides a considerable advantage for folding proteins that cannot fold with either osmolytes or chaperonins alone. This technique rapidly identifies superior folding solution conditions for a broad array of proteins that are difficult or impossible to fold by other methods. While testing the broad applicability of this technique, we have discovered that osmolytes greatly simplify the chaperonin reaction by eliminating the requirement for the co-chaperonin GroES which is normally involved in encapsulating folding proteins within the GroEL–GroES cavity. Therefore, combinations of soluble or immobilized GroEL, osmolytes and ATP or even ADP are sufficient to refold the test proteins. The first step in the chaperonin/osmolyte process is to form a stable long-lived chaperonin–substrate protein complex in the absence of nucleotide. In the second step, different osmolyte solutions are added along with nucleotides, thus forming a ‘folding array’ to identify superior folding conditions. The stable chaperonin–substrate protein complex can be concentrated or immobilized prior to osmolyte addition. This procedure prevents-off pathway aggregation during folding/refolding reactions and more importantly allows one to refold proteins at concentrations (~mg/ml) that are substantially higher than the critical aggregation concentration for given protein. This technique can be used for successful refolding of proteins from purified inclusion bodies. Recently, other investigators have used our chaperonin/osmolyte method to demonstrate that a mutant protein that misfolds in human disease can be rescued by GroEL/osmolyte system. Soluble or immobilized GroEL can be easily removed from the released folded protein using simple separation techniques. The method allows for isolation of folded monomeric or oligomeric proteins in quantities sufficient for X-ray crystallography or NMR structural determinations.     

Protein folding and protein refolding.

The functional three-dimensional structure of proteins is determined solely by their amino acid sequences. Protein folding occurs spontaneously beginning with the formation of local secondary structure concomitant with a compaction of the molecule. Secondary structure elements subsequently interact to form subdomains and domains stabilized by tertiary interactions. Disulfide bond formation, and cis-trans isomerization of X-Pro peptide bonds, as the rate-limiting folding reactions, are enzymatically catalyzed during protein folding in the cell. Although folding of domains is fast enough to occur cotranslationally in vivo, such vectorial folding on the ribosome is not essential for attainment of the native structure of a protein. Slow steps on the pathway to the functional protein structure are docking reactions of domains, association of subunits, or reshuffling reactions at the oligomer level. Aggregation as a competing side reaction is prevented, and the kinetic partition between competing polypeptide folding and translocation reactions is regulated by chaperone proteins binding to incompletely folded polypeptides.     

Folding of a model three-helix bundle protein: a thermodynamic and kinetic analysis.

The kinetics and thermodynamics of an off-lattice model for a three-helix bundle protein are investigated as a function of a bias gap parameter that determines the energy difference between native and non-native contacts. A simple dihedral potential is used to introduce the tendency to form right-handed helices. For each value of the bias parameter, 100 trajectories of up to one microsecond are performed. Such statistically valid sampling of the kinetics is made possible by the use of the discrete molecular dynamics method with square-well interactions. This permits much faster simulations for off-lattice models than do continuous potentials. It is found that major folding pathways can be defined, although ensembles with considerable structural variation are involved. The large gap models generally fold faster than those with a smaller gap. For the large gap models, the kinetic intermediates are non-obligatory, while both obligatory and non-obligatory intermediates are present for small gap models. Certain large gap intermediates have a two-helix microdomain with one helix extended outward (as in domain-swapped dimers); the small gap intermediates have more diverse structures. The importance of studying the kinetic, as well as the thermodynamics, of folding for an understanding of the mechanism is discussed and the relation between kinetic and equilibrium intermediates is examined. It is found that the behavior of this model system has aspects that encompass both the "new" view and the "old" view of protein folding.     

Comparison between long-range interactions and contact order in determining the folding rate of two-state proteins: application of long-range order to folding rate prediction.

The contact order is believed to be an important factor for understanding protein folding mechanisms. In our earlier work, we have shown that the long-range interactions play a vital role in protein folding. In this work, we analyzed the contribution of long-range contacts to determine the folding rate of two-state proteins. We found that the residues that are close in space and are separated by at least ten to 15 residues in sequence are important determinants of folding rates, suggesting the presence of a folding nucleus at an interval of approximately 25 residues. A novel parameter "long-range order" has been proposed to predict protein folding rates. This parameter shows as good a relationship with the folding rate of two-state proteins as contact order. Further, we examined the minimum limit of residue separation to determine the long-range contacts for different structural classes. We observed an excellent correlation between long-range order and folding rate for all classes of globular proteins. We suggest that in mixed-class proteins, a larger number of residues can serve as folding nuclei compared to all-alpha and all-beta proteins. A simple statistical method has been developed to predict the folding rates of two-state proteins using the long-range order that produces an agreement with experimental results that is better or comparable to other methods in the literature.     

An automated in vitro protein folding screen applied to a human dynactin subunit

The preparation of proteins for structural and functional analysis using the Escherichia coli expression system is often hampered by the formation of insoluble intracellular protein aggregates (inclusion bodies). Transferring those proteins into their native states by in vitro protein folding requires screening for the best buffer conditions and suitable additives. However, it is difficult to assess the success of such a screen if no biological assay is available. We established a fully automated folding screen and a system to detect folded protein that is based on analytical hydrophobic interaction chromatography and tryptophan fluorescence spectroscopy. The system was evaluated with two model enzymes (carbonic anhydrase II and malate dehydrogenase), and was successfully applied to the folding of the p22 subunit of human dynactin, which is expressed in inclusion bodies in E. coli. The described screen allows for high-throughput folding analysis of inclusion body proteins for structural and functional analyses.     

Thermodynamic and kinetic characterization of the association of triosephosphate isomerase: the role of diffusion

It is known that diffusion plays a central role in the folding of small monomeric proteins and in the rigid-body association of proteins, however, the role of diffusion in the association of the folding intermediates of oligomeric proteins has been scarcely explored. In this work, catalytic activity and fluorescence measurements were used to study the effect of viscosity in the unfolding and refolding of the homodimeric enzyme triosephosphate isomerase from Saccharamyces cerevisiae. Two transitions were found by equilibrium and kinetic experiments, suggesting a three-state model with a monomeric intermediate. Glycerol barely affects DeltaG(0)(fold) whereas DeltaG(0)(assoc) becomes more favourable in the presence of the cosolvent. From 0 to 60% (v/v) glycerol, the association rate constant showed a near unitary dependence on solvent viscosity. However, at higher glycerol concentrations deviations from Kramers theory were observed. The dissociation rate constant showed a viscosity effect much higher than one. This may be related to secondary effects such as short-range glycerol-induced repulsion between monomers. Nevertheless, after comparison under isostability conditions, a slope near one was also observed for the dissociation rate. These results strongly suggest that the bimolecular association producing the native dimer is limited by diffusional events of the polypeptide chains through the solvent.     

Convenient and efficient in vitro folding of disulfide-containing globular protein from crude bacterial inclusion bodies

We investigated how the folding yield of disulfide-containing globular proteins having positive net charges from crude bacterial inclusion bodies was affected by additives in the folding buffer. In screening folding conditions for human ribonucleases and its derivative, we found that addition of salt (about 0.4 M) to a folding buffer increased the folding yield. This suggested that electrostatic interaction between polyanionic impurities such as nucleic acids and cationic unfolded protein led to the formation of aggregates under the low-salt conditions. Since inclusion bodies were found to contain nucleic acids regardless of the electrostatic nature of the expressed protein, the electrostatic interaction between phosphate moieties of nucleic acids and basic amino acid residues of a denatured protein may be large enough to cause aggregation, and therefore the addition of salt in a folding buffer may generally be useful for promotion of protein folding from crude inclusion bodies. We further systematically investigated additives such as glycerol, guanidium chloride, and urea that are known to act as chemical chaperons, and found that these additives, together with salt, synergistically improved folding yield. This study, suggesting that the addition of salt into the folding buffer is one of the crucial points to be considered, may pave the way for a systematic investigation of the folding conditions of disulfide-containing foreign proteins from crude bacterial inclusion bodies.     

An Overlapping Region between the Two Terminal Folding Units of the Outer Surface Protein A (OspA) Controls Its Folding Behavior

Although many naturally occurring proteins consist of multiple domains, most studies on protein folding to date deal with single-domain proteins or isolated domains of multi-domain proteins. Studies of multi-domain protein folding are required for further advancing our understanding of protein folding mechanisms. Borrelia outer surface protein A (OspA) is a β-rich two-domain protein, in which two globular domains are connected by a rigid and stable single-layer β-sheet. Thus, OspA is particularly suited as a model system for studying the interplays of domains in protein folding. Here, we studied the equilibria and kinetics of the urea-induced folding–unfolding reactions of OspA probed with tryptophan fluorescence and ultraviolet circular dichroism. Global analysis of the experimental data revealed compelling lines of evidence for accumulation of an on-pathway intermediate during kinetic refolding and for the identity between the kinetic intermediate and a previously described equilibrium unfolding intermediate. The results suggest that the intermediate has the fully native structure in the N-terminal domain and the single layer β-sheet, with the C-terminal domain still unfolded. The observation of the productive on-pathway folding intermediate clearly indicates substantial interactions between the two domains mediated by the single-layer β-sheet. We propose that a rigid and stable intervening region between two domains creates an overlap between two folding units and can energetically couple their folding reactions.     

Cysteine-cysteine contact preference leads to target-focusing in protein folding

We perform a statistical analysis of amino-acid contacts to investigate possible preferences of amino-acid interactions. We include in the analysis only tertiary contacts, because they are less constrained--compared to secondary contacts--by proteins' backbone rigidity. Using proteins from the protein data bank, our analysis reveals an unusually high frequency of cysteine pairings relative to that expected from random. To elucidate the possible effects of cysteine interactions in folding, we perform molecular simulations on three cysteine-rich proteins. In particular, we investigate the difference in folding dynamics between a Gō-like model (where attraction only occurs between amino acids forming a native contact) and a variant model (where attraction between any two cysteines is introduced to mimic the formation/dissociation of native/nonnative disulfide bonds). We find that when attraction among cysteines is nonspecific and comparable to a solvent-averaged interaction, they produce a target-focusing effect that expedites folding of cysteine-rich proteins as a result of a reduction of conformational search space. In addition, the target-focusing effect also helps reduce glassiness by lowering activation energy barriers and kinetic frustration in the system. The concept of target-focusing also provides a qualitative understanding of a correlation between the rates of protein folding and parameters such as contact order and total contact distance.