EPR redox study of cytochromec3fromDesulfovibrio vulgarisMiyazaki

Very Low-Frequency Broadband Electron Paramagnetic Resonance Spectroscopy of Metalloproteins

A previously developed spectrometer for broadband electron paramagnetic resonance (EPR) spectroscopy of dilute randomly oriented systems has been considerably modified to extend the frequency reach down to the hundred MHz range and to boost concentration sensitivity by 1 to 2 orders of magnitude. The instrument is now suitable for the study of biological systems in particular metalloproteins. As a proof of concept, examples from the class of low-spin ferric hemoproteins are studied in terms of frequency-dependent changes in their EPR spectra. Mono-heme cytochrome c EPR is determined by g-strain over a wide frequency range, whereas a combination of unresolved ligand hyperfine interaction and concentration-dependent intermolecular dipolar interaction becomes dominant at very low frequencies. In the four heme containing cytochrome c₃, g-strain combines with intramolecular dipolar interaction over the full-studied frequency range of 0.23–12.0 GHz. It is concluded that the point-dipole approach is inappropriate to describe magnetic interactions between low-spin ferric heme systems and that a body of literature on redox interactions in multi-heme proteins will be affected by this conclusion.

Exploring Coupled Redox and pH Processes with a Force-Field-Based Approach: Applications to Five Different Systems

Coupled redox and pH-driven processes are at the core of many important biological mechanisms. As the distribution of protonation and redox states in a system is associated with the pH and redox potential of the solution, having efficient computational tools that can simulate under these conditions becomes very important. Such tools have the potential to provide information that complement and drive experiments. In previous publications we have presented the implementation of the constant pH and redox potential molecular dynamics (C(pH,E)MD) method in AMBER and we have shown how multidimensional replica exchange can be used to significantly enhance the convergence efficiency of our simulations. In the current work, after an improvement in our C(pH,E)MD approach that allows a given residue to be simultaneously pH- and redox-active, we have employed our methodologies to study five different systems of interest in the literature. We present results for capped tyrosine dipeptide, two maquette systems containing one pH- and redox-active tyrosine (α₃Y and peptide A), and two proteins that contain multiple heme groups (diheme cytochrome c from Rhodobacter sphaeroides and Desulfovibrio vulgaris Hildenborough cytochrome c₃). We show that our results can provide new insights into previous theoretical and experimental findings by using a fully force-field-based and GPU-accelerated approach, which allows the simulations to be executed with high computational performance.

Redox-coupled conformational alternations in cytochrome c(3) from D. vulgaris Miyazaki F on the basis of its reduced solution structure

Heteronuclear NMR spectroscopy was performed to determine the solution structure of (15)N-labeled ferrocytochrome c(3) from Desulfovibrio vulgaris Miyazaki F (DvMF). Although the folding of the reduced cytochrome c(3) in solution was similar to that of the oxidized one in the crystal structure, the region involving hemes 1 and 2 was different. The redox-coupled conformational change is consistent with the reported solution structure of D. vulgaris Hildenborough ferrocytochrome c(3), but is different from those of other cytochromes c(3). The former is homologous with DvMF cytochrome c(3) in amino acid sequence. Small displacements of hemes 1 and 2 relative to hemes 3 and 4 were observed. This observation is consistent with the unusual behavior of the 2(1)CH(3) signal of heme 3 reported previously. As shown by the (15)N relaxation parameters of the backbone, a region between hemes 1 and 2 has more flexibility than the other regions. The results of this work strongly suggest that the cooperative reduction of hemes 1 and 2 is based on the conformational changes of the C-13 propionate of heme 1 and the aromatic ring of Tyr43, and the interaction between His34 and His 35 through covalent and coordination bonds.

Further characterization of the two tetraheme cytochromes c3 from Desulfovibiro africanus: nucleotide sequences, EPR spectroscopy and biological activity

