Early events during folding of wild-type staphylococcal nuclease and a single-tryptophan variant studied by ultrarapid mixing

A continuous-flow mixing device with a dead time of 100 micros coupled with intrinsic tryptophan and 1-anilinonaphthalene-8-sulfonate (ANS) fluorescence was used to monitor structure formation during early stages of the folding of staphylococcal nuclease (SNase). A variant with a unique tryptophan fluorophore in the N-terminal beta-barrel domain (Trp76 SNase) was obtained by replacing the single Trp140 in wild-type SNase with His in combination with Trp substitution of Phe76. A common background of P47G, P117G and H124L mutations was chosen in order to stabilize the protein and prevent accumulation of cis proline isomers under native conditions. In contrast to WT(*) SNase, which shows no changes in tryptophan fluorescence prior to the rate-limiting folding step ( approximately 100 ms), the F76W/W140H variant shows additional changes (enhancement) during an early folding phase with a time constant of 75 micros. Both proteins exhibit a major increase in ANS fluorescence and identical rates for this early folding event. These findings are consistent with the rapid accumulation of an ensemble of states containing a loosely packed hydrophobic core involving primarily the beta-barrel domain while the specific interactions in the alpha-helical domain involving Trp140 are formed only during the final stages of folding. The fact that both variants exhibit the same number of kinetic phases with very similar rates confirms that the folding mechanism is not perturbed by the F76W/W140H mutations. However, the Trp at position 76 reports on the rapid formation of a hydrophobic cluster in the N-terminal beta-sheet region while the wild-type Trp140 is silent during this early stage of folding. Quantitative modeling of the (un)folding kinetics and thermodynamics of these two proteins versus urea concentration revealed that the F76W/W140H mutation selectively destabilizes the native state relative to WT(*) SNase while the stability of transient intermediates remains unchanged, leading to accumulation of intermediates under equilibrium conditions at moderate denaturant concentrations.     

Role of C-terminal region of Staphylococcal nuclease for foldability,stability, and activity

The role of the C-terminal region of Staphylococcal nuclease (SNase) was examined by deletion mutation. Deletions up to eight residues do not affect the structure and function. The structure and enzymatic activity were partially lost by deleting Ser141-Asn149 (Delta141-149), and deletion of Trp140-Asn149 (Delta140-149) resulted in further loss of structure and activity. A 13-residue deletion showed the same effect as the 10-residue deletion. Both Ser141Gln and Ser141Ala mutations for an eight-residue deletion mutant did not alter properties as well as Ser141A1a for full-length SNase. In contrast, Trp140Ala mutation for Delta141-149 shows the same effect as the deletion of Trp140. Trp140Ala mutation for full-length SNase causes the loss of native structure. These observations indicate the significance of the 140th and the 141st residues. The side-chain of the 140th residue is required to be tryptophan; however, the backbone of the 141st residue is solely critical for foldability, but the side-chain information is not crucial. All of the mutants that take a non-native conformation show enzymatic activity and inhibitor-induced folding, suggesting that foldability is required for the activity.     

Folding stability and cooperativity of the three forms of 1-110 residues fragment of staphylococcal nuclease

Folding stability and cooperativity of the three forms of 1-110 residues fragment of staphylococcal nuclease (SNase110) have been studied by various biophysical and NMR methods. Samples of G-88W- and V-66W-mutant SNase110, namely G-88W110 and V-66W110, in aqueous solution and SNase110 in 2.0 M TMAO are adopted in this study. The unfolding transitions and folded conformations of the three SNase fragments were detected by far- and near-ultraviolet circular dichroism and intrinsic tryptophan fluorescence measurements. The tertiary structures and internal motions of the fragments were determined by NMR spectroscopy. Both G-88W and V-66W single mutations as well as a small organic osmolyte (Trimethylamine N-oxide, TMAO) can fold the fragment into a native-like conformation. However, the tertiary structures of the three fragments exhibit different degrees of folding stability and compactness. G-88W110 adopts a relatively rigid structure representing a most stable native-like beta-subdomain conformation of the three fragments. V-66W110- and TMAO-stabilized SNase110 produce less compact structures having a less stable "beta-barrel" structural region. The different folding status accounts for the different backbone dynamic and urea-unfolding transition features of the three fragments. The G-20I/G-29I-mutant variants of the three fragments have provided the evidence that the folding status is correlated closely to the packing of the beta-strands in the beta-barrel of the fragments. The native-like beta-barrel structural region acts as a nonlocal nucleus for folding the fragment. The tertiary folding of the three fragments is initiated by formation of the local nucleation sites at two beta-turn regions, I-18-D-21 and Y-27-Q-30, and developed by the formation of a nonlocal nucleation site at the beta-barrel region. The formation of beta-barrel and overall structure is concerted, but the level of cooperativity is different for the three 1-110 residues SNase fragments.     

