c-di-GMP regulates activity of the PlzA RNA chaperone from the Lyme disease spirochete.

PlzA is a c-di-GMP-binding protein crucial for adaptation of the Lyme disease spirochete Borrelia (Borreliella) burgdorferi during its enzootic life cycle. Unliganded apo-PlzA is important for vertebrate infection, while liganded holo-PlzA is important for survival in the tick; however, the biological function of PlzA has remained enigmatic. Here, we report that PlzA has RNA chaperone activity that is inhibited by c-di-GMP binding. Holo- and apo-PlzA bind RNA and accelerate RNA annealing, while only apo-PlzA can strand displace and unwind double-stranded RNA. Guided by the crystal structure of PlzA, we identified several key aromatic amino acids protruding from the N- and C-terminal domains that are required for RNA-binding and unwinding activity. Our findings illuminate c-di-GMP as a switch controlling the RNA chaperone activity of PlzA, and we propose that complex RNA-mediated modulatory mechanisms allow PlzA to regulate gene expression during both the vector and host phases of the B. burgdorferi life cycle.

High-Accuracy Prediction of Stabilizing Surface Mutations to the Three-Helix Bundle, UBA(1), with EmCAST.

The accurate modeling of energetic contributions to protein structure is a fundamental challenge in computational approaches to protein analysis and design. We describe a general computational method, EmCAST (empirical Cα stabilization), to score and optimize the sequence to the structure in proteins. The method relies on an empirical potential derived from the database of the Cα dihedral angle preferences for all possible four-residue sequences, using the data available in the Protein Data Bank. Our method produces stability predictions that naturally correlate one-to-one with the experimental results for solvent-exposed mutation sites. EmCAST predicted four mutations that increased the stability of a three-helix bundle, UBA(1), from 2.4 to 4.8 kcal/mol by optimizing residues in both helices and turns. For a set of eight variants, the predicted and experimental stabilizations correlate very well (R2 = 0.97) with a slope near 1 and with a 0.16 kcal/mol standard error for EmCAST predictions. Tests against literature data for the stability effects of surface-exposed mutations show that EmCAST outperforms the existing stability prediction methods. UBA(1) variants were crystallized to verify and analyze their structures at an atomic resolution. Thermodynamic and kinetic folding experiments were performed to determine the magnitude and mechanism of stabilization. Our method has the potential to enable the rapid, rational optimization of natural proteins, expand the analysis of the sequence/structure relationship, and supplement the existing protein design strategies.

Residual Structure in the Denatured State of the Fast-Folding UBA(1) Domain from the Human DNA Excision Repair Protein HHR23A.

The structure of the first ubiquitin-associated domain from HHR23A, UBA(1), was determined by X-ray crystallography at a 1.60 Å resolution, and its stability, folding kinetics, and residual structure under denaturing conditions have been investigated. The concentration dependence of thermal denaturation and size-exclusion chromatography indicate that UBA(1) is monomeric. Guanidine hydrochloride (GdnHCl) denaturation experiments reveal that the unfolding free energy, ΔGu°'(H2O), of UBA(1) is 2.4 kcal mol-1. Stopped-flow folding kinetics indicates sub-millisecond folding with only proline isomerization phases detectable at 25 °C. The full folding kinetics are observable at 4 °C, yielding a folding rate constant, kf, in the absence of a denaturant of 13,000 s-1 and a Tanford ß-value of 0.80, consistent with a compact transition state. Evaluation of the secondary structure via circular dichroism shows that the residual helical structure in the denatured state is replaced by polyproline II structure as the GdnHCl concentration increases. Analysis of NMR secondary chemical shifts for backbone 15NH, 13CO, and 13Cα atoms between 4 and 7 M GdnHCl shows three islands of residual helical secondary structure that align in sequence with the three native-state helices. Extrapolation of the NMR data to 0 M GdnHCl demonstrates that helical structure would populate to 17-33% in the denatured state under folding conditions. Comparison with NMR data for a peptide corresponding to helix 1 indicates that this helix is stabilized by transient tertiary interactions in the denatured state of UBA(1). The high helical content in the denatured state, which is enhanced by transient tertiary interactions, suggests a diffusion-collision folding mechanism.

Alkaline State of the Domain-Swapped Dimer of Human Cytochrome c: A Conformational Switch for Apoptotic Peroxidase Activity.

