Kinetic refolding of beta-lactoglobulin. Studies by synchrotron X-ray scattering, and circular dichroism, absorption and fluorescence spectroscopy

A surprisingly simple three-state generic process for reversible protein denaturation by trifluoroethanol

Despite the rich knowledge of the influence of 2,2,2-trifluoroethanol (TFE) on the structure and conformation of peptides and proteins, the mode(s) of TFE-protein interactions and the mechanism by which TFE reversibly denatures a globular protein remain elusive. This study systematically examines TFE-induced equilibrium transition curves for six paradigmatic globular proteins by using basic fluorescence and circular dichroism measurements under neutral pH conditions. The results are remarkably simple. Low TFE invariably unfolds the tertiary structure of all proteins to produce the obligate intermediate (I) which retains nearly all of native-state secondary structure, but enables the formation of extra α-helices as the level of TFE is raised higher. Inspection of the transitions at once reveals that the tertiary structure unfolding is always a distinct process, necessitating the inclusion of at least one obligate intermediate in the TFE-induced protein denaturation. It appears that the intermediate in the minimal unfolding mechanism N⇌I⇌D somehow acquires higher α-helical propensity to generate α-helices in excess of that in the native state to produce the denatured state (D), also called the TFE state. The low TFE-populated intermediate I may be called a universal intermediate by virtue of its α-helical propensity. Contrary to many earlier suggestions, this study dismisses molten globule (MG)-like attribute of I or D.

The Molten Globule, and Two-State vs. Non-Two-State Folding of Globular Proteins

From experimental studies of protein folding, it is now clear that there are two types of folding behavior, i.e., two-state folding and non-two-state folding, and understanding the relationships between these apparently different folding behaviors is essential for fully elucidating the molecular mechanisms of protein folding. This article describes how the presence of the two types of folding behavior has been confirmed experimentally, and discusses the relationships between the two-state and the non-two-state folding reactions, on the basis of available data on the correlations of the folding rate constant with various structure-based properties, which are determined primarily by the backbone topology of proteins. Finally, a two-stage hierarchical model is proposed as a general mechanism of protein folding. In this model, protein folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the molten globule state as a folding intermediate. The two-state folding is thus merely a simplified version of the hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.

A complete picture of protein unfolding and refolding in surfactants

Interactions between proteins and surfactants are of relevance in many applications including food, washing powder formulations, and drug formulation. The anionic surfactant sodium dodecyl sulfate (SDS) is known to unfold globular proteins, while the non-ionic surfactant octaethyleneglycol monododecyl ether (C₁₂E₈) can be used to refold proteins from their SDS-denatured state. While unfolding have been studied in detail at the protein level, a complete picture of the interplay between protein and surfactant in these processes is lacking. This gap in our knowledge is addressed in the current work, using the β-sheet-rich globular protein β-lactoglobulin (bLG). We combined stopped-flow time-resolved SAXS, fluorescence, and circular dichroism, respectively, to provide an unprecedented in-depth picture of the different steps involved in both protein unfolding and refolding in the presence of SDS and C₁₂E₈. During unfolding, core-shell bLG-SDS complexes were formed within ∼10 ms. This involved an initial rapid process where protein and SDS formed aggregates, followed by two slower processes, where the complexes first disaggregated into single protein structures situated asymmetrically on the SDS micelles, followed by isotropic redistribution of the protein. Refolding kinetics (>100 s) were slower than unfolding (<30 s), and involved rearrangements within the mixing deadtime (∼5 ms) and transient accumulation of unfolded monomeric protein, differing in structure from the original bLG-SDS structure. Refolding of bLG involved two steps: extraction of most of the SDS from the complexes followed by protein refolding. These results reveal that surfactant-mediated unfolding and refolding of proteins are complex processes with rearrangements occurring on time scales from sub-milliseconds to minutes.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Non-Native α-Helices in the Initial Folding Intermediate Facilitate the Ordered Assembly of the β-Barrel in β-Lactoglobulin

The roles of non-native α-helices frequently observed in the initial folding stage of β-sheet proteins have been examined for many years. We herein investigated the residue-level structures of several mutants of bovine β-lactoglobulin (βLG) in quenched-flow pH-pulse labeling experiments. βLG assumes a collapsed intermediate with a non-native α-helical structure (I₀) in the early stage of folding, although its native form is predominantly composed of β-structures. The protection profile in I₀ of pseudo-wild type (WT*) βLG was found to deviate from the pattern of the "average area buried upon folding" (AABUF). In particular, the level of protection at the region of strand A, at which non-native α-helices form in the I₀ state, was significantly low compared to AABUF. G17E, the mutant with an increased helical propensity, showed a similar protection pattern. In contrast, the protection pattern for I₀ of E44L, the mutant with an increased β-sheet propensity, was distinct from that of WT* and resembled the AABUF pattern. Transverse relaxation measurements demonstrated that the positions of the residual structures in the unfolded states of these mutants were consistent with those of the protected residues in the respective I₀ states. On the basis of the slower conversion of I₀ to the native state for E44L to that for WT*, non-native α-helices facilitate the ordered assembly of the β-barrel by preventing interactions that trap folding.

Structural characterization of MG and pre-MG states of proteins by MD simulations, NMR, and other techniques

Almost all proteins fold via a number of partially structured intermediates such as molten globule (MG) and pre-molten globule states. Understanding the structure of these intermediates at atomic level is often a challenge, as these states are observed under extreme conditions of pH, temperature, and chemical denaturants. Furthermore, several other processes such as chemical modification, site-directed mutagenesis (or point mutation), and cleavage of covalent bond of natural proteins often lead to MG like partially unfolded conformation. However, the dynamic nature of proteins in these states makes them unsuitable for most structure determination at atomic level. Intermediate states studied so far have been characterized mostly by circular dichroism, fluorescence, viscosity, dynamic light scattering measurements, dye binding, infrared techniques, molecular dynamics simulations, etc. There is a limited amount of structural data available on these intermediate states by nuclear magnetic resonance (NMR) and hence there is a need to characterize these states at the molecular level. In this review, we present characterization of equilibrium intermediates by biophysical techniques with special reference to NMR.

