The sequential unfolding of ribonuclease A: detection of a fast initial phase in the kinetics of unfolding

Characterization of protein unfolding by fast cross-linking mass spectrometry using di-ortho-phthalaldehyde cross-linkers

Chemical cross-linking of proteins coupled with mass spectrometry is widely used in protein structural analysis. In this study we develop a class of non-hydrolyzable amine-selective di-ortho-phthalaldehyde (DOPA) cross-linkers, one of which is called DOPA2. Cross-linking of proteins with DOPA2 is 60-120 times faster than that with the N-hydroxysuccinimide ester cross-linker DSS. Compared with DSS cross-links, DOPA2 cross-links show better agreement with the crystal structures of tested proteins. More importantly, DOPA2 has unique advantages when working at low pH, low temperature, or in the presence of denaturants. Using staphylococcal nuclease, bovine serum albumin, and bovine pancreatic ribonuclease A, we demonstrate that DOPA2 cross-linking provides abundant spatial information about the conformations of progressively denatured forms of these proteins. Furthermore, DOPA2 cross-linking allows time-course analysis of protein conformational changes during denaturant-induced unfolding.

Network representation of conformational transitions between hidden intermediates of Rd-apocytochrome b562

The folding kinetics of Rd-apocytochrome b562 is two-state, but native-state hydrogen exchange experiments show that there are discrete partially unfolded (PUF) structures in equilibrium with the native state. These PUF structures are called hidden intermediates because they are not detected in kinetic experiments and they exist after the rate-limiting step. Structures of the mimics of hidden intermediates of Rd-apocytochrome b562 are resolved by NMR. Based upon their relative stability and structural features, the folding mechanism was proposed to follow a specific pathway (unfolded → rate-limiting transition state → PUF1 → PUF2 → native). Investigating the roles of equilibrium PUF structures in folding kinetics and their interrelationship not only deepens our understanding of the details of folding mechanism but also provides guides in protein design and prevention of misfolding. We performed molecular dynamics simulations starting from a hidden intermediate and the native state of Rd-apocytochrome b562 in explicit solvent, for a total of 37.18 μs mainly with Anton. We validated our simulations by detailed comparison with experimental data and other computations. We have verified that we sampled the post rate-limiting transition state region only. Markov state model was used to analyze the simulation results. We replace the specific pathway model with a network model. Transition-path theory was employed to calculate the net effective flux from the most unfolded state towards the most folded state in the network. The proposed sequential folding pathway via PUF1 then more stable, more native-like PUF2 is one of the routes in our network, but it is not dominant. The dominant path visits PUF2 without going through PUF1. There is also a route from PUF1 directly to the most folded state in the network without visiting PUF2. Our results indicate that the PUF states are not necessarily sequential in the folding. The major routes predicted in our network are testable by future experiments such as single molecule experiment.

Native contact density and nonnative hydrophobic effects in the folding of bacterial immunity proteins

The bacterial colicin-immunity proteins Im7 and Im9 fold by different mechanisms. Experimentally, at pH 7.0 and 10°C, Im7 folds in a three-state manner via an intermediate but Im9 folding is two-state-like. Accordingly, Im7 exhibits a chevron rollover, whereas the chevron arm for Im9 folding is linear. Here we address the biophysical basis of their different behaviors by using native-centric models with and without additional transferrable, sequence-dependent energies. The Im7 chevron rollover is not captured by either a pure native-centric model or a model augmented by nonnative hydrophobic interactions with a uniform strength irrespective of residue type. By contrast, a more realistic nonnative interaction scheme that accounts for the difference in hydrophobicity among residues leads simultaneously to a chevron rollover for Im7 and an essentially linear folding chevron arm for Im9. Hydrophobic residues identified by published experiments to be involved in nonnative interactions during Im7 folding are found to participate in the strongest nonnative contacts in this model. Thus our observations support the experimental perspective that the Im7 folding intermediate is largely underpinned by nonnative interactions involving large hydrophobics. Our simulation suggests further that nonnative effects in Im7 are facilitated by a lower local native contact density relative to that of Im9. In a one-dimensional diffusion picture of Im7 folding with a coordinate- and stability-dependent diffusion coefficient, a significant chevron rollover is consistent with a diffusion coefficient that depends strongly on native stability at the conformational position of the folding intermediate.

