Structural and dynamical properties of a partially unfolded Fe4S4 protein: role of the cofactor in protein folding

Heteronuclear multidimensional NMR spectroscopy was used to investigate in detail the structural and dynamical properties of a partially unfolded intermediate of the reduced high-potential iron-sulfur protein (HiPIP) from Chromatium vinosum present in 4 M guanidinium chloride solution. After an extensive assignment of 15N and 1H resonances, NOE data, proton longitudinal relaxation times, and 3JHNHalpha coupling constants as well as 15N relaxation parameters (T1, T2, T1rho, and 1H-15N NOE) were obtained and used to build a structural model of the intermediate. The Fe4S4 cluster of the HiPIP plays a decisive role in determining the resulting structure, which is random in the N-terminal half of the protein and partially organized in the loops between the cysteines bound to the cluster. Consistent with the structural data, the backbone mobility is typical of folded proteins in the regions where there are elements of structure and increases with the structural indetermination.     

Two partially unfolded states of Torpedo californica acetylcholinesterase.

Chemical modification with sulfhydryl reagents of the single, nonconserved cysteine residue Cys231 in each subunit of a disulfide-linked dimer of Torpedo californica acetylcholinesterase produces a partially unfolded inactive state. Another partially unfolded state can be obtained by exposure of the enzyme to 1-2 M guanidine hydrochloride. Both these states display several important features of a molten globule, but differ in their spectroscopic (CD, intrinsic fluorescence) and hydrodynamic (Stokes radii) characteristics. With reversal of chemical modification of the former state or removal of denaturant from the latter, both states retain their physiochemical characteristics. Thus, acetylcholinesterase can exist in two molten globule states, both of which are long-lived under physiologic conditions without aggregating, and without either intraconverting or reverting to the native state. Both states undergo spontaneous intramolecular thioldisulfide exchange, implying that they are flexible. As revealed by differential scanning calorimetry, the state produced by chemical modification lacks any heat capacity peak, presumably due to aggregation during scanning, whereas the state produced by guanidine hydrochloride unfolds as a single cooperative unit, thermal transition being completely reversible. Sucrose gradient centrifugation reveals that reduction of the interchain disulfide of the native acetylcholinesterase dimer converts it to monomers, whereas, after such reduction, the two subunits remain completely associated in the partially unfolded state generated by guanidine hydrochloride, and partially associated in that produced by chemical modification. It is suggested that a novel hydrophobic core, generated across the subunit interfaces, is responsible for this noncovalent association. Transition from the unfolded state generated by chemical modification to that produced by guanidine hydrochloride is observed only in the presence of the denaturant, yielding, on extrapolation to zero guanidine hydrochloride, a high free energy barrier (ca. 23.8 kcal/mol) separating these two flexible, partially unfolded states.     

The role of hydration in protein stability: comparison of the cold and heat unfolded states of Yfh1

Protein unfolding occurs at both low and high temperatures, although in most cases, only the high-temperature transition can be experimentally studied. A pressing question is how much the low- and high-temperature denatured states, although thermodynamically equivalent, are structurally and kinetically similar. We have combined experimental and computational approaches to compare the high- and low-temperature unfolded states of Yfh1, a natural protein that, at physiologic pH, undergoes cold and heat denaturation around 0 °C and 40 °C without the help of ad hoc destabilization. We observe that the two denatured states have similar but not identical residual secondary structures, different kinetics and compactness and a remarkably different degree of hydration. We use molecular dynamics simulations to rationalize the role of solvation and its effect on protein stability.     

Identification of rare partially unfolded states in equilibrium with the native conformation in an all beta-barrel protein

