Exploration of partially unfolded states of human alpha-lactalbumin by molecular dynamics simulation

Molecular dynamics simulations are used to probe the properties of non-native states of the protein human alpha-lactalbumin (human alpha-LA) with a detailed atomistic model in an implicit aqueous solvent environment. To sample the conformational space, a biasing force is introduced that increases the radius of gyration relative to the native state and generates a large number of low-energy conformers that differ in terms of their root-mean-square deviation, for a given radius of gyration. The resulting structures are relaxed by unbiased simulations and used as models of the molten globule and partly denatured states of human alpha-LA, based on measured radii of gyration obtained from nuclear magnetic resonance experiments. The ensembles of structures agree in their overall properties with experimental data available for the human alpha-LA molten globule and its more denatured states. In particular, the simulation results show that the native-like fold of the alpha-domain is preserved in the molten globule. Further, a considerable proportion of the antiparallel beta-strand in the beta-domain are present. This indicates that the lack of hydrogen exchange protection found experimentally for the beta-domain is due to rearrangement of the beta-sheet involving transient populations of non-native beta-structures. The simulations also provide details concerning the ensemble of structures that contribute as the molten globule unfolds and shows, in accord with experimental data, that unfolding is not cooperative; i.e. the various structural elements do not unfold simultaneously.     

Energetic basis of structural stability in the molten globule state: alpha-lactalbumin

The denatured states of alpha-lactalbumin, which have features of a molten globule state, have been studied to elucidate the energetics of the molten globule state and its contribution to the stability of the native conformation. Analysis of calorimetric and CD data shows that the heat capacity increment of alpha-lactalbumin denaturation highly correlates with the degree of disorder of the residual structure of the state. As a result, the denaturational transition of alpha-lactalbumin from the native to a highly ordered compact denatured state, and from the native to the disordered unfolded state are described by different thermodynamic functions. The enthalpy and entropy of the denaturation of alpha-lactalbumin to compact denatured state are always greater than the enthalpy and entropy of its unfolding. This difference represents the unfolding of the molten globule state. Calorimetric measurements of the heat effect associated with the unfolding of the molten globule state reveal that it is negative in sign over the temperature range of molten globule stability. This observation demonstrates the energetic specificity of the molten globule state, which, in contrast to a protein with unique tertiary structure, is stabilized by the dominance of negative entropy and enthalpy of hydration over the positive conformational entropy and enthalpy of internal interactions. It is concluded that at physiological temperatures the entropy of dehydration is the dominant factor providing stability for the compact intermediate state on the folding pathway, while for the stability of the native state, the conformational enthalpy is the dominant factor.     

Characterization of the molten globule state of retinol-binding protein using a molecular dynamics simulation approach

Retinol-binding protein transports retinol, and circulates in the plasma as a macromolecular complex with the protein transthyretin. Under acidic conditions retinol-binding protein undergoes a transition to the molten globule state, and releases the bound retinol ligand. A biased molecular dynamics simulation method has been used to generate models for the ensemble of conformers populated within this molten globule state. Simulation conformers, with a radius of gyration at least 1.1 A greater than that of the native state, contain on average 37%beta-sheet secondary structure. In these conformers the central regions of the two orthogonal beta-sheets that make up the beta-barrel in the native protein are highly persistent. However, there are sizable fluctuations for residues in the outer regions of the beta-sheets, and large variations in side chain packing even in the protein core. Significant conformational changes are seen in the simulation conformers for residues 85-104 (beta-strands E and F and the E-F loop). These changes give an opening of the retinol-binding site. Comparisons with experimental data suggest that the unfolding in this region may provide a mechanism by which the complex of retinol-binding protein and transthyretin dissociates, and retinol is released at the cell surface.     

Structural characterization of partially folded intermediates of apomyoglobin H64F

