Conformational plasticity in folding of the split beta-alpha-beta protein S6: evidence for burst-phase disruption of the native state

An increasing number of folding studies of two-state proteins shows that point mutations sometimes change the kinetic m-values, leading to kinks and curves in the chevron plots. The molecular origin of these changes is yet unclear although it is speculated that they are linked to structural rearrangement of the transition state or to accumulation of meta-stable intermediates. To shed more light on this issue, we present here a combined m and phi-value analysis of the split beta-alpha-beta protein S6. Wild-type S6 displays classical two-state kinetics with v-shaped chevron plot, but a majority of its mutants display distinct m-value changes or curved chevrons. We observe that this kinetic aberration of S6 is linked to mutations that are clustered in distinct regions of the native structure. The most pronounced changes, i.e. decrease in the m-value for the unfolding rate constant, are seen upon truncation of interactions between the N and C termini, whereas mutations in the centre of the hydrophobic core show smaller or even opposed effects. As a consequence, the calculated phi-values display a systematic increase upon addition of denaturant. In the case of S6, the phenomenon seems to arise from a general plasticity of the different species on the folding pathway. That is, the structure of the denatured ensemble, the transition state, and the native ground-state for unfolding seem to change upon mutation. From these changes, it is concluded that interactions spanning the centre of the hydrophobic core form early in folding, whereas the entropically disfavoured interactions linking the N and C termini consolidate very late, mainly on the down-hill-side of the folding barrier.     

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.     

Simulation and experiment conspire to reveal cryptic intermediates and a slide from the nucleation-condensation to framework mechanism of folding

There is a change from three-state to two-state kinetics of folding across the homeodomain superfamily of proteins as the mechanism slides from framework to nucleation-condensation. The tendency for framework folding in this family correlates with inherent helical propensity. The cellular myeloblastis protein (c-Myb) falls in the mechanistic transition region. An earlier, preliminary report of protein engineering experiments and molecular dynamics simulations (MD) showed that the folding mechanism for this protein has aspects of both the nucleation-condensation and framework models. In the more in-depth analysis of the MD trajectories presented here, we find that folding may be attributed to both of these mechanisms in different regions of the protein. The folding of the loop, middle helix, and turn is best described by nucleation-condensation, whereas folding of the N and C-terminal helices may be described by the framework model. Experimentally, c-Myb folds by apparent two-state kinetics, but the MD simulations predict that the kinetics hide a high-energy intermediate. We stabilized this hypothetical folding intermediate by deleting a residue (P174) in the loop between its second and third helices, and the mutant intermediate is long-lived in the simulations. Equilibrium and kinetic experiments demonstrate that folding of the DeltaP174 mutant is indeed three-state. The presence and shape of the intermediate observed in the simulations were confirmed by small angle X-ray scattering experiments.     

Thermal unfolding molecular dynamics simulation of Escherichia coli dihydrofolate reductase: thermal stability of protein domains and unfolding pathway

Temperature induced unfolding of Escherichia coli dihydrofolate reductase was carried out by using molecular dynamic simulations. The simulations show that the unfolding generally involves an initial end-to-end collapse of the adenine binding domain into partially extended loops, followed by a gradual breakdown of the remaining beta sheet core structure. The core, which consists of beta strands 5-7, was observed to be the most resistant to thermal unfolding. This region, which is made up of part of the N terminus domain and part of the large domain of the E. coli dihydrofolate reductase, may constitute the nucleation site for protein folding and may be important for the eventual formation of both domains. The unfolding of different domains at different stages of the unfolding process suggests that protein domains vary in stability and that the rate at which they unfold can affect the overall outcome of the unfolding pathway. This observation is compared with the recently proposed hierarchical folding model. Finally, the results of the simulation were found to be consistent with a previous experimental study (Frieden, Proc Natl Acad Sci USA 1990;87:4413-4416) which showed that the folding process of E. coli dihydrofolate reductase involves sequential formation of the substrate binding sites.     

Topological determinants of protein unfolding rates

