The core of apomyoglobin E-form folds at the diffusion limit

Heterogeneity in the Folding of Villin Headpiece Subdomain HP36

Small single domain proteins that fold on the microsecond time scale have been the subject of intense interest as models for probing the complexity of folding energy landscapes. The villin headpiece subdomain (HP36) has been extensively studied because of its simple three helix structure, ultrafast folding lifetime of a few microseconds, and stable native fold. We have previously shown that folding as measured by a single ¹³C═¹⁸O isotopic label on residue A57 in helix 2 occurs at a different rate than that measured by global probes of folding, indicating noncooperative complexity in the folding of HP36. In order to determine whether this complexity reflects intermediates or parallel pathways over a small activation barrier, ¹³C═¹⁸O labels were individually incorporated at six different positions in HP36, including into all 3 helices. The equilibrium thermal unfolding transitions and the folding/unfolding dynamics were monitored using the unique IR signature of the ¹³C═¹⁸O label by temperature dependent FTIR and temperature jump IR spectroscopy, respectively. Equilibrium experiments reveal that the ¹³C═¹⁸O labels at different positions in HP36 show drastic differences in the midpoint of their transitions ( T_m), ranging from 45 to 67 °C. Heterogeneity is also observed in the relaxation kinetics; there are differences in the microsecond phase when different labeled positions are probed. At a final temperature of 45 °C, the relaxation rate for ¹³C═¹⁸O A57 is 2.4e + 05 s^-1 whereas for ¹³C═¹⁸O L69 HP36 the relaxation rate is 5.1e + 05 s^-1, two times faster. The observation of site-dependent midpoints for the equilibrium unfolding transitions and differences in the relaxation rates of the labeled positions enables us to probe the progressive accumulation of the folded structure, providing insight into the microscopic details of the folding mechanism.

Microsecond and nanosecond polyproline II helix formation in aqueous nanodrops measured by mass spectrometry

The 1.5 μs and <400 ns time constants for the formation of polyproline II helix structures in 21 and 16 residue peptides, respectively, are measured using rapid mixing from theta-glass emitters coupled with mass spectrometry. Results from these studies should serve as useful benchmarks for comparison with computational simulation results.

Ultrafast (1 μs) Mixing and Fast Protein Folding in Nanodrops Monitored by Mass Spectrometry

The use of theta-glass emitters and mass spectrometry to monitor reactions that occur as fast as one μs is demonstrated. Acidified aqueous solutions containing unfolded proteins are mixed with aqueous ammonium acetate solutions to increase the solution pH and induce protein folding during nanoelectrospray ionization. Protein charge-state distributions show the extent to which folding occurs, and reaction times are obtained from known protein folding time constants. Shorter reaction times are obtained by decreasing the solution flow rate, and reaction times between 1.0 and 22 μs are obtained using flow rates between 48 and 2880 pL/s, respectively. Remarkably similar reaction times are obtained for three different proteins (Trp-cage, myoglobin, and cytochrome c) with folding time constants that differ by more than an order of magnitude (4.1, 7, and 57 μs, respectively), indicating that the reaction times obtained using rapid mixing from theta-glass emitters are independent of protein identity. A folding time constant of 2.2 μs is obtained for the formation of a β-hairpin structure of renin substrate tetradecapeptide, which is the fastest folding event measured using a rapid mixing technique. The 1.0 μs reaction time obtained here is about an order of magnitude lower than the shortest reaction time probed using a conventional mixer (8 μs). Moreover, this fast reaction time is obtained with a 48 pL/s flow rate, which is 2000-times less than the flow rate required to obtained the 8 μs reaction time using a conventional mixer. These results indicate that rapid mixing with theta-glass emitters can be used to access significantly faster reaction times while consuming substantially less sample than in conventional mixing apparatus.

Dimer domain swapping versus monomer folding in apo-myoglobin studied by molecular simulations