Human acidic fibroblast growth factor 1 (hFGF-1) is an all beta-barrel protein, and the secondary structural elements in the protein include 12 antiparallel beta-strands arranged into a beta-trefoil fold. In the present study, we investigate the stability of hFGF-1 by hydrogen-deuterium exchange as a function of urea concentration. Urea-induced equilibrium unfolding of hFGF-1 monitored by fluorescence and CD spectroscopy suggests that the protein unfolds by a two-state (native to denatured) mechanism. Hydrogen exchange in hFGF-1, under the experimental conditions used, occurs by the EX2 mechanism. In contrast to the equilibrium unfolding events monitored by optical probes, native state hydrogen exchange data show that the beta-trefoil architecture of hFGF-1 does not behave as a single cooperative unit. There are at least two structurally independent units with differing stabilities in hFGF-1. Beta-strands I, II, III, VI, VII, X, XI, and XII fit into the global unfolding isotherm. By contrast, residues in beta-strands IV, V, VIII, and IX exchange by the subfolding isotherm and could be responsible for the occurrence of high-energy partially unfolded state(s) in hFGF-1. There appears to be a broad continuum of stabilities among the four beta-strands (beta-strands IV, V, VIII, and IX) constituting the subglobal folding unit. The slow exchanging residues in hFGF-1 do not represent the folding nucleus of the protein.     

Congo red populates partially unfolded states of an amyloidogenic protein to enhance aggregation and amyloid fibril formation

Congo red (CR) has been reported to inhibit or enhance amyloid fibril formation by several proteins. To gain insight into the mechanism(s) for these apparently paradoxical effects, we studied as a model amyloidogenic protein, a dimeric immunoglobulin light chain variable domain. With a range of molar ratios of CR, i.e. r = [CR]/[protein dimer], we investigated the aggregation kinetics, conformation, hydrogen-deuterium exchange, and thermal stability of the protein. In addition, we used isothermal titration calorimetry to characterize the thermodynamics of CR binding to the protein. During incubation at 37 degrees C or during thermal scanning, with CR at r = 0.3, 1.3, and 4.8, protein aggregation was greatly accelerated compared with that measured in the absence of the dye. In contrast, with CR at r = 8.8, protein unfolding was favored over aggregation. The aggregates formed with CR at r = 0 or 0.3 were typical amyloid fibrils, but mixtures of amyloid fibrils and amorphous aggregates were formed at r = 1.3 and 4.8. CR decreased the apparent thermal unfolding temperature of the protein. Furthermore, CR perturbed the tertiary structure of the protein without significantly altering its secondary structure. Consistent with this result, CR also increased the rate of hydrogen-deuterium exchange by the protein. Isothermal titration calorimetry showed that CR binding to the protein was enthalpically driven, indicating that binding was mainly the result of electrostatic interactions. Overall, these results demonstrate that at low concentrations, CR binding to the protein favors a structurally perturbed, aggregation-competent species, resulting in acceleration of fibril formation. At high CR concentration, protein unfolding is favored over aggregation, and fibril formation is inhibited. Because low concentrations of CR can promote amyloid fibril formation, the therapeutic utility of this compound or its analogs to inhibit amyloidoses is questionable.     

Side-chain conformational entropy in protein unfolded states

The largest force disfavoring the folding of a protein is the loss of conformational entropy. A large contribution to this entropy loss is due to the side-chains, which are restricted, although not immobilized, in the folded protein. In order to accurately estimate the loss of side-chain conformational entropy that occurs upon folding it is necessary to have accurate estimates of the amount of entropy possessed by side-chains in the ensemble of unfolded states. A new scale of side-chain conformational entropies is presented here. This scale was derived from Monte Carlo computer simulations of small peptide models. It is demonstrated that the entropies are independent of host peptide length. This new scale has the advantage over previous scales of being more precise with low standard errors. Better estimates are obtained for long (e.g., Arg and Lys) and rare (e.g., Trp and Met) side-chains. Excellent agreement with previous side-chain entropy scales is achieved, indicating that further advancements in accuracy are likely to be small at best. Strikingly, longer side-chains are found to possess a smaller fraction of the theoretical maximum entropy available than short side-chains. This indicates that rotations about torsions after chi(2) are significantly affected by side-chain interactions with the polypeptide backbone. This finding invalidates previous assumptions about side-chain-backbone interactions. Proteins 2000;40:443-450.     

