Role of Hydrogen Bond Networks and Dynamics in Positive and Negative Cooperative Stabilization of a Protein

Probing the versatility of cytochrome c by spectroscopic means: A Laudatio on resonance Raman spectroscopy

Over the last 50 years resonance Raman spectroscopy has become an invaluable tool for the exploration of chromophores in biological macromolecules. Among them, heme proteins and metal complexes have attracted considerable attention. This interest results from the fact that resonance Raman spectroscopy probes the vibrational dynamics of these chromophores without direct interference from the surrounding. However, the indirect influence via through-bond and through-space chromophore-protein interactions can be conveniently probed and analyzed. This review article illustrates this point by focusing on class 1 cytochrome c, a comparatively simple heme protein generally known as electron carrier in mitochondria. The article demonstrates how through selective excitation of resonance Raman active modes information about the ligation, the redox state and the spin state of the heme iron can be obtained from band positions in the Raman spectra. The investigation of intensities and depolarization ratios emerged as tools for the analysis of in-plane and out-of-plane deformations of the heme macrocycle. The article further shows how resonance Raman spectroscopy was used to characterize partially unfolded states of oxidized cytochrome c. Finally, it describes its use for exploring structural changes due to the protein's binding to anionic surfaces like cardiolipin containing membranes.

C-Graphs Tool with Graphical User Interface to Dissect Conserved Hydrogen-Bond Networks: Applications to Visual Rhodopsins

Dynamic hydrogen-bond networks provide proteins with structural plasticity required to translate signals such as ligand binding into a cellular response or to transport ions and larger solutes across membranes and, thus, are of central interest to understand protein reaction mechanisms. Here, we present C-Graphs, an efficient tool with graphical user interface that analyzes data sets of static protein structures or of independent numerical simulations to identify conserved, vs unique, hydrogen bonds and hydrogen-bond networks. For static structures, which may belong to the same protein or to proteins with different sequences, C-Graphs uses a clustering algorithm to identify sites of the hydrogen-bond network where waters are conserved among the structures. Using C-Graphs, we identify an internal protein-water hydrogen-bond network common to static structures of visual rhodopsins and adenosine A2A G protein-coupled receptors (GPCRs). Molecular dynamics simulations of a visual rhodopsin indicate that the conserved hydrogen-bond network from static structure can recruit dynamic hydrogen bonds and extend throughout most of the receptor. We release with this work the code for C-Graphs and its graphical user interface.

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.

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.

A Heme Propionate Staples the Structure of Cytochrome c for Methionine Ligation to the Heme Iron

Ligand-switch reactions at the heme iron are common in biological systems, but their mechanisms and the features of the polypeptide fold that support dual ligation are not well understood. In cytochrome c (cyt c), two low-stability loops (Ω-loop C and Ω-loop D) are connected by the heme propionate HP6. At alkaline pH, the native Met80 ligand from Ω-loop D switches to a Lys residue from the same loop. Deprotonation of an as yet unknown group triggers the alkaline transition. We have created the two cyt c variants T49V/K79G and T78V/K79G with altered connections of these two loops to HP6. Electronic absorption, NMR, and EPR studies demonstrate that at pH 7.4 ferric forms of these variants are Lys-ligated, whereas ferrous forms maintain the native Met80 ligation. Measurements of protein stability, cyclic voltammetry, pH-jump and gated electron-transfer kinetics have revealed that these Thr to Val substitutions greatly affect the alkaline transition in both ferric and ferrous proteins. The substitutions modify the stability of the Met-ligated species and reduction potentials of the heme iron. The kinetics of ligand-switch processes are also altered, and analyses of these effects implicate redox-dependent differences in metal-ligand interactions and the role of the protein dynamics, including cross-talk between the two Ω-loops. With the two destabilized variants, it is possible to map energy levels for the Met- and Lys-ligated species in both ferric and ferrous proteins and assess the role of the protein scaffold in redox-dependent preferences for these two ligands. The estimated shift in the heme iron reduction potential upon deprotonation of the "trigger" group is consistent with those associated with deprotonation of an HP, suggesting that HP6, on its own or as a part of a hydrogen-bonded cluster, is a likely "trigger" for the Met to Lys ligand switch.

Rapid Sampling of Hydrogen Bond Networks for Computational Protein Design

Hydrogen bond networks play a critical role in determining the stability and specificity of biomolecular complexes, and the ability to design such networks is important for engineering novel structures, interactions, and enzymes. One key feature of hydrogen bond networks that makes them difficult to rationally engineer is that they are highly cooperative and are not energetically favorable until the hydrogen bonding potential has been satisfied for all buried polar groups in the network. Existing computational methods for protein design are ill-equipped for creating these highly cooperative networks because they rely on energy functions and sampling strategies that are focused on pairwise interactions. To enable the design of complex hydrogen bond networks, we have developed a new sampling protocol in the molecular modeling program Rosetta that explicitly searches for sets of amino acid mutations that can form self-contained hydrogen bond networks. For a given set of designable residues, the protocol often identifies many alternative sets of mutations/networks, and we show that it can readily be applied to large sets of residues at protein-protein interfaces or in the interior of proteins. The protocol builds on a recently developed method in Rosetta for designing hydrogen bond networks that has been experimentally validated for small symmetric systems but was not extensible to many larger protein structures and complexes. The sampling protocol we describe here not only recapitulates previously validated designs with performance improvements but also yields viable hydrogen bond networks for cases where the previous method fails, such as the design of large, asymmetric interfaces relevant to engineering protein-based therapeutics.

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.