The residual pro-part of cathepsin C fulfills the criteria required for an intramolecular chaperone in folding and stabilizing the human proenzyme

The 13.5 kDa N-terminal part of the propeptide remains associated with mature cathepsin C after proteolytic activation and excision of the activation peptide. This residual pro-part, isolated from the recombinant enzyme, folds spontaneously and rapidly to a stable, compact monomer with secondary structure and stable tertiary interactions. Folding and unfolding kinetics of the residual pro-part with intact disulfides are complex, and accumulation of transient intermediates is observed. The cleaved form of the pro-part isolated from natural human cathepsin C also folds, suggesting that the intact form comprises two folding domains. The linkages of the two disulfide bridges have been established as 30-118 and 54-136 for the native enzyme. The native disulfide bonds can be re-formed from the fully reduced and denatured state by oxidative refolding, resulting in a domain that is spectroscopically indistinguishable from the original refolded residual pro-part. Both disulfides are solvent-exposed and can be reduced in the absence of denaturant. The reduced form retains most or all of the native tertiary structure and is only approximately 2 kcal.mol(-1) less stable than the oxidized form. It folds fast relative to the rate of biosynthesis, to the same conformation as the oxidized form. Folding and disulfide formation are sequential. These results indicate that the proenzyme folds sequentially in vivo and that the residual pro-part constitutes a rapidly and independently folding domain that stabilizes the mature enzyme. It thus fulfills the criteria required of an intramolecular chaperone. It may also be involved in stabilizing the tetrameric structure of the mature enzyme.     

Is protein folding hierarchic? II. Folding intermediates and transition states

The folding reactions of some small proteins show clear evidence of a hierarchic process, whereas others, lacking detectable intermediates, do not. Evidence from folding intermediates and transition states suggests that folding begins locally, and that the formation of native secondary structure precedes the formation of tertiary interactions, not the reverse. Some notable examples in the literature have been interpreted to the contrary. For these examples, we have simulated the local structures that form when folding begins by using the LINUS program with nonlocal interactions turned off. Our results support a hierarchic model of protein folding.     

Tight complex formation between Cosmc chaperone and its specific client non-native T-synthase leads to enzyme activity and client-driven dissociation

The interaction of the endoplasmic reticulum molecular chaperone Cosmc with its specific client T-synthase (Core 1 β1-3-galactosyltransferase) is required for folding of the enzyme and eventual movement of the T-synthase to the Golgi, but the mechanism of interaction is unclear. Here we show that the lumenal domain of recombinant Cosmc directly interacts specifically in either free form or covalently bound to solid supports with denatured T-synthase but not with the active dimeric form of the enzyme. This leads to formation of a relatively stable complex of Cosmc and denatured T-synthase accompanied by formation of reactivated enzyme in an ATP-independent fashion that is not regulated by redox, calcium, pH, or intermolecular disulfide bond formation. The partly refolded and active T-synthase remains tightly bound noncovalently to Cosmc. Dissociation of T-synthase from the complex is promoted by further interactions of the complex with free forms of either native or non-native T-synthase. Taken together, these results demonstrate a novel mechanism in which Cosmc cycles to bind non-native T-synthase, leading to enzyme activity and release in a client-driven process.     

