Molten globules are compact, partially-folded proteins postulated to be general intermediates in protein folding. Human α-lactalbumin (α-LA) is a Ca2+-binding, four-disulphide protein whose native structure is divided into two lobes, one is largely helical, the α-domain, and the other has a significant β-sheet content, the β-domain. α-LA forms a “classical” molten globule at low pH which has been studied widely as a model system of a partially-folded protein. The α-LA molten globule is compact and has a native-like helical secondary structure content. All-Ala α-LA, which has all eight native cysteines mutated to alanine, also adopts a partially-folded molten globule conformation and gives a high-quality 1H-15N HSQC spectrum at pH 2 and 40 °C. The lack of cysteine residues makes all-Ala α-LA a suitable template for spin-labelling studies. In this report we present 1H, 13C and 15N assignments for human all-Ala α-LA in its molten globule and 8 M urea-denatured states. Analysis of the chemical shift data for the molten globule state shows they are consistent with high populations of conformations in the α region of φ,ψ space for residues in the α domain of the protein. In contrast, the data for the urea-denatured state are closely similar to those expected for a random coil.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12104-026-10260-x.

Effect of urea concentration on instant refolding of Nuclear Export Protein (NEP) from Influenza-A virus H1N1: A solution NMR based investigation

Protein structural changes characterized by high-pressure, pulsed field gradient diffusion NMR spectroscopy

A method is introduced that permits direct observation of the rates at which backbone amide hydrogens become protected from solvent exchange after rapidly dropping the hydrostatic pressure inside the NMR sample cell from denaturing (2.5 kbar) to native (1 bar) conditions. The method is demonstrated for a pressure-sensitized ubiquitin variant that contains two Val to Ala mutations. Increased protection against hydrogen exchange with solvent is monitored as a function of time during the folding process. Results for 53 backbone amides show narrow clustering with protection occurring with a time constant of ca. 85 ms, but slower protection is observed around a reverse turn near the C-terminus of the protein. Remarkably, the native NMR spectrum returns with this slower time constant of ca. 150 ms, indicating that the almost fully folded protein retains molten globule characteristics with severe NMR line broadening until the final hydrogen bonds are formed. Prior to crossing the transition state barrier, hydrogen exchange protection factors are close to unity, but with slightly elevated values in the β1-β2 hairpin, previously shown to be already lowly populated in the urea-denatured state.

The factors defining the correct folding and stability of integral membrane proteins are poorly understood. Folding of only a few select membrane proteins has been scrutinised, leaving considerable deficiencies in knowledge for large protein families, such as G protein coupled receptors (GPCRs). Complete reversible folding, which is problematic for any membrane protein, has eluded this dominant receptor family. Moreover, attempts to recover receptors from denatured states are inefficient, yielding at best 40–70% functional protein. We present a method for the reversible unfolding of an archetypal family member, the β1-adrenergic receptor, and attain 100% recovery of the folded, functional state, in terms of ligand binding, compared to receptor which has not been subject to any unfolding and retains its original, folded structure. We exploit refolding on a solid support, which could avoid unwanted interactions and aggregation that occur in bulk solution. We determine the changes in structure and function upon unfolding and refolding. Additionally, we employ a method that is relatively new to membrane protein folding; pulse proteolysis. Complete refolding of β1-adrenergic receptor occurs in n-decyl-β-D-maltoside (DM) micelles from a urea-denatured state, as shown by regain of its original helical structure, ligand binding and protein fluorescence. The successful refolding strategy on a solid support offers a defined method for the controlled refolding and recovery of functional GPCRs and other membrane proteins that suffer from instability and irreversible denaturation once isolated from their native membranes.

