Electrostatic analysis of the interaction of cytochrome c with native and dimethyl ester heme substituted cytochrome b5

Contributions to cytochrome c inner- and outer-sphere reorganization energy

Cytochromes c are small water-soluble proteins that catalyze electron transfer in metabolism and energy conversion processes. Hydrogenobacter thermophilus cytochrome c₅₅₂ presents a curious case in displaying fluxionality of its heme axial methionine ligand; this behavior is altered by single point mutation of the Q64 residue to N64 or V64, which fixes the ligand in a single configuration. The reorganization energy (λ) of these cytochrome c₅₅₂ variants is experimentally determined using a combination of rotating disc electrochemistry, chronoamperometry and cyclic voltammetry. The differences between the λ determined from these complementary techniques helps to deconvolute the contribution of the active site and its immediate environment to the overall λ (λ_Total). The experimentally determined λ values in conjunction with DFT calculations indicate that the differences in λ among the protein variants are mainly due to the differences in contributions from the protein environment and not just inner-sphere λ. DFT calculations indicate that the position of residue 64, responsible for the orientation of the axial methionine, determines the geometric relaxation of the redox active molecular orbital (RAMO). The orientation of the RAMO with respect to the heme is key to determining electron transfer coupling (H_AB) which results in higher ET rates in the wild-type protein relative to the Q64V mutant despite a 150 mV higher λ_Total in the former.

Efficient delocalization of the redox-active molecular orbital (RAMO) in HtWT results in an increase in H_AB value which in turn accelerates the electron transfer (ET) rate in spite of the higher reorganization energy (λ) than the HtQ64V mutant.

Binding Mechanisms of Electron Transport Proteins with Cyanobacterial Photosystem I: An Integrated Computational and Experimental Model

The stromal domain (PsaC, D, and E) of photosystem I (PSI) in cyanobacteria accepts electrons from PsaA and PsaB of photosystem I (PSI). These electrons are then used in the reduction of transiently bound ferredoxin (Fd) or flavodoxin. Experimental X-ray and NMR structures are known for all of these protein partners separately, yet to date, there is no known experimental structure of the PSI/Fd complexes in the published literature. Computational models of Fd docked with the stromal domain of cyanobacterial PSI were assembled here starting from X-ray and NMR structures of PSI and Fd. Predicted models of specific regions of protein-protein interactions were built on the basis of energetic frustration, residue conservation and evolutionary importance, as well as from experimental site-directed mutagenesis and cross-linking studies. Microsecond time-scale molecular dynamics simulations of the PSI/Fd complexes suggest, rather than a single complex structure, the possible existence of multiple transient complexes of Fd bound to PSI.

Histidine-Lysine Axial Ligand Switching in a Hemoglobin: A Role for Heme Propionates

The hemoglobin of Synechococcus sp. PCC 7002, GlbN, is a monomeric group I truncated protein (TrHb1) that coordinates the heme iron with two histidine ligands at neutral pH. One of these is the distal histidine (His46), a residue that can be displaced by dioxygen and other small molecules. Here, we show with mutagenesis, electronic absorption spectroscopy, and nuclear magnetic resonance (NMR) spectroscopy that at high pH and exclusively in the ferrous state, Lys42 competes with His46 for the iron coordination site. When b heme is originally present, the population of the lysine-bound species remains too small for detailed characterization; however, the population can be increased significantly by using dimethyl-esterified heme. Electronic absorption and NMR spectroscopies showed that the reversible ligand switching process occurs with an apparent pK_a of 9.3 and a Lys-ligated population of ∼60% at the basic pH limit in the modified holoprotein. The switching rate, which is slow on the chemical shift time scale, was estimated to be 20-30 s^-1 by NMR exchange spectroscopy. Lys42-His46 competition and attendant conformational rearrangement appeared to be related to weakened bis-histidine ligation and enhanced backbone dynamics in the ferrous protein. The pH- and redox-dependent ligand exchange process observed in GlbN illustrates the structural plasticity allowed by the TrHb1 fold and demonstrates the importance of electrostatic interactions at the heme periphery for achieving axial ligand selection. An analogy is drawn to the alkaline transition of cytochrome c, in which Lys-Met competition is detected at alkaline pH, but, in contrast to GlbN, in the ferric state only.

Stabilization of cytochrome b5 by a conserved tyrosine in the secondary sphere of heme active site: A spectroscopic and computational study

Heme proteins perform a large array of biological functions, with the heme group bound non-covalently or covalently. To probe the stabilization role of conserved tyrosine residue in the secondary sphere of heme site in heme proteins, we herein used cytochrome b₅ (Cyt b₅) as a model protein, and mutated Tyr30 to Phe or His by removal of Tyr30 associated H-bond network and hydrophobic interaction. We performed thermal-induced unfolding studies for the two mutants, Y30F Cyt b₅ and Y30H Cyt b₅, as monitored by both UV-Vis and CD spectroscopy, as well as heme transfer studies from these proteins to apo-myoglobin, with wild-type Cyt b₅ under the same conditions for comparison. The reduced stability of both mutants indicates that both the H-bonding and hydrophobic interactions associated with Tyr30 contribute to the protein stability. Moreover, we performed molecular modeling studies, which revealed that the hydrophobic interaction in the local region of Y30F Cyt b₅ was well-remained, whereas Y30H Cyt b₅ formed an H-bond network. These observations suggest that the conserved Tyr30 in Cyt b₅ is not replaceable due to the presence of both the H-bond network and hydrophobic interaction in the secondary sphere of the heme active site. As demonstrated here for Cyt b₅, it may be of practical importance for design of artificial heme proteins by engineering a Tyr in the secondary sphere with improved properties and functions.

STEREODYNAMICS, MASS EFFECT AND DOCKING FROM O(3P) + HD → OH + D AT Ecol = 0.4–1.0eV AND O(3P) + HD → OD + H AT Ecol = 0.5–1.0eV

Quasiclassical Trajectory (QCT) calculation for O(3P) + HD → OH + D and O(3P) + HD → OD + H at E col = 0.4–1.0 eV and 0.5–1.0 eV, respectively, on the lowest PES 1 3A″ of Kuppermann et al. has been done. Distribution p(ϑr) of azimuthal angles between the relative velocity k of the reactants and rotational angular momentum vector j′ of either OH or OD , p(φr) of polar as well as dihedral angles correlating k - k′ -j′, p(ϑr, φr), and PDDCS dependent upon the scattering angle ϑt of either OH , or OD between k and k′ of the reactants and products, respectively, are presented and discussed. The stereodynamics and isotopic mass effects at the smallest possible collision energies 0.4 eV and 0.5 eV for OH and OD , respectively, are significantly different. The significant mass effect with quotient 1/2 of H/D, at the corresponding collision threshold may be applied for the investigation of docking mechanism, drug modification and delivery.

Phospholipid membranes affect tertiary structure of the soluble cytochrome b5 heme-binding domain

The influence of charged phospholipid membranes on the conformational state of the water-soluble fragment of cytochrome b5 has been investigated by a variety of techniques at neutral pH. The results of this work provide the first evidence that aqueous solutions with high phospholipid/protein molar ratios (pH 7.2) induce the cytochrome to undergo a structural transition from the native conformation to an intermediate state with molten-globule like properties that occur in the presence of an artificial membrane surface and that leads to binding of the protein to the membrane. At other phospholipid/protein ratios, equilibrium was observed between cytochrome free in solution and cytochrome bound to the surface of vesicles. Inhibition of protein binding to the vesicles with increasing ionic strength indicated for the most part an electrostatic contribution to the stability of cytochrome b5-vesicle interactions at pH 7.2. The possible physiological role of membrane-induced conformational change in the structure of cytochrome b5 upon the interaction with its redox partners is discussed.

