Urea-dependent unfolding of murine adenosine deaminase: sequential destabilization as measured by 19F NMR

Murine adenosine deaminase (mADA) is a 40 kDa (beta/alpha)(8)-barrel protein consisting of eight central beta-strands and eight peripheral alpha-helices containing four tryptophan residues. In this study, we investigated the urea-dependent behavior of the protein labeled with 6-fluorotryptophan (6-(19)F-Trp). The (19)F NMR spectrum of 6-(19)F-Trp-labeled mADA reveals four distinct resonances in the native state and three partly overlapped resonances in the unfolded state. The resonances were assigned unambiguously by site-directed mutagenesis. Equilibrium unfolding of 6-(19)F-Trp-labeled mADA was monitored using (19)F NMR based on these assignments. The changes in intensity of folded and unfolded resonances as a function of urea concentration show transition midpoints consistent with data observed by far-UV CD and fluorescence spectroscopy, indicating that conformational changes in mADA during urea unfolding can be followed by (19)F NMR. Chemical shifts of the (19)F resonances exhibited different changes between 1.0 and 6.0 M urea, indicating that local structures around 6-(19)F-Trp residues change differently. The urea-induced changes in local structure around four 6-(19)F-Trp residues of mADA were analyzed on the basis of the tertiary structure and chemical shifts of folded resonances. The results reveal that different local regions in mADA have different urea-dependent behavior, and that local regions of mADA change sequentially from native to intermediate topologies on the unfolding pathway.     

Temperature dependence of the folding and unfolding kinetics of the GCN4 leucine zipper via 13C(alpha)-NMR

Studies by one-dimensional NMR are reported on the interconversion of folded and unfolded forms of the GCN4 leucine zipper in neutral saline buffer. The peptide bears 99% 13C(alpha) labels at three sites: V9, L12, and G31. Time-domain 13C(alpha)-NMR spectra are interpreted by global Bayesian lineshape analysis to extract the rate constants for both unfolding and folding as functions of temperature in the range 47-71 degrees C. The data are well fit by the assumption that the same rate constants apply at each labeled site, confirming that only two conformational states need be considered. Results show that 1) both processes require a free energy of activation; 2) unfolding is kinetically enthalpy-opposed and entropy-driven, while folding is the opposite; and 3) the transition state dimer ensemble averages approximately 40% helical. The activation parameters for unfolding, derived from NMR data at the elevated temperatures where both conformations are populated, lead to estimates of the rate constant at low temperatures (5-15 degrees C) that agree with extant values determined by stopped-flow CD via dilution from denaturing media. However, the corresponding estimated values for the folding rate constant are larger by two to three orders of magnitude than those obtained by stopped flow. We propose that this apparent disagreement is caused by the necessity, in the stopped-flow experiment, for initiation of new helices as the highly denaturant-unfolded molecule adjusts to the newly created benign solvent conditions. This must reduce the success rate of collisions in producing the folded molecule. In the NMR determinations, however, the unfolded chains always have a small, but essential, helix content that makes such initiation unnecessary. Support for this hypothesis is adduced from recent extant experiments on the helix-coil transition in single-chain helical peptides and from demonstration that the folding rate constants for coiled coils, as obtained by stopped flow, are influenced by the nature of the denaturant used.     

No cofactor effect on equilibrium unfolding of Desulfovibrio desulfuricans flavodoxin

Flavodoxins are proteins with an alpha/beta doubly wound topology that mediate electron transfer through a non-covalently bound flavin mononucleotide (FMN). The FMN moiety binds strongly to folded flavodoxin (K(D)=0.1 nM, oxidized FMN). To study the effect of this organic cofactor on the conformational stability, we have characterized apo and holo forms of Desulfovibrio desulfuricans flavodoxin by GuHCl-induced denaturation. The unfolding reactions for both holo- and apo-flavodoxin are reversible. However, the unfolding curves monitored by far-UV circular dichroism and fluorescence spectroscopy do not coincide. For both apo- and holo-flavodoxin, a native-like intermediate (with altered tryptophan fluorescence but secondary structure as the folded form) is present at low GuHCl concentrations. There is no effect on the flavodoxin stability imposed by the presence of the FMN cofactor (DeltaG=20(+/-2) and 19(+/-1) kJ/mol for holo- and apo-flavodoxin, respectively). A thermodynamic cycle, connecting FMN binding to folded and unfolded flavodoxin with the unfolding free energies for apo- and holo-flavodoxin, suggests that the binding strength of FMN to unfolded flavodoxin must be very high (K(D)=0.2 nM). In agreement, we discovered that the FMN remains coordinated to the polypeptide upon unfolding.     

