Effect of a molecular dipole on the ionic strength dependence of a biomolecular rate constant. Identification of the site of reaction

Does Ionic Screening Lower Activation Barriers for Conformational Transitions in Proteins?

In this work, we investigated the kinetics of binding of hen egg-white lysozyme with tri- N-acetylglucosamine in aqueous solutions, at two values of pH, 3.2 and 11, as a function of ionic strength, by a stopped-flow method with tryptophyl fluorescence observation of the transients. We analyzed registered reaction progress curves by employing numerical integration of appropriate chemical master equations. We discriminated between several binding models and established that the process observed in experiments follows a two-step mechanism, composed of four elementary stages: diffusional formation of an encounter complex, dissociation of the encounter complex, conformational transition of the encounter complex to the final complex, and the reverse transformation, i.e., from the final complex to the encounter complex. We evaluated rate constants of these elementary stages and determined their dependencies on solution ionic strength. Regardless of solution pH, rate constants of both forward and reverse conformational transitions increase with an increasing ionic strength. This suggests that ionic screening of intramolecular electrostatic interactions may act to lower the activation barrier for conformational transition in proteins.

Ionic strength effects on the photodegradation reactions of riboflavin in aqueous solution

A study of the effect of ionic strength on the photodegradation reactions (photoreduction and photoaddition) of riboflavin (RF) in phosphate buffer (pH7.0) has been carried out using a specific multicomponent spectrometric method. It has been found that the rates of photodegradation reactions of RF are dependent upon the ionic strength of the solutions at different buffer concentrations. The apparent first-order rate constants (kobs) for the photodegradation of riboflavin at ionic strengths of 0.1-0.5 (0.5M phosphate) lie in the range of 7.35-30.32 × 10(-3) min(-1). Under these conditions, the rate constants for the formation of the major products, lumichrome (LC) by photoreduction pathway, and cyclodehydroriboflavin (CDRF) by photoaddition pathway, are in the range of 3.80-16.03 and 1.70-6.07 × 10(-3) min(-1), respectively. A linear relationship has been observed between log kobs and √μ/1+√μ. A similar plot of log k/ko against √μ yields a straight line with a value of ~+1 for ZAZB showing the involvement of a charged species in the rate determining step. NaCl appears to promote the photodegradation reactions of RF probably by an excited state interaction. The implications of ionic strength on RF photodegradation by different pathways and flavin-protein interactions have been discussed.

Structural investigation on the electrostatic loop of native and mutated SOD1 and their interaction with therapeutic compounds

The electrostatic loop of the native and mutated SOD1 protein with single point mutation in the loop is subjected to MD simulation. The structure and electrostatic properties of the native and mutated loops before/after interacting with small compounds are compared.

Heparin-induced cis- and trans-dimerization modes of the thrombospondin-1 N-terminal domain

Through its interactions with proteins and proteoglycans, thrombospondin-1 (TSP-1) functions at the interface of the cell membrane and the extracellular matrix to regulate matrix structure and cellular phenotype. We have previously determined the structure of the high affinity heparin-binding domain of TSP-1, designated TSPN-1, in association with the synthetic heparin, Arixtra. To establish that the binding of TSPN-1 to Arixtra is representative of the association with naturally occurring heparins, we have determined the structures of TSPN-1 in complex with heparin oligosaccharides containing eight (dp8) and ten (dp10) subunits, by x-ray crystallography. We have found that dp8 and dp10 bind to TSPN-1 in a manner similar to Arixtra and that dp8 and dp10 induce the formation of trans and cis TSPN-1 dimers, respectively. In silico docking calculations partnered with our crystal structures support the importance of arginine residues in positions 29, 42, and 77 in binding sulfate groups of the dp8 and dp10 forms of heparin. The ability of several TSPN-1 domains to bind to glycosaminoglycans simultaneously probably increases the affinity of binding through multivalent interactions. The formation of cis and trans dimers of the TSPN-1 domain with relatively short segments of heparin further enhances the ability of TSP-1 to participate in high affinity binding to glycosaminoglycans. Dimer formation may also involve TSPN-1 domains from two separate TSP-1 molecules. This association would enable glycosaminoglycans to cluster TSP-1.

Kinetics properties of Cu,Zn-superoxide dismutase as a function of metal content

The kinetics of bovine Cu,Zn superoxide dismutase were studied by pulse radiolysis. To ensure the absence of catalytically active free copper, commercially obtained holo-superoxide dismutase was demetallated, and the apo-superoxide dismutase concentrations were determined by isothermal titration calorimetry prior to reconstitution with defined amounts of copper and zinc. The catalytic rate constant was determined as a function of ionic strength over the range of 4-154 mM, and of the copper and zinc content. The catalytic rate constant increases with ionic strength up to (1.5 +/- 0.2) x 10(9) M(-1) s(-1) at an ionic strength of 15 mM, and then decreases. At pH 7 and 50 mM ionic strength, k = (1.2 +/- 0.2) x 10(9) M(-1) s(-1), and at a physiologically relevant ionic strength of 150 mM, it is (0.7 +/- 0.1) x 10 (9) M(-1) s(-1). The effect of ionic strength is ascribed to the inhomogeneous electric field generated by the surface charges of superoxide dismutase. The value of the catalytic rate constant at 50 mM is ca. 2-fold smaller than earlier values reported in the literature. The relationship between copper content and the catalytic rate constant shows that addition of more than a stoichiometric amount of copper cannot be masked efficiently by EDTA. The possibility exists that earlier reported values were based on experiments contaminated with trace amounts of copper.