In order to fold correctly, a globular protein must avoid being trapped in wrong, i.e., nonnative conformations. Thus a biophysical account of how attractive nonnative interactions are bypassed by some amino acid sequences but not others is key to deciphering protein structure and function. We examine two closely related bacterial immunity proteins, Im7 and Im9, that are experimentally known to fold very differently: Whereas Im9 folds directly, Im7 folds through a mispacked conformational intermediate. A simple model we developed accounts for their intriguingly different folding kinetics in terms of a balance between the density of native-promoting contacts and the hydrophobicity of local amino acid sequences. This emergent principle is extensible to other biomolecular recognition processes.

An IMS-IMS threshold method for semi-quantitative determination of activation barriers: Interconversion of proline cis↔trans forms in triply protonated bradykinin

Collisional activation of selected conformations by multidimensional ion mobility spectrometry (IMS-IMS), combined with mass spectrometry (MS), is described as a method to determine semi-quantitative activation energies for interconversion of different structures of the nonapeptide bradykinin (BK, Arg-Pro-Pro-Gly-Phe-Ser-Pro-Phe-Arg). This analysis is based on a calibration involving collision-induced dissociation measurements of ions with known dissociation energies (i.e., "thermometer" ions) such as leucine enkephalin, BK, and amino acid-metal cation systems. The energetic barriers between six conformations of [BK+3H]3+ range from 0.23 ±0.01 to 0.55 ±0.03 eV. Prior results indicate that the major peaks in the IMS distributions correspond to specific combinations of cis and trans configurations of the three proline residues in the peptide sequence. The analysis allows us to directly assess pathways for specific transitions. The combination of structural assignments, experimentally determined barrier heights, onset of the quasi-equilibrium region, and dissociation threshold are used to derive a semi-quantitative potential energy surface for main features of [BK+3H]3+.

Transition paths, diffusive processes, and preequilibria of protein folding

Fundamental relationships between the thermodynamics and kinetics of protein folding were investigated using chain models of natural proteins with diverse folding rates by extensive comparisons between the distribution of conformations in thermodynamic equilibrium and the distribution of conformations sampled along folding trajectories. Consistent with theory and single-molecule experiment, duration of the folding transition paths exhibits only a weak correlation with overall folding time. Conformational distributions of folding trajectories near the overall thermodynamic folding/unfolding barrier show significant deviations from preequilibrium. These deviations, the distribution of transition path times, and the variation of mean transition path time for different proteins can all be rationalized by a diffusive process that we modeled using simple Monte Carlo algorithms with an effective coordinate-independent diffusion coefficient. Conformations in the initial stages of transition paths tend to form more nonlocal contacts than typical conformations with the same number of native contacts. This statistical bias, which is indicative of preferred folding pathways, should be amenable to future single-molecule measurements. We found that the preexponential factor defined in the transition state theory of folding varies from protein to protein and that this variation can be rationalized by our Monte Carlo diffusion model. Thus, protein folding physics is different in certain fundamental respects from the physics envisioned by a simple transition-state picture. Nonetheless, transition state theory can be a useful approximate predictor of cooperative folding speed, because the height of the overall folding barrier is apparently a proxy for related rate-determining physical properties.

Protein folding: then and now

Over the past three decades the protein folding field has undergone monumental changes. Originally a purely academic question, how a protein folds has now become vital in understanding diseases and our abilities to rationally manipulate cellular life by engineering protein folding pathways. We review and contrast past and recent developments in the protein folding field. Specifically, we discuss the progress in our understanding of protein folding thermodynamics and kinetics, the properties of evasive intermediates, and unfolded states. We also discuss how some abnormalities in protein folding lead to protein aggregation and human diseases.

An emerging concept of prion infections as a form of transmissible cerebral amyloidosis

