A lysine 73-->histidine variant of yeast iso-1-cytochrome c: evidence for a native-like intermediate in the unfolding pathway and implications for m value effects

Effect of proline content and histidine ligation on the dynamics of Ω-loop D and the peroxidase activity of iso-1-cytochrome c

To study how proline residues affect the dynamics of Ω-loop D (residues 70 to 85) of cytochrome c, we prepared G83P and G83A variants of yeast iso-1-cytochrome c (iso-1-Cytc) in the presence and absence of a K73H mutation. Ω-loop D is important in controlling both the electron transfer function of Cytc and the peroxidase activity of Cytc used in apoptosis because it provides the Met80 heme ligand. The G83P and G83A mutations have no effect on the global stability of iso-1-Cytc in presence or absence of the K73H mutation. However, both mutations destabilize the His73-mediated alkaline conformer relative to the native state. pH jump stopped-flow experiments show that the dynamics of the His73-mediated alkaline transition are significantly enhanced by the G83P mutation. Gated electron transfer studies show that the enhanced dynamics result from an increased rate of return to the native state, whereas the rate of loss of Met80 ligation is unchanged by the G83P mutation. Thus, the G83P substitution does not stiffen the conformation of the native state. Because bis-His heme ligation occurs when Cytc binds to cardiolipin-containing membranes, we studied the effect of His73 ligation on the peroxidase activity of Cytc, which acts as an early signal in apoptosis by causing oxygenation of cardiolipin. We find that the His73 alkaline conformer suppresses the peroxidase activity of Cytc. Thus, the bis-His ligated state of Cytc formed upon binding to cardiolipin is a negative effector for the peroxidase activity of Cytc early in apoptosis.

How to Turn an Electron Transfer Protein into a Redox Enzyme for Biosensing

Cytochrome c is a small globular protein whose main physiological role is to shuttle electrons within the mitochondrial electron transport chain. This protein has been widely investigated, especially as a paradigmatic system for understanding the fundamental aspects of biological electron transfer and protein folding. Nevertheless, cytochrome c can also be endowed with a non-native catalytic activity and be immobilized on an electrode surface for the development of third generation biosensors. Here, an overview is offered of the most significant examples of such a functional transformation, carried out by either point mutation(s) or controlled unfolding. The latter can be induced chemically or upon protein immobilization on hydrophobic self-assembled monolayers. We critically discuss the potential held by these systems as core constituents of amperometric biosensors, along with the issues that need to be addressed to optimize their applicability and response.

Met80 and Tyr67 affect the chemical unfolding of yeast cytochromec: comparing the solutionvs.immobilized state

The motional regime affects the unfolding propensity and axial heme coordination of the Met80Ala and Met80Ala/Tyr67Ala variants of yeast iso-1 cytochromec.

Naturally Occurring A51V Variant of Human Cytochrome c Destabilizes the Native State and Enhances Peroxidase Activity

The A51V variant of human cytochrome c is linked to thrombocytopenia 4 (THC4), a condition that causes decreased blood platelet counts. A 1.82 Å structure of the A51V variant shows only minor changes in tertiary structure relative to the wild-type (WT) protein. Guanidine hydrochloride denaturation demonstrates that the global stability of the A51V variant is 1.3 kcal/mol less than that of the WT protein. The midpoint pH, pH_1/2, of the alkaline transition of the A51V variant is 1 unit less than that of the WT protein. Stopped-flow pH jump experiments show that the A51V substitution affects the triggering ionization for one of two kinetically distinguishable alkaline conformers and enhances the accessibility of a high-spin heme transient. The pH_1/2 for acid unfolding of the A51V variant is 0.7 units higher than for that of the WT protein. Consistent with the greater accessibility of non-native conformers for the A51V variant, the k_cat values for its peroxidase activity increase by 6- to 15-fold in the pH range of 5-8 versus those of the WT protein. These data along with previously reported data for the other THC4-linked variants, G41S and Y48H, underscore the role of Ω-loop C (residues 40-57) in modulating the peroxidase activity of cytochrome c early in apoptosis.

Heightened Dynamics of the Oxidized Y48H Variant of Human Cytochrome c Increases Its Peroxidatic Activity