Determination and Quantification of Molecular Interactions in Protein Films: A Review

Protein based films are nowadays also prepared with the aim of replacing expensive, crude oil-based polymers as environmentally friendly and renewable alternatives. The protein structure determines the ability of protein chains to form intra- and intermolecular bonds, whereas the degree of cross-linking depends on the amino acid composition and molecular weight of the protein, besides the conditions used in film preparation and processing. The functionality varies significantly depending on the type of protein and affects the resulting film quality and properties. This paper reviews the methods used in examination of molecular interactions in protein films and discusses how these intermolecular interactions can be quantified. The qualitative determination methods can be distinguished by structural analysis of solutions (electrophoretic analysis, size exclusion chromatography) and analysis of solid films (spectroscopy techniques, X-ray scattering methods). To quantify molecular interactions involved, two methods were found to be the most suitable: protein film swelling and solubility. The importance of non-covalent and covalent interactions in protein films can be investigated using different solvents. The research was focused on whey protein, whereas soy protein and wheat gluten were included as further examples of proteins.

Transient non-native helix formation during the folding of β-lactoglobulin

In ideal proteins, only native interactions are stabilized step-by-step in a smooth funnel-like energy landscape. In real proteins, however, the transient formation of non-native structures is frequently observed. In this review, the transient formation of non-native structures is described using the non-native helix formation during the folding of β-lactoglobulin as a prominent example. Although β-lactoglobulin is a predominantly β-sheet protein, it has been shown to form non-native helices during the early stage of folding. The location of non-native helices, their stabilization mechanism, and their role in the folding reaction are discussed.

Transient helical structure during PI3K and Fyn SH3 domain folding

A growing list of proteins, including the β-sheet-rich SH3 domain, is known to transiently populate a compact α-helical intermediate before settling into the native structure. Examples have been discovered in cryogenic solvent as well as by pressure jumps. Earlier studies of λ repressor mutants showed that transient states with excess helix are robust in an all-α protein. Here we extend a previous study of src SH3 domain to two new SH3 sequences, phosphatidylinositol 3-kinase (PI3K) and a Fyn mutant, to see how robust such helix-rich transients are to sequence variations in this β-sheet fold. We quantify helical structure by circular dichroism (CD), protein compactness by small-angle X-ray scattering (SAXS), and transient helical populations by cryo-stopped-flow CD. Our results show that transient compact helix-rich intermediates are easily accessible on the folding landscape of different SH3 domains. In molecular dynamics simulations, force field errors are often blamed for transient non-native structure. We suggest that experimental examples of very fast α-rich transient misfolding could become a more subtle test for further force field improvements than observation of the native state alone.

Slow unfolding pathway of hyperthermophilic Tk-RNase H2 examined by pulse proteolysis using the stable protease Tk-subtilisin

The unfolding speed of some hyperthermophilic proteins is significantly slower than those of their mesostable homologues. Ribonuclease H2 from the hyperthermophilic archaeon Thermococcus kodakarensis (Tk-RNase H2) is stabilized by its remarkably slow unfolding rate. In this work, we examined the slow unfolding pathway of Tk-RNase H2 by pulse proteolysis using a superstable subtilisin-like serine protease from T. kodakarensis (Tk-subtilisin). Tk-subtilisin has enzymatic activity in highly concentrated guanidine hydrochloride (GdnHCl), in which Tk-RNase H2 unfolds slowly. The native state of Tk-RNase H2 was completely resistant to Tk-subtilisin, whereas the unfolded state (induced by 4 M GdnHCl) was degraded by Tk-subtilisin. Degradation products of Tk-RNase H2 created from pulse proteolysis during its unfolding were detected by tricine-sodium dodecyl sulfate-polyacrylamide gel electrophoresis. We identified the cleavage sites in Tk-RNase H2 by N-terminal sequencing and mass spectrometry and constructed mimics of the unfolding intermediate of Tk-RNase H2 by protein engineering. The mimics were biophysically characterized. We found that the native state of Tk-RNase H2 (N-state) changed to the I(A)-state that was digested by Tk-subtilisin in the early stage of unfolding. In the slow unfolding pathway, the I(A)-state shifted to two intermediate forms, I(B)-state and I(C)-state. The I(B)-state was digested by Tk-subtilisin in the C-terminal region, but the I(C)-state was a Tk-subtilisin resistant form. These states gradually unfolded through the I(D)-state, in which the N-terminal region was digested. The results indicate that pulse proteolysis, by a superstable protease, was a suitable strategy and an effective tool for analyzing intermediate structures of proteins with slow unfolding properties. We also showed that the N-terminal region contributes to the slow unfolding of Tk-RNase H2, and the C-terminal region is important for folding and stability.

Dynamic rheological properties of native and cross-linked gliadin proteins

A comparison of cross-linked and native gliadin suspensions, with respect to the state of protein globular structure was carried out using small-angle X-ray scattering (SAXS), dynamic light scattering (DLS) and rheological analysis. Gliadin suspensions were also analyzed in the presence and absence of glycerol. DLS analysis showed that R(h) increased only with gliadin/EDC/NHS suspensions. However, Kratky plots revealed that gliadin and gliadin/L-cysteine maintained their globular shape even in absence or presence of glycerol. Rheological experiments revealed that gliadin and gliadin/L-cysteine suspension exhibited a similar profile with three main domains, and a sol-gel transition. Gliadin/EDC/NHS did not present any sol-gel transition, and this fact corroborates with DLS results and the hypothesis of lower protein-protein interaction, which are in agreement with G″ > G'.

