Folding of cytochrome c initiated by submillisecond mixing

Design and evaluation of a passive alcove-based microfluidic mixer

A novel passive microfluidic silicon mixer has been designed, optimized and fabricated. The architecture of the mixer consists of a simple "T" junction, made up by a 20 microm wide by 82 microm deep channel, followed by three repeats of an alcove, each with a triangular obstruction, arranged in a zigzag fashion. Numerical simulations were employed to optimize the geometry, particularly the dimensions of the alcoves, the relative orientation and the spacing between them, and the degree of intrusion associated with them. The simulation results demonstrate that chaotic flow due to recirculation within the alcoves results in transverse velocity that promotes effective fluid mixing. The microfluidic mixer with the simulation-optimized geometry was fabricated with photolithographic techniques and characterized by optical imaging, fluorescence, and Raman microscope spectroscopy. At a sample flow rate of 20 microL/s, the mixer exhibits a short mixing deadtime of approximately 22 micros and a high mixing efficiency under both low and high viscosity conditions. The alcove-based microfluidic silicon mixer offers unique advantages for its short deadtime and slow sample consumption rate. In addition, it provides a valuable component for laboratory-on-a-chip applications for its ease of development into multiple networks for massively parallel analytical processes.

Stability of the heme Fe-N-terminal amino group coordination bond in denatured cytochrome c

In the denatured states of Hydrogenobacter thermophilus cytochrome c(552) (HT) and Pseudomonas aeruginosa cytochrome c(551) (PA), and their mutants, the N-terminal amino group of the polypeptide chain is coordinated to heme Fe in place of the axial Met, the His-N(term) form being formed. The coordination of the N-terminal amino group to heme Fe leads to loop formation by the N-terminal stretch preceding the first Cys residue bound to the heme, and the N-terminal stretches of HT and PA are different from each other in terms of both the sequence and the number of constituent amino acid residues. The His-N(term) form was shown to be rather stable, and hence it can influence the stability of the denatured state. We have investigated the heme Fe coordination structures and stabilities of the His-N(term) forms emerging upon guanidine hydrochloric acid-induced unfolding of the oxidized forms of the proteins. The Fe-N(term) coordination bond in the His-N(term) form with a 9-residue N-terminal stretch of HT proteins was found to be tilted to some extent away from the heme normal, as reflected by the great heme methyl proton shift spread. On the other hand, the small heme methyl proton shift spread of the His-N(term) form with an 11-residue stretch of PA proteins indicated that its Fe-N(term) bond is nearly parallel with the heme normal. The stability of the His-N(term) form was found to be affected by the structural properties of the N-terminal stretch, such as its length and the N-terminal residue. With a given N-terminal residue, the stability of the His-N(term) form is higher for a 9-residue N-terminal stretch than an 11-residue one. In addition, with a given length of the N-terminal stretch, the His-N(term) form with an N-terminal Glu is stabilized by a few kJ mol(-1) relative to that with an N-terminal Asn. These results provide a novel insight into the stabilizing interactions in the denatured cyts c that will facilitate elucidation of the folding/unfolding mechanisms of the proteins.

Three-pulse photon echo peak shift spectroscopy as a probe of flexibility and conformational heterogeneity in protein folding

We investigate the equilibrium unfolding of Zn-cytochrome c in guanidine hydrochloride by three-pulse photon echo peak shift (3PEPS) spectroscopy. Unexpectedly, the measurements reveal that inhomogeneous broadening of the sample at the midpoint of the denaturation is larger than that of either native or unfolded states. To interpret this finding, we present simulations of the peak shift for both two-state and three-state unfolding models. Both the denaturant concentration dependence of the asymptotic peak shift (APS) and the wavelength dependence of the APS at the midpoint of the denaturation are different for the two models. Our data are consistent with two-state unfolding.

Intramolecular vibrational excitation of unfolding reactions in ZnII-substituted and metal-free cytochromes c: activation enthalpies from integrated fluorescence stokes shift and line shape excitation profiles

We have employed continuous-wave fluorescence spectroscopy to observe the light-induced formation of partially unfolded states of Zn(II)-substituted and metal-free (or free-base) cytochrome c (ZnCytc and fbCytc, respectively). In these experiments, the intrinsic porphyrin chromophore provides a vibrational excitation to the protein structure via intramolecular vibrational redistribution of the excess vibronic energy above the first excited singlet state. As the excitation light source is tuned, the fluorescence spectrum of both systems exhibits steplike transitions of the integrated Stokes shift, vibronic structure, and line width that mark apparent activation enthalpy barriers for structural transitions of the protein from the native state to a set of at least three partially unfolded states. The vibronic structure of the ZnCytc spectrum reports the exchange of the Zn(II) ion's native H18 and M80 axial ligands with non-native ligands as the excitation wavenumber is scanned through the three barriers. The metal ion's axial ligands contribute substantially to the stability of ZnCytc; the activation enthalpies for the corresponding transitions in fbCytc are one-third of those in ZnCytc. A comparison of the present results from ZnCytc with those obtained previously with picosecond time-resolved methods [Lampa-Pastirk and Beck, J. Phys. Chem. B 2006, 110, 22971-22974] indicates that the vibrationally excited protein structure propagates along an unfolding pathway from the native state that specifically populates the three states in order of their activation enthalpies. The excitation-wavenumber profile of the fluorescence line width is markedly inconsistent with a Maxwell-Boltzmann distribution over the three states. These results contrast with the general expectation of the protein-folding funnel hypothesis that a distribution of intermediate structures should result from the diffusive propagation of a nonequilibrium protein structure.

Reliable protein folding on complex energy landscapes: the free energy reaction path

