pH dependent unfolding characteristics of DLC8 dimer: Residue level details from NMR

Environment dependence of folding and unfolding of a protein is central to its function. In the same vein, knowledge of pH dependence of stability and folding/unfolding is crucial for many biophysical equilibrium and kinetic studies designed to understand protein folding mechanisms. In the present study we investigated the guanidine induced unfolding transition of dynein light chain protein (DLC8), a cargo adaptor of the dynein complex in the pH range 7-10. It is observed that while the protein remains a dimer in the entire pH range, its stability is somewhat reduced at alkaline pH. Global unfolding features monitored using fluorescence spectroscopy revealed that the unfolding transition of DLC8 at pH 7 is best described by a three-state model, whereas, that at pH 10 is best described by a two-state model. Chemical shift perturbations due to pH change provided insights into the corresponding residue level structural perturbations in the DLC8 dimer. Likewise, backbone (15)N relaxation measurements threw light on the corresponding motional changes in the dimeric protein. These observations have been rationalized on the basis of expected changes with increasing pH in the protonation states of the titratable residues on the structure of the protein. These, in turn provide an explanation for the change from three-state to two-state guanidine induced unfolding transition as the pH is increased from 7 to 10. All these results exemplify and highlight the role of environment vis-à-vis the sequence and structure of a given protein in dictating its folding/unfolding characteristics.     

Multifrequency calorimetry of the folding/unfolding transition of cytochrome c.

The folding-unfolding transition of Fe(III) cytochrome c has been studied with the new technique of multifrequency calorimetry. Multifrequency calorimetry is aimed at measuring directly the dynamics of the energetic events that take place during a thermally induced transition by measuring the frequency dispersion of the heat capacity. This is done by modulating the folding/unfolding equilibrium using a variable frequency, small oscillatory temperature perturbation (approximately 0.05-0.1 degrees C) centered at the equilibrium temperature of the system. Fe(III) cytochrome c at pH 4 undergoes a fully reversible folding/unfolding transition centered at 67.7 degrees C and characterized by an enthalpy change of 81 kcal/mol and heat capacity difference between unfolded and folded states of 0.9 kcal/K*mol. By measuring the temperature dependence of the frequency dispersion of the heat capacity in the frequency range of 0.1-1 Hz it has been possible to examine the time regime of the enthalpic events associated with the transition. The multifrequency calorimetry results indicate that approximately 85% of the excess heat capacity associated with the folding/unfolding transition relaxes with a single relaxation time of 326 +/- 68 ms at the midpoint of the transition region. This is the first time that the time regime in which heat is absorbed and released during protein folding/unfolding has been measured.     

pH and cation effects on the properties of parallel pyrimidine motif DNA triplexes

The effects of cytosine protonation and various cations on the properties of parallel pyrimidine motif DNA triplexes were intensively investigated and characterized by several different techniques, such as circular dichroism (CD) conformation, ultraviolet (UV) melting, differential scanning calorimetry (DSC) thermal denaturation, and surface plasmon resonance (SPR) real-time dynamics. The comparative CD spectra of the triplex and the corresponding homoduplexes showed that the negative peak at approximately 218 nm would be the eigenpeak of the Hoogsteen paired strand, and moreover, the formation pathway of a triplex was significantly pH-dependent and fell into three groups: under acidic conditions, the triplex is formed by a one-step docking, under near physiological conditions, the Watson-Crick duplex is first structured and then accepts the Hoogsteen third strand into its major groove, and under basic conditions, the triplex is not formed. The pH-dependent thermodynamics of the global triplex, the Watson-Crick antiparallel duplex, and the Crick-Hoogsteen parallel duplex were comparatively discussed for the first time. These data revealed that the thermodynamic stabilities of the Watson-Crick-Hoogsteen triplex and the Crick-Hoogsteen duplex would be strongly dependent on cytosine protonation, but a low-pH environment somewhat destabilized the Watson-Crick duplex. The binding energy of triplex formation would be different from the unfolding energy of triplex melting under acidic conditions due to the disparity in the pathway between the formation and unfolding of a triplex. Real-time dynamic measurements showed that the association and dissociation rate constants of a duplex-to-triplex formation are (1.98 +/- 0.24) x 10(3) M(-1) s(-1) and (4.09 +/- 0.96) x 10(-4) s(-1) at 20 degrees C and pH 6.0, respectively. The formation energy of the duplex-to-triplex transition derived from SPR measurements was in agreement with the unfolding energy of the free Hoogsteen paired duplex derived from UV measurements. The calorimetric enthalpies of the triplex-to-duplex-to-single transition were 39.3 and 75.3 kcal/mol under near physiological conditions (pH 7.0), respectively, which were underestimated relative to the van't Hoff enthalpies. In addition, the effects of various cations, ionic strength, mixed-valent cations, and the position of the C(+)xG.C triplets on the thermodynamics of the triplexes were addressed under near physiological conditions. The interaction of metal ions with the triplexes clearly depended on the type and ionic strength of the cations, and the efficiency with which the cations stabilized the global triplex was in the order Mg(2+) > Mn(2+) > Ca(2+) > Ba(2+) > Na(+). These observations would be useful for the design of triplex-forming oligonucleotides for antigene drugs and therapeutic purposes.     