A 2.08 Å structure of an alkaline conformer of the domain-swapped dimer of K72A human cytochrome c (Cytc) crystallized at pH 9.9 is presented. In the structure, Lys79 is ligated to the heme. All other domain-swapped dimer structures of Cytc have water bound to this coordination site. Part of Ω-loop D (residues 70-85) forms a flexible linker between the subunits in other Cytc domain-swapped dimer structures but instead converts to a helix in the alkaline conformer of the dimer combining with the C-terminal helix to form two 26-residue helices that bracket both sides of the dimer. The alkaline transition of the K72A human dimer monitored at both 625 nm (high spin heme) and 695 nm (Met80 ligation) yields midpoint pH values of 6.6 and 7.6, respectively, showing that the Met80 → Lys79 and high spin to low spin transitions are distinct. The dimer peroxidase activity increases rapidly below pH 7, suggesting that population of the high spin form of the heme is what promotes peroxidase activity. Comparison of the structures of the alkaline dimer and the neutral pH dimer shows that the neutral pH conformer has a better electrostatic surface for binding to a cardiolipin-containing membrane and provides better access for small molecules to the heme iron. Given that the pH of mitochondrial cristae ranges from 6.9 to 7.2, the alkaline transition of the Cytc dimer could provide a conformational switch to tune the peroxidase activity of Cytc that oxygenates cardiolipin in the early stages of apoptosis.

Denatured State Conformational Biases in Three-Helix Bundles Containing Divergent Sequences Localize near Turns and Helix Capping Residues.

Rhodopseudomonas palustris cytochrome c', a four-helix bundle, and the second ubiquitin-associated domain, UBA(2), a three-helix bundle from the human homologue of yeast Rad23, HHR23A, deviate from random coil behavior under denaturing conditions in a fold-specific manner. The random coil deviations in each of these folds occur near interhelical turns and loops in their tertiary structures. Here, we examine an additional three-helix bundle with an identical fold to UBA(2), but a highly divergent sequence, the first ubiquitin-associated domain, UBA(1), of HHR23A. We use histidine-heme loop formation methods, employing eight single histidine variants, to probe for denatured state conformational bias of a UBA(1) domain fused to the N-terminus of iso-1-cytochrome c (iso-1-Cytc). Guanidine hydrochloride (GuHCl) denaturation shows that the iso-1-Cytc domain unfolds first, followed by the UBA(1) domain. Denatured state (4 and 6 M GuHCl) histidine-heme loop formation studies show that as the size of the histidine-heme loop increases, loop stability decreases, as expected for the Jacobson-Stockmayer relationship. However, loops formed with His35, His31, and His15, of UBA(1), are 0.6-1.1 kcal/mol more stable than expected from the Jacobson-Stockmayer relationship, confirming the importance of deviations of the denatured state from random coil behavior near interhelical turns of helical domains for facilitating folding to the correct topology. For UBA(1) and UBA(2), hydrophobic clusters on either side of the turns partially explain deviations from random coil behavior; however, helix capping also appears to be important.

The Human Cytochrome c Domain-Swapped Dimer Facilitates Tight Regulation of Intrinsic Apoptosis.

Oxidation of cardiolipin (CL) by cytochrome c (cytc) has been proposed to initiate the intrinsic pathway of apoptosis. Domain-swapped dimer (DSD) conformations of cytc have been reported both by our laboratory and by others. The DSD is an alternate conformer of cytc that could oxygenate CL early in apoptosis. We demonstrate here that the cytc DSD has a set of properties that would provide tighter regulation of the intrinsic pathway. We show that the human DSD is kinetically more stable than horse and yeast DSDs. Circular dichroism data indicate that the DSD has a less asymmetric heme environment, similar to that seen when the monomeric protein binds to CL vesicles at high lipid:protein ratios. The dimer undergoes the alkaline conformational transition near pH 7.0, 2.5 pH units lower than that of the monomer. Data from fluorescence correlation spectroscopy and fluorescence anisotropy suggest that the alkaline transition of the DSD may act as a switch from a high affinity for CL nanodiscs at pH 7.4 to a much lower affinity at pH 8.0. Additionally, the peroxidase activity of the human DSD increases 7-fold compared to that of the monomer at pH 7 and 8, but by 14-fold at pH 6 when mixed Met80/H2O ligation replaces the lysine ligation of the alkaline state. We also present data that indicate that cytc binding shows a cooperative effect as the concentration of cytc is increased. The DSD appears to have evolved into a pH-inducible switch that provides a means to control activation of apoptosis near pH 7.0.

Curvature-Dependent Binding of Cytochrome c to Cardiolipin.

