Dynamic docking and electron transfer between Zn-myoglobin and cytochrome b(5)

We present a broad study of the effect of neutralizing the two negative charges of the Mb propionates on the interaction and electron transfer (ET) between horse Mb and bovine cyt b(5), through use of Zn-substituted Mb (ZnMb, 1) to study the photoinitiated reaction, ((3)ZnP)Mb + Fe(3+)cyt b(5) --> (ZnP)(+)Mb + Fe(2+)cyt b(5). The charge neutralization has been carried out both by replacing the Mb heme with zinc-deuteroporphyrin dimethylester (ZnMb(dme), 2), which replaces the charges by small neutral hydrophobic patches, and also by replacement with the newly prepared zinc-deuteroporphyrin diamide (ZnMb(diamide), 3), which converts the charged groups to neutral, hydrophilic ones. The effect of propionate neutralization on the conformation of the zinc-porphyrin in the Mb heme pocket has been studied by multinuclear NMR with an (15)N labeled zinc porphyrin derivative (ZnMb((15)N-diamide), 4). The rates of photoinitiated ET between the Mb's (1-3) and cyt b(5) have been measured over a range of pH values and ionic strengths. Isothermal titration calorimetry (ITC) and NMR methods have been used to independently investigate the effect of charge neutralization on Mb/b(5) binding. The neutralization of the two heme propionates of ZnMb by formation of the heme diester or, for the first time, the diamide increases the second-order rate constant of the ET reaction between ZnMb and cyt b(5) by as much as several 100-fold, depending on pH and ionic strength, while causing negligible changes in binding affinity. Brownian dynamic (BD) simulations and ET pathway calculations provide insight into the protein docking and ET process. The results support a new "dynamic docking" paradigm for protein-protein reactions in which numerous weakly bound conformations of the docked complex contribute to the binding of cyt b(5) to Mb and Hb, but only a very small subset of these are ET active, and this subset does not include the conformations most favorable for binding; the Mb surface is a large "target" with a small "bullseye" for the cyt b(5) "arrow". This paradigm differs sharply from the more familiar, "simple" docking within a single, or narrow range of conformations, where binding strength and ET reactivity increase in parallel. Likewise, it is distinct from, although complementary to, the well-known picture of conformational control of ET through "gating", or a related picture of "conformational coupling". The new model describes situations in which tight binding does not correlate with efficient ET reactivity, and explains how it is possible to modulate reactivity without changing affinity. Such "decoupling" of reactivity from binding clearly is of physiological relevance for the reduction of met-Mb in muscle and of met-Hb in a red cell, where tight binding of cyt b(5) to the high concentration of ferrous-Mb/Hb would prevent the cytochrome from finding and reducing the oxidized proteins; it likely is of physiological relevance in other situations, as well.     

Dynamic docking and electron-transfer between cytochrome b5 and a suite of myoglobin surface-charge mutants. Introduction of a functional-docking algorithm for protein-protein complexes

Horse myoglobin (Mb) provides a convenient "workbench" for probing the effects of electrostatics on binding and reactivity in the dynamic [Mb, cytochrome b(5)] electron-transfer (ET) complex. We have combined mutagenesis and heme neutralization to prepare a suite of six Mb surface-charge variants: the [S92D]Mb and [V67R]Mb mutants introduce additional charges on the "front" face, and incorporation of the heme di-ester into each of these neutralizes the charge on the heme propionates which further increases the positive charge on the "front" face. For this set of mutants, the nominal charge of Mb changes by -1 to +3 units relative to that for native Mb. For each member of this set, we have measured the bimolecular quenching rate constant (k(2)) for the photoinitiated (3)ZnDMb --> Fe(3+)b(5) ET reaction as a function of ionic strength. We find: (i) a dramatic decoupling of binding and reactivity, in which k(2) varies approximately 10(3)-fold within the suite of Mbs without a significant change in binding affinity; (ii) the ET reaction occurs within the "thermodynamic" or "rapid exchange" limit of the "Dynamic Docking" model, in which a large ensemble of weakly bound protein-protein configurations contribute to binding, but only a few are reactive, as shown by the fact that the zero-ionic-strength bimolecular rate constant varies exponentially with the net charge on Mb; (iii) Brownian dynamic docking profiles allow us to visualize the microscopic basis of dynamic docking. To describe these results we present a new theoretical approach which mathematically combines PATHWAY donor/acceptor coupling calculations with Poisson-Boltzmann-based electrostatics estimates of the docking energetics in a Monte Carlo (MC) sampling framework that is thus specially tailored to the intermolecular ET problem. This procedure is extremely efficient because it targets only the functionally active complex geometries by introducing a "reactivity filter" into the computations themselves, rather than as a subsequent step. This efficiency allows us to employ more computationally expensive and accurate methods to describe the relevant intermolecular interaction energies and the protein-mediated donor/acceptor coupling interactions. It is employed here to compute the changes in the bimolecular rate constant for ET between Mb and cyt b(5) upon variations in the myoglobin surface charge, pH, and ionic strength.     