The genes encoding the basic and acidic tetraheme cytochromes c3 from Desulfovibrio africanus have been sequenced. The corresponding amino acid sequences of the basic and acidic cytochromes c3 indicate that the mature proteins consist of a single polypeptide chain of 117 and 103 residues, respectively. Their molecular masses, 15102 and 13742 Da, respectively, determined by mass spectrometry, are in perfect agreement with those calculated from their amino acid sequences. Both D. africanus cytochromes c3 are synthesized as precursor proteins with signal peptides of 23 and 24 residues for the basic and acidic cytochromes, respectively. These cytochromes c3 exhibit the main structural features of the cytochrome c3 family and contain the 16 strictly conserved cysteine + histidine residues directly involved in the heme binding sites. The D. africanus acidic cytochrome c3 differs from all the other homologous cytochromes by its low content of basic residues and its distribution of charged residues in the amino acid sequence. The presence of four hemes per molecule was confirmed by EPR spectroscopy in both cytochromes c3. The g-value analysis suggests that in both cytochromes, the angle between imidazole planes of the axial histidine ligands is close to 90 degrees for one heme and much lower for the three others. Moreover, an unusually high exchange interaction (approximately 10[-2] cm[-1]) was evidenced between the highest potential heme (-90 mV) and one of the low potential hemes in the basic cytochrome c3. The reactivity of D. africanus cytochromes c3 with heterologous [NiFe] and [Fe] hydrogenases was investigated. Only the basic one interacts with the two types of hydrogenase to achieve efficient electron transfer, whereas the acidic cytochrome c3 exchanges electrons specifically with the basic cytochrome c3. The difference in the specificity of the two D. africanus cytochromes c3 has been correlated with their highly different content of basic and acidic residues.

Regulation of the redox order of four hemes by pH in cytochrome c3 from D. vulgaris Miyazaki F

The assignment of 1H-NMR signals of the heme methyl and propionate groups of cytochrome c3 of D. vulgaris Miyazaki F was performed. The heme assignment was revised for hemes 2 and 3 (sequential heme numbering). Namely, heme 4 is mainly reduced at first with hemes 1, 2 and 3 following it in this order. The p2H titration of heme methyl signals in four macroscopic oxidation states was performed in the p2H range of 5.2 to 9.0. While the heme methyl resonances in the fully oxidized state showed just small changes with p2H, most resonances in the intermediate oxidation states revealed clear p2H dependence. In particular, the methyl resonances of heme 1 shifted significantly in the acidic region. Then, the chemical shifts of beta-CH2 (next to the carboxyl group) of all propionate groups in the fully oxidized state were observed at various p2H in the range of 4.5 to 9.0. Only the propionate group at C-13 (IUPAC-IUB nomenclature) of heme 1 showed a clear change in this p2H range, its titration curve being similar to those of the methyl resonances of heme 1 in the intermediate oxidation states. pKa of the propionate group was 5.95 +/- 0.05. Analysis of the microscopic formal redox potentials was carried out for the observations at p2H 5.2, 7.1 and 9.0. The redox potentials of heme 1 showed the most remarkable p2H dependence, resulting in the change of the order of the redox potentials of four hemes. A significant change was also found in the interacting potential between hemes 1 and 2. In the light of the p2H-titration experiments, the propionate at C-13 of heme 1 was identified as the most plausible ionizable group responsible for the p2H dependence of microscopic redox potentials of heme 1 in the acidic region.

Catalytic properties of adenylylsulfate reductase from Desulfovibrio vulgaris Miyazaki

Adenylylsulfate reductase (EC 1.8.99.2) isolated from Desulfovibrio vulgaris Miyazaki catalyzes electron transfer from dihydroflavin coenzyme (FADH2, FMNH2, or dihydroriboflavin) to adenylyl sulfate (APS), and catalyzes flavin-mediated oxidation of ferrocytochrome c3 with APS. The reaction with FAD as an electron mediator was markedly stimulated in the presence of menadione. Km of the enzyme was about 0.015 mM for riboflavin and FAD in the presence of menadione. Free flavin coenzyme was found to be the normal cellular constituent. These observations suggested that free flavin coenzyme may be a physiological electron carrier for APS reductase, and the enzyme may be called AMP, sulfite:flavin oxidoreductase. Km (APS) of this enzyme is lower than 1 microM. The enzyme is not inhibited by ATP and GTP, but was inhibited by AMP and sulfite. Its extremely low Km (APS) enables this enzyme to reduce any traces of cytosolic APS which is present only at micromolar concentration, and inhibition by sulfite makes this organism utilize an energetically favorable electron acceptor, sulfite, preferentially over APS which is produced from sulfate at the cost of ATP.