A theoretical framework is developed to study the dynamics of protein folding. The key insight is that the search for the native protein conformation is influenced by the rate r at which external parameters, such as temperature, chemical denaturant, or pH, are adjusted to induce folding. A theory based on this insight predicts that 1), proteins with complex energy landscapes can fold reliably to their native state; 2), reliable folding can occur as an equilibrium or out-of-equilibrium process; and 3), reliable folding only occurs when the rate r is below a limiting value, which can be calculated from measurements of the free energy. We test these predictions against numerical simulations of model proteins with a single energy scale.

Hierarchical folding mechanism of apomyoglobin revealed by ultra-fast H/D exchange coupled with 2D NMR

The earliest steps in the folding of proteins are complete on an extremely rapid time scale that is difficult to access experimentally. We have used rapid-mixing quench-flow methods to extend the time resolution of folding studies on apomyoglobin and elucidate the structural and dynamic features of members of the ensemble of intermediate states that are populated on a submillisecond time scale during this process. The picture that emerges is of a continuum of rapidly interconverting states. Even after only 0.4 ms of refolding time a compact state is formed that contains major parts of the A, G, and H helices, which are sufficiently well folded to protect amides from exchange. The B, C, and E helix regions fold more slowly and fluctuate rapidly between open and closed states as they search docking sites on this core; the secondary structure in these regions becomes stabilized as the refolding time is increased from 0.4 to 6 ms. No further stabilization occurs in the A, G, H core at 6 ms of folding time. These studies begin to time-resolve a progression of compact states between the fully unfolded and native folded states and confirm the presence an ensemble of intermediates that interconvert in a hierarchical sequence as the protein searches conformational space on its folding trajectory.

Native and unfolded cytochrome c--comparison of dynamics using 2D-IR vibrational echo spectroscopy

Unfolded vs native CO-coordinated horse heart cytochrome c (h-cyt c) and a heme axial methionine mutant cyt c552 from Hydrogenobacter thermophilus ( Ht-M61A) are studied by IR absorption spectroscopy and ultrafast 2D-IR vibrational echo spectroscopy of the CO stretching mode. The unfolding is induced by guanidinium hydrochloride (GuHCl). The CO IR absorption spectra for both h-cyt c and Ht-M61A shift to the red as the GuHCl concentration is increased through the concentration region over which unfolding occurs. The spectra for the unfolded state are substantially broader than the spectra for the native proteins. A plot of the CO peak position vs GuHCl concentration produces a sigmoidal curve that overlays the concentration-dependent circular dichroism (CD) data of the CO-coordinated forms of both Ht-M61A and h-cyt c within experimental error. The coincidence of the CO peak shift curve with the CD curves demonstrates that the CO vibrational frequency is sensitive to the structural changes induced by the denaturant. 2D-IR vibrational echo experiments are performed on native Ht-M61A and on the protein in low- and high-concentration GuHCl solutions. The 2D-IR vibrational echo is sensitive to the global protein structural dynamics on time scales from subpicosecond to greater than 100 ps through the change in the shape of the 2D spectrum with time (spectral diffusion). At the high GuHCl concentration (5.1 M), at which Ht-M61A is essentially fully denatured as judged by CD, a very large reduction in dynamics is observed compared to the native protein within the approximately 100 ps time window of the experiment. The results suggest the denatured protein may be in a glassy-like state involving hydrophobic collapse around the heme.

Applications of terahertz spectroscopy in biosystems

Terahertz (THz) spectroscopic investigations of condensed-phase biological samples are reviewed ranging from the simple crystalline forms of amino acids, carbohydrates and polypeptides to the more complex aqueous forms of small proteins, DNA and RNA. Vibrationally resolved studies of crystalline samples have revealed the exquisite sensitivity of THz modes to crystalline order, temperature, conformational form, peptide sequence and local solvate environment and have given unprecedented measures of the binding force constants and anharmonic character of the force fields, properties necessary to improve predictability but not readily obtainable using any other method. These studies have provided benchmark vibrational data on extended periodic structures for direct comparisons with classical (CHARMm) and quantum chemical (density functional theory) theories. For the larger amorphous and/or aqueous phase samples, the THz modes form a continuum-like absorption that arises because of the full accessibility to conformational space and/or the rapid time scale for inter-conversion in these environments. Despite severe absorption by liquid water, detailed investigations have uncovered the photo- and hydration-induced conformational flexibility of proteins, the solvent shell depth of the water/biomolecule boundary layers and the solvent reorientation dynamics occurring in these interfacial layers that occur on sub-picosecond time scales. As such, THz spectroscopy has enhanced and extended the accessibility to intermolecular forces, length- and timescales important in biological structure and activity.

Characterization of N-terminal amino group–heme ligation emerging upon guanidine hydrochloric acid induced unfolding of Hydrogenobacter thermophilus ferricytochrome c 552

Nonnative heme coordination structures emerging upon guanidine hydrochloric acid (GdnHCl) induced unfolding of Hydrogenobacter thermophilus ferricytochrome c552 were characterized by means of paramagnetic NMR. The heme coordination structure possessing the N-terminal amino group of the peptide chain in place of axial Met (His-Nterm form) was determined in the presence of GdnHCl concentrations in excess of 1.5 M at neutral pH. The stability of the His-Nterm form at pH 7.0 was found to be comparable with that of the bis-His form which has been recognized as a major nonnative heme coordination structure in cytochrome c folding/unfolding. Consequently, in addition to the bis-His form, the His-Nterm form is a substantial intermediate which affects the pathway and kinetics of the folding/unfolding of cytochromes c, of which the N-terminal amino groups are not acetylated.

Folding kinetics of staphylococcal nuclease studied by tryptophan engineering and rapid mixing methods