Proteins performing multiple biochemical functions are called "moonlighting proteins" or extreme multifunctional (EMF) proteins. Mitochondrial cytochrome c is an EMF protein that binds multiple partner proteins to act as a signaling molecule, transfers electrons in the respiratory chain, and acts as a peroxidase in apoptosis. Mutations in the cytochrome c gene lead to the disease thrombocytopenia, which is accompanied by enhanced apoptotic activity. The Y48H variant arises from one such mutation and is found in the 40-57 Ω-loop, the lowest-unfolding free energy substructure of the cytochrome c fold. A 1.36 Å resolution X-ray structure of the Y48H variant reveals minimal structural changes compared to the wild-type structure, with the axial Met80 ligand coordinated to the heme iron. Despite this, the intrinsic peroxidase activity is enhanced, implying that a pentacoordinate heme state is more prevalent in the Y48H variant, corroborated through determination of a Met80 "off rate" of >125 s^-1 compared to a rate of ∼6 s^-1 for the wild-type protein. Heteronuclear nuclear magnetic resonance measurements with the oxidized Y48H variant reveal heightened dynamics in the 40-57 Ω-loop and the Met80-containing 71-85 Ω-loop relative to the wild-type protein, illustrating communication between these substructures. Placed into context with the G41S cytochrome c variant, also implicated in thrombocytopenia, a dynamic picture associated with this disease relative to cytochrome c is emerging whereby increasing dynamics in substructures of the cytochrome c fold serve to facilitate an increased population of the peroxidatic pentacoordinate heme state in the following order: wild type < G41S < Y48H.

The response of Ω-loop D dynamics to truncation of trimethyllysine 72 of yeast iso-1-cytochrome c depends on the nature of loop deformation

Trimethyllysine 72 (tmK72) has been suggested to play a role in sterically constraining the heme crevice dynamics of yeast iso-1-cytochrome c mediated by the Ω-loop D cooperative substructure (residues 70-85). A tmK72A mutation causes a gain in peroxidase activity, a function of cytochrome c that is important early in apoptosis. More than one higher energy state is accessible for the Ω-loop D substructure via tier 0 dynamics. Two of these are alkaline conformers mediated by Lys73 and Lys79. In the current work, the effect of the tmK72A mutation on the thermodynamic and kinetic properties of wild-type iso-1-cytochrome c (yWT versus WT*) and on variants carrying a K73H mutation (yWT/K73H versus WT*/K73H) is studied. Whereas the tmK72A mutation confers increased peroxidase activity in wild-type yeast iso-1-cytochrome c and increased dynamics for formation of a previously studied His79-heme alkaline conformer, the tmK72A mutation speeds return of the His73-heme alkaline conformer to the native state through destabilization of the His73-heme alkaline conformer relative to the native conformer. These opposing behaviors demonstrate that the response of the dynamics of a protein substructure to mutation depends on the nature of the perturbation to the substructure. For a protein substructure which mediates more than one function of a protein through multiple non-native structures, a mutation could change the partitioning between these functions. The current results suggest that the tier 0 dynamics of Ω-loop D that mediates peroxidase activity has similarities to the tier 0 dynamics required to form the His79-heme alkaline conformer.

Native state dynamics affects the folding transition of porcine pancreatic phospholipase A2

Porcine pancreatic phospholipase A2, a small and disulfide rich protein, is extremely resistant against chemically or thermally induced unfolding. Despite this marked resistance, the protein displays broad unfolding transitions resulting in comparatively low apparent thermodynamic stability. Broad unfolding transitions may result from undetected folding intermediates, residual structures in the unfolded state or an inhomogeneity of the native state. Using circular dichroism, fluorescence, and NMR spectroscopy, we ruled out the existence of stably populated folding intermediates, whereas UV absorbance measurements hinted at stable residual structures in the unfolded state. These residual structures proved, however, to have no impact on the folding parameters. Studies by limited proteolysis, CD, and NMR spectroscopy under non-denaturing conditions suggested pronounced dynamics of the protein in the native state, which as long as unrestrained by acidic pH or bound Ca(2+) ions exert considerable influence on the unfolding transition.

Effect of an Ala81His mutation on the Met80 loop dynamics of iso-1-cytochrome c