For proteins that fold by two-state kinetics, the folding and unfolding processes are believed to be closely related to their native structures. In particular, folding and unfolding rates are influenced by the native structures of proteins. Thus, we focus on finding important topological quantities from a protein structure that determine its unfolding rate. After constructing graphs from protein native structures, we investigate the relationships between unfolding rates and various topological quantities of the graphs. First, we find that the correlation between the unfolding rate and the contact order is not as prominent as in the case of the folding rate and the contact order. Next, we investigate the correlation between the unfolding rate and the clustering coefficient of the graph of a protein native structure, and observe no correlation between them. Finally, we find that a newly introduced quantity, the impact of edge removal per residue, has a good overall correlation with protein unfolding rates. The impact of edge removal is defined as the ratio of the change of the average path length to the edge removal probability. From these facts, we conclude that the protein unfolding process is closely related to the protein native structure.     

A stable intermediate in the equilibrium unfolding of Escherichia coli citrate synthase.

Urea-induced unfolding of Escherichia coli citrate synthase occurs in two phases, as monitored by circular dichroism at 222 nm (measuring secondary structure) or by tryptophan fluorescence. In this paper we characterize the intermediate state, which retains about 40% of the ellipticity of the native state, and is stable between 2.5 M and 5.5 M urea, approximately. This intermediate binds significant amounts of the probe for hydrophobic surfaces, anilinonaphthalene sulfonate, but forms aggregates at least as high as an octamer, as shown by transverse urea gradient polyacrylamide electrophoresis. Thermal denaturation of E. coli citrate synthase also produces an intermediate at temperatures near 60 degrees C, which also retains about 40% of the native ellipticity and forms aggregates, as measured by electrospray-ionization/time-of-flight mass spectrometry. We have used a collection of "cavity-forming" mutant proteins, in which bulky buried hydrophobic residues are replaced by alanines, to explore the nature of the intermediate state further. A certain amount of these mutant proteins shows a destabilized intermediate, as measured by the urea concentration range in which the intermediate is observed. These mutants are found in parts of the citrate synthase sequence that, in a native state, form helices G, M, N, Q, R, and S. From this and other evidence, it is argued that the intermediate state is an aggregated state in which these six helices, or parts of them, remain folded, and that formation of this intermediate is also likely to be a key step in the folding of E. coli citrate synthase.     

Amyloidogenic unfolding intermediates differentiate sheep prion protein variants

Sheep is a unique example among mammalian species to present a strong correlation between genotype and prion disease susceptibility phenotype. Indeed a well-defined set of PrP polymorphisms at positions 136, 154 and 171 (sheep numbering) govern scrapie susceptibility, ranging from very high susceptibility for V136-R154-Q171 variant (VRQ) to resistance for A136-R154-R171 variant (ARR).To get better insight into the molecular mechanisms of scrapie susceptibility/resistance, the unfolding pathways of the different full-length recombinant sheep prion protein variants were analysed by differential scanning calorimetry in a wide range of pH. In the pH range 4.5-6.0, thermal unfolding occurs through a reversible one-step process while at pH <4.5 and >6.0 unfolding intermediates are formed, which are stable in the temperature range 65-80 degrees C. While these general behaviours are shared by all variants, VRQ and ARQ (susceptibility variants) show higher thermal stability than AHQ and ARR (resistance variants) and the formation of their unfolding intermediates requires higher activation energy than in the case of AHQ and ARR. Furthermore, secondary structures of the unfolding intermediates differentiate variants: ARR unfolding intermediate exhibits random coil structure, contrasting with the beta-sheet structure of VRQ and ARQ unfolding intermediates. The rate of the unfolding intermediate formation allows us to classify genetic variants along increasing scrapie susceptibility at pH 4.0, VRQ and ARQ rates being the highest. Rather poor correlation is observed at pH 7.2. Upon cooling, these intermediates refold into stable species, which are rich in beta-type secondary structures and, as revealed by thioflavin T fluorescence and electron microscopy, share amyloid characteristics. These results highlight the prion protein plasticity genetically modulated in sheep, and might provide a molecular basis for sheep predisposition to scrapie taking into account both thermodynamic stability and transconformation rate of prion protein.     