Cellular tolerance of prion protein PrP in yeast involves proteolysis and the unfolded protein response

Secretory proteins undergo a stringent quality control process in the endoplasmic reticulum (ER). Misfolded ER proteins are returned to the cytosol and destroyed by the proteasome. Prion protein PrP is degraded by the proteasome in mammalian cells. However, the significance of proteolysis on PrP-induced cell death is controversial. Moreover, the specific pathway involved in PrP degradation remains unknown. Here, we demonstrate that the unglycosylated form of human PrP is subjected to the ER-associated protein degradation (ERAD) process in the yeast Saccharomyces cerevisiae. We also show that unglycosylated PrP is degraded by the Hrd1-Hrd3 pathway. Accumulation of misfolded proteins triggers the unfolded protein response (UPR), which promotes substrate refolding. Interestingly, we find that the expression of PrP leads to growth impairment in cells deficient in UPR and ERAD. These findings raise the possibility that decreased UPR activity and proteolysis may contribute to the pathogenesis of some prion-related diseases.     

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.     

Populating partially unfolded forms by hydrogen exchange-directed protein engineering

The native-state hydrogen exchange of a redesigned apocytochrome b(562) suggests that at least two partially unfolded forms (PUFs) exist for this four-helix bundle protein under native conditions. The more stable PUF has the N-terminal helix unfolded. To verify the conclusion further and obtain more detailed structural information about this PUF, five hydrophobic core residues in the N-terminal helix were mutated to Gly and Asp to destabilize the native state selectively and populate the PUF for structural studies. The secondary structure and the backbone dynamics of this mutant were characterized using multidimensional NMR. Consistent with the prediction, the N-terminal region of the mutant was found to be unfolded while other parts of the proteins remained folded. These results suggest that native-state hydrogen exchange-directed protein engineering can be a useful approach to populating partially unfolded forms for detailed structural studies.     

Selective stabilization of a partially unfolded protein by a metabolite

When proteins fold in vivo, the intermediates that exist transiently on their folding pathways are exposed to the potential interactions with a plethora of metabolites within the cell. However, these potential interactions are commonly ignored. Here, we report a case in which a ubiquitous metabolite interacts selectively with a nonnative conformation of a protein and facilitates protein folding and unfolding process. From our previous proteomics study, we have discovered that Escherichia coli glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which is not known to bind ATP under native conditions, is apparently destabilized in the presence of a physiological concentration of ATP. To decipher the origin of this surprising effect, we investigated the thermodynamics and kinetics of folding and unfolding of GAPDH in the presence of ATP. Equilibrium unfolding of the protein in urea showed that a partially unfolded equilibrium intermediate accumulates in the presence of ATP. This intermediate has a quaternary structure distinct from the native protein. Also, ATP significantly accelerates the unfolding of GAPDH by selectively stabilizing a transition state that is distinct from the native state of the protein. Moreover, ATP also significantly accelerates the folding of GAPDH. These results demonstrate that ATP interacts specifically with a partially unfolded form of GAPDH and affects the kinetics of folding and unfolding of this protein. This unusual effect of ATP on the folding of GAPDH implies that endogenous metabolites may facilitate protein folding in vivo by interacting with partially unfolded intermediates.     

Partial molar volumes and adiabatic compressibilities of unfolded protein states

We determined the partial molar volumes, V degrees , and adiabatic compressibilities, K degrees (S), of N-acetyl amino acids with neutralized carboxyl termini, N-acetyl amino acid amides, and N-acetyl amino acid methylamides between 18 and 55 degrees C. The individual compounds in the three classes have been selected so as to collectively cover the 20 naturally occurring amino acid side chains. We interpret our experimental results in terms of the volumetric contributions and hydration properties of individual amino acid side chains and their constituent atomic groups. We also conducted pH-dependent densimetric and acoustic measurements to determine changes in volume and compressibility accompanying protonation of the aspartic acid, glutamic acid, histidine, lysine, and arginine side chains. We use our resulting data to develop an additive scheme for calculating the partial molar (specific) volume and adiabatic compressibility of fully extended polypeptide chains as a function of pH and temperature. We discuss the differences and similarities between our proposed scheme and the reported additive approaches. We compare our calculated volumetric characteristics of the fully extended conformations of apocytochrome c and apomyoglobin with the experimental values measured in water (for apocytochrome c) or acidic pH (for apomyoglobin). At these respective experimental conditions, the two proteins are unfolded. However, the comparison between the calculated and experimental volumetric characteristics suggests that neither apocytochrome c nor apomyoglobin are fully unfolded and retain a sizeable core of solvent-inaccessible groups.     