Individual redox characteristics and kinetic properties of the hemes in cytochromes c3: new methods of investigation

The elucidation of the role of the four hemes in cytochromes c3 requires several complementary approaches. The measurements and the assignment of the redox potentials resort to magnetic spectroscopies, EPR and NMR, which are able to discriminate the hemes. The origin of the differences between the redox properties of the hemes can be studied by comparing their thermodynamic parameters delta S and delta H, as measured by the temperature dependence of their individual potentials. Lastly, the available data concerning the electron exchange between cytochromes c3 and their redox partners can be analysed through a detailed kinetic model which provides important information on the role of the different hemes.

Reduction kinetics of the four hemes of cytochrome c3 from Desulfovibrio vulgaris by flash photolysis

The reduction of the tetraheme cytochrome c3 (from Desulfovibrio vulgaris, strains Miyazaki F and Hildenbourough) by flavin semiquinone and reduced methyl viologen follows a monophasic kinetic profile, even though the four hemes do not have equivalent reduction potentials. Rate constants for reduction of the individual hemes are obtained subsequent to incrementally reducing the cytochrome by phototitration. The dependence of each rate constant on the reduction potential difference between the heme and the reductant can be described by outer sphere electron transfer theroy. Thus, the very low reduction potentials of the cytochrome c3 hemes compensate for the very large solvent accessibility of the hemes. The relative rate constants for electron transfer to the four hemes of cytochrome c3 are consistent with the assignments of reduction potential to hemes previously made by Park et al. (Park, J.-S., Kano, K., Niki, S. and Akutsu, H. (1991) FEBS Lett. 285, 149-151) using NMR techniques. The ionic strength dependence of the observed rate constant for reduction by the methyl viologen radical cation indicates that ionic strength substantially alters the structure and/or the heme reduction potentials of the cytochrome. This result is confirmed by reduction with a neutral flavin species (5-deazariboflavin semiquinone) in which the reactivity of the highest potential heme decreases and the reactivity of the lowest potential heme increases at high (500 mM) ionic strength, and by the sensitivity of heme methyl resonances to ionic strength as observed by 1H-NMR. These unusual ionic strength-dependent effects may be due to a combination of structural changes in the cytochrome and alterations of the electrostatic fields at elevated ionic strengths.

Mössbauer characterization of the tetraheme cytochrome c3 from Desulfovibrio baculatus (DSM 1743). Spectral deconvolution of the heme components

Mössbauer spectroscopy was used to study the tetraheme cytochrome c3 from Desulfovibrio baculatus (DSM 1743). Samples with different degrees of reduction were prepared using a redoxtitration technique. In the reduced cytochrome c3, all four hemes are reduced and exhibit diamagnetic Mössbauer spectra typical for low-spin ferrous hemes (S = 0). In the oxidized protein, the hemes are low-spin ferric (S = 1/2) and exhibit overlapping magnetic Mössbauer spectra. A method of differential spectroscopy was applied to deconvolute the four overlapping heme spectra and a crystal-field model was used for data analysis. Characteristic Mössbauer spectral components for each heme group are obtained. Hyperfine and crystal-field parameters for all four hemes are determined from these deconvoluted spectra.

Full assignment of heme redox potentials of cytochrome c3 of D. vulgaris Miyazaki F by 1H-NMR

Site-specific heme assignment of the 1H-NMR spectrum of cytochrome c3 of D. vulgaris Miyazaki F, a tetraheme protein, was established. The major reduction of the heme turned out to take place in the order of hemes I, III, IV and II (numbering in the crystal structure). The hemes with the smallest and greatest solvent accessibility were reduced at the highest and lowest potentials on average, respectively. A cooperative interheme interaction was attributed to a pair of the closest hemes, namely, hemes III and IV. This assignment can provide the physiochemical basis for the elucidation of electron transfer of this protein.