An A81H variant of yeast iso-1-cytochrome c is prepared to test the hypothesis that the steric size of the amino acid at sequence position 81 of cytochrome c, which has evolved from Ala in yeast to Ile in mammals, slows the dynamics of the opening of the heme crevice. The A81H mutation is used both to increase steric size and to provide a probe of the dynamics of the heme crevice through measurement of the thermodynamics and kinetics of the His81-mediated alkaline conformational transition of A81H iso-1-cytochrome c. Thermodynamic measurements show that the native conformer is more stable than the His81-heme alkaline conformer for A81H iso-1-cytochrome c. ΔGu°(H2O) is approximately 1.9 kcal/mol for formation of the His81-heme alkaline conformer. By contrast, for K79H iso-1-cytochrome c, the native conformer is less stable than the His79-heme alkaline conformer. ΔGu°(H2O) is approximately -0.34 kcal/mol for formation of the His79-heme alkaline conformer. pH jump and gated electron transfer kinetics demonstrate that this stabilization of the native conformer in A81H iso-1-cytochrome c arises primarily from a decrease in the rate constant for formation of the His81-heme alkaline conformer, kf,His81, relative to kf,His79 for formation of the His79-heme alkaline conformer, which forms by a mechanism similar to that observed for the His81-heme alkaline conformer. The result is discussed in terms of the effect of global protein stability on protein dynamics and in terms of optimization of the sequence of cytochrome c for its role as a peroxidase in the early stages of apoptosis in higher eukaryotes.

Mutation of trimethyllysine 72 to alanine enhances His79-heme-mediated dynamics of iso-1-cytochrome c

Trimethyllysine 72 (Tml72) of yeast iso-1-cytochrome c lies across the surface of the heme crevice loop (Ω-loop D, residues 70-85) like a brace. Lys72 is oriented similarly in horse cytochrome c (Cytc). To determine whether this residue affects the dynamics of opening the heme crevice loop, we have studied the effect of a Tml72 to Ala substitution on the formation of the His79-heme alkaline conformer near neutral pH using a variant of iso-1-Cytc including K72A and K79H mutations. Guanidine hydrochloride denaturation shows that the Tml72 to Ala substitution within error does not affect the global stability of the protein. The effect of the Tml72 to Ala substitution on the thermodynamics of the His79-heme alkaline transition is also small. However, pH-jump kinetic studies of the His79-heme alkaline transition show that both the forward and backward rates of conformational change are increased by the Tml72 to Ala substitution. The barrier for opening the heme crevice is reduced by 0.5 kcal/mol and for closing the heme crevice by 0.3 kcal/mol. The ability of Tml72 to modulate the heme crevice dynamics may indicate a crucial role in regulating function, such as in the peroxidase activity seen in the early stages of apoptosis.

Thermodynamics of loop formation in the denatured state of rhodopseudomonas palustris cytochrome c': scaling exponents and the reconciliation problem

The observation that denatured proteins yield scaling exponents, nu, consistent with random-coil behavior and yet can also have pockets of residual or nonrandom structure has been termed the "reconciliation problem". To provide greater insight into the denatured state of a foldable sequence, we have measured histidine-heme loop formation equilibria in the denatured state of a class II c-type cytochrome, cytochrome c' from Rhodopseudomonas palustris. We have prepared a series of variants that provide His-heme loop stabilities, pK(loop)(His), for loop sizes ranging from 10 to 111 residues at intervals of 7 to 11 residues along the sequence of the protein. We observe a scaling exponent for loop formation, nu(3), of 2.5+/-0.3. Theoretical values for nu(3) range from 1.8 to 2.4; thus, the observed nu(3) is consistent with random-coil behavior. However, in contrast to data for loop formation as a function of loop size obtained with peptides of homogeneous sequence, we observe considerable scatter about the linear dependence of loop stability on loop size. Thus, foldable sequences behave very differently from homogeneous peptide sequences. The observed scatter suggests that there is considerable variation in the conformational properties along the backbone of a foldable sequence, consistent with alternating compact and extended regions. With regard to the reconciliation problem, it is evident that a scaling exponent consistent with a random coil is necessary but not sufficient to demonstrate random-coil behavior.

Compressing the free energy range of substructure stabilities in iso-1-cytochrome c