Cytochrome c binds cardiolipin on the concave surface of the inner mitochondrial membrane, before oxidizing the lipid and initiating the apoptotic pathway. This interaction has been studied in vitro, where mimicking the membrane curvature of the binding environment is difficult. Here we report binding to concave, cardiolipin-containing, membrane surfaces and compare findings to convex binding under the same conditions. For binding to the convex outer surface of cardiolipin-containing vesicles, a two-step structural rearrangement is observed with a small rearrangement detectable by Soret circular dichroism (CD) occurring at an exposed lipid-to-protein ratio (LPR) near 10 and partial unfolding detectable by Trp59 fluorescence occurring at an exposed LPR near 23. On the concave inner surface of cardiolipin-containing vesicles, the structural transitions monitored by Soret CD and Trp59 fluorescence are coincident and occur at an exposed LPR near 58. On the concave inner surface of mitochondrial cristae, we estimate the LPR of cardiolipin to cytochrome c is between 50 and 100. Thus, cytochrome c may have adapted to its native environment so that it can undergo a conformational change that switches on its peroxidase activity when it binds to CL-containing membranes in the cristae early in apoptosis. Our results show that membrane curvature qualitatively affects peripheral protein-lipid interactions and also highlights the disparity between in vitro binding studies and their physiological counterparts where cone-shaped lipids, like cardiolipin, are involved.

Effect of V83G and I81A Substitutions to Human Cytochrome c on Acid Unfolding and Peroxidase Activity below a Neutral pH.

Mitochondrial cytochrome c is a highly conserved protein in eukaryotes. Certain functions of cytochrome c have been tuned during evolution. For instance, the intrinsic peroxidase activity of human cytochrome c is much lower than that of the yeast counterpart. Structural studies on K72A yeast iso-1-cytochrome c [McClelland, L. J., et al. (2014) Proc. Natl. Acad. Sci. USA, 111, 6648-6653] revealed that residues 81 and 83 in Ω-loop D (residues 70-85) may be gatekeeper residues for the peroxidase activity linked to intrinsic apoptosis. Amino acids at both positions evolve to more sterically demanding amino acids. We hypothesized that residues 81 and 83 evolved such that steric constraints at these positions tune down the peroxidase activity of human cytochrome c. To test this hypothesis, I81A and V83G variants of human cytochrome c were prepared. Our results show that the I81A substitution significantly influences the thermodynamics and kinetics of access to alternate conformers of human cytochrome c, while the V83G substitution has a more modest effect on these properties. The I81A variant also shows a significant enhancement in peroxidase activity, particularly below pH 7, whereas the V83G variant has a similar peroxidase activity to wild-type human cytochrome c. However, neither variant increases the peroxidase activity of human cytochrome c to the level of yeast iso-1-cytochrome c, indicating that other substructures of cytochrome c are also involved in tuning the intrinsic peroxidase activity of mitochondrial cytochrome c.

Helical Propensity Affects the Conformational Properties of the Denatured State of Cytochrome c'.

Changing the helical propensity of a polypeptide sequence might be expected to affect the conformational properties of the denatured state of a protein. To test this hypothesis, alanines at positions 83 and 87 near the center of helix 3 of cytochrome c' from Rhodopseudomonas palustris were mutated to serine to decrease the stability of this helix. A set of 13 single histidine variants in the A83S/A87S background were prepared to permit assessment of the conformational properties of the denatured state using histidine-loop formation in 3 M guanidine hydrochloride. The data are compared with previous histidine-heme loop formation data for wild-type cytochrome c'. As expected, destabilization of helix 3 decreases the global stabilities of the histidine variants in the A83S/A87S background relative to the wild-type background. Loop stability versus loop size data yields a scaling exponent of 2.1 ± 0.2, similar to the value of 2.3 ± 0.2 obtained for wild-type cytochrome c'. However, the stabilities of all histidine-heme loops, which contain the helix 3 sequence segment, are increased in the A83S/A87S background compared to the wild-type background. Rate constants for histidine-heme loop breakage are similar for the wild-type and A83S/A87S variants. However, for histidine-heme loops that contain the helix 3 sequence segment, the rate constants for loop formation increase in the A83S/A87S background compared to the wild-type background. Thus, residual helical structure appears to stiffen the polypeptide chain slowing loop formation in the denatured state. The implications of these results for protein folding mechanisms are discussed.

Electrostatic Constituents of the Interaction of Cardiolipin with Site A of Cytochrome c.