Im7 folds via an on-pathway intermediate that contains three of the four native α-helices. The missing helix, helix III, is the shortest and its failure to be formed until late in the pathway is related to frustration in the structure. Im7H3M3, a 94-residue variant of the 87-residue Im7 in which helix III is the longest of the four native helices, also folds via an intermediate. To investigate the structural basis for this we calculated the frustration in the structure of Im7H3M3 and used NMR to investigate its dynamics. We found that the native state of Im7H3M3 is highly frustrated and in equilibrium with an intermediate state that lacks helix III, similar to Im7. Model-free analysis identified residues with chemical exchange contributions to their relaxation that aligned with the residues predicted to have highly frustrated interactions, also like Im7. Finally, we determined properties of urea-denatured Im7H3M3 and identified four clusters of interacting residues that corresponded to the α-helices of the native protein. In Im7 the cluster sizes were related to the lengths of the α-helices with cluster III being the smallest but in Im7H3M3 cluster III was also the smallest, despite this region forming the longest helix in the native state. These results suggest that the conformational properties of the urea-denatured states promote formation of a three-helix intermediate in which the residues that form helix III remain non-helical. Thus it appears that features of the native structure are formed early in folding linked to collapse of the unfolded state.

Erratum to: Backbone 1H, 13C and 15N assignments of YibK and a variant containing a unique cysteine residue at C-terminus in 8 M urea-denatured states

We have analyzed the thermodynamic properties of the von Willebrand factor (VWF) A3 domain using urea-induced unfolding at variable temperature and thermal unfolding at variable urea concentrations to generate a phase diagram that quantitatively describes the equilibrium between native and denatured states. From this analysis, we were able to determine consistent thermodynamic parameters with various spectroscopic and calorimetric methods that define the urea-temperature parameter plane from cold denaturation to heat denaturation. Urea and thermal denaturation are experimentally reversible and independent of the thermal scan rate indicating that all transitions are at equilibrium and the van't Hoff and calorimetric enthalpies obtained from analysis of individual thermal transitions are equivalent demonstrating two-state character. Global analysis of the urea-temperature phase diagram results in a significantly higher enthalpy of unfolding than obtained from analysis of individual thermal transitions and significant cross correlations describing the urea dependence of ΔH0 and ΔCP0 that define a complex temperature dependence of the m-value. Circular dichroism (CD) spectroscopy illustrates a large increase in secondary structure content of the urea-denatured state as temperature increases and a loss of secondary structure in the thermally denatured state upon addition of urea. These structural changes in the denatured ensemble make up ∼40% of the total ellipticity change indicating a highly compact thermally denatured state. The difference between the thermodynamic parameters obtained from phase diagram analysis and those obtained from analysis of individual thermal transitions illustrates that phase diagrams capture both contributions to unfolding and denatured state expansion and by comparison are able to decipher these contributions.

YibK is a tRNA methyltransferase from Haemophilus influenzae, which forms a stable homodimer in solution and contains a deep trefoil 31 knot encompassing the C-terminal helix that threads through a long loop. It has been a model system for investigating knotted protein folding pathways. Recent data have shown that the polypeptide chain of YibK remains loosely knotted under highly denaturing conditions. Here, we report (1)H, (13)C and (15)N chemical shift assignments for YibK and its variant in the presence of 8 M urea. This work forms the basis for further analysis using NMR techniques such as paramagnetic relaxation enhancement, residual dipolar couplings and spin-relaxation dynamics analysis.

YbeA is a 3-methylpseudoridine methyltransferase from Escherichia coli that forms a stable homodimer in solution. It is one of the deeply trefoil 31 knotted proteins, of which the knot encompasses the C-terminal helix that threads through a long loop. Recent studies on the knotted protein folding pathways using YbeA have suggested that the protein knot remains present under chemically denaturing conditions. Here, we report (1)H, (13)C and (15)N chemical shift assignments for urea-denatured YbeA, which will serve as the basis for further structural characterisations using solution state NMR spectroscopy with paramagnetic spin labeled and partial alignment media.