Evolutionary conservation of substructure architecture between yeast iso-1-cytochrome c and the well-characterized horse cytochrome c is studied with limited proteolysis, the alkaline conformational transition and global unfolding with guanidine-HCl. Mass spectral analysis of limited proteolysis cleavage products for iso-1-cytochrome c show that its least stable substructure is the same as horse cytochrome c. The limited proteolysis data yield a free energy of 3.8 +/- 0.4 kcal mol(-1) to unfold the least stable substructure compared with 5.05 +/- 0.30 kcal mol(-1) for global unfolding of iso-1-cytochrome c. Thus, substructure stabilities of iso-1-cytochrome c span only approximately 1.2 kcal mol(-1) compared with approximately 8 kcal mol(-1) for horse cytochrome c. Consistent with the less cooperative folding thus expected for the horse protein, the guanidine-HCl m-values are approximately 3 kcal mol(-1)M(-1) versus approximately 4.5 kcal mol(-1)M(-1) for horse versus yeast cytochrome c. The tight free energy spacing of the yeast cytochrome c substructures suggests that its folding has more branch points than for horse cytochrome c. Studies on a variant of iso-1-cytochrome c with an H26N mutation indicate that the least and most stable substructures unfold sequentially and the two least stable substructures unfold independently as for horse cytochrome c. Thus, important aspects of the substructure architecture of horse cytochrome c, albeit compressed energetically, are preserved evolutionally in yeast iso-1-cytochrome c.

ATP acts as a regulatory effector in modulating structural transitions of cytochrome c: implications for apoptotic activity

The binding of lipids (free fatty acids as well as acidic phospholipids) to cytochrome c (cyt c) induces conformational changes and partial unfolding of the protein, strongly influencing cyt c oxidase/peroxidase activity. ATP is unique among the nucleotides in being able to turn non-native states of cyt c back to the native conformation. The peroxidase activity acquired by lipid-bound cyt c turns out to be very critical in the early stages of apoptosis. Nucleotide specificity is observed for apoptosome formation and caspase activation, the cleavage occurring only in the presence of dATP or ATP. In this study, we demonstrate the connection between peroxidase activity and oleic acid-induced conformational transitions of cyt c and show how ATP is capable of modulating such interplay. By NMR measurement, we have demonstrated that ATP interacts with a site (S1) formed by K88, R91, and E62 and such interaction was weakened by mutation of E62, suggesting the selective role in the interaction played by the base moiety. Interestingly, the interactions of ATP and GTP with cyt c are significantly different at low nucleotide concentrations, with GTP being less effective in perturbing the S1 site and in eliciting apoptotic activity. To gain insights into the structural features of cyt c required for its pro-apoptotic activity and to demonstrate a regulatory role for ATP (compared to the effect of GTP), we have performed experiments on cell lysates by using cyt c proteins mutated on amino acid residues that, as suggested by NMR measurements, belong to S1. Thus, we provide evidence that ATP acts as an allosteric effector, regulating structural transitions among different conformations and different oxidation states of cyt c, which are endowed with apoptotic activity or not. On this basis, we suggest a previously unrecognized role for ATP binding to cyt c at low millimolar concentrations in the cytosol, beyond the known regulatory role during the oxidative phosphorylation in mitochondria.

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.

Zinc porphyrin: a fluorescent acceptor in studies of Zn-cytochrome c unfolding by fluorescence resonance energy transfer

FRET between the zinc porphyrin (ZnP) chromophore in zinc-substituted cytochrome c (Zn-cyt c) and an Alexa Fluor dye attached to specific surface sites was used to characterize Zn-cyt c unfolding. The use of ZnP as a fluorescent acceptor eliminates the need to doubly label the protein with exogenous dyes to perform FRET experiments in which both donor and acceptor fluorescence is monitored. The requirement for attachment of only one dye also minimizes perturbation to the protein. This sensitive technique allowed for the determination of distances between the label placed at six different sites and ZnP through a range of denaturant concentrations. Fitting of the data to a three-state model provides distances in the unfolding intermediate. The use of ZnP as a fluorescent acceptor of energy in FRET has a significant potential for application to a range of other systems including heme-binding proteins and proteins to which a covalently attached heme tag may be added.

Probing the bottom of a folding funnel using conformationally gated electron transfer reactions

The effect of global stability on the kinetics of interconversion between the native (N) and a compact, partially unfolded form (I) of iso-1-cytochrome c stabilized by His73-heme ligation is investigated using a novel conformationally gated ET method. For the K73H variant and the 2-fold less stable AcH73 variant, the N and I conformers are of nearly equal stability at pH 7.5. The pH jump kinetic data yield kobs = kNI + kIN of 35-40 s-1 at final pH values from 6 to 8 for the AcH73 variant, about 3-fold faster than for the more stable K73H variant. Gated ET measurements give kNI = 28 s-1 and kIN = 13 s-1 for the AcH73 variant, 10- and 2-fold greater than that for the more stable K73H variant. Thus, funneled landscapes have evolved such that loss of global stability lowers barriers at the bottom of a folding funnel, still allowing for efficient folding.

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.