Kinetics of the Triplex-Duplex Transition in DNA

The kinetics of triplex folding/unfolding is investigated by the single-molecule fluorescence resonance energy transfer (FRET) technique. In neutral pH conditions, the average dwell times in both high-FRET (folded) and low-FRET (unfolded) states are comparable, meaning that the triplex is marginally stable. The dwell-time distributions are qualitatively different: while the dwell-time distribution of the high-FRET state should be fit with at least a double-exponential function, the dwell-time distribution of the low-FRET state can be fit with a single-exponential function. We propose a model where the folding can be trapped in metastable states, which is consistent with the FRET data. Our model also accounts for the fact that the relevant timescales of triplex folding/unfolding are macroscopic.     

Noninteracting local-structure model of folding and unfolding transition in globular proteins. II. Application to two-dimensional lattice proteins

The noninteracting local-structure model of the folding and unfolding transition in globular proteins, the formulation of which was given in the preceding paper, is applied to the analysis of the two-dimensional lattice model of proteins. The lattice model of proteins is a theoretical tool designed to study the statistical-mechanical aspect of the folding and unfolding transition. Its dynamics have been studied by a method of Monte Carlo simulation. The noninteracting local-structure model reproduces the equilibrium properties of the lattice model obtained previously by computer simulation remarkably well, when the specificity of the long-range interactions is strong. This observation indicates that the basic assumption of the noninteracting local-structure model is equivalent to the assumption of strong specificity of intramolecular interactions. It is argued that by assuming this strong specificity, we can emphasize the correct main paths of folding and unfolding transition. The way local structures grow and/or merge along the most probable path of folding in the lattice model is discussed by the noninteracting local-structure model.     

Protein dynamics: From the native to the unfolded state and back again

Simulations to study protein unfolding and folding were performed. The unfolding simulations make use of molecular dynamics and treat an atomic model of barnase in aqueous solvent. The cooperative nature of the unfolding transition and the important role of water are described. The folding simulations are based on a bead model of the protein on a cubic lattice. It is shown for the 27-mer model that a large energy gap between the lowest energy (native) state and the excited states is a necessary and sufficient condition for fast folding.     

Spectroscopic and calorimetric investigation on the DNA triplex formed by d(CTCTTCTTTCTTTTCTTTCTTCTC) and d(GAGAAGAAAGA) at acidic pH.

The equimolar mixture of d(CTCTTCTTTCTTTTCTTTCTTCTC) (dY24) and d(GAGAAGAAAGA) (dR11) [designated (dY24).(dR11)], forms at pH = 5 a DNA triplex, which mimicks the H-DNA structure. The DNA triplex was identified by the following criteria: (i) dY24 and dR11 co-migrate in a poly-acrylamide gel, with a mobility and a retardation coefficient comparable to those observed for an 11-triad DNA triplex, previously characterized in our laboratories (1); (ii) the intercalator ethidium bromide shows a poor affinity for (dR11).(dY24) at pH = 5, and a high affinity at pH = 8; (iii) the (dR11).(dY24) mixture is not a substrate for DNase I at pH = 5; (iv) the CD spectrum of (dR11).(dY24), at pH = 5, is consistent with those previously reported for triple-stranded DNA. The (dR11).(dY24) mixture exhibits a thermally induced co-operative transition, which appears to be monophasic, reversible and concentration dependent. This transition is attributed to the disruption of the DNA triplex into single strands. The enthalpy change of the triplex-coil transition was measured by DSC (delta Hcal = 129 +/- 6 kcal/mol) and, assuming a two-state model, by analysis of UV-denaturation curves (average of two methods delta HUV = 137 +/- 13 kcal/mol). Subtracting from delta Hcal of triplex formation the contributions due to the Watson-Crick helix and to the protonation of the C-residues, we found that each pyrimidine binding into the major groove of the duplex, through a Hoogsteen base pair, is accompanied by an average delta H = -5.8 +/- 0.6 kcal/mol. The effect on the stability of the (dR11).(dY24) triplex due to the substitution of a T:A:T triad with a T:T:T one was also investigated.     