We present a detailed investigation of unfolded and partially folded states of a mutant apomyoglobin (apoMb) where the distal histidine has been replaced by phenylalanine (H64F). Previous studies have shown that substitution of His64, located in the E helix of the native protein, stabilizes the equilibrium molten globule and native states and leads to an increase in folding rate and a change in the folding pathway. Analysis of changes in chemical shift and in backbone flexibility, detected via [1H]-15N heteronuclear nuclear Overhauser effect measurements, indicates that the phenylalanine substitution has only minor effects on the conformational ensemble in the acid- and urea-unfolded states, but has a substantial effect on the structure, dynamics, and stability of the equilibrium molten globule intermediate formed near pH 4. In H64F apomyoglobin, additional regions of the polypeptide chain are recruited into the compact core of the molten globule. Since the phenylalanine substitution has negligible effect on the unfolded ensemble, its influence on folding rate and stability comes entirely from interactions within the compact folded or partly folded states. Replacement of His64 with Phe leads to favorable hydrophobic packing between the helix E region and the molten globule core and leads to stabilization of helix E secondary structure and overall thermodynamic stabilization of the molten globule. The secondary structure of the equilibrium molten globule parallels that of the burst phase kinetic intermediate; both intermediates contain significant helical structure in regions of the polypeptide that comprise the A, B, E, G, and H helices of the fully folded protein.     

Molecular dynamics model structures for the molten globule state of alpha-lactalbumin: aromatic residue clusters I and II

To model the molten globule structure of alpha-lactalbumin, molecular dynamics (MD) simulations were carried out for the protein in explicit water at high temperature. In these simulations, long-range Coulomb interactions were evaluated explicitly with an original method (particle-particle and particle-cell: PPPC) to avoid artifacts caused by the cut-off. The MD simulations were started from two initial conditions to verify that similar results would be obtained. From the last 150 ps trajectories of the two MD simulations, two partially unfolded average structures were obtained. These structures had the following common structural features which are characteristic of the molten globule state. The radii of gyration for these conformations were 7.4 and 9.6% larger than that of the native state. These values were almost the same as the experimental value (9.6%) observed recently by small-angle X-ray scattering (Kataoka,M., Kuwajima,K., Tokunaga,F. and Goto,Y., 1997, Protein Sci., 6, 422-430). Furthermore, aromatic residues of clusters I and II in these structures were far apart from each other except for Try103-Trp104. This result is in good agreement with NMR experimental results for the acid-denatured molten globule state (Alexandrescu et al., 1992, 1993); that is, NOE signals between the aromatic residues were not observed, except for that of Try103-Trp104 in the molten globule state. Other structural features of these models for the molten globule state are discussed with reference to native state structures.     

Illuminating the Off-Pathway Nature of the Molten Globule Folding Intermediate of an α-β Parallel Protein

Partially folded protein species transiently form during folding of most proteins. Often, these species are molten globules, which may be on- or off-pathway to the native state. Molten globules are ensembles of interconverting protein conformers that have a substantial amount of secondary structure, but lack virtually all tertiary side-chain packing characteristics of natively folded proteins. Due to solvent-exposed hydrophobic groups, molten globules are prone to aggregation, which can have detrimental effects on organisms. The molten globule observed during folding of the 179-residue apoflavodoxin from Azotobacter vinelandii is off-pathway, as it has to unfold before native protein can form. Here, we study folding of apoflavodoxin and characterize its molten globule using fluorescence spectroscopy and Förster Resonance Energy Transfer (FRET). Apoflavodoxin is site-specifically labeled with fluorescent donor and acceptor dyes, utilizing dye-inaccessibility of Cys69 in cofactor-bound protein. Donor (i.e., Alexa Fluor 488) is covalently attached to Cys69 in all apoflavodoxin variants used. Acceptor (i.e., Alexa Fluor 568) is coupled to Cys1, Cys131 and Cys178, respectively. Our FRET data show that apoflavodoxin’s molten globule forms in a non-cooperative manner and that its N-terminal 69 residues fold last. In addition, striking conformational differences between molten globule and native protein are revealed, because the inter-label distances sampled in the 111-residue C-terminal segment of the molten globule are shorter than observed for native apoflavodoxin. Thus, FRET sheds light on the off-pathway nature of the molten globule during folding of an α-β parallel protein.     

Structural characterization of the molten globule of alpha-lactalbumin by solution X-ray scattering.