Proteins are a major constituent of cells with specific biological functions. Besides the primary structure that is simply the sequence of amino acids that comprise a protein, the secondary structure represents the first step of folding defining its general conformation. The biological functions of proteins are directly dependent on the acquisition of their conformation. The same protein can have different stable states, which may participate with different functions in the cell. The amyloid diseases comprise Alzheimer's and Parkinson's diseases, type II diabetes mellitus and systemic amyloidosis. Amyloid fibers are insoluble, resistant to proteolysis and show an extremely high content of beta-sheet, in a very similar structure to the one observed among prion rods, associated to the transmissible spongiform encephalopathies. All these diseases are "infectious" in the sense that misfolded beta-sheeted conformers formed in a nucleation process in which preformed metastable oligomer acts as a seed to convert a normal isoform into an abnormal protein with a misfolded conformation. Only prion infections have a proven infectivity in a microbiological sense; some recent observations, however, detected the transmissibility of systemic amyloidosis by a prion-like mechanism among mice. Prions diseases and amyloidosis present many similar aspects of the so-called conformational diseases; according to this interpretation the prion infections could be considered as a form of transmissible cerebral amyloidosis.

Tuning lambda6-85 towards downhill folding at its melting temperature

The five-helix bundle lambda6-85* is a fast two-state folder. Several stabilized mutants have been reported to fold kinetically near-downhill or downhill. These mutants undergo a transition to two-state folding kinetics when heated. It has been suggested that this transition is caused by increased hydrophobicity at higher temperature. Here we investigate two histidine-containing mutants of lambda6-85* to see if a weaker hydrophobic core can extend the temperature range of downhill folding. The very stable lambdaHA is the fastest-folding lambda repressor to date (k(f)(-1) approximately k(obs)(-1)=2.3 micros at 44 degrees C). It folds downhill at low temperature, but transits back to two-state folding at its unfolding midpoint. lambdaHG has a weakened hydrophobic core. It is less stable than some slower folding mutants of lambda6-85*, and it has more exposed hydrophobic surface area in the folded state. This mutant nonetheless folds very rapidly, and has the non-exponential folding kinetics of an incipient downhill folder even at the unfolding midpoint (k(m)(-1) approximately 2 micros, k(a)(-1)=15 micros at 56 degrees C). We also compare the thermodynamic melting transition of lambdaHG with the nominal two-state folding mutant lambdaQG, which has a similar melting temperature. Unlike lambdaQG, lambdaHG yields fluorescence wavelength-dependent cooperativities and probe-dependent melting temperatures. This result combined with previous work shows that the energy landscapes of lambda repressor mutants support all standard folding mechanisms.

Denaturation of jack-bean urease by sodium n-dodecyl sulphate: A kinetic study below the critical micelle concentration

Kinetics of urease denaturation by anionic surfactant (sodium n-dodecyl sulphate, SDS) at concentrations below the critical micelle concentration (CMC) is investigated spectrophotometrically at neutral pH and the corresponding two-phase kinetic parameters of the process are estimated from a three-state reversible process using a binomial exponential relation based on the relaxation time method as: Using a prepared computer program, the experimental data are properly fitted into a binomial exponential relation, considering a two-phase denaturation pathway including a kinetically stable folded intermediate formed at SDS concentration of 1.1 mM. Forward and backward rate constants are estimated as: k(1)=0.2141+/-4.5 x 10(-3), k(2)=5.173 x 10(-3)+/-8.3 x 10(-5), k(-1)=0.09432+/-3.6 x 10(-4) and k(-2)=2.079 x 10(-3)+/-5.6 x 10(-5)s(-1) for the proposed mechanism. The rate-limiting step as well as the reaction coordinates in the denaturation mechanism are established. The mechanism involves formation of a kinetically stable folded native like intermediate through the electrostatic interactions. The intermediate was found to be more stable even than the native form (by about 9 kJmol(-1)) and still hexamer, because no loss of amplitude was observed. Electrophoresis experiments on the native and surfactant/urease complexes indicated a higher mobility for the kinetically folded native like intermediate.

Cooperativity in two-state protein folding kinetics

We present a solvable model that predicts the folding kinetics of two-state proteins from their native structures. The model is based on conditional chain entropies. It assumes that folding processes are dominated by small-loop closure events that can be inferred from native structures. For CI2, the src SH3 domain, TNfn3, and protein L, the model reproduces two-state kinetics, and it predicts well the average Phi-values for secondary structures. The barrier to folding is the formation of predominantly local structures such as helices and hairpins, which are needed to bring nonlocal pairs of amino acids into contact.

Fast-folding protein kinetics, hidden intermediates, and the sequential stabilization model