Investigation of the electron-transfer mechanism by cross-linking between Zn-substituted myoglobin and cytochrome b(5)

We have investigated the photoinduced electron transfer (ET) in the 1:1 cross-linked complex (CL-ZnMb/b(5)) formed by a cross-linking reagent, EDC, between Zn-substituted myoglobin (ZnMb) and cytochrome b(5) (Cytb(5)) to reveal the mechanism of the inter-protein ET reactions under the condition of multiple encounter complexes. A variety of the ZnMb-Cytb(5) orientations was suggested because of failure to identify the single and specific cross-linking site on Cytb(5) by the peptide-mapping analysis using mass spectrometry. In CL-ZnMb/b(5), a laser pulse generates the triplet excited state of the ZnMb domain ((3)ZnMb()), which can transfer one electron to the Cytb(5) domain. The decay kinetics of (3)ZnMb() in CL-ZnMb/b(5) consists of a facile power-law ET phase to Cytb(5) domain ( approximately 30%) and a slower single-exponential phase ( approximately 70%). The application of the Marcus equation to this power-law phase indicates that CL-ZnMb/b(5) has a variety of ZnMb-Cytb(5) orientations for the facile ET in which the distance between the redox centers (D-A distance) is distributed over 13-20 A. The single-exponential phase in the (3)ZnMb() decay kinetics of CL-ZnMb/b(5) is similar to the intrinsic decay of (3)ZnMb() in its rate constant, 65 s(-)(1). This implies that the ET is impeded in about 70% of the total ZnMb-Cytb(5) orientations due to the D-A distance larger than 20 A. Combined with the results of the Brownian dynamics simulations for the encounter complexes, the overall bimolecular ET rate, k(app), can be reproduced by the sum of the ET rates for the minor encounter complexes of which D-A distance is less than 20 A. On the other hand, the encounter complexes with longer D-A distance, which are the majority of the encounter complexes between ZnMb and Cytb(5), have little contribution to the overall bimolecular ET rate. These observations experimentally demonstrate that ZnMb forms a variety of encounter complexes with Cytb(5), among which a minor set of the complexes with the shorter D-A distance (< approximately 20 A) regulates the overall bimolecular ET between the proteins.     

Unequivocal determination of metal atom oxidation state in naked heme proteins: Fe(III)myoglobin, Fe(III)cytochrome c, Fe(III)cytochrome b5, and Fe(III)cytochrome b5 L47R

Unambiguous determination of metal atom oxidation state in an intact metalloprotein is achieved by matching experimental (electrospray ionization 9.4 tesla Fourier transform ion cyclotron resonance) and theoretical isotopic abundance mass distributions for one or more holoprotein charge states. The ion atom oxidation state is determined unequivocally as Fe(III) for each of four gas-phase unhydrated heme proteins electrosprayed from H2O: myoglobin, cytochrome c, cytochrome b5, and cytochrome b5 L47R (i.e., the solution-phase oxidation state is conserved following electrospray to produce gas-phase ions). However, the same Fe(III) oxidation state in all four heme proteins is observed after prior reduction by sodium dithionite to produce Fe(II) heme proteins in solution: thus proving that oxygen was present during the electrospray process. Those results bear directly on the issue of similarity (or lack thereof) of solution-phase and gas-phase protein conformations. Finally, infrared multiphoton irradiation of the gas-phase Fe(III)holoproteins releases Fe(III)heme from each of the noncovalently bound Fe(III)heme proteins (myoglobin, cytochrome b5 and cytochrome b5 L47R), but yields Fe(II)heme from the covalently bound heme in cytochrome c.     