A comparative study of the unfolding of the endoglucanase Cel45 from Humicola insolens in denaturant and surfactant.

Cellulases are increasingly being used for industrial purposes, particularly in washing powders, yet little is known of the factors governing the stability of proteins in detergent solutions. We present a comparative analysis of the behavior of the cellulase Cel45 from Humicola insolens in the presence of the denaturant guanidinium chloride and the anionic detergent C12-LAS. Although Cel45 unfolds in GdmCl according to a simple two-state model under equilibrium conditions, it accumulates a transient intermediate during refolding. The four disulfide bonds do not contribute detectably to the stability of the native state. Cel45 is unfolded by very low concentrations of C12-LAS (1-4 mM). An analysis of 16 mutants of Cel45 shows a very weak correlation between unfolding rates in denaturant and detergent; mutants that have the same unfolding rate in GdmCl (within a factor of 1.5) vary 1,000-fold in their unfolding rates in C12-LAS. The data support a simple model for unfolding by detergent, in which the introduction of positive charges or removal of negative charges greatly increases detergent sensitivity, while interactions with the hydrophobic detergent tail contribute to a smaller extent. This implies that different detergent-mediated unfolding pathways exist, whose accessibilities depend on individual residues. Double-mutant cycles reveal that mutations in two proximal residues lead to repulsion and a destabilization greater than the sum of the individual mutations as measured by GdmCl denaturation, but they also reduce the affinity for LAS and therefore actually stabilize the protein relative to wild-type. Ligands that interact strongly with the denatured state may therefore alter the unfolding process.     

Structural basis for unfolding pathway‐dependent stability of proteins: Vectorial unfolding versus global unfolding

Point mutations in proteins can have different effects on protein stability depending on the mechanism of unfolding. In the most interesting case of I27, the Ig‐like module of the muscle protein titin, one point mutation (Y9P) yields opposite effects on protein stability during denaturant‐induced “global unfolding” versus “vectorial unfolding” by mechanical pulling force or cellular unfolding systems. Here, we assessed the reason for the different effects of the Y9P mutation of I27 on the overall molecular stability and N‐terminal unraveling by NMR. We found that the Y9P mutation causes a conformational change that is transmitted through β‐sheet structures to reach the central hydrophobic core in the interior and alters its accessibility to bulk solvent, which leads to destabilization of the hydrophobic core. On the other hand, the Y9P mutation causes a bend in the backbone structure, which leads to the formation of a more stable N‐terminal structure probably through enhanced hydrophobic interactions.     

Kinetic folding and unfolding of staphylococcal nuclease and its six mutants studied by stopped-flow circular dichroism

Kinetics of refolding and unfolding of staphylococcal nuclease and its six mutants, each carrying single or double amino acid substitutions, are studied by stopped-flow circular dichroism measurements. A transient kinetic intermediate formed within 10 ms after refolding starts possesses a substantial part of the N-domain core β-structure, whereas helices are formed at the later stages. The structure of the kinetic intermediate is less organized than the structure that is known to be formed by a nuclease 1-136 fragment. Only the refolding kinetics are affected by the mutations in all the mutants except two in which the mutations have changed the native structure. From this result and also from the locations of the mutation sites, the major N-terminal domain of the nuclease in the transition state of folding has a structure nearly identical to the native one. On the other hand, the minor C-terminal domain has previously been shown to be still disorganized in the transition state. The effects of the amino acid substitutions on the stability of the native and the transition states are in good agreement with the changes in the hydration free energy, expected for the corresponding amino acid replacements in the unfolded polypeptide. Since side chains of all the mutated residues are not accessible to solvent in the native structure, the result suggests that it is the unfolded state that is mainly affected by the mutations. © 1995 Wiley-Liss, Inc.     

Structure of stable protein folding intermediates by equilibrium phi-analysis: the apoflavodoxin thermal intermediate