Using a coarse-grained symmetrized Go model, we performed a series of folding simulations of two apo-myoglobin molecules restrained at a high density, addressing competition of formation of a domain-swapped dimer with folding to two monomer structures.

Investigating protein folding and unfolding in electrospray nanodrops upon rapid mixing using theta-glass emitters

Theta-glass emitters are used to rapidly mix two solutions to induce either protein folding or unfolding during nanoelectrospray (nanoESI). Mixing acid-denatured myoglobin with an aqueous ammonium acetate solution to increase solution pH results in protein folding during nanoESI. A reaction time and upper limit to the droplet lifetime of 9 ± 2 μs is obtained from the relative abundance of the folded conformer in these rapid mixing experiments compared to that obtained from solutions at equilibrium and a folding time constant of 7 μs. Heme reincorporation does not occur, consistent with the short droplet lifetime and the much longer time constant for this process. Similar mixing experiments with acid-denatured cytochrome c and the resulting folding during nanoESI indicate a reaction time of between 7 and 25 μs depending on the solution composition. The extent of unfolding of holo-myoglobin upon rapid mixing with theta-glass emitters is less than that reported previously (Fisher et al. Anal. Chem.2014, 86, 4581−458824702054), a result that is attributed to the much smaller, ∼1.5 μm, average o.d. tips used here. These results indicate that the time frame during which protein folding or unfolding can occur during nanoESI depends both on the initial droplet size, which can be varied by changing the emitter tip diameter, and on the solution composition. This study demonstrates that protein folding or unfolding processes that occur on the ∼10 μs time scale can be readily investigated using rapid mixing with theta-glass emitters combined with mass spectrometry.

A simple three-dimensional-focusing, continuous-flow mixer for the study of fast protein dynamics

We present a simple, yet flexible microfluidic mixer with a demonstrated mixing time as short as 80 μs that is widely accessible because it is made of commercially available parts. To simplify the study of fast protein dynamics, we have developed an inexpensive continuous-flow microfluidic mixer, requiring no specialized equipment or techniques. The mixer uses three-dimensional, hydrodynamic focusing of a protein sample stream by a surrounding sheath solution to achieve rapid diffusional mixing between the sample and sheath. Mixing initiates the reaction of interest. Reactions can be spatially observed by fluorescence or absorbance spectroscopy. We characterized the pixel-to-time calibration and diffusional mixing experimentally. We achieved a mixing time as short as 80 μs. We studied the kinetics of horse apomyoglobin (apoMb) unfolding from the intermediate (I) state to its completely unfolded (U) state, induced by a pH jump from the initial pH of 4.5 in the sample stream to a final pH of 2.0 in the sheath solution. The reaction time was probed using the fluorescence of 1-anilinonaphthalene-8-sulfonate (1,8-ANS) bound to the folded protein. We observed unfolding of apoMb within 760 μs, without populating additional intermediate states under these conditions. We also studied the reaction kinetics of the conversion of pyruvate to lactate catalyzed by lactate dehydrogenase using the intrinsic tryptophan emission of the enzyme. We observe sub-millisecond kinetics that we attribute to Michaelis complex formation and loop domain closure. These results demonstrate the utility of the three-dimensional focusing mixer for biophysical studies of protein dynamics.

An instrument for fast acquisition of fluorescence decay curves at picosecond resolution designed for "double kinetics" experiments: application to fluorescence resonance excitation energy transfer study of protein folding