Kinetics of low pH-induced lamellar to bicontinuous cubic phase transition in dioleoylphosphatidylserine∕monoolein

Recently, it has been well recognized that the modulation of electrostatic interactions due to surface charges can induce transitions between lamellar liquid-crystalline (L(α)) and inverse bicontinuous double-diamond cubic (Q(II)(D)) phases in biological lipids. To reveal their kinetic pathway and mechanism, we investigated the low pH-induced L(α) to Q(II)(D) phase transitions in 20%-dioleoylphosphatidylserine (DOPS)/80%-monoolein (MO) using time-resolved small-angle x-ray scattering and a rapid mixing method. At a final pH of 2.6-2.9, the L(α) phase was transformed completely into the hexagonal II (H(II)) phase within 2-10 s after mixing a low pH buffer with a suspension of multilamellar vesicles of 20%-DOPS∕80%-MO (the initial step). Subsequently, the H(II) phase slowly converted into the Q(II)(D) phase and completely disappeared within 15-30 min (the second step). The rate constants of the second step were obtained using the singular value decomposition analysis. On the basis of these data, we discuss the underlying mechanism of the kinetic pathway of the low pH-induced L(α) to Q(II)(D) phase transitions.

A circumventing role for the non-native intermediate in the folding of β-lactoglobulin

Folding experiments have suggested that some proteins have kinetic intermediates with a non-native structure. A simple G ̅o model does not explain such non-native intermediates. Therefore, the folding energy landscape of proteins with non-native intermediates should have characteristic properties. To identify such properties, we investigated the folding of bovine β-lactoglobulin (βLG). This protein has an intermediate with a non-native α-helical structure, although its native form is predominantly composed of β-structure. In this study, we prepared mutants whose α-helical and β-sheet propensities are modified and observed their folding using a stopped-flow circular dichroism apparatus. One interesting finding was that E44L, whose β-sheet propensity was increased, showed a folding intermediate with an amount of β-structure similar to that of the wild type, though its folding took longer. Thus, the intermediate seems to be a trapped intermediate. The high α-helical propensity of the wild-type sequence likely causes the folding pathway to circumvent such time-consuming intermediates. We propose that the role of the non-native intermediate is to control the pathway at the beginning of the folding reaction.

Time resolved SAXS and RNA folding

Small angle X-ray scattering provides low resolution structural information about macromolecules in solution. When coupled with rapid mixing methods, SAXS reports time-dependent conformational changes of RNA induced by the addition of Mg(2+) to trigger folding. Thus time-resolved SAXS provides unique information about the global or overall structures of transient intermediates populated during folding. Notably, SAXS provides information about the earliest folding events, which can evade detection by other methods.

Structural analysis of an equilibrium folding intermediate in the apoflavodoxin native ensemble by small-angle X-ray scattering

Intermediate conformations are crucial to our understanding of how proteins fold into their native structures and become functional. Conventional spectroscopic measurements of thermal denaturation transitions allow the detection of equilibrium intermediates but often provide little structural detail; thus, application of more informative techniques is required. Here we used small-angle X-ray scattering (SAXS) to study the thermal denaturation of four variants of Anabaena PCC 7119 flavodoxin, including the wild-type apo and holo forms, and two mutants, E20K/E72K and F98N. Denaturation was monitored from changes in SAXS descriptors. Although the starting and final points of the denaturation were similar for the flavodoxin variants tested, substantial differences in the unfolding pathway were apparent between them. In agreement with calorimetric data, analysis of the SAXS data sets indicated a three-state unfolding equilibrium for wild-type apoflavodoxin, a two-state equilibrium for the F98N mutant, and increased thermostability of the E20K/E72K mutant and holoflavodoxin. Although the apoflavodoxin intermediate consistently appeared mixed with significant amounts of either native or unfolded conformations, its SAXS profile was derived from the deconvolution of the temperature-dependent SAXS data set. The apoflavodoxin thermal intermediate was structurally close to the native state but less compact, thereby indicating incipient unfolding. The residues that foster denaturation were explored by an ensemble of equilibrium ϕ-value restrained molecular dynamics. These simulations pointed to residues located in the cofactor and partner-protein recognition regions as the initial sites of denaturation and suggest a conformational adaptation as the mechanism of action in apoflavodoxin.

Protein stability and structure in HIC: hydrogen exchange experiments and COREX calculations

Hydrogen exchange mass spectrometry (HXMS) coupled to proteolytic digestion has been used to probe the conformation of bovine β-lactoglobulin (BLG), bovine α-lactalbumin (BLA), and human serum albumin (HSA) in solution and while adsorbed to the hydrophobic interaction chromatography media Phenyl Sepharose 6FF. All three proteins show evidence of EX1 exchange kinetics, indicating a loss of stability on the surface. HX protection patterns for all three proteins also indicate that the unfolded form is only partially solvent exposed. The hydrogen-deuterium exchange patterns of BLG and BLA on the surface suggest a structure that resembles each protein's respective solution phase molten globule state. The low stability of Domain II of HSA observed on Phenyl Sepharose 6FF also suggests a link to solution stability because Domain II is frequently cited as the least stable domain in solution unfolding pathways. COREX, an algorithm used to compute protein folding stabilities, correctly predicts solution hydrogen-deuterium exchange patterns for BLG and offers insight into its adsorbed phase stabilities but is unreliable for BLA predictions. The results of this work demonstrate a link between solution-phase local stability patterns and the nature of partially unfolded states that proteins can adopt on HIC surfaces.

