The low-affinity complex of cytochrome c and its peroxidase

Mechanistic insights into the chemistry of compound I formation in heme peroxidases: quantum chemical investigations of cytochrome c peroxidase

Peroxidases are heme containing enzymes that catalyze peroxide-dependant oxidation of a variety of substrates through forming key ferryl intermediates, compounds I and II. Cytochrome c peroxidase (Ccp1) has served for decades as a chemical model toward understanding the chemical biology of this heme family of enzymes. It is known to feature a distinctive electronic behaviour for its compound I despite significant structural similarity to other peroxidases. A water-assisted mechanism has been proposed over a dry one for the formation of compound I in similar peroxidases. To better identify the viability of these mechanisms, we employed quantum chemistry calculations for the heme pocket of Ccp1 in three different spin states. We provided comparative energetic and structural results for the six possible pathways that suggest the preference of the dry mechanism energetically and structurally. The doublet state is found to be the most preferable spin state for the mechanism to proceed and for the formation of the Cpd I ferryl-intermediate irrespective of the considered dielectric constant used to represent the solvent environment. The nature of the spin state has negligible effects on the calculated structures but great impact on the energetics. Our analysis was also expanded to explain the major contribution of key residues to the peroxidase activity of Ccp1 through exploring the mechanism at various in silico generated Ccp1 variants. Overall, we provide valuable findings toward solving the current ambiguity of the exact mechanism in Ccp1, which could be applied to peroxidases with similar heme pockets.

Discerning the feasibility of a no-water peroxidase mechanism in the doublet spin state irrespective of the environment surrounding the heme pocket.

Enhancing the population of the encounter complex affects protein complex formation efficiency

Optimal charge distribution is considered to be important for efficient formation of protein complexes. Electrostatic interactions guide encounter complex formation that precedes the formation of an active protein complex. However, disturbing the optimized distribution by introduction of extra charged patches on cytochrome c peroxidase does not lead to a reduction in productive encounters with its partner cytochrome c. To test whether a complex with a high population of encounter complex is more easily affected by suboptimal charge distribution, the interactions of cytochrome c mutant R13A with wild‐type cytochrome c peroxidase and a variant with an additional negative patch were studied. The complex of the peroxidase and cytochrome c R13A was reported to have an encounter state population of 80%, compared to 30% for the wild‐type cytochrome c. NMR analysis confirms the dynamic nature of the interaction and demonstrates that the mutant cytochrome c samples the introduced negative patch. Kinetic experiments show that productive complex formation is fivefold to sevenfold slower at moderate and high ionic strength values for cytochrome c R13A but the association rate is not affected by the additional negative patch on cytochrome c peroxidase, showing that the total charge on the protein surface can compensate for less optimal charge distribution. At low ionic strength (44 mm), the association with the mutant cytochrome c reaches the same high rates as found for wild‐type cytochrome c, approaching the diffusion limit.

Erv1 and Cytochrome c Mediate Rapid Electron Transfer via A Collision-Type Interaction

Being essential for oxidative protein folding in the mitochondrial intermembrane space, the mitochondrial disulfide relay relies on the electron transfer (ET) from the sulfhydryl oxidase Erv1 to cytochrome c (Cc). Using solution NMR spectroscopy, we demonstrate that while the yeast Cc-Erv1 system is functionally active, no observable binding of the protein partners takes place. The transient interaction between Erv1 and Cc can be rationalized by molecular modeling, suggesting that a large surface area of Erv1 can sustain a fast ET to Cc via a collision-type mechanism, without the need for a canonical protein complex formation. We suggest that, by preventing the direct ET to molecular oxygen (O₂), the collision-type Cc-Erv1 interaction plays a role in protecting the organism against reactive oxygen species.

The Charge Distribution on a Protein Surface Determines Whether Productive or Futile Encounter Complexes Are Formed

Protein complex formation depends strongly on electrostatic interactions. The distribution of charges on the surface of redox proteins is often optimized by evolution to guide recognition and binding. To test the degree to which the electrostatic interactions between cytochrome c peroxidase (CcP) and cytochrome c (Cc) are optimized, we produced five CcP variants, each with a different charge distribution on the surface. Monte Carlo simulations show that the addition of negative charges attracts Cc to the new patches, and the neutralization of the charges in the regular, stereospecific binding site for Cc abolishes the electrostatic interactions in that region entirely. For CcP variants with the charges in the regular binding site intact, additional negative patches slightly enhance productive complex formation, despite disrupting the optimized charge distribution. Removal of the charges in the regular binding site results in a dramatic decrease in the complex formation rate, even in the presence of highly negative patches elsewhere on the surface. We conclude that additional charge patches can result in either productive or futile encounter complexes, depending on whether negative residues are located also in the regular binding site.