Protein intermediates in equilibrium with native states may play important roles in protein dynamics but, in cases, can initiate harmful aggregation events. Investigating equilibrium protein intermediates is thus important for understanding protein behaviour (useful or pernicious) but it is hampered by difficulties in gathering structural information. We show here that the phi-analysis techniques developed to investigate transition states of protein folding can be extended to determine low-resolution three-dimensional structures of protein equilibrium intermediates. The analysis proposed is based solely on equilibrium data and is illustrated by determination of the structure of the apoflavodoxin thermal unfolding intermediate. In this conformation, a large part of the protein remains close to natively folded, but a 40 residue region is clearly unfolded. This structure is fully consistent with the NMR data gathered on an apoflavodoxin mutant designed specifically to stabilise the intermediate. The structure shows that the folded region of the intermediate is much larger than the proton slow-exchange core at 25 degrees C. It also reveals that the unfolded region is made of elements whose packing surface is more polar than average. In addition, it constitutes a useful guide to rationally stabilise the native state relative to the intermediate state, a far from trivial task.     

Increasing temperature accelerates protein unfolding without changing the pathway of unfolding

We have traditionally relied on extremely elevated temperatures (498K, 225 degrees C) to investigate the unfolding process of proteins within the timescale available to molecular dynamics simulations with explicit solvent. However, recent advances in computer hardware have allowed us to extend our thermal denaturation studies to much lower temperatures. Here we describe the results of simulations of chymotrypsin inhibitor 2 at seven temperatures, ranging from 298K to 498K. The simulation lengths vary from 94ns to 20ns, for a total simulation time of 344ns, or 0.34 micros. At 298K, the protein is very stable over the full 50ns simulation. At 348K, corresponding to the experimentally observed melting temperature of CI2, the protein unfolds over the first 25ns, explores partially unfolded conformations for 20ns, and then refolds over the last 35ns. Above its melting temperature, complete thermal denaturation occurs in an activated process. Early unfolding is characterized by sliding or breathing motions in the protein core, leading to an unfolding transition state with a weakened core and some loss of secondary structure. After the unfolding transition, the core contacts are rapidly lost as the protein passes on to the fully denatured ensemble. While the overall character and order of events in the unfolding process are well conserved across temperatures, there are substantial differences in the timescales over which these events take place. We conclude that 498K simulations are suitable for elucidating the details of protein unfolding at a minimum of computational expense.     

Autophosphorylation of brain microtubule protein: evidence for endogenous protein kinase/phosphoprotein phosphatase cycling and multiple phosphorylation of a microtubule associated protein

Brain microtubule protein, prepared by two types of recycling methods, undergoes “flash” phosphorylation in the presence of [γ-³²p]ATP through sequential action of protein kinase and phosphoprotein phosphatase present in microtubule protein. SDS electrophoretic analysis indicates that MAP₁, tau protein, and tubulin are poorly phosphorylated, and MAP₂ is the major site of phosphorylation. To improve [³²P]phosphoprotein stability in the presence of the kinase/phosphatase cycle, 3′,5′-cyclicAMP, orthophosphate, or fluoride ion may be added. After separation from tubulin by phosphocellulose chromatography, the MAP fraction exhibits autophosphorylation. Finally, the maximal extent of autophosphorylation is observed with an ATP regenerating system using ADP, [³²P]acetyl-P, and bacterial acetate kinase; this results in the incorporation of 3–4 phosphoryl groups per MAP₂ subunit.     

ADP-ribosylation of nuclear protein A24

Nuclear protein A24, which is composed of histone H2A and ubiquitin, a nonhistone protein, joined by an isopeptide linkage [Goldknopf and Busch (1977) $\text{Proc}$. $\text{Natl}$. $\text{Acad}$. $\text{Sci}$. $\text{USA}$$\text{74}$, 864–868], is found to be ADP-ribosylated in isolated rat liver nuclei.     

Inactivation and unfolding of aminoacylase during denaturation in sodium dodecyl sulfate solutions