Multistage collapse of a bacterial ribozyme observed by time-resolved small-angle X-ray scattering

Ribozymes must fold into compact, native structures to function properly in the cell. The first step in forming the RNA tertiary structure is the neutralization of the phosphate charge by cations, followed by collapse of the unfolded molecules into more compact structures. The specificity of the collapse transition determines the structures of the folding intermediates and the folding time to the native state. However, the forces that enable specific collapse in RNA are not understood. Using time-resolved SAXS, we report that upon addition of 5 mM Mg(2+) to the Azoarcus group I ribozyme up to 80% of chains form compact structures in less than 1 ms. In 1 mM Mg(2+), the collapse transition produces extended structures that slowly approach the folded state, while > or = 1.5 mM Mg(2+) leads to an ensemble of random coils that fold with multistage kinetics. Increased flexibility of molecules in the intermediate ensemble correlates with a Mg(2+)-dependent increase in the fast folding population and a previously unobserved crossover in the collapse kinetics. Partial denaturation of the unfolded RNA with urea also increases the fraction of chains following the fast-folding pathway. These results demonstrate that the preferred collapse mechanism depends on the extent of Mg(2+)-dependent charge neutralization and that non-native interactions within the unfolded ensemble contribute to the heterogeneity of the ribozyme folding pathways at the very earliest stages of tertiary structure formation.

alpha-Helical burst on the folding pathway of FHA domains from Rad53 and Ki67

We investigated refolding processes of beta-sheeted protein FHA domains (FHA1 domain of Rad53 and Ki67 FHA domain) by cryo-stopped-flow (SF) method combined with far-ultraviolet (far-UV) circular dichroism (CD, the average secondary structure content) and small angle X-ray scattering (SAXS, measuring the radius of gyration). In case of FHA1 domain of Rad53, no detectable time course was observed except the initial burst on its refolding process at 4 degrees C, suggesting that the FHA1 domain of Rad53 was already refolded to its native state within the dead time of the SF apparatus and the rate of the refolding is too fast to be observed at this temperature. In contrast, there was an observable alpha-helical burst at -15 degrees C and -20 degrees C in the presence of 45% ethylene glycol (EGOH) by CD-SF. Besides, the radius of gyration (Rg) of the burst phase intermediate at -20 degrees C shows the intermediate is already compact, and the compaction process was accompanied with the decrease of alpha-helical content at the same temperature. In case of Ki67 FHA domain, ellipticity change at 222 nm was observed on its refolding pathway at -28 degrees C in the presence of 45% EGOH and 2 mM DTT, indicating that Ki67 FHA domain also takes non-native alpha-helix-rich intermediate on its folding pathway. Time-resolved SAXS experiment was done. As the signal/noise ratio is low, we could not observe the time-dependent signal change through the time course. However, the initial Rg value was obtained as 18.2 +/- 0.5 A, which is much smaller than the unfolded Rg value (26.5 +/- 1.2 A), and is slightly larger than the native one (15.9 +/- 1.8 A). These results suggest that Ki67 FHA domain also forms compact non-native alpha-helix-rich intermediate before refolding to its native beta-structure on the refolding pathway. These results are in good agreement with other beta-proteins, such as bovine beta-lactoglobulin (BLG), src SH3 domain proteins. It seems the alpha-helical burst phases appear on the folding pathway of beta-sandwiched proteins.

Insights into protein aggregation by NMR characterization of insoluble SH3 mutants solubilized in salt-free water

Protein aggregation in vivo has been extensively associated with a large spectrum of human diseases. On the other hand, mechanistic insights into protein aggregation in vitro were incomplete due to the inability in solubilizing insoluble proteins for high-resolution biophysical investigations. However, a new avenue may be opened up by our recent discovery that previously-thought insoluble proteins can in fact be solubilized in salt-free water. Here we use this approach to study the NMR structural and dynamic properties of an insoluble SH3 mutant with a naturally-occurring insertion of Val22 at the tip of the diverging turn. The obtained results reveal: 1) regardless of whether the residue is Val, Ala, Asp or Arg, the insertion will render the first hNck2 SH3 domain to be insoluble in buffers. Nevertheless, all four mutants could be solubilized in salt-free water and appear to be largely unfolded as evident from their CD and NMR HSQC spectra. 2) Comparison of the chemical shift deviations reveals that while in V22-SH3 the second helical region is similarly populated as in the wild-type SH3 at pH 2.0, the first helical region is largely unformed. 3) In V22-SH3, many non-native medium-range NOEs manifest to define non-native helical conformations. In the meanwhile a small group of native-like long-range NOEs still persists, indicating the existence of a rudimentary native-like tertiary topology. 4) Although overall, V22-SH3 has significantly increased backbone motions on the ps-ns time scale, some regions still own restricted backbone motions as revealed by analyzing ¹⁵N relaxation data. Our study not only leads to the establishment of the first high-resolution structural and dynamic picture for an insoluble protein, but also shed more light on the molecular events for the nonhierarchical folding mechanism. Furthermore, a general mechanism is also proposed for in vivo protein aggregation triggered by the genetic mutation and posttranslational modification.

Slowing down downhill folding: a three-probe study

The mutant Tyr22Trp/Glu33Tyr/Gly46Ala/Gly48Ala of lambda repressor fragment lambda(6-85) was previously assigned as an incipient downhill folder. We slow down its folding in a cryogenic water-ethylene-glycol solvent (-18 to -28 degrees C). The refolding kinetics are probed by small-angle x-ray scattering, circular dichroism, and fluorescence to measure the radius of gyration, the average secondary structure content, and the native packing around the single tryptophan residue. The main resolved kinetic phase of the mutant is probe independent and faster than the main phase observed for the pseudo-wild-type. Excess helical structure formed early on by the mutant may reduce the formation of turns and prevent the formation of compact misfolded states, speeding up the overall folding process. Extrapolation of our main cryogenic folding phase and previous T-jump measurements to 37 degrees C yields nearly the same refolding rate as extrapolated by Oas and co-workers from NMR line-shape data. Taken together, all the data consistently indicate a folding speed limit of approximately 4.5 micros for this fast folder.

