Fast and slow tracks in lysozyme folding: insight into the role of domains in the folding process

Molecular insights into the phase transition of lysozyme into amyloid nanostructures: Implications of therapeutic strategies in diverse pathological conditions

Lysozyme, a well-known bacteriolytic enzyme, exhibits a fascinating yet complex behavior when it comes to protein aggregation. Under certain conditions, this enzyme undergoes flexible transformation, transitioning from partially unfolded intermediate units of native conformers into complex cross-β-rich nano fibrillar amyloid architectures. Formation of such lysozyme amyloids has been implicated in a multitude of pathological and medical severities, like hepatic dysfunction, hepatomegaly, splenic rupture as well as spleen dysfunction, nephropathy, sicca syndrome, renal dysfunction, renal amyloidosis, and systemic amyloidosis. In this comprehensive review, we have attempted to provide in-depth insights into the aggregating behavior of lysozyme across a spectrum of variables, including concentrations, temperatures, pH levels, and mutations. Our objective is to elucidate the underlying mechanisms that govern lysozyme's aggregation process and to unravel the complex interplay between its structural attributes. Moreover, this work has critically examined the latest advancements in the field, focusing specifically on novel strategies and systems, that have been implemented to delay or inhibit the lysozyme amyloidogenesis. Apart from this, we have tried to explore and advance our fundamental understanding of the complex processes involved in lysozyme aggregation. This will help the research community to lay a robust foundation for screening, designing, and formulating targeted anti-amyloid therapeutics offering improved treatment modalities and interventions not only for lysozyme-linked amyloidopathy but for a wide range of amyloid-related disorders.

Heterotrimeric collagen helix with high specificity of assembly results in a rapid rate of folding

The most abundant natural collagens form heterotrimeric triple helices. Synthetic mimics of collagen heterotrimers have been found to fold slowly, even compared to the already slow rates of homotrimeric helices. These prolonged folding rates are not understood. Here we compare the stabilities, specificities and folding rates of three heterotrimeric collagen mimics designed through a computationally assisted approach. The crystal structure of one ABC-type heterotrimer verified a well-controlled composition and register and elucidated the geometry of pairwise cation-π and axial and lateral salt bridges in the assembly. This collagen heterotrimer folds much faster (hours versus days) than comparable, well-designed systems. Circular dichroism and NMR data suggest the folding is frustrated by unproductive, competing heterotrimer species and these species must unwind before refolding into the thermodynamically favoured assembly. The heterotrimeric collagen folding rate is inhibited by the introduction of preformed competing triple-helical assemblies, which suggests that slow heterotrimer folding kinetics are dominated by the frustration of the energy landscape caused by competing triple helices.

Structural Lessons From the Mutant Proinsulin Syndrome

Insight into folding mechanisms of proinsulin has been provided by analysis of dominant diabetes-associated mutations in the human insulin gene (INS). Such mutations cause pancreatic β-cell dysfunction due to toxic misfolding of a mutant proinsulin and impairment in trans of wild-type insulin secretion. Anticipated by the "Akita" mouse (a classical model of monogenic diabetes mellitus; DM), this syndrome illustrates the paradigm endoreticulum (ER) stress leading to intracellular proteotoxicity. Diverse clinical mutations directly or indirectly perturb native disulfide pairing leading to protein misfolding and aberrant aggregation. Although most introduce or remove a cysteine (Cys; leading in either case to an unpaired thiol group), non-Cys-related mutations identify key determinants of folding efficiency. Studies of such mutations suggest that the hormone's evolution has been constrained not only by structure-function relationships, but also by the susceptibility of its single-chain precursor to impaired foldability. An intriguing hypothesis posits that INS overexpression in response to peripheral insulin resistance likewise leads to chronic ER stress and β-cell dysfunction in the natural history of non-syndromic Type 2 DM. Cryptic contributions of conserved residues to folding efficiency, as uncovered by rare genetic variants, define molecular links between biophysical principles and the emerging paradigm of Darwinian medicine: Biosynthesis of proinsulin at the edge of non-foldability provides a key determinant of "diabesity" as a pandemic disease of civilization.

Diabetes mellitus due to toxic misfolding of proinsulin variants

BACKGROUND

Dominant mutations in the human insulin gene (INS) lead to pancreatic β-cell dysfunction and diabetes mellitus (DM) due to toxic misfolding of a mutant proinsulin. Analogous to a classical mouse model of monogenic DM ("Akita"), this syndrome highlights the susceptibility of β-cells to endoreticulum (ER) stress due to protein misfolding and aberrant aggregation.

SCOPE OF REVIEW

Diverse clinical mutations directly or indirectly perturb native disulfide pairing. Whereas most introduce or remove a cysteine (Cys; leading in either case to an unpaired thiol group), non-Cys-related mutations identify key determinants of folding efficiency. Studies of such mutations suggest that the hormone's evolution has been constrained not only by structure-function relationships but also by the susceptibility of its single-chain precursor to impaired foldability. An intriguing hypothesis posits that INS overexpression in response to peripheral insulin resistance likewise leads to chronic ER stress and β-cell dysfunction in the natural history of nonsyndromic Type 2 DM.

MAJOR CONCLUSIONS

Cryptic contributions of conserved residues to folding efficiency, as uncovered by rare genetic variants, define molecular links between biophysical principles and the emerging paradigm of Darwinian medicine: Biosynthesis of proinsulin at the edge of nonfoldability provides a key determinant of "diabesity" as a pandemic disease of civilization.

Effects of Protein Unfolding on Aggregation and Gelation in Lysozyme Solutions

In this work, we investigate the role of folding/unfolding equilibrium in protein aggregation and formation of a gel network. Near the neutral pH and at a low buffer ionic strength, the formation of the gel network around unfolding conditions prevents investigations of protein aggregation. In this study, by deploying the fact that in lysozyme solutions the time of folding/unfolding is much shorter than the characteristic time of gelation, we have prevented gelation by rapidly heating the solution up to the unfolding temperature (~80 °C) for a short time (~30 min.) followed by fast cooling to the room temperature. Dynamic light scattering measurements show that if the gelation is prevented, nanosized irreversible aggregates (about 10–15 nm radius) form over a time scale of 10 days. These small aggregates persist and aggregate further into larger aggregates over several weeks. If gelation is not prevented, the nanosized aggregates become the building blocks for the gel network and define its mesh length scale. These results support our previously published conclusion on the nature of mesoscopic aggregates commonly observed in solutions of lysozyme, namely that aggregates do not form from lysozyme monomers in their native folded state. Only with the emergence of a small fraction of unfolded proteins molecules will the aggregates start to appear and grow.