4-nitroimidazole binding to horse metmyoglobin: evidence for preferential anion binding

The ionization of 4-nitroimidazole to 4-nitroimidazolate was investigated as a function of ionic strength. The apparent pKa varies from 8.99 to 9.50 between 0.001 and 1.0 M ionic strength, respectively, at 25 degrees C. The ionic strength dependence of this ionization is anomalous. The binding of 4-nitroimidazole by horse metmyoglobin was studied between pH 5.0 and 11.5 and as a function of ionic strength between 0.01 and 1.0 M. The association rate constant is pH-dependent, varying from 24 M(-1)s(-1) at pH 5 to a maximum value of 280 M(-1)s(-1) at pH 9.5 and then decreasing to 10 M(-1)s(-1) at pH 11.5 in 0.1 M ionic strength buffers. The dissociation rate constant has a much smaller pH dependence, varying from 0.082 s(-1) at low pH to 0.035 s(-1) at high pH, with an apparent pKa of 6.5. The binding affinity of 4-nitroimidazole to horse metmyoglobin is about 2.5 orders of magnitude stronger than that for imidazole and this increased affinity is attributed to the much slower dissociation rate for 4-nitroimidazole compared to that of imidazole. Although the ionic strength dependence of the binding rate is small and secondary kinetic salt effects can account for the ionic strength dependence of the association rate constant, the pH dependence of the rate constants and microscopic reversibility arguments indicate that the anionic form of the ligand binds more rapidly to all forms of metmyoglobin than does the neutral form of the ligand. However, the spectrum of the complex is similar to model complexes involving neutral imidazole and not imidazolate. The latter observation suggests that the initial metmyoglobin/4-nitroimidazolate complex rapidly binds a proton and the neutral form of the bound ligand is stabilized, probably through hydrogen binding with the distal histidine.

Folding and assembly of dimeric human glutathione transferase A1-1

Glutathione transferases function as detoxification enzymes and ligand-binding proteins for many hydrophobic endogenous and xenobiotic compounds. The molecular mechanism of folding of urea-denatured homodimeric human glutathione transferase A1-1 (hGSTA1-1) was investigated. The kinetics of change were investigated using far-UV CD, Trp20 fluorescence, fluorescence-detected ANS binding, acrylamide quenching of Trp20 fluorescence, and catalytic reactivation. The very early stages of refolding (millisecond time range) involve the formation of structured monomers with native-like secondary structure and exposed hydrophobic surfaces that have a high binding capacity for the amphipathic dye ANS. Dimerization of the monomeric intermediates was detected using Trp fluorescence and occurs as fast and intermediate events. The intermediate event was distinguished from the fast event because it is limited by a preceding slow trans-to-cis isomerization reaction (optically silent in this study). At high concentrations of hFKBP, dimerization is not limited by the isomerization reaction, and only the fast event was detected. The fast (tau = 200 ms) and intermediate (tau = 2.5 s) events show similar urea-, temperature-, and ionic strength-dependent properties. The dimeric intermediate has a partially functional active site ( approximately 20%). Final reorganization to form the native tertiary and quaternary structures occurs during a slow, unimolecular, urea- and ionic strength-independent event. During this slow event (tau = 250 s), structural rearrangements at the domain interface occur at/near Trp20 and result in burial of Trp20. The slow event results in the regain of the fully functional dimer. The role of the C-terminus helix 9 (residues 210-221) as a structural determinant for this final event is proposed.

Electrostatic mutations in iberiotoxin as a unique tool for probing the electrostatic structure of the maxi-K channel outer vestibule

Iberiotoxin (IbTX or alpha-KTx 1.3), a selective, high-affinity blocker of the large-conductance, calcium-activated (maxi-K) channel, exhibits a unique, asymmetric distribution of charge. To test how these charges control kinetics of IbTX binding, we generated five mutants at two positions, K27 and R34, that are highly conserved among other isotoxins. The dissociation and association rate constants, koff and kon, were determined from toxin-blocked and -unblocked durations of single maxi-K channels incorporated into planar lipid bilayers. Equilibrium dissociation constant (Kd) values were calculated from koff/kon. The IbTX mutants K27N, K27Q, and R34N caused large increases in Kd values compared to wild-type, suggesting that the IbTX interaction surface encompasses these residues. A well-established pore-blocking mechanism for IbTX predicts a voltage dependence of toxin-blocked times following occupancy of a potassium binding site in the channel pore. Time constants for block by K27R were approximately 5-fold slower at -20 mV versus +40 mV, while neutralization of K27 relieved the voltage dependence of block. This suggests that K27 in IbTX interacts with a potassium binding site in the pore. Neutralized mutants of K27 and R34, with zero net charge, displayed toxin association rate constants approximately 10-fold slower than wild-type. Association rates for R34N diminished approximately 19-fold when external potassium was increased from 30 to 300 mM. These findings suggest that simple net charge and diffusional processes do not control ingress of IbTX into the channel vestibule.