Thermodynamics of protein denatured states

Recent work on the thermodynamics of protein denatured states is providing insight into the stability of residual structure and the conformational constraints that affect the disordered states of proteins. Current data from native state hydrogen exchange and the pH dependence of protein stability indicate that residual structure can modulate the stability of the denatured state by up to 4 kcal mol(-1). NMR structural data have emphasized the role of hydrophobic clusters in stabilizing denatured state residual structures, however recent results indicate that electrostatic interactions, both favorable and unfavorable, are also important modulators of the stability of the denatured state. Thermodynamics methods that take advantage of histidine-heme ligation chemistry have also been developed to probe the conformational constraints that act on denatured states. These methods have provided insights into the role of excluded volume, chain stiffness, and loop persistence in modulating the conformational preferences of highly disordered proteins. New insights into protein folding and novel methods to manipulate protein stability are emerging from this work.

Tuning the Rate and pH Accessibility of a Conformational Electron Transfer Gate

Methods to fine-tune the rate of a fast conformational electron transfer (ET) gate involving a His-heme alkaline conformer of iso-1-cytochrome c (iso-1-Cytc) and to adjust the pH accessibility of a slow ET gate involving a Lys-heme alkaline conformer are described. Fine-tuning the fast ET gate employs a strategy of making surface mutations in a substructure unfolded in the alkaline conformer. To make the slow ET gate accessible at neutral pH, the strategy involves mutations at buried sequence positions which are expected to more strongly perturb the stability of native versus alkaline iso-1-Cytc. To fine-tune the rate of the fast His 73-heme ET gate, we mutate the surface-exposed Lys 79 to Ala (A79H73 variant). This mutation also simplifies ET gating by removing Lys 79, which can serve as a ligand in the alkaline conformer of iso-1-Cytc. To adjust the pH accessibility of the slow Lys 73-heme ET gate, we convert the buried side chain Asn 52 to Gly and also mutate Lys 79 to Ala to simplify ET gating (A79G52 variant). ET kinetics is studied as a function of pH using hexaammineruthenium(II) chloride (a6Ru2+) to reduce the variants. Both variants show fast direct ET reactions dependent on [a6Ru2+] and slower gated ET reactions that are independent of [a6Ru2+]. The observed gated ET rates correlate well with rates for the alkaline-to-native state conformational change measured independently. Together with the previously reported H73 variant (Baddam, S.; Bowler, B. E. J. Am. Chem. Soc. 2005, 127, 9702-9703), the A79H73 variant allows His 73-heme-mediated ET gating to be fine-tuned from 75 to 200 ms. The slower Lys 73-heme (15-20 s time scale) ET gate for the A79G52 variant is now accessible over the pH range 6-8.

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.

Conformationally gated electron transfer in iso-1-cytochrome c: engineering the rate of a conformational switch

An engineered form of iso-1-cytochrome c with lysine 73 mutated to histidine is shown to increase by nearly 500-fold the rate of a conformational gate that modulates the rate of electron transfer into this protein. This result demonstrates the potential of protein engineering to provide electron transfer gates with tailored properties. The pH dependence of the rate of the conformational electron transfer gate correlates well with the pH dependence of the conformational change from a His 73-ligated heme to a Met 80-ligated heme, determined independently by pH jump methods, allowing unambiguous assignment of the conformational electron transfer gating step. The rate of the electron transfer gate is also modulated by a cis to trans proline isomerization, indicating that both amino acid sequence and the nature of the heme ligand provide avenues for rational design of electron transfer gates which open at different rates.

Proton-mediated dynamics of the alkaline conformational transition of yeast iso-1-cytochrome c

The kinetics of the alkaline conformational transition of a Lys 73-->His variant of iso-1-cytochrome c have been investigated using pH jump stopped-flow methods to probe the nature of the ionizable "trigger" group for this conformational change. This mutation moves the pK(a) of the ligand replacing Met 80 from about 10.5 to approximately 6.6 and has unmasked two other ionizable groups, besides the ligand replacing Met 80, that modulate the kinetics of this process. The results are discussed in terms of the impact of ionization equilibria on protein folding mechanisms.

Redesigning the Folding Energetics of a Model Three-helix Bundle Protein by Site-directed Mutagenesis