We present here the characterization of the structural, dynamics, and energetics of properties of the urea-denatured state of ubiquitin, a small prototypical soluble protein. By combining state-of-the-art molecular dynamics simulations with NMR and small-angle X-ray scattering data, we were able to: (i) define the unfolded state ensemble, (ii) understand the energetics stabilizing unfolded structures in urea, (iii) describe the dedifferential nature of the interactions of the fully unfolded proteins with urea and water, and (iv) characterize the early stages of protein refolding when chemically denatured proteins are transferred to native conditions. The results presented herein are unique in providing a complete picture of the chemically unfolded state of proteins and contribute to deciphering the mechanisms that stabilize the native state of proteins, as well as those that maintain them unfolded in the presence of urea.

Thermomyces lanuginosus lipase (TlL) is a kinetically stable protein, resistant toward both denaturation and refolding in the presence of the ionic surfactant sodium dodecyl sulfate (SDS) and the nonionic surfactant decyl maltoside (DecM). We investigate the pH dependence of this kinetic stability. At pH 8, TlL remains folded and enzymatically active at multimillimolar surfactant concentrations but fails to refold from the acid urea-denatured state at submillimolar concentrations of SDS and DecM, indicating a broad concentration range of kinetic trapping or hysteresis. At pH 8, very few SDS molecules bind to TlL. The hysteresis SDS concentration range shrinks when moving to pH 4-6; in this pH range, SDS binds as micellelike clusters. Although hysteresis can be eliminated by reducing disulfide bonds, destabilizing the native state, and lowering the unfolding activation barrier, SDS sensitivity is not directly linked to intrinsic kinetic stability [its resistance to the general chemical denaturant guanidinium chloride (GdmCl)], because TlL unfolds more slowly in GdmCl at pH 6.0 than at pH 8.0. However, the estimated net charge drops from approximately -12 to approximately -5 between pH 8 and 6. SDS denatures TlL at pH 6.0 by nucleating via a critical number of bound SDS molecules on the surface of native TlL to form clusters. These results imply that SDS sensitivity is connected to the availability of appropriately charged regions on the protein. We suggest that conformational rigidity is a necessary but not sufficient feature of SDS resistance, because this has to be combined with sufficient negative electrostatic potential to avoid extensive SDS binding.

The outer membrane protein OmpA from Escherichia coli can fold into lipid vesicles and surfactant micelles from the urea-denatured state. However, a complete kinetic description of the folding and unfolding of OmpA, which can provide the basis for subsequent protein engineering studies of the protein's folding pathway, is lacking. Here we use two different denaturants to probe the unfolding mechanism of OmpA in the presence of the surfactant octyl maltoside (OM). Unfolding of OmpA in the presence of micelles, achieved with the potent denaturant guanidinium chloride (GdmCl), leads to single-phase unfolding. In contrast, OmpA unfolds in urea only below OM's critical micelle concentration, and this occurs in different phases, which we attribute to the existence of states that have bound different amounts of surfactant, from completely "naked" to partly covered by surfactant. Multiple parallel refolding phases are attributed to different levels of collapse prior to folding. Kinetic results used to derive the stability of OmpA in surfactant, using either urea or GdmCl as the denaturing agent, give comparable results and indicate a minimalist three-state folding scheme involving denatured state D, folding intermediate I, and native state N. N and I are stabilized by 15.6 and 2.6 kcal/mol, respectively, relative to D. The periplasmic domain of OmpA does not contribute to stability in surfactant micelles. However, BBP, a minimalist transmembrane β-barrel version of OmpA with shortened loops, is destabilized by ~10 kcal/mol compared to OmpA, highlighting loop contributions to OmpA stability.