Role of the electrostatic loop charged residues in Cu,Zn superoxide dismutase

We have expressed and characterized a mutant of Xenopus laevis Cu,Zn superoxide dismutase in which four highly conserved charged residues belonging to the electrostatic loop have been replaced by neutral side chains: Lys120 --> Leu, Asp130 --> Gln, Glu131 --> Gln, and Lys134 --> Thr. At low ionic strength, the mutant enzyme is one of the fastest superoxide dismutases ever assayed (k = 6.7 x 10(9) M(-1) s(-1), at pH 7 and mu = 0.02 M). Brownian dynamics simulations give rise to identical enzyme-substrate association rates for both wild-type and mutant enzymes, ruling out the possibility that enhancement of the activity is due to pure electrostatic factors. Comparative analysis of the experimental catalytic rate of the quadruple and single mutants reveals the nonadditivity of the mutation effects, indicating that the hyperefficiency of the mutant is due to a decrease of the energy barrier and/or to an alternative pathway for the diffusion of superoxide within the active site channel. At physiological ionic strength the catalytic rate of the mutant at neutral pH is similar to that of the wild-type enzyme as it is to the catalytic rate pH dependence. Moreover, mutation effects are additive. These results show that, at physiological salt conditions, electrostatic loop charged residues do not influence the diffusion pathway of the substrate and, if concomitantly neutralized, are not essential for high catalytic efficiency of the enzyme, pointing out the role of the metal cluster and of the invariant Arg141 in determining the local electrostatic forces facilitating the diffusion of the substrate towards the active site.

Effects of single and double mutations in plastocyanin on the rate constant and activation parameters for the rearrangement gating the electron-transfer reaction between the triplet state of zinc cytochrome c and cupriplastocyanin

The unimolecular rate constant for the photoinduced electron-transfer reaction 3Zncyt/pc(II) --> Zncyt+/pc(I) within the electrostatic complex of zinc cytochrome c and spinach cupriplastocyanin is kF. We report the effects on kF of the following factors, all at pH 7.0: 12 single mutations on the plastocyanin surface (Leu12Asn, Leu12Glu, Leu12Lys, Asp42Asn, Asp42Lys, Glu43Asn, Glu59Gln, Glu59Lys, Glu60Gln, Glu60Lys, Gln88Glu, and Gln88Lys), the double mutation Glu59Lys/Glu60Gln, temperature (in the range 273.3-302.9 K), and solution viscosity (in the range 1. 00-116.0 cP) at 283.2 and 293.2 K. We also report the effects of the plastocyanin mutations on the association constant (Ka) and the corresponding free energy of association (DeltaGa) with zinc cytochrome c at 298.2 K. Dependence of kF on temperature yielded the activation parameters DeltaH, DeltaS, and DeltaG. Dependence of kF on solution viscosity yielded the protein friction and confirmed the DeltaG values determined from the temperature dependence. The aforementioned intracomplex reaction is not a simple electron-transfer reaction because donor-acceptor electronic coupling (HAB) and reorganizational energy (lambda), obtained by fitting of the temperature dependence of kF to the Marcus equation, deviate from the expectations based on precedents and because kF greatly depends on viscosity. This last dependence and the fact that certain mutations affect Ka but not kF are two lines of evidence against the mechanism in which the electron-transfer step is coupled with the faster, but thermodynamically unfavorable, rearrangement step. The electron-transfer reaction is gated by the slower, and thus rate determining, structural rearrangement of the diprotein complex; the rate constant kF corresponds to this rearrangement. Isokinetic correlation of DeltaH and DeltaS parameters and Coulombic energies of the various configurations of the Zncyt/pc(II) complex consistently show that the rearrangement is a facile configurational fluctuation of the associated proteins, qualitatively the same process regardless of the mutations in plastocyanin. Correlation of kF with the orientation of the cupriplastocyanin dipole moment indicates that the reactive configuration of the diprotein complex involves the area near the residue 59, between the upper acidic cluster and the hydrophobic patch. Kinetic effects and noneffects of plastocyanin mutations show that the rearrangement from the initial (docking) configuration, which involves both acidic clusters, to the reactive configuration does not involve the lower acidic cluster and the hydrophobic patch but involves the upper acidic cluster and the area near the residue 88.

Imidazole binding to horse metmyoglobin: dependence upon pH and ionic strength