Role of Disulfide Bonds and Topological Frustration in the Kinetic Partitioning of Lysozyme Folding Pathways

Disulfide bonds in proteins can strongly influence the folding pathways by constraining the conformational space. Lysozyme has four disulfide bonds and is widely studied for its antibacterial properties. Experiments on lysozyme infer that the protein folds through a fast and a slow pathway. However, the reasons for the kinetic partitioning in the folding pathways are not completely clear. Using a coarse-grained protein model and simulations, we show that two out of the four disulfide bonds, which are present in the α-domain of lysozyme, are responsible for the slow folding pathway. In this pathway, a kinetically trapped intermediate state, which is close to the native state, is populated. In this state, the orientations of α-helices present in the α-domain are misaligned relative to each other. The protein in this state has to partially unfold by breaking down the interhelical contacts between the misaligned helices to fold to the native state. However, the topological constraints due to the two disulfide bonds present in the α-domain make the protein less flexible, and it is trapped in this conformation for hundreds of milliseconds. On disabling these disulfide bonds, we find that the kinetically trapped intermediate state and the slow folding pathway disappear. Simulations mimicking the folding of protein without disulfide bonds under oxidative conditions show that the native disulfide bonds are formed as the protein folds, indicating that folding guides the formation of disulfide bonds. The sequence of formation of the disulfide bonds is Cys64-Cys80 → Cys76-Cys94 → Cys30-Cys115 → Cys6-Cys127. Any disulfide bond that forms before its precursor in the sequence has to break and follow the sequence for the protein to fold. These results show that lysozyme also serves as a very good model system to probe the role of disulfide bonds and topological frustration in protein folding. The predictions from the simulations can be verified by single-molecule fluorescence resonance energy transfer or single-molecule pulling experiments, which can probe heterogeneity in the folding pathways.

Direct observation of parallel folding pathways revealed using a symmetric repeat protein system

Although progress has been made to determine the native fold of a polypeptide from its primary structure, the diversity of pathways that connect the unfolded and folded states has not been adequately explored. Theoretical and computational studies predict that proteins fold through parallel pathways on funneled energy landscapes, although experimental detection of pathway diversity has been challenging. Here, we exploit the high translational symmetry and the direct length variation afforded by linear repeat proteins to directly detect folding through parallel pathways. By comparing folding rates of consensus ankyrin repeat proteins (CARPs), we find a clear increase in folding rates with increasing size and repeat number, although the size of the transition states (estimated from denaturant sensitivity) remains unchanged. The increase in folding rate with chain length, as opposed to a decrease expected from typical models for globular proteins, is a clear demonstration of parallel pathways. This conclusion is not dependent on extensive curve-fitting or structural perturbation of protein structure. By globally fitting a simple parallel-Ising pathway model, we have directly measured nucleation and propagation rates in protein folding, and have quantified the fluxes along each path, providing a detailed energy landscape for folding. This finding of parallel pathways differs from results from kinetic studies of repeat-proteins composed of sequence-variable repeats, where modest repeat-to-repeat energy variation coalesces folding into a single, dominant channel. Thus, for globular proteins, which have much higher variation in local structure and topology, parallel pathways are expected to be the exception rather than the rule.

Diabetes mellitus due to the toxic misfolding of proinsulin variants

Dominant mutations in the human insulin gene can lead to pancreatic β-cell dysfunction and diabetes mellitus due to toxic folding of a mutant proinsulin. Analogous to a classical mouse model (the Akita mouse), this monogenic syndrome highlights the susceptibility of human β-cells to endoreticular stress due to protein misfolding and aberrant aggregation. The clinical mutations directly or indirectly perturb native disulfide pairing. Whereas the majority of mutations introduce or remove a cysteine (leading in either case to an unpaired residue), non-cysteine-related mutations identify key determinants of folding efficiency. Studies of such mutations suggest that the evolution of insulin has been constrained not only by its structure and function, but also by the susceptibility of its single-chain precursor to impaired foldability.

Constraint Network Analysis (CNA): a Python software package for efficiently linking biomacromolecular structure, flexibility, (thermo-)stability, and function

For deriving maximal advantage from information on biomacromolecular flexibility and rigidity, results from rigidity analyses must be linked to biologically relevant characteristics of a structure. Here, we describe the Python-based software package Constraint Network Analysis (CNA) developed for this task. CNA functions as a front- and backend to the graph-based rigidity analysis software FIRST. CNA goes beyond the mere identification of flexible and rigid regions in a biomacromolecule in that it (I) provides a refined modeling of thermal unfolding simulations that also considers the temperature-dependence of hydrophobic tethers, (II) allows performing rigidity analyses on ensembles of network topologies, either generated from structural ensembles or by using the concept of fuzzy noncovalent constraints, and (III) computes a set of global and local indices for quantifying biomacromolecular stability. This leads to more robust results from rigidity analyses and extends the application domain of rigidity analyses in that phase transition points ("melting points") and unfolding nuclei ("structural weak spots") are determined automatically. Furthermore, CNA robustly handles small-molecule ligands in general. Such advancements are important for applying rigidity analysis to data-driven protein engineering and for estimating the influence of ligand molecules on biomacromolecular stability. CNA maintains the efficiency of FIRST such that the analysis of a single protein structure takes a few seconds for systems of several hundred residues on a single core. These features make CNA an interesting tool for linking biomacromolecular structure, flexibility, (thermo-)stability, and function. CNA is available from http://cpclab.uni-duesseldorf.de/software for nonprofit organizations.

Optimal use of data in parallel tempering simulations for the construction of discrete-state Markov models of biomolecular dynamics

Parallel tempering (PT) molecular dynamics simulations have been extensively investigated as a means of efficient sampling of the configurations of biomolecular systems. Recent work has demonstrated how the short physical trajectories generated in PT simulations of biomolecules can be used to construct the Markov models describing biomolecular dynamics at each simulated temperature. While this approach describes the temperature-dependent kinetics, it does not make optimal use of all available PT data, instead estimating the rates at a given temperature using only data from that temperature. This can be problematic, as some relevant transitions or states may not be sufficiently sampled at the temperature of interest, but might be readily sampled at nearby temperatures. Further, the comparison of temperature-dependent properties can suffer from the false assumption that data collected from different temperatures are uncorrelated. We propose here a strategy in which, by a simple modification of the PT protocol, the harvested trajectories can be reweighted, permitting data from all temperatures to contribute to the estimated kinetic model. The method reduces the statistical uncertainty in the kinetic model relative to the single temperature approach and provides estimates of transition probabilities even for transitions not observed at the temperature of interest. Further, the method allows the kinetics to be estimated at temperatures other than those at which simulations were run. We illustrate this method by applying it to the generation of a Markov model of the conformational dynamics of the solvated terminally blocked alanine peptide.