Probing the folding capacity and residual structures in 1-79 residues fragment of staphylococcal nuclease by biophysical and NMR methods

1-79 residues SNase fragment (SNase79) has chain length containing a sequence for helix alpha(1), omega-loop, beta(I)-sheet, and partial beta(II)-sheet of native SNase. The incomplete "beta-barrel" structural region of SNase79 makes this fragment to be interested in investigation of its conformation. For this study, we use CD, fluorescence, and NMR spectroscopy to probe the folding capacity and the residual structures in SNase79. The optical spectra obtained for SNase79 and its mutants reveal the presence of retained capacity for folding of the fragment. The NMR derived (13)C(alpha) secondary chemical shifts, (3)J(NH-Halpha) coupling constants, amide-proton temperature coefficients, interresidue NOEs, and (15)N relaxation data determine the intrinsic propensities for helix- and turn- or beta-sheet-like conformations of SNase79, which is not the result of stabilizing inter-molecular interactions by oligomerization effects. The residual turn- and helix-like structures may serve as potential local nucleation sites, whereas the residual beta(I)-sheet-like structure can be regarded as a potential non-local nucleation site in the folding of SNase79. The intrinsic local and non-local interactions in these potential initiation sites are insufficient to stabilize the folding of SNase79 due to the shortage of relevant long-range interactions from other part of the fragment. The conformational ensemble of SNase79 is a highly heterogeneous collection of interconverting conformations having transiently populated helix- and beta-sheet- or turn-like structures.     

Local interactions and the optimization of protein folding

The role of local interactions in protein folding has recently been the subject of some controversy. Here we investigate an extension of Zwanzig's simple and general model of folding in which local and nonlocal interactions are represented by functions of single and multiple conformational degrees of freedom, respectively. The kinetics and thermodynamics of folding are studied for a series of energy functions in which the energy of the native structure is fixed, but the relative contributions of local and nonlocal interactions to this energy are varied over a broad range. For funnel shaped energy landscapes, we find that 1) the rate of folding increases, but the stability of the folded state decreases, as the contribution of local interactions to the energy of the native structure increases, and 2) the amount of native structure in the unfolded state and the transition state vary considerably with the local interaction strength. Simple exponential kinetics and a well-defined free energy barrier separating folded and unfolded states are observed when nonlocal interactions make an appreciable contribution to the energy of the native structure; in such cases a transition state theory type approximation yields reasonably accurate estimates of the folding rate. Bumps in the folding funnel near the native state, which could result from desolvation effects, side chain freezing, or the breaking of nonnative contacts, significantly alter the dependence of the folding rate on the local interaction strength: the rate of folding decreases when the local interaction strength is increased beyond a certain point. A survey of the distribution of strong contacts in the protein structure database suggests that evolutionary optimization has involved both kinetics and thermodynamics: strong contacts are enriched at both very short and very long sequence separations. Proteins 29:282–291, 1997. © 1997 Wiley-Liss, Inc.     

Compact dimension of denatured states of staphylococcal nuclease

