Mapping the interactions present in the transition state for unfolding/folding of FKBP12.

The structure of the transition state for folding/unfolding of the immunophilin FKBP12 has been characterised using a combination of protein engineering techniques, unfolding kinetics, and molecular dynamics simulations. A total of 34 mutations were made at sites throughout the protein to probe the extent of secondary and tertiary structure in the transition state. The transition state for folding is compact compared with the unfolded state, with an approximately 30 % increase in the native solvent-accessible surface area. All of the interactions are substantially weaker in the transition state, as probed by both experiment and molecular dynamics simulations. In contrast to some other proteins of this size, no element of structure is fully formed in the transition state; instead, the transition state is similar to that found for smaller, single-domain proteins, such as chymotrypsin inhibitor 2 and the SH3 domain from alpha-spectrin. For FKBP12, the central three strands of the beta-sheet, beta-strand 2, beta-strand 4 and beta-strand 5, comprise the most structured region of the transition state. In particular Val101, which is one of the most highly buried residues and located in the middle of the central beta-strand, makes approximately 60 % of its native interactions. The outer beta-strands and the ends of the central beta-strands are formed to a lesser degree. The short alpha-helix is largely unstructured in the transition state, as are the loops. The data are consistent with a nucleation-condensation model of folding, the nucleus of which is formed by side-chains within beta-strands 2, 4 and 5, and the C terminus of the alpha-helix. The precise residues involved in the nucleus differ in the two simulated transition state ensembles, but the interacting regions of the protein are conserved. These residues are distant in the primary sequence, demonstrating the importance of tertiary interactions in the transition state. The two independently derived transition state ensembles are structurally similar, which is consistent with a Bronsted analysis confirming that the transition state is an ensemble of states close in structure.     

Stabilisation of alpha-helices by site-directed mutagenesis reveals the importance of secondary structure in the transition state for acylphosphatase folding

The effects of stabilising mutations on the folding process of common-type acylphosphatase have been investigated. The mutations were designed to increase the helical propensity of the regions of the polypeptide chain corresponding to the two alpha-helices of the native protein. Various synthetic peptides incorporating the designed mutations were produced and their helical content estimated by circular dichroism. The most substantial increase in helical content is found for the peptide carrying five mutations in the second alpha-helix. Acylphosphatase variants containing the corresponding mutations display, to different extents, enhanced conformational stabilities as indicated by equilibrium urea denaturation experiments monitored by changes of intrinsic fluorescence. All the protein variants studied here refold with apparent two-state kinetics. Mutations in the first alpha-helix are responsible for a small increase in the refolding rate, accompanied by a marked decrease in the unfolding rate. On the other hand, multiple mutations in the second helix result in a considerable increase in the refolding rate without any significant effect on the unfolding rate. Addition of trifluoroethanol was found to accelerate the folding of the acylphosphatase variants, the extent of the acceleration being inversely proportional to the intrinsic rate of folding of the corresponding mutant. The trifluoroethanol-induced acceleration is far less marked for those variants whose alpha-helical structure is efficiently stabilised by amino acid replacements. This observation suggests that trifluoroethanol acts in a similar manner to the stabilising mutations in promoting native-like secondary structure. Analysis of the kinetic data indicates that the second helix is fully consolidated in the transition state for folding of acylphosphatase, whereas the first helix is only partially formed. These data suggest that the second helix is an important element in the folding process of the protein.     

Microsecond folding of the cold shock protein measured by a pressure-jump technique

A pressure-jump apparatus was employed in investigating the kinetics of protein unfolding and refolding. In the reaction cell, the pressure can be increased or decreased by 100-160 bar within 50-100 microseconds and then held constant. Thus, unfolding and refolding reactions in the time range from 70 microseconds to 70 s can be followed with this technique. Measurements are possible in the transition regions of thermally or denaturant-induced folding in a wide range of temperatures and solvent conditions. We used this pressure-jump method to determine the temperature dependence of the rate constants of unfolding and refolding of the cold shock protein of Bacillus subtilis and of three variants thereof with Phe --> Ala substitutions in the central beta-sheet region. For all variants, the change in heat capacity occurred in refolding between the unfolded and activated states, suggesting that the overall native-like character of the activated state of folding was not changed by the deletion of individual Phe side chains. The Phe27Ala mutation affected the rate of unfolding only; the Phe15Ala and Phe17Ala mutations changed the kinetics of both unfolding and refolding. Although the activated state of folding of the cold shock protein is overall native-like, individual side chains are still in a non-native environment.     