A compact denatured state is often observed under a mild denaturation condition for various proteins. A typical example is the alpha-lactalbumin molten globule. Although the molecular compactness and shape are the essential properties for defining the molten globule, there have been ambiguities of these properties for the molten globule of alpha-lactalbumin. Using solution X-ray scattering, we have examined the structural properties of two types of molten globule of alpha-lactalbumin, the apo-protein at neutral pH and the acid molten globule. The radius of gyration for the native holo-protein was 15.7 A, but the two different molten globules both had a radius of gyration of 17.2 A. The maximum dimension of the molecule was also increased from 50 A for the native state to 60 A for the molten globule. These values clearly indicate that the molten globule is not as compact as the native state. The increment in the radius of gyration was less than 10% for the alpha-lactalbumin molten globule, compared with up to 30% for the molten globules of other globular proteins. Intramolecular disulfide bonds restrict the molecular expansion of the molten globule. The distance distribution function of the alpha-lactalbumin molten globule is composed of a single peak suggesting a globular shape, which is simply swollen from the native state. The scattering profile in the high Q region of the molten globule indicates the presence of a significant amount of tertiary fold. Based on the structural properties obtained by solution X-ray scattering, general and conceptual structural images for the molten globules of various proteins are described and compared with the individual, detailed structural model obtained by nuclear magnetic resonance.     

Comparison of the molten globule states of thermophilic and mesophilic alpha-amylases

In recent years great interest has been generated in the process of protein folding, and the formation of intermediates during the folding process has been proven with new experimental strategies. In the present work, we have examined the molten globule state of Bacillus licheniformis alpha-amylase (BLA) by intrinsic fluorescence and circular dichroism spectra, 1-anilino naphthalene-8-sulfonate (ANS) binding and proteolytic digestion by pepsin, for comparison to its mesophilic counterpart, Bacillus amyloliquefaciens alpha-amylase (BAA). At pH 4.0, both enzymes acquire partially folded state which show characteristics of molten globule state. They unfold in such a way that their hydrophobic surfaces are exposed to a greater extent compared to the native forms. Chemical denaturation studies by guanidine hydrochloride and proteolytic digestion with pepsin show that molten globule state of BLA is more stable than from BAA. Results from gel filtration indicate that BAA has the same compactness at pH 4.0 and 7.5. However, molten globule state of BLA is less compact than its native state. The effects of polyols such as trehalose, sorbitol and glycerol on refolding of enzymes from molten globule to native state were also studied. These polyols are effective on refolding of mesophilic alpha-amylase but only slightly effect on BLA refolding. In addition, the folding pathway and stability of intermediate state of the thermophilic and the mesophilic alpha-amylases are discussed.     

High pressure NMR reveals that apomyoglobin is an equilibrium mixture from the native to the unfolded

Pressure-induced reversible conformational changes of sperm whale apomyoglobin have been studied between 30 bar and 3000 bar on individual residue basis by utilizing 1H/15N hetero nuclear single-quantum coherence two-dimensional NMR spectroscopy at pH 6.0 and 35 degrees C. Apomyoglobin showed a series of pressure-dependent NMR spectra as a function of pressure, assignable to the native (N), intermediates (I), molten globule (MG) and unfolded (U) conformers. At 30 bar, the native fold (N) shows disorder only in the F helix. Between 500 bar and 1200 bar, a series of locally disordered conformers I are produced, in which local disorder occurs in the C helix, the CD loop, the G helix and part of the H helix. At 2000 bar, most cross-peaks exhibit severe line-broadening, suggesting the formation of a molten globule, but at 3000 bar all the cross-peaks reappear, showing that the molten globule turns into a well-hydrated, mobile unfolded conformation U. Since all the spectral changes were reversible with pressure, apomyoglobin is considered to exist as an equilibrium mixture of the N, I, MG and U conformers at all pressures. MG is situated at 2.4+/-(0.1) kcal/mol above N at 1 bar and the unfolding transition from the combined N-I state to MG is accompanied by a loss of partial molar volume by 75+/-(3) ml/mol. On the basis of these observations, we postulate a theorem that the partial molar volume of a protein decreases in parallel with the loss of its conformational order.     

Hydrogen exchange in a large 29 kD protein and characterization of molten globule aggregation by NMR

