Bruce E. Bowler

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.

Residual structure in unfolded proteins.

The denatured state ensemble (DSE) of unfolded proteins, once considered to be well-modeled by an energetically featureless random coil, is now well-known to contain flickering elements of residual structure. The position and nature of DSE residual structure may provide clues toward deciphering the protein folding code. This review focuses on recent advances in our understanding of the nature of DSE collapse under folding conditions, the quantification of the stability of residual structure in the DSE, the determination of the location and types of residues involved in thermodynamically significant residual structure and advances in detection of long-range interactions in the DSE.

Structure of a mitochondrial cytochrome c conformer competent for peroxidase activity

Significance Cytochrome c is essential to two important biochemical pathways, the electron transport chain and the intrinsic pathway of apoptosis. The heme crevice loop, which provides the Met80 ligand to the heme cofactor, is the most highly conserved segment of the cytochrome c sequence. The dynamics of this loop are likely important for both functions. Cytochrome c-mediated peroxidation of cardiolipin in the inner mitochondrial membrane is an early signal in apoptosis. We show that mutation of trimethyllysine 72 to alanine in yeast iso-1-cytochrome c allows formation of a conformer of the protein with Met80 displaced from the heme and enhances peroxidase activity. Thus, this residue is likely an important modulator of the peroxidase function of cytochrome c. 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.

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.

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 ∼1.2 kcal mol−1 compared with ∼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 ∼3 kcal mol−1M−1 versus ∼4.5 kcal mol−1M−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.

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.

Probing the dynamics of a His73-heme alkaline transition in a destabilized variant of yeast iso-1-cytochrome c with conformationally gated electron transfer methods.

The alkaline transition of cytochrome c involves substitution of the Met80 heme ligand of the native state with a lysine ligand from a surface Ω-loop (residues 70 to 85). The standard mechanism for the alkaline transition involves a rapid deprotonation equilibrium followed by the conformational change. However, recent work implicates multiple ionization equilibria and stable intermediates. In previous work, we showed that the kinetics of formation of a His73-heme alkaline conformer of yeast iso-1-cytochrome c requires ionization of the histidine ligand (pK(HL) ~ 6.5). Furthermore, the forward and backward rate constants, k(f) and k(b), respectively, for the conformational change are modulated by two auxiliary ionizations (pK(H1) ~ 5.5, and pK(H2) ~ 9). A possible candidate for pK(H1) is His26, which has a strongly shifted pK(a) in native cytochrome c. Here, we use the AcH73 iso-1-cytochrome c variant, which contains an H26N mutation, to test this hypothesis. pH jump experiments on the AcH73 variant show no change in k(obs) for the His73-heme alkaline transition from pH 5 to 8, suggesting that pK(H1) has disappeared. However, direct measurement of k(f) and k(b) using conformationally gated electron transfer methods shows that the pH independence of k(obs) results from coincidental compensation between the decrease in k(b) due to pK(H1) and the increase in k(f) due to pK(HL). Thus, His26 is not the source of pK(H1). The data also show that the H26N mutation enhances the dynamics of this conformational transition from pH 5 to 10, likely as a result of destabilization of the protein.

Denatured states of low-complexity polypeptide sequences differ dramatically from those of foldable sequences.

How the primary sequence of a protein encodes conformational preferences that operate early in folding to promote efficient formation of the correct native topology is still poorly understood. To address this issue, we have prepared a set of yeast iso-1-cytochrome c variants that contain polyalanine inserts ranging from 6 to 30 residues in length near the N terminus of the protein. We study the thermodynamics and kinetics of His-heme loop formation in the denatured state at 3 and 6 M guanidine-HCl concentration. We find that polyalanine closely approximates a random coil with excluded volume giving scaling exponents, ν3, for equilibrium loop formation of 2.26 ± 0.13 and 1.97 ± 0.04 in 3 and 6 M guanidine-HCl, respectively. The rate of loop breakage initially decreases and then becomes independent of loop size as would be expected for a random coil. Comparison with previously reported data for denatured state His-heme loop formation for iso-1-cytochrome c and Rhodopseudomonas palustris cytochrome c′, shows that foldable sequences deviate significantly from random coil behavior and that the deviation is fold-dependent.

Disruption of a hydrogen bond network in human versus spider monkey cytochrome c affects heme crevice stability.

The hypothesis that the recent rapid evolution of primate cytochromes c, which primarily involves residues in the least stable Ω-loop (Ω-loop C, residues 40-57), stabilizes the heme crevice of cytochrome c relative to other mammals, is tested. To accomplish this goal, we have compared the properties of human and spider monkey cytochrome c and a set of four variants produced in the process of converting human cytochrome c into spider monkey cytochrome c. The global stability of all variants has been measured by guanidine hydrochloride denaturation. The stability of the heme crevice has been assessed with the alkaline conformational transition. Structural insight into the effects of the five amino acid substitutions needed to convert human cytochrome c into spider monkey cytochrome c is provided by a 1.15Å resolution structure of spider monkey cytochrome c. The global stability for all variants is near 9.0kcal/mol at 25°C and pH7, which is higher than that observed for other mammalian cytochromes c. The heme crevice stability is more sensitive to the substitutions required to produce spider monkey cytochrome c with decreases of up to 0.5 units in the apparent pKa of the alkaline conformational transition relative to human cytochrome c. The structure of spider monkey cytochrome c indicates that the Y46F substitution destabilizes the heme crevice by disrupting an extensive hydrogen bond network that connects three surface loops including Ω-loop D (residues 70-85), which contains the Met80 heme ligand.

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.