Characterization of a slow folding reaction for the alpha subunit of tryptophan synthase

The equilibria and kinetics of urea-induced unfolding and refolding of the alpha subunit of tryptophan synthase of E. coli have been examined for their dependences on viscosity, pH, and temperature in order to investigate the properties of one of the rate-limiting steps, domain association. A viscosity enhancer, 0.58 M sucrose, was found to slow unfolding and accelerate refolding. This apparently anomalous result was shown to be due to the stabilizing effect of sucrose on the folding reaction. After accounting for this stabilization effect by using linear free-energy plots, the unfolding and refolding kinetics were found to have a viscosity dependence. A decrease in pH was found to stabilize the domain association reaction by increasing the refolding rate and decreasing the unfolding rate. This effect was accounted for by protonation of a single residue with a pK value of 8.8 in the native state and 7.1 in the intermediate, in which the two domains are not yet associated. The activation energy of unfolding is 4.8 kcal/mol, close to the diffusion limit. The negative activation entropy of unfolding, -47 cal/deg-mol, which controls this reaction, may result from ordering of solvent about the newly exposed domain interface of the transition state. These results may provide information on the types of noncovalent interactions involved in domain association and improve the ability to interpret the folding of mutants with single amino-acid substitutions at the interface.     

Thermodynamics and kinetics of the hairpin ribozyme from atomistic folding/unfolding simulations

We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active “docked” conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.     

High temperature unfolding simulations of the TRPZ1 peptide

We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 μs with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.     

Fluorescence studies on denaturation and stability of recombinant human interferon-gamma

Unfolding/folding transitions of recombinant human interferon-gamma (hIFNgamma) in urea and guanidine chloride (Gn.HCl) solutions were studied by fluorescence spectroscopy. At pH 7.4 Gn.HCl was a much more efficient denaturant (midpoint of unfolding C* = 1.1 M and deltaG0 = 13.4 kJ/mol) than urea (C* = 2.8 M and deltaG0 = 11.7 kJ/mol). The close deltaG0 values indicate that the contribution of electrostatic interactions to the stability of hIFNgamma is insignificant. Both the pH dependence of the fluorescence intensity and the unfolding experiments in urea at variable pH showed that hIFNgamma remains native in the pH range of 4.8-9.5. Using two quenchers, iodide and acrylamide, and applying the Stern-Volmer equation, a cluster of acidic groups situated in close proximity to the single tryptophan residue was identified. Based on the denaturation experiments at different pH values and on our earlier calculations of the electrostatic interactions in hIFNgamma, we assume that the protonation of Asp63 causes conformational changes having a substantial impact on the stability of hIFNgamma.     

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.     

Pressure denaturation of staphylococcal nuclease studied by neutron small-angle scattering and molecular simulation

We studied the pressure-induced folding/unfolding transition of staphylococcal nuclease (SN) over a pressure range of approximately 1-3 kilobars at 25 degrees C by small-angle neutron scattering and molecular dynamics simulations. We find that applying pressure leads to a twofold increase in the radius of gyration derived from the small-angle neutron scattering spectra, and P(r), the pair distance distribution function, broadens and shows a transition from a unimodal to a bimodal distribution as the protein unfolds. The results indicate that the globular structure of SN is retained across the folding/unfolding transition although this structure is less compact and elongated relative to the native structure. Pressure-induced unfolding is initiated in the molecular dynamics simulations by inserting water molecules into the protein interior and applying pressure. The P(r) calculated from these simulations likewise broadens and shows a similar unimodal-to-bimodal transition with increasing pressure. The simulations also reveal that the bimodal P(r) for the pressure-unfolded state arises as the protein expands and forms two subdomains that effectively diffuse apart during initial stages of unfolding. Hydrophobic contact maps derived from the simulations show that water insertions into the protein interior and the application of pressure together destabilize hydrophobic contacts between these two subdomains. The findings support a mechanism for the pressure-induced unfolding of SN in which water penetration into the hydrophobic core plays a central role.     