Because of their limited size and complexity, de novo designed proteins are particularly useful for the detailed investigation of folding thermodynamics and mechanisms. Here, we describe how subtle changes in the hydrophobic core of a model three-helix bundle protein (GM-0) alter its folding energetics. To explore the folding tolerance of GM-0 toward amino acid sequence variability, two mutant proteins (GM-1 and GM-2) were generated. In the mutants, cavities were created in the hydrophobic core of the protein by either singly (GM-1; L35A variant) or doubly (GM-2; L35A/I39A variant) replacing large hydrophobic side chains by smaller Ala residues. The folding of GM-0 is characterized by two partially folded intermediate states exhibiting characteristics of molten globules, as evidenced by pressure-unfolding and pressure-assisted cold denaturation experiments. In contrast, the folding energetics of both mutants, GM-1 and GM-2, exhibit only one folding intermediate. Our results support the view that simple but biologically important folding motifs such as the three-helix bundle can reveal complex folding plasticity, and they point to the role of hydrophobic packing as a determinant of the overall stability and folding thermodynamic of the helix bundle.

Extended application of gel-permeation chromatography by spin column

Separation of small volumes of proteins from unbound ligands or reequilibration with buffer by passing through a 1-ml Sephadex G-50 column under mild centrifugal force is a popular technique. Here it has been demonstrated that other Sephadex matrix could similarly be used for complete or partial separation of protein molecules. Proteins to be eluted at void volume are recovered near quantitatively, while others are partly or almost completely retained depending on molecular size. Calibration curves using standard proteins of Mw 12.5 to 440 kDa with Sephadex G-50-G-200 representing recovery versus molecular weight show profiles as expected from the fractionation ranges of the column matrix. The procedure may be applied to follow protein association-dissociation reactions if the molecular weights of the species concerned are known and a proper matrix exists for separating them. Equilibrium unfolding transitions constructed with model proteins in presence of 0-8M urea using recovery as an index correspond to profiles obtained from other physical measurements. This may be a convenient approach to follow change of protein hydrodynamic volume quickly when a parallel methodology is not readily available.

Denatured state thermodynamics: residual structure, chain stiffness and scaling factors

A set of nine variants of yeast iso-1-cytochrome c with zero or one surface histidine have been engineered such that the N-terminal amino group is acetylated in vivo. N-terminal acetylation has been confirmed by mass spectral analysis of intact and proteolytically digested protein. The histidine-heme loop-forming equilibrium, under denaturing conditions (3 M guanidine hydrochloride), has been measured by pH titration providing an observed pK(a), pK(a)(obs), for each variant. N-terminal acetylation prevents the N-terminal amino group-heme binding equilibrium from interfering with measurements of histidine-heme affinity. Significant deviation is observed from the linear dependence of pK(a)(obs) on the log of the number of monomers in the loop formed, expected for a random coil denatured state. The maximum histidine-heme affinity occurs for a loop size of 37 monomers. For loop sizes of 37-83 monomers, histidine-heme pK(a)(obs) values are consistent with a scaling factor of -4.2+/-0.3. This value is much larger than the scaling factor of -1.5 for a freely jointed random coil, which is commonly used to represent the conformational properties of protein denatured states. For loop sizes of nine to 22 monomers, chain stiffness is likely responsible for the decreases in histidine-heme affinity relative to a loop size of 37. The results are discussed in terms of residual structure and sequence composition effects on the conformational properties of the denatured states of proteins.

Cofactor and tryptophan accessibility and unfolding of brain glutamate decarboxylase

