Snapshots of protein folding. A study on the multiple transition state pathway of cytochrome c(551) from Pseudomonas aeruginosa

A short structural extension dictates the early stages of folding of a PDZ domain

PDZ domains are highly abundant protein-protein interaction modules in human. One of the most extensively characterized PDZ domain, the third PDZ domain from PSD-95 (PDZ3), contains an α-helical C-terminal extension that has a key role in the function of the domain. Here we compared the folding of PDZ3 with a truncated variant (PDZ3Δα3), lacking the additional helix, by means of the so-called Φ-value analysis, an experimental technique that allows inferring the structure of folding transition states. Experiments reveal subtle but detectable differences in the folding of PDZ3Δα3 versus PDZ3, as probed by structural characterization of the folding transition states. These differences appear more remarkable in the early stages of folding, with a detectable shift of the folding nucleus. The presented results allow demonstrating that the native state exerts a weak bias at the early stages of folding, which appear to be characterized by alternative pathways.

Studies of cytochrome c-551 unfolding using fluorescence correlation spectroscopy and other biophysical techniques

Fluorescence correlation spectroscopy studies with a bacterial cytochrome c labeled at different positions complement NMR hydrogen exchange results.

Domain-swapped dimer of Pseudomonas aeruginosa cytochrome c551: structural insights into domain swapping of cytochrome c family proteins

Cytochrome c (cyt c) family proteins, such as horse cyt c, Pseudomonas aeruginosa cytochrome c₅₅₁ (PA cyt c₅₅₁), and Hydrogenobacter thermophilus cytochrome c₅₅₂ (HT cyt c₅₅₂), have been used as model proteins to study the relationship between the protein structure and folding process. We have shown in the past that horse cyt c forms oligomers by domain swapping its C-terminal helix, perturbing the Met–heme coordination significantly compared to the monomer. HT cyt c₅₅₂ forms dimers by domain swapping the region containing the N-terminal α-helix and heme, where the heme axial His and Met ligands belong to different protomers. Herein, we show that PA cyt c₅₅₁ also forms domain-swapped dimers by swapping the region containing the N-terminal α-helix and heme. The secondary structures of the M61A mutant of PA cyt c₅₅₁ were perturbed slightly and its oligomer formation ability decreased compared to that of the wild-type protein, showing that the stability of the protein secondary structures is important for domain swapping. The hinge loop of domain swapping for cyt c family proteins corresponded to the unstable region specified by hydrogen exchange NMR measurements for the monomer, although the swapping region differed among proteins. These results show that the unstable loop region has a tendency to become a hinge loop in domain-swapped proteins.

The role of key residues in structure, function, and stability of cytochrome-c

Cytochrome-c (cyt-c), a multi-functional protein, plays a significant role in the electron transport chain, and thus is indispensable in the energy-production process. Besides being an important component in apoptosis, it detoxifies reactive oxygen species. Two hundred and eighty-five complete amino acid sequences of cyt-c from different species are known. Sequence analysis suggests that the number of amino acid residues in most mitochondrial cyts-c is in the range 104 ± 10, and amino acid residues at only few positions are highly conserved throughout evolution. These highly conserved residues are Cys14, Cys17, His18, Gly29, Pro30, Gly41, Asn52, Trp59, Tyr67, Leu68, Pro71, Pro76, Thr78, Met80, and Phe82. These are also known as "key residues", which contribute significantly to the structure, function, folding, and stability of cyt-c. The three-dimensional structure of cyt-c from ten eukaryotic species have been determined using X-ray diffraction studies. Structure analysis suggests that the tertiary structure of cyt-c is almost preserved along the evolutionary scale. Furthermore, residues of N/C-terminal helices Gly6, Phe10, Leu94, and Tyr97 interact with each other in a specific manner, forming an evolutionary conserved interface. To understand the role of evolutionary conserved residues on structure, stability, and function, numerous studies have been performed in which these residues were substituted with different amino acids. In these studies, structure deals with the effect of mutation on secondary and tertiary structure measured by spectroscopic techniques; stability deals with the effect of mutation on T m (midpoint of heat denaturation), ∆G D (Gibbs free energy change on denaturation) and folding; and function deals with the effect of mutation on electron transport, apoptosis, cell growth, and protein expression. In this review, we have compiled all these studies at one place. This compilation will be useful to biochemists and biophysicists interested in understanding the importance of conservation of certain residues throughout the evolution in preserving the structure, function, and stability in proteins.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

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.

Submolecular unfolding units of Pseudomonas aeruginosa cytochrome c-551