The nature of denatured ensembles of the enzyme human carbonic anhydrase (HCA) has been extensively studied by various methods in the past. The protein constitutes an interesting model for folding studies that does not unfold by a simple two-state transition, instead a molten globule intermediate is highly populated at 1.5 M GuHCl. In this work, NMR and H/D exchange studies have been conducted on one of the isozymes, HCA I. The H/D exchange studies, which were enabled by the previously obtained resonance assignment of HCA I, have been used to identify unfolded forms that are accessible from the native state. In addition, the GuHCl-induced unfolded states of HCA I have also been characterized by NMR at GuHCl concentrations in the 0-5 M range. The most important findings in this work are as follows: (1) Amide protons located in the center of the beta-sheet require global unfolding events for efficient H/D exchange. (2) The molten globule and the native state give similar protection against H/D exchange for all of the observable amide protons (i.e., water seems not to efficiently penetrate the interior of the molten globule). (3) At high protein concentrations, the molten globule can form large aggregates, which are not detectable by solution-state NMR methods. (4) The unfolded state (U), present at GuHCl concentrations above 2 M, is composed of an ensemble of conformations having residual structures with different stabilities.     

Side chain accessibility and dynamics in the molten globule state of alpha-lactalbumin: a (19)F-NMR study

The molten globule state of alpha-lactalbumin (alpha-LA) has been considered a prototype of partially folded proteins. Despite the importance of molten globules in understanding the mechanisms of protein folding and its relevance to some biological phenomena, site-specific information on the structure and dynamics of a molten globule is limited, largely because of the high conformational flexibility and heterogeneity. Here, we use selective isotope labeling and (19)F NMR to investigate the solvent accessibility and side-chain dynamics of aromatic residues in the molten globule of alpha-LA. Comparison of these properties with those of the native and unfolded protein indicates that the alpha-LA molten globule is highly heterogeneous; each residue has its unique solvent accessibility and motional environment. Many aromatic residues normally buried in the interior of native alpha-LA remain significantly buried in the molten globule and the side-chain dynamics of these residues are highly restricted. Our results suggest that hydrophobic and van der Waals interactions mediated by the inaccessible surface area could be sufficient to account for all the stability of the alpha-LA molten globule, which is approximately 50% of the value for the native protein.     

Local and global cooperativity in the human alpha-lactalbumin molten globule

NMR spectroscopy has been used to follow the urea-induced unfolding of the low pH molten globule states of a single-disulfide variant of human alpha-lactalbumin ([28-111] alpha-LA) and of two mutants, each with a single proline substitution in a helix. [28-111] alpha-LA forms a molten globule very similar to that formed by the wild-type four-disulfide protein, and this variant has been used as a model for the alpha-lactalbumin (alpha-LA) molten globule in a number of studies. The urea-induced unfolding behavior of [28-111] alpha-LA is similar to that of the four-disulfide form of the protein, except that [28-111] alpha-LA is less stable and has greater cooperativity in the loss of different elements of structure. For one mutant, L11P, the helix containing the mutation is highly destabilized such that it is completely unfolded even in the absence of urea. By contrast, for the other mutant, Q117P, the helix containing the mutation retains its compact structure. Both mutations, however, show significant long-range destabilization of the overall fold showing that the molten globule state has a degree of global cooperativity. The results reveal that different permutations of three of the four major alpha-helices of the protein can form a stable, locally cooperative, compact structural core. Taken together, these findings demonstrate that the molten globule state of alpha-LA is an ensemble of conformations, with different subsets of structures linked by a range of long-range interactions.     

Molecular basis of cooperativity in protein folding. V. Thermodynamic and structural conditions for the stabilization of compact denatured states