The information obtained by studying fluorescence decay of labeled biopolymers is a major resource for understanding the dynamics of their conformations and interactions. The lifetime of the excited states of probes attached to macromolecules is in the nanosecond time regime, and hence, a series of snapshot decay curves of such probes might - in principle - yield details of fast changes of ensembles of labeled molecules down to sub-microsecond time resolution. Hence, a major current challenge is the development of instruments for the low noise detection of fluorescence decay curves within the shortest possible time intervals. Here, we report the development of an instrument, picosecond double kinetics apparatus, that enables recording of multiple fluorescence decay curves with picosecond excitation pulses over wide spectral range during microsecond data collection for each curve. The design is based on recording and averaging multiphoton pulses of fluorescence decay using a fast 13 GHz oscilloscope during microsecond time intervals at selected time points over the course of a chemical reaction or conformational transition. We tested this instrument in a double kinetics experiment using reference probes (N-acetyl-tryptophanamide). Very low stochastic noise level was attained, and reliable multi-parameter analysis such as derivation of distance distributions from time resolved FRET (fluorescence resonance excitation energy transfer) measurements was achieved. The advantage of the pulse recording and averaging approach used here relative to double kinetics methods based on the established time correlated single photon counting method, is that in the pulse recording approach, averaging of substantially fewer kinetic experiments is sufficient for obtaining the data. This results in a major reduction in the consumption of labeled samples, which in many cases, enables the performance of important experiments that were not previously feasible.

Negative regulation of fibrin polymerization and clot formation by nanoparticles of silver

Thrombolytic therapy in acute stroke has reduced ischemia; however, it is also associated with increased incidence of intracerebral hemorrhage and expanding stroke. Platelets and fibrin are the major components of thrombi. Since fibrin is available in large concentration at lesion sites and in all types of thrombi, it is an obvious target for majority of antithrombotic therapies. Previously we have demonstrated innate antiplatelet properties with nanosilver. It can effectively prevent platelet activation in response to physiological agonists, under both in vitro as well as ex vivo conditions, and immobilize and stabilize proteins. Here we report for the first time that nanosilver can significantly retard fibrin polymerization kinetics both in pure and plasma-incorporated systems and hence can impede thrombus formation. We also discuss the conformational changes ensued upon fibrinogen following interaction with nanosilver. Together with its inherent antiplatelet and antibacterial properties, capacity to inhibit fibrin polymerization can open up possibilities of newer biomedical application and research potential involving silver nanoparticles.

Solvation and desolvation dynamics in apomyoglobin folding monitored by time-resolved infrared spectroscopy

Solvation and desolvation dynamics around helices during the kinetic folding process of apomyoglobin (apoMb) were investigated by using time-resolved infrared (IR) spectroscopy based on continuous-flow rapid mixing devices and an IR microscope. The folding of apoMb can be described by the collapse and search mechanism, in which the initial collapse occurring within several hundreds of microseconds is followed by the search for the correct secondary and tertiary structures. The time-resolved IR measurements showed a significant increase in solvated helix possessing a component of amide I' at 1633 cm(-1) within 100 mus after initiating the folding by a pD jump from pD2.2 to 6.0. In contrast, there was a minor increase in buried helices having amide I' at 1652 cm(-1) in this time domain. The observations demonstrate that the initially collapsed conformation of apoMb possesses a large amount of solvated helices, and suggest that much water is retained inside the collapsed domain. The contents of solvated and buried helices decrease and increase, respectively, in the time domain after the collapse, showing that the stepwise desolvation around helices is associated with the conformational search process. Interestingly, the largest changes in solvated and buried helices were observed at the final rate-limiting step of the apoMb folding. The persistence of the solvated helix until the final stage of apoMb folding suggests that the dissociation of hydrogen bonds between water and main-chain amides contributes to the energy barrier in the rate-determining step of the folding.

Analysis of heterogeneous fluorescence decays in proteins. Using fluorescence lifetime of 8-anilino-1-naphthalenesulfonate to probe apomyoglobin unfolding at equilibrium

