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

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.

Accurate prediction of protein folding mechanisms by simple structure-based statistical mechanical models

Recent breakthroughs in highly accurate protein structure prediction using deep neural networks have made considerable progress in solving the structure prediction component of the 'protein folding problem'. However, predicting detailed mechanisms of how proteins fold into specific native structures remains challenging, especially for multidomain proteins constituting most of the proteomes. Here, we develop a simple structure-based statistical mechanical model that introduces nonlocal interactions driving the folding of multidomain proteins. Our model successfully predicts protein folding processes consistent with experiments, without the limitations of protein size and shape. Furthermore, slight modifications of the model allow prediction of disulfide-oxidative and disulfide-intact protein folding. These predictions depict details of the folding processes beyond reproducing experimental results and provide a rationale for the folding mechanisms. Thus, our physics-based models enable accurate prediction of protein folding mechanisms with low computational complexity, paving the way for solving the folding process component of the 'protein folding problem'.

The Right-Handed Parallel β-Helix Topology of Erwinia chrysanthemi Pectin Methylesterase Is Intimately Associated with Both Sequential Folding and Resistance to High Pressure

The complex topologies of large multi-domain globular proteins make the study of their folding and assembly particularly demanding. It is often characterized by complex kinetics and undesired side reactions, such as aggregation. The structural simplicity of tandem-repeat proteins, which are characterized by the repetition of a basic structural motif and are stabilized exclusively by sequentially localized contacts, has provided opportunities for dissecting their folding landscapes. In this study, we focus on the Erwinia chrysanthemi pectin methylesterase (342 residues), an all-β pectinolytic enzyme with a right-handed parallel β-helix structure. Chemicals and pressure were chosen as denaturants and a variety of optical techniques were used in conjunction with stopped-flow equipment to investigate the folding mechanism of the enzyme at 25 °C. Under equilibrium conditions, both chemical- and pressure-induced unfolding show two-state transitions, with average conformational stability (ΔG° = 35 ± 5 kJ·mol⁻¹) but exceptionally high resistance to pressure (P_m = 800 ± 7 MPa). Stopped-flow kinetic experiments revealed a very rapid (τ < 1 ms) hydrophobic collapse accompanied by the formation of an extended secondary structure but did not reveal stable tertiary contacts. This is followed by three distinct cooperative phases and the significant population of two intermediate species. The kinetics followed by intrinsic fluorescence shows a lag phase, strongly indicating that these intermediates are productive species on a sequential folding pathway, for which we propose a plausible model. These combined data demonstrate that even a large repeat protein can fold in a highly cooperative manner.

Major conformational changes in the structure of lysozyme obtained from a crystal with a very low solvent content

A new crystal form of lysozyme with a very low solvent content (26.35%) has been obtained in the orthorhombic space group P212121 (with unit-cell parameters a = 30.04, b = 51.68, c = 61.53 Å). The lysozyme structure obtained from these crystals does not show the typical overall fold. Instead, major conformational changes take place in some elements of the secondary structure and in the hydrophobic core of the protein. At the end of the central α-helix (α2), Glu35 is usually buried in the catalytic site and shows an abnormally high pKa value, which is key to the activity of the enzyme. The high pKa value of this glutamate residue is favoured by the hydrophobic environment, particularly by its neighbour Trp108, which is important for structural stability and saccharide binding. In this new structure, Trp108 shows a 90° rotation of its side chain, which results in the rearrangement of the hydrophobic core. Conformational changes also result in the exposure of Glu35 to the solvent, which impairs the catalytic site by increasing the distance between Glu35 and Asp52 and lowering the pKa value of the glutamate. Altogether, this new lysozyme structure reveals major conformational changes in the hydrophobic core and catalytic site that might play a role in the folding and bactericidal function of the protein.

Helix-strand interaction regulates stability and aggregation of the human mitochondrial membrane protein channel VDAC3

Human mitochondrial VDACs bind amyloidogenic proteins, but do not intrinsically aggregate. Gupta and Mahalakshmi find that an interaction between the N-terminal α-helix and strands β7–β9 regulates VDAC aggregation and stability, providing a plausible mechanism for VDAC coaggregation in cells.