Striking stabilization of Arc repressor by an engineered disulfide bond

A solvent-exposed Cys11-Cys11' disulfide bond was designed to link the antiparallel strands of the beta sheet both in the Arc repressor dimer and in a single-chain variant in which the Arc subunits are connected by a 15-residue peptide tether. In both proteins, the presence of the disulfide bond increased the T(m) by approximately 40 degrees C. In the single-chain background, the disulfide bond stabilized Arc by 8.5 kcal/mol relative to the reduced form, a significantly larger degree of stabilization than caused by other engineered disulfides and most natural disulfides. This exceptional stabilization arises from a modest effective concentration of the Cys11-Cys11' disulfide in the native state (71 M) and an anomalously low effective concentration in the denatured state (40 microM). Disulfide cross-linking of the two beta strands in the single-chain Arc background accelerated refolding by a factor of 170 into the sub-microsecond time scale. However, the major energetic effect of the disulfide occurs after the transition state for Arc refolding, slowing unfolding by 200 000-fold.     

Thermal unfolding simulations of a multimeric protein--transition state and unfolding pathways

The folding of an oligomeric protein poses an extra challenge to the folding problem because the protein not only has to fold correctly; it has to avoid nonproductive aggregation. We have carried out over 100 molecular dynamics simulations using an implicit solvation model at different temperatures to study the unfolding of one of the smallest known tetramers, p53 tetramerization domain (p53tet). We found that unfolding started with disruption of the native tetrameric hydrophobic core. The transition state for the tetramer to dimer transition was characterized as a diverse ensemble of different structures using Phi value analysis in quantitative agreement with experimental data. Despite the diversity, the ensemble was still native-like with common features such as partially exposed tetramer hydrophobic core and shifts in the dimer-dimer arrangements. After passing the transition state, the secondary and tertiary structures continued to unfold until the primary dimers broke free. The free dimer had little secondary structure left and the final free monomers were random-coil like. Both the transition states and the unfolding pathways from these trajectories were very diverse, in agreement with the new view of protein folding. The multiple simulations showed that the folding of p53tet is a mixture of the framework and nucleation-condensation mechanisms and the folding is coupled to the complex formation. We have also calculated the entropy and effective energy for the different states along the unfolding pathway and found that the tetramerization is stabilized by hydrophobic interactions.     

Conformational stability of dimeric and monomeric forms of the C-terminal domain of human immunodeficiency virus-1 capsid protein

The unfolding equilibrium of the C-terminal domain of human immunodeficiency virus-1 (HIV-1) capsid protein has been analyzed by circular dichroism and fluorescence spectroscopy. The results for the dimeric, natural domain are consistent with a three-state model (N(2)<-->2I<-->2U). The dimer (N(2)) dissociates and partially unfolds in a coupled cooperative process, into a monomeric intermediate (I) of very low conformational stability. This intermediate, which is the only significantly populated form at low (1 microM) protein concentrations, fully preserves the secondary structure but has lost part of the tertiary (intramonomer) interactions found in the dimer. In a second transition, the intermediate cooperatively unfolds into denatured monomer (U). The latter process is the equivalent of a two-state unfolding transition observed for a monomeric domain in which Trp184 at the dimer interface had been truncated to Ala. A highly conserved, disulfide-bonded cysteine, but not the disulfide bond itself, and three conserved residues within the major homology region of the retroviral capsid are important for the conformational stability of the monomer. All these residues are involved also in the association process, despite being located far away from the dimerization interface. It is proposed that dimerization of the C-terminal domain of the HIV-1 capsid protein involves induced-fit recognition, and the conformational reorganization also improves substantially the low intrinsic stability of each monomeric half.     

An on-pathway hidden intermediate and the early rate-limiting transition state of Rd-apocytochrome b562 characterized by protein engineering