The heat-denatured state of proteins has been usually assumed to be a fully hydrated random coil. It is now evident that under certain solvent conditions or after chemical or genetic modifications, the protein molecule may exhibit a hydrophobic core and residual secondary structure after thermal denaturation. This state of the protein has been called the “compact denatured” or “molten globule” state. Recently is has been shown that α-lactalbumin at pH < 5 denatures into a molten globule state upon increasing the temperature (Griko, Y., Freire, E., Privalov, P. L. Biochemistry 33:1889–1899, 1994). This state has a lower heat capacity and a higher enthalpy at low temperatures than the unfolded state. At those temperatures the stabilization of the molten globule state is of an entropic origin since the enthalpy contributes unfavorably to the Gibbs free energy. Since the molten globule is more structured than the unfolded state and, therefore, is expected to have a lower configurational entropy, the net entropic gain must originate primarily from solvent related entropy arising from the hydrophobic effect, and to a lesser extent from protonation or electrostatic effects. In this work, we have examined a large ensemble of partly folded states derived from the native structure of α-lactalbumin in order to identify those states that satisfy the energetic criteria of the molten globule. It was found that only few states satisfied the experimental constraints and that, furthermore, those states were part of the same structural family. In particular, the regions corresponding to the A, B, and C helices were found to be folded, while the β sheet and the D helix were found to be unfolded. At temperatures below 45°C the states exhibiting those structural characteristics are enthalpically higher than the unfolded state in agreement with the experimental data. Interestingly, those states have a heat capacity close to that observed for the acid pH compact denatured state of α-lactalbumin [980 cal (mol.K)^?l]. In addition, the folded regions of these states include those residues found to be highly protected by NMR hydrogen exchange experiments. This work represents an initial attempt to model the structural origin of the thermodynamic properties of partly folded states. The results suggest a number of structural features that are consistent with experimental data. © 1994 Wiley-Liss, Inc.     

Validation of NMR side-chain conformations by packing calculations

The precision and accuracy of protein structures determined by nuclear magnetic resonance (NMR) spectroscopy depend on the completeness of input experimental data set. Typically, rather than a single structure, an ensemble of up to 20 equally representative conformers is generated and routinely deposited in the Protein Database. There are substantially more experimentally derived restraints available to define the main-chain coordinates than those of the side chains. Consequently, the side-chain conformations among the conformers are more variable and less well defined than those of the backbone. Even when a side chain is determined with high precision and is found to adopt very similar orientations among all the conformers in the ensemble, it is possible that its orientation might still be incorrect. Thus, it would be helpful if there were a method to assess independently the side-chain orientations determined by NMR. Recently, homology modeling by side-chain packing algorithms has been shown to be successful in predicting the side-chain conformations of the buried residues for a protein when the main-chain coordinates and sequence information are given. Since the main-chain coordinates determined by NMR are consistently more reliable than those of the side-chains, we have applied the side-chain packing algorithms to predict side-chain conformations that are compatible with the NMR-derived backbone. Using four test cases where the NMR solution structures and the X-ray crystal structure of the same protein are available, we demonstrate that the side-chain packing method can provide independent validation for the side-chain conformations of NMR structures. Comparison of the side-chain conformations derived by side-chain packing prediction and by NMR spectroscopy demonstrates that when there is agreement between the NMR model and the predicted model, on average 78% of the time the X-ray structure also concurs. While the side-chain packing method can confirm the reliable residue conformations in NMR models, more importantly, it can also identify the questionable residue conformations with an accuracy of 60%. This validation method can serve to increase the confidence level for potential users of structural models determined by NMR.     

Increased exposure of hydrophobic surface in molten globule state of alpha-lactalbumin. Fluorescence and hydrophobic photolabeling studies.

The involvement of molten globule state as a distinct intermediate in the denaturation process in proteins is well documented. However, the structural characterization of such an intermediate is far from complete. We have, using fluorescence and fluorescence quenching, studied the molten globule state of bovine alpha-lactalbumin. Unlike the native state, where all the 4 tryptophans are buried in the protein, 2 tryptophans are exposed in the molten globule state. Using the hydrophobic photoactivable reagent [3H]diazofluorene, we observe an increased hydrophobic exposure in the molten globule state. These structural characteristics conform to the current views on the molten globule state, i.e. it has similar secondary structure but a poorly defined tertiary structure. Our fluorescence studies indicate the involvement of a premolten globule state in the native to molten globule state transition. This premolten globule state exists at pH 5.0 and has a very compact structure involving increased hydrophobic interactions in the protein interior. These results are also supported by circular dichroism studies.     