Do two-state proteins fold by pathways or funnels? Native-state hydrogen exchange experiments show discrete nonnative structures in equilibrium with the native state. These could be called hidden intermediates (HI) because their populations are small at equilibrium, and they are not detected in kinetic experiments. HIs have been invoked as disproof of funnel models, because funnel pictures appear to indicate (1) no specific sequences of events in folding; (2) a continuum, rather than a discrete ladder, of structures; and (3) smooth landscapes. In the present study, we solve the exact dynamics of a simple model. We find, instead, that the present microscopic model is indeed consistent with HIs and transition states, but such states occur in parallel, rather than along the single pathway predicted by the sequential stabilization model. At the microscopic level, we observe a huge multiplicity of trajectories. But at the macroscopic level, we observe two pathways of specific sequences of events that are relatively traditional except that they are in parallel, so there is not a single reaction coordinate. Using singular value decomposition, we show an accurate representation of the shapes of the model energy landscapes. They are highly complex funnels.

Quantitative study of polymer conformation and dynamics by single-particle tracking

We present a new method for analyzing the dynamics of conformational fluctuations of individual flexible polymer molecules. In single-particle tracking (SPT), one end of the polymer molecule is tethered to an immobile substratum. A microsphere attached to the other end serves as an optical marker. The conformational fluctuations of the polymer molecule can be measured by optical microscopy via the motion of the microsphere. The bead-and-spring theory for polymer dynamics is further developed to account for the microsphere, and together the measurement and the theory yield quantitative information about molecular conformations and dynamics under nonperturbing conditions. Applying the method to measurements carried out on DNA molecules provides information complementary to recent studies of single DNA molecules under extensional force. Combining high precision measurements with the theoretical analysis presented here creates a powerful tool for studying conformational dynamics of biological and synthetic macromolecules at the single-molecule level.

Protein folding: matching theory and experiment

The impact of folding funnels and folding simulations on the way experimentalists interpret results is examined. The image of the transition state has changed from a unique species that has a strained configuration, with a correspondingly high free energy, to a more ordinary folding intermediate, whose balance between limited conformational entropy and stabilizing contacts places it at the top of the free energy barrier. Evidence for a broad transition barrier comes from studies showing that mutations can change the position of the barrier. The main controversial issue now is whether populated folding intermediates are productive on-pathway intermediates or dead-end traps. Direct experimental evidence is needed. Theories suggesting that populated intermediates are trapped in a glasslike state are usually based on mechanisms which imply that trapping would only be extremely short-lived (e.g., nanoseconds) in water at 25 degrees C. There seems to be little experimental evidence for long-lived trapping in monomers, if folding aggregates are excluded. On the other hand, there is good evidence for kinetic trapping in dimers. alpha-Helix formation is currently the fastest known process in protein folding, and incipient helices are present at the start of folding. Fast helix formation has the effect of narrowing drastically the choice of folding routes. Thus helix formation can direct folding. It changes the folding metaphor from pouring liquid down a folding funnel to a train leaving a switchyard with only a few choices of exit tracks.

Protein folding in the landscape perspective: chevron plots and non-Arrhenius kinetics

We use two simple models and the energy landscape perspective to study protein folding kinetics. A major challenge has been to use the landscape perspective to interpret experimental data, which requires ensemble averaging over the microscopic trajectories usually observed in such models. Here, because of the simplicity of the model, this can be achieved. The kinetics of protein folding falls into two classes: multiple-exponential and two-state (single-exponential) kinetics. Experiments show that two-state relaxation times have "chevron plot" dependences on denaturant and non-Arrhenius dependences on temperature. We find that HP and HP+ models can account for these behaviors. The HP model often gives bumpy landscapes with many kinetic traps and multiple-exponential behavior, whereas the HP+ model gives more smooth funnels and two-state behavior. Multiple-exponential kinetics often involves fast collapse into kinetic traps and slower barrier climbing out of the traps. Two-state kinetics often involves entropic barriers where conformational searching limits the folding speed. Transition states and activation barriers need not define a single conformation; they can involve a broad ensemble of the conformations searched on the way to the native state. We find that unfolding is not always a direct reversal of the folding process.

From Levinthal to pathways to funnels