To monitor the development of tertiary structural contacts during folding, a unique tryptophan residue was introduced at seven partially buried locations (residues 15, 27, 61, 76, 91, 102 and 121) of a tryptophan-free variant of staphylococcal nuclease (P47G/P117G/H124L/W140H). Thermal unfolding measurements by circular dichroism indicate that the variants are destabilized, but maintain the ability to fold into a native-like structure. For the variants with Trp at positions 15, 27 and 61, the intrinsic fluorescence is significantly quenched in the native state due to close contact with polar side-chains that act as intramolecular quenchers. All other variants exhibit enhanced fluorescence under native conditions consistent with burial of the tryptophan residues in an apolar environment. The kinetics of folding was observed by continuous and stopped-flow fluorescence measurements over refolding times ranging from 100 micros to 10 s. The folding kinetics of all variants is quantitatively described by a mechanism involving a major pathway with a series of intermediate states and a minor parallel channel. The engineered tryptophan residues in the beta-barrel and the N-terminal part of the alpha-helical domain become partially shielded from the solvent at an early stage (<1 ms), indicating that this region undergoes a rapid collapse. For some variants, a major increase in fluorescence coincides with the rate-limiting step of folding on the 100 ms time scale, indicating that these tryptophan residues are buried only during the late stages of folding. Other variants exhibit a transient increase in fluorescence during the 10 ms phase followed by a decrease during the rate-limiting phase. These observations are consistent with burial of these probes in a collapsed, but loosely packed intermediate, followed by the rate-limiting formation of the densely packed native core, which brings the tryptophan residues into close contact with intramolecular quenchers.

Heme coordination states of unfolded ferrous cytochrome C

The structural changes of ferrous Cyt-c that are induced by binding to SDS micelles, phospholipid vesicles, DeTAB, and GuHCl as well as by high temperatures and changes in the pH have been studied by RR and UV-Vis absorption spectroscopies. Four species have been identified in which the native methionine-80 ligand is removed from the heme iron. This coordination site is either occupied by a histidine (His-33 or His-26) to form a 6cLS configuration, which is the prevailing species in GuHCl at pH 7.0 and ambient temperature, or remains vacant to yield a 5cHS configuration. The three identified 5cHS species differ with respect to the hydrogen-bond interactions of the proximal histidine ligand (His-18) and include a nonhydrogen-bonded, a hydrogen-bonded, and a deprotonated imidazole ring. These structural motifs have been found irrespective of the unfolding conditions used. An unambiguous spectroscopic distinction of these 5cHS species is possible on the basis of the Fe-N(imidazole) stretching vibrations, the RR bands in the region between 1300 and 1650 cm(-1), and the electronic transitions in the Soret- and Q-band regions. In acid and neutral solutions, the species with a hydrogen-bonded and a nonhydrogen-bonded His-18 prevail, whereas in alkaline solutions a configuration with a deprotonated His-18 ligand is also observed. Upon lowering the pH or increasing the temperature in GuHCl solutions, the structure on the proximal side of the heme is perturbed, resulting in a loss of the hydrogen-bond interactions of the His-18 ligand. Conversely, the hydrogen-bonded His-18 of ferrous Cyt-c is stabilized by electrostatic interactions which increase in strength from phospholipid vesicles to SDS micelles. The results here suggest that unfolding of Cyt-c is initiated by the rupture of the Fe-Met-80 bond and structural reorganizations on the distal side of the heme pocket, whereas the proximal part is only affected in a later stage of the denaturation process.

Temperature-jump apparatus with Raman detection based on a solid-state tunable (1.80-2.05 microm) kHz optical parametric oscillator laser

The operating characteristics of a pulsed (10 ns) tunable near-infrared (NIR) laser source are described for temperature-jump (T-jump) applications. A Q-switched Nd:YLF laser (approximately 10 ns pulses) with a 1 kHz repetition rate is used to pump a potassium titanyl arsenate (KTA) crystal-based optical parametric oscillator (OPO), producing approximately 1 mJ NIR pulses that are tunable (1.80-2.05 microm) across the 1.9 microm vibrational overtone band of water. This T-jump source has been coupled to a deep ultraviolet (UV) probe laser for Raman studies of protein dynamics. T-jumps of up to 30 degrees C, as measured via the O-H stretching Raman band of water, are readily achieved. Application to cytochrome c unfolding is demonstrated.

A Weak Fe–O Bond in the Oxygenated Complex of the Nitric-oxide Synthase of Staphylococcus aureus

Little is known about the intermediates formed during catalysis by nitric-oxide synthase (NOS). We report here the characterization by resonance Raman spectroscopy of the oxygenated complex of the NOS from Staphylococcus aureus (saNOS) as well as the kinetics of formation and decay of the complex. An oxygenated complex transiently formed after mixing reduced saNOS with oxygen and decayed to the ferric enzyme with kinetics that were dependent on the substrate L-arginine and the cofactor H(4)B. The oxygenated complex displayed a Soret absorption band centered at 430 nm. Resonance Raman spectroscopy revealed that it can be described as a ferric superoxide form (Fe(III)O(2)(-)) with a single nu(O-O) mode at 1135 cm(-1). In the presence of L-arginine, an additional nu(O-O) mode at 1123 cm(-1) was observed, indicating an increased pi back-bonding electron donation to the bound oxygen induced by the substrate. With saNOS, this is the first time that the nu(Fe-O) mode of a NOS has been observed. The low frequency of this mode, at 517 cm(-1), points to an oxygenated complex that differs from that of P450(cam). The electronic structure of the oxygenated complex and the effect of L-arginine are discussed in relation to the kinetic properties of saNOS and other NOS.

Time-resolved Small-angle X-ray Scattering Investigation of the Folding Dynamics of Heme Oxygenase: Implication of the Scaling Relationship for the Submillisecond Intermediates of Protein Folding