An engineered leucine zipper a position mutant with an unusual three-state unfolding pathway

The leucine zipper is a dimeric coiled-coil structural motif consisting of four to six heptad repeats, designated (abcdefg)(n). In the GCN4 leucine zipper, a position 16 in the third heptad is occupied by an Asn residue whereas the other a positions are Val residues. Recently, we have constructed variants of the GCN4 leucine zipper in which the a position Val residues were replaced by Ile. The folding and unfolding of the wild-type GCN4 leucine zipper and the Val to Ile variant both adhere to a simple two-state mechanism. In this study, another variant of the GCN4 leucine zipper was constructed by moving the single Asn residue from a position 16 to a position 9. This switch causes the thermal unfolding of the GCN4 leucine zipper to become three state. The unfolding pathway of this variant was determined by thermal denaturation, limited proteinase K digestion, and sedimentation equilibrium analysis. Our data are consistent with a model in which the variant first unfolds from its N terminus and changes the oligomerization specificity from a native dimer to a partially unfolded intermediate containing a mixture of dimers and trimers and then completely unfolds to unstructured monomers.     

Conformational specificity of the lac repressor coiled-coil tetramerization domain

Predictive understanding of how the folded, functional shape of a native protein is encoded in the linear sequence of its amino acid residues remains an unsolved challenge in modern structural biology. Antiparallel four-stranded coiled coils are relatively simple protein structures that embody a heptad sequence repeat and rich diversity for tertiary packing of alpha-helices. To explore specific sequence determinants of the lac repressor coiled-coil tetramerization domain, we have engineered a set of buried nonpolar side chains at the a-, d-, and e-positions into the hydrophobic interior of the dimeric GCN4 leucine zipper. Circular dichroism and equilibrium ultracentrifugation studies show that this core variant (GCN4-pAeLV) forms a stable tetrameric structure with a reversible and highly cooperative thermal unfolding transition. The X-ray crystal structure at 1.9 A reveals that GCN4-pAeLV is an antiparallel four-stranded coiled coil of the lac repressor type in which the a, d, and e side chains associate by means of combined knobs-against-knobs and knobs-into-holes packing with a characteristic interhelical offset of 0.25 heptad. Comparison of the side chain shape and packing in the antiparallel tetramers shows that the burial of alanine residues at the e positions between the neighboring helices of GCN4-pAeLV dictates both the antiparallel orientation and helix offset. This study fills in a gap in our knowledge of the determinants of structural specificity in antiparallel coiled coils and improves our understanding of how specific side chain packing forms the teritiary structure of a functional protein.     

CD and (13)C(alpha)-NMR studies of folding equilibria in a two-stranded coiled coil formed by residues 190-254 of alpha-tropomyosin