Cytochrome c binds to cardiolipin (CL) on the inner mitochondrial membrane during the initial stages of apoptosis where it oxidizes CL, promoting its release into the cytoplasm where it initiates apoptosis. Previous work has identified interaction sites on cytochrome c involved in the cytochrome c-CL interaction. The contributions of the lysines attributed to site A, the anionic site, are studied here to elucidate the relative importance of each for electrostatic interaction of cytochrome c with CL at pH 8, conditions where site A is dominant. A set of single, double, and quadruple lysine to alanine variants of yeast iso-1-cytochrome c, at sequence positions 72, 73, 86, and 87, show that all contribute to the site A-mediated interaction with CL. All variants experience two sequential structural rearrangements as the lipid to protein ratio (LPR) increases. At a low LPR near 10, all variants undergo a small heme-centered structural change detected by Soret circular dichroism. At higher LPRs ranging from 22 to 34, all variants partially unfold as detected by Trp59 emission. The robustness of the mechanism of interaction to sequential neutralization of the four lysines assigned to site A demonstrates that site A is more extensive than previously supposed. The nature of both structural rearrangements also depends on which lysines constitute site A. The peroxidase activity of cytochrome c in the early stages of apoptosis depends on the nature of structural rearrangement near the heme. Thus, the lysines that comprise site A may have evolved to optimize the peroxidase signaling switch.

Probing Denatured State Conformational Bias in a Three-Helix Bundle, UBA(2), Using a Cytochrome c Fusion Protein.

Previous work with the four-helix-bundle protein cytochrome c' from Rhodopseudomonas palustris using histidine-heme loop formation methods revealed fold-specific deviations from random coil behavior in its denatured state ensemble. To examine the generality of this finding, we extend this work to a three-helix-bundle polypeptide, the second ubiquitin-associated domain, UBA(2), of the human DNA excision repair protein. We use yeast iso-1-cytochrome c as a scaffold, fusing the UBA(2) domain at the N-terminus of iso-1-cytochrome c. We have engineered histidine into highly solvent accessible positions of UBA(2), creating six single histidine variants. Guanidine hydrochloride denaturation studies show that the UBA(2)-cytochrome c fusion protein unfolds in a three-state process with iso-1-cytochrome c unfolding first. Furthermore, engineered histidine residues in UBA(2) strongly destabilize the iso-1-cytochrome c domain. Equilibrium and kinetic histidine-heme loop formation measurements in the denatured state at 4 and 6 M guanidine hydrochloride show that loop stability decreases as the size of the histidine-heme loop increases, in accord with the Jacobson-Stockmayer equation. However, we observe that the His27-heme loop is both more stable than expected from the Jacobson-Stockmayer relationship and breaks more slowly than expected. These results show that the sequence near His27, which is in the reverse turn between helices 2 and 3 of UBA(2), is prone to persistent interactions in the denatured state. Therefore, consistent with our results for cytochrome c', this reverse turn sequence may help to establish the topology of this fold by biasing the conformational distribution of the denatured state.

Rapid quantification of vesicle concentration for DOPG/DOPC and Cardiolipin/DOPC mixed lipid systems of variable composition.

A novel approach to quantify mixed lipid systems is described. Traditional approaches to lipid vesicle quantification are time consuming, require large amounts of material and are destructive. We extend our recently described method for quantification of pure lipid systems to mixed lipid systems. The method only requires a UV-Vis spectrometer and does not destroy sample. Mie scattering data from absorbance measurements are used as input into a Matlab program to calculate the total vesicle concentration and the concentrations of each lipid in the mixed lipid system. The technique is fast and accurate, which is essential for analytical lipid binding experiments.

Site A-Mediated Partial Unfolding of Cytochrome c on Cardiolipin Vesicles Is Species-Dependent and Does Not Require Lys72.