Voltage-dependent anion channels (VDACs) are β-sheet–rich transmembrane β-barrels that are vital for metabolite transport across the mitochondrial membrane. Under cellular stress, human VDACs hetero-oligomerize and coaggregate with proteins that can form amyloidogenic and neurodegenerative deposits, implicating a role for VDACs in proteotoxicity. However, whether VDACs possess intrinsic interaction sites that can lead to protein aggregation is not known. Here, we couple a systematic thiol replacement strategy with far-UV circular dichroism spectropolarimetry and UV scattering spectroscopy to map aggregation-prone regions of human VDACs, using isoform 3 as our model VDAC. We show that the region comprising strands β7–β9 is highly aggregation prone. Further, we find that an α1–β7–β9 interaction (involving the hVDAC3 N-terminal α1 helix) can lower protein aggregation, whereas perturbations of this interaction promote VDAC aggregation. We also show that hVDAC3 aggregation proceeds via a partially unfolded structure. Our findings allow us to propose a plausible mechanism for the role of human VDACs in forming proteotoxic aggregates in the cell. The key target sites on VDACs—strands β7–β9—may be useful for developing VDAC aggregation inhibitors.

Rutin attenuates negatively charged surfactant (SDS)-induced lysozyme aggregation/amyloid formation and its cytotoxicity

Amyloid fibrils are highly ordered protein assemblies known to contribute to the pathology of a variety of genetic and aging-associated diseases. Here, we have investigated the aggregation propensity of lysozyme in the presence of a negatively charged surfactant (SDS) and evaluated the anti-aggregation activity of rutin. Multiple approaches such as turbidity measurements, dye binding assays, intrinsic fluorescence, circular dichroism (CD), transmission electron microscopy (TEM), MTT and comet assays have been used for this purpose. We inferred that SDS induces aggregation of lysozyme in 0.2-0.6 mM concentration range while at higher concentration range (0.8-1.0 mM), it leads to solubilization/stabilization of protein. Intrinsic/extrinsic fluorescence and CD analysis confirmed significant conformational changes in lysozyme at 0.2 mM SDS. Thioflavin T (ThT), congo red binding and TEM analysis further reaffirmed the formation of lysozyme fibrils. Moreover, MTT assay demonstrated cytotoxicity of these fibrils towards neuroblastoma cell lines (SH-SY5Y) and their attenuation by rutin. Comet assay supported the cytotoxicity mechanism via DNA damage. Molecular docking results also advocate a strong interaction between lysozyme and rutin. The current study indicates a mechanistic approach assuming structural constraints and specific aromatic interactions of rutin with HEWL aggregates.

Glutathionylation of the Bacterial Hsp70 Chaperone DnaK Provides a Link between Oxidative Stress and the Heat Shock Response

DnaK is the major bacterial Hsp70, participating in DNA replication, protein folding, and the stress response. DnaK cooperates with the Hsp40 co-chaperone DnaJ and the nucleotide exchange factor GrpE. Under non-stress conditions, DnaK binds to the heat shock transcription factor σ(32)and facilitates its degradation. Oxidative stress results in temporary inactivation of DnaK due to depletion of cellular ATP and thiol modifications such as glutathionylation until normal cellular ATP levels and a reducing environment are restored. However, the biological significance of DnaK glutathionylation remains unknown, and the mechanisms by which glutathionylation may regulate the activity of DnaK are also unclear. We investigated the conditions under which Escherichia coli DnaK undergoesS-glutathionylation. We observed glutathionylation of DnaK in lysates of E. coli cells that had been subjected to oxidative stress. We also obtained homogeneously glutathionylated DnaK using purified DnaK in the apo state. We found that glutathionylation of DnaK reversibly changes the secondary structure and tertiary conformation, leading to reduced nucleotide and peptide binding ability. The chaperone activity of DnaK was reversibly down-regulated by glutathionylation, accompanying the structural changes. We found that interaction of DnaK with DnaJ, GrpE, or σ(32)becomes weaker when DnaK is glutathionylated, and the interaction is restored upon deglutathionylation. This study confirms that glutathionylation down-regulates the functions of DnaK under oxidizing conditions, and this down-regulation may facilitate release of σ(32)from its interaction with DnaK, thus triggering the heat shock response. Such a mechanism provides a link between oxidative stress and the heat shock response in bacteria.

Improving Convergence of Replica-Exchange Simulations through Coupling to a High-Temperature Structure Reservoir

Parallel tempering or replica-exchange molecular dynamics (REMD) significantly increases the efficiency of conformational sampling for complex molecular systems. However, obtaining converged data with REMD remains challenging, especially for large systems with complex topologies. We propose a new variant to REMD where the replicas are also permitted to exchange with an ensemble of structures that have been generated in advance using high-temperature MD simulations, similar in spirit to J-walking methods. We tested this approach on two non-trivial model systems, a β-hairpin and a 3-stranded β-sheet and compared the results to those obtained from very long (>100 ns) standard REMD simulations. The resulting ensembles were indistinguishable, including relative populations of different conformations on the unfolded state. The use of the reservoir is shown to significantly reduce the time required for convergence.