NMR evidence for forming highly populated helical conformations in the partially folded hNck2 SH3 domain

Recent studies of several proteins implied that the folding of beta-proteins may follow a nonhierarchical mechanism in which two major transitions are essential, i.e., the collapse of a random coil to form a nonnative helical intermediate, followed by a transformation into the native beta-structure. We report that the first hNck2 SH3 domain, assuming an all-beta barrel in the native form, can be reversibly transformed into a stable and nonnative helical state by acid-unfolding. We also conducted extensive NMR and mutagenesis studies that led to two striking findings: 1), NMR analysis reveals that in the helical state formed at pH 2.0, the first and last beta-strands in the native form become unstructured, whereas the rest is surprisingly converted into two highly populated helices with a significantly limited backbone motion; and 2), a conserved four-residue sequence is identified on the second beta-strand, a mutation of which suddenly renders the SH3 domain into a helical state even at pH 6.5, with NMR conformational and dynamic properties highly similar to those of the wild-type at pH 2.0. This observation implies that the region might contribute key interactions to disrupt the helical state, and to facilitate a further transformation into the native SH3 fold in the second transition.

Helix-rich transient and equilibrium intermediates of equine beta-lactoglobulin in alkaline buffer

Acidic buffer conditions are known to stabilize helix-rich states of even those proteins with a predominantly beta-sheet native secondary structure. Here we investigated whether such states also exist under alkaline buffer conditions. The guanidine hydrochloride (GuHCl)-induced unfolding transition and kinetic refolding of equine beta-lactoglobulin (ELG) by GuHCl-jump were investigated at pH 8.7 by far-ultraviolet circular dichroism. We found that an equilibrium intermediate appeared in 45% ethylene glycol (EGOH) buffer with 1.5 M GuHCl. The intermediate is rich in non-native alpha-helix, which is similar to the helix-rich state of ELG at pH 4.0. A kinetic study was done on the folding rate of ELG and compared with bovine beta-lactoglobulin (BLG). Transient intermediates, which were observed as the burst phase of the refolding reaction, were also rich in alpha-helix. The activation enthalpy of ELG was calculated to be c.a. 80 kJ/mol, whereas that of BLG was c.a. 70 kJ/mol in the presence of 45% EGOH. The ellipticities of the transient intermediate of ELG show temperature dependence in the presence of 45% EGOH, whereas that of BLG did not show significant dependence. This study therefore extends the existence of helix-rich equilibrium and transient intermediates of predominantly beta-sheet proteins to alkaline buffer conditions.

A stable alpha-helix-rich intermediate is formed by a single mutation of the beta-sheet protein, src SH3, at pH 3

Recently, we have found a transient intermediate on the folding pathway of src SH3. Intending to investigate the structure of the transient intermediate, we tested a mutant of src SH3, named A45G, using circular dichroism, fluorescence and X-ray solution scattering, and incidentally found that it forms a stable alpha-helix-rich intermediate (I(eq)) (different from the native beta-sheet-based secondary structure) at pH 3.0, but contains only beta-sheets at pH 6.0, whereas wild-type SH3 forms only beta-sheets at both pH 3.0 and pH 6.0. The intermediate I(eq) shows a circular dichroism measured at theta(222)=-10,300 deg.cm(2) dmol(-1), indicating a 31% alpha-helix proportion, as estimated by the CONTIN program. X-ray scattering gave the radius of gyration for I(eq) as 19.1 A at pH 3.0 and 15.4 A at pH 6.0, and Kratky plots showed a clear peak at pH 3.0, 4.0 and 6.0, indicating that I(eq) too is compact. In these parameters, I(eq) closely resembles the kinetically-obtained intermediate I(kin) which we found on the folding pathway of wild-type SH3 at pH 3.0 (radius of gyration 18.7 A and theta(222)=-8700 deg.cm(2)dmol(-1)), indicating a 26% alpha-helix proportion in our previous paper. Refolding experiments with A45G were done at pH 6.0 by stopped-flow apparatus monitored by circular dichroism, and compared to kinetic experiments with wild-type SH3 at pH 6.0. The result showed an alpha-helix-rich intermediate at the same dichroism amplitude, but nine times slower in formation-rate. A pH-jump experiment from pH 3.0 to pH 5.9 on A45G was also performed. This showed no bursts, and the rate of conformation-change was almost as fast as the refolding rate of A45G at pH 6.0. These kinetic experiment data would be consistent with I(eq) being nearly identical to the I(kin), which appeared on the folding pathways of both wild-type SH3 and A45G at pH 3.

Microsecond Hydrophobic Collapse in the Folding of Escherichia coli Dihydrofolate Reductase, an α/β-Type Protein

Using small-angle X-ray scattering combined with a continuous-flow mixing device, we monitored the microsecond compaction dynamics in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein. A significant collapse of the radius of gyration from 30 A to 23.2 A occurs within 300 micros after the initiation of refolding by a urea dilution jump. The subsequent folding after the major chain collapse occurs on a considerably longer time-scale. The protein folding trajectories constructed by comparing the development of the compactness and the secondary structure suggest that the specific hydrophobic collapse model rather than the framework model better explains the experimental observations. The folding trajectory of this alpha/beta-type protein is located between those of alpha-helical and beta-sheet proteins, suggesting that native structure determines the folding landscape.