Cofactor and tryptophan accessibility of the 65-kDa form of rat brain glutamate decarboxylase (GAD) was investigated by fluorescence quenching measurements using acrylamide, I-, and Cs+ as the quenchers. Trp residues were partially exposed to solvent. I- was less able and Cs+ was more able to quench the fluorescence of Trp residues in the holoenzyme of GAD (holoGAD) than the apoenzyme (apoGAD). The fraction of exposed Trp residues were in the range of 30-49%. In contrast, pyridoxal-P bound to the active site of GAD was exposed to solvent. I- was more able and Cs+ was less able to quench the fluorescence of pyridoxal-P in holoGAD. The cofactor was present in a positively charged microenvironment, making it accessible for interactions with anions. A difference in the exposure of Trp residues and pyridoxal-P to these charged quenchers suggested that the exposed Trp residues were essentially located outside of the active site. Changes in the accessibility of Trp residues upon pyridoxal-P binding strongly supported a significant conformational change in GAD. Fluorescence intensity measurements were also carried out to investigate the unfolding of GAD using guanidine hydrochloride (GdnHCl) as the denaturant. At 0.8-1.5 M GdnHCl, an intermediate step was observed during the unfolding of GAD from the native to the denatured state, and was not found during the refolding of GAD from the denatured to native state, indicating that this intermediate step was not a reversible process. However, at >1.5 M GdnHCl for holoGAD and >2.0 M GdnHCl for apoGAD, the transition leading to the denatured state was reversible. It was suggested that the intermediate step involved the dissociation of native dimer of GAD into monomers and the change in the secondary structure of the protein. Circular dichroism revealed a decrease in the alpha-helix content of GAD from 36 to 28%. The unfolding pattern suggested that GAD may consist of at least two unfolding domains. Unfolding of the lower GdnHCl-resisting domain occurred at a similar concentration of denaturant for apoGAD and holoGAD, while unfolding of the higher GdnHCl-resisting domain occurred at a higher concentration of GdnHCl for apoGAD than holoGAD.

Core and surface mutations affect folding kinetics, stability and cooperativity in IL-1 beta: does alteration in buried water play a role?

Interleukin-1 beta (IL-1 beta) is a cytokine and a member of the beta-trefoil superfamily of protein structures. An interesting feature in the folding of IL-1 beta, shared with some other members of the same topological family, is the existence of a slow step in folding to the native conformation from a discrete intermediate. Wanting to probe the nature of this slow step in the folding of WT IL-1 beta (tau(1)=45 seconds), we made ten sequence variants of IL-1 beta (L10A, T9Q, T9G, C8S, C8A, N7G, N7D, L6A, R4P, and R4Q), where all mutations are located along strand 1. This strand is not protected from hydrogen exchange until late in folding. Most of the mutations showed little effect on the kinetics of folding for IL-1 beta. However, C8 is clearly involved in both the late and the early steps in folding, while sequence variants at L10 and L6 affect only late events in folding. The value of the slowest relaxation time, tau(1), which is associated with the rate of native protein formation, increased for the refolding of C8S, while C8A, L6A, and L10A showed smaller but systematic increases in the value of tau(1.)For both C8S and C8A, the value of the step associated with formation of the intermediate, tau(2), was independent of denaturant concentration. In addition, mutations in the hydrophobic core (L10A, C8A, C8S, and L6A) and, surprisingly, along the surface (T9G, T9Q, and N7G) alter the stability. The most destabilizing mutations show changes in equilibrium unfolding cooperativity, which is atypical for destabilizing mutations in IL-1 beta. Crystallographic studies indicate that mutations along strand 1 may alter the number of ordered water molecules within the core. Thus, side-chain replacement in this region can disrupt essential main-chain interactions mediated by ordered water contacts in a highly cooperative network of hydrogen bonding.

NMR investigation of ferricytochrome c unfolding: detection of an equilibrium unfolding intermediate and residual structure in the denatured state

Horse ferricytochrome c (cyt c) undergoes exchange of one of its axial heme ligands (Met-80) for one or more non-native ligands under denaturing conditions. We have used (1)H NMR spectroscopy to detect two conformations of paramagnetic cyt c with non-native heme ligation through a range of urea concentrations. One non-native form is an equilibrium unfolding intermediate observed under partially denaturing conditions and is attributed to replacement of Met-80 with one or more Lys side chains. The second non-native form, in which the native Met ligand is replaced by a His, is observed under strongly denaturing conditions. Thermodynamic analysis of these data indicates a relatively small DeltaG (17 kJ/mol) for the transition from native to the Lys-ligated intermediate and a significantly larger DeltaG (47 kJ/mol) for the transition from native to the His-ligated species. Although CD and fluorescence data indicate that the equilibrium unfolding of cyt c is a two-state process, these NMR results implicate an intermediate with His-Lys ligation.

Measuring denatured state energetics: deviations from random coil behavior and implications for the folding of iso-1-cytochrome c