Adaptive Steered Molecular Dynamics of the Long-Distance Unfolding of Neuropeptide Y

Neuropeptide Y (NPY) has been found to adopt two stable conformations in vivo: (1) a monomeric form called the PP-fold in which a polyproline tail is folded onto an α-helix via a β-turn and (2) a dimeric form of the unfolded proteins in which the α-helices interact with each other via side chains. The transition pathway and rates between the two conformations remain unknown and are important to the nature of the binding of the protein. Toward addressing this question, the present work suggests that the unfolding of the PP-fold is too slow to play a role in NPY monomeric binding unless the receptor catalyzes it to do so. Specifically, the dynamics and structural changes of the unfolding of a monomeric NPY protein have been investigated in this work. Temperature accelerated molecular dynamics (MD) simulations at 500 K under constant (N,V,E) conditions suggests a hinge-like unraveling of the tail rather than a random unfolding. The free energetics of the proposed unfolding pathway have been described using an adaptive steered MD (SMD) approach at various temperatures. This approach generalizes the use of Jarzynski's equality through a series of stages that allows for better convergence along nonlinear and long-distance pathways. Results acquired using this approach provide a potential of mean force (PMF) with narrower error bars and are consistent with some of the earlier reports on the qualitative behavior of NPY binding.

GroEL/ES chaperonin modulates the mechanism and accelerates the rate of TIM-barrel domain folding

The GroEL/ES chaperonin system functions as a protein folding cage. Many obligate substrates of GroEL share the (βα)8 TIM-barrel fold, but how the chaperonin promotes folding of these proteins is not known. Here, we analyzed the folding of DapA at peptide resolution using hydrogen/deuterium exchange and mass spectrometry. During spontaneous folding, all elements of the DapA TIM barrel acquire structure simultaneously in a process associated with a long search time. In contrast, GroEL/ES accelerates folding more than 30-fold by catalyzing segmental structure formation in the TIM barrel. Segmental structure formation is also observed during the fast spontaneous folding of a structural homolog of DapA from a bacterium that lacks GroEL/ES. Thus, chaperonin independence correlates with folding properties otherwise enforced by protein confinement in the GroEL/ES cage. We suggest that folding catalysis by GroEL/ES is required by a set of proteins to reach native state at a biologically relevant timescale, avoiding aggregation or degradation.

Non-Arrhenius protein aggregation

Protein aggregation presents one of the key challenges in the development of protein biotherapeutics. It affects not only product quality but also potentially impacts safety, as protein aggregates have been shown to be linked with cytotoxicity and patient immunogenicity. Therefore, investigations of protein aggregation remain a major focus in pharmaceutical companies and academic institutions. Due to the complexity of the aggregation process and temperature-dependent conformational stability, temperature-induced protein aggregation is often non-Arrhenius over even relatively small temperature windows relevant for product development, and this makes low-temperature extrapolation difficult based simply on accelerated stability studies at high temperatures. This review discusses the non-Arrhenius nature of the temperature dependence of protein aggregation, explores possible causes, and considers inherent hurdles for accurately extrapolating aggregation rates from conventional industrial approaches for selecting accelerated conditions and from conventional or more advanced methods of analyzing the resulting rate data.

Sequential four-state folding/unfolding of goat α-lactalbumin and its N-terminal variants

Equilibria and kinetics of folding/unfolding of α-lactalbumin and its two N-terminal variants were studied by circular dichroism spectroscopy. The two variants were wild-type recombinant and Glu1-deletion (E1M) variants expressed in Escherichia coli. The presence of an extra methionine at the N terminus in recombinant α-lactalbumin destabilized the protein by 2 kcal/mol, while the stability was recovered in the E1M variant in which Glu1 was replaced by Met1. Kinetic folding/unfolding reactions of the proteins, induced by stopped-flow concentration jumps of guanidine hydrochloride, indicated the presence of a burst-phase in refolding, and gave chevron plots with significant curvatures in both the folding and unfolding limbs. The folding-limb curvature was interpreted in terms of accumulation of the burst-phase intermediate. However, there was no burst phase observed in the unfolding kinetics to interpret the unfolding-limb curvature. We thus assumed a sequential four-state mechanism, in which the folding from the burst-phase intermediate takes place via two transition states separated by a high-energy intermediate. We estimated changes in the free energies of the burst-phase intermediate and two transition states, caused by the N-terminal variations and also by the presence of stabilizing calcium ions. The Φ values at the N terminus and at the Ca(2+)-binding site thus obtained increased successively during folding, demonstrating the validity of the sequential mechanism. The stability and the folding behavior of the E1M variant were essentially identical to those of the authentic protein, allowing us to use this variant as a pseudo-wild-type α-lactalbumin in future studies.