A tightly packed hydrophobic cluster directs the formation of an off-pathway sub-millisecond folding intermediate in the alpha subunit of tryptophan synthase, a TIM barrel protein

Protein misfolding is now recognized as playing a crucial role in both normal and pathogenic folding reactions. An interesting example of misfolding at the earliest state of a natural folding reaction is provided by the alpha-subunit of tryptophan synthase, a (beta/alpha)(8) TIM barrel protein. The molecular basis for the formation of this off-pathway misfolded intermediate, I(BP), and a subsequent on-pathway intermediate, I1, was probed by mutational analysis of 20 branched aliphatic side-chains distributed throughout the sequence. The elimination of I(BP) and the substantial destabilization of I1 by replacement of a selective set of the isoleucine, leucine or valine residues (ILV) with alanine in a large ILV cluster external-to-the-barrel and spanning the N and C termini (cluster 2) implies tight-packing at most sites in both intermediates. Differential effects on I(BP) and I1 for replacements in alpha3, beta4 and alpha8 at the boundaries of cluster 2 suggest that their incorporation into I1 but not I(BP) reflects non-native folds at the edges of the crucial (beta/alpha)(1-2)beta(3) core in I(BP). The retention of I(BP) and the smaller and consistent destabilization of both I(BP) and I1 by similar replacements in an internal-to-the-barrel ILV cluster (cluster 1) and a second external-to-the-barrel ILV cluster (cluster 3) imply molten globule-like packing. The tight packing inferred, in part, for I(BP) or for all of I1 in cluster 2, but not in clusters 1 and 3, may reflect the larger size of cluster 2 and/or the enhanced number of isoleucine, leucine and valine self-contacts in and between contiguous elements of secondary structure. Tightly packed ILV-dominated hydrophobic clusters could serve as an important driving force for the earliest events in the folding and misfolding of the TIM barrel and other members of the (beta/alpha)(n) class of proteins.

Solvent-tuning the collapse and helix formation time scales of lambda(6-85)*

The lambda(6-85)(*) pseudo-wild type of lambda repressor fragment is a fast two-state folder (k(f) approximately 35 microsec(-1) at 58 degrees C). Previously, highly stable lambda(6-85)(*) mutants with k(f) > 30 microsec(-1) have been engineered to fold nearly or fully downhill. Stabilization of the native state by solvent tuning might also tune lambda(6-85)(*) away from two-state folding. We test this prediction by examining the folding thermodynamics and kinetics of lambda(6-85)(*) in a stabilizing solvent, 45% by weight aqueous ethylene glycol at -28 degrees C. Detection of kinetics by circular dichroism at 222 nm (sensitive to helix content) and small angle X-ray scattering (measuring the radius of gyration) shows that refolding from guanidine hydrochloride denatured conditions exhibits very different time scales for collapse and secondary structure formation: the two processes become decoupled. Collapse remains a low-barrier activated process, while the fastest of several secondary structure formation time scales approaches the downhill folding limit. Two-state folding of lambda(6-85)(*) is not a robust process.

Conformational changes of beta-lactoglobulin induced by anionic phospholipid

Conformational changes of beta-lactoglobulin (beta-LG) induced by anionic phospholipid (dimyristoylphosphatidylglycerol, DMPG) at physiological conditions (pH 7.0) have been investigated by UV-VIS, circular dichroism (CD) and fluorescence spectra. The experimental results suggest that beta-LG-DMPG interactions cause beta-LG a structural reorganization of the secondary structure elements accompanied by an increase in alpha-helical content, and a loosening of the protein tertiary structure. The interaction forces between beta-LG and DMPG are further evaluated by fluorescence spectra. The fluorescence spectral data show that conformational changes in the protein are driven by electrostatic interaction at first, then by hydrophobic interaction between a protein with a negative net charge and a negatively charged phospholipid.

Comparative analysis of refolding of chemically denatured β-lactoglobulin types A and B using the dilution additive mode

The kinetic refolding of beta-lactoglobulin (BLG), types A and B, by beta-cyclodextrin, glucose and sorbitol has been investigated in aqueous solution using fluorescence, far UV-CD and UV-spectrophotometric techniques. A new Pd-complex has been used to denature the protein. CD and fluorescence studies indicated that when incubated with sugar, the denatured BLG is refolded into the native-like structure through the dilution additive mode resulting in a higher yield of active protein than without sugar. CD studies show that these sugars can induce a non-native alpha-helical structure in denatured BLG-A and -B, then aid in the refolding of the protein. Based on the present study, these sugars have a different effect on BLG-A than BLG-B because of their differences in protein thermal stability. BLG-A has a higher thermal stability than BLG-B due to differences in the amino acid sequences.

Carbon nanotube transistors for biosensing applications

Electronic detection of biomolecules, although still in its early stages, is gradually emerging as an effective alternative to optical detection methods. We describe field effect transistor devices with carbon nanotube conducting channels that have been developed and used for biosensing and biodetection. Both transistors with single carbon nanotube conducting channels and devices with nanotube network conducting channels have been fabricated and their electronic characteristics examined. The devices readily respond to changes in the environment, and such effects have been examined using gas molecules and coatings with specific properties. Device operation in (conducting) buffer and in a dry environment--after buffer removal--is also discussed. Applications in the biosensing area are illustrated with three examples: the investigation of the interaction between devices and biomolecules, the electronic monitoring of biomolecular processes, and attempts to integrate cell membranes with active electronic devices.

Helical and Expanded Conformation of Equine β-Lactoglobulin in the Cold-denatured State