Humanlike substitutions to Ω-loop D of yeast iso-1-cytochrome c only modestly affect dynamics and peroxidase activity

Structural studies of yeast iso-1-cytochrome c (L.J. McClelland, T.-C. Mou, M.E. Jeakins-Cooley, S.R. Sprang, B.E. Bowler, Proc. Natl. Acad. Sci. U.S.A. 111 (2014) 6648-6653) show that modest movement of Ω-loop D (residues 70-85, average RMSD versus the native structure: 0.81 Å) permits loss of Met80-heme ligation creating an available coordination site to catalyze the peroxidase activity mediated by cytochrome c early in apoptosis. However, Ala81 and Gly83 move significantly (RMSDs of 2.18 and 1.26 Å, respectively). Ala81 and Gly83 evolve to Ile and Val, respectively, in human cytochrome c and peroxidase activity decreases 25-fold relative to the yeast protein at pH 7. To test the hypothesis that these residues evolved to restrict the peroxidase activity of cytochrome c, A81I and G83V variants of yeast iso-1-cytochrome c were prepared. For both variants, the apparent pK_a of the alkaline transition increases by 0.2 to 0.3 relative to the wild type (WT) protein and the rate of opening the heme crevice is slowed. The cooperativity of acid unfolding is decreased for the G83V variant. At pH 7 and 8, the catalytic rate constant, k_cat, for the peroxidase activity of both variants decreases relative to WT, consistent with the effects on alkaline isomerization. Below pH 7, the loss in the cooperativity of acid unfolding causes k_cat for peroxidase activity to increase for the G83V variant relative to WT. Neither variant decreases k_cat to the level of the human protein, indicating that other residues also contribute to the low peroxidase activity of human cytochrome c.

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.

Remote Perturbations in Tertiary Contacts Trigger Ligation of Lysine to the Heme Iron in Cytochrome c

Perturbations in protein structure define the mechanism of allosteric regulation and biological information transfer. In cytochrome c (cyt c), ligation of Met80 to the heme iron is critical for the protein's electron-transfer (ET) function in oxidative phosphorylation and for suppressing its peroxidase activity in apoptosis. The hard base Lys is a better match for the hard ferric iron than the soft base Met is, suggesting the key role of the protein scaffold in favoring Met ligation. To probe the role of the protein structure in the maintenance of Met ligation, mutations T49V and Y67R/M80A were designed to disrupt hydrogen bonding and packing of the heme coordination loop, respectively. Electronic absorption, nuclear magnetic resonance, and electron paramagnetic resonance spectra reveal that ferric forms of both variants are Lys-ligated at neutral pH. A minor change in the tertiary contacts in T49V, away from the heme coordination loop, appears to be sufficient to execute a change in ligation, suggesting a cross-talk between the different regions of the protein structure and a possibility of built-in conformational switches in cyt c. Analyses of thermodynamic stability, kinetics of Lys binding and dissociation, and the pH-dependent changes in ligation provide a detailed characterization of the Lys coordination in these variants and relate these properties to the extent of structural perturbations. The findings emphasize the importance of the hydrogen-bonding network in controlling ligation of the native Met80 to the heme iron.

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.

Structure-function relationships in human cytochrome c: The role of tyrosine 67

Spectroscopic and functional properties of human cytochrome c and its Tyr67 residue mutants (i.e., Tyr67His and Tyr67Arg) have been investigated. In the case of the Tyr67His mutant, we have observed only a very limited structural alteration of the heme pocket and of the Ω-loop involving, among others, the residue Met80 and its bond with the heme iron. Conversely, in the Tyr67Arg mutant the Fe-Met80 bond is cleaved; consequently, a much more extensive structural alteration of the Ω-loop can be envisaged. The structural, and thus the functional modifications, of the Tyr67Arg mutant are present in both the ferric [Fe(III)] and the ferrous [Fe(II)] forms, indicating that the structural changes are independent of the heme iron oxidation state, depending instead on the type of substituting residue. Furthermore, a significant peroxidase activity is evident for the Tyr67Arg mutant, highlighting the role of Arg as a basic, positively charged residue at pH7.0, located in the heme distal pocket, which may act as an acid to cleave the O-O bond in H2O2. As a whole, our results indicate that a delicate equilibrium is associated with the spatial arrangement of the Ω-loop. Clearly, Arg, but not His, is able to stabilize and polarize the negative charge on the Fe(III)-OOH complex during the formation of Compound I, with important consequences on cytochrome peroxidation activity and its role in the apoptotic process, which is somewhat different in yeast and mammals.

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.