Direct quantification of the attempt frequency determining the mechanical unfolding of ubiquitin protein

Understanding protein dynamics requires a comprehensive knowledge of the underlying potential energy surface that governs the motion of each individual protein molecule. Single molecule mechanical studies have provided the unprecedented opportunity to study the individual unfolding pathways along a well defined coordinate, the end-to-end length of the protein. In these experiments, unfolding requires surmounting an energy barrier that separates the native from the extended state. The calculation of the absolute value of the barrier height has traditionally relied on the assumption of an attempt frequency, υ(‡). Here we used single molecule force-clamp spectroscopy to directly determine the value of υ(‡) for mechanical unfolding by measuring the unfolding rate of the small protein ubiquitin at varying temperatures. Our experiments demonstrate a significant effect of the temperature on the mechanical rate of unfolding. By extrapolating the unfolding rate in the absence of force for different temperatures, varying within the range spanning from 5 to 45 °C, we measured a value for the activation barrier of ΔG(‡) = 71 ± 5 kJ/mol and an exponential prefactor υ(‡) ∼4 × 10(9) s(-1). Although the measured prefactor value is 3 orders of magnitude smaller than the value predicted by the transition state theory (∼6 × 10(12) s(-1)), it is 400-fold higher than that encountered in analogous experiments studying the effect of temperature on the reactivity of a protein-embedded disulfide bond (∼10(7) M(-1) s(-1)). This approach will allow quantitative characterization of the complete energy landscape of a folding polypeptide from highly extended states, of capital importance for proteins with elastic function.

Existence of metastable intermediate lysozyme conformation highlights the role of alcohols in altering protein stability

Alcohols have a manifold effect on the conformational and thermodynamic stability of native proteins. Here, we study the effect of moderate concentrations of trifluoroethanol (TFE) on the thermal stability of hen egg-white lysozyme (HEWL), by far-UV circular dichroism and by steady-state and time-resolved photoluminescence of intrinsic tryptophans. Our results highlight that TFE affects lysozyme stability by preferential solvation of the protein molecule. Furthermore, we discovered the existence at 20% TFE of an equilibrium partially folded state of lysozyme, intermediate between the native and the unfolded state. A three-state model is therefore used to interpolate the thermal denaturation data. Our analysis explains how the stabilization of the intermediate conformation enhances the entropic contribution to unfolding, and thus decreases the unfolding temperature, while, at the same time, TFE enhances the conformational stability of the native fold at room temperature. Eventually, we challenged the ability of these intermediate structures to form supramolecular aggregates by heating experiments at different TFE concentrations.

Evaluation of a non-Arrhenius model for therapeutic monoclonal antibody aggregation

Understanding antibody aggregation is of great significance for the pharmaceutical industry. We studied the aggregation of five different therapeutic monoclonal antibodies (mAbs) with size-exclusion chromatography-high-performance liquid chromatography (SEC-HPLC), fluorescence spectroscopy, electron microscopy, and light scattering methods at various temperatures with the aim of gaining insight into the aggregation process and developing models of it. In particular, we find that the kinetics can be described by a second-order model and are non-Arrhenius. Thus, we develop a non-Arrhenius model to connect accelerated aggregation experiments at high temperature to long-term storage experiments at low temperature. We evaluate our model by predicting mAb aggregation and comparing it with long-term behavior. Our results suggest that the number of monomers and mAb conformations within aggregates vary with the size and age of the aggregates, and that only certain sizes of aggregates are populated in the solution. We also propose a kinetic model based on conformational changes of proteins and monomer peak loss kinetics from SEC-HPLC. This model could be employed for a detail analysis of mAb aggregation kinetics.

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)).

Backbone 1H, 13C, and 15N resonance assignments for lysozyme from bacteriophage lambda