The solvatochromic fluorescent dye 8-anilino-1-naphthalenesulfonate (ANS) is one of the popular probes of protein folding. Folding kinetics is tracked with ANS fluorescence intensity, usually interpreted as a reflection of protein structure-the hydrophobicity of the binding environments. Such simplistic view overlooks the complicated nature of ANS-protein complexes: the fluorescence characteristics are convoluted results of the ground state populational distribution of the probe-protein complex, the structural changes in the protein and the excited state photophysics of the probe. Understanding of the interplay of these aspects is crucial in accurate interpretation of the protein dynamics. In this work, the fluorescence decay of ANS complexed with apomyoglobin at different conformations denatured by pH is modeled. The fluorescence decay of the ANS-apomyoglobin complex contains information on not only apomyoglobin structure but also molecular populational distributions. The challenge in modeling fluorescence decay profiles originates from the convolution of heterogeneous binding and excited-state relaxation of the fluorescent probe. We analyzed frequency-domain fluorescence lifetime data of ANS-apomyoglobin with both maximum entropy methods (MEM) and nonlinear least squares methods (NLLS). MEM recovers a model of two expanding-and-merging lifetime distributions for ANS-apomyoglobin in the equilibrium transition from the native (N) through an intermediate (I-1) to the acid-unfolded state U(A). At pH 6.5 and above, when apomyoglobin is mostly populated at the N-state, ANS-apomyoglobin emits a predominant long-lifetime fluorescence from a relaxed charge transfer state S(1,CT) of ANS, and a short-lifetime fluorescence that is mainly from a nascent excited-state S(1,np) of ANS stabilized by the strong ANS-apomyoglobin interaction. Lowering the pH diminishes the contribution from the S(1,np) state. Meanwhile, more protein molecules become populated at the U(A) state, which exhibits a short lifetime that is not distinguishable from the S(1,np) state. At pH 3.4, when the population of the U(A) becomes significant, the short-lifetime fluorescence comes predominantly from ANS binding to the U(A). Further lowering the pH leads to more exposure of the bound ANS. The long lifetime shifts toward and finally merges with the short lifetime and becomes one broad distribution that stands for ANS binding to the U(A) below pH 2.4. The above expanding-and-merging model is consistent with F-statistic analysis of NLLS models. The consistency of this model with the knowledge from the literature, as well as the continuity of the decay parameters changing upon experimental conditions are also crucial in drawing the conclusions.

Steady-state and transient ultraviolet resonance Raman spectrometer for the 193-270 nm spectral region

We describe a state-of-the-art tunable ultraviolet (UV) Raman spectrometer for the 193-270 nm spectral region. This instrument allows for steady-state and transient UV Raman measurements. We utilize a 5 kHz Ti-sapphire continuously tunable laser (approximately 20 ns pulse width) between 193 nm and 240 nm for steady-state measurements. For transient Raman measurements we utilize one Coherent Infinity YAG laser to generate nanosecond infrared (IR) pump laser pulses to generate a temperature jump (T-jump) and a second Coherent Infinity YAG laser that is frequency tripled and Raman shifted into the deep UV (204 nm) for transient UV Raman excitation. Numerous other UV excitation frequencies can be utilized for selective excitation of chromophoric groups for transient Raman measurements. We constructed a subtractive dispersion double monochromator to minimize stray light. We utilize a new charge-coupled device (CCD) camera that responds efficiently to UV light, as opposed to the previous CCD and photodiode detectors, which required intensifiers for detecting UV light. For the T-jump measurements we use a second camera to simultaneously acquire the Raman spectra of the water stretching bands (2500-4000 cm(-1)) whose band-shape and frequency report the sample temperature.

Collapse and search dynamics of apomyoglobin folding revealed by submillisecond observations of alpha-helical content and compactness

The characterization of protein folding dynamics in terms of secondary and tertiary structures is important in elucidating the features of intraprotein interactions that lead to specific folded structures. Apomyoglobin (apoMb), possessing seven helices termed A-E, G, and H in the native state, has a folding intermediate composed of the A, G, and H helices, whose formation in the submillisecond time domain has not been clearly characterized. In this study, we used a rapid-mixing device combined with circular dichroism and small-angle x-ray scattering to observe the submillisecond folding dynamics of apoMb in terms of helical content (f(H)) and radius of gyration (R(g)), respectively. The folding of apoMb from the acid-unfolded state at pH 2.2 was initiated by a pH jump to 6.0. A significant collapse, corresponding to approximately 50% of the overall change in R(g) from the unfolded to native conformation, was observed within 300 micros after the pH jump. The collapsed intermediate has a f(H) of 33% and a globular shape that involves >80% of all its atoms. Subsequently, a stepwise helix formation was detected, which was interpreted to be associated with a conformational search for the correct tertiary contacts. The characterized folding dynamics of apoMb indicates the importance of the initial collapse event, which is suggested to facilitate the subsequent conformational search and the helix formation leading to the native structure.