While the classical view of protein folding kinetics relies on phenomenological models, and regards folding intermediates in a structural way, the new view emphasizes the ensemble nature of protein conformations. Although folding has sometimes been regarded as a linear sequence of events, the new view sees folding as parallel microscopic multi-pathway diffusion-like processes. While the classical view invoked pathways to solve the problem of searching for the needle in the haystack, the pathway idea was then seen as conflicting with Anfinsen's experiments showing that folding is pathway-independent (Levinthal's paradox). In contrast, the new view sees no inherent paradox because it eliminates the pathway idea: folding can funnel to a single stable state by multiple routes in conformational space. The general energy landscape picture provides a conceptual framework for understanding both two-state and multi-state folding kinetics. Better tests of these ideas will come when new experiments become available for measuring not just averages of structural observables, but also correlations among their fluctuations. At that point we hope to learn much more about the real shapes of protein folding landscapes.

Ribonucleases from the extreme thermophilic archaebacterium S. solfataricus

A purification procedure consisting of DEAE-Sephacel chromatography, heparin-Sepharose CL-6B chromatography and Mono-S chromatography led to the isolation of three proteins endowed with RNase activity from the thermoacidophilic archaebacterium Sulfolobus solfataricus. They were referred to as p1, p2 and p3, according to their elution order from the Mono-S column. Complete amino acid sequence of p2 and partial sequence of p3 displayed high sequence similarity to the 7-kDa DNA-binding proteins previously isolated in Sulfolobus strains [Choli, T., Wittman-Liebold, B. & Reinhardt, R. (1988) J. Biol. Chem. 263, 7087-7093]. The molecular mass of p2, calculated from sequence data, was 7.02 kDa, which compares fairly well with the value of 7.4 kDa determined by SDS/PAGE. Gel filtration of the molecule under native conditions displayed, however, a largely prevailing form with an assessed molecular mass of 13.0 kDa, which points to a dimeric structure. Kinetic characterization of protein p2 showed a broad pH optimum in the range 6.7-7.6 using yeast RNA as substrate; also, it was shown that activity was unaffected by EDTA, Mg2+ and phosphate. The enzyme did not accept as substrate any homopolyribonucleotide, which points to a rather narrow substrate specificity. This was also confirmed by incubating p2 with tRNA(fMet)Met (fMet, N-formylmethionine) from Escherichia coli: the hydrolysis products were thus identified as 3'-phosphooligonucleotides.

Kinetics of reversible protein denaturation

The kinetics of denaturation and renaturation of Aplysia Myoglobin by temperature-jump and stopped-flow are reported. The time course of the unfolding and refolding reveals the presence of significant amounts of intermediates both in absence and in the presence of alcohols. A linear three states reaction mechanism is proposed and discussed.

A quantitative treatment of the kinetics of the folding transition of ribonuclease A

New experimental data and a quantitative theoretical treatment are given for the kinetics of the thermal folding transition of ribonuclease A at pH 3.0. A three-species mechanism is used as a starting point for the analysis: U1 (slow) in equilibrium U2(fast) in equilibrium N, where U1 and U2 are two forms of the unfolded enzyme with markedly different rates of refolding and N is the native enzyme. This mechanism is based on certain facts established in previous studies of refolding. The kinetics of unfolding and refolding show two phases a fast phase and a slow phase, over a range of temperatures extending above the transition midpoint, Tm. The three-species mechanism can be used in this range. At higher temperatures a new much faster kinetic phase is also observed corresponding to the transient formation of a new intermediate (I). Although the general solution for a four-species mechanism is complex it is not difficult to extend the three-species analysis for the special case found here, in which the fast reaction (I in equilibrium N) is well separated from the other two reactions. At temperatures below the transition zone the slow phase of refolding becomes kinetically complex. No attempt has been made to extend the analysis to include this effect. The basic test of the three-state analysis is the prediction as a function of temperature of alpha2, the relative amplitude of the fast phase, both for unfolding and refolding. At temperatures above Tm for which the three-state analysis must be extended to include the new intermediate I, a crresponding quanitity alpha2(cor) is predicted and compared with measured values. Data used in the three-state prediction are values of tau2 and tau1, the time constants of the fast and slow kinetic phases, plus a single value of alpha2 measured when tau2 and tau1 are well separated. The observed and predicted values of alpha2 agree within experimental error. The analysis predicts correctly that, for these experiments, alpha2 should have the same value in unfolding as in refolding in the final conditions. The analysis also predicts satisfactorily the equilibrium transition curve from kinetic data alone. Four striking properties of the kinetics are explained or correlated by the analysis: (a) the drop in alpha2 to a minimum near Tm as well as the delayed rise in alpha2 above Tm;(b) the vanishing of alpha1 above the transition zone; (c) the sharp drop in tau1 inside the transition zone followed by a partial leveling off outside this zone; and (d) the passage of tau2 through a maximum near Tm. Through a comparison of observed and predicted values of alpha2, the analysis also rules out the alternative three-species mechanism U1 (slow) in equilibrium N (fast) in equilibrium U2. Finally, the temperature dependence of the amplitude for the fast reaction (I in equilibrium N) is discussed; the behavior of I is like that of U2 and I may be an unfolded species populated at equilibrium...