Lysozyme from lambda bacteriophage (lambda lysozyme) is an 18 kDa globular protein displaying some of the structural features common to all lysozymes; in particular, lambda lysozyme consists of two structural domains connected by a helix, and has its catalytic residues located at the interface between these two domains. An interesting feature of lambda lysozyme, when compared to the well-characterised hen egg-white lysozyme, is its lack of disulfide bridges; this makes lambda lysozyme an interesting system for studies of protein folding. A comparison of the folding properties of lambda lysozyme and hen lysozyme will provide important insights into the role that disulfide bonds play in the refolding pathway of the latter protein. Here we report the (1)H, (13)C and (15)N backbone resonance assignments for lambda lysozyme by heteronuclear multidimensional NMR spectroscopy. These assignments provide the starting point for detailed investigation of the refolding pathway using pulse-labelling hydrogen/deuterium exchange experiments monitored by NMR.

Comparison of the thermodynamic landscapes of unfolding and formation of the energy dissipative state in the isolated light harvesting complex II

In biochemistry and cell biology, understanding the molecular mechanisms by which physiological processes are regulated is regarded as an ultimate goal. In higher plants, one of the most widely investigated regulatory processes occurs in the light harvesting complexes (LHCII) of the chloroplast thylakoid membranes. Under limiting photon flux densities, LHCII harvests sunlight with high efficiency. When the intensity of incident radiation reaches levels close to the saturation of the photosynthesis, the efficiency of light harvesting is decreased by a process referred to as nonphotochemical quenching (NPQ), which enhances the singlet-excited state deactivation via nonradiative dissipative processes. Conformational rearrangements in LHCII are known to be crucial in promoting and controlling NPQ in vitro and in vivo. In this article, we address the thermodynamic nature of the conformational rearrangements promoting and controlling NPQ in isolated LHCII. A combined, linear reaction scheme in which the folded, quenched state represents a stable intermediate on the unfolding pathway was employed to describe the temperature dependence of the spectroscopic signatures associated with the chlorophyll fluorescence quenching and the loss of secondary structure motifs in LHCII. The thermodynamic model requires considering the temperature dependence of Gibbs free energy difference between the quenched and the unquenched states, as well as the unfolded and quenched states, of LHCII. Even though the same reaction scheme is adequate to describe the quenching and the unfolding processes in LHCII monomers and trimers, their thermodynamic characteristics were found to be markedly different. The results of the thermodynamic analysis shed light on the physiological importance of the trimeric state of LHCII in stabilizing the efficient light harvesting mode as well as preventing the quenched conformation of the protein from unfolding. Moreover, the transition to the quenched conformation in trimers reveals a larger degree of cooperativity than in monomers, explained by a small characteristic entropy (DeltaH(q) = 85 +/- 3 kJ mol(-1) compared to 125 +/- 5 kJ mol(-1) in monomers), which enables the fine-tuning of nonphotochemical quenching in vivo.

Hydrogen/deuterium exchange mass spectrometry

Amide hydrogen/deuterium (H/D) exchange of proteins monitored by mass spectrometry has established itself as a powerful method for probing protein conformational dynamics and protein interactions. The method uses isotope labeling to probe the rate at which protein backbone amide hydrogens undergo exchange. Backbone amide hydrogen exchange rates are particularly sensitive to hydrogen bonding; hydrogen bonding slows the exchange rates dramatically. Exchange rates reflect on the conformational mobility, hydrogen bonding strength, and solvent accessibility in protein structure. Mass spectrometric techniques are used to monitor the exchange events as mass shifts that arise through the incorporation of deuterium into the protein. Global conformational information can be deduced by monitoring the exchange profiles over time. Combining the labeling experiment with proteolysis under conditions that preserve the exchange information allows for localizing exchange events to distinct regions of the protein backbone and thus, the study of protein conformation with medium spatial resolution. Over the past decade, H/D exchange mass spectrometry has evolved into a versatile technique for investigating conformational dynamics and interactions in proteins, protein-ligand and protein-protein complexes.

Distinct unfolding and refolding pathways of ribonuclease a revealed by heating and cooling temperature jumps