Kinetic analysis of the acid and the alkaline unfolded states of staphylococcal nuclease

Thermodynamic analysis by differential scanning calorimetry shows that the folding/unfolding transition of staphylococcal nuclease is consistent with the two-state process. Stopped-flow kinetic measurements, monitoring the Trp140 fluorescence and covering five decades in time (2 ms to 300 s), indicate that the unfolding from pH 7.0 to 3.1 is monophasic (time constant 1.15 s) and from pH 7.0 to 12.2 is biphasic (time constants: one less than 2 ms and the other 0.6 s). However, the folding, either from pH 3.1 to 7.0 or from pH 12.2 to 7.0, is triphasic (time constants 150 ms, 850 ms and 30 s from acid, 90 ms, 565 ms and 33 s from alkaline). A simple sequential model, which agrees with the above observations for acidic folding/unfolding is, D3 in equilibrium D2 in equilibrium D1 in equilibrium N. The three Ds denote three sub-states of the unfolded state and N denotes the native state. These sub-states of D have similar enthalpy and tryptophan fluorescence, and their equilibrium cannot be shifted by temperature changes. However, they are kinetically distinctive. Data do not favor alternative mechanisms assuming parallel transitions of the three Ds to N, or complexity of the N state, or parallel transitions of sub-states of N1, N2 and N3 to D. Other more complex, branched or cyclic, kinetics are not considered because of the lack of evidence, pH dependence of the unfolding kinetics suggests that the unfolding is triggered by protonation of 0.8(+/- 0.3) ionizable groups, with a pKa of 3.9 or by deprotonation of 1.6(+/- 0.4) ionizable groups with pKa values near 10.5. Circular dichroisms indicate that these three D states retain nonrandom chain conformation. Possible role of these "chain conformation" in the protein folding is discussed.     

Calorimetric unfolding of the bimolecular and i-motif complexes of the human telomere complementary strand, d(C(3)TA(2))(4)

A combination of spectroscopic and calorimetric techniques is used to determine the unfolding thermodynamics of the complexes formed by the complementary sequence of the human telomere, d(C(3)TA(2))(4), in the pH range of 4.2 to 6. Calorimetric melting curves show biphasic transitions; both transitions are shifted to higher temperatures as the pH is decreased, indicative of cytosine protonation, which favors the formation of C*C(+) base pairs. Furthermore, the transition temperature, T(M), of the lower transition depends on strand concentration, while the T(M) of the higher transition is independent of strand concentration, indicating the following sequential melting: bimolecular complex(s)-->intramolecular complex-->random coil. The thermodynamic profiles for the formation of each complex, bimolecular and i-motif reveals small favorable free energy terms resulting from favorable enthalpy-unfavorable entropy compensations, uptake of protons, marginal uptake of counterions (i-motif) and marginal release of water molecules (i-motif). Furthermore, an enthalpy of 3.2 kcal/mol (bimolecular complex) and 5.0 kcal/mol (i-motif) is estimated for a single C*C(+)/C*C(+) base-pair stack.     

Implications from a network-based topological analysis of ubiquitin unfolding simulations

Background

The architectural organization of protein structures has been the focus of intense research since it can hopefully lead to an understanding of how proteins fold. In earlier works we had attempted to identify the inherent structural organization in proteins through a study of protein topology. We obtained a modular partitioning of protein structures with the modules correlating well with experimental evidence of early folding units or “foldons”. Residues that connect different modules were shown to be those that were protected during the transition phase of folding.

Methodology/Principal Findings

In this work, we follow the topological path of ubiquitin through molecular dynamics unfolding simulations. We observed that the use of recurrence quantification analysis (RQA) could lead to the identification of the transition state during unfolding. Additionally, our earlier contention that the modules uncovered through our graph partitioning approach correlated well with early folding units was vindicated through our simulations. Moreover, residues identified from native structure as connector hubs and which had been shown to be those that were protected during the transition phase of folding were indeed more stable (less flexible) well beyond the transition state. Further analysis of the topological pathway suggests that the all pairs shortest path in a protein is minimized during folding.

Conclusions