The thermal unfolding transition of equine beta-lactoglobulin (ELG) was investigated by circular dichroism (CD) over a temperature range of -15 degrees C to 85 degrees C. In the presence of 2 M urea, a cooperative unfolding transition was observed both with increasing and decreasing temperature. The CD spectrum indicated that the heat and cold-denatured states of ELG have substantial secondary structures but lack persistent tertiary packing of the side-chains. In order to clarify the relation between the heat or cold-denatured state and the acid-denatured (A) state characterized previously, we have attempted to observe the temperature dependence of the CD spectrum at pH 1.5. The CD spectrum in the heat-denatured state is similar to that in the A state. The CD spectrum in the A state does not change cooperatively with increasing temperature. These results indicate that the heat-denatured state and the A state are the same structural state. On the other hand, the CD intensity at acid pH cooperatively increased with decreasing temperature. The CD spectrum at low temperature and acid pH is consistent with that in the cold-denatured state. Therefore, the cold-denatured state is distinguished from the heat-denatured state or the A state, and ELG assumes a larger amount of non-native alpha-helices in the cold-denatured state. Small angle X-ray scattering and analytical ultracentrifugation have indicated that ELG assumes an expanded chain-like conformation in the cold-denatured state in contrast to the compact globular conformation in the A state. The relation between the molecular size and the helical content in the partially folded states is discussed.

Protein folding pathways from replica exchange simulations and a kinetic network model

We present an approach to the study of protein folding that uses the combined power of replica exchange simulations and a network model for the kinetics. We carry out replica exchange simulations to generate a large ( approximately 10(6)) set of states with an all-atom effective potential function and construct a kinetic model for folding, using an ansatz that allows kinetic transitions between states based on structural similarity. We use this network to perform random walks in the state space and examine the overall network structure. Results are presented for the C-terminal peptide from the B1 domain of protein G. The kinetics is two-state after small temperature perturbations. However, the coil-to-hairpin folding is dominated by pathways that visit metastable helical conformations. We propose possible mechanisms for the alpha-helix/beta-hairpin interconversion.

Temperature Weighted Histogram Analysis Method, Replica Exchange, and Transition Paths†

We analyzed the data from a replica exchange molecular dynamics simulation using the weighted histogram analysis method to combine data from all of the temperature replicas (T-WHAM) to obtain the room-temperature potential of mean force of the G-peptide (the C-terminal beta-hairpin of the B1 domain of protein G) in regions of conformational space not sampled at room temperature. We were able to determine the potential of mean force in the transition region between a minor alpha-helical population and the major beta-hairpin population and identify a possible transition path between them along which the peptide retains a significant amount of secondary structure. This observation provides new insights into a possible mechanism of formation of beta-sheet secondary structures in proteins. We developed a novel Bayesian statistical uncertainty estimation method for any quantity derived from WHAM and used it to validate the calculated potential of mean force. The feasibility of estimating regions of the potential of mean force with unfavorable free energy at room temperature by T-WHAM analysis of replica exchange simulations was further tested on a system that can be solved analytically and presented some of the same challenges found in more complex chemical systems.

Porcine beta-lactoglobulin chemical unfolding: Identification of a non-native α-helical intermediate

The chemical unfolding behavior of porcine beta-lactoglobulin (PLG) has been followed at pH 2 and 6 in the presence of guanidinium hydrochloride. The PLG unfolding transition, monitored by tryptophan fluorescence, far and near UV circular dichroism and 1D-NMR, can be described by a three-state transition suggesting the presence of at least one intermediate state that appears to display an excess of non-native alpha-helical structures. The thermodynamic parameters, as determined through a global analysis fitting procedure, give estimates of the free energy differences of the transitions connecting the native, the intermediate and the unfolded state: DeltaG(NI) (0) = 2.8 +/- 0.7 kcal mol(-1) (pH 2) and 4.2 +/- 0.5 kcal mol(-1) (pH 6) and DeltaG(NU) (0) = 7.2 +/- 0.6 kcal mol(-1) (pH 2) and 6.9 +/- 0.6 kcal mol(-1) (pH 6). CD unfolding data of the bovine species (BLG) have been collected here under the same experimental conditions of PLG to allow a careful comparison of the two beta-lactoglobulins. Intermediates with different characteristics have been identified for BLG and PLG, and their nature has been discussed on a structural analysis basis. The thermodynamic data reported here for PLG and BLG and the comparative analysis with data reported for equine beta lactoglobulin, show that homologous beta-barrel proteins, belonging to the same family and displaying high sequence identity (52-64%) populate unfolding intermediates to different extents, even though a common tendency to the formation of non-native alpha-helical intermediates, can be envisaged. The present results provide a prerequisite foundation of knowledge for the design and interpretation of future folding kinetic studies.

Localized nature of the transition-state structure in goat alpha-lactalbumin folding

To investigate whether the structure partially formed in the molten globule folding intermediate of goat alpha-lactalbumin is further organized in the transition state of folding, we constructed a number of mutant proteins and performed Phi-value analysis on them. For this purpose, we measured the equilibrium unfolding transitions and kinetic refolding and unfolding reactions of the mutants using equilibrium and stopped-flow kinetic circular dichroism techniques. The results show that the mutants with mutations located in the A-helix (V8A, L12A), the B-helix (V27A), the beta-domain (L52A, W60A), the C-helix (K93A, L96A), the C-D loop (Y103F), the D-helix (L105A, L110A), and the C-terminal 3(10)-helix (W118F), have low Phi-values, less than 0.2. On the other hand, D87N, which is located on the Ca(2+)-binding site, has a high Phi-value, 0.91, indicating that tight packing of the side-chain around Asp87 occurs in the transition state. One beta-domain mutant (I55V) and three C-helix mutants (I89V, V90A, and I95V) demonstrated intermediate Phi-values, between 0.4 and 0.7. These results indicate that the folding nucleus in the transition state of goat alpha-LA is not extensively distributed over the alpha-domain of the protein, but very localized in a region that contains the Ca(2+)-binding site and the interface between the C-helix and the beta-domain. This is apparently in contrast with the fact that the molten globule state of alpha-lactalbumin has a partially formed structure inside the alpha-domain. It is concluded that the specific docking of the alpha and beta-domains at a domain interface is necessary for this protein to organize its native structure from the molten globule intermediate.