Polypeptide collapse is generally observed as the initial folding dynamics of proteins with more than 100 residues, and is suggested to be caused by the coil-globule transition explained by Flory's theory of polymers. To support the suggestion by establishing a scaling behavior between radius of gyration (Rg) and chain length for the initial folding intermediates, the folding dynamics of heme oxygenase (HO) was characterized by time-resolved, small-angle X-ray scattering. HO is a highly helical protein without disulfide bridges, and is the largest protein (263 residues) characterized by the method. The folding process of HO was found to contain a transient oligomerization; however, the conformation within 10 ms was demonstrated to be monomeric and to possess Rg of 26.1(+/-1.1) A. Together with the corresponding data for proteins with different chain lengths, the seven Rg values demonstrated the scaling relationship to chain length with a scaling exponent of 0.35+/-0.11, which is close to the theoretical value of 1/3 predicted for globules in solutions where monomer-monomer interactions are favored over monomer-solvent interactions (poor solvent). The finding indicated that the initial folding dynamics of proteins bears the signature of the coil-globule transition, and offers a clue to explain the folding mechanisms of proteins with different chain lengths.

Structural Characterization of an Equilibrium Unfolding Intermediate in Cytochrome c

Although the denaturant-induced unfolding transition of cytochrome c was initially thought to be a cooperative process, recent spectroscopic studies have shown deviations from two-state behavior consistent with accumulation of an equilibrium intermediate. However, little is known about the structural and thermodynamic properties of this state, and whether it is stabilized by the presence of non-native heme ligands. We monitored the reversible denaturant-induced unfolding equilibrium of oxidized horse cytochrome c using various spectroscopic probes, including fluorescence, near and far-UV CD, heme absorbance bands in the Soret, visible and near-IR regions of the spectrum, as well as 2D NMR. Global fitting techniques were used for a quantitative interpretation of the results in terms of a three-state model, which enabled us to determine the intrinsic spectroscopic properties of the intermediate. A well-populated intermediate was observed in equilibrium experiments at pH 5 using either guanidine-HCl or urea as a denaturant, both for wild-type cytochrome c as well as an H33N mutant chosen to prevent formation of non-native His-heme ligation. For a more detailed structural characterization of the intermediate, we used 2D 1H-15N correlation spectroscopy to follow the changes in peak intensity for individual backbone amide groups. The equilibrium state observed in our optical and NMR studies contains many native-like structural features, including a well-structured alpha-helical sub-domain, a short Trp59-heme distance and solvent-shielded heme environment, but lacks the native Met80 sulfur-iron linkage and shows major perturbations in side-chain packing and other tertiary interactions. These structural properties are reminiscent of the A-state of cytochrome c, a compact denatured form found under acidic high-salt conditions, as well as a kinetic intermediate populated at a late stage of folding. The denaturant-induced intermediate also resembles alkaline forms of cytochrome c with altered heme ligation, suggesting that disruption of the native methionine ligand favors accumulation of structurally analogous states both in the presence and absence of non-native ligands.

Rapid binding of copper(I) to folded aporusticyanin

Kinetics of copper uptake in both oxidation states by the folded and unfolded forms of the type 1 copper protein rusticyanin have been studied. The speed of the binding of copper(I) to the folded rusticyanin is fast, and of the same order of magnitude as copper(I) uptake by the unfolded form. Thus, the binding of copper can be subsequent to the protein folding, contrary to previous proposals. Implications for the mechanism of the formation of the active holoprotein in vivo are discussed.

Effect of sol–gel encapsulation on the unfolding of ferric horse heart cytochrome c

Electronic absorption and resonance Raman spectra of ferric cytochrome c embedded in wet silica gels, in the presence of guanidine HCl as unfolding agent, between pH 0.35 and 7.0 are presented. The data clearly show that the ferric form of the protein encapsulated in sol-gel preserves its active site conformation. However, the spectra of the unfolded embedded protein are different from the corresponding spectra in solution suggesting that a strong interaction between the protein and the sol-gel takes place upon unfolding. The unfolding process mainly depends on the interaction between the exposed positive charges of the unfolded protein and the negatively charged functional groups of the silica surfaces. While this interaction partially stabilizes the protein in its native structure even at very acidic pH, in the presence of denaturants it has the opposite effect, causing mainly the weakening of both the heme-protein and the heme-ligand interactions.

Hydrogen bonding dynamics during protein folding of reduced cytochrome c: temperature and denaturant concentration dependence

Folding dynamics of reduced cytochrome c triggered by the laser-induced reduction method is investigated from a viewpoint of the intermolecular interaction change. Change of the diffusion coefficient of cytochrome c during the refolding process is traced in the time domain from the unfolded value to the native value continuously at various denaturant (guanidine hydrochloride (GdnHCl)) concentrations and temperatures. In the temperature range of 288 K-308 K and GdnHCl concentration range of 2.5 M-4.25 M, the diffusion change can be analyzed well by the two-state model consistently. It was found that the m(double dagger)-value and the activation energy of the transition state from the unfolded state for the hydrogen bonding network change are surprisingly similar to that for the local structural change around the heme group monitored by the fluorescence quenching experiment. This agreement suggests the existence of common or similar fundamental dynamics including water molecular movement to control the refolding dynamics. The nature of the transition state is discussed.

Specifically collapsed intermediate in the early stage of the folding of ribonuclease A

Nature of the burst-phase signals of protein folding has been the subject of much debate as to whether the signals represent the formation of early intermediates or the non-specific collapse of unfolded polypeptides. To distinguish the two possibilities, the submillisecond folding dynamics of ribonuclease A (RNase A) was examined, and compared with those of the disulfide bond-ruptured analog of RNase A (r-RNase A). The circular dichroism measurements on RNase A showed the burst-phase signal within 320 micros after the initiation of the folding reaction, which was identical to that observed for r-RNase A. In contrast, the burst phase increase in the extrinsic fluorescence from 1-anilino-8-naphthalene sulfonate (ANS) was observed for RNase A but not for r-RNase A. The kinetic titration experiment of the ANS fluorescence intensity showed the presence of a specific binding site for ANS in the fast-refolding component of RNase A. The small-angle X-ray scattering measurements at approximately 22 ms after initiating the folding reaction demonstrated that the burst phase conformations of the medium and slow-refolding components of RNase A were distinctly smaller than that of r-RNase A. These results indicated the difference in the burst phase conformations of RNase A and r-RNase A. Since r-RNase A is denatured in the physiological solution condition, the burst-phase signal of RNase A was interpreted as the formation of the folding intermediate with specific conformations.