Efficient Encounter Complex Formation and Electron Transfer to Cytochrome c Peroxidase with an Additional, Distant Electrostatic Binding Site

Electrostatic interactions can strongly increase the efficiency of protein complex formation. The charge distribution in redox proteins is often optimized to steer a redox partner to the electron transfer active binding site. To test whether the optimized distribution is more important than the strength of the electrostatic interactions, an additional negative patch was introduced on the surface of cytochrome c peroxidase, away from the stereospecific binding site, and its effect on the encounter complex as well as the rate of complex formation was determined. Monte Carlo simulations and paramagnetic relaxation enhancement NMR experiments indicate that the partner, cytochrome c, interacts with the new patch. Unexpectedly, the rate of the active complex formation was not reduced, but rather slightly increased. The findings support the idea that for efficient protein complex formation the strength of the electrostatic interaction is more critical than an optimized charge distribution.

The Transient Complex of Cytochrome c and Cytochrome c Peroxidase: Insights into the Encounter Complex from Multifrequency EPR and NMR Spectroscopy

We present a novel approach to study transient protein-protein complexes with standard, 9 GHz, and high-field, 95 GHz, electron paramagnetic resonance (EPR) and paramagnetic NMR at ambient temperatures and in solution. We apply it to the complex of yeast mitochondrial iso-1-cytochrome c (Cc) with cytochrome c peroxidase (CcP) with the spin label [1-oxyl-2,2,5,5-tetramethyl-Δ3-pyrroline-3-methyl)-methanethiosulfonate] attached at position 81 of Cc (SL-Cc). A dissociation constant KD of 20±4×10-6  M (EPR and NMR) and an equal amount of stereo-specific and encounter complex (NMR) are found. The EPR spectrum of the fully bound complex reveals that the encounter complex has a significant population (60 %) that shares important features, such as the Cc-interaction surface, with the stereo-specific complex.

Tuning Radical Relay Residues by Proton Management Rescues Protein Electron Hopping

Transient tyrosine and tryptophan radicals play key roles in the electron transfer (ET) reactions of photosystem (PS) II, ribonucleotide reductase (RNR), photolyase, and many other proteins. However, Tyr and Trp are not functionally interchangeable, and the factors controlling their reactivity are often unclear. Cytochrome c peroxidase (CcP) employs a Trp191^•+ radical to oxidize reduced cytochrome c (Cc). Although a Tyr191 replacement also forms a stable radical, it does not support rapid ET from Cc. Here we probe the redox properties of CcP Y191 by non-natural amino acid substitution, altering the ET driving force and manipulating the protic environment of Y191. Higher potential fluorotyrosine residues increase ET rates marginally, but only addition of a hydrogen bond donor to Tyr191^• (via Leu232His or Glu) substantially alters activity by increasing the ET rate by nearly 30-fold. ESR and ESEEM spectroscopies, crystallography, and pH-dependent ET kinetics provide strong evidence for hydrogen bond formation to Y191^• by His232/Glu232. Rate measurements and rapid freeze quench ESR spectroscopy further reveal differences in radical propagation and Cc oxidation that support an increased Y191^• formal potential of ∼200 mV in the presence of E232. Hence, Y191 inactivity results from a potential drop owing to Y191^•+ deprotonation. Incorporation of a well-positioned base to accept and donate back a hydrogen bond upshifts the Tyr^• potential into a range where it can effectively oxidize Cc. These findings have implications for the Y_Z/Y_D radicals of PS II, hole-hopping in RNR and cryptochrome, and engineering proteins for long-range ET reactions.

Constraints on the Radical Cation Center of Cytochrome c Peroxidase for Electron Transfer from Cytochrome c