Heating and cooling temperature jumps (T-jumps) were performed using a newly developed technique to trigger unfolding and refolding of wild-type ribonuclease A and a tryptophan-containing variant (Y115W). From the linear Arrhenius plots of the microscopic folding and unfolding rate constants, activation enthalpy (DeltaH(#)), and activation entropy (DeltaS(#)) were determined to characterize the kinetic transition states (TS) for the unfolding and refolding reactions. The single TS of the wild-type protein was split into three for the Y115W variant. Two of these transition states, TS1 and TS2, characterize a slow kinetic phase, and one, TS3, a fast phase. Heating T-jumps induced protein unfolding via TS2 and TS3; cooling T-jumps induced refolding via TS1 and TS3. The observed speed of the fast phase increased at lower temperature, due to a strongly negative DeltaH(#) of the folding-rate constant. The results are consistent with a path-dependent protein folding/unfolding mechanism. TS1 and TS2 are likely to reflect X-Pro(114) isomerization in the folded and unfolded protein, respectively, and TS3 the local conformational change of the beta-hairpin comprising Trp(115). A very fast protein folding/unfolding phase appears to precede both processes. The path dependence of the observed kinetics is suggestive of a rugged energy protein folding funnel.

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.

The first step of hen egg white lysozyme fibrillation, irreversible partial unfolding, is a two-state transition

Amyloid fibril depositions are associated with many neurodegenerative diseases as well as amyloidosis. The detailed molecular mechanism of fibrillation is still far from complete understanding. In our previous study of in vitro fibrillation of hen egg white lysozyme, an irreversible partially unfolded intermediate was characterized. A similarity of unfolding kinetics found for the secondary and tertiary structure of lysozyme using deep UV resonance Raman (DUVRR) and tryptophan fluorescence spectroscopy leads to a hypothesis that the unfolding might be a two-state transition. In this study, chemometric analysis, including abstract factor analysis (AFA), target factor analysis (TFA), evolving factor analysis (EFA), multivariate curve resolution-alternating least squares (ALS), and genetic algorithm, was employed to verify that only two principal components contribute to the DUVRR and fluorescence spectra of soluble fraction of lysozyme during the fibrillation process. However, a definite conclusion on the number of conformers cannot be made based solely on the above spectroscopic data although chemometric analysis suggested the existence of two principal components. Therefore, electrospray ionization mass spectrometry (ESI-MS) was also utilized to address the hypothesis. The protein ion charge state distribution (CSD) envelopes of the incubated lysozyme were well fitted with two principal components. Based on the above analysis, the partial unfolding of lysozyme during in vitro fibrillation was characterized quantitatively and proven to be a two-state transition. The combination of ESI-MS and Raman and fluorescence spectroscopies with advanced statistical analysis was demonstrated to be a powerful methodology for studying protein structural transformations.

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.

Aqueous two-phase system formation kinetics for elastin-like polypeptides of varying chain length

The kinetics of aqueous two-phase system (ATPS) formation for elastin-like polypeptides (ELP) with defined chemical composition and chain length was investigated by dark field microscopy in an on-chip format with a linear temperature gradient. Scattering intensities from peptide solutions in the presence and absence of sodium dodecyl sulfate (SDS) were recorded as a function of temperature and time, simultaneously. It was found that the formation of the ATPS for three ELPs of different molecular weights (36 075, 59 422, and 129 856 Da) in the absence of SDS followed a coalescence mechanism, and the rate constant and activation energy were independent of chain length. With the introduction of SDS into the ELP solutions, the rate constants were attenuated more strongly with increasing chain length. Moreover, the coalescence process in the presence of SDS showed non-Arrhenius kinetics as a function of temperature. For the two shorter ELPs, ATPS formation occurred via coalescence at all SDS concentrations and temperatures investigated. On the other hand, the coalescence process was greatly suppressed for the longest ELP at elevated temperatures and higher SDS concentrations. Under these circumstances, ATPS formation was forced to proceed via a mixed Ostwald ripening and coalescence mechanism.

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.

Sorbitol prevents the self-aggregation of unfolded lysozyme leading to and up to 13 degrees C stabilisation of the folded form

We present a calorimetric investigation of stabilisation of hen egg-white lysozyme with sorbitol in the pH range 3.8-10.5. Differential scanning calorimetry and steady-state fluorescence were used to determine the denaturation temperatures of lysozyme as a function of sorbitol concentration. The fluorescence data were collected in the presence of 2M urea to lower the melting point of the protein to an observable range of the instrument. The effect of sorbitol on the activation energy of unfolding was investigated by scanrate studies. The effect of sorbitol lysozyme interaction was investigated using isothermal titration calorimetry. The titration experiments were performed with folded as well as unfolded lysozyme to investigate in more detail the nature of the interaction. The data obtained in those experiments show a remarkable stabilisation effect of sorbitol. We observed a 4.0 degrees C increase in the Tm for 1 M sorbitol in the pH range 3.8-8.5 by scanning calorimetry. The effect increases dramatically at pH 9.5 where we observe a 9.5 degrees C stabilisation. An increase in the sorbitol concentration to 2 M stabilises lysozyme by 11.3-13.4 degrees C in the pH range 9.5-10.5. In the absence of urea, no significant effects of sorbitol were observed on the activation energy for unfolding for lysozyme at pH 4.5. This indicates together with the results from the titration experiments that sorbitol may stabilise the folded form of lysozyme by destabilising the unfolded form of lysozyme. At pH values at and above lysozyme's pI (approximately 9.3), the unfolding of the protein is accompanied with a substantial amount of self-aggregation seen in the calorimetry experiments in the ratio of DeltaH(cal)/DeltaH(vH). In the presence of sorbitol, the self-aggregation was counterbalanced by higher sorbitol concentrations. These results strongly suggest a negative influence of sorbitol on the unfolded form of lysozyme and thereby stabilising the native form.

Desolvation is a likely origin of robust enthalpic barriers to protein folding

Experimental data from global analyses of temperature (T) and denaturant dependence of the folding rates of small proteins led to an intrinsic enthalpic folding barrier hypothesis: to a good approximation, the T-dependence of folding rate under constant native stability conditions is Arrhenius. Furthermore, for a given protein, the slope of isostability folding rate versus 1/T is essentially independent of native stability. This hypothesis implies a simple relationship between chevron and Eyring plots of folding that is easily discernible when both sets of rates are expressed as functions of native stability. Using experimental data in the literature, we verify the predicted chevron-Eyring relationship for 14 proteins and determine their intrinsic enthalpic folding barriers, which vary approximately from 15 kcal/mol to 40 kcal/mol for different proteins. These enthalpic barriers do not appear to correlate with folding rates, but they exhibit correlation with equilibrium unfolding enthalpy at room temperature. Intrinsic enthalpic barriers with similarly high magnitudes apply as well to at least two cases of peptide-peptide and peptide-protein association, suggesting that these barriers are a hallmark of certain general and fundamental kinetic processes during folding and binding. Using a class of explicit-chain C(alpha) protein models with constant elementary enthalpic desolvation barriers between C(alpha) positions, we show that small microscopic pairwise desolvation barriers, which are a direct consequence of the particulate nature of water, can act cooperatively to give rise to a significant overall enthalpic barrier to folding. This theoretical finding provides a physical rationalization for the high intrinsic enthalpic barriers in protein folding energetics. Ramifications of entropy-enthalpy compensation in hydrophobic association for the height of enthalpic desolvation barrier are discussed.

Energy landscapes for cooperative processes: nearly ideal glass transitions, liquid-liquid transitions and folding transitions

We describe basic phenomenology in the physics of supercooling liquids at constant volume (most simulations), and at constant pressure (most laboratory experiments) before focusing attention on the exceptional cases that exhibit liquid-liquid phase transitions on constant-pressure cooling. We give evidence for point defects in glasses and liquids near T(g). Models based on defects predict transitions with density gaps in constant-pressure systems. We describe the energy landscape representation of such systems. Water, in these terms, is post-critical, and its nearly ideal glass formation can be related to nucleation-free protein 'funnel-folding'. For nucleated folding of proteins, a pseudo-gap should be present. Experimental methods of distinguishing between alternative folding scenarios are described.

Denaturant sensitive regions in creatine kinase identified by hydrogen/deuterium exchange

The GdmHCl-induced unfolding of creatine kinase (CK) has been studied by hydrogen/deuterium (H/D) exchange combined with mass spectrometry. MM-CK unfolded for various periods in different denaturant concentrations was pulsed-labeled with deuterium to identify different conformational intermediate states. For all denaturation times or GdmHCl concentrations, we observed variable proportions of only two species. The low-mass envelope of isotope peaks corresponds to a species that has gained about 10 deuteriums more than native CK, and the high-mass envelope to a completely deuterated species. To localize precisely the unfolded regions in the states highly populated during denaturation, the protein was digested with two proteases (pepsin and type XIII protease) after H/D exchange and rapid quenching of the reaction. The two sets of fragments obtained were analyzed by liquid chromatography coupled to mass spectrometry to determine the deuterium level in each fragment. Bimodal distributions of deuterium were found for most peptides, indicating that these regions were either folded or unfolded. This behavior is consistent with cooperative, localized unfolding. However, we observed a monomodal distribution of deuterium in two regions (1-12 and 162-186). We conclude that the increment of mass observed in the low-mass species of the intact protein (+10 Da) has its origin in these two segments. These regions, which are very sensitive to low GdmHCl concentrations, are involved in the monomer-monomer interface of CK and their perturbation is likely to weaken the dimeric structure. At higher denaturant concentration, this would induce dissociation of the dimer.

Unusually Slow Denaturation and Refolding Processes of Pyrrolidone Carboxyl Peptidase from a Hyperthermophile Are Highly Cooperative: Real-Time NMR Studies

The refolding rate of heat-denatured cysteine-free pyrrolidone carboxyl peptidase (PCP-0SH) from Pyrococcus furiosus has been reported to be unusually slow under some conditions. To elucidate the structural basis of the unusually slow kinetics of the protein, the denaturation and refolding processes of the PCP-0SH were investigated using a real-time 2D (1)H-(15)N HSQC and CD experiments. At 2 M urea denaturation of the PCP-0SH in the acidic region, all of the native peaks in the 2D HSQC spectrum completely disappeared. The conformation of the PCP-0SH just after removal of 6 M GuHCl could be observed as a stable intermediate (D(1) state) in 2D HSQC and CD experiments, which is similar to a molten globule structure. The D(1) state of the PCP-0SH, which is the initial state of refolding, corresponded to the state at 2 M urea and seemed to be the denatured state in equilibrium with the native state under the physiological conditions. The refolding of PCP-0SH from the D(1) state to the native state could be observed to be highly cooperative without any intermediates between them, even if the refolding rate was quite slow. In the higher concentration of denaturants, PCP-0SH showed HSQC and CD spectra characteristic of completely unfolded proteins called the D(2) state. The unusually slow refolding rate was discussed as originating in the conformations in the transition state and/or the retardation of reorganization in an ensemble of nonrandom denatured structures in the D(1) state.

Conformational Dynamics of the GdmHCl-Induced Molten Globule State of Creatine Kinase Monitored by Hydrogen Exchange and Mass Spectrometry

Our understanding of the mechanism of protein folding can be improved by the characterization of folding intermediate states. Intrinsic tryptophan fluorescence measurements of equilibrium GdmHCl-induced unfolding of MM-CK allow for the construction of a "phase diagram", which shows the presence of five different conformational states, including three partially folded intermediates. However, only three states are detected by using pulsed-labeled H-D exchange analyzed by electrospray ionization mass spectrometry. One of them is the native state, and the two other species are present in proportions strongly dependent on the GdmHCl concentration and denaturation time. The low-mass peak is due to a largely exchange-incompetent state, which has gained only approximately 10 deuteriums more than the native protein. This population of MM-CK molecules has undergone a small conformational change induced by low GdmHCl concentrations. However, this limited change is in itself not sufficient to inactivate the enzyme or is easily reversible. The high-mass peak corresponds to a population of MM-CK that is fully deuterated. The comparison of fluorescence, activity, and H-D exchange measurements shows that the maximally populated intermediate at 0.8 M GdmHCl has the characteristics of a molten globule. It has no activity; it has 55% of its native alpha-helices and a maximum fluorescence emission wavelength of approximately 341 nm, and it binds ANS strongly. However, no protection against exchange is detected under the conditions used in this work. This paradox, the presence of significant residual secondary and tertiary structures detected by optical probes and the total deuteration of its amide protons detected by H-D exchange and mass spectrometry, could be explained by a highly dynamic MM-CK molten globule.

Temperature-induced dissociation of protein aggregates: accessing the denatured state

The thermal denaturation of lysozyme and ribonuclease A (RNase A) under reducing and nonreducing conditions at neutral pH has been monitored by Fourier transform infrared spectroscopy. In the absence of the reductant, lysozyme and RNase A undergo apparent three- and two-state denaturation, respectively, as observed from the conformation-sensitive amide I' band. For both proteins the hydrogen-deuterium exchange takes place at lower temperatures than the main denaturation temperatures, suggesting that a transient denaturation mechanism occurs. The observed transition at 51.2 degrees C during the denaturation of lysozyme is attributed to this transient effect, rather than to the loss of tertiary structure. Under reducing conditions lysozyme aggregates during the heating phase, whereas RNase A shows only a minor aggregation, which further increases during the cooling step. The reduced stability of both proteins can be correlated with the transient denaturation behavior, which is also suggested to be involved in protein aggregation at physiologically relevant temperatures. In addition, it is shown that when the temperature is further increased, the amorphous aggregates dissociate. Comparison of the dissociated states with the denatured states obtained under nonreducing conditions indicates that these states have the same conformation. By using a two-dimensional correlation analysis we were able to show that the dissociation is preceded by a conformational change. It is argued that this extends to other types of perturbation.