An acid induced conformational transition of denatured cytochrome c in urea and guanidine hydrochloride solutions

Previous work has shown that at neutral pH ferricytochrome c (horse heart) retains certain residual structures in concentrated solutions of urea or guanidine hydrochloride (Tsong, T. Y. (1974), J. Biol. Chem. 249, 1988). Present studies reveal that cooperative unfolding of these residual structures can be achieved by acidification of the protein to pH 4 in 9 M urea but can only be partially achieved in a 6 M guanidine hydrochloride solution. The evidence that the residual structures unfold in 9 M urea upon acidification is twofold. (1) Further uncoupling of the Trp-59-heme interaction occurs; this is reflected in the intensification of the tryptophan fluorescence from 55 to 90 percent relative to that of free tryptophan in the same solvent. (2) The intrinsic viscosity of the protein solution increases from 15.0 to 21 ml/g. The acidification also induces a spin-state transformation of the heme group at pH 5 both in urea and in guanidine hydrochloride. Acidic titration of the protein in urea and guanidine hydrochloride indicates that the unfolding involves the absorption of a single proton. However, the kinetics of the spin-state transformation are triphasic. These results suggest that the displacement of the ligand His-18 by a solvent molecule and the subsequent disintegration of the residual structures are complex processes and involve at least three kinetic steps. The ineffectiveness of guanidine hydrochloride as a denaturant for ferricytochrome c is shown to be due to the presence of the high concentration of Cl minus which can stabilize certain elements of the protein structure.

Nuclear magnetic resonance study of the thermal denaturation of ribonuclease A: implications for multistate behavior at low pH

The thermal denaturation of ribonuclease A has been studied by use of Fourier transform nuclear magnetic resonance by monitoring the imidazole C-2 proton resonances of the histidine residues as a function of temperature at pH 1.3. As the temperature is raised, a slow chemical exchange process results in the disappearance of the peaks corresponding to the native conformation and the appearance of a single peak corresponding to histidine in the denatured state. The disappearance of the native peaks is not simultaneous, implying that at least two regions of the molecule denature at different temperatures. Also, fast chemical exchange processes result in small chemical shifts that appear to be related to local conformational changes. The observed phenomena have been shown to be reversible by the measurement of absorbance at 278 nm, enzyme activity, and nuclear magnetic resonance spectroscopy. The results of this equilibrium study support a multistate denaturation mechanism for ribonuclease A at pH 1.3.

Properties of the refolding and unfolding reactions of ribonuclease A

Both the refolding kinetics and unfolding kinetics of ribonuclease A have been measured at the same final conditions, as a function of temperature at pH 3.9, by stopped-flow (pH-jump) experiments; absorbance changes at 240 and 286.5 nm were measured. Refolding follows first-order kinetics in the upper two-thirds of the thermal transition zone. Under the same conditions, the unfolding kinetics are biphasic; the terminal phase has the same rate constant as refolding. The biphasic kinetics of unfolding demonstrate the presence of intermediate states. Since both the refolding and unfolding kinetics are consistent with a simple sequential model, the intermediates satisfy kinetic criteria for being on the direct pathway of unfolding. At temperatures just above the transition zone, the fast phase of unfolding becomes the major kinetic phase. The rate of the slow unfolding reaction increases rapidly with temperature, and approaches the average rate of the fast phase at temperatures just above the transition zone. The entire set of kinetic results can be reproduced semiquantitatively by assignment of values to four parameters in a cooperative sequential model. However, reasons are given for the belief that this simple model will have to be generalized before it can give a realistic description of the kinetics of the unfolding reaction.