A model for the misfolded bis-His intermediate of cytochrome c: the 1-56 N-fragment

We have characterized the ferric and ferrous forms of the heme-containing (1-56 residues) N-fragment of horse heart cytochrome c (cyt c) at different pH values and low ionic strength by UV-visible absorption and resonance Raman (RR) scattering. The results are compared with native cyt c in the same experimental conditions as this may provide a deeper insight into the cyt c unfolding-folding process. Folding of cyt c leads to a state having the heme iron coordinated to a histidine (His18) and a methionine (Met80) as axial ligands. At neutral pH the N-fragment (which lacks Met80) shows absorption and RR spectra that are consistent with the presence of a bis-His low spin heme, like several non-native forms of the parental protein. In particular, the optical spectra are identical to those of cyt c in the presence of a high concentration of denaturants; this renders the N-fragment a suitable model to study the heme pocket microenvironment of the misfolded (His-His) intermediate formed during folding of cyt c. Acid pH affects the ligation state in both cyt c and the N-fragment. Data obtained as a function of pH allow a correlation between the structural properties in the heme pocket of the N-fragment and those of non-native forms of cyt c. The results underline that the (57-104 residues) segment under native-like conditions imparts structural stability to the protein by impeding solvent access into the heme pocket.

Kinetics of intermolecular interaction during protein folding of reduced cytochrome c

Kinetics of intermolecular interaction between reduced cytochrome c (Cyt c) protein and solvent during the protein-refolding process is studied by monitoring the time dependence of apparent diffusion coefficient (D) using the pulsed-laser-induced transient grating technique. The refolding was triggered by photoinduced reduction of unfolded Fe(III) Cyt c in 3.5 M guanidine hydrochloride (GdnHCl) solution and the change in the diffusion coefficient was monitored in time domain. The relationship between D and the protein conformations under equilibrium condition were investigated at various GdnHCl concentrations using a photolabeling reagent. The time dependence of the observed transient grating signal was analyzed using these data and two models: a continuous change model of the intermolecular interaction and a two-state model. It was found that the TG signals in various time ranges can be consistently reproduced well by the two-state model. The dynamics of D is expressed well by a single exponential function with a rate constant of 22 +/- 7 s(-1) in a whole time range. The folding process of Cyt c is discussed based on these observations.

Specific collapse followed by slow hydrogen-bond formation of beta-sheet in the folding of single-chain monellin

Characterization of the conformational landscapes for proteins with different secondary structures is important in elucidating the mechanism of protein folding. The folding trajectory of single-chain monellin composed of a five-stranded beta-sheet and a helix was investigated by using a pH-jump from the alkaline unfolded to native state. The kinetic changes in the secondary structures and in the overall size and shape were measured by circular dichroism spectroscopy and small-angle x-ray scattering, respectively. The formation of the tertiary structure was monitored by intrinsic and extrinsic fluorescence. A significant collapse was observed within 300 micros after the pH-jump, leading to the intermediate with a small amount of secondary and tertiary structures but with an overall oblate shape. Subsequently, the stepwise formation of secondary and tertiary structures was detected. The current observation was consistent with the theoretical prediction that a more significant collapse precedes the formation of secondary structures in the folding of beta-sheet proteins than that of helical proteins [Shea, J. E., Onuchic, J. N. & Brooks, C. L., III (2002) Proc. Natl. Acad. Sci. USA 99, 16064-16068]. Furthermore, it was implied that the initial collapse was promoted by the formation of some specific structural elements, such as tight turns, to form the oblate shape.

The heme iron coordination of unfolded ferric and ferrous cytochrome c in neutral and acidic urea solutions. Spectroscopic and electrochemical studies

The heme iron coordination of unfolded ferric and ferrous cytochrome c in the presence of 7-9 M urea at different pH values has been probed by several spectroscopic techniques including magnetic and natural circular dichroism (CD), electrochemistry, UV-visible (UV-vis) absorption and resonance Raman (RR). In 7-9 M urea at neutral pH, ferric cytochrome c is found to be predominantly a low spin bis-His-ligated heme center. In acidic 9 M urea solutions the UV-vis and near-infrared (NIR) magnetic circular dichroism (MCD) measurements have for the first time revealed the formation of a high spin His/H(2)O complex. The pK(a) for the neutral to acidic conversion is 5.2. In 9 M urea, ferrous cytochrome c is shown to retain its native ligation structure at pH 7. Formation of a five-coordinate high spin complex in equilibrium with the native form of ferrous cytochrome c takes place below the pK(a) 4.8. The formal redox potential of the His/H(2)O complex of cytochrome c in 9 M urea at pH 3 was estimated to be -0.13 V, ca. 100 mV more positive than E degrees ' estimated for the bis-His complex of cytochrome c in urea solution at pH 7.

The 40s ?-loop plays a critical role in the stability and the alkaline conformational transition of cytochrome c