Lysozyme adsorption on polyethylene surfaces: why are long simulations needed?

The adsorption of lysozyme onto a polyethylene (PE) surface in an aqueous environment was investigated via molecular dynamics (MD) simulation. The adsorption can be divided into three processes: diffusion to the surface, dehydration induced by hydrophobic surface-protein interactions, and denaturation. The dehydration process is very long, around 70 ns. Structural deformations start soon after the protein reaches the surface and continue during the whole trajectory. The hydrophobic residues are slowly driven toward the surface, inducing changes in the protein's secondary structure. The protein's secondary structural components near the surface are more disturbed than those farther away from the surface. The lysozyme is adsorbed with its long axis parallel to the surface and displays an anisotropic mobility on the surface that is probably due to the intrinsic structure of the PE surface. Our study demonstrates the need for long-time atomistic simulation in order to gain a complete understanding of the adsorption process.

Probing molecular kinetics with Markov models: metastable states, transition pathways and spectroscopic observables

Markov (state) models (MSMs) have attracted a lot of interest recently as they (1) can probe long-term molecular kinetics based on short-time simulations, (2) offer a way to analyze great amounts of simulation data with relatively little subjectivity of the analyst, (3) provide insight into microscopic quantities such as the ensemble of transition pathways, and (4) allow simulation data to be reconciled with measurement data in a rigorous and explicit way. Here we sketch our current perspective of Markov models and explain in short their theoretical basis and assumptions. We describe transition path theory which allows the entire ensemble of protein folding pathways to be investigated and that combines naturally with Markov models. Experimental observations can be naturally linked to Markov models with the dynamical fingerprint theory, by which experimentally observable timescales can be equipped with an understanding of the structural rearrangement processes that take place at these timescales. The concepts of this paper are illustrated by a simple kinetic model of protein folding.

Rapid collapse into a molten globule is followed by simple two-state kinetics in the folding of lysozyme from bacteriophage λ

Stopped-flow fluorescence and circular dichroism spectroscopy have been used in combination with quenched-flow hydrogen exchange labeling, monitored by two-dimensional NMR and electrospray ionization mass spectrometry, to investigate the folding kinetics of lysozyme from bacteriophage λ (λ lysozyme) at pH 5.6, 20 °C. The first step in the folding of λ lysozyme occurs very rapidly (τ < 1 ms) after refolding is initiated and involves both hydrophobic collapse and formation of a high content of secondary structure but only weak protection from (1)H/(2)H exchange and no fixed tertiary structure organization. This early folding step is reflected in the dead-time events observed in the far-UV CD and ANS fluorescence experiments. Following accumulation of this kinetic molten globule species, the secondary structural elements are stabilized and the majority (ca. 88%) of refolding molecules acquire native-like properties in a highly cooperative two-state process, with τ = 0.15 ± 0.03 s. This is accompanied by the acquisition of substantial native-like protection from hydrogen exchange. A double-mixing experiment and the absence of a denaturant effect reveal that slow (τ = 5 ± 1 s) folding of the remaining (ca. 12%) molecules is rate limited by the cis/trans isomerization of prolines that are trans in the folded enzyme. In addition, native state hydrogen exchange and classical denaturant unfolding experiments have been used to characterize the thermodynamic properties of the enzyme. In good agreement with previous crystallographic evidence, our results show that λ lysozyme is a highly dynamic protein, with relatively low conformational stability (ΔG°(N-U) = 25 ± 2 kJ·mol(-1)).

Structural heterogeneity of 6 M GdmCl-denatured proteins: implications for the mechanism of protein folding

An in vitro experiment with protein folding is typically initiated with 6 M GdmCl-denatured proteins, which are generally considered fully unfolded. However, studies conducted by various laboratories have shown that many 6 M GdmCl-denatured proteins are structurally heterogeneous and still retain nativelike residual structures. The extent of conformational heterogeneity of the 6 M GdmCl-denatured protein has significant implications for the folding landscape as well as the interpretation of the observed early stage folding mechanism. Using the method of disulfide scrambling, we are able to gain rough insight into the diverse structural properties of 6 M GdmCl-denatured proteins. It demonstrates that most 6 M GdmCl-denatured proteins are approximately fully denatured, but partially unfolded. Most of them comprise diverse conformational isomers. We review here the cumulative evidence obtained from various laboratories and also provide experimental data obtained in our laboratory.

Constructing the equilibrium ensemble of folding pathways from short off-equilibrium simulations

Characterizing the equilibrium ensemble of folding pathways, including their relative probability, is one of the major challenges in protein folding theory today. Although this information is in principle accessible via all-atom molecular dynamics simulations, it is difficult to compute in practice because protein folding is a rare event and the affordable simulation length is typically not sufficient to observe an appreciable number of folding events, unless very simplified protein models are used. Here we present an approach that allows for the reconstruction of the full ensemble of folding pathways from simulations that are much shorter than the folding time. This approach can be applied to all-atom protein simulations in explicit solvent. It does not use a predefined reaction coordinate but is based on partitioning the state space into small conformational states and constructing a Markov model between them. A theory is presented that allows for the extraction of the full ensemble of transition pathways from the unfolded to the folded configurations. The approach is applied to the folding of a PinWW domain in explicit solvent where the folding time is two orders of magnitude larger than the length of individual simulations. The results are in good agreement with kinetic experimental data and give detailed insights about the nature of the folding process which is shown to be surprisingly complex and parallel. The analysis reveals the existence of misfolded trap states outside the network of efficient folding intermediates that significantly reduce the folding speed.

What have we learned from the studies of two-state folders, and what are the unanswered questions about two-state protein folding?

Small proteins with globular structures often fold by simple all-or-none mechanisms, both in an equilibrium and a kinetic sense, despite the very large number of partly folded conformations available. This type of 'two-state' folding will be discussed in terms of experimental tests, underlying molecular mechanisms, and limits to two-state behavior. Factors that appear to be important for two-state folding include topology (sequence distance of contacts in the native structure), molecular cooperativity and local energy distribution. Because their local stability distributions and cooperativities can be dissected and analyzed separately from topological features, recent studies of the folding of symmetric proteins will be discussed as a means to better understand the origins of two-state folding.

Kinetic model for the coupling between allosteric transitions in GroEL and substrate protein folding and aggregation