The tryptophan 191 cation radical of cytochrome c peroxidase (CcP) compound I (Cpd I) mediates long-range electron transfer (ET) to cytochrome c (Cc). Here we test the effects of chemical substitution at position 191. CcP W191Y forms a stable tyrosyl radical upon reaction with peroxide and produces spectral properties similar to those of Cpd I but has low reactivity toward reduced Cc. CcP W191G and W191F variants also have low activity, as do redox ligands that bind within the W191G cavity. Crystal structures of complexes between Cc and CcP W191X (X = Y, F, or G), as well as W191G with four bound ligands reveal similar 1:1 association modes and heme pocket conformations. The ligands display structural disorder in the pocket and do not hydrogen bond to Asp235, as does Trp191. Well-ordered Tyr191 directs its hydroxyl group toward the porphyrin ring, with no basic residue in the range of interaction. CcP W191X (X = Y, F, or G) variants substituted with zinc-porphyrin (ZnP) undergo photoinduced ET with Cc(III). Their slow charge recombination kinetics that result from loss of the radical center allow resolution of difference spectra for the charge-separated state [ZnP(+), Cc(II)]. The change from a phenyl moiety at position 191 in W191F to a water-filled cavity in W191G produces effects on ET rates much weaker than the effects of the change from Trp to Phe. Low net reactivity of W191Y toward Cc(II) derives either from the inability of ZnP(+) or the Fe-CcP ferryl to oxidize Tyr or from the low potential of the resulting neutral Tyr radical.

Contribution of Electrostatics to the Kinetics of Electron Transfer from NADH-Cytochrome b5 Reductase to Fe(III)-Cytochrome b5

Brownian dynamics (BD) simulations provide here a theoretical atomic-level treatment of the reduction of human ferric cytochrome b5 (cyt b5) by NADH-cytochrome b5 reductaste (cyt b5r) and several of its mutants. BD is used to calculate the second-order rate constant of electron transfer (ET) between the proteins for direct correlation with experiments. Interestingly, the inclusion of electrostatic forces dramatically increases the reaction rate of the native proteins despite the overall negative charge of both proteins. The role played by electrostatic charge distribution in stabilizing the ET complexes and the role of mutations of several amino acid residues in stabilizing or destabilizing the complexes are analyzed. The complex with the shortest ET reaction distance (d = 6.58 Å) from rigid body BD is further subjected to 1 ns of molecular dynamics (MD) in a periodic box of TIP3P water to produce a more stable complex allowed by flexibility and with a shorter average reaction distance d = 6.02 Å. We predict a docking model in which the following ion-ion interactions are dominant (cyt b5r/cyt b5): Lys162-Heme O1D/Lys163-Asp64/Arg91-Heme O1A/Lys125-Asp70.

Structure and Function of Transient Encounters of Redox Proteins

Many biomolecular interactions proceed via lowly populated, transient intermediates. Believed to facilitate formation of a productive complex, these short-lived species are inaccessible to conventional biophysical and structural techniques and, until recently, could only be studied by theoretical simulations. Recent development of experimental approaches sensitive to the presence of minor species--in particular paramagnetic relaxation enhancement (PRE) NMR spectroscopy--has enabled direct visualization and detailed characterization of such lowly populated states. Collectively referred to as an encounter complex, the binding intermediates are particularly important in transient protein interactions, such as those orchestrating signaling cascades or energy-generating electron transfer (ET) chains. Here I discuss encounter complexes of redox proteins mediating biological ET reactions, which are essential for many vital cellular activities including oxidative phosphorylation and photosynthesis. In particular, this Account focuses on the complex of cytochrome c (Cc) and cytochrome c peroxidase (CcP), which is a paradigm of biomolecular ET and an attractive system for studying protein binding and enzymatic catalysis. The Cc-CcP complex formation proceeds via an encounter state, consisting of multiple protein-protein orientations sampled in the search of the dominant, functionally active bound form and exhibiting a broad spatial distribution, in striking agreement with earlier theoretical simulations. At low ionic strength, CcP binds another Cc molecule to form a weak ternary complex, initially inferred from kinetics experiments and postulated to account for the measured ET activity. Despite strenuous efforts, the ternary complex could not be observed directly and remained eagerly sought for the past two decades. Very recently, we have solved its structure in solution and shown that it consists of two binding forms: the dominant, ET-inactive geometry and an ensemble of lowly populated species with short separations between Cc and CcP cofactors, which summarily account for the measured ET rate. Unlike most protein complexes, which require accurate alignment of the binding surfaces in a single, well-defined orientation to carry out their function, redox proteins can form multiple productive complexes. As fast ET will occur any time the redox centers of the binding partners are close enough to ensure efficient electron tunneling across the interface, many protein-protein orientations are expected to be ET active. The present analysis confirms that the low-occupancy states can support the functional ET activity and contribute to the stability of redox protein complexes. As illustrated here, boundaries between the dominant and the encounter forms become blurred for many dynamic ET systems, which are more aptly described by ensembles of functionally and structurally heterogeneous bound forms.