Hydrogen exchange rates for backbone amide protons of oxidized Pseudomonas aeruginosa cytochrome c-551 (P. aeruginosa cytochrome c) have been measured in the presence of low concentrations of the denaturant guanidine hydrochloride. Analysis of the data has allowed identification of submolecular unfolding units known as foldons. The highest-energy foldon bears similarity to the proposed folding intermediate for P. aeruginosa cytochrome c. Parallels are seen to the foldons of the structurally homologous horse cytochrome c, although the heme axial methionine-bearing loop has greater local stability in P. aeruginosa cytochrome c, in accord with previous folding studies. Regions of low local stability are observed to correspond with regions that interact with redox partners, providing a link between foldon properties and function.

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.

Plasticity of the protein folding landscape: switching between on- and off-pathway intermediates

Proteins may fold via parallel routes partitioned by the relative effect of solvent conditions on the relevant transition states. Thus, intermediates may or may not necessarily be obligatory species accumulated during the folding process, but rather kinetic traps due to the ruggedness of the folding landscape. Implicit in this view is the notion of plasticity of the folding pathway: proteins can be rerouted through the energy landscape by mutational, topological or solvent perturbations. Our work was specifically aimed to the experimental identification of a switch in the folding mechanism of a c-type cytochrome from the thermophilic bacterium Hydrogenobacter thermophilus (HT cyt c(552)) induced by acidic conditions. We present evidence that, by destabilizing the relevant transition state, the native state of HT cyt c(552) can be reached along alternative folding routes, which may involve an off-pathway intermediate.

Unveiling a hidden folding intermediate in c-type cytochromes by protein engineering

Several investigators have highlighted a correlation between the basic features of the folding process of a protein and its topology, which dictates the folding pathway. Within this conceptual framework we proposed that different members of the cytochrome c (cyt c) family share the same folding mechanism, involving a consensus partially structured state. Pseudomonas aeruginosa cyt c(551) (Pa cyt c(551)) folds via an apparent two-state mechanism through a high energy intermediate. Here we present kinetic evidence demonstrating that it is possible to switch its folding mechanism from two to three state, stabilizing the high energy intermediate by rational mutagenesis. Characterization of the folding kinetics of one single-site mutant of the Pa cyt c(551) (Phe(7) to Ala) indeed reveals an additional refolding phase and a fast unfolding process which are explained by the accumulation of a partially folded species. Further kinetic analysis highlights the presence of two parallel processes both leading to the native state, suggesting that the above mentioned species is a non obligatory on-pathway intermediate. Determination of the crystallographic structure of F7A shows the presence of an extended internal cavity, which hosts three "bound" water molecules and a H-bond in the N-terminal helix, which is shorter than in the wild type protein. These two features allow us to propose a detailed structural interpretation for the stabilization of the native and especially the intermediate states induced by a single crucial mutation. These results show how protein engineering, x-ray crystallography and state-of-the-art kinetics concur to unveil a folding intermediate and the structural determinants of its stability.

The transition state for folding of a peripheral subunit-binding domain contains robust and ionic-strength dependent characteristics

The denaturant dependencies of the folding and unfolding kinetics were used to characterize the structure of the transition state for folding of E3BD, a peripheral subunit-binding domain. For the majority of E3BD mutants, the Phi-values calculated at 298 K from the analysis of chevron plots were in good agreement with those previously determined at 325 K using Arrhenius analysis. This agreement further demonstrates the general robustness of Phi-value analyses, since different experiments, methods of denaturation and thermodynamic assumptions were used to determine each set of Phi(F) values. The structure of the transition state for folding was grossly conserved at 298 K and 325 K, with residues in Helix I playing a lesser role in folding than those located in the 3(10) helix, disordered loop and Helix II. However, the energetic contributions of a cluster of basic residues close to the N-terminus and Helix I, which are an integral part of the ligand-binding site, were susceptible to ionic strength effects because of electrostatic strain in native and transition states of E3BD at low ionic strength. We found no evidence of the downhill folding previously proposed for E3BD, even though the conditions employed in this study significantly increased the energetic bias towards the native state.

Kinetic folding mechanism of PDZ2 from PTP-BL

PDZ domains represent a large family of protein-interaction modules associated with a variety of unrelated proteins with different functions. We report a complete characterization of the kinetic folding mechanism of a fluorescent variant of PDZ2 from PTP-BL, investigated under a variety of different experimental conditions. For this purpose, we engineered a fluorescent variant of this protein Y43W (called pseudo-wild-type, pWT43). The results suggest the presence of a high-energy intermediate in the folding of PDZ2, as revealed by a pronounced non-linear dependence of the unfolding rate constant on denaturant concentration. Such an intermediate may or may not be detectable depending on the experimental conditions, giving rise to apparent two-state folding under stabilizing conditions (e.g. in the presence of sodium sulfate). Interestingly, even under these conditions, three-state folding can be restored by selectively destabilizing the native-like rate-limiting barrier by one specific mutation (V44A). Finally, we show that data taken on pWT43 under different experimental conditions (e.g. different pH values from 2.1 to 8.0 or in the presence of a stabilizing salt) and also data on a site-directed conservative mutant can be rationalized in terms of a simple reaction scheme involving a single set of intermediates and transition states.