Scoring a diverse set of high-quality docked conformations: A metascore based on electrostatic and desolvation interactions

Predicting protein-protein interactions involves sampling and scoring docked conformations. Barring some large structural rearrangement, rapidly sampling the space of docked conformations is now a real possibility, and the limiting step for the successful prediction of protein interactions is the scoring function used to reduce the space of conformations from billions to a few, and eventually one high affinity complex. An atomic level free-energy scoring function that estimates in units of kcal/mol both electrostatic and desolvation interactions (plus van der Waals if appropriate) of protein-protein docked conformations is used to rerank the blind predictions (860 in total) submitted for six targets to the community-wide Critical Assessment of PRediction of Interactions (CAPRI; http://capri.ebi.ac.uk). We found that native-like models often have varying intermolecular contacts and atom clashes, making unlikely that one can construct a universal function that would rank all these models as native-like. Nevertheless, our scoring function is able to consistently identify the native-like complexes as those with the lowest free energy for the individual models of 16 (out of 17) human predictors for five of the targets, while at the same time the modelers failed to do so in more than half of the cases. The scoring of high-quality models developed by a wide variety of methods and force fields confirms that electrostatic and desolvation forces are the dominant interactions determining the bound structure. The CAPRI experiment has shown that modelers can predict valuable models of protein-protein complexes, and improvements in scoring functions should soon solve the docking problem for complexes whose backbones do not change much upon binding. A scoring server and programs are available at http://structure.pitt.edu.

Characterization and calculation of a cytochrome c-cytochrome b5 complex using NMR data

To identify the binding site for bovine cytochrome b(5) (cyt b(5)) on horse cytochrome c (cyt c), cross-saturation transfer NMR experiments were performed with (2)H- and (15)N-enriched cyt c and unlabeled cyt b(5). In addition, chemical shift changes of the cyt c backbone amide and side chain methyl resonances were monitored as a function of cyt b(5) concentration. The chemical shift changes indicate that the complex is in fast exchange, and are consistent with a 1:1 stoichiometry. A K(a) of (4 +/- 3) x 10(5) M(-)(1) was obtained with a lower limit of 855 s(-)(1) for the dissociation rate of the complex. Mapping of the chemical shift variations and intensity changes upon cross-saturation NMR experiments in the complex reveals a single, contiguous interaction interface on cyt c. Using NMR data as constraints, a protein docking program was used to calculate two low-energy model complex clusters. Independent calculations of the effect of the cyt b(5) heme ring current-induced magnetic dipole on cyt c were used to discriminate between the different models. The interaction surface of horse cyt c in the current experimentally constrained model of the cyt c-cyt b(5) complex is similar but not identical to the interface predicted in yeast cyt c by Brownian dynamics and docking calculations. The occurrence of different amino acids at the protein-protein interface and the dissimilar assumptions employed in the calculations can largely account for the nonidentical interfaces.

The orientations of cytochrome c in the highly dynamic complex with cytochrome b5 visualized by NMR and docking using HADDOCK

The interaction of bovine microsomal ferricytochrome b5 with yeast iso-1-ferri and ferrocytochrome c has been investigated using heteronuclear NMR techniques. Chemical-shift perturbations for 1H and 15N nuclei of both cytochromes, arising from the interactions with the unlabeled partner proteins, were used for mapping the interacting surfaces on both proteins. The similarity of the binding shifts observed for oxidized and reduced cytochrome c indicates that the complex formation is not influenced by the oxidation state of the cytochrome c. Protein-protein docking simulations have been performed for the binary cytochrome b5-cytochrome c and ternary (cytochrome b5)-(cytochrome c)2 complexes using a novel HADDOCK approach. The docking procedure, which makes use of the experimental data to drive the docking, identified a range of orientations assumed by the proteins in the complex. It is demonstrated that cytochrome c uses a confined surface patch for interaction with a much more extensive surface area of cytochrome b5. Taken together, the experimental data suggest the presence of a dynamic ensemble of conformations assumed by the proteins in the complex.

A further investigation of the cytochrome b5-cytochrome c complex

The interaction of reduced rabbit cytochrome b(5) with reduced yeast iso-1 cytochrome c has been studied through the analysis of (1)H-(15)N HSQC spectra, of (15)N longitudinal ( R(1)) and transverse ( R(2)) relaxation rates, and of the solvent exchange rates of protein backbone amides. For the first time, the adduct has been investigated also from the cytochrome c side. The analysis of the NMR data was integrated with docking calculations. The result is that cytochrome b(5) has two negative patches capable of interacting with a single positive surface area of cytochrome c. At low protein concentrations and in equimolar mixture, two different 1:1 adducts are formed. At high concentration and/or with excess cytochrome c, a 2:1 adduct is formed. All the species are in fast exchange on the scale of differences in chemical shift. By comparison with literature data, it appears that the structure of one 1:1 adduct changes with the origin or primary sequence of cytochrome b(5).

X-ray crystallography, CD and kinetic studies revealed the essence of the abnormal behaviors of the cytochrome b5 Phe35-->Tyr mutant

Conserved phenylalanine 35 is one of the hydrophobic patch residues on the surface of cytochrome b5 (cyt b5). This patch is partially exposed on the surface of cyt b5 while its buried face is in direct van der Waals' contact with heme b. Residues Phe35 and Phe/Tyr74 also form an aromatic channel with His39, which is one of the axial ligands of heme b. By site-directed mutagenesis we have produced three mutants of cyt b5: Phe35-->Tyr, Phe35-->Leu, and Phe35-->His. We found that of these three mutants, the Phe35-->Tyr mutant displays abnormal properties. The redox potential of the Phe35-->Tyr mutant is 66 mV more negative than that of the wild-type cyt b5 and the oxidized Phe35-->Tyr mutant is more stable towards thermal and chemical denaturation than wild-type cyt b5. In this study we studied the most interesting mutant, Phe35-->Tyr, by X-ray crystallography, thermal denaturation, CD and kinetic studies of heme dissociation to explore the origin of its unusual behaviors. Analysis of crystal structure of the Phe35-->Tyr mutant shows that the overall structure of the mutant is basically the same as that of the wild-type protein. However, the introduction of a hydroxyl group in the heme pocket, and the increased van der Waals' and electrostatic interactions between the side chain of Tyr35 and the heme probably result in enhancement of stability of the Phe35-->Tyr mutant. The kinetic difference of the heme trapped by the heme pocket also supports this conclusion. The detailed conformational changes of the proteins in response to heat have been studied by CD for the first time, revealing the existence of the folding intermediate.

Dimethyl propionate ester heme-containing cytochrome b5: structure and stability

A derivative of rat microsomal cytochrome b5, obtained by substitution of the native heme moiety with protoporphyrin IX dimethyl ester, has been characterized by 1H and 15N NMR spectroscopy. Besides the two usual A and B forms, which depend on the orientation of the heme in the prostethic group cavity, two other minor forms have been detected which presumably indicate different conformations of the vinyl side chains. The shifts of the heme methyls, as well as the directions of the rhombic axes of the magnetic susceptibility tensor, indicate a small difference in the orientation of the imidazole planes of the histidine axial ligands. The solution structure was determined by using 1,303 meaningful NOEs and 241 pseudocontact shifts, the latter being derived from the native reduced protein. A family of 40 energy-minimized conformers was obtained with average RMSD of 0.56+/-0.09 A and 1.04+/-0.12 A for backbone and heavy atoms, respectively, and distance and pseudocontact shift penalty functions of 0.50+/-0.07 A2 and 0.51+/-0.02 ppm2. The structure shows some changes around the cavity and in particular a movement of the 60-70 backbone segment owing to the absence of two hydrogen bonds between the Ser64 backbone NH and side-chain OH and the carboxylate oxygen of propionate-7, present in the native protein. The analysis of the NMR spectra in the presence of unfolding agents indicates that this protein is less stable than the native form. The decrease in stability may be the result of the loss of the two hydrogen bonds connecting propionate-7 to Ser64 in the native protein. The available data on the reduction potential and the electron transfer rates are discussed on the basis of the present structural data.

Structure, interaction and electron transfer between cytochrome b5, its E44A and/or E56A mutants and cytochrome c

Site-directed mutagenesis has been used to produce variants of a tryptic fragment of bovine liver cytochrome b5 in which Glu44 and Glu56 are mutated to alanine. The reduction potentials measured by spectroelectrochemical titration (in the presence of 1 mM (Ru(NH3)6)3+, pH 7.0 and I=0.1 M) are 4.5, 6.0, 6.0 and 7.5 mV versus the standard hydrogen electrode (SHE) for the wild-type and E44A, E56A and E44/56A mutants of cytochrome b5, respectively. A comparative two-dimensional NMR study of cytochrome b5 and its E44/56A mutant in water solution has been achieved. Resonance assignments of side-chains have been completed successfully. The NMR results suggest that the secondary structures and global folding of the E44/56A mutant remain unchanged, but the mutation of both Glu44 and Glu56 to hydrophobic alanine may lead to the two helices containing mutated residues contracting towards the heme center. The inner mobility of the Gly42 approximately Glu44 segment in cytochrome b5 may be responsible for the difference of the binding mode between Glu44 and Glu56 with cytochrome c. The binding between cytochrome c and cytochrome b5 was studied by optical difference spectra of cytochrome c and variants of cytochrome b5. The association constants (KA) for the wild-type, E44A, E56A, and E44/56A mutants of cytochrome b5 with cytochrome c, are 4.70(+/-0. 10)x10(6) M-1, 1.88(+/-0.03)x10(6) M-1, 2.70(+/-0.13)x10(6) M-1, and 1.14(+/-0.05)x10(6) M-1, respectively. This is indicative that both Glu44 and Glu56 are involved in the complex formation between cytochrome b5 and cytochrome c. The reduction of horse heart ferricytochrome c by recombinant ferrocytochrome b5 and its mutants has been studied. The rate constant of the electron transfer reaction between ferricytochrome c and wild-type ferrocytochrome b5 (1.074(+/-0.49)x10(7) M-1 s-1) is higher than those of the mutant protein E44A (8.98(+/-0.20)x10(6) M-1 s-1), E56A (8.76(+/-0. 39)x10(6) M-1 s-1), and E44/56A (8.02(+/-0.38)x10(6) M-1 s-1) at 15 degreesC, pH 7.0, I=0.35 M. The rate constants are strongly dependent on ionic strength and temperature. These studies, by means of a series of techniques, provide conclusive results that the interaction between cytochrome b5 and cytochrome c is electrostatically guided, and, more importantly, that both Glu44 and Glu56 participate in the electron transfer reaction.

N epsilon,N epsilon-dimethyl-lysine cytochrome c as an NMR probe for lysine involvement in protein-protein complex formation

The reductively dimethylated derivatives of horse and yeast iso-1-ferricytochromes c have been prepared and characterized for use as NMR probes of the complexes formed by cytochrome c with bovine liver cytochrome b5 and yeast cytochrome c peroxidase. The electrostatic properties and structures of the derivatized cytochromes are not significantly perturbed by the modifications; neither are the electrostatics of protein-protein complex formation or rates of interprotein electron transfer. Two-dimensional 1H-13C NMR spectroscopy of the complexes formed by the derivatized cytochromes with cytochrome b5 and cytochrome c peroxidase has been used to investigate the number and identity of lysine residues of cytochrome c that are involved in interprotein interactions of cytochrome c. The NMR data are incompatible with simple static models proposed previously for the complexes formed by these proteins, but are consistent with the presence of multiple, interconverting complexes of comparable stability, consistent with studies employing Brownian dynamics to model the complexes. The NMR characteristics of the Nepsilon,Nepsilon-dimethyl-lysine groups, their chemical shift dispersion, oxidation state and temperature dependences and the possibility of chemical exchange phenomena are discussed with relevance to the utility of Nepsilon, Nepsilon-dimethyl-lysine's being a generally useful derivative for characterizing protein-protein complexes.

Electrostatic interaction between NADH-cytochrome b5 reductase and cytochrome b5 studied by site-directed mutagenesis

Electrostatic interaction between NADH-cytochrome b5 reductase and cytochrome b5 was studied by site-directed mutagenesis. The target residues for mutagenesis were selected on the basis of the previously reported chemical cross-linking study of these two proteins, which implicated possible charge-pair interactions between Lys-41, Lys-125, Lys-162, and Lys-163 of the enzyme, and Glu-47, Glu-48, Glu-52, Glu-60, Asp-64 (group A), and heme propionate of cytochrome b5. Mutant reductases that lost one of the above-listed Lys residues showed higher K(m) values for cytochrome b5 and lower kcat values than those of the wild type, suggesting that all of the examined Lys residues participate in binding with cytochrome b5 as reported previously. In contrast, a removal of one of (or even all of) the group A residues from cytochrome b5 by mutagenesis caused no significant effect on the catalytic properties of cytochrome b5. Additional elimination of another set of negative residues (Glu-41, Glu-42, Asp-57, and Glu-63 (Group B)), which are also located close to heme, elevated the K(m) value by more than five folds. These results suggest that there should be other acidic residue(s) than group A in cytochrome b5 which participate in binding with NADH-cytochrome b5 reductase.

Pseudocontact shifts used in the restraint of the solution structures of electron transfer complexes

The geometry of the ferricytochrome b5-ferricytochrome c complex has been analysed using long-range interprotein paramagnetic dipolar shifts. Heteronuclear filtered NMR spectra of samples containing 15N-labelled cytochrome b5 in complex with unlabelled cytochrome c allowed unambiguous assessment of pseudocontact shifts relative to diamagnetic reference states. Because pseudocontact shifts can be observed for protons as much as 20 A from the paramagnetic centre, this approach allows study of electron transfer proteins in fast exchange. Our findings provide the first physical evidence confirming hypotheses presented in previous theoretical studies. This absence of certain predicted shifts that are expected based on the best fit to a static model of the complex suggests that cytochrome b5 is more dynamic in solution than in the crystal, in agreement with molecular dynamics simulations.

Comparison of low oxidoreduction potential cytochrome c553 from Desulfovibrio vulgaris with the class I cytochrome c family

The cytochrome c553 from Desulfovibrio vulgaris (DvH c553) is of importance in the understanding of the relationship of structure and function of cytochrome c due to its lack of sequence homology with other cytochromes, and its abnormally low oxido-reduction potential. In evolutionary terms, this protein also represents an important reference point for the understanding of both bacterial and mitochondrial cytochromes c. Using the recently determined nuclear magnetic resonance (NMR) structure of the reduced protein we compare the structural, dynamic, and functional characteristics of DvH c553 with members of both the mitochondrial and bacterial cytochromes c to characterize the protein in the context of the cytochrome c family, and to understand better the control of oxide-reduction potential in electron transfer proteins. Despite the low sequence homology, striking structural similarities between this protein and representatives of both eukaryotic [cytochrome c from tuna (tuna c)] and prokaryotic [Pseudomonas aeruginosa c551 (Psa c551)] cytochromes c have been recognized. The previously observed helical core is also found in the DvH c553. The structural framework and hydrogen bonding network of the DvH c553 is most similar to that of the tuna c, with the exception of an insertion loop of 24 residues closing the heme pocket and protecting the propionates, which is absent in the DvH c553. In contrast, the Psa c551 protects the propionates from the solvent principally by extending the methionine ligand arm. The electrostatic distribution at the recognized encounter surface around the heme in the mitochondrial cytochrome is reproduced in the DvH c553, and corresponding hydrogen bonding networks, particularly in the vicinity of the heme cleft, exist in both molecules. Thus, although the cytochrome DvH c553 exhibits higher primary sequence homology to other bacterial cytochromes c, the structural and physical homology is significantly greater with respect to the mitochondrial cytochrome c. The major structural and functional difference is the absence of solvent protection for the heme, differentiating this cytochrome from both reference cytochromes, which have evolved different mechanisms to cover the propionates. This suggests that the abnormal redox potential of the DvH c553 is linked to the raised accessibility of the heme and supports the theory that redox potential in cytochromes is controlled by heme propionate solvent accessibility.

Beta-thiopropionyl cytochromes c modified at lysyl residues: preparation and characterization of the monosubstituted horse cytochromes c

beta-Thiopropionyl derivatives of horse cytochrome c singly modified at each of 18 different lysine epsilon-amino groups have been prepared using sulfosuccinimidyl-2-(biotinamido)ethyl-1,3-dithiopropionate and purified to homogeneity by high-pressure liquid chromatography. These derivatives were characterized by determination of: (i) the location of the modification; (ii) reduction potentials; (iii) visible and NMR spectra: and by (iv) measurement of electron transfer activity with cytochrome-c oxidase. No significant changes in structure were indicated, except for the ferric forms of the derivatives modified at lysines 72, 73, and 79 which are discussed separately. The electron transfer activity of the beta-thiopropionyl cytochromes c with bovine heart cytochrome-c oxidase was decreased to extents dependent on the position of the modification. Aminoethylation, a secondary modification which reverses the charge change, restored the electron transfer rate to that observed with the unmodified cytochrome c, irrespective of the location of the primary modification. These results afford a direct experimental demonstration that alterations in kinetics with physiological electron transfer partners resulting from modifications which cause a change of the charge of surface side chains are solely due to the electrostatic effects. Of the many chemically modified cytochromes c prepared to date, the singly substituted beta-thiopropionyl cytochromes c are likely to be particularly useful as the thiol allows covalent linkage of any sulfhydryl-reactive reagent to a well-defined location on the protein surface by a simple procedure, even when the secondary modifier is relatively unstable, a crucial advantage not otherwise readily achieved.

Molecular dynamics simulation of cytochrome b5: implications for protein-protein recognition

Cytochrome b5 participates in electron-transfer reactions with a variety of different proteins. To explore how this protein might discern between structurally varied proteins, we have performed a molecular dynamics simulation focusing on its structural stability and dynamic behavior in solution. The protein was simulated in water at 298 K and pH 6.9 for 2.5 ns. The protein deviated significantly from the crystal structure midway through the simulation, but ultimately the crystalline conformation was regained. The simulation was at all times well behaved as judged by comparison to structural NMR data obtained in solution. One region of the protein backbone that deviated from the crystal conformation contains acidic residues implicated in electrostatic-based protein-protein recognition. The mobility in this region caused the protein to display different patterns of residues at the surface with time, as well as the formation of a large cleft partially exposing the hydrophobic core lining the heme pocket. Furthermore, the position and cyclical formation of this cleft suggest that hydrophobic interactions may be important in protein-protein recognition events and possibly even electron transfer, as the cleft allows for easy access to the heme group. These results indicate that thermal motion could provide a low-energy mechanism for controlling recognition events. Thus, the dynamical behavior observed through the varying solution conformations sampled may be important in influencing the diverse range of protein-protein interactions in which cytochrome b5 participates.

Experimental and theoretical analysis of the interaction between cytochrome c and cytochrome b5

Experimental and theoretical investigation of the interaction of cytochrome c and cytochrome b5 performed over nearly twenty years has produced considerable insight into the manner in which these proteins recognize and bind to each other. The results of these studies and the experimental and theoretical strategies that have been developed to achieve these results have significant implications for understanding the behavior of similar complexes formed by more complex and less-well characterized electron transfer proteins. The current review provides a comprehensive summary and critical evaluation of the literature on which the current status of our understanding of the interaction of cytochrome c and cytochrome b5 is based. The general issues related to the study of electron transfer complexes of this type are discussed and some new directions for future investigation of such systems are considered.

Electron transfer from cytochrome b5 to cytochrome c

The reaction of cytochrome b5 with cytochrome c has become a very prominent system for investigating fundamental questions regarding interprotein electron transfer. One of the first computer modeling studies of electron transfer and protein/protein interaction was reported using this system. Subsequently, numerous studies focused on the experimental determination of the features which control protein/protein interactions. Kinetic measurements of the intracomplex electron transfer reaction have only appeared in the last 10 years. The current review will provide a summary of the kinetic measurements and a critical assessment of the interpretation of these experiments.

Cytochrome b5, its functions, structure and membrane topology

The first part of the present communication reviews recent advances in our understanding of the known physiological functions of cytochrome b5. In addition, one section is devoted to a description of a recently discovered function of cytochrome b5, namely its involvement in the synthesis of the oncofetal antigen N-glycolylneuraminic acid. The second part of the article summarizes site-directed mutagenesis studies, primarily conducted in the author's laboratory, in both the catalytic heme-binding and membrane-binding domain of cytochrome b5. These studies have shown that: 1) the membrane binding domain of cytochrome b5 spans the bilayer; 2) cytochrome b5 lacking 19 COOH-terminal amino acids does not bind to membrane bilayers; and 3) specific amino acids in the membrane binding domain have been mutated and shown not to be essential for the function of cytochrome b5 with its redox partners.

Thermodynamics of ferredoxin binding to ferredoxin:NADP+ reductase and the role of water at the complex interface

The association of ferredoxin with ferredoxin:NADP+ reductase (both proteins from spinach chloroplasts) was characterized by isothermal titration calorimetry and fluorescence quenching titration. The formation of the complex is mainly driven by a positive entropy change (delta S = 125 +/- 8 J mol-1 K-1). The calorimetric enthalpy of binding is small between 10 and 37 degrees C and either negative or positive, with an inversion temperature near 25 degrees C. The pH dependence of the association constant [Batie, C. J., & Kamin, H. (1981) J. Biol. Chem. 256, 7756-7763] was shown to correlate with the uptake of a single proton by a group exhibiting a heat of protonation of -26 kJ mol-1. This value agrees with the protonation of an imidazole group. Possible residues to become protonated in the complex are His-19 or His-90 of ferredoxin:NADP+ reductase. The temperature dependence of the free energy of binding, delta G, is weak because of the enthalpy-entropy compensation caused by a heat capacity change, delta Cp, of -680 +/- 44 J mol-1 K-1. The favorable binding entropy and the negative delta Cp indicate a large contribution to binding from hydrophobic effects, which seem to originate from dehydration of the protein-protein interface. Dehydration was demonstrated by osmotic stress experiments in which the association constant was found to increase by 2-4-fold in the presence of 52% (w/w) glycerol. The increase in the association constant with osmotic pressure points to the release of several water molecules from the complex interface.

A "parallel plate" electrostatic model for bimolecular rate constants applied to electron transfer proteins

A "parallel plate" model describing the electrostatic potential energy of protein-protein interactions is presented that provides an analytical representation of the effect of ionic strength on a biomolecular rate constant. The model takes into account the asymmetric distribution of charge on the surface of the protein and localized charges at the site of electron transfer that are modeled as elements of a parallel plate condenser. Both monopolar and dipolar interactions are included. Examples of simple (monophasic) and complex (biphasic) ionic strength dependencies obtained from experiments with several electron transfer protein systems are presented, all of which can be accommodated by the model. The simple cases do not require the use of both monopolar and dipolar terms (i.e., they can be fit well by either alone). The biphasic dependencies can be fit only by using dipolar and monopolar terms of opposite sign, which is physically unreasonable for the molecules considered. Alternatively, the high ionic strength portion of the complex dependencies can be fit using either the monopolar term alone or the complete equation; this assumes a model in which such behavior is a consequence of electron transfer mechanisms involving changes in orientation or site of reaction as the ionic strength is varied. Based on these analyses, we conclude that the principal applications of the model presented here are to provide information about the structural properties of intermediate electron transfer complexes and to quantify comparisons between related proteins or site-specific mutants. We also conclude that the relative contributions of monopolar and dipolar effects to protein electron transfer kinetics cannot be evaluated from experimental data by present approximations.

Structural water in oxidized and reduced horse heart cytochrome c

The existence of structural water in the interior of both oxidized and reduced horse-heart cytochrome c in solution is demonstrated using nuclear magnetic resonance spectroscopy. Six water molecules have been located in ferrocytochrome c and five in ferricytochrome c, with residence times greater than a few hundred picoseconds. Two water molecules are located in the haem crevice, one of which is found to undergo a large change in position with a change of oxidation state. Both of these observations indicate that buried structural waters in the haem crevice have, by microscopic dielectric effects, significant roles in the setting of the solvent reorganization energy associated with electron transfer.

Tertiary structure of the heme-binding domain of rat cytochrome b5 based on homology modeling

The in vitro complexes formed between cytochrome b5 and other proteins (e.g. cytochrome c) have served as a useful means to probe electrostatic contributions to macromolecular recognition. Extensive experimentation has been carried out to test the specificity and stability of these complexes, including site-directed mutagenesis based on the heterologous expression of rat cytochrome b5 in E. coli. Despite this interest, there has not been a determination of the complete structure of cytochrome b5. Here we report coordinates for the complete tertiary structure of the heme-binding domain of rat cytochrome b5 based on homology modeling. Protein Data Bank (PDB) coordinates derived from the crystal structure of the highly homologous bovine cytochrome b5 were used for main chain scaffolding. Secondary structures for the termini missing in the bovine structure were generated using homologous sequences derived from an exhaustive search of the PDB database. The model structure was solvated and further refined using energy minimization techniques. The N-terminal residues of the model appear to be in a beta sheet conformation while the carboxy terminus is in a helical conformation. The rest of the rat model is folded virtually identically to the bovine x-ray crystal structure (r.m.s. deviation 1.28 A), despite six sequence differences between the two cores. This homology-based structure should be useful for structure-function analyses of molecular recognition involving cytochrome b5.

Analysis of the bimolecular reduction of ferricytochrome c by ferrocytochrome b5 through mutagenesis and molecular modelling

Site-directed mutagenesis has been used to produce variants of cytochrome c in which selected structural or functional properties of this protein are altered that have been implicated previously in contributing to the rate at which ferricytochrome c is reduced by ferrocytochrome b5. In total, 18 variants have been studied by kinetics and electrochemical methods to assess the contributions of thermodynamic driving force, surface charge and hydrophobic interactions, and redox-linked structural reorganization of the protein to the rate of electron transfer between these two proteins under conditions where the reaction is bimolecular. While some variants (those at position-38) appear to affect primarily the driving force of the reaction, others appear to influence the rearrangement barrier to electron transfer (those at positions-67 and -52) while the interface between electron donor and acceptor centers is the principal effect of substitutions for a conserved aromatic heme contact residue at the surface of the protein (position-82). Interpretation of these results has been facilitated through the use of energy minimization calculations to refine the hypothetical models previously suggested for the cytochrome c- cytochrome b5 precursor complex on the basis of Brownian dynamics simulations of the bimolecular encounter event.

The cytochrome b5-fold: an adaptable module

The family of b5-like cytochromes encompasses, besides cytochrome b5 itself, hemoprotein domains covalently associated with other redox proteins, in flavocytochrome b2 (L-lactate dehydrogenase), sulfite oxidase and assimilatory nitrate reductase. A comparison of about 40 amino acid sequences deposited in data banks shows that eight residues are invariant and about 15 positions carry strongly conservative substitutions. Examination of the location of these invariant and conserved positions in the light of the three-dimensional structures of beef cytochrome b5 and S cerevisiae flavocytochrome b2 suggests a strongly conserved protein structure for the b5-like heme-binding domain throughout evolution. Numerous NMR studies have demonstrated the existence of a positional isomerism for the heme, which involves both a 180 degree-rotation around the heme alpha,gamma-meso carbon atoms and a rotation through an axis normal to the heme plane at the iron. NMR studies did not detect significant differences in protein structure between reduced and oxidized states, or between species. The role of a number of side chains was probed by site-directed mutagenesis. Studies of complex formation and of electron transfer rates between cytochrome b5 and redox partners have led to the idea that complexation is driven by electrostatic forces, that it is generally the exposed heme edge which makes contact with electron donors and acceptors, but that there are multiple overlapping sites within this general area. For the bi- and trifunctional members of the family, extrapolation of available data would suggest a mobile heme-binding domain within a complex structure. In these cases the existence of a single interaction area for both electron donor and acceptor, or of two different ones, remains open to discussion.

Responses of two protein-protein complexes to solvent stress: does water play a role at the interface?

We have analyzed the stability of the cytochrome c-cytochrome b5 and cytochrome c-cytochrome c oxidase complexes as a function of solvent stress. High concentrations of glycerol were used to displace the two equilibria. Glycerol promotes complex formation between cytochrome c and cytochrome b5 but inhibits that between cytochrome c and cytochrome c oxidase. The results with cytochrome b5 and cytochrome c were expected; the association of this complex is largely entropy driven. Our interpretation is that the cytochrome c-cytochrome b5 complex excludes water. The results with the cytochrome c oxidase and cytochrome c couple were not expected. We interpret them to mean that either glycerol is binding to the oxidase, thereby displacing the cytochrome c, or that water is required at this protein-protein interface. A requirement for substantial quantities of water at the interface of some protein complexes is logical but has been reported only once.

Species-specific differences in covalently crosslinked complexes of yeast cytochrome c peroxidase with horse and yeast iso-1 ferricytochromes c

1. The results of chemically crosslinking yeast cytochrome c peroxidase with both horse and yeast iso-1 ferricytochromes c have been studied by a combination of gel electrophoresis and proton NMR spectroscopy. 2. The complexes were formed at a variety of potassium phosphate concentrations ranging from 10 to 300 mM using the water soluble crosslinking agent, EDC (1-ethyl-3-[3-(dimethylamino)propyl]-carbodiimide). 3. The primary crosslinking product in both cases is the 1:1 covalent complex, but, for each pair of partner proteins the yield of the 1:1 crosslinked complex varies with the salt concentration. 4. Furthermore, at low salt concentrations the yield of the 1:1 covalent complex involving horse cytochrome c is much larger than the yield of the 1:1 covalent complex formed with yeast iso-1 cytochrome c, whereas at high salt concentrations the situation is reversed. 5. Proton NMR spectroscopy, in combination with gel electrophoresis, provides evidence for the formation of different types of 1:1 complexes for the peroxidase/yeast cytochrome c pair and has been used to study the effect of changes in the solution ionic strength upon both the peroxidases/horse cytochrome c and the peroxidase/yeast cytochrome c complexes. 6. This work indicates that electrostatic interactions between proteins play a dominant role in formation of complexes between cytochrome c peroxidase and horse ferricytochrome c, whereas the hydrophobic effect plays a comparatively larger role in stabilizing complexes between cytochrome c peroxidase and yeast iso-1 ferricytochrome c.

Intracomplex electron transfer between ruthenium-65-cytochrome b5 and position-82 variants of yeast iso-1-cytochrome c

We tested the idea that the aromatic ring on the invariant residue Phe-82 in cytochrome c acts as an electron-transfer bridge between cytochrome c and cytochrome b5. Ru-65-cyt b5 was prepared by labeling the single sulfhydryl group on T65C cytochrome b5 with [4-(bromomethyl)-4'-methylbipyridine][bis(bipyridine)]ruthenium 2+ as previously described [Willie, A., Stayton, P.S., Sligar, S.G., Durham, B., & Millett, F. (1992) Biochemistry 31, 7237-7242]. Laser excitation of the complex formed between Ru-65-cyt b5 and Saccharomyces cerevisiae iso-1-cytochrome c at low ionic strength results in rapid electron transfer from the excited-state Ru(II*) to the heme group of Ru-65-cyt b5 followed by biphasic electron transfer to the heme group of cytochrome c with rate constants of (1.0 +/- 0.2) x 10(5) s-1 and (2.0 +/- 0.04) x 10(4) s-1. Variants of iso-1-cytochrome c substituted at Phe-82 with Tyr, Gly, Leu, and Ile have fast-phase rate constants of 0.4, 1.9, 2.1, and 2.0 x 10(5) s-1 and slow-phase rate constants of 5.3, 3.5, 2.4, and 2.0 x 10(3) s-1, respectively. Increasing the ionic strength to 50 mM results in single-phase intracomplex electron transfer with rate constants of 3.8, 3.1, 3.0, 5.0, and 4.5 x 10(4) s-1 for the wild-type, Tyr, Gly, Leu, and Ile variants, respectively. These results demonstrate that an aromatic side chain at residue 82 is not needed for rapid electron transfer with cytochrome b5. Furthermore, two conformational forms of the complex are present at low ionic strength with fast and slow electron-transfer rates.(ABSTRACT TRUNCATED AT 250 WORDS)

Effects of charged amino acid mutations on the bimolecular kinetics of reduction of yeast iso-1-ferricytochrome c by bovine ferrocytochrome b5

The reduction of wild-type yeast iso-1-ferricytochrome c (ycytc) and several mutants by trypsin-solubilized bovine liver ferrocytochrome b5 (cytb5) has been studied under conditions in which the electron-transfer reaction is bimolecular. The effect of electrostatic charge modifications and steric changes on the kinetics has been determined by experimental and theoretical observations of the electron-transfer rates of ycytc mutants K79A, K'72A, K79A/K'72A, and R38A (K' is used to signify trimethyllysine (Tml)). A structurally robust Brownian dynamics (BD) method simulating diffusional docking and electron transfer was employed to predict the mutation effect on the rate constants. A realistic model of the electron-transfer event embodied in an intrinsic unimolecular rate constant is used which varies exponentially with donor-acceptor distance. The BD method quantitatively predicts rate constants over a considerable range of ionic strengths. Semiquantitative agreement is obtained in predicting the perturbing influence of the mutations on the rate constants. Both the experimentally observed rate constants and those predicted by BD descend in the following order: native ycytc > K79A > K'72A > K79A/K'72A. Variant R38A was studied at a different ionic strength than this series of mutations, and the theory agreed with experiment in predicting a smaller rate constant for the mutant. In all cases the predicted effect of mutation was in the correct direction, but not as large as that observed. The BD simulations predict that the two proteins dock through essentially a single domain, with a distance of closest approach of the two heme groups in rigid body docking typically around 12 A. Two predominant classes of complexes were calculated, the most frequent involving the quartet of cytb5/ycytc interactions, Glu48-Arg13, Glu56-Lys87, Asp60-Lys86, and heme-Tml72, having an average electrostatic energy of -13.0 kcal/mol. The second most important complexes were of the type previously postulated (Salemme, 1976; Mauk et al., 1986; Rodgers et al., 1988) with interactions Glu44-Lys27, Glu48-Arg13, Asp60-Tml72, and heme-Lys79 and having an energy of -6.4 kcal/mol. The ionic strength dependence of the bimolecular reaction rate was well reproduced using a discontinuous dielectric model, but poorly so for a uniform dielectric model.

Electron-transfer kinetics and electrostatic properties of the Rhodobacter sphaeroides reaction center and soluble c-cytochromes

The kinetics of electron transfer between the Rhodobacter sphaeroides R-26 reaction center and nine soluble c-cytochromes have been analyzed and compared to the patterns of the surface electrostatic potentials for each of the proteins. Characteristic first-order electron-transfer rates for 1:1 complexes formed at low ionic strength between the reaction center and the different c-cytochromes were identified and found to vary by a factor of almost 100, while second-order rates were found to differ by greater than 10(6). A correlation was found between the location of likely electrostatic interaction domains on each cytochrome and its characteristic rate of electron transfer. The interaction domains were identified by mapping electrostatic potentials, calculated from the Poisson-Boltzmann equation, onto simulated "encounter surfaces" for each of the cytochromes and the reaction center. For the reaction center, the c-cytochrome binding domain was found to have almost exclusively net negative potential (< -3 kT) and to be shifted slightly toward the M-subunit side of the reaction center. The location of interaction domains of complementary, positive potential (> 3 kT) differed for each cytochrome. The correspondence between electrostatic, structural, and kinetic properties of 1:1 reaction center-cytochrome complexes leads to a proposed mechanism for formation of reaction center-cytochrome electron-transfer complexes that is primarily driven by the juxtaposition of regions of delocalized complementary potential. In this mechanism the clustering of charged residues is of primary importance and not the location of specific residues. A consequence of this mechanism is that many different sets of charge distributions are predicted to be capable of stabilizing a specific configuration for a reaction center-cytochrome complex. This mechanism for reaction center association with water-soluble c-cytochromes fits molecular recognition mechanisms proposed for c-cytochromes in nonphotosynthetic systems. In general, the kinetic scheme for reaction center driven cytochrome oxidation was found to vary between a simple two-state model, involving cytochrome in free and reaction center bound states, and a three-state model, that includes cytochrome binding in kinetically competent ("proximal") and incompetent ("distal") modes. The kinetically incompetent mode of cytochrome binding is suggested not to be an intrinsic feature of the reaction center-cytochrome association but is likely to be due to variation in the physical state of the reaction center.

Motional dynamics of functional cytochrome c delivered by low pH fusion into the intermembrane space of intact mitochondria

We have investigated the motional dynamics of cytochrome c in the intact, functional rat liver mitochondrion. To do this, functional, FITC-cytochrome c (fluorescein isothiocyanate monoderivatized cytochrome c) was incorporated into the intermembrane space (IMS) of intact mitochondria through encapsulation of cytochrome c into asolectin liposomes followed by low pH-induced fusion of the liposomes with the outer membranes of the mitochondria. A cytochrome c controlled enrichment of between 15%-50% (1800-7200 molecules incorporated per mitochondrion) was obtained. All cytochrome c incorporated, regardless of the quantity, participated in the function of electron transport, indicative of a functional, independent random diffusant. Resonance energy transfer was determined from the IMS-entrapped functional FITC-cytochrome c to octadecylrhodamine B incorporated into the mitochondrial membranes. Resonance energy transfer from FITC-cytochrome c to octadecylrhodamine B in isolated inner or outer mitochondrial membranes (IMM and OMM, respectively) was also measured. We found substantial differences in the effects of ionic strength (I) on the proximity of cytochrome c to isolated IMM and OMM. Interactions with isolated IMM were very dynamic, i.e., very I-dependent, and cytochrome c binding to IMM was significant only at very low I. I-dependent interactions of cytochrome c with isolated OMM were less I-dependent than those for the IMM. However, FITC-cytochrome c was essentially released from IMM and OMM at physiological I. The proximity of FITC-cytochrome c to each mitochondrial membrane after its incorporation into the IMS of intact mitochondria in the condensed configuration was estimated at different external, bulk I using: (a) resonance energy transfer from IMS-entrapped FITC-cytochrome c to octadecylrhodamine B-label evenly distributed in both mitochondrial membranes; and (b) resonance energy transfer from IMS-entrapped FITC-cytochrome c to octadecylrhodamine B-label concentrated in the OMM. Resonance energy transfer showed that the average distance between cytochrome c and the two IMS-membrane surfaces increased with increasing IMS-I, approaching a maximal measurable distance of 85 A at 150 mM I. This result is consistent with a dissociation of FITC-cytochrome c and both membranes of intact mitochondria at physiological I, i.e., when the activity of cytochrome c in electron transport is highest. Our findings reveal a primarily three-dimensional diffusion mode for IMS-cytochrome c during its function in electron transport in intact mitochondria at physiological I, and offer further evidence that mitochondrial electron transport is a process driven by random collisions between its independently diffusing electron transferring, redox components.

Laser flash photolysis studies of electron transfer to the cytochrome b5-cytochrome c complex

Rate constants for electron transfer in the complex between recombinant rat mitochondrial outer membrane cytochrome b5 or the tryptic fragment of bovine liver cytochrome b5 and horse mitochondrial cytochrome c were measured by laser flash photolysis of 5-deazariboflavin-EDTA solutions. When an excess of cytochrome b5 was titrated with increasing amounts of cytochrome c at low ionic strength and electron transfer was initiated by a laser flash, both proteins were rapidly reduced by deazariboflavin semiquinone. The initial photoreduction was followed by a slower second-order reduction of b5 complexed oxidized cytochrome c by free reduced cytochrome b5. At an 8:1 ratio of cytochromes b5 to c, the pseudo-first-order rate constant for reduction of complexed cytochrome c increased 3-5-fold between ionic strengths of 5 and 40 mM, and then dropped precipitously at higher ionic strengths. The ionic strength dependent increase in rate constant is likely to be due to relief of steric hindrance via rearrangement of cytochrome c in the complex. The reaction rate showed no sign of saturation at any ionic strength, indicating a first-order rate constant greater than 10(4) s-1 within a transient ternary protein complex; i.e., interprotein electron transfer approaches the largest values previously reported for the stable binary protein complex (approximately 4 x 10(5) s-1). Our results emphasize the flexibility of electron-transfer protein complexes, which had previously been modeled in a single conformation with specific salt bridges. It appears that a variety of orientations can exist within such protein-protein complexes and that the population of conformations changes with ionic strength.(ABSTRACT TRUNCATED AT 250 WORDS)

Cytochrome C (Fe2+) as a competitive inhibitor of NADPH-dependent reduction of cytochrome P450 LM2: locating protein-protein interaction sites in microsomal electron carriers

The kinetics of NADPH-dependent reduction of cytochrome P450 LM2 in the soluble monomeric reconstituted system in the absence of any substrate is shown to be monophasic. We show that ferrous cytochrome c acts as a competitive inhibitor of the reduction. In the presence of 1 mM benzphetamine an additional extremely fast phase was observed. Under these conditions ferrous cytochrome c was found to be a competitive inhibitor of the slow phase of the reduction process, which accounted for 80% of the total reduction amplitude. Inhibition experiments yield a dissociation constant for the LM2-reductase complex of 3.0 +/- 1.5 microM. This constant was the same both in the presence and in the absence of benzphetamine. Based on these data we conclude that cytochromes P450 and c bind to the same center on the NADPH-cytochrome P450 reductase molecule. Comparative analysis of the amino acid sequences reveals a detectable similarity between cytochrome c and cytochrome P450 LM2 at positions 68-87 and 121-145, respectively. In addition, a substantial similarity was shown for sequence fragments 204-224 of NADPH-cytochrome P450 reductase and 40-60 of cytochrome b5. Based on these findings a hypothesis for the location of the centers of intermolecular interactions on the molecules of cytochrome P450 LM2 and NADPH-cytochrome P450 reductase is proposed.

Spectrophotometric detection of the interaction between cytochrome c and heparin

Heparin inhibits transport of electrons from reduced cytochrome c to cytochrome c oxidase. The effect is due to the interaction of heparin with cytochrome c. It has been observed that binding of heparin to the reduced or oxidized cytochrome c changes the spectrum of cytochrome c at the Soret region. Affinity chromatography of heparin in cytochrome c immobilized to thiol-Sepharose shows that commercial heparin is eluted in the low-affinity and high-affinity fractions. Both participate in the interaction with cytochrome c. Polylysine induces decay of the cytochrome c-heparin complex.

The identification of cation-binding domains on the surface of microsomal cytochrome b5 using 1H-NMR paramagnetic difference spectroscopy

One-dimensional and two-dimensional 1H-NMR methods and paramagnetic difference spectroscopy have defined cation binding domains on the surface of the tryptic fragment of microsomal cytochrome b5. The addition of tris(ethylenediamine) chromium(III) [Cr(en)3(3+)] to solutions of ferricytochrome b5 reveals at least three distinct sites on the surface of the protein to which highly charged cations may bind (20 mM phosphate pH 7.0, T = 300 K). Surprisingly only one of these sites is located close to the haem edge region of the protein, whilst the remaining two sites are more remote. Site I contains the exposed haem C13 propionate and a series of carboxylate residues that includes glutamates 37, 38, 43, 44, and 48. Sites II and III are located away from the haem edge region and are delineated by the broadening of aromatic resonances of histidines 26 and 80. Further investigation of the interaction between Cr(en)3(3+) and cytochrome b5 using two-dimensional double-quantum-filtered correlated spectroscopy shows that resonances assigned to Glu59, Asp60, Glu79, Asp82 and Asp83 are broadened with the distribution of these charged side chains correlating with the relaxation broadening observed from one-dimensional experiments. In a binary complex with ferricytochrome c, Cr(en3(3+) broadens many cytochrome b45 resonances including the haem propionates, His26, Ala54, Thr55 and His80. Although the pattern of line-broadening of resonances at sites II and III is unaltered by complex formation, cytochrome c shields residues at site I, the haem edge site. The results indicate that the interaction between cytochrome b5 and c in a binary complex involves multiple protein configurations.

Mapping electrostatic interactions in macromolecular associations

In the association of electron transfer proteins, electrostatics has been proposed to play a role in maintaining the stability and specificity of the biomolecular complexes formed. An excellent model system is the interaction between mammalian cytochrome b5 and cytochrome c, in which the X-ray structures of the individual components reveal a complementary asymmetry of charges surrounding their respective redox centers. Determining the exact extent of the electrostatic interactions and identifying the specific residues involved in the formation of the electron transfer complex has proved more elusive. We report herein the utilization of high-pressure techniques, together with site-directed mutagenesis, to provide a map of the interaction domains in biomolecular complex formation. The application of high pressure disrupts macromolecular associations since dissociation of the complex results in a decreased volume of the system due to the solvation of charges that had been previously sequestered in the interface region and force solvation of hydrophobic surfaces. Site-directed mutagenesis of a totally synthetic gene for rat liver cytochrome b5, which expresses this mammalian protein in Escherichia coli as a hemecontaining soluble component, was used to selectively alter negatively charged residues of cytochrome b5 to neutral amide side-chains. We have demonstrated that the interaction domain of cytochrome b5 with cytochrome c can be mapped from a comparison of dissociation volumes of these modified cytochrome b5-cytochrome c complexes with the native complex. Using these techniques we can specifically investigate the role of particular residues in the equilibrium association of these two electron transfer proteins. Single-point mutations in the interaction domain give nearly identical effects on the measured dissociation volumes, yet removal of acidic residues outside the recognition surface yield volumes similar to wild-type protein. Multiple mutations in the proposed protein-protein interaction site are found to allow greater solvent-accessibility of the interface as reflected in a diminution in the volume changes on subsequent charge removal. This is indicative that the interprotein salt-bridges in this complex provide a mechanism for a greater exclusion of solvent from the interfacial domain of the complex, resulting in a more stable association.

Proton linkage of complex formation between cytochrome c and cytochrome b5: electrostatic consequences of protein-protein interactions

Two potentiometric methods have been used to study the pH-dependent changes in proton binding that accompany complex formation between cytochrome c and cytochrome b5. With one method, the number of protons bound or released upon addition of one cytochrome to the other has been measured as a function of pH. The results from these studies are correlated with the complexation-induced difference titration curve calculated from the titration curves of the preformed complex and of the individual proteins. Both methods demonstrate that complex formation at acid pH is accompanied by proton release, that complex formation at basic pH is accompanied by proton uptake, and that the change in proton binding at neutral pH, where stability of complex formation is maximal, is relatively small. Under all conditions studied, the stoichiometry of cytochrome c-cytochrome b5 complex formation is 1:1 with no evidence of higher order complex formation. Although the dependence of complex formation on pH for interaction between different species of cytochrome c and cytochrome b5 are qualitatively similar, they are quantitatively different. In particular, complex formation between yeast iso-1-cytochrome c and lipase-solubilized bovine cytochrome b5 occurs with a stability constant that is 10-fold greater than observed for the other two pairs of proteins under all conditions studied. Interaction between these two proteins is also significantly less dependent on ionic strength than observed for complexes formed by horse heart cytochrome c with either form of cytochrome b5.(ABSTRACT TRUNCATED AT 250 WORDS)

1H NMR study of the role of heme carboxylate side chains in modulating heme pocket structure and the mechanism of reconstitution of cytochrome b5

1H nuclear magnetic resonance spectroscopy was used to assign the hyperfine-shifted resonances and determine the position of a side chain in the heme cavity of wild-type rat apocytochrome b5 reconstituted with a series of synthetic hemins possessing systematically perturbed carboxylate side chains. The hemins included protohemin derivatives with individually removed or pairwise shortened and lengthened carboxylate side chains, as well as (propionate)n(methyl)8-nporphine-iron(III) isomers with n = 1-3 designed to force occupation of nonnative propionate sites. The resonance assignments were effected on the basis of available empirical heme contact shift correlations and steady-state nuclear Overhauser effect measurements in the low-spin oxidized proteins. The failure to detect holoproteins with certain hemins dictates that the stable holoproteins, unlike the case of myoglobin, demand the axial iron-His bonds and cannot accommodate carboxylate side chains at interior positions in the binding pocket. Hence, the heme pocket interior in cytochrome b5 is judged much less polar and less sterically accommodating than that of myoglobin. The propionate occupational preference was greatest as the native 7-propionate site, but also possible at the nonnative crystallographic 5-methyl or 8-methyl positions. Only for a propionate at the crystallographic 8-methyl position was a significant perturbation of the native molecular/electronic structure observed, and this was attributed to an alternative propionate-protein hydrogen bond at the crystallographic 8-methyl position. The structures of the transient protein complexes detected only shortly after reconstitution reveal that the initial encounter complexes during assembly of holoprotein from apoprotein and hemin involve one of the two alternate propionate-protein links at either the 7-propionate or native 8-methyl position. In a monopropionate hemin, this leads to the characterization of a new type of heme orientational disorder involving rotation about a N-Fe-N axis.

The formation of protein complexes between ferricytochrome b5 and ferricytochrome c studied using high-resolution 1H-NMR spectroscopy

The association of the tryptic fragment of bovine microsomal cytochrome b5 with cytochrome c has been studied by one- and two-dimensional 1H-NMR spectroscopy. The association of cytochromes to form protein complexes is apparent from the increase in linewidths for resonances of ferricytochrome b5 as well as small perturbations in their chemical shifts that occur upon increasing the cytochrome c/b5 molar ratio. The changes in the chemical shifts of hyperfine shifted resonances of ferricytochrome b5 with increasing ratios of ferricytochrome c indicate the formation of binary 1:1 complexes and ternary 1:2 complexes. Similarly, titrations of the linewidth of resolved resonances of ferricytochrome b5 are consistent with stoichiometries of 1:1 and 1:2 for complexes formed between cytochromes b5 and c. Surprisingly, in the 1:1 complex, mobility is shown to be a function of ionic strength. Two-dimensional correlated spectroscopy (COSY) and nuclear Overhauser enhancement spectroscopy (NOESY) of the binary complex formed between ferricytochrome b5 and c indicate that the positions of many resonances attributable to amino acids are unaltered by protein association, although distinctive chemical shift changes are detected in the alpha-CH of the haem C17 propionate. The protein complex detected by NMR is discussed with respect to the model for the binary complex proposed by Salemme and possible mechanisms of electron transfer.

Visualization of molecular flexibility and its effects on electrostatic recognition

To study the effect of protein flexibility on electrostatic recognition, we have devised two novel computer graphic representations of the changes in the electrostatic field of a protein resulting from its internal motions. The atomic structure of Cu, Zn superoxide dismutase was minimized, and the 200 lowest frequency normal modes of the enzyme were determined. Individual and combined normal-mode vibrations were visualized interactively with the program Flex. Normal-mode motions are fast enough (approximately 10(-11) s cycle-1) to evade solvent damping, thus allowing long-range electrostatic interactions to dominate. The changing electrostatic environment of the protein was examined by animating precalculated frames of electrostatic field vectors with GRAMPS. With Vu, changes in electrostatic potential were displayed as variations in the color-coding of dots lying on a consensus surface that maintains the protein's shape. The consensus surface was calculated with the program Sphinx, and was derived from spherical harmonic approximations of expanded molecular surfaces. The ability to view the effects of molecular motions interactively should be useful in understanding the relationships of protein structure to function.

Direct electrochemistry of protein-protein complexes involving cytochrome c, cytochrome b5, and plastocyanin

The direct electrochemistry of the cytochrome c/cytochrome b5 and cytochrome c/plastocyanin complexes has been investigated at edge-plane graphite and modified gold electrode surfaces, which are selective for one of the two components of the complex. Electrochemical response of one protein at an otherwise electrostatically unfavorable electrode surface was achieved in the presence of the other protein, and the calculated heterogeneous electron-transfer rate constant and diffusion coefficient were found to be in good agreement with the values determined previously from the electrochemistry of the individual proteins [Armstrong, F. A., Hill, H. A. O., & Walton, N. J. (1988) Acc. Chem. Res. 21, 407 and references therein]. A dynamic model of the protein-protein-electrode ternary complex is proposed to explain the promotion effect, and this model is supported by a study comparing the electrochemical responses of covalent and electrostatic cytochrome c/plastocyanin complexes. It is also suggested that the behavior of protein-protein complexes at electrode surfaces could be related to that of the complexes associated with biological membranes.

Structural studies with the uveopathogenic peptide M derived from retinal S-antigen

The 18-residue fragment of bovine S-antigen, corresponding to amino acid positions 303-320, is highly immunogenic and is known to induce experimental autoimmune uveitis. The solution conformation of this immunogenic peptide, known as peptide M, was studied by Fourier-transform infrared spectroscopy and by circular dichroism. In the pH range between approximately 4 and 9.5, peptide M has a strong tendency to form macromolecular assemblies in which it adopts an intermolecular beta-sheet structure. The intermolecular beta-sheets are stabilized by ionic interactions ("salt bridges") between the carboxylate groups and basic residues of the neighboring peptide molecules. These interactions can be disrupted by neutralization of either acidic (pH range below 4) or basic residues (pH range above 9.5) or by elevated hydrostatic pressure. The secondary structure of the peptide under conditions favoring the monomeric state appears to be a mixture of unordered structure and beta-sheets. The present data are consistent with a recently proposed model [Sette, A., Buns, S., Colon, S., Smith, J. A., Miles, C., & Grey, H. M. (1987) Nature 328, 395-399], which assumes that certain immunogenic peptides adopt an extended beta-type conformation in which they are "sandwiched" between the major histocompatibility complex and the T-cell receptor.