The changes in the free energy of the denatured state of a set of yeast iso-1-cytochrome c variants with single surface histidine residues have been measured in 3 M guanidine hydrochloride. The thermodynamics of unfolding by guanidine hydrochloride is also reported. All variants have decreased stability relative to the wild-type protein. The free energy of the denatured state was determined in 3 M guanidine hydrochloride by evaluating the strength of heme-histidine ligation through determination of the pK(a) for loss of histidine binding to the heme. The data are corrected for the presence of the N-terminal amino group which also ligates to the heme under similar solution conditions. Significant deviations from random coil behavior are observed. Relative to a variant with a single histidine at position 26, residual structure of the order of -1.0 to -2.5 kcal/mol is seen for the other variants studied. The data explain the slower folding of yeast iso-1-cytochrome c relative to the horse protein. The greater number of histidines and the greater strength of ligation are expected to slow conversion of the histidine-misligated forms to the obligatory aquo-heme intermediate during the ligand exchange phase of folding. The particularly strong association of histidine residues at positions 54 and 89 may indicate regions of the protein with strong energetic propensities to collapse against the heme during early folding events, consistent with available data in the literature on early folding events for cytochrome c.

The magnitude of changes in guanidine-HCl unfolding m-values in the protein, iso-1-cytochrome c, depends upon the substructure containing the mutation

Hydrophilic to hydrophobic mutations have been made at 11 solvent exposed sites on the surface of iso-1-cytochrome c. Most of these mutations involve the replacement of lysine with methionine, which is nearly isosteric with lysine. Minimal perturbation to the native structure is expected, and this expectation is confirmed by infrared amide I spectroscopy. Guanidine hydrochloride denaturation studies demonstrate that these variants affect the magnitude of the m-value, the rate of change of free energy with respect to denaturant concentration, to different degrees. Changes in m-values are indicative of changes in the equilibrium folding mechanism of a protein. Decreases in m-values are normally thought to result either from an increased population of intermediates during unfolding or from a more compact denatured state. When cytochrome c is considered in terms of its thermodynamic substructures, the changes in the m-value for a given variant appear to depend upon the substructure in which the mutation is made. These data indicate that the relative stabilities and physical properties of substructures of cytochrome c play an important determining role in the equilibrium folding mechanism of this protein.

Cytochrome c folding traps are not due solely to histidine-heme ligation: direct demonstration of a role for N-terminal amino group-heme ligation

In previous work, heme ligation effects on the folding of cytochrome c have been attributed to histidine side-chains. A variant of yeast iso-1-cytochrome c designated TM, which lacks all histidine residues except His18, still shows evidence of denatured state heme ligation in the pH range between 5 and 6 where normally only histidine ligation is expected. Conversion of the N-terminal amino group of TM to a carbonyl group through a transamination reaction with glyoxylate produced a protein (ModTM) with no terminal amino group. The midpoint pH (pH1/2) for loss of heme ligation in 3 M guanidine-HCl shifts from 5.9 to 7.4 as a result of this modification, providing direct evidence for N-terminal amino group-heme ligation under these conditions. The N-terminal amino group thus competes with histidine for misligation of iso-1-cytochrome c under denaturing conditions. To assess the effect of denatured state N-terminal amino group-heme ligation on the folding of iso-1-cytochrome c, stopped-flow kinetics experiments were conducted. At pH 6.2, the major refolding lifetimes (3 M-->0.27 M guanidine-HCl) for ModTM, TM and the wild-type protein are 11.6 ms, 30 ms and 1.3 seconds, respectively. Denatured state ligation of the N-terminal amino group thus slows folding 2.6-fold.

A histidine variant of yeast iso-1-cytochrome c that strongly affects the energetics of the denatured state

Iso-1-cytochrome c has been engineered to remove all histidine residues not involved in heme ligation in the native state to produce a variant designated TM. Single histidine residues were then introduced at positions 26, TM + His26, and 54, TM + His54. Since histidine residues not involved in native state heme ligation are known to replace the methionine 80 heme ligand in denatured cytochrome c, these variants were expected to affect the structure of the denatured state. Guanidine hydrochloride denaturations were performed to assess the stability of these proteins relative to the wild-type protein. The free energy difference for heme ligation in the denatured state was assessed by pH titration. The experimentally observed mutation-induced change (delta deltaG(D-state)) in the free energy of heme ligation for unfolded TM + His54 versus TM + His26 is -0.4 kcal/mol. The expected mutation-induced change in delta deltaG(D-state) calculated for a random coil unfolded state is +2 kcal/mol. Thus, unfolded TM + His54 has residual structure stabilizing its denatured state by -2.4 kcal/mol relative to TM + His26. The results imply that the denatured state can contribute significantly to mutation-induced changes in the free energy of unfolding of a protein.