Structural characterization of unfolded states of apomyoglobin using residual dipolar couplings

The conformational propensities of unfolded states of apomyoglobin have been investigated by measurement of residual dipolar couplings between (15)N and (1)H in backbone amide groups. Weak alignment of apomyoglobin in acid and urea-unfolded states was induced with both stretched and compressed polyacrylamide gels. In 8 M urea solution at pH 2.3, conditions under which apomyoglobin contains no detectable secondary or tertiary structure, significant residual dipolar couplings of uniform sign were observed for all residues. At pH 2.3 in the absence of urea, a change in the magnitude and/or sign of the residual dipolar couplings occurs in local regions of the polypeptide where there is a high propensity for helical secondary structure. These results are interpreted on the basis of the statistical properties of the unfolded polypeptide chain, viewed as a polymer of statistical segments. For a folded protein, the magnitude and sign of the residual dipolar couplings depend on the orientation of each bond vector relative to the alignment tensor of the entire molecule, which reorients as a single entity. For unfolded proteins, there is no global alignment tensor; instead, residual dipolar couplings are attributed to alignment of the statistical segments or of transient elements of secondary structure. For apomyoglobin in 8 M urea, the backbone is highly extended, with phi and psi dihedral angles favoring the beta or P(II) regions. Each statistical segment has a highly anisotropic shape, with the N-H bond vectors approximately perpendicular to the long axis, and becomes weakly aligned in the anisotropic environment of the strained acrylamide gels. Local regions of enhanced flexibility or chain compaction are characterized by a decrease in the magnitude of the residual dipolar couplings. The formation of a small population of helical structure in the acid-denatured state of apomyoglobin leads to a change in sign of the residual dipolar couplings in local regions of the polypeptide; the population of helix estimated from the residual dipolar couplings is in excellent agreement with that determined from chemical shifts. The alignment model described here for apomyoglobin can also explain the pattern of residual dipolar couplings reported previously for denatured states of staphylococcal nuclease and other proteins. In conjunction with other NMR experiments, residual dipolar couplings can provide valuable insights into the dynamic conformational propensities of unfolded and partly folded states of proteins and thereby help to chart the upper reaches of the folding landscape.     

Definable equilibrium states in the folding of human prion protein

The role of conformational intermediates in the conversion of prion protein from its normal cellular form (PrP(C)) to the disease-associated scrapie form (PrP(Sc)) remains unknown. To look for such intermediates in equilibrium conditions, we have examined the unfolding transitions of PrP(C), primarily using the chemical denaturant guanidine hydrochloride (GuHCl). When the protein conformation is assessed by NMR, there is a gradual shift of NMR signals in the regions between residues 125-146 and 186-196. The denaturant dependence of these shifts shows that in aqueous solution the native and locally unfolded conformations are both significantly populated. Following this shift, there is the major unfolding transition to generate a substantially unfolded population. However, analysis of NMR chemical shift and intensity changes shows that there is persistent structure in the molecule well beyond this major cooperative unfolding transition. Residual structure within this state is extensive and encompasses the majority of the secondary structure elements found in the native state of the protein.     

Probing protein hydration and conformational states in solution.