We observed that treating a protein native structure as a network by having amino acid residues as nodes and the non-covalent interactions among them as links allows for the rationalization of many aspects of the folding process. The possibility to derive this information directly from 3D structure opens the way to the prediction of important residues in proteins, while the confirmation of the minimization of APSP for folding allows for the establishment of a potentially useful proxy for kinetic optimality in the validation of sequence-structure predictions.     

Triplex formation at physiological pH by 5-Me-dC-N4-(spermine) [X] oligodeoxynucleotides: non protonation of N3 in X of X*G:C triad and effect of base mismatch/ionic strength on triplex stabilities.

Oligodeoxynucleotide (ODN) directed triplex formation has therapeutic importance and depends on Hoogsteen hydrogen bonds between a duplex DNA and a third DNA strand. T*A:T triplets are formed at neutral pH and C+*G:C are favoured at acidic pH. It is demonstrated that spermine conjugation at N4 of 5-Me-dC in ODNs 1-5 (sp-ODNs) imparts zwitterionic character, thus reducing the net negative charge of ODNs 1-5. sp-ODNs form triplexes with complementary 24mer duplex 8:9 show foremost stability at neutral pH 7.3 and decrease in stability towards lower pH, unlike the normal ODNs where optimal stability is found at an acidic pH 5.5. At pH 7.3, control ODNs 6 and 7 carrying dC or 5-Me-dC, respectively, do not show any triple helix formation. The stability order of triplex containing 5-Me-dC-N4-(spermine) with normal and mismatched duplex was found to be X*G:C approximately X*A:T > X*C:G > X*T:A. The hysteresis curve of sp-ODN triplex 3*8:9 indicated a better association with complementary duplex 8:9 as compared to unmodified ODN 6 in triplex 6*8:9. pH-dependent UV difference spectra suggest that N3 protonation is not a requirement for triplex formation by sp-ODN and interstrand interaction of conjugated spermine more than compensates for loss in stability due to absence of a single Hoogsteen hydrogen bond. These results may have importance in designing oligonucleotides for antigene applications.     

Kinetics of folding and unfolding of goat alpha-lactalbumin

The equilibrium unfolding and the kinetic folding and unfolding of goat alpha-lactalbumin (GLA) were studied by near- and far-ultraviolet circular dichroism (CD) and by stopped-flow fluorescence spectroscopy. Specifically, the influence of environmental conditions such as pH and Ca2+ binding was examined. Compared to the apo-form, the Ca2+-bound form was found to be strongly stabilized in equilibrium conditions at pH 7.5 and 25 degrees C. The kinetics of the refolding of apo-GLA show a major change of fluorescence intensity during the experimental dead-time, but this unresolved effect is strongly diminished in holo-GLA. In both cases, however, the chevron plots can adequately be fitted to a three-state model. Moreover, double-mix stopped-flow experiments showed that the native state (N) is reached through one major pathway without the occurrence of alternative tracks. In contrast to the homologous bovine alpha-lactalbumin (BLA), the compactness of GLA is strongly influenced by the presence of Ca2+ ions. Unlike the two-state transition observed in guanidine hydrochloride (GdnHCl)-induced equilibrium denaturation experiments at higher pH, an equilibrium intermediate state (I) is involved in denaturation at pH 4.5. In the latter case, analysis of the kinetic data makes clear that the intermediate and the unfolded states (U) show practically no Gibbs free energy difference and that they are in rapid equilibrium with each other. A possible explanation for these variations in stability and in folding characteristics with pH could be the degree of protonation of His107 that directly influences non-native interactions. Variation of environmental conditions and even small differences in sequence, therefore, can result in important effects on thermodynamic and folding parameters.     

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

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

The folding pathway of spectrin R17 from experiment and simulation: using experimentally validated MD simulations to characterize States hinted at by experiment

We present an experimental and computational analysis of the folding pathway of the 17th domain of chicken brain alpha-spectrin, R17. Wild-type R17 folds in a two-state manner and the chevron plot (plot of the logarithm of the observed rate constant against concentration of urea) shows essentially linear folding and unfolding arms. A number of mutant proteins, however, show a change in slope of the unfolding arm at high concentration of denaturant, hinting at complexity in the folding landscape. Through a combination of mutational studies and high temperature molecular dynamics simulations we show that the folding of R17 can be described by a model with two sequential transition states separated by an intermediate species. The rate limiting transition state for folding in water has been characterized both through experimental Phi-value analysis and by simulation. In contrast, a detailed analysis of the transition state predicted to dominate under highly denaturing conditions is only possible by simulation.