Synthesis and CD and (13)C(alpha)-NMR studies in a near-neutral saline buffer are reported for a 65-residue peptide ((190)Tm(254)) comprising residues 190-254 of the alpha-tropomyosin chain. CD on a version disulfide cross-linked via the N-terminal cysteine side chains indicates that this dimer is highly helical and melts near 48 degrees C. The CD is independent of peptide concentration, showing that association of (190)Tm(254) stops at the two-strand level. Similar studies on the reduced version show much lower helix content at low temperature, melting points below room temperature, and the expected concentration dependence. The observed melting temperature of the reduced peptide is far below (by 27 degrees C) that expected from an extant analysis of calorimetry data on parent tropomyosin that designates (190)Tm(254) as an independently melting "cooperative block." This disagreement and the pronounced nonadditivity seen when data for (190)Tm(254) are combined with extant data for other subsequences argue decisively against the concept of specific independently melting blocks within the tropomyosin chain. The data for (190)Tm(254) also serve to test recent ideas on the sequence determinants of structure and stability in coiled coils. Analysis shows that some ideas, such as the stabilizing effect of leucine in the d heptad position, find support, but others--such as the destabilizing effect of alanine in d, the dimer-disfavoring effect of beta-branching in d and its dimer-favoring effect in a, and the dimer-directing effect of asparagine in a--are more questionable in tropomyosin than in the leucine zipper coiled coils. (13)C(alpha)-NMR data at two labeled sites, L228(d) and V246(a), of (190)Tm(254) display well-separated resonances for folded and unfolded forms at each site, indicating that the transition is slow on the NMR time scale and thus demonstrating the possibility of obtaining thermodynamic and kinetic information on the transition at the residue level.     

A measure of conformational entropy change during thermal protein unfolding using neutron spectroscopy

Thermal unfolding of proteins at high temperatures is caused by a strong increase of the entropy change which lowers Gibbs free energy change of the unfolding transition (DeltaG(unf) = DeltaH - TDeltaS). The main contributions to entropy are the conformational entropy of the polypeptide chain itself and ordering of water molecules around hydrophobic side chains of the protein. To elucidate the role of conformational entropy upon thermal unfolding in more detail, conformational dynamics in the time regime of picoseconds was investigated with neutron spectroscopy. Confined internal structural fluctuations were analyzed for alpha-amylase in the folded and the unfolded state as a function of temperature. A strong difference in structural fluctuations between the folded and the unfolded state was observed at 30 degrees C, which increased even more with rising temperatures. A simple analytical model was used to quantify the differences of the conformational space explored by the observed protein dynamics for the folded and unfolded state. Conformational entropy changes, calculated on the basis of the applied model, show a significant increase upon heating. In contrast to indirect estimates, which proposed a temperature independent conformational entropy change, the measurements presented here, demonstrated that the conformational entropy change increases with rising temperature and therefore contributes to thermal unfolding.     

Thermodynamics of alligator metmyoglobin unfolding

The conformational free energy of alligator metmyoglobin was examined over a pH range of 4.4-8.0, a temperature range of 18-48 degrees C, and a guanidinium chloride concentration of 0-2.3 M. For isothermal unfolding at 25 degrees C essentially the same value was obtained for the conformational free energy from all the data: 7.0 +/- 0.5 kcal/mol. The cooperativity of the unfolding with respect to denaturant is considerably less than for mammalian myoglobins. On unfolding three to four side chains with a pKa of 6.3 in the unfolded protein are protonated instead of the six expected. The temperature at which delta H (unfolding) is zero is much lower than for previously characterized myoglobins. Alligator metmyoglobin, considerably less stable than other previously characterized myoglobins, may not be as compactly folded.     

Localized unfolding of collagen explains collagenase cleavage near imino-poor sites

Collagenases cleave all three chains of type III collagen at specific sites characterized by a Gly-Leu or a Gly-Ile bond that is upstream from an imino acid-poor region. Molecular dynamics trajectories were used to calculate the free energy of unfolding for collagen-like model peptides. The free energy profiles suggest that such imino-poor regions can adopt a low-energy, partially unfolded state where one of the peptide chains forms a solvent-exposed loop. The results are consistent with a model for collagenase cleavage where partial unfolding of imino-poor regions enables collagenases to gain access to their cleavage sites.     

Local unfolding of Cu, Zn superoxide dismutase monomer determines the morphology of fibrillar aggregates