Mechanism of formation of a productive molten globule form of barstar

Structural analysis of the initial steps in protein folding is difficult because of the swiftness with which these steps occur. Hence, the link between initial polypeptide chain collapse and formation of secondary and other specific structures remains poorly understood. Here, an equilibrium model has been developed for characterizing the initial steps of folding of the small protein barstar, which lead to the formation of a productive molten globule in the folding pathway. In this model, the high-pH-unfolded form (D form) of barstar, which is shown to be as unstructured as the urea-denatured form, is transformed progressively into a molten globule B form by incremental addition of the salt Na(2)SO(4) at pH 12. At very low concentrations of Na(2)SO(4), the D form collapses into a pre-molten globule (P) form, whose volume exceeds that of the native (N) state by only 20%, and which lacks any specific structure as determined by far- and near-UV circular dichroism. At higher concentrations of Na(2)SO(4), the P form transforms into the molten globule (B) form in a highly noncooperative transition populated by an ensemble of at least two intermediates. The B form is a dry molten globule in which water is excluded from the core, and in which secondary structure develops to 65% and tertiary contacts develop to 40%, relative to that of the native protein. Kinetic refolding experiments carried out at pH 7 and at high Na(2)SO(4) concentrations, in which the rate of folding of the D form to the N state is compared to that of the B form to the N state, indicate conclusively that the B form is a productive intermediate that forms on the direct pathway of folding from the D form to the N state.     

Compactness of the kinetic molten globule of bovine alpha-lactalbumin: a dynamic light scattering study.

During folding of globular proteins, the molten globule state was observed as an equilibrium intermediate under mildly denaturing conditions as well as a transient intermediate in kinetic refolding experiments. While the high compactness of the equilibrium intermediate of alpha-lactalbumin has been verified, direct measurements of the compactness of the kinetic intermediate have not been reported until now. Our dynamic light scattering measurements provide a complete set of the hydrodynamic dimensions of bovine alpha-lactalbumin in different conformational states, particularly in the kinetic molten globule state. The Stokes radii for the native, kinetic molten globule, equilibrium molten globule, and unfolded states are 1.91, 1.99, 2.08, and 2.46 nm, respectively. Therefore, the kinetic intermediate appears to be even more compact than its equilibrium counterpart. Remarkable differences in the concentration dependence of the Stokes radius exist revealing strong attractive but repulsive intermolecular interactions in the kinetic and equilibrium molten globule states, respectively. This underlines the importance of extrapolation to zero protein concentration in measurements of the molecular compactness.     

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

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

Solution structure of native proteins with irregular folds from Raman optical activity

Raman optical activity (ROA) spectra have been measured for the proteins hen phosvitin, yeast invertase, bovine alpha-casein, soybean Bowman-Birk protease inhibitor, and rabbit Cd(7)-metallothionein, all of which have irregular folds in the native state. The results show that ROA is able to distinguish between two types of disorder. Specifically, invertase, alpha-casein, the Bowman-Birk inhibitor, and metallothionein appear to possess a "static" type of disorder similar to that in disordered states of poly(L-lysine) and poly(L-glutamic acid); whereas phosvitin appears to possess a more "dynamic" type of disorder similar to that in reduced (unfolded) lysozyme and ribonuclease A and also in molten globule protein states. In the delimiting cases, static disorder corresponds to that found in loops and turns within native proteins with well-defined tertiary folds that contain sequences of residues with fixed but nonrepetitive phi,psi angles; and dynamic disorder corresponds to that envisaged for the model random coil in which there is a distribution of Ramachandran phi,psi angles for each amino acid residue, giving rise to an ensemble of interconverting conformers. In both cases there is a propensity for the phi,psi angles to correspond to the alpha, beta and poly(L-proline) II (PPII) regions of the Ramachandran surface, as in native proteins with well-defined tertiary folds. Our results suggest that, with the exception of invertase and metallothionein, an important conformational element present in the polypeptide and protein states supporting the static type of disorder is that of the PPII helix. Long sequences of relatively unconstrained PPII helix, as in alpha-casein, may impart a plastic (rheomorphic) character to the structure.