Atomic-level analyses of non-native protein ensembles constitute an important aspect of protein folding studies to reach a more complete understanding of how proteins attain their native form exhibiting biological activity. Previously, formation of hydrophobic clusters in the 6 M urea-denatured state of an ultrafast folding mini-protein known as TC5b from both photo-CIDNP NOE transfer studies and FCS measurements was observed. Here, we elucidate the structural properties of this mini-protein denatured in 6 M urea performing 15N NMR relaxation studies together with a thorough NOE analysis. Even though our results demonstrate that no elements of secondary structure persist in the denatured state, the heterogeneous distribution of R 2 rate constants together with observing pronounced heteronuclear NOEs along the peptide backbone reveals specific regions of urea-denatured TC5b exhibiting a high degree of structural rigidity more frequently observed for native proteins. The data are complemented with studies on two TC5b point mutants to verify the importance of hydrophobic interactions for fast folding. Our results corroborate earlier findings of a hydrophobic cluster present in urea-denatured TC5b comprising both native and non-native contacts underscoring their importance for ultra rapid folding. The data assist in finding ways of interpreting the effects of pre-existing native and/or non-native interactions on the ultrafast folding of proteins; a fact, which might have to be considered when defining the starting conditions for molecular dynamics simulation studies of protein folding.

The unfolded ensemble in aqueous solution represents the starting point of protein folding. Characterisation of this species is often difficult since the native state is usually predominantly populated at equilibrium. Previous work has shown that the four-helix protein, Im7 (immunity protein 7), folds via an on-pathway intermediate. While the transition states and folding intermediate have been characterised in atomistic detail, knowledge of the unfolded ensemble under the same ambient conditions remained sparse. Here, we introduce destabilising amino acid substitutions into the sequence of Im7, such that the unfolded state becomes predominantly populated at equilibrium in the absence of denaturant. Using far- and near-UV CD, fluorescence, urea titration and heteronuclear NMR experiments, we show that three amino acid substitutions (L18A–L19A–L37A) are sufficient to prevent Im7 folding, such that the unfolded state is predominantly populated at equilibrium. Using measurement of chemical shifts, 15N transverse relaxation rates and sedimentation coefficients, we show that the unfolded species of L18A–L19A–L37A deviates significantly from random-coil behaviour. Specifically, we demonstrate that this unfolded species is compact (Rh=25 Å) relative to the urea-denatured state (Rh≥30 Å) and contains local clusters of hydrophobic residues in regions that correspond to the four helices in the native state. Despite these interactions, there is no evidence for long-range stabilising tertiary interactions or persistent helical structure. The results reveal an unfolded ensemble that is conformationally restricted in regions of the polypeptide chain that ultimately form helices I, II and IV in the native state.

Investigating the Refolding Pathway of Human Acidic Fibroblast Growth Factor (hFGF-1) from the Residual Structure(s) Obtained by Denatured-State Hydrogen/Deuterium Exchange

The interior of cells is highly crowded with macromolecules, which impacts all physiological processes. To explore how macromolecular crowding may influence cellular protein folding, we interrogated the folding landscape of a model beta-rich protein, cellular retinoic acid-binding protein I (CRABP I), in the presence of an inert crowding agent (Ficoll 70). Urea titrations revealed a crowding-induced change in the water-accessible polar amide surface of its denatured state, based on an observed ca. 15% decrease in the change in unfolding free energy with respect to urea concentration (the m-value), and the effect of crowding on the equilibrium stability of CRABP I was less than our experimental error (i.e., < or = 1.2 kcal/mol). Consequently, we directly probed the effect of crowding on the denatured state of CRABP I by measuring side-chain accessibility using iodide quenching of tryptophan fluorescence and chemical modification of cysteines. We observed that the urea-denatured state is more compact under crowded conditions, and the observed extent of reduction of the m-value by crowding agent is fully consistent with the extent of reduction of the accessibility of the Trp and Cys probes, suggesting a random and nonspecific compaction of the unfolded state. The thermodynamic consequences of crowding-induced compaction are discussed. In addition, over a wide range of Ficoll concentration, crowding significantly retarded the unfolding kinetics of CRABP I without influencing the urea dependence of the unfolding rate, arguing for no appreciable change in the nature of the transition state. Our results demonstrate how macromolecular crowding may influence protein folding by effects on both the unfolded state ensemble and unfolding kinetics.