The addition of polyethylene glycol (PEG), of various molecular weights, to solutions bathing yeast hexokinase increases the affinity of the enzyme for its substrate glucose. The results can be interpreted on the basis that PEG acts directly on the protein or indirectly through water activity. The nature of the effects suggests to us that PEG's action is indirect. Interpretation of the results as an osmotic effect yields a decrease in the number of water molecules, delta Nw, associated with the glucose binding reaction. delta Nw is the difference in the number of PEG-inaccessible water molecules between the glucose-bound and glucose-free conformations of hexokinase. At low PEG concentrations, delta Nw increases from 50 to 326 with increasing MW of the PEG from 300 to 1000, and then remains constant for MW-PEG up to 10,000. This suggests that up to MW 1000, solutes of increasing size are excluded from ever larger aqueous compartments around the protein. Three hundred and twenty-six waters is larger than is estimated from modeling solvent volumes around the crystal structures of the two hexokinase conformations. For PEGs of MW > 1000, delta Nw falls from 326 to about 25 waters with increasing PEG concentration, i.e., PEG alone appears to "dehydrate" the unbound conformation of hexokinase in solution. Remarkably, the osmotic work of this dehydration would be on the order of only one k T per hexokinase molecule. We conclude that under thermal fluctuations, hexokinase in solution has a conformational flexibility that explores a wide range of hydration states not seen in the crystal structure.     

Structural requirements for efficient prion protein conversion

To understand why cross species infection of prion disease often results in inefficient transmission and reduced protein conversion, most research has focussed on defining the effect of variations in PrP primary structures, including sequence compatibility of substrate and seed. By contrast, little research has been aimed at investigating structural differences between different variants of PrP^C and secondary structural requirements for efficient conversion. This is despite a clear role for molecular chaperones in formation of prions in non-mammalian systems, indicating the importance of secondary/tertiary structure during the conversion process. Recent data from our laboratory on the cellular location of disease-specific prion cofactors supports the critical role of specific secondary structural motifs and the stability of these motifs in determining the efficiency of disease-specific prion protein conversion. In this paper we summarise our recent results and build on the hypothesis previously suggested by Wuthrich and colleagues, that stability of certain regions of the prion protein is crucial for protein conversion to abnormal isoforms in vivo. It is suggested that one role for molecular co-factors in the conversion process is to stabilise PrP^C structure in a form that is amenable for conversion to PrP^Sc.     

Structural properties of prion protein protofibrils and fibrils: an experimental assessment of atomic models

Decades after the prion protein was implicated in transmissible spongiform encephalopathies, the structure of its toxic isoform and its mechanism of toxicity remain unknown. By gathering available experimental data, albeit low resolution, a few pieces of the prion puzzle can be put in place. Currently, there are two fundamentally different models of a prion protofibril. One has its building blocks derived from a molecular dynamics simulation of the prion protein under amyloidogenic conditions, termed the spiral model. The other model was constructed by threading a portion of the prion sequence through a beta-helical structure from the Protein Data Bank. Here we compare and contrast these models with respect to all of the available experimental information, including electron micrographs, symmetries, secondary structure, oligomerization interfaces, enzymatic digestion, epitope exposure, and disaggregation profiles. Much of this information was not available when the two models were introduced. Overall, we find that the spiral model is consistent with all of the experimental results. In contrast, it is difficult to reconcile several of the experimental observables with the beta-helix model. While the experimental constraints are of low resolution, in bringing together the previously disconnected experiments, we have developed a clearer picture of prion aggregates. Both the improved characterization of prion aggregates and the existing atomic models can be used to devise further experiments to better elucidate the misfolding pathway and the structure of prion protofibrils.     

Ligand binding to a high‐energy partially unfolded protein

The conformational energy landscape of a protein determines populations of all possible conformations of the protein and also determines the kinetics of the conversion between the conformations. Interaction with ligands influences the conformational energy landscapes of proteins and shifts populations of proteins in different conformational states. To investigate the effect of ligand binding on partial unfolding of a protein, we use Escherichia coli dihydrofolate reductase (DHFR) and its functional ligand NADP⁺ as a model system. We previously identified a partially unfolded form of DHFR that is populated under native conditions. In this report, we determined the free energy for partial unfolding of DHFR at varying concentrations of NADP⁺ and found that NADP⁺ binds to the partially unfolded form as well as the native form. DHFR unfolds partially without releasing the ligand, though the binding affinity for NADP⁺ is diminished upon partial unfolding. Based on known crystallographic structures of NADP⁺‐bound DHFR and the model of the partially unfolded protein we previously determined, we propose that the adenosine‐binding domain of DHFR remains folded in the partially unfolded form and interacts with the adenosine moiety of NADP⁺. Our result demonstrates that ligand binding may affect the conformational free energy of not only native forms but also high‐energy non‐native forms.     