Measurements at pH 8 allow evaluation of binding of 100% cardiolipin vesicles to site A of cytochrome c without interference from other known binding sites. Site A encompasses Lys72, Lys73, Lys86, and Lys87, located in or adjacent to Ω-loop D (residues 70-85), which positions Met80 for binding to the heme. Binding of cytochrome c to cardiolipin disrupts Met80 heme binding, permitting peroxidase activity. Binding of cardiolipin to yeast iso-1-cytochrome c versus human cytochrome c is compared to assess how binding of cardiolipin to site A has evolved for cytochrome c from species that do not have a complete intrinsic apoptotic pathway to species that do. Using a nondestructive method of quantifying cardiolipin concentration, highly reproducible binding curves are obtained. The results indicate two sequential structural rearrangements on the surface of 100% cardiolipin vesicles. The first, more modest, structural rearrangement occurs at an exposed (outer leaflet) lipid:protein ratio of 8-10 for both cytochromes c. The second, occurring at higher lipid:protein ratios, causes significant unfolding of cytochrome c and requires a much higher lipid:protein ratio for human versus yeast cytochrome c. Higher lipid:protein ratios enhance the peroxidase activity of cytochrome c, suggesting that human cytochrome c has evolved a more stringent on/off switch for cardiolipin peroxidation in the early stages of apoptosis. For both human and yeast cytochrome c, the K72A mutation has only minor effects on binding to site A, suggesting that other nearby lysines can compensate for the lack of Lys72.

Effect of a K72A Mutation on the Structure, Stability, Dynamics, and Peroxidase Activity of Human Cytochrome c.

We test the hypothesis that Lys72 suppresses the intrinsic peroxidase activity of human cytochrome c, as observed previously for yeast iso-1-cytochrome c [McClelland, L. J., et al. (2014) Proc. Natl. Acad. Sci. U. S. A. 111, 6648-6653]. A 1.25 Å X-ray structure of K72A human cytochrome c shows that the mutation minimally affects structure. Guanidine hydrochloride denaturation demonstrates that the K72A mutation increases global stability by 0.5 kcal/mol. The K72A mutation also increases the apparent pKa of the alkaline transition, a measure of the stability of the heme crevice, by 0.5 unit. Consistent with the increase in the apparent pKa, the rate of formation of the dominant alkaline conformer decreases, and this conformer is no longer stabilized by proline isomerization. Peroxidase activity measurements show that the K72A mutation increases kcat by 1.6-4-fold at pH 7-10, an effect larger than that seen for the yeast protein. X-ray structures of wild type and K72A human cytochrome c indicate that direct interactions of Lys72 with the far side of Ω-loop D, which are seen in X-ray structures of horse and yeast cytochrome c and could suppress peroxidase activity, are lacking. Instead, we propose that the stronger effect of the K72A mutation on the peroxidase activity of human versus yeast cytochrome c results from relief of steric interactions between the side chains at positions 72 and 81 (Ile in human vs Ala in yeast), which suppress the dynamics of Ω-loop D necessary for the intrinsic peroxidase activity of cytochrome c.

Effect of an Imposed Contact on Secondary Structure in the Denatured State of Yeast Iso-1-cytochrome c.

There is considerable evidence that long-range interactions stabilize residual protein structure under denaturing conditions. However, evaluation of the effect of a specific contact on structure in the denatured state has been difficult. Iso-1-cytochrome c variants with a Lys54 → His mutation form a particularly stable His-heme loop in the denatured state, suggestive of loop-induced residual structure. We have used multidimensional nuclear magnetic resonance methods to assign 1H and 15N backbone amide and 13C backbone and side chain chemical shifts in the denatured state of iso-1-cytochrome c carrying the Lys54 → His mutation in 3 and 6 M guanidine hydrochloride and at both pH 6.4, where the His54-heme loop is formed, and pH 3.6, where the His54-heme loop is broken. Using the secondary structure propensity score, with the 6 M guanidine hydrochloride chemical shift data as a random coil reference state for data collected in 3 M guanidine hydrochloride, we found residual helical structure in the denatured state for the 60s helix and the C-terminal helix, but not in the N-terminal helix in the presence or absence of the His54-heme loop. Non-native helical structure is observed in two regions that form Ω-loops in the native state. There is more residual helical structure in the C-terminal helix at pH 6.4 when the loop is formed. Loop formation also appears to stabilize helical structure near His54, consistent with induction of helical structure observed when His-heme bonds form in heme-peptide model systems. The results are discussed in the context of the folding mechanism of cytochrome c.

Rapid quantification of cardiolipin and DOPC lipid and vesicle concentration.

A novel approach to quantification of cardiolipin and DOPC lipid and vesicle concentration that is rapid and inexpensive is described. Traditional approaches to quantifying vesicle concentration destroy sample and are often time consuming. Using common laboratory equipment and software, lipid vesicles were reliably quantified allowing for immediate use without significant sample loss. Once calibrated, only absorbance measurements with a UV-Vis spectrophotometer are necessary as input into a Matlab program, which calculates the corresponding vesicle and lipid concentration. Fast and accurate concentration determination for preparations of vesicles is essential for analytical titration experiments necessary for protein/vesicle binding curves.