The sequence specific (1)H, (13)C and (15)N resonance assignments of hahellin in 8M urea-denatured state have been accomplished by NMR spectroscopy. Secondary chemical shift analysis reveals the native-like propensities for β-rich conformation in the denatured state.

The Escherichia coli outer membrane protein X (OmpX) contains two polypeptide segments that present nonrandom residual structure in 8 M aqueous urea, whereas the remainder of the protein is in a flexibly disordered conformation (Tafer et al. in Biochemistry 43:860-869, 2004). In the present study, the results of two long-timescale (0.4 micros) unrestrained explicit-solvent molecular dynamics (MD) simulations of a tetradecapeptide representative of one of these two segments in 8 M aqueous urea are reported and analyzed. The two simulations were initiated either from the conformation of the corresponding segment in an NMR model structure of the unfolded protein or from an entirely extended configuration. The sampled conformational ensembles agree qualitatively with the experimentally observed NOEs, but not quantitatively, suggesting that a number of relevant configurations were not visited on the 2 x 0.4 micros timescale. Major conformational transitions occur on the 0.1 micros timescale, and the ensembles corresponding to the two independent simulations overlap only to a limited extent. However, both simulations show in multiple events the reversible formation and disruption of alpha-helical secondary structure (characteristic of the urea-denatured state) and beta-turn secondary structure (characteristic of the native state). Events of helix formation are correlated with the appearance of hydrogen bonds between two side chains (Asp75-Ser78) and of a persistent hydrophobic contact (Trp76-Tyr80). They also evidence a peculiar helix stabilization and N-terminal capping role for a negatively charged residue (Asp75). These features are in good qualitative agreement with the NMR model for the structured state of the corresponding segment in the urea-denatured protein. The analysis of the simulations provides a detailed picture of the structural and dynamic features of the considered peptide at atomic resolution that is of high relevance in the understanding of the OmpX folding process.

Current NMR information on side-chain conformations of unfolded protein states is sparse due to the poor dispersion particularly of side-chain proton resonances. We present here optimized schemes for the detection of (3)J(HalphaHbeta), (3)J(NHbeta), and (3)J(C'Hbeta) scalar and (1)D(CbetaHbeta) residual dipolar couplings (RDCs) in unfolded proteins. For urea-denatured ubiquitin and protein G, up to six (3)J-couplings to (1)H(beta) are detected, which define the chi(1) angle at very high precision. Interpretation of the (3)J couplings by a model of mixed staggered chi(1) rotamers yields excellent agreement and also provides stereoassignments for (1)H(beta) methylene protons. For all observed amino acids with the exception of leucine, the chemical shift of (1)H(beta3) protons was found downfield from (1)H(beta2). For most residues, the precision of individual chi(1) rotamer populations is better than 2%. The experimental chi(1) rotamer populations are in the vicinity of averages obtained from coil regions in folded protein structures. However, individual variations from these averages of up to 40% are highly significant and indicate sequence- and residue-specific interactions. Particularly strong deviations from the coil average are found for serine and threonine residues, an effect that may be explained by a weakening of side-chain to backbone hydrogen bonds in the urea-denatured state. The measured (1)D(CbetaHbeta) RDCs correlate well with predicted RDCs that were calculated from a sterically aligned coil model ensemble and the (3)J-derived chi(1) rotamer populations. This agreement supports the coil model as a good first approximation of the unfolded state. Deviations between measured and predicted values at certain sequence locations indicate that the description of the local backbone conformations can be improved by incorporation of the RDC information. The ease of detection of a large number of highly precise side-chain RDCs opens the possibility for a more rigorous characterization of both side-chain and backbone conformations in unfolded proteins.