Influences of the Hydrophobicity of the Heme－binding Pocket on the Properties and Functions of Cytochrome b5 Mutants

The mutation sites of the four mutants F35Y, P40V, V45E and V45Y of cytochrome b5 are located at the edge of the heme-binding pocket. The solvent accessible areas of the “pocket inte-rior“ of the four mutants and the wild-type cytochrome b5 have been calculated based on their crystal structures at high resolu-tion. The change in the hydrophobicity of the heme-binding pocket resulting from the mutation can be quantitatively de-scribed using the difference of the solvent accessible area of the “pocket interior“ of each mutant from that of the wild-type cy-tochrome b5. The influences of the hydrophobicity of the heme-binding pocket on the protein stability and redox potential are discussed.     

Fe vibrational spectroscopy of myoglobin and cytochrome f

The Fe vibrational density of states (VDOS) has been determined for the heme proteins deoxymyoglobin, metmyoglobin, and cytochrome f in the oxidized and reduced states, using nuclear resonance vibrational spectroscopy (NRVS). For cytochrome f in particular, the NRVS spectrum is compared with multiwavelength resonance Raman spectra to identify those Raman modes with significant Fe displacement. Modes not seen by Raman due to optical selection rules appear in the NRVS spectrum. The mean Fe force constant extracted from the VDOS illustrates how Fe dynamics varies among these four monoheme proteins, and is correlated with oxidation and spin state trends seen in model heme compounds. The protein's contribution to Fe motion is dominant at low frequencies, where coupling to the backbone tightly constrains Fe displacements in cytochrome f, in contrast to enhanced heme flexibility in myoglobin.     

Ligand binding properties of myoglobin reconstituted with iron porphycene: unusual O2 binding selectivity against CO binding