An obligatory intermediate in the folding pathway of cytochrome c552 from Hydrogenobacter thermophilus

The folding mechanism of many proteins involves the population of partially organized structures en route to the native state. Identification and characterization of these intermediates is particularly difficult, as they are often only transiently populated and may play different mechanistic roles, being either on-pathway productive species or off-pathway kinetic traps. Following different spectroscopic probes, and employing state-of-the-art kinetic analysis, we present evidence that the folding mechanism of the thermostable cytochrome c552 from Hydrogenobacter thermophilus does involve the presence of an elusive, yet compact, on-pathway intermediate. Characterization of the folding mechanism of this cytochrome c is particularly interesting for the purpose of comparative folding studies, because H. thermophilus cytochrome c552 shares high sequence identity and structural homology with its homologue from the mesophilic bacterium Pseudomonas aeruginosa cytochrome c551, which refolds through a broad energy barrier without the accumulation of intermediates. Analysis of the folding kinetics and correlation with the three-dimensional structure add new evidence for the validity of a consensus folding mechanism in the cytochrome c family.

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.

Unifying features in protein-folding mechanisms

We compare the folding of representative members of a protein superfamily by experiment and simulation to investigate common features in folding mechanisms. The homeodomain superfamily of three-helical, single-domain proteins exhibits a spectrum of folding processes that spans the complete transition from concurrent secondary and tertiary structure formation (nucleation-condensation mechanism) to sequential secondary and tertiary formation (framework mechanism). The unifying factor in their mechanisms is that the transition state for (un)folding is expanded and very native-like, with the proportion and degree of formation of secondary and tertiary interactions varying. There is a transition, or slide, from the framework to nucleation-condensation mechanism with decreasing stability of the secondary structure. Thus, framework and nucleation-condensation are different manifestations of an underlying common mechanism.

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.

Hydrogen exchange in a bacterial cytochrome c: a fingerprint of the cytochrome c fold

The hydrogen exchange rates of backbone amides in a minimal (71 amino acid long) monoheme cytochrome c were determined as a function of pH in the absence and in the presence of guanidinium chloride. These data permitted the identification of units undergoing the opening reaction that precedes hydrogen exchange through a common mechanism. The opening units broadly correlate with the secondary structure elements of the protein. It is found that, despite the significant difference in primary sequence, the distribution of the opening units within the three-dimensional structure of the cytochrome studied here closely resembles that determined in mitochondrial c-type cytochromes. It is proposed that the observed distribution represents a fingerprint of the cytochrome c fold and has a role in directing the folding/unfolding of the protein.

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.

Peroxidase activity as a tool for studying the folding of c-type cytochromes

The peroxidase activity of c-type cytochromes increases substantially by unfolding. This phenomenon was used to study the equilibrium unfolding of ferricytochrome c. The peroxidase activity is already enhanced at low denaturant concentrations. The lowest free energy folding intermediate is easily detected by this method, while it is invisible using fluorescence or optical spectroscopy. The free energy difference between this folding intermediate and the native state depends on the strength of the sixth ligand of the heme-iron and the increase in peroxidase activity upon unfolding is shown to be a sensitive indicator of the strength of this ligand. Under fully denaturing conditions, the peroxidase activity is inhibited by protein-based ligands. It is shown that at least three different ligand groups can be responsible for this inhibition, and that at neutral or alkaline pH, the predominant ligand is not histidine. The use of peroxidase activity assays as a method to study the unfolding of cytochrome c is evaluated.

Refolding kinetics of cytochrome c(551) reveals a mechanistic difference between urea and guanidine

The energetic parameters for the folding of small globular proteins can be very different if derived from guanidine hydrochloride (GdnHCl) or urea denaturation experiments. A study of the equilibrium and kinetics of the refolding of wild-type (wt) cytochrome c(551) (cyt c(551)) from Pseudomonas aeruginosa and of two site-directed mutants (E70Q and E70V) shows that the nonionic nature of urea reveals the role of a salt bridge between residues E70 and K10 on the transition state, which is otherwise completely masked in GdnHCl experiments. Mixed denaturant refolding experiments allow us to conclude that the masking effect of GdnHCl is complete at fairly low GdnHCl concentrations ( congruent with 0.1 M). The fact that potassium chloride is unable to reproduce this quenching effect, together with the results obtained on the mutants, suggests a specific binding of the Gdn(+) cation, which involves the E70-K10 ion pair in wt cyt c(551). We propose, therefore, a simple kinetic test to obtain a mechanistic interpretation of nonlinear dependences of DeltaG(w) on GdnHCl concentration on the basis of kinetic refolding experiments in the presence of both denaturants.