The folding pathway of Rd-apocytochrome b562, a four-helix bundle protein, was characterized using Trp and Ala/Gly pair mutations. We found that the Trp mutants (F65W) of both the fully folded Rd-apocytochrome b562 and a partially unfolded intermediate with the N-terminal helix (helix I) unfolded, fold with identical folding rates, providing direct evidence for the conclusion that the rate-limiting transition state folds before the partially unfolded intermediate; and that this hidden intermediate is an on-pathway intermediate. We further characterized the helical structures formed in the rate-limiting transition state by measuring the folding/unfolding rates for Ala/Gly pair mutations at solvent-exposed positions. Little change in folding rates occurred for the Ala/Gly pair mutations at positions in helix I and the C-terminal regions of helix II and IV. In contrast, a significant difference in folding rates was observed for the Ala/Gly pair mutations in helix III and the N-terminal regions of helix II and IV, suggesting that helix III and the N-terminal regions of helix II and IV are formed in the rate-limiting transition state. These results complement those obtained from earlier studies and help to define the folding pathway of Rd-apocytochrome b562 in more detail.     

The Folding Free-Energy Surface of HIV-1 Protease: Insights into the Thermodynamic Basis for Resistance to Inhibitors

Spontaneous mutations at numerous sites distant from the active site of human immunodeficiency virus type 1 protease enable resistance to inhibitors while retaining enzymatic activity. As a benchmark for probing the effects of these mutations on the conformational adaptability of this dimeric β-barrel protein, the folding free-energy surface of a pseudo-wild-type variant, HIV-PR^?, was determined by a combination of equilibrium and kinetic experiments on the urea-induced unfolding/refolding reactions. The equilibrium unfolding reaction was well described by a two-state model involving only the native dimeric form and the unfolded monomer. The global analysis of the kinetic folding mechanism reveals the presence of a fully folded monomeric intermediate that associates to form the native dimeric structure. Independent analysis of a stable monomeric version of the protease demonstrated that a small-amplitude fluorescence phase in refolding and unfolding, not included in the global analysis of the dimeric protein, reflects the presence of a transient intermediate in the monomer folding reaction. The partially folded and fully folded monomers are only marginally stable with respect to the unfolded state, and the dimerization reaction provides a modest driving force at micromolar concentrations of protein. The thermodynamic properties of this system are such that mutations can readily shift the equilibrium from the dimeric native state towards weakly folded states that have a lower affinity for inhibitors but that could be induced to bind to their target proteolytic sites. Presumably, subsequent secondary mutations increase the stability of the native dimeric state in these variants and, thereby, optimize the catalytic properties of the resistant human immunodeficiency virus type 1 protease.     

Intermonomer Hydrogen Bonds Enhance GxxxG-Driven Dimerization of the BNIP3 Transmembrane Domain: Roles for Sequence Context in Helix-Helix Association in Membranes

We determined the sequence dependence of human BNIP3 transmembrane domain dimerization using the biological assay TOXCAT. Mutants in which intermonomer hydrogen bonds between Ser172 and His173 are abolished show moderate interaction, indicating that side-chain hydrogen bonds contribute to dimer stability but are not essential to dimerization. Mutants in which a GxxxG motif composed of Gly180 and Gly184 has been abolished show little or no interaction, demonstrating the critical nature of the GxxxG motif to BNIP3 dimerization. These findings show that side-chain hydrogen bonds can enhance the intrinsic dimerization of a GxxxG motif and that sequence context can control how hydrogen bonds influence helix-helix interactions in membranes. The dimer interface mapped by TOXCAT mutagenesis agrees closely with the interfaces observed in the NMR structure and inferred from mutational analysis of dimerization on SDS-PAGE, showing that the native dimer structure is retained in detergents. We show that TOXCAT and SDS-PAGE give complementary and consistent information about BNIP3 transmembrane domain dimerization: TOXCAT is insensitive to mutations that have modest effects on self-association in detergents but readily discriminates among mutations that completely disrupt detergent-resistant dimerization. The close agreement between conclusions reached from TOXCAT and SDS-PAGE data for BNIP3 suggests that accurate estimates of the relative effects of mutations on native-state protein-protein interactions can be obtained even when the detergent environment is strongly disruptive.     

Roles of dimerization in folding and stability of ketosteroid isomerase from Pseudomonas putida biotype B