Structural changes in a hydrophobic domain of the prion protein induced by hydration and by ala-->Val and pro-->Leu substitutions

X-ray diffraction was used to study the structure of assemblies formed by synthetic peptide fragments of the prion protein (PrP) that include the hydrophobic domain implicated in the Gerstmann-Str?ussler-Scheinker (GSS) mutation (P102L). The effects of hydration on polypeptide assembly and of Ala-->Val substitutions in the hydrophobic domain were characterized. Synthetic peptides included: (i) Syrian hamster (SHa) hydrophobic core, SHa106-122 (KTNMKHMAGAAAAGAVV); (ii) SHa104-122(3A-V), with A-->V mutations at 113, 115 and 118 (KPKTNMKHMVGVAAVGAVV); (iii) mouse (Mo) wild-type sequence of the N-terminal hydrophobic domain, Mo89-143WT; and (iv) the same mouse sequence with leucine substitution for proline at residue number 101, Mo89-143(P101L). Samples of SHa106-122 that formed assemblies while drying under ambient conditions showed X-ray patterns indicative of 33 A thick slab-like structures having extensive H-bonding and intersheet stacking. By contrast, lyophilized peptide that was equilibrated against 100 % relative humidity showed assemblies with only a few layers of beta-sheets. The Ala-->Val substitutions in SHa104-122 and Mo89-143(P101L) resulted in the formation of 40 A wide, cross-beta fibrils. Observation of similar size beta-sheet fibrils formed by peptides SHa104-122(3A-V) and the longer Mo89-143(P101L) supports the notion that the hydrophobic sequence forms a template or core that promotes the beta-folding of the longer peptide. The substitution of amino acids in the mutants, e.g. 3A-->V and P101L, enhances the folding of the peptide into compact structural units, significantly enhancing the formation of the extensive beta-sheet fibrils.     

Conformational plasticity of cryptolepain: accumulation of partially unfolded states in denaturants induced equilibrium unfolding

pH and chemical denaturant dependent conformational changes of a serine protease cryptolepain from Cryptolepis buchanani are presented in this paper. Activity measurements, near UV, far UV CD, fluorescence emission spectroscopy, and ANS binding studies have been carried out to understand the folding mechanism of the protein in the presence of denaturants. pH and chemical denaturants have a marked effect on the stability, structure, and function of many globular proteins due to their ability to influence the electrostatic interactions. The preliminary biophysical study on cryptolepain shows that major elements of secondary structure are beta-sheets. Under neutral conditions the enzyme was stable in urea while GuHCl-induced equilibrium unfolding was cooperative. Cryptolepain shows little ANS binding even under neutral conditions due to more hydrophobicity of beta-sheets. Multiple intermediates were populated during the pH-induced unfolding of cryptolepain. Temperature-induced denaturation of cryptolepain in the molten globule like state is non-cooperative, contrary to the cooperativity seen with the native protein, suggesting the presence of two parts, possibly domains, in the molecular structure of cryptolepain, with different stability that unfolds in steps. Interestingly, the GuHCl-induced unfolding of A state (molten globule state) of cryptolepain is unique, as lower concentration of denaturant, not only induces structure but also facilitate transition from one molten globule like state (MG(1)) into another (MG(2)). The increase of pH drives the protein into alkaline denatured state characterized by the absence of any ANS binding. GuHCl- and urea-induced unfolding transition curves at pH 12.0 were non-coincidental indicating the presence of an intermediate in the unfolding pathway.