Aggregation of Cu, Zn superoxide dismutase (SOD1) is often found in amyotrophic lateral sclerosis patients. The fibrillar aggregates formed by wild type and various disease-associated mutants have recently been found to have distinct cores and morphologies. Previous computational and experimental studies of wild-type SOD1 suggest that the apo-monomer, highly aggregation prone, displays substantial local unfolding dynamics. The residual folded structure of locally unfolded apoSOD1 corresponds to peptide segments forming the aggregation core as identified by a combination of proteolysis and mass spectroscopy. Therefore, we hypothesize that the destabilization of apoSOD1 caused by various mutations leads to distinct local unfolding dynamics. The partially unfolded structure, exposing the hydrophobic core and backbone hydrogen bond donors and acceptors, is prone to aggregate. The peptide segments in the residual folded structures form the "building block" for aggregation, which in turn determines the morphology of the aggregates. To test this hypothesis, we apply a multiscale simulation approach to study the aggregation of three typical SOD1 variants: wild type, G37R, and I149T. Each of these SOD1 variants has distinct peptide segments forming the core structure and features different aggregate morphologies. We perform atomistic molecular dynamics simulations to study the conformational dynamics of apoSOD1 monomer and coarse-grained molecular dynamics simulations to study the aggregation of partially unfolded SOD1 monomers. Our computational studies of monomer local unfolding and the aggregation of different SOD1 variants are consistent with experiments, supporting the hypothesis of the formation of aggregation "building blocks" via apo-monomer local unfolding as the mechanism of SOD1 fibrillar aggregation.     

Methanol-stabilized intermediates in the thermal unfolding of ribonuclease A: Characterization by 1H nuclear magnetic resonance

The thermal unfolding of ribonuclease A has been studied in solutions of 25, 35 and 50% methanol ($\text{v}\text{v}$), using 360 MHz proton magnetic resonance spectroscopy. Several observations indicate that the native structure of the protein in methanol cryosolvents is very similar to that in aqueous solution. A detailed analysis of the unfolding process has been made using the C-2 protons of the imidazole side-chains of the four histidine residues. As denaturation proceeds new resonances appear, whose chemical shifts correspond to neither native nor fully unfolded species. These have been assigned to particular His residues by selective deuteration studies. The thermal denaturation transitions reveal a multiphasic process in each of the solvents, and become less co-operative with increasing concentrations of methanol. The denaturation is fully reversible with no evidence of hysteresis.The new resonances that appear during the unfolding process are attributed to partially folded species, which are stabilized by the presence of the relatively hydrophobic methanol. Based on the temperature dependence of the chemical shifts and the relative areas of the various resonances, a detailed sequence of events has been proposed to describe the unfolding process. Key features include the initial general loosening of the two domains, the subsequent movement of the upper S-peptide region (residues 13 to 25) away from the main body of the protein, followed by partial separation of the sheet structure and full exposure of the N-terminal helix, leading to complete separation of the “winged domains”, and ultimately the loss of the residual sheet and helix structure.     

Kinetics of nucleosome unfolding at low ionic strength

The kinetics of a conformational change which occurs in nucleosome core particles at about 1 mM ionic strength have been studied by observing changes in the fluorescence of labeled histone H3. The unfolding reaction is intramolecular since no concentration dependence is observed. However, the kinetics are unexpectedly complicated and reveal evidence of at least three relaxation times. It is possible to fit the kinetics observed under several conditions to a consistent four-state cyclic mechanism in which folded and unfolded forms can inter-convert by two parallel pathways, each involving a distinct intermediate. While the data are not sufficient to establish this mechanism as a unique choice, they exclude many simpler possibilities. The cyclic mechanism is quite reasonable in view of what is currently known about the structures of the folded and unfolded forms.     

The equilibrium unfolding pathway of a (beta/alpha)8 barrel

The (beta/alpha)(8) barrel is the most commonly occurring fold among enzymes. A key step towards rationally engineering (beta/alpha)(8) barrel proteins is to understand their underlying structural organization and folding energetics. Using misincorporation proton-alkyl exchange (MPAX), a new tool for solution structural studies of large proteins, we have performed a native-state exchange analysis of the prototypical (beta/alpha)(8) barrel triosephosphate isomerase. Three cooperatively unfolding subdomains within the structure are identified, as well as two partially unfolded forms of the protein. The C-terminal domain coincides with domains reported to exist in four other (beta/alpha)(8) barrels, but the two N-terminal domains have not been observed previously. These partially unfolded forms may represent sequential intermediates on the folding pathway of triosephosphate isomerase. The methods reported here should be applicable to a variety of other biological problems involving protein conformational changes.     