The structural and redox properties of a non-covalent complex reconstituted upon mixing two non-contiguous fragments of horse cytochrome c, the residues 1-38 heme-containing N-fragment with the residues 57-104 C-fragment, have been investigated. With respect to native cyt c, the complex lacks a segment of 18 residues, corresponding, in the native protein, to an omega (Omega)-loop region. The fragment complex shows compact structure, native-like alpha-helix content but a less rigid atomic packing and reduced stability with respect to the native protein. Structural heterogeneity is observed at pH 7.0, involving formation of an axially misligated low-spin species and consequent partial displacement of Met80 from the sixth coordination position of the heme-iron. Spectroscopic data suggest that a lysine (located in the Met80-containing loop, namely Lys72, Lys73, or Lys79) replaces the methionine residue. The residues 1-38/57-104 fragment complex shows an unusual biphasic alkaline titration characterized by a low (p K(a1)=6.72) and a high p K(a)-associated state transition (p K(a2)=8.56); this behavior differs from that of native cyt c, which shows a monophasic alkaline transition (p K(a)=8.9). The data indicate that the 40s Omega-loop plays an important role in the stability of cyt c and in ensuring a correct alkaline conformational transition of the protein.

Optical spectroscopic differentiation of various equilibrium denatured states of horse cytochrome c

Thermal unfolding of cytochrome c (cyt c) from several states has been studied using equilibrium spectroscopic techniques. CD in the uv, vibrational circular dichroism, infrared, and uv-vis absorption spectra measured at various temperatures, pHs, salt concentrations, and GuHCl concentrations are used to show the conformational as well as heme structural differences between native and various denatured states. The difference in thermal denaturation behaviors of cyt c starting from acid denatured, molten globule (MG), and the A and native states are explored. Different final high temperature states were observed for cytochrome c unfolding from four different initial states (native, MG, A, and acid denatured state) by electronic CD, Fourier transform infrared (FTIR), and vibrational CD (VCD). Consistent with this, different thermal unfolding pathways for the MG and A states are suggested by the FTIR and VCD data for this process.

Refolding processes of cytochrome P450cam from ferric and ferrous acid forms to the native conformation. Formations of folding intermediates with non-native heme coordination state

Changes in heme coordination state and protein conformation of cytochrome P450(cam) (P450(cam)), a b-type heme protein, were investigated by employing pH jump experiments coupled with time-resolved optical absorption, fluorescence, circular dichroism, and resonance Raman techniques. We found a partially unfolded form (acid form) of ferric P450(cam) at pH 2.5, in which a Cys(-)-heme coordination bond in the native conformation was ruptured. When the pH was raised to pH 7.5, the acid form refolded to the native conformation through a distinctive intermediate. Formations of similar acid and intermediate forms were also observed for ferrous P450(cam). Both the ferric and ferrous forms of the intermediate were found to have an unidentified axial ligand of the heme at the 6th coordination sphere, which is vacant in the high spin ferric and ferrous forms at the native conformation. For the ferrous form, it was also indicated that the 5th axial ligand is different from the native cysteinate. The folding intermediates identified in this study demonstrate occurrences of non-native coordination state of heme during the refolding processes of the large b-type heme protein, being akin to the well known folding intermediates of cytochromes c, in which c-type heme is covalently attached to a smaller protein.

Early Events During Folding of Wild-type Staphylococcal Nuclease and a Single-tryptophan Variant Studied by Ultrarapid Mixing

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

Crystal structure of the dioxygen-bound heme oxygenase from Corynebacterium diphtheriae: implications for heme oxygenase function

HmuO, a heme oxygenase of Corynebacterium diphtheriae, catalyzes degradation of heme using the same mechanism as the mammalian enzyme. The oxy form of HmuO, the precursor of the catalytically active ferric hydroperoxo species, has been characterized by ligand binding kinetics, resonance Raman spectroscopy, and x-ray crystallography. The oxygen association and dissociation rate constants are 5 microm(-1) s(-1) and 0.22 s(-1), respectively, yielding an O(2) affinity of 21 microm(-1), which is approximately 20 times greater than that of mammalian myoglobins. However, the affinity of HmuO for CO is only 3-4-fold greater than that for mammalian myoglobins, implying the presence of strong hydrogen bonding interactions in the distal pocket of HmuO that preferentially favor O(2) binding. Resonance Raman spectra show that the Fe-O(2) vibrations are tightly coupled to porphyrin vibrations, indicating the highly bent Fe-O-O geometry that is characteristic of the oxy forms of heme oxygenases. In the crystal structure of the oxy form the Fe-O-O angle is 110 degrees, the O-O bond is pointed toward the heme alpha-meso-carbon by direct steric interactions with Gly-135 and Gly-139, and hydrogen bonds occur between the bound O(2) and the amide nitrogen of Gly-139 and a distal pocket water molecule, which is a part of an extended hydrogen bonding network that provides the solvent protons required for oxygen activation. In addition, the O-O bond is orthogonal to the plane of the proximal imidazole side chain, which facilitates hydroxylation of the porphyrin alpha-meso-carbon by preventing premature O-O bond cleavage.

Hydrophobic effect of trityrosine on heme ligand exchange during folding of cytochrome c

Effect of a hydrophobic peptide on folding of oxidized cytochrome c (cyt c) is studied with trityrosine. Folding of cyt c was initiated by pH jump from 2.3 (acid-unfolded) to 4.2 (folded). The Soret band of the 2-ms transient absorption spectrum during folding decreased its intensity and red-shifted from 397 to 400 nm by interaction with trityrosine, whereas tyrosinol caused no significant effect. The change in the transient absorption spectrum by interaction with trityrosine was similar to that obtained with 100 mM imidazole, which showed that the population of the intermediate His/His coordinated species increased during folding of cyt c by interaction with trityrosine. The absorption change was biphasic, the fast phase (82+/-9s(-1)) corresponding to the transition from the His/H(2)O to the His/Met coordinated species, whereas the slow phase (24+/-3s(-1)) from His/His to His/Met. By addition of trityrosine, the relative ratio of the slow phase increased, due to increase of the His/His species at the initial stage of folding. According to the resonance Raman spectra of cyt c, the high-spin 6-coordinate and low-spin 6-coordinate species were dominated at pH 2.3 and 4.2, respectively, and these species were not affected by addition of trityrosine. These results demonstrated that the His/His species increased by interaction with trityrosine at the initial stage of cyt c folding, whereas the heme coordination structure was not affected by trityrosine when the protein was completely unfolded or folded. Hydrophobic peptides thus may be useful to study the effects of hydrophobic interactions on protein folding.

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.