Equilibrium and kinetic analyses have been performed to elucidate the roles of dimerization in folding and stability of KSI from Pseudomonas putida biotype B. Folding was reversible in secondary and tertiary structures as well as in activity. Equilibrium unfolding transition, as monitored by fluorescence and ellipticity measurements, could be modeled by a two-state mechanism without thermodynamically stable intermediates. Consistent with the two-state model, one dimensional (1D) NMR spectra and gel-filtration chromatography analysis did not show any evidence for a folded monomeric intermediate. Interestingly enough, Cys 81 located at the dimeric interface was modified by DTNB before unfolding. This inconsistent result might be explained by increased dynamic motion of the interface residues in the presence of urea to expose Cys 81 more frequently without the dimer dissociation. The refolding process, as monitored by fluorescence change, could best be described by five kinetic phases, in which the second phase was a bimolecular step. Because <30% of the total fluorescence change occurred during the first step, most of the native tertiary structure may be driven to form by the bimolecular step. During the refolding process, negative ellipticity at 225 nm increased very fast within 80 msec to account for >80% of the total amplitude. This result suggests that the protein folds into a monomer containing most of the alpha-helical structures before dimerization. Monitoring the enzyme activity during the refolding process could estimate the activity of the monomer that is not fully active. Together, these results stress the importance of dimerization in the formation and maintenance of the functional native tertiary structure.     

Identification of equilibrium and kinetic intermediates involved in folding of urea-denatured creatine kinase

The unfolding transition and kinetic refolding of dimeric creatine kinase after urea denaturation were monitored by intrinsic fluorescence and far ultraviolet circular dichroism. An equilibrium intermediate and a kinetic folding intermediate were identified and characterized. The fluorescence intensity of the equilibrium intermediate is close to that of the unfolded state, whereas its ellipticity at 222 nm is about 50% of the native state. The transition curves measured by these two methods are therefore non-coincident. The kinetic folding intermediate, formed during the burst phase of refolding under native-like conditions, possesses 75% of the native secondary structure, but is mostly lacking in native tertiary structure. In moderate concentrations of urea, only the initial, rapid change in fluorescence intensity or negative ellipticity is observed, and the final state values do not reach the equivalent unfolding values. The unfolding and refolding transition curves measured under identical conditions are non-coincident within the transition from intermediate to fully unfolded state. It is observed by SDS-PAGE that disulfide bond-linked dimeric or oligomeric intermediates are formed in moderate urea concentrations, especially in the refolding reaction. These rapidly formed, soluble intermediates represent an off-pathway event that leads to the hysteresis in the refolding transition curves.     

Equilibrium unfolding of dimeric and engineered monomeric forms of lambda Cro (F58W) repressor and the effect of added salts: evidence for the formation of folded monomer induced by sodium perchlorate

The equilibrium unfolding transitions of Cro repressor variants, dimeric variant Cro F58W and monomer Cro K56[DGEVK]F58W, have been studied by urea and guanidine hydrochloride to probe the folding mechanism. The unfolding transitions of a dimeric variant are well described by a two state process involving native dimer and unfolded monomer with a free energy of unfolding, DeltaG(0,un)(0), of approximately 10-11 kcal/mol. The midpoint of transition curves is dependent on total protein concentration and DeltaG(0,un)(0) is independent of protein concentration, as expected for this model. Unfolding of Cro monomer is well described by the standard two state model. The stability of both forms of protein increases in the presence of salt but decreases with the decrease in pH. Because of the suggested importance of a N2<-->2F dimerization process in DNA binding, we have also studied the effect of sodium perchlorate, containing the chaotropic perchlorate anion, on the conformational transition of Cro dimer by CD, fluorescence and NMR (in addition to urea and guanidine hydrochloride) in an attempt both to characterize the thermodynamics of the process and to identify conditions that lead to an increase in the population of the folded monomers. Data suggest that sodium perchlorate stabilizes the protein at low concentration (<1.5 M) and destabilizes the protein at higher perchlorate concentration with the formation of a "significantly folded" monomer. The tryptophan residue in the "significantly folded" monomer induced by perchlorate is more exposed to the solvent than in native dimer.     

Folding pathway of FKBP12 and characterisation of the transition state.