The reaction between metmyoglobin and imidazole has been studied as a function of pH between pH 4.2 and 11.5 and as a function of ionic strength at integral pH values (5 to 10) between 0.001 and 1.0 M ionic strength. The reaction between metmyoglobin and 1-methylimidazole has also been investigated as a function of pH. Comparison of the pH dependence of the association rate constants for the two ligands indicates that the negatively charged imidazolate ion does not contribute to the observed rate of imidazole binding at pH < or = 11.5. At all pH values between pH 4.2 and pH 11.5 the initial complex formed involves the neutral form of bound imidazole. At pH 11.5, the neutral imidazole complex is converted slowly (t1/2 approximately 10 s) into an imidazolate complex. The kinetic data were analyzed according to two mechanisms, one involving the binding of neutral imidazole only and one involving the direct binding of both imidazole and the imidazolium ion to metmyoglobin. Although secondary kinetic salt effects account for the ionic strength dependence of the association rate constant, evidence which indicates that metmyoglobin reacts with imidazole and with the imidazolium ion with similar rates is provided. A self-consistent analysis indicates that the rate constants for imidazole and imidazolium ion binding to metmyoglobin are 350 and 230 M-1 s-1, respectively, at neutral pH and 0.1 M ionic strength. Imidazole can react directly with hydroxymetmyoglobin with a rate of 56 M-1 s-1 at 0.1 M ionic strength, about sixfold slower than binding to aquometmyoglobin. Protonation of a second heme-linked group, thought to be His-97, has little influence on the binding of imidazole but does decrease the rate of imidazolium binding by about eightfold to 29 M-1 s-1 at 0.1 M ionic strength.

Site-directed mutagenesis of rubredoxin reveals the molecular basis of its electron transfer properties

Rubredoxins contain a single non-heme iron atom coordinated by four cysteines. This iron is redox active and confers a role to these proteins in electron transfer chains. The structural features responsible for setting the values of the reduction potential and of the electron self-exchange rate constant have been probed by site-directed mutagenesis. Replacements of the highly conserved residues in positions 8, 10, and 11 (valine, glycine, and tyrosine, respectively) all lead to shifts of the reduction potential, up to 75 mV. These cannot be explained by simple considerations about the physicochemical properties of the substituting side chains but rather indicate that the value of the reduction potential is finely tuned by a variety of interactions. In contrast, the electron self exchange rate constant measured by nuclear magnetic resonance does not vary much, except when a charged residue is included in position 8 or 10, at the surface of the protein closest to the iron atom. Analysis of the data with a model for electrostatic interactions, including both monopolar and dipolar terms, indicates that the presence of a charge in this region not only increases the repulsion between molecules but also affects the electron transfer efficiency of the bimolecular complexes formed. The studies presented constitute a first step toward probing the structural elements modulating the reactivity of the FeS4 unit in a protein and defining the electron transfer active site(s) of rubredoxin.

Very rapid, ionic strength-dependent association and folding of a heterodimeric leucine zipper