Quaternary structure of aldolase leads to differences in its folding and unfolding intermediates

Pulsed hydrogen exchange mass spectrometry has been used to investigate folding of rabbit muscle aldolase, an alpha/beta-barrel protein exhibiting the classic TIM structure. Aldolase unfolded in GdHCl refolded as the denaturant concentration was reduced by dialysis. Samples withdrawn during dialysis were pulse-labeled with deuterium to identify unfolded regions in structural forms highly populated during the folding process. Intact, labeled aldolase was digested into fragments, which were analyzed by HPLC electrospray ionization mass spectrometry to detect the H/D exchange along the aldolase backbone. For some concentrations of GdHCl, bimodal distributions of deuterium were found for most peptic fragments, indicating that regions represented by these fragments were either unfolded or folded in the intact polypeptide prior to labeling. The extent of folding was determined from these mass spectra, as well as by CD (220 nm) and enzymatic activity. These results show that folding to the active form involves three domains and two intermediates. Approximately 110 residues fold to highly compact forms in each step. These results also show that each folding domain includes widely separated regions of the backbone. When compared with the results of a previous study of aldolase unfolding, these results show that the folding and unfolding domains include most of the same residues. However, three short segments change domains depending on whether the process is folding or unfolding. These changes are attributed to the very stable quaternary structure of rabbit muscle aldolase.     

Site-specific experiments on folding/unfolding of Jun coiled coils: thermodynamic and kinetic parameters from spin inversion transfer nuclear magnetic resonance at leucine-18

The 32-residue leucine zipper subsequence, called here Jun-lz, associates in benign media to form a parallel two-stranded coiled coil. Studies are reported of its thermal unfolding/folding transition by circular dichroism (CD) on samples of natural isotopic abundance and by both equilibrium and spin inversion transfer (SIT) nuclear magnetic resonance (NMR) on samples labeled at the leucine-18 alpha-carbon with 99% 13C. The data cover a wide range of temperature and concentration, and show that Jun-lz unfolds below room temperature, being far less stable than some other leucine zippers such as GCN4. 13C-NMR shows two well-separated resonances. We ascribe the upfield one to 13C spins on unfolded single chains and the downfield one to 13C spins on coiled-coil dimers. Their relative intensities provide a measure of the unfolding equilibrium constant. In SIT NMR, the recovery of the equilibrium magnetization after one resonance is inverted is modulated in part by the unfolding and folding rate constants, which are accessible from the data. Global Bayesian analysis of the equilibrium and SIT NMR data provide values for the standard enthalpy, entropy, and heat capacity of unfolding, and show the latter to be unusually large. The CD results are compatible with the NMR findings. Global Bayesian analysis of the SIT NMR data yields the corresponding activation parameters for unfolding and folding. The results show that both reaction directions are activated processes. Activation for unfolding is entropy driven, enthalpy opposed. Activation for folding is strongly enthalpy opposed and somewhat entropy opposed, falsifying the idea that the barrier for folding is solely due to a purely entropic search for properly registered partners. The activation heat capacity is much larger for folding, so almost the entire overall change is due to the folding direction. This latter finding, if it applies to GCN4 leucine zippers, clears up an extant apparent disagreement between folding rate constants for GCN4 as determined by chevron analysis and NMR in differing temperature regimes.     

pH-induced conformation changes and equilibrium unfolding in yeast iso-2 cytochrome c