The folding pathway of human FKBP12, a 12 kDa FK506-binding protein (immunophilin), has been characterised. Unfolding and refolding rate constants have been determined over a wide range of denaturant concentrations and data are shown to fit to a two-state model of folding in which only the denatured and native states are significantly populated, even in the absence of denaturant. This simple model for folding, in which no intermediate states are significantly populated, is further supported from stopped-flow circular dichroism experiments in which no fast "burst" phases are observed. FKBP12, with 107 residues, is the largest protein to date which folds with simple two-state kinetics in water (kF=4 s(-1)at 25 degrees C). The topological crossing of two loops in FKBP12, a structural element suggested to cause kinetic traps during folding, seems to have little effect on the folding pathway.The transition state for folding has been characterised by a series of experiments on wild-type FKBP12. Information on the thermodynamic nature of, the solvent accessibility of, and secondary structure in, the transition state was obtained from experiments measuring the unfolding and refolding rate constants as a function of temperature, denaturant concentration and trifluoroethanol concentration. In addition, unfolding and refolding studies in the presence of ligand provided information on the structure of the ligand-binding pocket in the transition state. The data suggest a compact transition state relative to the unfolded state with some 70 % of the surface area buried. The ligand-binding site, which is formed mainly by two loops, is largely unstructured in the transition state. The trifluoroethanol experiments suggest that the alpha-helix may be formed in the transition state. These results are compared with results from protein engineering studies and molecular dynamics simulations (see the accompanying paper).     

Acceleration of the refolding of Arc repressor by nucleic acids and other polyanions.

The refolding rate of the Arc repressor dimer can be accelerated 30-fold or more by negatively charged polymers including single-stranded and double-stranded DNA, RNA, and polyvinylsulfate but not by neutral or positively charged polymers. The salt-dependence of the polyanion-mediated process and mutant studies indicate that electrostatic interactions are important in the rate acceleration. Urea-dependence studies suggest that Arc is relatively unstructured in the transition state for polyanion-stimulated refolding. At low ionic strength, the observed kinetics of refolding are consistent with a model in which denatured Arc monomers bind rapidly and nonspecifically to the polyanion and complete folding in the bound state.     

Dimerization of the p53 oligomerization domain: identification of a folding nucleus by molecular dynamics simulations

Dimerization of the p53 oligomerization domain involves coupled folding and binding of monomers. To examine the dimerization, we have performed molecular dynamics (MD) simulations of dimer folding from the rate-limiting transition state ensemble (TSE). Among 799 putative transition state structures that were selected from a large ensemble of high-temperature unfolding trajectories, 129 were identified as members of the TSE via calculation of a 50% transmission coefficient from at least 20 room-temperature simulations. This study is the first to examine the refolding of a protein dimer using MD simulations in explicit water, revealing a folding nucleus for dimerization. Our atomistic simulations are consistent with experiment and offer insight that was previously unobtainable.     

The thermal denaturation of ribonuclease A in aqueous-methanol solvents.

Circular dichroism was used to monitor the thermal unfolding of ribonuclease A in 50% aqueous methanol. The spectrum of the protein at temperatures below -10 degrees C (pH* 3.0) was essentially identical to that of native ribonuclease A in aqueous solution. The spectrum of the thermally denatured material above 70 degrees C revealed some residual secondary structure in comparison to protein unfolded by 5 M Gdn.HCl at 70 degrees C in the presence or absence of methanol. The spectra as a function of temperature were deconvoluted to determine the contributions of different types of secondary structure. The position of the thermal unfolding transition as monitored by alpha-helix, with a midpoint at 38 degrees C, was at a much higher temperature than that monitored by beta-sheet, 26 degrees C, which also corresponded to that observed by delta A286, tyrosine fluorescence and hydrodynamic radius (from light scattering measurements). Thus, the loss of beta-sheet structure is decoupled from that of alpha-helix, suggesting a step-wise unfolding of the protein. The transition observed for loss of alpha-helix coincides with the previously measured transition for His-12 by NMR from a partially folded state to the unfolded state, suggesting that the unfolding of the N-terminal helix in RNase A is lost after unfolding of the core beta-sheet during thermal denaturation. The thermally denatured protein was relatively compact, as measured by dynamic light scattering.     

Protein Folding Mechanism of the Dimeric AmphiphysinII/Bin1 N-BAR Domain