Fluorescence and circular dichroism stopped-flow have been widely used to determine the kinetics of protein folding including folding rates and possible folding pathways. Yet, these measurements are not able to provide spatial information of protein folding/unfolding. Especially, conformations of denatured states cannot be elaborated in detail. In this study, we apply the method of fluorescence energy transfer with a stopped-flow technique to study global structural changes of the staphylococcal nuclease (SNase) mutant K45C, where lysine 45 is replaced by cysteine, during folding and unfolding. By labeling the thiol group of cysteine with TNB (5,5'-dithiobis-2-nitrobenzoic acid) as an energy acceptor and the tryptophan at position 140 as a donor, distance changes between the acceptor and the donor during folding and unfolding are measured from the efficiency of energy transfer. Results indicate that the denatured states of SNase are highly compact regardless of how the denatured states (pH-induced or GdmCl-induced) are induced. The range of distance changes between two probes is between 25.6 and 25.4 A while it is 20.4 A for the native state. Furthermore, the folding process consists of three kinetic phases while the unfolding process is a single phase. These observations agree with our previous sequential model: N(0) left arrow over right arrow D(1) left arrow over right arrow D(2) left arrow over right arrow D(3) (Chen et al., J Mol Biol 1991;220:771-778). The efficiency of protein folding may be attributed to initiating the folding process from these compact denatured structures.     

Acceleration of disulfide-coupled protein folding using glutathione derivatives

Protein folding occurs simultaneously with disulfide bond formation. In general, the in vitro folding of proteins containing disulfide bond(s) is carried out in the presence of redox reagents, such as glutathione, to permit native disulfide pairing to occur. It is well known that the formation of a disulfide bond and the correct tertiary structure of a target protein are strongly affected by the redox reagent used. However, little is known concerning the role of each amino acid residue of the redox reagent, such as glutathione. Therefore, we prepared glutathione derivatives - glutamyl-cysteinyl-arginine (ECR) and arginyl-cysteinyl-glycine (RCG) - and examined their ability to facilitate protein folding using lysozyme and prouroguanylin as model proteins. When the reduced and oxidized forms of RCG were used, folding recovery was greater than that for a typical glutathione redox system. This was particularly true when high protein concentrations were employed, whereas folding recovery using ECR was similar to that of the glutathione redox system. Kinetic analyses of the oxidative folding of prouroguanylin revealed that the folding velocity (K(RCG) = 3.69 × 10(-3) s(-1)) using reduced RCG/oxidized RCG was approximately threefold higher than that using reduced glutathione/oxidized glutathione. In addition, folding experiments using only the oxidized form of RCG or glutathione indicated that prouroguanylin was converted to the native conformation more efficiently in the case of RCG, compared with glutathione. The findings indicate that a positively charged redox molecule is preferred to accelerate disulfide-exchange reactions and that the RCG system is effective in mediating the formation of native disulfide bonds in proteins.     

Effect of mutations involving charged residues on the stability of staphylococcal nuclease: a continuum electrostatics study

A continuum electrostatics model is used to calculate the relative stabilities of 117 mutants of staphylococcal nuclease (SNase) involving the mutation of a charged residue to an uncharged residue. The calculations are based on the crystallographic structure of the wild-type protein and attempt to take implicitly into account the effect of the mutations in the denatured state by assuming a linear relationship between the free energy changes caused by the mutation in the native and denatured states. A good correlation (linear correlation coefficient of approximately 0.8) is found with published experimental relative stabilities of these mutants. The results suggest that in the case of SNase (i) charged residues contribute to the stability of the native state mainly through electrostatic interactions, and (ii) native-like electrostatic interactions may persist in the denatured state. The continuum electrostatics method is only moderately sensitive to model parameters and leads to quasi-predictive results for the relative mutant stabilities (error of 2-3 kJ mol(-1) or of the order of k(B)T), except for mutants in which a charged residue is mutated to glycine.     

Perturbations of the denatured state ensemble: modeling their effects on protein stability and folding kinetics.

By considering the denatured state of a protein as an ensemble of conformations with varying numbers of sequence-specific interactions, the effects on stability, folding kinetics, and aggregation of perturbing these interactions can be predicted from changes in the molecular partition function. From general considerations, the following conclusions are drawn: (1) A perturbation that enhances a native interaction in denatured state conformations always increases the stability of the native state. (2) A perturbation that promotes a non-native interaction in the denatured state always decreases the stability of the native state. (3) A change in the denatured state ensemble can alter the kinetics of aggregation and folding. (4) The loss (or increase) in stability accompanying two mutations, each of which lowers (or raises) the free energy of the denatured state, will be less than the sum of the effects of the single mutations, except in cases where both mutations affect the same set of partially folded conformations. By modeling the denatured state as the ensemble of all non-native conformations of hydrophobic-polar (HP) chains configured on a square lattice, it can be shown that the stabilization obtained from enhancement of native interactions derives in large measure from the avoidance of non-native interactions in the D state. In addition, the kinetic effects of fixing single native contacts in the denatured state or imposing linear gradients in the HH contact probabilities are found, for some sequences, to significantly enhance the efficiency of folding by a simple hydrophobic zippering algorithm. Again, the dominant mechanism appears to be avoidance of non-native interactions. These results suggest stabilization of native interactions and imposition of gradients in the stability of local structure are two plausible mechanisms involving the denatured state that could play a role in the evolution of protein folding and stability.     

Two peptide fragments G55-I72 and K97-A109 from staphylococcal nuclease exhibit different behaviors in conformational preferences for helix formation

Two synthetic peptides, SNasealpha1 and SNasealpha2, corresponding to residues G55-I72 and K97-A109, respectively, of staphylococcal nuclease (SNase), are adopted for detecting the role of helix alpha1 (E57-A69) and helix alpha2 (M98-Q106) in the initiation of folding of SNase. The helix-forming tendencies of the two SNase peptide fragments are investigated using circular dichroism (CD) and two-dimensional (2D) nuclear magnetic resonance (NMR) methods in water and 40% trifluoroethanol (TFE) solutions. The coil-helix conformational transitions of the two peptides in the TFE-H2O mixture are different from each other. SNasealpha1 adopts a low population of localized helical conformation in water, and shows a gradual transition to helical conformation with increasing concentrations of TFE. SNasealpha2 is essentially unstructured in water, but undergoes a cooperative transition to a predominantly helical conformation at high TFE concentrations. Using the NMR data obtained in the presence of 40% TFE, an ensemble of alpha-helical structures has been calculated for both peptides in the absence of tertiary interactions. Analysis of all the experimental data available indicates that formation of ordered alpha-helical structures in the segments E57-A69 and M98-Q106 of SNase may require nonlocal interactions through transient contact with hydrophobic residues in other parts of the protein to stabilize the helical conformations in the folding. The folding of helix alpha1 is supposed to be effective in initiating protein folding. The formation of helix alpha2 depends strongly on the hydrophobic environment created in the protein folding, and is more important in the stabilization of the tertiary conformation of SNase.     

Effect of salts on the stability and folding of staphylococcal nuclease

The stability and folding kinetics of wild-type and a mutant staphylococcal nuclease (SNase) at neutral pH are significantly perturbed by the presence of moderate to high concentrations of salts. Very substantial increases in stability toward thermal and urea denaturation were observed; for example, 0.4 M sodium sulfate increased the free energy of wild-type SNase by more than 2 kcal/mol. For the NCA SNase mutant, the presence of the salts abolished the cold denaturation observed at neutral pH with this variant, and substantially increased its stability. Significant effects of salts on the kinetics of refolding were also observed. For NCA SNase, the presence of the salts markedly increased the folding rates (up to 5-fold). On the other hand, chloride, in particular, substantially decreased the rate of folding of the wild-type protein. Since the rates of the slow phases due to proline isomerization were increased by salt, these steps must be coupled to conformational processes. Fluorescence energy transfer between the lone tryptophan (Trp140) and an engineered fluorescent acceptor at residue 64 revealed that the addition of a high concentration of KCl led to the formation of a transient folding intermediate not observed at lower salt concentrations, and in which residues 140 and 64 were much closer than in the native state. The salt-induced effects on the kinetics of folding are attributed to the enhanced stability of the transient folding intermediates. It is likely that the combination of the high net charge, due to the high isoelectric point, and the relatively low intrinsic hydrophobicity, leads to staphylococcal nuclease having only marginal stability at neutral pH. The salt-induced effects on the structure, stability, and kinetics of staphylococcal nuclease are attributed to the binding of counterions, namely, anions, resulting in minimization of intramolecular electrostatic repulsion. This leads to increased stability, more structure, and greater compactness, as observed. Consequently, localized electrostatic repulsion is present at neutral pH in SNase, probably contributing to its marginal stability. The results suggest that, in general, marginally stable globular proteins will be significantly stabilized by salts under conditions where they have a substantial net charge.     

Stability and folding kinetics of a ubiquitin mutant with a strong propensity for nonnative beta-hairpin conformation in the unfolded state

A F45W mutant of yeast ubiquitin has been used as a model system to examine the effects of nonnative local interactions on protein folding and stability. Mutating the native TLTGK G-bulged type I turn in the N-terminal beta-hairpin to NPDG stabilizes a nonnative beta-strand alignment in the isolated peptide fragment. However, NMR structural analysis of the native and mutant proteins shows that the NPDG mutant is forced to adopt the native beta-strand alignment and an unfavorable type I NPDG turn. The mutant is significantly less stable (approximately 9 kJ mol(-1)) and folds 30 times slower than the native sequence, demonstrating that local interactions can modulate protein stability and that attainment of a nativelike beta-hairpin conformation in the transition state ensemble is frustrated by the turn mutations. Surprising, alcoholic cosolvents [5-10% (v/v) TFE] are shown to accelerate the folding rate of the NPDG mutant. We conclude, backed-up by NMR data on the peptide fragments, that even though nonnative states in the denatured ensemble are highly populated and their stability further enhanced in the presence of cosolvents, the simultaneous increase in the proportion of nativelike secondary structure (hairpin or helix), in rapid equilibrium with nonnative states, is sufficient to accelerate the folding process. It is evident that modulating local interactions and increasing nonnative secondary structure propensities can change protein stability and folding kinetics. However, nonlocal contacts formed in the global cooperative folding event appear to determine structural specificity.     

Searching for folding initiation sites of staphylococcal nuclease: a study of N-terminal short fragments

The N-terminal short fragments of staphylococcal nuclease (SNase), SNase20, SNase28, and SNase36, corresponding to the sequence regions, Ala1-Gly20, Ala1-Lys28, and Ala1-Leu36, respectively, as well as an 8-residue peptide (Ala17-Ile18-Asp19-Gly20-Asp21-Thr22-Val23-Lys24) have been synthesized. The conformational states of these fragments were investigated using CD and NMR spectroscopy in aqueous solution and in trifluoroethanol (TFE)-H(2)O mixture. SNase20 containing a sequence corresponding to a bent peptide in native SNase shows a transient population of bend-like conformation around Ala12-Thr13-Leu14 in TFE-H(2)O mixture. The sequence region of Ala17-Thr22 of SNase28 displays a localized propensity for turn-like conformation in both aqueous solution and TFE-H(2)O mixture. The conformational ensemble of SNase36 in aqueous solution includes populated turn-like conformations localized in sequence regions Ala17-Thr22 and Tyr27-Gln30. The analysis suggests that these sequence regions, which form the regular secondary structures in native protein, may serve as the folding nucleation sites of SNase fragments of different chain lengths starting from the N-terminal end. Thus, the formation of bend- and turn-like conformations of these sequence regions may be involved in the early folding events of the SNase polypeptide chain in vitro.     

Atomic-level structure characterization of an ultrafast folding mini-protein denatured state

Atomic-level analyses of non-native protein ensembles constitute an important aspect of protein folding studies to reach a more complete understanding of how proteins attain their native form exhibiting biological activity. Previously, formation of hydrophobic clusters in the 6 M urea-denatured state of an ultrafast folding mini-protein known as TC5b from both photo-CIDNP NOE transfer studies and FCS measurements was observed. Here, we elucidate the structural properties of this mini-protein denatured in 6 M urea performing (15)N NMR relaxation studies together with a thorough NOE analysis. Even though our results demonstrate that no elements of secondary structure persist in the denatured state, the heterogeneous distribution of R(2) rate constants together with observing pronounced heteronuclear NOEs along the peptide backbone reveals specific regions of urea-denatured TC5b exhibiting a high degree of structural rigidity more frequently observed for native proteins. The data are complemented with studies on two TC5b point mutants to verify the importance of hydrophobic interactions for fast folding. Our results corroborate earlier findings of a hydrophobic cluster present in urea-denatured TC5b comprising both native and non-native contacts underscoring their importance for ultra rapid folding. The data assist in finding ways of interpreting the effects of pre-existing native and/or non-native interactions on the ultrafast folding of proteins; a fact, which might have to be considered when defining the starting conditions for molecular dynamics simulation studies of protein folding.     

Role of kinetic intermediates in the folding of leech carboxypeptidase inhibitor

The oxidative folding and reductive unfolding pathways of leech carboxypeptidase inhibitor (LCI; four disulfides) have been characterized in this work by structural and kinetic analysis of the acid-trapped folding intermediates. The oxidative folding of reduced and denatured LCI proceeds rapidly through a sequential flow of 1-, 2-, 3-, and 4-disulfide (scrambled) species to reach the native form. Folding intermediates of LCI comprise two predominant 3-disulfide species (designated as III-A and III-B) and a heterogeneous population of scrambled isomers that consecutively accumulate along the folding reaction. Our study reveals that forms III-A and III-B exclusively contain native disulfide bonds and correspond to stable and partially structured species that interconvert, reaching an equilibrium prior to the formation of the scrambled isomers. Given that these intermediates act as kinetic traps during the oxidative folding, their accumulation is prevented when they are destabilized, thus leading to a significant acceleration of the folding kinetics. III-A and III-B forms appear to have both native disulfides bonds and free thiols similarly protected from the solvent; major structural rearrangements through the formation of scrambled isomers are required to render native LCI. The reductive unfolding pathway of LCI undergoes an apparent all-or-none mechanism, although low amounts of intermediates III-A and III-B can be detected, suggesting differences in protection against reduction among the disulfide bonds. The characterization of III-A and III-B forms shows that the former intermediate structurally and functionally resembles native LCI, whereas the III-B form bears more resemblance to scrambled isomers.     

Using motion planning to study protein folding pathways.

We present a framework for studying protein folding pathways and potential landscapes which is based on techniques recently developed in the robotics motion planning community. Our focus in this work is to study the protein folding mechanism assuming we know the native fold. That is, instead of performing fold prediction, we aim to study issues related to the folding process, such as the formation of secondary and tertiary structure, and the dependence of the folding pathway on the initial denatured conformation. Our work uses probabilistic roadmap (PRM) motion planning techniques which have proven successful for problems involving high-dimensional configuration spaces. A strength of these methods is their efficiency in rapidly covering the planning space without becoming trapped in local minima. We have applied our PRM technique to several small proteins (~60 residues) and validated the pathways computed by comparing the secondary structure formation order on our paths to known hydrogen exchange experimental results. An advantage of the PRM framework over other simulation methods is that it enables one to easily and efficiently compute folding pathways from any denatured starting state to the (known) native fold. This aspect makes our approach ideal for studying global properties of the protein's potential landscape, most of which are difficult to simulate and study with other methods. For example, in the proteins we study, the folding pathways starting from different denatured states sometimes share common portions when they are close to the native fold, and moreover, the formation order of the secondary structure appears largely independent of the starting denatured conformation. Another feature of our technique is that the distribution of the sampled conformations is correlated with the formation of secondary structure and, in particular, appears to differentiate situations in which secondary structure clearly forms first and those in which the tertiary structure is obtained more directly. Overall, our results applying PRM techniques are very encouraging and indicate the promise of our approach for studying proteins for which experimental results are not available.     

Local stability identification and the role of a key aromatic amino acid residue in staphylococcal nuclease refolding

Staphylococcal nuclease (SNase) is a model protein that contains one domain and no disulfide bonds. Its stability in the native state may be maintained mainly by key amino acids. In this study, two point-mutated proteins each with a single base substitution [alanine for tryptophan (W140A) and alanine for lysine (K133A)] and two truncated fragment proteins (positions 1-139 [SNase(1-139) or W140O] and positions 1-141 [SNase(1-141) or E142O]) were generated. Differential scanning microcalorimetry in thermal denaturation experiments showed that K133A and E142O have nearly unchanged DeltaH(cal) relative to the wild-type, whereas W140A and W140O display zero enthalpy change (DeltaH(cal) approximately 0). Far-UV CD measurements indicate secondary structure in W140A but not W140O, and near-UV CD measurements indicate no tertiary structure in either W140 mutant. These observations indicate an unusually large contribution of W140 to the stability and structural integrity of SNase.