Structural and functional approach toward a classification of the complex cytochrome c system found in sulfate-reducing bacteria

Following the discovery of the tetraheme cytochrome c3 in the strict anaerobic sulfate-reducing bacteria (Postgate, J.R. (1954) Biochem. J. 59, xi; Ishimoto et al. (1954) Bull. Chem. Soc. Japan 27, 564-565), a variety of c-type cytochromes (and others) have been reported, indicating that the array of heme proteins in these bacteria is complex. We are proposing here a tentative classification of sulfate- (and sulfur-) reducing bacteria cytochromes c based on: number of hemes per monomer, heme axial ligation, heme spin state and primary structures (whole or fragmentary). Different and complementary spectroscopic tools have been used to reveal the structural features of the heme sites.

Single-crystal electron paramagnetic resonance study of cytochrome c3 from Desulfovibrio desulfuricans Norway Strain. Assignment of the heme midpoint redox potentials

A single crystal of cytochrome c3 from Desulfovibrio desulfuricans Norway is studied by electron paramagnetic resonance at low temperature. The orientation of the principal axis corresponding to the largest g value is determined for the 12 heme groups in the crystal unit cell. The comparison of these directions to the normals to the heme planes, determined from the crystallographic data at 2.5 A resolution, gives strong evidence for the following assignment of the midpoint redox potentials to the heme groups H1 to H4, defined in the three-dimensional structure: -150 mV is assigned to H3, -300 mV to H4, -330 mV to H1 and -355 mV to H2. This assignment is in agreement with a partial correspondence previously established from an independent study performed on cytochrome c3 in solution.

Estimation of microscopic redox potentials of a tetraheme protein, cytochrome c3 of Desulfovibrio vulgaris, Miyazaki F, and partial assignments of heme groups

The microscopic formal redox potentials of a tetraheme protein, cytochrome c3 from Desulfovibrio vulgaris, Miyazaki F, were estimated from the chemical shifts of the heme methyl signals in its 1H NMR spectrum. All chemical shifts in the five macroscopic oxidation states were determined for eight of the heme methyl protons by the saturation-transfer method. The electron-distribution probability at each heme in each oxidation state was estimated directly from the chemical shifts. To minimize errors due to interheme pseudocontact contributions, the average electron-distribution probability was used for calculation of the microscopic formal redox potentials. By introducing interacting potentials, 32 parameters were reduced to 10. The 10 parameters were determined analytically from the 9 independent electron-distribution probabilities and 2 macroscopic formal redox potentials. The results showed the presence of a strong positive interaction between a pair of particular hemes. The microscopic formal redox potential changes dramatically with the extent of reduction because of the intramolecular interheme interactions. NMR signals of two hemes were assigned to particular hemes in the crystal structures by nuclear Overhauser effect experiments. The results showed that the hemes with the highest and lowest redox potentials in the one-electron reduction process correspond to hemes I and IV in the crystal structure.

EPR study of the redox interactions in cytochrome c3 from Desulfovibrio vulgaris Miyazaki

We present a new examination of the EPR redox titration data for the tetraheme cytochrome c3 from Desulfovibrio vulgaris Miyazaki. Our analysis includes the contribution of the interaction potentials between the four redox sites and is based on the model previously developed for the study of cytochrome c3 from Desulfovibrio desulfuricans Norway. We observed, as for D. desulfuricans Norway cytochrome c3, that the conformation of the heme with the lowest redox potential, heme 4, is sensitive to the redox state of the heme with the highest potential, heme 1. However in D. vulgaris Miyazaki cytochrome c3 spectral simulations show that heme 4 is present in two conformational states which interconvert partially when heme 1 is reduced. The sets of redox parameters which satisfy the fitting procedure of the titration curves are in the following domain: -250 mV less than or equal to e41 less than or equal to -220 mV, -325 mV less than or equal to e2 less than or equal to -320 mV, -335 mV less than or equal to e3 less than or equal to -330 mV, -360 mV less than or equal to e4 less than or equal to -355 mV, -5 mV less than or equal to I12 less than or equal to 20 mV, -10 mV less than or equal to I13 less than or equal to 5 mV, -15 mV less than or equal to I23 less than or equal to -5 mV, -15 mV less than or equal to I24 less than or equal to -10 mV, -25 mV, less than or equal to I34 less than or equal to -15 mV. As in D. desulfuricans Norway cytochrome c3 the interactions are moderate. Simple electrostatic considerations suggest that these moderate values could be related to the large accessibility of the hemes to the solvent. Our work does not confirm the existence of a cooperative interaction between heme 2 and heme 3 which has been proposed on the basis of electrochemical measurements.