During denaturation by sodium dodecyl sulfate (SDS), aminoacylase shows a rapid decrease in activity with increasing concentration of the detergent to reach complete inactivation at 1.0 mM SDS. The denatured minus native-enzyme difference spectrum showed two negative peaks at 287 and 295 nm. With the increase of concentration of SDS, both negative peaks increased in magnitude to reach maximal values at 5.0 mM SDS. The fluorescence emission intensity of the enzyme decreased, whereas there was no red shift of emission maximum in SDS solutions of increasing concentration. In the SDS concentration regions employed in the present study, no marked changes of secondary structure of the enzyme have been observed by following the changes in far-ultraviolet CD spectra. The inactivation of this enzyme has been followed and compared with the unfolding observed during denaturation in SDS solutions. A marked inactivation is already evident at low SDS concentration before significant conformational changes can be detected by ultraviolet absorbance and fluorescence changes. The inactivation rate constants of free enzyme and substrate-enzyme complex were determined by the kinetics method of the substrate reaction in the presence of inactivator previously described by Tsou [Tsou (1988),Adv. Enzymol. Related Areas Mol. Biol. 61, 381–436]. It was found that substrate protects against inactivation and at the same SDS concentrations, the inactivation rate of the free enzyme is much higher than the unfolding rate. The above results show that the active sites of metal enzyme containing Zn²⁺ are also situated in a limited and flexible region of the enzyme molecule that is more fragile to denaturants than the protein as a whole.     

The synthesis of freezing-point-depressing protein of the winter flounder Pseudopleuronectusamericanus in Xenpuslaevis oocytes

The serum of the winter flounder $\text{Pseudopleuronectus}$$\text{americanus}$ contains a freezing-point-depressing protein of a molecular weight approximately 10,000 with 60% alanine in its composition. When injected into Xenopus o?cyte, a 6–10 S, poly A-rich RNA preparation isolated from the fish liver polysomes stimulated 3–4 fold the incorporation of [³H] alanine into 10% trichloroacetic acid-soluble, non-dialysable proteins. Analysis of the protein fractions showed a translation product similar in molecular weight and electrophoretic mobility to flounder freezing-point-depressing protein. These observations indicated that the 6–10 S RNA from the flounder contained mRNA for the synthesis of flounder's freezing-point-depressing protein.     

Thermal versus mechanical unfolding of ubiquitin

The authors studied the temperature-induced unfolding of ubiquitin by all-atom Monte Carlo simulations. The unfolding behavior is compared with that seen in previous simulations of the mechanical unfolding of this protein, based on the same model. In mechanical unfolding, secondary-structure elements were found to break in a quite well-defined order. In thermal unfolding, the authors saw somewhat larger event-to-event fluctuations, but the unfolding pathway was still far from random. Two long-lived secondary-structure elements could be identified in the simulations. These two elements have been found experimentally to be the thermally most stable ones. Interestingly, one of these long-lived elements, the first beta-hairpin, was found to break early in the mechanical unfolding simulations. Their combined simulation results thus enable the authors to predict in detail important differences between the thermal and mechanical unfolding behaviors of ubiquitin.     

Investigating protein unfolding kinetics by pulse proteolysis

Investigation of protein unfolding kinetics of proteins in crude samples may provide many exciting opportunities to study protein energetics under unconventional conditions. As an effort to develop a method with this capability, we employed “pulse proteolysis” to investigate protein unfolding kinetics. Pulse proteolysis has been shown to be an effective and facile method to determine global stability of proteins by exploiting the difference in proteolytic susceptibilities between folded and unfolded proteins. Electrophoretic separation after proteolysis allows monitoring protein unfolding without protein purification. We employed pulse proteolysis to determine unfolding kinetics of E. coli maltose binding protein (MBP) and E. coli ribonuclease H (RNase H). The unfolding kinetic constants determined by pulse proteolysis are in good agreement with those determined by circular dichroism. We then determined an unfolding kinetic constant of overexpressed MBP in a cell lysate. An accurate unfolding kinetic constant was successfully determined with the unpurified MBP. Also, we investigated the effect of ligand binding on unfolding kinetics of MBP using pulse proteolysis. On the basis of a kinetic model for unfolding of MBP•maltose complex, we have determined the dissociation equilibrium constant (K_d) of the complex from unfolding kinetic constants, which is also in good agreement with known K_d values of the complex. These results clearly demonstrate the feasibility and the accuracy of pulse proteolysis as a quantitative probe to investigate protein unfolding kinetics.