Structure of a mitochondrial cytochrome c conformer competent for peroxidase activity.

At the onset of apoptosis, the peroxidation of cardiolipin at the inner mitochondrial membrane by cytochrome c requires an open coordination site on the heme. We report a 1.45-Å resolution structure of yeast iso-1-cytochrome c with the Met80 heme ligand swung out of the heme crevice and replaced by a water molecule. This conformational change requires modest adjustments to the main chain of the heme crevice loop and is facilitated by a trimethyllysine 72-to-alanine mutation. This mutation also enhances the peroxidase activity of iso-1-cytochrome c. The structure shows a buried water channel capable of facilitating peroxide access to the active site and of moving protons produced during peroxidase activity to the protein surface. Alternate positions of the side chain of Arg38 appear to mediate opening and closing of the buried water channel. In addition, two buried water molecules can adopt alternate positions that change the network of hydrogen bonds in the buried water channel. Taken together, these observations suggest that low and high proton conductivity states may mediate peroxidase function. Comparison of yeast and mammalian cytochrome c sequences, in the context of the steric factors that permit opening of the heme crevice, suggests that higher organisms have evolved to inhibit peroxidase activity, providing a more stringent barrier to the onset of apoptosis.

Lower Protein Stability Does Not Necessarily Increase Local Dynamics.

Overall protein stability is thought to have an important impact on the millisecond time scale dynamics modulating enzyme function. In order to better understand the effects of overall stability on the substructure dynamics of mitochondrial cytochrome c, we test the effect of a destabilizing L85A mutation on the kinetics and equilibrium thermodynamics of the alkaline conformational transition. The alkaline conformational transition replaces the Met80 ligand of the heme with a lysine residue from Ω-loop D, the heme crevice loop, consisting of residues 70-85. Residues 67-87 are the most conserved portion of the sequence of mitochondrial cytochrome c, suggesting that this region is of prime importance for function. Mutations to Ω-loop D affect the stability of the heme crevice directly, modulating the pKapp of the alkaline transition. Two variants of yeast iso-1-cytochrome c, WT*/L85A and WT*/K73H/L85A, were prepared for these studies. Guanidine-HCl unfolding monitored by circular dichroism and pH titrations at 695 nm, respectively, were used to study the thermodynamics of global and local unfolding of these variants. The kinetics of the alkaline transition were measured by pH-jump stopped-flow methods. Gated electron transfer techniques using bis(2,2',2″-terpyridine)cobalt(II) as a reducing reagent were implemented to measure the heme crevice dynamics for the WT*/K73H/L85A variant. Contrary to the expectation that dynamics around the heme crevice would be faster for the less stable WT*/K73H/L85A variant, based on the behavior of psychrophilic versus mesophilic enzymes, they were similar to those for a variant without the L85A mutation. In fact, below pH 7, the dynamics of the WT*/K73H/L85A variant were slower.

Cytochrome c Can Form a Well-Defined Binding Pocket for Hydrocarbons.

Cytochrome c can acquire peroxidase activity when it binds to cardiolipin in mitochondrial membranes. The resulting oxygenation of cardiolipin by cytochrome c provides an early signal for the onset of apoptosis. The structure of this enzyme-substrate complex is a matter of considerable debate. We present three structures at 1.7-2.0 Å resolution of a domain-swapped dimer of yeast iso-1-cytochrome c with the detergents, CYMAL-5, CYMAL-6, and ω-undecylenyl-ß-d-maltopyranoside, bound in a channel that places the hydrocarbon moieties of these detergents next to the heme. The heme is poised for peroxidase activity with water bound in place of Met80, which serves as the axial heme ligand when cytochrome c functions as an electron carrier. The hydroxyl group of Tyr67 sits 3.6-4.0 Å from the nearest carbon of the detergents, positioned to act as a relay in radical abstraction during peroxidase activity. Docking studies with linoleic acid, the most common fatty acid component of cardiolipin, show that C11 of linoleic acid can sit adjacent to Tyr67 and the heme, consistent with the oxygenation pattern observed in lipidomics studies. The well-defined hydrocarbon binding pocket provides atomic resolution evidence for the extended lipid anchorage model for cytochrome c/cardiolipin binding. Dimer dissociation/association kinetics for yeast versus equine cytochrome c indicate that formation of mammalian cytochrome c dimers in vivo would require catalysis. However, the dimer structure shows that only a modest deformation of monomeric cytochrome c would suffice to form the hydrocarbon binding site occupied by these detergents.

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.