Primary folding dynamics of sperm whale apomyoglobin: core formation

The structure, thermodynamics, and kinetics of heat-induced unfolding of sperm whale apomyoglobin core formation have been studied. The most rudimentary core is formed at pH(*) 3.0 and up to 60 mM NaCl. Steady state for ultraviolet circular dichroism and fluorescence melting studies indicate that the core in this acid-destabilized state consists of a heterogeneous composition of structures of approximately 26 residues, two-thirds of the number involved for horse heart apomyoglobin under these conditions. Fluorescence temperature-jump relaxation studies show that there is only one process involved in Trp burial. This occurs in 20 micro s for a 7 degrees jump to 52 degrees C, which is close to the limits placed by diffusion on folding reactions. However, infrared temperature jump studies monitoring native helix burial are biexponential with times of 5 micro s and 56 micro s for a similar temperature jump. Both fluorescence and infrared fast phases are energetically favorable but the slow infrared absorbance phase is highly temperature-dependent, indicating a substantial enthalpic barrier for this process. The kinetics are best understood by a multiple-pathway kinetics model. The rapid phases likely represent direct burial of one or both of the Trp residues and parts of the G- and H-helices. We attribute the slow phase to burial and subsequent rearrangement of a misformed core or to a collapse having a high energy barrier wherein both Trps are solvent-exposed.

Two-state expansion and collapse of a polypeptide

The initial phase of folding for many proteins is presumed to be the collapse of the polypeptide chain from expanded to compact, but still denatured, conformations. Theory and simulations suggest that this collapse may be a two-state transition, characterized by barrier-crossing kinetics, while the collapse of homopolymers is continuous and multi-phasic. We have used a laser temperature-jump with fluorescence spectroscopy to measure the complete time-course of the collapse of denatured cytochrome c with nanosecond time resolution. We find the process to be exponential in time and thermally activated, with an apparent activation energy approximately 9 k(B)T (after correction for solvent viscosity). These results indicate that polypeptide collapse is kinetically a two-state transition. Because of the observed free energy barrier, the time scale of polypeptide collapse is dramatically slower than is predicted by Langevin models for homopolymer collapse.

A statistical appraisal of native state hydrogen exchange data: evidence for a burst phase continuum?

For a number of proteins, folding occurs via the rapid accumulation of secondary and tertiary structural features in a so-called burst phase, preceding the relatively slow, highly activated transition leading to the native state. A fundamental question is: do these burst phase reactions comprise two phase-separated thermodynamic states or a continuum of states? Ribonuclease HI (RNase H) from Escherichia coli and phage T4 lysozyme (T4L) both exhibit such a phenomenon. Native-state hydrogen exchange (NHX) data have been collected for these proteins, providing residue-specific free energies and m-values (a measure of hydrocarbon solvation) for the manifold of partially unfolded, exchange-competent forms that are accessible from the native state (DeltaG(sg) and m(sg), where the sg subscript denotes sub-global). There is good evidence that these parameters pertain to exchange-competent species comprising the burst phase observed in the global folding kinetics. We combine the results from the global folding kinetics of these proteins with a statistical analysis of their NHX parameters to determine if the distribution of experimental (m(sg), DeltaG(sg)) values derive from a mechanism where the burst phase is two-state. For RNase H, this analysis demonstrates that the burst phase of this protein is not two-state; the results imply a distribution of states, m and DeltaG exhibiting a linear functional relationship consistent with the global folding parameters. For T4L, it is difficult to distinguish the observed distribution of m(sg), DeltaG(sg) values from that expected for a mechanism where the burst phase is two-state. The results for RNase H* lend support for the idea that the burst phase reaction of this protein comprises a continuum of states. This has important implications for how we model the process of structural acquisition in protein folding reactions.