Leucine zippers (coiled coils) are dimerization motifs found in several DNA-binding transcription factors. A parallel leucine zipper composed of the acidic chain X1-EYQALEKEVAQLEAENX2-ALEKEVAQLEHEG-amide and the basic chain X1-EYQALKKKVAQLKAKNX2ALKKKVAQLKHKG-amide was designed to study the kinetics of folding of a heterodimeric leucine zipper and to investigate the role of electrostatic attraction between oppositely charged peptide chains to the folding reaction. Each peptide alone did not form a leucine zipper at ionic strength (mu) < 1 M because of electrostatic repulsion between like charges in a homodimer. Therefore, the formation of the heterodimeric leucine zipper could be investigated by simple mixing of acidic and basic chains. To monitor folding, a fluorescent label was located either at the N-terminus (X1 = fluorescein-GGG, X2 = Q) or in the center of the coiled coil (X1 = acetyl, X2 = W). Folding could be described by a simple two-state reaction involving the disordered monomers and the folded heterodimer. The same bimolecular rate constant (k(on)) was observed independent of the location of the fluorescent label, indicating that both fluorescence probes monitored the same reaction. Lowering of the ionic strength increased k(on) from 4 x 10(6) M-1 s-1 (mu = 525 mM) to about 5 x 10(7) M-1 s-1 (mu = 74 mM). When extrapolated to mu = O, k(on) was approximately 10(9) M-1 s-1, which is near the diffusion limit. In contrast, the rate of dissociation depended very weakly on ionic strength; k(off) decreased only by about 2-fold when mu was lowered from 525 to 74 mM. Equilibrium association constants (Ka) of the heterodimeric zippers measured directly and calculated from kinetic constants (Ka = k(on)/k(off) were in good agreement. The observed two-state mechanism, the strong dependence on ionic strength of k(on) but not of k(off), and the nearly diffusion-limited association rate at very low ionic strength point to a folding pathway in which the formation of an electrostatically stabilized dimeric intermediate may be rate-limiting and the subsequent folding to the final dimer is very rapid and follows a "down-hill" free energy landscape. The small increase of k(off) at increasing ionic strength indicates a minor contribution of electrostatics to the stability of the folded leucine zipper.

Interaction of cytochrome c with flavocytochrome b2

Flavocytochrome b2 from Saccharomyces cerevisiae couples L-lactate dehydrogenation to cytochrome c reduction. At 25 degrees C, 0.10 M ionic strength, and saturating L-lactate concentration, the turnover rate is 207 s-1 [per cytochrome c reduced; Miles, C. S., Rouviere, N., Lederer, F., Mathews, F. S., Reid, G. A., Black, M. T., & Chapman, S. K. (1992) Biochem. J. 285, 187-192]. The second-order rate constant for cytochrome c reduction in the pre-steady-state has been determined by stopped-flow spectrophotometry to be 34.8 (+/- 0.9) muM-1 s-1 in the presence of 10 mM L-lactate. This rate constant has been found to be dependent entirely on the rate of complex formation, the electron-transfer rate in the pre-formed complex being in excess of 1000 s-1. Inhibition of the pre-steady-state reduction of cytochrome c by either zinc-substituted cytochrome c or ferrocytochrome c has led to the estimation of a Kd for the catalytically competent complex of 8 microM, and from this the dissociation rate constant of 280 s-1, a value much less than the actual electron-transfer rate. The inhibition observed is only partial which indicates that electron transfer from the 1:1 complex to another cytochrome c can occur and that alternative electron transfer sites exist. The cytochrome c binding site proposed by Tegoni et al. [Tegoni, M., White, S. A., Roussel, A., Mathews, F. S. & Cambillau, C. (1993) Proteins 16, 408-422] has been tested using site-directed mutagenesis. Mutations designed to affect the complex stability and putative electron-transfer pathway had little effect, suggesting that the primary cytochrome c binding site on flavocytochrome b2 lies elsewhere. The combination of tight binding and multiple electron-transfer sites gives flavocytochrome b2 a low K(m) and a high kcat, maximizing its catalytic efficiency. In the steady-state, the turnover rate is therefore largely limited by other steps in the catalytic cycle, a conclusion which is discussed in the preceding paper in this issue [Daff, S., Ingledew, W. J., Reid, G. A., & Chapman, S. K. (1996) Biochemistry 35, 6345-6350].

Calorimetric studies on the interaction of horse ferricytochrome c and yeast cytochrome c peroxidase

The binding of horse ferricytochrome c to yeast cytochrome c peroxidase at pH 6.0 in 8.7 mM phosphate buffer (0.0100 M ionic strength) is characterized by a small, unfavorable enthalpy change (+1.91 +/- 0.16 kcal mol-1) and a large, positive entropy change (+37 +/- 1 eu). The free energy of binding depends strongly upon ionic strength, increasing from -9.01 to -4.51 kcal mol-1 between 0.0100 and 0.200 M ionic strength. The increase in free energy is due solely to the change in entropy over this ionic strength range, with the entropy change decreasing from 37 +/- 1 to 22 +/- 3 eu between 0.0100 and 0.200 M ionic strength. The observed enthalpy change remains constant over the same ionic strength range. At 0.0100 M ionic strength, complex formation is accompanied by the release of 0.54 +/- 0.11 proton, causing a variation in the observed enthalpy of reaction depending upon the buffer. After accounting for proton binding to the buffer, the corrected values for the enthalpy and entropy of binding are +2.84 +/- 0.26 kcal mol-1 and +21 +/- 3 eu, respectively. At 0.05 M ionic strength, the observed change in heat capacity, delta Cp, for the reaction between horse ferricytochrome c and cytochrome c peroxidase is essentially zero, 1.6 +/- 9.6 cal mol-1 K-1. The corrected delta Cp for binding is -28 +/- 10 cal mol-1 K-1 after accounting for proton binding to the buffer. No evidence for formation of a 2:1 horse ferricytochrome c/cytochrome c peroxidase complex was obtained in this study.

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.

Strategies for the study of cytochrome c structure and function by site-directed mutagenesis

The class I cytochromes c have been extensively studied by biochemical and biophysical methods; however, many questions remain concerning the roles of specific amino acids in electron transfer and stability properties. The method of site-directed mutagenesis, which substitutes specific amino acid residues by genetic methods, is ideal for addressing these questions of cytochrome c structure and function. Practical considerations of mutational effects on protein processing and stability will be addressed. The criteria for the selection of mutation sites will be discussed. Examples of site-directed mutagenesis studies, which were designed to elucidate the factors controlling biological electron transfer, protein processing, and protein stability, are given.

A study of the interaction between two proteins, one containing a flavin mononucleotide

Interaction between cytochrome c and flavocytochrome b2 has been studied in presence of 2-p-toluidinylnaphthalene-6-sulfonate (TNS). Affinity of the probe to flavocytochrome b2 increases when the complex between the two proteins is obtained. Binding of TNS increases the fluorescence of flavocytochrome b2 FMN. When the stoichiometry of the complex between the two proteins is reached, TNS looses its affinity and stops binding on the flavocytochrome b2; consequently, FMN fluorescence increase is no more observed. The dissociation constant of the complex was found equal to 0.1 microM. A similar result was obtained for the interaction between cytochrome c and flavodehydrogenase domain. The latter was obtained by proteolysis of flavocytochrome b2.

Ionic-strength-dependence of the oxidation of native and pyridoxal 5'-phosphate-modified cytochromes c by cytochrome c oxidase

The ionic-strength-dependences of the rate constants (log k plotted versus square root of 1) for oxidation of native and pyridoxal 5'-phosphate-modified cytochromes c by three different preparations of cytochrome c oxidase have complex non-linear character, which may be explained on the basis of present knowledge of the structure of the oxidase and the monomer-dimer equilibrium of the enzyme. The wave-type curve (with a minimum and a maximum) for oxidation of native cytochrome c by purified cytochrome c oxidase depleted of phospholipids may reflect consecutively inhibition of oxidase monomers (initial descending part), competition between this inhibition and dimer formation, resulting in increased activity (second part with positive slope), and finally inhibition of oxidase dimers (last descending part of the curve). The dependence of oxidation of native cytochrome c by cytochrome c oxidase reconstituted into phospholipid vesicles is a curve with a maximum, without the initial descending part described above. This may reflect the lack of pure monomers in the vesicles, where equilibrium is shifted to dimers even at low ionic strength. Subunit-III-depleted cytochrome c oxidase does not exhibit the maximum seen with the other two enzyme preparations. This may mean that removal of subunit III hinders dimer formation. The charge interactions of each of the cytochromes c (native or modified) with the three cytochrome c oxidase preparations are similar, as judged by the similar slopes of the linear dependences at I values above the optimal one. This shows that subunit III and the phospholipid membrane do not seem to be involved in the specific charge interaction of cytochrome c oxidase with cytochrome c.

Electrostatic and steric control of electron self-exchange in cytochromes c, c551, and b5

The ionic strength dependence of the electron self-exchange rate constants of cytochromes c, c551, and b5 has been analyzed in terms of a monopole-dipole formalism (van Leeuwen, J.W. 1983. Biochim. Biophys. Acta. 743:408-421). The dipole moments of the reduced and oxidized forms of Ps. aeruginosa cytochrome c551 are 190 and 210 D, respectively (calculated from the crystal structure). The projections of these on the vector from the center of mass through the exposed heme edge are 120 and 150 D. For cytochrome b5, the dipole moments calculated from the crystal structure are 500 and 460 D for the reduced and oxidized protein; the projections of these dipole moments through the exposed heme edge are -330 and -280 D. A fit of the ionic strength dependence of the electron self-exchange rate constants gives -280 (reduced) and -250 (oxidized) D for the center of mass to heme edge vector. The self-exchange rate constants extrapolated to infinite ionic strength of cytochrome c, c551, and b5 are 5.1 x 10(5), 2 x 10(7), and 3.7 x 10(5) M-1 s-1, respectively. The extension of the monopole-dipole approach to other cytochrome-cytochrome electron transfer reactions is discussed. The control of electron transfer by the size and shape of the protein is investigated using a model which accounts for the distance of the heme from each of the surface atoms of the protein. These calculations indicate that the difference between the electrostatically corrected self-exchange rate constants of cytochromes c and c551 is due only in part to the different sizes and heme exposures of the two proteins.

Electrostatic interactions of 4-carboxy-2,6-dinitrophenyllysine-modified cytochromes c with physiological and non-physiological redox partners

An analysis of the effect of electrostatic properties of 4-carboxy-2,6-dinitrophenyllysine (CDNP-lysine) cytochromes c on their reactions with strongly and weakly binding redox partners is given. For strongly binding systems (cytochrome-c oxidase, cytochrome-c reductase, sulphite oxidase and yeast cytochrome-c peroxidase) the magnitude of the dipole moments of the CDNP cytochromes c determines their relative reactivities. For weakly binding redox agents, such as hexacyanoferrate(III), cobalt(III)tris(1,10-phenanthroline), azurin and plastocyanin, the electrostatic potential at the haem edge accounts for the greater part of the relative activities. Relative rate data were obtained from the literature. It is concluded that the dipole moment of native cytochromes c may account for an approx. 50-fold increase in the efficiency of its physiological activity towards membrane-bound enzymes. A correction on a formula to describe the contribution of a molecular dipole moment to the ionic strength dependence of a bimolecular rate constant (Koppenol, W. H. (1980) Biophys. J. 29, 493-508) leads to an equation nearly identical to that obtained by Van Leeuwen et al. (Van Leeuwen, J.W., Mofers, F.J.M. and Verrman, E.C.I. (1981) Biochim. Biophys. Acta 635, 434-439).

Electron self-exchange in Pseudomonas cytochromes

Electron self-exchange has been measured by an NMR technique for cytochromes c551 from Pseudomonas aeruginosa and Pseudomonas stutzeri. The rate for P. aeruginosa cyt c551 is 1.2 x 10(7) M-1 s-1 at 40 degrees C in 50 mM phosphate at pH 7. For P. stutzeri, under the same conditions, the rate is 4 x 10(7) M-1 s-1. For both cytochromes, the rate was independent of ionic strength up to 0.5 M in added NaC1, the enthalpy of activation was 20 +/- 4 kcal mol-1, and the entropy of activation was 38 +/- 10 cal mol-1 deg-1.

Cytochrome b2, an electron carrier between flavocytochrome b2 and cytochrome c. Rapid kinetic characterization of the electron-transfer parameters with ionic-strength-dependence

The oxidation-reduction properties of free cytochrome b2 isolated by controlled proteolysis from flavocytochrome b2, i.e. the flavodehydrogenase-bound cytochrome b2, were investigated by using stopped-flow spectrophotometry. The rapid kinetics of the reduction of cytochrome b2 by flavocytochrome b2 in the presence of L-lactate are reported. The self-exchange rate constant between reduced cytochrome b2 bound to the flavodehydrogenase and free cytochrome b2 was determined to be 10(5) M-1 X S-1 at 5 degrees C, I 0.2 and pH 7.0. The specific electron-transfer reaction between reduced cytochrome b2 and cytochrome c was also studied, giving an apparent second-order rate constant of 10(7) M-1 X S-1 at 5 degrees C, I 0.2 and pH 7.0. This electron-exchange rate is slightly modulated by ionic strength, following the Debye-Hückel relationship with a charge factor Z1Z2 = -1.9. Comparison of these data with those for the reduction of cytochrome c by flavodehydrogenase-bound cytochrome b2 [Capeillère-Blandin (1982) Eur. J. Biochem. 128, 533-542] leads to the conclusion that the intramolecular electron exchange between haem b2 and haem c within the reaction complex occurs at a rate very similar to that determined experimentally in presence of the flavodehydrogenase domain. The low reaction rate observed with free cytochrome b2 is ascribed to the low stability of the reaction complex formed between free cytochrome b2 and cytochrome c.

Influence of buffer composition, membrane lipids and proteases on the kinetics of reconstituted cytochrome-c oxidase from bovine liver and heart

Isolated cytochrome-c oxidases from bovine heart and liver were reconstituted in liposomes with asolectin and the kinetics of cytochrome c oxidation were measured under various uncoupled conditions. With 40 mM KCl, 10 mM Hepes, pH 7.4, the liver enzyme showed a higher Vmax in the polarographic but a lower Vmax in the photometric assay. With 125 mM phosphate buffer at pH 6.0 both enzymes revealed identical kinetics. Reconstitution with pure phosphatidylcholine leads to a low activity, which is specifically stimulated for the heart enzyme by inclusion of 10% cardiolipin. Proteoliposomes of both enzymes prepared with asolectin have a high activity, which is unaffected by cardiolipin. Exchanging the intraliposomal buffer, Hepes, for phosphate causes an opposite change of the Vmax and a similar change of the Km for both enzymes suggesting a conformational change of the extraliposomal binding domain for cytochrome c through the membrane. Proteases change the kinetics of both enzymes, but to a different degree. The data indicate a complex and tissue-specific influence of nucleus-coded subunits on the catalytic activity of cytochrome-c-oxidase.

Is cytochrome c reactivity determined by dipole moment or by local charges?

A criticism of a recent paper by M. Fragata and F. Bellemare (Biophys. Chem. 15 (1982) 111) is presented. These authors developed a model of polarity-dependent ferrocytochrome c oxidation which is shown to be incorrect. It fails to show that the use of the 'overal dipole moment' is likely to be unreliable, and that reactivity is best explained by a polarity effect on the dipole of the haem of cytochrome c.

Toward a multistep mechanism of cytochrome c reactivity. Answer to a comment

Koppenol's rejection (Biophys. Chem. 18 (1983) 203) of a model of polarity-dependent ferrocytochrome c oxidation (M. Fragata and F. Bellemare, Biophys. Chem. 15 (1982) 111) places emphasis on the role of the protein surface charges in reactivity but is at the same time too restrictive as it neglects largely the polarity (dielectric constant) of the aqueous and hydrophobic interfaces of the exposed heme edge and the inner cleft (heme crevice) of cytochrome c which appear to be the oxidation-reduction sites. It is suggested that a more general model should take into account (i) a recognition (or diffusion) step where the distance travelled by cytochrome c at the membrane surface and/or the Brownian displacements in the bulk solution are greatly influenced by ionic strength, and (ii) a redox step where low polarity effects prevail with concomitant weakening of ionic activity.

Transient kinetics of the one-electron transfer reaction between reduced flavocytochrome b2 and oxidized cytochrome c. Evidence for the existence of a protein complex in the reaction

The one-electron transfer reaction from reduced flavocytochrome b2 (fully reduced by three electron equivalents) to ferricytochrome c, both purified from the yeast Hansenula anomala, has been studied using stopped-flow spectrophotometry in the course of a single turnover, for reactants initially mixed in a heme molar ratio equal to one. The cytochrome c reduction proceeded to completion through an apparently first-order process. Depending on the experimental conditions (concentrations and or ionic strength), the reduction is of second-order or first-order character. To interpret these kinetic results computer simulation studies have been performed based on a kinetic scheme involving, besides the formation of a complex before the electron transfer step, intramolecular electron transfer steps within flavocytochrome b2 to maintain the concentration of the specific electron donor center, the reduced cytochrome b2. As far as the cytochrome c reduction rate constant, ka, and its variations were concerned the simulated data showed that this complicated scheme could approximate a mechanism which is by far the simplest, involving only the two former steps. Such a scheme accounts firstly for the hyperbolic dependence of the rate of reduction of cytochrome c, ka, upon reductant concentrations which had provided clear evidence for the kinetic existence of a complex in the reaction pathway. At 5 degrees C the rate constant for the electron transfer is 380 s-1 with an activation energy of 13.8kJ mol-1 (3.3 kcal mol-1). Secondly it predicts the observed variations of ka with ionic strength and provides estimates of the rate constants of the binding step.

The reaction of cytochrome aa3 with (porphyrin) cytochrome c as studied by pulse radiolysis

(1) Using the pulse-radiolysis and stopped-flow techniques, the reactions of iron-free (porphyrin) cytochrome c and native cytochrome c with cytochrome aa3 were investigated. The porphyrin cytochrome c anion radical (generated by reduction of porphyrin cytochrome c by the hydrated electron) can transfer its electron to cytochrome aa3. The bimolecular rate constant for this reaction is 2 x 10(7) M-1 . s-1 (5 mM potassium phosphate, 0.5% Tween 20, pH 7.0, 20 degrees C). (2) The ionic strength dependence of the cytochrome c-cytochrome aa3 interaction was measured in the ionic strength range between 40 and 120 mM. At ionic strengths below 30 mM, a cytochrome c-cytochrome aa3 complex is formed in which cytochrome c is no longer reducible by the hydrated electron. A method is described by which the contributions of electrostatic forces to the reaction rate can be determined. (3) Using the stopped-flow technique, the effect of the dielectric constant (epsilon) of the reaction medium on the reaction of cytochrome C with cytochrome aa3 was investigated. With increasing epsilon the second-order rate constant decreased.

Dielectric constant dependence of biological oxidation-reduction. 1. A model of polarity-dependent ferrocytochrome c oxidation

A theoretical model for the effect of the dielectric constant (c) of the solvent medium on ferrocytochrome c oxidation by ferricyanide is developed to account for the observed variations of the rate constant (k) of reactions in aqueous binary mixtures with alcohols (less than 5-10 mol% ethanol and propranolol). A correlation between k and c is found if ln k is expressed as a function of the Kirkwood parameter (c-1)(2c+1). The results of calculations indicate that the use of the 'overall dipole moment' of cytochrome c in oxidoreduction studies is likely to be unreliable. Instead, the decrease in k in alcohol/water mixtures is best explained--in conformity with Onsager's theory of the reaction field--by a polarity effect on the dipole moment of the cytochrome c heme upon diffusion of the polar solvent molecules into the low dielectric constant heme crevice.

Molecular interpretation of kinetic-ionic strength effects

A recent and important approach to investigating electron transfer mechanisms of redox proteins has been through kinetic-ionic strength studies. There is, however, significant controversy as to whether such studies (1) yield information regarding the charge (or location) of the electron transfer site or (2) more simply reflect the influence of net or overall protein charge on the electrostatic interactions. A critical analysis using different theoretical approaches is made of our recent work and of the bulk of the published non-physiological small molecule-protein and protein-protein kinetic ionic strength studies; it is concluded that (1) the approximated Bronsted-Debye-Huckel equation can not be used at all for protein redox reactions, (2) irrespective of the theoretical approaches discussed, such studies do not provide information regarding the charge of the electron transfer site, (3) it is the net charge of the reactants that control the electrostatic interactions, (4) both the equation derived by Wherland and Gray and the full Bronsted-Debye-Huckel equation provide reasonably good approximations of net protein charge, (5) pH changes quantitatively modulate net protein charge, and (6) thus, protein redox rates need to be electrostatically corrected if relevant interpretations of kinetic-ionic strength experiments are to be made.

Laser flash photolysis studies of electron transfer between semiquinone and fully reduced free flavins and horse heart cytochrome c

Laser flash photolysis has been used to determine the second-order rate constants for the reduction of horse heart cytochrome c by the semiquinone and fully reduced forms of various flavin analogs. We find that substitution in the dimethylbenzene ring of the flavin causes appreciable changes in the rate constants, whereas substitutions at the N-10 position do not. Placing a charged phosphate group in the N-10 ribityl side chain leads to only small ionic strength effects on the rate constants, whereas a charged group attached to the dimethylbenzene ring produces a large ionic strength effect. These results can be accounted for by assuming that a productive collision between flavin and cytochrome involves an orientation that positions the aromatic ring--N-5 region of the flavin toward the heme crevice and the N-10--pyrimidine ring region away from it. Our observations have implications for mechanistic understanding of biological electron transfer reactions and are discussed in this context.