Cytochrome c(551) as a model system for protein folding

Considerable progress was made over the last few years in understanding the mechanism of folding of cytochrome c(551), a small acidic hemeprotein from Pseudomonas aeruginosa. Comparison of our results with those obtained by others on horse heart cytochrome c allows to draw some general conclusions on the structural features that are common determinants in the folding of members of the cytochrome c family.

Exploring the cytochrome c folding mechanism: cytochrome c552 from thermus thermophilus folds through an on-pathway intermediate

Understanding the role of partially folded intermediate states in the folding mechanism of a protein is a crucial yet very difficult problem. We exploited a kinetic approach to demonstrate that a transient intermediate of a thermostable member of the widely studied cytochrome c family (cytochrome c552 from Thermus thermophilus) is indeed on-pathway. This is the first clear indication of an obligatory intermediate in the folding mechanism of a cytochrome c. The fluorescence properties of this intermediate demonstrate that the relative position of the heme and of the only tryptophan residue cannot correspond to their native orientation. Based on an analysis of the three-dimensional structure of cytochrome c552, we propose an interpretation of the data which explains the residual fluorescence of the intermediate and is consistent with the established role played by some conserved interhelical interactions in the folding of other members of this family. A limited set of topologically conserved contacts may guide the folding of evolutionary distant cytochromes c through the same partially structured state, which, however, can play different kinetic roles, acting either as an intermediate or a transition state.

A novel method for study of protein folding kinetics by monitoring diffusion coefficient in time domain

Molecular diffusion process after the photo-induced electron injection to ferric cytochrome c (Fe(III) cyt c) in guanidine hydrochloride (GdnHCl) 3.5 M buffer solution is studied by the time-resolved transient grating technique. Circular dichroism studies have revealed that Fe(III) cyt c is unfolded under this condition but the reduced form, Fe(II) cyt c, is folded. Hence, this pulsed laser-induced reduction should initiate the folding process of cyt c. The observed transient grating signal shows prominent features, which have never been observed before. Based on several characteristic points, we concluded that the apparent diffusion coefficient (D) of Fe(II) cyt c after the reduction is time dependent, which must be associated with the protein folding dynamics. This time-dependent apparent D should reflect either the continuous time development of the hydrodynamic radius or population change of the unfolded and folded states during the folding dynamics. This is the first observation of the time-dependent apparent D during any chemical reaction, and this time-dependent measurement of D should be a unique and powerful way to study the protein folding kinetics from a viewpoint of the protein's shape or the protein-water intermolecular interaction.

Parallel pathways in cytochrome c(551) folding

The folding of cytochrome c(551) from Pseudomonas aeruginosa was previously thought to follow a simple sequential mechanism, consistent with the lack of histidine residues, other than the native His16 heme ligand, that can give rise to mis-coordinated species. However, further kinetic analysis reveals complexities indicative of a folding mechanism involving parallel pathways. Double-jump interrupted refolding experiments at low pH indicate that approximately 50% of the unfolded cytochrome c(551) population can reach the native state via a fast (10 ms) folding track, while the rest follows a slower folding path with populated intermediates. Stopped-flow experiments using absorbance at 695 nm to monitor refolding confirm the presence of a rapidly folding species containing the native methionine-iron bond while measurements on carboxymethylated cytochrome c(551) (which lacks the Met-Fe coordination bond) indicate that methionine ligation occurs late during folding along the fast folding track, which appears to be dominant at physiological pH. Continuous-flow measurements of tryptophan-heme energy transfer, using a capillary mixer with a dead time of about 60 micros, show evidence for a rapid chain collapse within 100 micros preceding the rate-limiting folding phase on the milliseconds time scale. A third process with a time constant in the 10-50 ms time range is consistent with a minor population of molecules folding along a parallel channel, as confirmed by quantitative kinetic modeling. These findings indicate the presence of two or more slowly inter-converting ensembles of denatured states that give rise to pH-dependent partitioning among fast and slow-folding pathways.

Anion concentration modulates the conformation and stability of the molten globule of cytochrome c

Anions induce collapse of acid-denatured cytochrome c into a compact state, the A-state, showing molten globule character. Since structural information on partially folded forms of proteins is important for a deeper understanding of folding mechanisms and of the factors affecting protein stabilization, in this paper we have investigated in detail the effects of anions on the tertiary conformation of the A-state. We have found that the salt-induced collapse of acid-denatured cytochrome c leads to a number of equilibria between high-spin and low-spin heme states and between two types of low-spin states. The two latter states are characterized by conformations leading to a native-like Met-Fe-His axial coordination and a bis-His configuration. The equilibrium between these two A-states is dependent on the concentration and/or size of the anions (i.e. the bigger the anion, the greater its effect). Further, on the basis of fast kinetic data, a kinetic model of the folding process from the acid-unfolded protein to the A-state (at low and high anion concentration) is described.

Kinetics of cytochrome C folding: atomically detailed simulations