A Compact Structure of Cytochrome c Trapped in a Lysine-Ligated State: Loop Refolding and Functional Implications of a Conformational Switch

It has been suggested that the alkaline form of cytochrome c (cyt c) regulates function of this protein as an electron carrier in oxidative phosphorylation and as a peroxidase that reacts with cardiolipin (CL) during apoptosis. In this form, Met80, the native ligand to the heme iron, is replaced by a Lys. While it has become clear that the structure of cyt c changes, the extent and sequence of conformational rearrangements associated with this ligand replacement remain a subject of debate. Herein we report a high-resolution crystal structure of a Lys73-ligated cyt c conformation that reveals intricate change in the heme environment upon this switch in the heme iron ligation. The structure is surprisingly compact, and the heme coordination loop refolds into a β-hairpin with a turn formed by the highly conserved residues Pro76 and Gly77. Repositioning of residue 78 modifies the intraprotein hydrogen-bonding network and, together with adjustments of residues 52 and 74, increases the volume of the heme pocket to allow for insertion of one of the CL acyl moieties next to Asn52. Derivatization of Cys78 with maleimide creates a solution mimic of the Lys-ligated cyt c that has enhanced peroxidase activity, adding support for a role of the Lys-ligated cyt c in the apoptotic mechanism. Experiments with the heme peptide microperoxidase-8 and engineered model proteins provide a thermodynamic rationale for the switch to Lys ligation upon perturbations in the protein scaffold.

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.

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.

Dynamics of the His79-heme alkaline transition of yeast iso-1-cytochrome c probed by conformationally gated electron transfer with Co(II)bis(terpyridine)

Alkaline conformers of cytochrome c may be involved in both its electron transport and apoptotic functions. We use cobalt(II)bis(terpyridine), Co(terpy)2(2+), as a reagent for conformationally gated electron-transfer (gated ET) experiments to study the alkaline conformational transition of K79H variants of yeast iso-1-cytochrome c expressed in Escherichia coli , WT*K79H, with alanine at position 72 and Saccharomyces cerevisiae , yK79H, with trimethyllysine (Tml) at position 72. Co(terpy)2(2+) is well-suited to the 100 ms to 1 s time scale of the His79-mediated alkaline conformational transition of these variants. Reduction of the His79-heme alkaline conformer by Co(terpy)2(2+) occurs primarily by gated ET, which involves conversion to the native state followed by reduction, with a small fraction of the His79-heme alkaline conformer directly reduced by Co(terpy)2(2+). The gated ET experiments show that the mechanism of formation of the His79-heme alkaline conformer involves only two ionizable groups. In previous work, we showed that the mechanism of the His73-mediated alkaline conformational transition requires three ionizable groups. Thus, the mechanism of heme crevice opening depends upon the position of the ligand mediating the process. The microscopic rate constants provided by gated ET studies show that mutation of Tml72 (yK79H variant) in the heme crevice loop to Ala72 (WT*K79H variant) affects the dynamics of heme crevice opening through a small destabilization of both the native conformer and the transition state relative to the His79-heme alkaline conformer. Previous pH jump data had indicated that the Tml72→Ala mutation primarily stabilized the transition state for the His79-mediated alkaline conformational transition.

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.

Chemical double-mutant cycles: dissecting non-covalent interactions

Thermodynamic double-mutant cycles and triple-mutant boxes are widely employed for the experimental quantification of non-covalent interactions and cooperative effects in proteins. This review describes the application of these powerful methodologies to the study of non-covalent interactions in synthetic systems.

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.

Order of steps in the cytochrome C folding pathway: evidence for a sequential stabilization mechanism

Previous work used hydrogen exchange (HX) experiments in kinetic and equilibrium modes to study the reversible unfolding and refolding of cytochrome c (Cyt c) under native conditions. Accumulated results now show that Cyt c is composed of five individually cooperative folding units, called foldons, which unfold and refold as concerted units in a stepwise pathway sequence. The first three steps of the folding pathway are linear and sequential. The ordering of the last two steps has been unclear because the fast HX of the amino acid residues in these foldons has made measurement difficult. New HX experiments done under slower exchange conditions show that the final two foldons do not unfold and refold in an obligatory sequence. They unfold separately and neither unfolding obligately contains the other, as indicated by their similar unfolding surface exposure and the specific effects of destabilizing and stabilizing mutations, pH change, and oxidation state. These results taken together support a sequential stabilization mechanism in which folding occurs in the native context with prior native-like structure serving to template the stepwise formation of subsequent native-like foldon units. Where the native structure of Cyt c requires sequential folding, in the first three steps, this is found. Where structural determination is ambiguous, in the final two steps, alternative parallel folding is found.