A new folding intermediate of apomyoglobin from Aplysia limacina: stepwise formation of a molten globule

Apomyoglobin from Aplysia limacina (al-apoMb), despite having only 20 % sequence identity with the more commonly studied mammalian globins (m-apoMbs), properties which result in an increased number of hydrophobic contacts and a loss of most internal salt bridges, shares a number of features of their folding profiles. We show here that it contains an unusually stable core which resists unfolding even at 70 degrees C. The equilibrium intermediate (I(T)) at this high temperature is distinct from the acid unfolded state I(A) which has many properties in common with the acid intermediate observed for the mammalian apoproteins (I(AGH)). It contains a smaller amount of secondary structure (27 % alpha-helical instead of 35 %) and is more highly solvated as evidenced from its fluorescence spectrum (lambda(max)=344 nm instead of 338 nm). Its stability is greatly increased (DeltaDeltaG(w)=-6.75 kcal mol(-1)) in the presence of high salt (2 M KCl), lending support to the view that hydrophobic interactions are responsible for its stability. Kinetic data show classical two-state kinetics between I(A) and the folded state both in the presence and absence of salt. Both I(A) and I(T) can be populated within the dead time of the stopped-flow apparatus, since initiating the refolding reaction from I(T) or I(A) rather than the completely unfolded state does not affect the observed refolding time-course. Our conclusion is that al-apoMb, as other "apo" proteins (including for example alpha-lactalbumin in the absence of Ca(2+)), may be described as "uncoupled" with an unusually high and exploitable tendency to populate partially folded states.

Two-state expansion and collapse of a polypeptide

The initial phase of folding for many proteins is presumed to be the collapse of the polypeptide chain from expanded to compact, but still denatured, conformations. Theory and simulations suggest that this collapse may be a two-state transition, characterized by barrier-crossing kinetics, while the collapse of homopolymers and random heteropolymers is continuous and multi-phasic. A new rapid-mixing flow technique has been used to resolve the late stages of polypeptide collapse, at time scales >/=45 microseconds. We have used a laser temperature-jump with fluorescence spectroscopy to resolve the complete time-course of the collapse of denatured cytochrome c with nanosecond time resolution. We find the process to be exponential in time and thermally activated, with an apparent activation energy approximately 9 k(B)T (after correction for solvent viscosity). These results indicate that polypeptide collapse is kinetically a two-state transition. Because of the observed free energy barrier, the time scale of polypeptide collapse is dramatically slower than is predicted by Langevin models for homopolymer collapse.

Localized, stereochemically sensitive hydrophobic packing in an early folding intermediate of dihydrofolate reductase from Escherichia coli

Mutational analysis was performed to probe the development of hydrophobic clusters during the early events in the folding of dihydrofolate reductase. Replacements were made in several hydrophobic subdomains to examine the roles of hydrophobicity and stereochemistry in the formation of kinetic intermediates. Amide protons in two of these clusters, including residues I91, I94, and I155, have been shown to be protected against solvent exchange within 13 ms of folding. Additional hydrophobic clusters were probed by substitutions at residues I2, I61, and L112; these residues are not protected from exchange until later in the folding reaction. Valine and leucine replacements at positions I91, I94, and I155 significantly diminish the formation of the burst phase kinetic intermediate, relative to the wild-type protein. In contrast, I2 and I61 are insensitive to these substitutions in the first 5 ms of the folding reaction, as is the replacement of L112 with either isoleucine or valine. These results demonstrate that the tightly packed core of dihydrofolate reductase is acquired in a non-uniform fashion, beginning in the submillisecond time frame. The progressive development of specific side-chain packing in localized hydrophobic clusters may be a common theme for complex protein folding reactions.

Fast kinetics and mechanisms in protein folding