The vast range of time scales (from nanoseconds to seconds) during protein folding is a challenge for experiments and computations. To make concrete predictions on folding mechanisms, atomically detailed simulations of protein folding, using potentials derived from chemical physics principles, are desired. However, due to their computational complexity, straightforward molecular dynamics simulations of protein folding are impossible today. An alternative algorithm is used that makes it possible to compute approximate atomically detailed long time trajectories (the Stochastic Difference Equation in Length). This algorithm is used to compute 26 atomically detailed folding trajectories of cytochrome c (a millisecond process). The early collapse of the protein chain (with marginal formation of secondary structure), and the earlier formation of the N and C helices (compare to the 60's helix) are consistent with the experiment. The existence of an energy barrier upon entry to the molten globule is examined as well. In addition to (favorable) comparison to experiments, we show that non-native contacts drive the formation of the molten globule. In contrast to popular folding models, the non-native contacts do not form off-pathway kinetic traps in cytochrome c.

Activation of hydrogen peroxide in horseradish peroxidase occurs within approximately 200 micro s observed by a new freeze-quench device

To observe the formation process of compound I in horseradish peroxidase (HRP), we developed a new freeze-quench device with approximately 200 micro s of the mixing-to-freezing time interval and observed the reaction between HRP and hydrogen peroxide (H(2)O(2)). The developed device consists of a submillisecond solution mixer and rotating copper or silver plates cooled at 77 K; it freezes the small droplets of mixed solution on the surface of the rotating plates. The ultraviolet-visible spectra of the sample quenched at approximately 1 ms after the mixing of HRP and H(2)O(2) suggest the formation of compound I. The electron paramagnetic resonance spectra of the same reaction quenched at approximately 200 micro s show a convex peak at g = 2.00, which is identified as compound I due to its microwave power and temperature dependencies. The absence of ferric signals in the electron paramagnetic resonance spectra of the quenched sample indicates that compound I is formed within approximately 200 micro s after mixing HRP and H(2)O(2). We conclude that the activation of H(2)O(2) in HRP at ambient temperature completes within approximately 200 micro s. The developed device can be generally applied to investigate the electronic structures of short-lived intermediates of metalloenzymes.

Inter- and intra-molecular electron transfer in the cytochrome bc(1) complex

In this review, we compare the intra-molecular and inter-molecular electron transfer rate constants of the high-potential branch of the cytochrome bc(1) complex. Several methods such as the conventional stopped-flow spectroscopy, pH-induced electron transfer, photoactivated ruthenium complex induced electron transfer and photoreleaseable caged quinol, have been used to determine reaction rates between redox centers in an attempt to elucidate the reaction mechanism of this vital energy conserving complex. Since the most active pure cytochrome bc(1) complex has a turnover number of 800 s(-1), any step with a rate constant much larger than this will not be rate-limiting. The most likely rate-limiting step is the cytochrome b redox state governed movement of the head domain of iron-sulfur protein from its electron-accepting site ("fixed" or "b-state" position) to its electron donating site ("c(1)-state" position).

Protein folding in the post-genomic era

Protein folding is a topic of fundamental interest since it concerns the mechanisms by which the genetic message is translated into the three-dimensional and functional structure of proteins. In these post-genomic times, the knowledge of the fundamental principles are required in the exploitation of the information contained in the increasing number of sequenced genomes. Protein folding also has practical applications in the understanding of different pathologies and the development of novel therapeutics to prevent diseases associated with protein misfolding and aggregation. Significant advances have been made ranging from the Anfinsen postulate to the "new view" which describes the folding process in terms of an energy landscape. These new insights arise from both theoretical and experimental studies. The problem of folding in the cellular environment is briefly discussed. The modern view of misfolding and aggregation processes that are involved in several pathologies such as prion and Alzheimer diseases. Several approaches of structure prediction, which is a very active field of research, are described.

Reaction of hydrogen peroxide and peroxidase activity in carboxymethylated cytochrome c: spectroscopic and kinetic studies

The peroxidase activity of carboxymethylated cytochrome c (Cmcytc) has been investigated by spectroscopic and kinetic techniques to examine the effect of carboxymethylation on the peroxidase activity of native cytochrome c (cytc). The optical spectrum suggests that the reaction of Cmcytc with H(2)O(2) proceeds through only one intermediate, compound I. The apparent rate constant (k(app)) for the reaction was found to be 17, 72 and 210 M(-1) s(-1) at pH 7.0, 5.0 and 3.5 respectively. These values are about 60 times larger than those reported for native cytc (0.236 M(-1) s(-1) at pH 7.0), and about five orders of magnitude lower than those for classical peroxidases. Cmcytc was found to catalyse oxidation of organic and inorganic substrates. The second order rate constant for the oxidation of 2,2'-azino-bis(3-ethylbenzthiazoline-6-sulphonic acid) (ABTS) by Cmcytc (205 [H(2)O(2)] s(-1)) is found to be larger than the corresponding value for native cytc (50 [H(2)O(2)] s(-1)) at pH 6.0. The carboxymethylation of cytc ruptures the Fe-S (Met 80) bond and increases the rate of its reaction with H(2)O(2), and its catalytic activity. The specific activity of Cmcytc was measured spectrophotometrically by the reported method using ABTS as substrate, and was found to be 288, 473 and 872 microM min(-1) mg(-1) at pH 7.0, 5.0 and 3.5 respectively. Resonance Raman studies indicated the presence of a bis-histidine coordinated form of Cmcytc at neutral pH, and the existence of a population distribution of different ligation states such as bis-histidine (HH), histidine-water (HW) and five coordinate (5C) forms at lower pH. The relative population of different species in Cmcytc was found to be HH (approximately 100%, approximately 50%, approximately 44%), HW (approximately 0%, approximately 44%, 41%) and 5C (approximately 0%, approximately 6%, 15%) at pH 7.0, 4.7 and 3.1 respectively. We have attempted to correlate the pH dependence of the reaction of Cmcytc with hydrogen peroxide and its peroxidase activity with the haem stereochemical structures observed for Cmcytc. Steady-state and time-resolved tryptophan fluorescence studies on Cmcytc were done to probe the conformational changes around the haem pocket of Cmcytc.