The bacterial chaperonin GroEL and the co-chaperonin GroES assist in the folding of a number of structurally unrelated substrate proteins (SPs). In the absence of chaperonins, SP folds by the kinetic partitioning mechanism (KPM), according to which a fraction of unfolded molecules reaches the native state directly, while the remaining fraction gets trapped in a potentially aggregation-prone misfolded state. During the catalytic reaction cycle, GroEL undergoes a series of allosteric transitions (T<-->R-->R"-->T) triggered by SP capture, ATP binding and hydrolysis, and GroES binding. We developed a general kinetic model that takes into account the coupling between the rates of the allosteric transitions and the folding and aggregation of the SP. Our model, in which the GroEL allosteric rates and SP-dependent folding and aggregation rates are independently varied without prior assumption, quantitatively fits the GroEL concentration-dependent data on the yield of native ribulose bisphosphate carboxylase/oxygenase (Rubisco) as a function of time. The extracted kinetic parameters for the GroEL reaction cycle are consistent with the available values from independent experiments. In addition, we also obtained physically reasonable parameters for the kinetic steps in the reaction cycle that are difficult to measure. If experimental values for GroEL allosteric rates are used, the time-dependent changes in native-state yield at eight GroEL concentrations can be quantitatively fit using only three SP-dependent parameters. The model predicts that the differences in the efficiencies (as measured by yields of the native state) of GroEL, single-ring mutant (SR1), and variants of SR1, in the rescue of mitochondrial malate dehydrogenase, citrate synthase, and Rubisco, are related to the large variations in the allosteric transition rates. We also show that GroEL/S mutants that efficiently fold one SP at the expense of all others are due to a decrease in the rate of a key step in the reaction cycle, which implies that wild-type GroEL has evolved as a compromise between generality and specificity. We predict that, under maximum loading conditions and saturating ATP concentration, the efficiency of GroEL (using parameters for Rubisco) depends predominantly on the rate of R-->R" transition, while the equilibrium constant of the T<-->R has a small effect only. Both under sub- and superstoichiometric GroEL concentrations, enhanced efficiency is achieved by rapid turnover of the reaction cycle, which is in accord with the predictions of the iterative annealing mechanism. The effects are most dramatic at substoichiometric conditions (most relevant for in vivo situations) when SP aggregation can outcompete capture of SP by chaperonins.

The protein folding problem

The "protein folding problem" consists of three closely related puzzles: (a) What is the folding code? (b) What is the folding mechanism? (c) Can we predict the native structure of a protein from its amino acid sequence? Once regarded as a grand challenge, protein folding has seen great progress in recent years. Now, foldable proteins and nonbiological polymers are being designed routinely and moving toward successful applications. The structures of small proteins are now often well predicted by computer methods. And, there is now a testable explanation for how a protein can fold so quickly: A protein solves its large global optimization problem as a series of smaller local optimization problems, growing and assembling the native structure from peptide fragments, local structures first.

Fibril formation of lysozyme upon interaction with sodium dodecyl sulfate at pH 9.2

Fibril formation seems to be a general property of all proteins. Its occurrence in hen or human lysozyme depends on certain conditions, namely acidic pHs or the presence of some additives. This paper studies the interaction of lysozyme with sodium dodecyl sulfate (SDS) at pH 9.2, using UV-visible spectrophotometry, circular dichroism (CD) spectropolarimetry, electron microscopy (EM) and chemometry. Based on observations such as the strange increase in absorbance at 650nm (pH 9.2) and the presence of intermediates, it is assumed that lysozyme fibrils have been formed at pH 9.2 in the presence of SDS as an anionic surfactant. Thioflavin T emission fluorescence and an EM image confirmed this assumption. beta-cyclodextrin was then used as a turbidity inhibitor to establish its effect on the distribution of intermediates that participate in fibril formation.

Folding of tandem-linked domains

One of the factors, which influences protein folding in vivo, is a linkage of protein domains into multidomain tandems. However, relatively little is known about the impact of domain connectivity on protein folding mechanisms. In this article, we use coarse grained models of proteins to explore folding of tandem-linked domains (TLD). We found TLD folding to follow two scenarios. In the first, the tandem connectivity produces relatively minor impact on folding and the mechanisms of folding of tandem-linked and single domains remain similar. The second scenario involves qualitative changes in folding mechanism because of tandem linkage. As a result, protein domains, which fold via two-state mechanism as single isolated domains, may form new stable intermediates when inserted into tandems. The new intermediates are created by topological constraints imposed by the linkers between domains. In both cases tandem linkage slows down folding. We propose that the impact of tandem connectivity can be minimized, if the terminal secondary structure elements (SSEs) are flexible. In particular, two factors appear to facilitate TLD folding: (1) the interactions between terminal SSE are poorly ordered in the folding transition state, whereas nonterminal SSE are better structured, (2) the interactions between terminal SSE are weak in the native state. We apply these findings to wild-type proteins by examining experimental phi-value data and by performing all-atom molecular dynamics simulations. We show that immunoglobulin-like domains appear to utilize the factors, which minimize the impact of tandem connectivity on their folding. Several single domain proteins, which are likely to misfold in tandems, are also identified.

A unified mechanism for protein folding: predetermined pathways with optional errors

There is a fundamental conflict between two different views of how proteins fold. Kinetic experiments and theoretical calculations are often interpreted in terms of different population fractions folding through different intermediates in independent unrelated pathways (IUP model). However, detailed structural information indicates that all of the protein population folds through a sequence of intermediates predetermined by the foldon substructure of the target protein and a sequential stabilization principle. These contrary views can be resolved by a predetermined pathway--optional error (PPOE) hypothesis. The hypothesis is that any pathway intermediate can incorporate a chance misfolding error that blocks folding and must be reversed for productive folding to continue. Different fractions of the protein population will then block at different steps, populate different intermediates, and fold at different rates, giving the appearance of multiple unrelated pathways. A test of the hypothesis matches the two models against extensive kinetic folding results for hen lysozyme which have been widely cited in support of independent parallel pathways. The PPOE model succeeds with fewer fitting constants. The fitted PPOE reaction scheme leads to known folding behavior, whereas the IUP properties are contradicted by experiment. The appearance of a conflict with multipath theoretical models seems to be due to their different focus, namely on multitrack microscopic behavior versus cooperative macroscopic behavior. The integration of three well-documented principles in the PPOE model (cooperative foldons, sequential stabilization, optional errors) provides a unifying explanation for how proteins fold and why they fold in that way.

Denaturation of jack-bean urease by sodium n-dodecyl sulphate: A kinetic study below the critical micelle concentration

Kinetics of urease denaturation by anionic surfactant (sodium n-dodecyl sulphate, SDS) at concentrations below the critical micelle concentration (CMC) is investigated spectrophotometrically at neutral pH and the corresponding two-phase kinetic parameters of the process are estimated from a three-state reversible process using a binomial exponential relation based on the relaxation time method as: Using a prepared computer program, the experimental data are properly fitted into a binomial exponential relation, considering a two-phase denaturation pathway including a kinetically stable folded intermediate formed at SDS concentration of 1.1 mM. Forward and backward rate constants are estimated as: k(1)=0.2141+/-4.5 x 10(-3), k(2)=5.173 x 10(-3)+/-8.3 x 10(-5), k(-1)=0.09432+/-3.6 x 10(-4) and k(-2)=2.079 x 10(-3)+/-5.6 x 10(-5)s(-1) for the proposed mechanism. The rate-limiting step as well as the reaction coordinates in the denaturation mechanism are established. The mechanism involves formation of a kinetically stable folded native like intermediate through the electrostatic interactions. The intermediate was found to be more stable even than the native form (by about 9 kJmol(-1)) and still hexamer, because no loss of amplitude was observed. Electrophoresis experiments on the native and surfactant/urease complexes indicated a higher mobility for the kinetically folded native like intermediate.

Chemometric studies of lysozyme upon interaction with sodium dodecyl sulfate and β-cyclodextrin

The interaction of hen egg-white lysozyme with sodium n-dodecyl sulfate (SDS) as an anionic surfactant was investigated by UV-vis spectrophotometry at different pHs at 25 degrees C using HCl/glycine and NaOH/glycine for acidic and basic pH ranges, respectively. Analysis of the spectral data using chemometric method gave the evidence for the existence of intermediate components during the cited interaction. Results also indicated a connection between turbidity of the protein solution upon interaction with SDS and distribution of our newly found intermediates. As intermediates are important in aggregation of proteins, beta-cyclodextrin was employed as an anti-aggregation agent and the results obtained for the lysozyme-SDS-beta-cyclodextrin ternary system were compared with those obtained in the absence of beta-cyclodextrin on distribution and mole fraction of intermediates with. It is also shown that as the distribution of intermediates broadens in a range of SDS concentrations, the turbidity and aggregation state of solution are reduced.

The folding nucleus of the insulin superfamily: a flexible peptide model foreshadows the native state

Oxidative folding of insulin-like growth factor I (IGF-I) and single-chain insulin analogs proceeds via one- and two-disulfide intermediates. A predominant one-disulfide intermediate in each case contains the canonical A20-B19 disulfide bridge (cystines 18-61 in IGF-I and 19-85 in human proinsulin). Here, we describe a disulfide-linked peptide model of this on-pathway intermediate. One peptide fragment (19 amino acids) spans IGF-I residues 7-25 (canonical positions B8-B26 in the insulin superfamily); the other (18 amino acids) spans IGF-I residues 53-70 (positions A12-A21 and D1-D8). Containing only half of the IGF-I sequence, the disulfide-linked polypeptide (designated IGF-p) is not well ordered. Nascent helical elements corresponding to native alpha-helices are nonetheless observed at 4 degrees C. Furthermore, (13)C-edited nuclear Overhauser effects establish transient formation of a native-like partial core; no non-native nuclear Overhauser effects are observed. Together, these observations suggest that early events in the folding of insulin-related polypeptides are nucleated by a native-like molten subdomain containing Cys(A20) and Cys(B19). We propose that nascent interactions within this subdomain orient the A20 and B19 thiolates for disulfide bond formation and stabilize the one-disulfide intermediate once formed. Substitutions in the corresponding region of insulin are associated with inefficient chain combination and impaired biosynthetic expression. The intrinsic conformational propensities of a flexible disulfide-linked peptide thus define a folding nucleus, foreshadowing the structure of the native state.

Refolding studies using pressure: The folding landscape of lysozyme in the pressure–temperature plane

Refolding of hen egg white lysozyme after pressure unfolding was measured by FTIR spectroscopy. The high-pressure treatment was found to be useful for unfolding/refolding studies because pressure acts against aggregation, and therefore no irreversible aggregation takes place during the pressure treatment. After the release of the pressure, folding intermediate structures were found which were formed during the decompression of the lysozyme. These were aggregation prone when heated, as indicated by their lower stability against aggregation. The intermediates were only formed if the protein was unfolded, subdenaturing pressures could not populate these intermediates. We introduced the notion of a superfunnel to describe the free energy landscape of interacting polypeptide chains. This can explain the propensity of folding intermediates to aggregate. A possible Gibbs-free energy landscape for lysozyme was constructed for the whole pressure-temperature plane.

Probing Lysozyme Conformation with Light Reveals a New Folding Intermediate

A means to control lysozyme conformation with light illumination has been developed using the interaction of the protein with a photoresponsive surfactant. Upon exposure to the appropriate wavelength of light, the azobenzene surfactant undergoes a reversible photoisomerization, with the visible-light (trans) form being more hydrophobic than the UV-light (cis) form. As a result, surfactant binding to the protein and, thus, protein unfolding, can be tuned with light. Small-angle neutron scattering (SANS) measurements were used to provide detailed information of the protein conformation in solution. Shape-reconstruction methods applied to the SANS data indicate that under visible light the protein exhibits a native-like form at low surfactant concentrations, a partially swollen form at intermediate concentrations, and a swollen/unfolded form at higher surfactant concentrations. Furthermore, the SANS data combined with FT-IR spectroscopic analysis of the protein secondary structure reveal that unfolding occurs primarily in the alpha domain of lysozyme, while the beta domain remains relatively intact. Thus, the surfactant-unfolded intermediate of lysozyme appears to be a separate structure than the well-known alpha-domain intermediate of lysozyme that contains a folded alpha domain and unfolded beta domain. Because the interactions between the photosurfactant and protein can be tuned with light, illumination with UV light returns the protein to a native-like conformation. Fluorescence emission data of the nonpolar probe Nile red indicate that hydrophobic domains become available for probe partitioning in surfactant-protein solutions under visible light, while the availability of these hydrophobic domains to the probe decrease under UV light. Dynamic light scattering and UV-vis spectroscopic measurements further confirm the shape-reconstruction findings and reveal three discrete conformations of lysozyme. The results clearly demonstrate that visible light causes a greater degree of lysozyme swelling than UV light, thus allowing for the protein conformation to be controlled with light.

RNA and protein folding: common themes and variations

Visualizing the navigation of an ensemble of unfolded molecules through the bumpy energy landscape in search of the native state gives a pictorial view of biomolecular folding. This picture, when combined with concepts in polymer theory, provides a unified theory of RNA and protein folding. Just as for proteins, the major folding free energy barrier for RNA scales sublinearly with the number of nucleotides, which allows us to extract the elusive prefactor for RNA folding. Several folding scenarios can be anticipated by considering variations in the energy landscape that depend on sequence, native topology, and external conditions. RNA and protein folding mechanism can be described by the kinetic partitioning mechanism (KPM) according to which a fraction (Phi) of molecules reaches the native state directly, whereas the remaining fraction gets kinetically trapped in metastable conformations. For two-state folders Phi approximately 1. Molecular chaperones are recruited to assist protein folding whenever Phi is small. We show that the iterative annealing mechanism, introduced to describe chaperonin-mediated folding, can be generalized to understand protein-assisted RNA folding. The major differences between the folding of proteins and RNA arise in the early stages of folding. For RNA, folding can only begin after the polyelectrolyte problem is solved, whereas protein collapse requires burial of hydrophobic residues. Cross-fertilization of ideas between the two fields should lead to an understanding of how RNA and proteins solve their folding problems.

Mixed macromolecular crowding accelerates the oxidative refolding of reduced, denatured lysozyme: implications for protein folding in intracellular environments

The oxidative refolding of reduced, denatured hen egg white lysozyme in the presence of a mixed macromolecular crowding agent containing both bovine serum albumin (BSA) and polysaccharide has been studied from a physiological point of view. When the total concentration of the mixed crowding agent is 100 g/liter, in which the weight ratio of BSA to dextran 70 is 1:9, the refolding yield of lysozyme after refolding for 4 h under this condition increases 24% compared with that in the presence of BSA and 16% compared with dextran 70. A remarkable increase in the refolding yield of lysozyme by a mixed crowding agent containing BSA and Ficoll 70 is also observed. Further folding kinetics analyses show that these two mixed crowding agents accelerate the oxidative refolding of lysozyme remarkably, compared with single crowding agents. These results suggest that the stabilization effects of mixed macromolecular crowding agents are stronger than those of single polysaccharide crowding agents such as dextran 70 and Ficoll 70, whereas the excluded volume effects of mixed macromolecular crowding agents are weaker than those of single protein crowding agents such as BSA. Both the refolding yield and the rate of the oxidative refolding of lysozyme in these two mixed crowded solutions with suitable weight ratios are higher than those in single crowded solutions, indicating that mixed macromolecular crowding agents are more favorable to lysozyme folding and can be used to simulate the intracellular environments more accurately than single crowding agents do.

Integrating folding kinetics and protein function: biphasic kinetics and dual binding specificity in a WW domain

Because of the association of beta-sheet formation with the initiation and propagation of amyloid diseases, model systems have been sought to further our understanding of this process. WW domains have been proposed as one such model system. Whereas the folding of the WW domains from human Yes-associated protein (YAP) and Pin have been shown to obey single-exponential kinetics, the folding of the WW domain from formin-binding protein (FBP) 28 has been shown to proceed via biphasic kinetics. From an analysis of free-energy landscapes from atomic-level molecular dynamics simulations, the biphasic folding kinetics observed in the FBP WW domain may be traced to the ability of this WW domain to adopt two slightly different forms of packing in its hydrophobic core. This conformational change is propagated along the peptide backbone and affects the position of a tryptophan residue shown in other WW domains to play a key role in binding. The WW domains of Pin and YAP do not support more than one type of packing each, leading to monophasic folding kinetics. The ability of the FBP WW domain to assume two different types of packing may, in turn, explain the capacity of this WW domain to bind two classes of ligand, a property that is not shared by other WW domains. These findings lead to the hypothesis that lability with respect to conformations separated by an observable barrier as a requirement for function is incompatible with the ability of a protein to fold via single-exponential kinetics.

The Nature of the Rate-limiting Steps in the Refolding of the Cofactor-dependent Protein Aspartate Aminotransferase

The refolding of mitochondrial aspartate aminotransferase (mAAT; EC 2.6.1.1) has been studied following unfolding in 6 m guanidine hydrochloride for different periods of time. Whereas reactivation of equilibrium-unfolded mAAT is sigmoidal, reactivation of the short term unfolded protein displays a double exponential behavior consistent with the presence of fast and slow refolding species. The amplitude of the fast phase decreases with increasing unfolding times (k approximately 0.75 min(-1) at 20 degrees C) and becomes undetectable at equilibrium unfolding. According to hydrogen exchange and stopped-flow intrinsic fluorescence data, unfolding of mAAT appears to be complete in less than 10 s, but hydrolysis of the Schiff base linking the coenzyme pyridoxal 5'-phosphate (PLP) to the polypeptide is much slower (k approximately 0.08 min(-1)). This implies the existence in short term unfolded samples of unfolded species with PLP still attached. However, since the disappearance of the fast refolding phase is about 10-fold faster than the release of PLP, the fast refolding phase does not correspond to folding of the coenzyme-containing molecules. The fast refolding phase disappears more rapidly in the pyridoxamine and apoenzyme forms of mAAT, both of which lack covalently attached cofactor. Thus, bound PLP increases the kinetic stability of the fast refolding unfolding intermediates. Conversion between fast and slow folding forms also takes place in an early folding intermediate. The presence of cyclophilin has no effect on the reactivation of either equilibrium or short term unfolded mAAT. These results suggest that proline isomerization may not be the only factor determining the slow refolding of this cofactor-dependent protein.

Is there a unifying mechanism for protein folding?

Proteins appear to fold by diverse pathways, but variations of a simple mechanism - nucleation-condensation - describe the overall features of folding of most domains. In general, secondary structure is inherently unstable and its stability is enhanced by tertiary interactions. Consequently, an extensive interplay of secondary and tertiary interactions determines the transition-state for folding, which is structurally similar to the native state, being formed in a general collapse (condensation) around a diffuse nucleus. As the propensity for stable secondary structure increases, folding becomes more hierarchical and eventually follows a framework mechanism where the transition state is assembled from pre-formed secondary structural elements.

Equilibrium and kinetic folding of hen egg-white lysozyme under acidic conditions

The equilibrium and kinetic folding of hen egg-white lysozyme was studied by means of circular dichroism spectra in the far- and near-ultraviolet (UV) regions at 25 degrees C under the acidic pH conditions. In equilibrium condition at pH 2.2, hen lysozyme shows a single cooperative transition in the GdnCl-induced unfolding experiment. However, in the GdnCl-induced unfolding process at lower pH 0.9, a distinct intermediate state with molten globule characteristics was observed. The time-dependent unfolding and refolding of the protein were induced by concentration jumps of the denaturant and measured by using stopped-flow circular dichroism at pH 2.2. Immediately after the dilution of denaturant, the kinetics of refolding shows evidence of a major unresolved far-UV CD change during the dead time (<10 ms) of the stopped-flow experiment (burst phase). The observed refolding and unfolding curves were both fitted well to a single-exponential function, and the rate constants obtained in the far- and near-UV regions coincided with each other. The dependence on denaturant concentration of amplitudes of burst phase and both rate constants was modeled quantitatively by a sequential three-state mechanism, U<-->I<-->N, in which the burst-phase intermediate (I) in rapid equilibrium with the unfolded state (U) precedes the rate-determining formation of the native state (N). The role of folding intermediate state of hen lysozyme was discussed.

Fast folding of a four-helical bundle protein

The FK506-FKBP12 binding-domain of the kinase FRAP (FRB) forms a classic up-down four-helical bundle. The folding pathway of this protein has been investigated using a combination of equilibrium and kinetic studies. The native state of the protein is stable with respect to the unfolded state by some 7 kcal mol(-1) at pH 6.0, 10 degrees C. A kinetic analysis of unfolding and refolding rate constants as a function of chemical denaturant concentration suggests that an intermediate state may be populated during folding at low concentrations of denaturant. The presence of this intermediate state is confirmed by refolding experiments performed in the presence of the hydrophobic dye 8-anilinonaphthalene-1 sulfonate (ANS). ANS binds to the partially folded intermediate state populated during the folding of FRB and undergoes a large change in fluorescence that can be detected using stopped-flow techniques. Analysis of the kinetic data suggests that the intermediate state is compact and it may even be a misfolded species that has to partially unfold before it can reach the transition state. Folding and unfolding rate constants in water are approximately 150-200 s(-1) and 0.005-0.06 s(-1), respectively, at neutral pH and 10 degrees C. The folding of FRB is somewhat slower than for other all-helical proteins, probably as a consequence of the formation of a metastable intermediate state. The folding rate constant in the absence of any populated intermediate can be estimated to be 8800 s(-1). Despite the presence of an intermediate state, which effectively slows folding, the protein still folds rapidly with a half-life of 5 ms at 10 degrees C. The dependence of the rate constants on denaturant concentration indicates that the transition state for folding is compact with some 80% of the surface area exposed in the unfolded state buried in the transition state. Data presented for FRB is compared with kinetic data obtained for other all-helical proteins.

Fast-folding protein kinetics, hidden intermediates, and the sequential stabilization model

Do two-state proteins fold by pathways or funnels? Native-state hydrogen exchange experiments show discrete nonnative structures in equilibrium with the native state. These could be called hidden intermediates (HI) because their populations are small at equilibrium, and they are not detected in kinetic experiments. HIs have been invoked as disproof of funnel models, because funnel pictures appear to indicate (1) no specific sequences of events in folding; (2) a continuum, rather than a discrete ladder, of structures; and (3) smooth landscapes. In the present study, we solve the exact dynamics of a simple model. We find, instead, that the present microscopic model is indeed consistent with HIs and transition states, but such states occur in parallel, rather than along the single pathway predicted by the sequential stabilization model. At the microscopic level, we observe a huge multiplicity of trajectories. But at the macroscopic level, we observe two pathways of specific sequences of events that are relatively traditional except that they are in parallel, so there is not a single reaction coordinate. Using singular value decomposition, we show an accurate representation of the shapes of the model energy landscapes. They are highly complex funnels.

Exploring protein aggregation and self-propagation using lattice models: phase diagram and kinetics

Many seemingly unrelated neurodegenerative disorders, such as amyloid and prion diseases, are associated with propagating fibrils whose structures are dramatically different from the native states of the corresponding monomers. This observation, along with the experimental demonstration that any protein can aggregate to form either fibrils or amorphous structures (inclusion bodies) under appropriate external conditions, suggest that there must be general principles that govern aggregation mechanisms. To probe generic aspects of prion-like behavior we use the model of Harrison, Chan, Prusiner, and Cohen. In this model, aggregation of a structure, that is conformationally distinct from the native state of the monomer, occurs by three parallel routes. Kinetic partitioning, which leads to parallel assembly pathways, occurs early in the aggregation process. In all pathways transient unfolding precedes polymerization and self-propagation. Chain polymerization is consistent with templated assembly, with the dimer being the minimal nucleus. The kinetic effciency of R(n-1) + G --> R(n) (R is the aggregation prone state and G is either U, the unfolded state, or N, the native state of the monomer) is increased when polymerization occurs in the presence of a "seed" (a dimer). These results support the seeded nucleated-polymerization model of fibril formation in amyloid peptides. To probe generic aspects of aggregation in two-state proteins, we use lattice models with side chains. The phase diagram in the (T,C) plane (T is the temperature and C is the polypeptide concentration) reveals a bewildering array of "phases" or structures. Explicit computations for dimers show that there are at least six phases including ordered structures and amorphous aggregates. In the ordered region of the phase diagram there are three distinct structures. We find ordered dimers (OD) in which each monomer is in the folded state and the interaction between the monomers occurs via a well-defined interface. In the domain-swapped structures a certain fraction of intrachain contacts are replaced by interchain contacts. In the parallel dimers the interface is stabilized by favorable intermolecular hydrophobic interactions. The kinetics of folding to OD shows that aggregation proceeds directly from U in a dynamically cooperative manner without populating partially structured intermediates. These results support the experimental observation that ordered aggregation in the two-state folders U1A and CI2 takes place from U. The contrasting aggregation processes in the two models suggest that there are several distinct mechanisms for polymerization that depend not only on the polypeptide sequence but also on external conditions (such as C, T, pH, and salt concentration).

Structural transitions in neutral and charged proteins in vacuo

In vacuo proteins provide a simple laboratory to explore the roles of sequence, temperature, charge state, and initial configuration in protein folding. Moreover, by the very absence of solvent, the study of anhydrous proteins in vacuo will also help us to understand specific environmental effects. From the experimental viewpoint, these systems are now beginning to be characterized at low resolution. Molecular dynamics (MD) simulations, in combination with tools for protein shape analysis, can complement experiments and provide further insights on the folding-unfolding transitions of these proteins. We review some aspects of this issue by using the results from a detailed MD study of hen egg-white lysozyme. For lysozyme ions, unfolding can be triggered by Coulombic repulsion. In neutral lysozyme, unfolding can be induced by centrifugal forces and also by weakening the monomer-monomer interaction. In both cases, the resulting unfolded transients can be used as initial configurations for relaxation dynamics. All trajectories are analyzed in terms of global molecular shape features of the backbone, including its anisometry and chain entanglement complexity. This strategy allows us to quantify separately the degree of polymer collapse and the evolution of large-scale folding features. Using these last two notions, we discuss some basic questions regarding the nature of the accessible paths associated with unfolding from, and refolding into, compact conformers.

Population analyses of kinetic partitioning in protein folding

A simple lattice model of protein folding is studied in order to analyze the kinetic partitioning phenomena in the energy landscape perspective. By restricting the area of conformational space, it becomes possible to follow many Monte Carlo trajectories until they reach equilibrium. Alteration of population of trajectories is monitored and the relations between the energy landscape and kinetics are examined. Kinetic partitioning phenomena are categorized into different types in terms of characteristic time constants and partitioning ratio. In a specific partitioning process, refolding proceeds along the parallel pathways; the time constants have a temperature dependence similar to that observed in hen lysozyme. High-energy conformations are classified into groups according to the probability that the trajectories starting from those conformations will reach each energy valley. The partitioning ratio is determined by the way in which the conformational space is organized into these groups.

Organic cosolvents and hen egg white lysozyme folding

Studies on the influence of organic cosolvents on lysozyme folding have been reported. As most of the researches are confined to a few specific molecules and focus on equilibrium states, less is known about the effect on folding dynamics. We have studied the influence of six soluble organic cosolvents on hen egg white lysozyme heat induced denaturation and refolding dynamics. It was found that trifluoroethanol (TFE) can change the folding pathway significantly. With the presence of TFE, the overshot phenomenon generally observed in lysozyme folding at 222 nm disappears. The common mechanism of how organic cosolvents influence folding is analyzed. The heat induced denaturation temperature was found to have a quantitative relationship with the slow phase rate constant during folding. We discuss this finding and hypothesize that it is due to the similar influence of organic cosolvent on the transition state of heat denaturation and refolding.

Macromolecular crowding perturbs protein refolding kinetics: implications for folding inside the cell

We have studied the effects of macromolecular crowding on protein folding kinetics by studying the oxidative refolding of hen lysozyme in the absence and presence of high concentrations of bovine serum albumin and Ficoll 70. The heterogeneity characteristic of the lysozyme refolding process is preserved under crowded conditions. This, together with the observation that the refolding intermediates that accumulate to significant levels are very similar in the absence and presence of Ficoll, suggests that crowding does not alter substantially the energetics of the protein folding reaction. However, the presence of high concentrations of macromolecules results in the acceleration of the fast track of the refolding process whereas the slow track is substantially retarded. The results can be explained by preferential excluded volume stabilization of compact states relative to more unfolded states, and suggest that, relative to dilute solutions, the rates of many protein folding processes are likely to be altered under conditions that more closely resemble the intracellular environment.

Understanding protein folding via free-energy surfaces from theory and experiment

The ability of protein molecules to fold into their highly structured functional states is one of the most remarkable evolutionary achievements of biology. In recent years, our understanding of the way in which this complex self-assembly process takes place has increased dramatically. Much of the reason for this advance has been the development of energy surfaces (landscapes), which allow the folding reaction to be described and visualized in a meaningful manner. Analysis of these surfaces, derived from the constructive interplay between theory and experiment, has led to the development of a unified mechanism for folding and a recognition of the underlying factors that control the rates and products of the folding process.

Thermal unfolding of an intermediate is associated with non-Arrhenius kinetics in the folding of hen lysozyme

A variety of techniques, including quenched-flow hydrogen exchange labelling monitored by electrospray ionization mass spectrometry, and stopped-flow absorbance, fluorescence and circular dichroism spectroscopy, has been used to investigate the refolding kinetics of hen lysozyme over a temperature range from 2 degrees C to 50 degrees C. Simple Arrhenius behaviour is not observed, and although the overall rate of folding increases from 2 to 40 degrees C, it decreases above 40 degrees C. In addition, the transient intermediate on the major folding pathway at 20 degrees C, in which the alpha-domain is persistently structured in the absence of a stable beta-domain, is thermally unfolded in a sigmoidal transition (T(m) approximately 40 degrees C) indicative of a cooperatively folded state. At all temperatures, however, there is evidence for fast ( approximately 25 %) and slow ( approximately 75 %) populations of refolding molecules. By using transition state theory, the kinetic data from various experiments were jointly fitted to a sequential three-state model for the slow folding pathway. Together with previous findings, these results indicate that the alpha-domain intermediate is a productive species on the folding route between the denatured and native states, and which accumulates as a consequence of its intrinsic stability. Our analysis suggests that the temperature dependence of the rate constant for lysozyme folding depends on both the total change in the heat capacity between the ground and transition states (the dominant factor at low temperatures) and the heat-induced destabilization of the alpha-domain intermediate (the dominant factor at high temperatures). Destabilization of such kinetically competent intermediate species is likely to be a determining factor in the non-Arrhenius temperature dependence of the folding rate of those proteins for which one or more intermediates are populated.

Mechanism of the antichaperone activity of protein disulfide isomerase: facilitated assembly of large, insoluble aggregates of denatured lysozyme and PDI

Protein disulfide isomerase (PDI), a folding catalyst and chaperone can, under certain conditions, facilitate the misfolding and aggregation of its substrates. This behavior, termed antichaperone activity [Puig, A., and Gilbert, H. F., (1994) J. Biol. Chem. 269, 25889] may provide a common mechanism for aggregate formation in the cell, both as a normal consequence of cell function or as a consequence of disease. When diluted from the denaturant, reduced, denatured lysozyme (10-50 microM) remains soluble, although it does aggregate to form an ensemble of species with an average sedimentation coefficient of 23 +/- 5 S (approximately 600 +/- 100 kDa). When low concentrations of PDI (1-5 microM) are present, the majority (80 +/- 8%) of lysozyme molecules precipitate in large, insoluble aggregates, together with 87 +/- 12% of the PDI. PDI-facilitated aggregation occurs even when disulfide formation is precluded by the presence of dithiothreitol (10 mM). Maximal lysozyme-PDI precipitation occurs at a constant lysozyme/PDI ratio of 10:1 over a range of lysozyme concentrations (10-50 microM). Concomitant resolubilization of PDI and lysozyme from these aggregates by increasing concentrations of urea suggests that PDI is an integral component of the mixed aggregate. PDI induces lysozyme aggregation by noncovalently cross-linking 23 S lysozyme species to form aggregates that become so large (approximately 38,000 S) that they are cleared from the analytical ultracentrifuge even at low speed (1500 rpm). The rate of insoluble aggregate formation increases with increasing PDI concentration (although a threshold PDI concentration is observed). However, increasing lysozyme concentration slows the rate of aggregation, presumably by depleting PDI from solution. A simple mechanism is proposed that accounts for these unusual aggregation kinetics as well as the switch between antichaperone and chaperone behavior observed at higher concentrations of PDI.