The human AmphyphisinII/Bin1 N-BAR domain belongs to the BAR domain superfamily, whose members sense and generate membrane curvatures. The N-BAR domain is a 57 kDa homodimeric protein comprising a six helix bundle. Here we report the protein folding mechanism of this protein as a representative of this protein superfamily. The concentration dependent thermodynamic stability was studied by urea equilibrium transition curves followed by fluorescence and far-UV CD spectroscopy. Kinetic unfolding and refolding experiments, including rapid double and triple mixing techniques, allowed to unravel the complex folding behavior of N-BAR. The equilibrium unfolding transition curve can be described by a two-state process, while the folding kinetics show four refolding phases, an additional burst reaction and two unfolding phases. All fast refolding phases show a rollover in the chevron plot but only one of these phases depends on the protein concentration reporting the dimerization step. Secondary structure formation occurs during the three fast refolding phases. The slowest phase can be assigned to a proline isomerization. All kinetic experiments were also followed by fluorescence anisotropy detection to verify the assignment of the dimerization step to the respective folding phase. Based on these experiments we propose for N-BAR two parallel folding pathways towards the homodimeric native state depending on the proline conformation in the unfolded state.     

Context-dependent effects of proline residues on the stability and folding pathway of ubiquitin.

Substitution of trans-proline at three positions in ubiquitin (residues 19, 37 and 38) produces significant context-dependent effects on protein stability (both stabilizing and destabilizing) that reflect changes to a combination of parameters including backbone flexibility, hydrophobic interactions, solvent accessibility to polar groups and intrinsic backbone conformational preferences. Kinetic analysis of the wild-type yeast protein reveals a predominant fast-folding phase which conforms to an apparent two-state folding model. Temperature-dependent studies of the refolding rate reveal thermodynamic details of the nature of the transition state for folding consistent with hydrophobic collapse providing the overall driving force. Br?nsted analysis of the refolding and unfolding rates of a family of mutants with a variety of side chain substitutions for P37 and P38 reveals that the two prolines, which are located in a surface loop adjacent to the C terminus of the main alpha-helix (residues 24-33), are not significantly structured in the transition state for folding and appear to be consolidated into the native structure only late in the folding process. We draw a similar conclusion regarding position 19 in the loop connecting the N-terminal beta-hairpin to the main alpha-helix. The proline residues of ubiquitin are passive spectators in the folding process, but influence protein stability in a variety of ways.     

Cofactor and tryptophan accessibility and unfolding of brain glutamate decarboxylase

Cofactor and tryptophan accessibility of the 65-kDa form of rat brain glutamate decarboxylase (GAD) was investigated by fluorescence quenching measurements using acrylamide, I-, and Cs+ as the quenchers. Trp residues were partially exposed to solvent. I- was less able and Cs+ was more able to quench the fluorescence of Trp residues in the holoenzyme of GAD (holoGAD) than the apoenzyme (apoGAD). The fraction of exposed Trp residues were in the range of 30-49%. In contrast, pyridoxal-P bound to the active site of GAD was exposed to solvent. I- was more able and Cs+ was less able to quench the fluorescence of pyridoxal-P in holoGAD. The cofactor was present in a positively charged microenvironment, making it accessible for interactions with anions. A difference in the exposure of Trp residues and pyridoxal-P to these charged quenchers suggested that the exposed Trp residues were essentially located outside of the active site. Changes in the accessibility of Trp residues upon pyridoxal-P binding strongly supported a significant conformational change in GAD. Fluorescence intensity measurements were also carried out to investigate the unfolding of GAD using guanidine hydrochloride (GdnHCl) as the denaturant. At 0.8-1.5 M GdnHCl, an intermediate step was observed during the unfolding of GAD from the native to the denatured state, and was not found during the refolding of GAD from the denatured to native state, indicating that this intermediate step was not a reversible process. However, at >1.5 M GdnHCl for holoGAD and >2.0 M GdnHCl for apoGAD, the transition leading to the denatured state was reversible. It was suggested that the intermediate step involved the dissociation of native dimer of GAD into monomers and the change in the secondary structure of the protein. Circular dichroism revealed a decrease in the alpha-helix content of GAD from 36 to 28%. The unfolding pattern suggested that GAD may consist of at least two unfolding domains. Unfolding of the lower GdnHCl-resisting domain occurred at a similar concentration of denaturant for apoGAD and holoGAD, while unfolding of the higher GdnHCl-resisting domain occurred at a higher concentration of GdnHCl for apoGAD than holoGAD.