The relationship between pH-induced conformational changes in iso-2 cytochrome c from Saccharomyces cerevisiae and the guanidine hydrochloride induced unfolding transition has been investigated. Comparison of equilibrium unfolding transitions at acid, neutral, and alkaline pH shows that stability toward guanidine hydrochloride denaturation is decreased at low pH but increased at high pH. In the acid range the decrease in stability of the folded protein is correlated with changes in the visible spectrum, which indicate conversion to a high-spin heme state--probably involving the loss of heme ligands. The increase in stability at high pH is correlated with a pH-induced conformational change with an apparent pK near 8. As in the case of homologous cytochromes c, this transition involves the loss of the 695-nm absorbance band with only minor changes in other optical parameters. For the unfolded protein, optical spectroscopy and 1H NMR spectroscopy are consistent with a random coil unfolded state in which amino acid side chains serve as (low-spin) heme ligands at both neutral and alkaline pH. However, the paramagnetic region of the proton NMR spectrum of unfolded iso-2 cytochrome c indicates a change in the (low-spin) heme-ligand complex at high pH. Apparently, the folded and unfolded states of the (inactive) alkaline form differ from the corresponding states of the less stable native protein.     

Site-specific thermodynamics and kinetics of a coiled-coil transition by spin inversion transfer NMR.

A 33-residue pseudo-wild-type GCN4 leucine zipper peptide is used to probe the equilibrium conformational population in proteins. 13Calpha-NMR shows that chain sites differ in structural content at a given temperature, and that two dimeric folded forms are evident at many sites. Spin inversion transfer experiments are reported bearing on the thermodynamics and kinetics of interconversion of the two dimeric folded forms (Fa <--> Fb) at the 13Calpha-labeled position L13. At each temperature, at conditions wherein the population of unfolded chains is quite small, inversion of the Fa spins via a tuned Gaussian pi-pulse is followed by a time interval (tau), interrogation, and recording of the free induction decay. Fifteen such inversions, with varying tau, provide the time course for recovery of equilibrium magnetization after inversion. Similar experiments follow inversion of the Fb spins. Re-equilibration is known to be modulated by four first-order rate constants: two (T1a(-1) and T1b(-1)) for spin-lattice relaxation intrinsic to the respective sites, and two (kab and kba) for the conformational change. All four follow from joint, Bayesian analysis of all the data at each temperature. The equilibrium constant at each temperature for this local transition, determined simply from the equilibrium relative magnetizations at Fa and Fb sites, agrees well with the kinetic ratio kab/kba. The standard Gibbs energies, enthalpy, and entropy follow. Activation parameters, both ways, are accessible from the rate constants and suggest a transition state with high Gibbs energy and enthalpy, but with entropy between those of Fa and Fb.     

Folding and unfolding of gammaTIM monomers and dimers

Kinetic simulations of the folding and unfolding of triosephosphate isomerase (TIM) from yeast were conducted using a single monomer gammaTIM polypeptide chain that folds as a monomer and two gammaTIM chains that fold to the native dimer structure. The basic protein model used was a minimalist Gō model using the native structure to determine attractive energies in the protein chain. For each simulation type--monomer unfolding, monomer refolding, dimer unfolding, and dimer refolding--thirty simulations were conducted, successfully capturing each reaction in full. Analysis of the simulations demonstrates four main conclusions. First, all four simulation types have a similar "folding order", i.e., they have similar structures in intermediate stages of folding between the unfolded and folded state. Second, despite this similarity, different intermediate stages are more or less populated in the four different simulations, with 1), no intermediates populated in monomer unfolding; 2), two intermediates populated with beta(2)-beta(4) and beta(1)-beta(5) regions folded in monomer refolding; 3), two intermediates populated with beta(2)-beta(3) and beta(2)-beta(4) regions folded in dimer unfolding; and 4), two intermediates populated with beta(1)-beta(5) and beta(1)-beta(5) + beta(6) + beta(7) + beta(8) regions folded in dimer refolding. Third, simulations demonstrate that dimer binding and unbinding can occur early in the folding process before complete monomer-chain folding. Fourth, excellent agreement is found between the simulations and MPAX (misincorporation proton alkyl exchange) experiments. In total, this agreement demonstrates that the computational Gō model is accurate for gammaTIM and that the energy landscape of gammaTIM appears funneled to the native state.