This review describes how kinetic experiments using techniques with dramatically improved time resolution have contributed to understanding mechanisms in protein folding. Optical triggering with nanosecond laser pulses has made it possible to study the fastest-folding proteins as well as fundamental processes in folding for the first time. These include formation of alpha-helices, beta-sheets, and contacts between residues distant in sequence, as well as overall collapse of the polypeptide chain. Improvements in the time resolution of mixing experiments and the use of dynamic nuclear magnetic resonance methods have also allowed kinetic studies of proteins that fold too fast (greater than approximately 10(3) s-1) to be observed by conventional methods. Simple statistical mechanical models have been extremely useful in interpreting the experimental results. One of the surprises is that models originally developed for explaining the fast kinetics of secondary structure formation in isolated peptides are also successful in calculating folding rates of single domain proteins from their native three-dimensional structure.

Fast events in protein folding: structural volume changes accompanying the early events in the N-->I transition of apomyoglobin induced by ultrafast pH jump

Ultrafast, laser-induced pH jump with time-resolved photoacoustic detection has been used to investigate the early protonation steps leading to the formation of the compact acid intermediate (I) of apomyoglobin (ApoMb). When ApoMb is in its native state (N) at pH 7.0, rapid acidification induced by a laser pulse leads to two parallel protonation processes. One reaction can be attributed to the binding of protons to the imidazole rings of His24 and His119. Reaction with imidazole leads to an unusually large contraction of -82 +/- 3 ml/mol, an enthalpy change of 8 +/- 1 kcal/mol, and an apparent bimolecular rate constant of (0.77 +/- 0.03) x 10(10) M(-1) s(-1). Our experiments evidence a rate-limiting step for this process at high ApoMb concentrations, characterized by a value of (0. 60 +/- 0.07) x 10(6) s(-1). The second protonation reaction at pH 7. 0 can be attributed to neutralization of carboxylate groups and is accompanied by an apparent expansion of 3.4 +/- 0.2 ml/mol, occurring with an apparent bimolecular rate constant of (1.25 +/- 0.02) x 10(11) M(-1) s(-1), and a reaction enthalpy of about 2 kcal/mol. The activation energy for the processes associated with the protonation of His24 and His119 is 16.2 +/- 0.9 kcal/mol, whereas that for the neutralization of carboxylates is 9.2 +/- 0.9 kcal/mol. At pH 4.5 ApoMb is in a partially unfolded state (I) and rapid acidification experiments evidence only the process assigned to carboxylate protonation. The unusually large contraction and the high energetic barrier observed at pH 7.0 for the protonation of the His residues suggests that the formation of the compact acid intermediate involves a rate-limiting step after protonation.

THEFASTPROTEINFOLDINGPROBLEM

During protein folding, many of the events leading to secondary and tertiary structure occur in milliseconds or faster. Modern nuclear magnetic resonance and laser detection techniques, coupled with fast initiation of the folding reaction, are probing these events in great detail. Theory, ranging from analytical models to molecular dynamics calculations, is beginning to match up with experiment. As a result, timescales, from such elementary steps as the addition of a residue to a helix to strange kinetics of collapsing protein backbones, can now be measured and interpreted.

Submillisecond unfolding kinetics of apomyoglobin and its pH 4 intermediate

Submillisecond mixing experiments and tryptophan fluorescence spectroscopy are used to address two questions raised in earlier stopped-flow studies of the folding and unfolding kinetics of sperm whale apomyoglobin. A study of the pH 4 folding intermediate (I) revealed, surprisingly, that its folding and unfolding kinetics are measurable and fit the two-state model except for a possible burst phase in unfolding. Submillisecond mixing experiments confirm the unfolding burst phase and show that its properties are consistent with the recently discovered interconversion between two forms of I, Ia equilibrium Ib. In urea-induced unfolding, Ib is converted to Ia before Ia unfolds, and the unfolding kinetics of Ia fit the two-state model when the burst phase is assigned to Ib-->Ia. The second question is whether the Ia, Ib intermediates accumulate transiently when the native protein (N) unfolds to the acid unfolded form (U). Earlier work showed that Ia and Ib accumulate when U refolds to N at pH 6.0 and the results fit the linear folding pathway U equilibrium Ia equilibrium Ib equilibrium N. We report here that either or both Ia and Ib accumulate transiently when N unfolds to U at pH 2.7 and that the position of the rate-limiting step in the pathway changes between unfolding at pH 2. 7 and refolding at pH 6.0. In unfolding as in refolding, we do not detect a fast track that bypasses the Ia, Ib intermediates.