Characterization and complete amino acid sequence of ferredoxin II from Desulfovibrio vulgaris Miyazaki

One of 2 ferredoxins, Fd II, purified from Desulfovibrio vulgaris Miyazaki (DvM) has been characterized and its complete amino acid sequence established. Fd II is composed of 63 amino acid residues and contains 7 cysteinyl residues but has only 4 iron atoms in an iron-sulfur cluster of a standard redox potential of -405 mV. The arrangement of cysteinyl residues in the protein suggests that some cysteinyl residues are not directly involved in ligation to the iron-sulfur cluster. Homology is recognized among Fd II (DvM), Fd I (D. desulfuricans Norway), Fd I (D. africanus Benghazi), and Fd (D. gigas). Although Fd I and Fd II in DvM are poorly homologous, the C-terminal half of Fd I is fairly homologous to the N-terminal half of Fd II. Fd II is 40% as effective as Fd I as an electron carrier for pyruvate dehydrogenase coupled with hydrogenase and cytochrome c3.

EPR determination of interaction redox potentials in a multiheme cytochrome: cytochrome c3 from Desulfovibrio desulfuricans Norway

In cytochromes c3 which contain four hemes per molecule, the redox properties of each heme may depend upon the redox state of the others. This effect can be described in terms of interaction redox potentials between the hemes and must be taken into account in the characterization of the redox properties of the molecule. We present here a method of measurement of these interactions based on the EPR study of the redox equilibria of the protein. The microscopic and macroscopic midpoint potentials and the interaction potentials are deduced from the analysis of the redox titration curves of the intensity and the amplitude of the EPR spectrum. This analysis includes a precise simulation of the spectrum of the protein in the oxidized state in order to determine the relative contribution of each heme to the spectral amplitude. Using our method on cytochrome c3 from D. desulfuricans Norway, we found evidence for the existence of weak interaction potentials between the hemes. The three interaction potentials which have been measured are characterized by absolute values lower than 20 mV in contrast with the values larger than 40-50 mV which have been reported for cytochrome c3 from D. gigas. Simulations of the spectra of samples poised at different potentials indicate a structural modification of the heme with the most negative potential during the first step of reduction. The correspondence between the redox sites as characterized by the EPR potentiometric titration and the hemes in the tridimensional structure is discussed.

Assignment of individual heme EPR signals of Desulfovibrio baculatus (strain 9974) tetraheme cytochrome c3. A redox equilibria study

An EPR redox titration was performed on the tetraheme cytochrome c3 isolated from Desulfovibrio baculatus (strain 9974), a sulfate-reducer. Using spectral differences at different poised redox states of the protein, it was possible to individualize the EPR g-values of each of the four hemes and also to determine the mid-point redox potentials of each individual heme: heme 4 (-70 mV) at gmax = 2.93, gmed = 2.26 and gmin = 1.51; heme 3 (-280 mV) at gmax = 3.41; heme 2 (-300 mV) at gmax = 3.05, gmed = 2.24 and gmin = 1.34; and heme 1 (-355 mV) at gmx = 3.18. A previously described multi-redox equilibria model used for the interpretation of NMR data of D. gigas cytochrome c3 [Santos, H., Moura, J.J.G., Moura, I., LeGall, J. & Xavier, A. V. (1984) Eur. J. Biochem. 141, 283-296] is discussed in terms of the EPR results.