Conformational changes of beta-lactoglobulin in sodium bis(2-ethylhexyl) sulfosuccinate reverse micelles. A fluorescence and CD study

The effect of beta-lactoglobulin encapsulation in sodium bis(2-ethylhexyl) sulfosuccinate reverse micelles on the environment of protein and on Trp was analysed at different water contents (omega0). CD data underlined the distortion of the beta-sheet and a less constrained tertiary structure as the omega0 increased, in agreement with a concomitant red shift and a decrease in the signal intensity obtained in steady-state fluorescence measurements. Fluorescence lifetimes, evaluated by biexponential analysis, were tau1 = 1.28 ns and tau2 = 3.36 ns in neutral water. In reverse micelles, decay-associated spectra indicated the occurrence of important environmental changes associated with omega0. Bimolecular fluorescence quenching by CCl4 and acrylamide was employed to analyse alterations in the accessibility of the two Trp residues in beta-lactoglobulin, induced by changes in omega0. The average bimolecular quenching constant <kq(CCl4)> was found not to depend on omega0, confirming the insolubility of this quencher in the aqueous interface, while <kq(arcylamide)> increases with omega0. The drastic decrease with omega0 of kq, associated with the longest lifetime kq2(CCl4), comparatively to the increase of kq2(acrylamide), emphasizes the location of beta-lactoglobulin in the aqueous interfacial region especially at omega0> or = 10. The fact that (omega0 = 30) >> kq2(acrylamide) (water) also confirms the important conformational changes of encapsulated beta-lactoglobulin.

Structural Stability and Solution Structure of Chaperonin GroES Heptamer Studied by Synchrotron Small-angle X-ray Scattering

The GroES protein from Escherichia coli is a well-known member of the molecular chaperones. GroES consists of seven identical 10 kDa subunits, and forms a dome-like oligomeric structure. In order to obtain information on the structural stability and unfolding-refolding mechanism of GroES protein, especially at protein concentrations (0.4-1.2 mM GroES monomer) that would mimic heat stress conditions in vivo, we have performed synchrotron small-angle X-ray scattering (SAXS) experiments. Surprisingly, in spite of the high protein concentration, reversibility in the unfolding-refolding reaction was confirmed by SAXS experiments structurally. Although the unfolding-refolding reaction showed an apparent single transition with a Cm of 1.1 M guanidium hydrochloride, a more detailed analysis of this transition demonstrated that the unfolding mechanism could be best explained by a sequential three-state model, which consists of native heptamer, dissociated monomer, and unfolded monomer. Together with our previous result that GroES unfolded completely via a partially folded monomer according to a three-state model at low protein concentration (5 microM monomer), the unfolding-refolding mechanism of GroES protein could be explained uniformly by the three-state model from low to high protein concentrations. Furthermore, to clarify an ambiguity of the native GroES structure in solution, especially mobile loop structures, we have estimated a solution structure of GroES using SAXS profiles obtained from experiments and simulation analysis. The result suggested that the native structure of GroES in solution was very similar to that seen in GroES-GroEL complex determined by crystallography.

Conserved signature proposed for folding in the lipocalin superfamily

We systematically identify a group of evolutionarily conserved residues proposed for folding in a model beta-barrel superfamily, the lipocalins. The nature of conservation at the structural level is defined and we show that the conserved residues are involved in a network of interactions that form the core of the fold. Exploratory kinetic studies are conducted with a model superfamily member, human serum retinol-binding protein, to examine their role. The present results, coupled with key experimental studies conducted with another lipocalin beta-lactoglobulin, suggest that the evolutionarily conserved regions fold on a faster folding time-scale than the non-conserved regions.

Refolding mechanism of ovalbumin: investigation by using a starting urea-denatured disulfide isomer with mispaired CYS367-CYS382

Ovalbumin, a member of the serpin superfamily, contains one cystine disulfide (Cys73-Cys120) and four cysteine sulfhydryls (Cys11, Cys30, Cys367, and Cys382) in the native state. To investigate the folding mechanism of ovalbumin, a urea-denatured disulfide isomer with a mispaired disulfide Cys367-Cys382 (D[367-382]) and its derivative (D[367-382/CM-73]) in which a native cystine counterpart of Cys73 is blocked by carboxymethylation were produced. Both the denatured isomers refolded within an instrumental dead time of 4 ms into an initial burst intermediate IN with partially folded conformation. After the initial burst phase, most of the D[367-382] molecules further refolded into the native form. In contrast, upon dilution of D[367-382/CM-73] with the refolding buffer, the protein stayed in the IN state as a stable form, which displayed a partial regain of the native secondary structure and a compact conformation with a similar Stokes radius to the native form. The structural characteristics of IN were clearly differentiated from those of an equilibrium intermediate IA that was produced by dilution with an acidic buffer of urea-denatured ovalbumin; IA showed much more hydrophobic dye binding and a larger Stokes radius than the IN state, despite their indistinguishable far-UV circular dichroic spectra. The non-productive nature of IA highlighted the importance of a compact conformation of the IN state for subsequent native refolding. These observations were consistent with a refolding model of ovalbumin that includes the regain of the partial secondary structure and of the compactness of overall conformation in an initial burst phase before the subsequent native refolding.