Polymer principles and protein folding

This paper surveys the emerging role of statistical mechanics and polymer theory in protein folding. In the polymer perspective, the folding code is more a solvation code than a code of local phipsi propensities. The polymer perspective resolves two classic puzzles: (1) the Blind Watchmaker's Paradox that biological proteins could not have originated from random sequences, and (2) Levinthal's Paradox that the folded state of a protein cannot be found by random search. Both paradoxes are traditionally framed in terms of random unguided searches through vast spaces, and vastness is equated with impossibility. But both processes are partly guided. The searches are more akin to balls rolling down funnels than balls rolling aimlessly on flat surfaces. In both cases, the vastness of the search is largely irrelevant to the search time and success. These ideas are captured by energy and fitness landscapes. Energy landscapes give a language for bridging between microscopics and macroscopics, for relating folding kinetics to equilibrium fluctuations, and for developing new and faster computational search strategies.

The compact and expanded denatured conformations of apomyoglobin in the methanol-water solvent

We have performed a detailed study of methanol-induced conformational transitions of horse heart apomyoglobin (apoMb) to investigate the existence of the compact and expanded denatured states. A combination of far- and near-ultraviolet circular dichroism, NMR spectroscopy, and small-angle X-ray scattering (SAXS) was used, allowing a phase diagram to be constructed as a function of pH and the methanol concentration. The phase diagram contains four conformational states, the native (N), acid-denatured (U(A)), compact denatured (I(M)), and expanded helical denatured (H) states, and indicates that the compact denatured state (I(M)) is stable under relatively mild denaturing conditions, whereas the expanded denatured states (U(A) and H) are realized under extreme conditions of pH (strong electric repulsion) or alcohol concentration (weak hydrophobic interaction). The results of this study, together with many previous studies in the literature, indicate the general existence of the compact denatured states not only in the salt-pH plane but also in the alcohol-pH plane. Furthermore, to determine the general feature of the H conformation we used several proteins including ubiquitin, ribonuclease A, alpha-lactalbumin, beta-lactoglobulin, and Streptomyces subtilisin inhibitor (SSI) in addition to apoMb. SAXS studies of these proteins in 60% methanol showed that the H states of these all proteins have expanded and nonglobular conformations. The qualitative agreement of the experimental data with computer-simulated Kratky profiles also supports this structural feature of the H state.

Changes in side chain packing during apomyoglobin folding characterized by pulsed thiol-disulfide exchange

It is clear that close-packed side chain interactions play a dominant role in stabilizing native proteins, but the extent to which they stabilize kinetic intermediates and shape the energetic landscape of folding is not known. A method for characterizing structural changes at the level of individual side chains is presented and applied to study the refolding of apomyoglobin mutants containing engineered cysteine residues at key helical packing interfaces. The formation of buried side chain structure at the probe sites is followed by the extent of thiol-disulfide exchange during a pulse of thiol labeling reagent (either methyl methanethiosulfonate or 5,5'-dithiobis (2-nitrobenzoic acid)) applied at various stages of folding. The results suggest that the eight helices pack in at least three distinct stages, involving formation of two intermediates with time constants of <2 ms and 50 ms. In some parts of the refolding protein, stable side chain structure can be attained very rapidly, possibly in advance of backbone hydrogen bond formation as detected by previous pulsed amide hydrogen exchange experiments.