Sperm whale myoglobin, an oxygen storage hemoprotein, was successfully reconstituted with the iron porphycene having two propionates, 2,7-diethyl-3,6,12,17-tetramethyl-13,16-bis(carboxyethyl)porphycenatoiron. The physicochemical properties and ligand bindings of the reconstituted myoglobin were investigated. The ferric reconstituted myoglobin shows the remarkable stability against acid denaturation and only a low-spin characteristic in its EPR spectrum. The Fe(III)/Fe(II) redox potential (-190 mV vs NHE) determined by the spectroelectrochemical measurements was much lower than that of the wild-type. These results can be attributed to the strong coordination of His93 to the porphycene iron, which is induced by the nature of the porphycene ring symmetry. The O2 affinity of the ferrous reconstituted myoglobin is 2600-fold higher than that of the wild-type, mainly due to the decrease in the O2 dissociation rate, whereas the CO affinity is not so significantly enhanced. As a result, the O2 affinity of the reconstituted myoglobin exceeds its CO affinity (M' = K(CO)/K(O2) < 1). The ligand binding studies on H64A mutants support the fact that the slow O2 dissociation of the reconstituted myoglobin is primarily caused by the stabilization of the Fe-O2 sigma-bonding. The IR spectra for the carbon monoxide (CO) complex of the reconstituted myoglobin suggest several structural and/or electrostatic conformations of the Fe-C-O bond, but this is not directly correlated with the CO dissociation rate. The high O2 affinity and the unique characteristics of the myoglobin with the iron porphycene indicate that reconstitution with a synthesized heme is a useful method not only to understand the physiological function of myoglobin but also to create a tailor-made function on the protein.     

Studies on heme binding in myoglobin,hemoglobin, and cytochrome c by ion spray mass spectrometry

The ion spray mass spectra of three representative heme-containing proteins were studied, with an emphasis on results obtained under neutral (pH 7) aqueous conditions. The noncovalently bound heme in myoglobin and hemoglobin may be readily distinguished from the covalently bound heme prosthetic group attached to cytochrome c by using collisioninduced dissociation in the free-jet expansion region of the mass spectrometer as well as in the collision quadrupole with premass selection. The charge state of iron in the expelled heme from myoglobin and hemoglobin appears to be 3+ but 2f for heme expelled from cytochrome c.     

The effectiveness of three multi-component binding models in describing the binary competitive equilibrium adsorption of two cytochrome b(5) mutants

Competitive adsorption isotherms for two conservative surface charge-neutralizing mutants of cytochrome b(5), E11Q and E44Q, previously measured with competitor concentration held constant over the range of the isotherm, were used to test three widely-used multi-component isotherm models. The extended Langmuir-Freundlich, Langmuir and Jovanovic-Freundlich models each adequately described the weaker infinite dilution adsorption of the E44Q protein in the presence of the strong binding E11Q. The extended Langmuir-Freundlich model generally gave the lowest errors at higher concentrations, and the Jovanovic-Freundlich model gave the best fits when using empirically optimized maximal loading values based on multi-component as well as pure-component isotherm data.     

Structural and oxygen binding properties of dimeric horse myoglobin

Myoglobin (Mb) stores dioxygen in muscles, and is a fundamental model protein widely used in molecular design. The presence of dimeric Mb has been known for more than forty years, but its structural and oxygen binding properties remain unknown. From an X-ray crystallographic analysis at 1.05 ? resolution, we found that dimeric metMb exhibits a domain-swapped structure with two extended α-helices. Each new long α-helix is formed by the E and F helices and the EF-loop of the original monomer, and as a result the proximal and distal histidines of the heme originate from different protomers. The heme orientation in the dimer was in the normal mode as in the monomer, but regulated faster from the reverse to normal orientation. The dimer possessed the oxygen binding property, although it exhibited a slightly higher oxygen binding affinity (～1.4 fold) compared to the monomer and showed no cooperativity for oxygen binding. The oxygen binding rate constant (k(on)) of the dimer ((14.0 ± 0.7) × 10(6) M(-1) s(-1)) was similar to that of the monomer, whereas the oxygen dissociation rate constant (k(off)) of the dimer (8 ± 1 s(-1)) was smaller than that of the monomer (12 ± 1 s(-1)). We attribute the similar k(on) values to their active site structures being similar, whereas the faster regulation of the heme orientation and the smaller k(off) in the dimer are presumably due to the slight change in the active site structure and/or more rigid structure compared to the monomer. These results show that domain swapping may be a new tool for protein engineering.     

Discover binding pathways using the sliding binding-box docking approach: application to binding pathways of oseltamivir to avian influenza H5N1 neuraminidase

Drug binding and unbinding are transient processes which are hardly observed by experiment and difficult to analyze by computational techniques. In this paper, we employed a cost-effective method called “pathway docking” in which molecular docking was used to screen ligand-receptor binding free energy surface to reveal possible paths of ligand approaching protein binding pocket. A case study was applied on oseltamivir, the key drug against influenza a virus. The equilibrium pathways identified by this method are found to be similar to those identified in prior studies using highly expensive computational approaches.     

A study of the electron transfer and oxygen binding reactions of myoglobin

Quasi-reversible electron transfer kinetics are reported for sperm whale myoglobin reacting at tin-doped indium oxide electrodes. This reaction was studied by three different electrochemical methods; cyclic voltammetry (CV), single potential step chronoabsorptometry (SPS/CA) and derivative cyclic voltabsorptometry (DCVA). Kinetic parameters were determined from experiments which followed the purification of the protein, the pretreatment of the electrode surface and the removal of dioxygen from the sample solution. A formal heterogeneous electron transfer rate constant, k^0'. of 2.6 (±0.5) × 10⁻⁵ cm s⁻¹ and a transfer coefficient, α, of 0.48 (±0.05) were obtained using SPS/CA. These results are shown to correlate well with those obtained using a second spectroelectrochemical method, DCVA. Anodic and cathodic responses for the heterogeneous electron transfer of myoglobin can be observed using cyclic voltammetry, but these responses are not as reproducible as those that are obtained using the optical methods described here. The advantages of these optical methods over those of cyclic voltammetry are clear from this study.