Exploring Electron Transfer Between Heme Proteins of Cytochrome c Super Family in Biosensors: A Molecular Mechanics Study

Abstract

Electron transfer between heme proteins with mediators plays an important role in the fabrication of sensitive bio-nano sensors. Heme protein Cytochrome c (pdb code—1HRC) was chosen as the mediator with Cytochrome c′ (pdb code—1A7V) as the probe protein for our investigation on the electron transfer process. We used the software GRAMM, HEX, and MACRODOX to build the protein complex with further evaluation by GROMACS potential. After molecular mechanics refinement by GROMACS the protein complexes were evaluated in terms of the following criteria: Hydrophobic packing, proximity of the hemes, hydrogen bonds, enthalpy and entropy of binding. The free energy was calculated for each complex to derive the feasible stable models. The combined electron transport of the chosen geometric models was evaluated to choose the possible models. Electrostatic potential was calculated using the program APBS around the heme in the presence and absence of other proteins. From our studies, we derived multiple feasible models and possible electronic path. These studies helped us to understand the relay mechanism between the two proteins and to design mutant proteins by rational site directed mutagenesis to enhance the redox potential and thereby improving the signal to noise ratio in amperometric bionano sensors.     

Yeast cytochrome c peroxidase: mechanistic studies via protein engineering

Cytochrome c peroxidase (CcP) is a yeast mitochondrial enzyme that catalyzes the reduction of hydrogen peroxide to water by ferrocytochrome c. It was the first heme enzyme to have its crystallographic structure determined and, as a consequence, has played a pivotal role in developing ideas about structural control of heme protein reactivity. Genetic engineering of the active site of CcP, along with structural, spectroscopic, and kinetic characterization of the mutant proteins has provided considerable insight into the mechanism of hydrogen peroxide activation, oxygen-oxygen bond cleavage, and formation of the higher-oxidation state intermediates in heme enzymes. The catalytic mechanism involves complex formation between cytochrome c and CcP. The cytochrome c/CcP system has been very useful in elucidating the complexities of long-range electron transfer in biological systems, including protein-protein recognition, complex formation, and intracomplex electron transfer processes.     

A needle in a haystack: the active site of the membrane-bound complex cytochrome c nitrite reductase

Cytochrome c nitrite reductase is a multicenter enzyme that uses a five-coordinated heme to perform the six-electron reduction of nitrite to ammonium. In the sulfate reducing bacterium Desulfovibrio desulfuricans ATCC 27774, the enzyme is purified as a NrfA2NrfH complex that houses 14 hemes. The number of closely-spaced hemes in this enzyme and the magnetic interactions between them make it very difficult to study the active site by using traditional spectroscopic approaches such as EPR or UV-Vis. Here, we use both catalytic and non-catalytic protein film voltammetry to simply and unambiguously determine the reduction potential of the catalytic heme over a wide range of pH and we demonstrate that proton transfer is coupled to electron transfer at the active site.     

Cytochrome bc1 complexes of microorganisms.

The cytochrome bc1 complex is the most widely occurring electron transfer complex capable of energy transduction. Cytochrome bc1 complexes are found in the plasma membranes of phylogenetically diverse photosynthetic and respiring bacteria, and in the inner mitochondrial membrane of all eucaryotic cells. In all of these species the bc1 complex transfers electrons from a low-potential quinol to a higher-potential c-type cytochrome and links this electron transfer to proton translocation. Most bacteria also possess alternative pathways of quinol oxidation capable of circumventing the bc1 complex, but these pathways generally lack the energy-transducing, protontranslocating activity of the bc1 complex. All cytochrome bc1 complexes contain three electron transfer proteins which contain four redox prosthetic groups. These are cytochrome b, which contains two b heme groups that differ in their optical and thermodynamic properties; cytochrome c1, which contains a covalently bound c-type heme; and a 2Fe-2S iron-sulfur protein. The mechanism which links proton translocation to electron transfer through these proteins is the proton motive Q cycle, and this mechanism appears to be universal to all bc1 complexes. Experimentation is currently focused on understanding selected structure-function relationships prerequisite for these redox proteins to participate in the Q-cycle mechanism. The cytochrome bc1 complexes of mitochondria differ from those of bacteria, in that the former contain six to eight supernumerary polypeptides, in addition to the three redox proteins common to bacteria and mitochondria. These extra polypeptides are encoded in the nucleus and do not contain redox prosthetic groups. The functions of the supernumerary polypeptides of the mitochondrial bc1 complexes are generally not known and are being actively explored by genetically manipulating these proteins in Saccharomyces cerevisiae.     

Structure of cytochrome c6A, a novel dithio-cytochrome of Arabidopsis thaliana, and its reactivity with plastocyanin: implications for function

Cytochrome c6A is a unique dithio-cytochrome present in land plants and some green algae. Its sequence and occurrence in the thylakoid lumen suggest that it is derived from cytochrome c6, which functions in photosynthetic electron transfer between the cytochrome b6f complex and photosystem I. Its known properties, however, and a strong indication that the disulfide group is not purely structural, indicate that it has a different, unidentified function. To help in the elucidation of this function the crystal structure of cytochrome c6A from Arabidopsis thaliana has been determined in the two redox states of the heme group, at resolutions of 1.2 A (ferric) and 1.4 A (ferrous). These two structures were virtually identical, leading to the functionally important conclusion that the heme and disulfide groups do not communicate by conformational change. They also show, however, that electron transfer between the reduced disulfide and the heme is feasible. We therefore suggest that the role of cytochrome c6A is to use its disulfide group to oxidize dithiol/disulfide groups of other proteins of the thylakoid lumen, followed by internal electron transfer from the dithiol to the heme, and re-oxidation of the heme by another thylakoid oxidant. Consistent with this model, we found a rapid electron transfer between ferro-cytochrome c6A and plastocyanin, with a second-order rate constant, k2=1.2 x 10(7) M(-1) s(-1).     

Dynamics of electron transfer pathways in cytochrome C oxidase

Cytochrome c oxidase mediates the final step of electron transfer reactions in the respiratory chain, catalyzing the transfer between cytochrome c and the molecular oxygen and concomitantly pumping protons across the inner mitochondrial membrane. We investigate the electron transfer reactions in cytochrome c oxidase, particularly the control of the effective electronic coupling by the nuclear thermal motion. The effective coupling is calculated using the Green's function technique with an extended Huckel level electronic Hamiltonian, combined with all-atom molecular dynamics of the protein in a native (membrane and solvent) environment. The effective coupling between Cu(A) and heme a is found to be dominated by the pathway that starts from His(B204). The coupling between heme a and heme a(3) is dominated by a through-space jump between the two heme rings rather than by covalent pathways. In the both steps, the effective electronic coupling is robust to the thermal nuclear vibrations, thereby providing fast and efficient electron transfer.     

Interaction of an anticancer ruthenium complex HInd[RuInd2Cl4] with cytochrome c

Cytochrome c is an important electron transfer protein in the respiratory chain, shuttling electrons from cytochrome c reductase to cytochrome c oxidase. Extensive chemical modification studies indicate significant electrostatic interactions between these proteins and show that all structural and conformational changes of cytochrome c can influence the electron transport. In the present work we examine the effect of an anticancer ruthenium complex, trans-Indazolium (bisindazole) tetrachlororuthenate(III) (HInd[RuInd(2)Cl(4)]), on the conformation of cytochrome c, the state of the heme moiety, formation of the protein dimer and on the folding state of apocytochrome c. For this purpose, gel-filtration chromatography, absorption second derivative spectroscopy, circular dichroism (CD) and inductively coupled plasma atomic emission spectroscopy (ICP(AES)) were used. The present data have revealed that binding of the potential anticancer drug HInd[RuInd(2)Cl(4)] complex to cytochrome c induces a conformation of the protein with less organized secondary and tertiary structure.     

Active site structure and dynamics of cytochrome c3 from Desulfovibrio gigas immobilized on electrodes

Cytochrome c3 from Desulfovibrio gigas is electrostatically adsorbed on Ag electrodes coated with self-assembled monolayers (SAMs) of 11-mercaptoundecanoic acid. The redox equilibria and electron transfer dynamics of the adsorbed four-heme protein are studied by surface enhanced resonance Raman spectroscopy. Immobilization on the coated electrodes does not cause any structural changes in the redox sites. The potential-dependent stationary experiments distinguish the redox potential of heme IV (-0.19 V versus normal hydrogen electrode) from those of the other hemes for which an average value of -0.3 V is determined. Taking into account the interfacial potential drops, these values are in good agreement with the redox potentials of the protein in solution. The heterogenous electron transfer between the electrode and heme IV of the adsorbed cytochrome c3 is analyzed on the basis of time-resolved experiments, leading to a formal electron transfer rate constant of 15 s(-1), which is a factor of 3 smaller than that of the monoheme protein cytochrome c.     

Catalysis and electron transfer in protein crystals: the binary and ternary complexes of methylamine dehydrogenase with electron acceptors

Polarized absorption microspectrophotometry has been used to detect catalysis and intermolecular electron transfer in single crystals of two multiprotein complexes: (1) the binary complex between Paracoccus denitrificans methylamine dehydrogenase, which contains tryptophan-tryptophylquinone (TTQ) as a cofactor, and its redox partner, the blue copper protein amicyanin; (2) the ternary complex between the same two proteins and cytochrome c-551i. Continuous wave electron paramagnetic resonance has been used to compare the state of copper in polycrystalline powders of the two systems. While catalysis and intermolecular electron transfer from reduced TTQ to copper are too fast to be accessible to our measurements, heme reduction occurs over a period of several minutes. The observed rate constant is about four orders of magnitude lower than in solution. The analysis of the temperature dependence of this apparent constant provides values for the parameters H(AB), related to electronic coupling between the two centers, and lambda, the reorganizational energy, that are compatible with electron transfer being the rate-determining step. From these parameters and the known distance between copper and heme, it is possible to calculate the parameter beta, which depends on the nature of the intervening medium, obtaining a value typical of electron transfer across a protein matrix. These findings suggest that the ternary complex in solution might achieve a higher efficiency than the rigid crystal structure thanks to an as yet unidentified role of protein dynamics.     

Photosynthetic cytochrome c550

Cytochrome c550 (cyt c550) is a membrane component of the PSII complex in cyanobacteria and some eukaryotic algae, such as red and brown algae. Cyt c550 presents a bis-histidine heme coordination which is very unusual for monoheme c-type cytochromes. In PSII, the cyt c550 with the other extrinsic proteins stabilizes the binding of Cl(-) and Ca(2+) ions to the oxygen evolving complex and protects the Mn(4)Ca cluster from attack by bulk reductants. The role (if there is one) of the heme of the cyt c550 is unknown. The low midpoint redox potential (E(m)) of the purified soluble form (from -250 to -314mV) is incompatible with a redox function in PSII. However, more positive values for the Em have been obtained for the cyt c550 bound to the PSII. A very recent work has shown an E(m) value of +200mV. These data open the possibility of a redox function for this protein in electron transfer in PSII. Despite the long distance (22?) between cyt c550 and the nearest redox cofactor (Mn(4)Ca cluster), an electron transfer reaction between these components is possible. Some kind of protective cycle involving a soluble redox component in the lumen has also been proposed. The aim of this article is to review previous studies done on cyt c550 and to consider its function in the light of the new results obtained in recent years. The emphasis is on the physical properties of the heme and its redox properties. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.     

Cytochrome P-450scc-adrenodoxin interactions. Ionic effects on binding, and regulation of cytochrome reduction by bound steroid substrates

The binding of adrenodoxin to cytochrome P-450scc and the intracomplex electron transfer from the iron-sulfur center to the heme have been studied. Salt sensitivity of the protein complex suggests the participation of electrostatic forces, as is also seen for the complex of adrenodoxin with NADPH-adrenodoxin reductase. Differences in ion specificities for the complexes of adrenodoxin with the other two proteins suggest some differences in binding requirements. Insensitivity of the heme reduction to solution conditions (salt, detergent) and kinetic analysis indicate that the protein complex is formed rapidly and that intracomplex electron transfer then occurs more slowly. Factors governing the rate of this electron transfer were investigated; binding of a series of cholesterol derivatives was used to perturb the spin state, midpoint potential, and reduction rate of the heme, and thus to test for relationships among these parameters. A linear free energy relationship between the substrate-induced midpoint potential and reduction rate is seen, but none of the other parameters (including the strength of substrate binding) are correlated. Data indicate that factors other than spin state (i.e. steric requirements and bonding groups within the steroid-binding site) regulate the strength of steroid binding. The bound steroid then modulates both midpoint potential/reduction rate and spin state but by independent mechanisms.     

Rubredoxin as a paramagnetic relaxation-inducing probe

The paramagnetic effect due to the presence of a metal center with unpaired electrons is no longer considered a hindrance in protein NMR spectroscopy. In the present work, the paramagnetic effect due to the presence of a metal center with unpaired electrons was used to map the interface of an electron transfer complex. Desulfovibrio gigas cytochrome c₃ was chosen as target to study the effect of the paramagnetic probe, Fe-rubredoxin, which produced specific line broadening in the heme IV methyl resonances M2¹ and M18¹. The rubredoxin binding surface in the complex with cytochrome c₃ was identified in a heteronuclear 2D NMR titration. The identified heme methyls on cytochrome c₃ are involved in the binding interface of the complex, a result that is in agreement with the predicted complexes obtained by restrained molecular docking, which shows a cluster of possible solutions near heme IV. The use of a paramagnetic probe in ¹HNMR titration and the mapping of the complex interface, in combination with a molecular simulation algorithm proved to be a valuable strategy to study electron transfer complexes involving non-heme iron proteins and cytochromes.     

Redox interaction of cytochrome c3 with [NiFe] hydrogenase from Desulfovibrio vulgaris Miyazaki F

Cytochrome c3 isolated from a sulfate-reducing bacterium, Desulfovibrio vulgaris Miyazaki F, is a tetraheme protein. Its physiological partner, [NiFe] hydrogenase, catalyzes the reversible oxidoreduction of molecular hydrogen. To elucidate the mechanism of electron transfer between cytochrome c3 and [NiFe] hydrogenase, the transient complex formation by these proteins was investigated by means of NMR. All NH signals of uniformly 15N-labeled ferric cytochrome c3 except N-terminus, Pro, and Gly73 were assigned. 1H-15N HSQC spectra were recorded for 15N-labeled ferric and ferrous cytochrome c3, in the absence and presence of hydrogenase. Chemical shift perturbations were observed in the region around heme 4 in both oxidation states. Additionally, the region between hemes 1 and 3 in ferrous cytochrome c3 was affected in the presence of hydrogenase, suggesting that the mode of interaction is different in each redox state. Heme 3 is probably the electron gate for ferrous cytochrome c3. To investigate the transient complex of cytochrome c3 and hydrogenase in detail, modeling of the complex was performed for the oxidized proteins using a docking program, ZDOCK 2.3, and NMR data. Furthermore, the roles of lysine residues of cytochrome c3 in the interaction with hydrogenase were investigated by site-directed mutagenesis. When the lysine residues around heme 4 were replaced by an uncharged residue, methionine, one by one, the Km of the electron-transfer kinetics increased. The results showed that the positive charges of Lys60, Lys72, Lys95, and Lys101 around heme 4 are important for formation of the transient complex with [NiFe] hydrogenase in the initial stage of the cytochrome c3 reduction. This finding is consistent with the most possible structure of the transient complex obtained by modeling.     

Conformational dynamics and molecular interaction reactions of recombinant cytochrome p450scc (CYP11A1) detected by fluorescence energy transfer.

Bovine adrenocortical cytochrome P450scc (P450scc) was expressed in Escherichia coli and purified as the substrate bound, high-spin complex (16.7 nmol of heme per mg of protein, expression level in E. coli about 400-700 nmol/l). The recombinant protein was characterized by comparison with native P450scc purified from adrenal cortex mitochondria. To study the interaction of the electron transfer proteins during the functioning of the heme protein, recombinant P450scc was selectively modified with fluorescein isothiocyanate (FITC). The present paper shows that modified P450scc, purified by affinity chromatography using adrenodoxin-Sepharose to remove non-covalently bound FITC, retains the functional activity of the unmodified enzyme, including its ability to bind adrenodoxin. Based on the efficiency of resonance fluorescence energy transfer in the donor-acceptor pair, FITC-heme, we calculated the distance between Lys(338), selectively labeled with the dye, and the heme of P450scc. The intensity of fluorescence from the label dramatically changes during: (a) denaturation of P450scc; (b) changing the spin state or redox potential of the heme protein; (c) formation of the carbon monoxide complex of reduced P450scc; (d) as well as during reactions of intermolecular interactions, such as changes of the state of aggregation, complex formation with the substrate, binding to the electron transfer partner adrenodoxin, or insertion of the protein into an artificial phospholipid membrane. Selective chemical modification of P450scc with FITC proved to be a very useful method to study the dynamics of conformational changes of the recombinant heme protein. The data obtained indicate that functionally important conformational changes of P450scc are large-scale ones, i.e. they are not limited only to changes in the dynamics of the protein active center. The results of the present study also indicate that chemical modification of Lys(338) of bovine adrenocortical P450scc does not dramatically alter the activity of the heme protein, but does result in a decrease of protein stability.     

Cytochrome P450 biosensors-a review

Cytochrome P450 (CYP) is a large family of enzymes containing heme as the active site. Since their discovery and the elucidation of their structure, they have attracted the interest of scientist for many years, particularly due to their catalytic abilities. Since the late 1970s attempts have concentrated on the construction and development of electrochemical sensors. Although sensors based on mediated electron transfer have also been constructed, the direct electron transfer approach has attracted most of the interest. This has enabled the investigation of the electrochemical properties of the various isoforms of CYP. Furthermore, CYP utilized to construct biosensors for the determination of substrates important in environmental monitoring, pharmaceutical industry and clinical practice.     

Contribution of leucine 85 to the structure and function of Saccharomyces cerevisiae iso-1 cytochrome c.

Cytochrome c is a small electron-transport protein whose major role is to transfer electrons between complex III (cytochrome reductase) and complex IV (cytochrome c oxidase) in the inner mitochondrial membrane of eukaryotes. Cytochrome c is used as a model for the examination of protein folding and structure and for the study of biological electron-transport processes. Amongst 96 cytochrome c sequences, residue 85 is generally conserved as either isoleucine or leucine. Spatially, the side chain is associated closely with that of the invariant residue Phe82, and this interaction may be important for optimal cytochrome c activity. The functional role of residue 85 has been examined using six site-directed mutants of Saccharomyces cerevisiae iso-1 cytochrome c, including, for the first time, kinetic data for electron transfer with the principle physiological partners. Results indicate two likely roles for the residue: first, heme crevice resistance to ligand exchange, sensitive to both the hydrophobicity and volume of the side chain; second, modulation of electron-transport activity through maintenance of the hydrophobic character of the protein in the vicinity of Phe82 and the exposed heme edge, and possibly of the ability of this region to facilitate redox-linked conformational change.     

Definition of the interaction domain for cytochrome c on cytochrome c oxidase. III. Prediction of the docked complex by a complete, systematic search

The electron transfer complex between bovine cytochrome c oxidase and horse cytochrome c has been predicted with the docking program DOT, which performs a complete, systematic search over all six rotational and translational degrees of freedom. Energies for over 36 billion configurations were calculated, providing a free-energy landscape showing guidance of positively charged cytochrome c to the negative region on the cytochrome c oxidase surface formed by subunit II. In a representative configuration, the solvent-exposed cytochrome c heme edge is within 4 A of the indole ring of subunit II residue Trp(104), indicating a likely electron transfer path. These two groups are surrounded by a small, hydrophobic contact region, which is surrounded by electrostatically complementary hydrophilic interactions. Cytochrome c/cytochrome c oxidase interactions of Lys(13) with Asp(119) and Lys(72) with Gln(103) and Asp(158) are the most critical polar interactions due to their proximity to the hydrophobic region and exclusion from bulk solvent. The predicted complex matches previous mutagenesis, binding, and time-resolved kinetics studies that implicate Trp(104) in electron transfer and show the importance of specific charged residues to protein affinity. Electrostatic forces not only enhance long range protein/protein association; they also predominate in short range alignment, creating the transient interaction needed for rapid turnover.     

Electron transfer between globular proteins. Dependence of the rate of transfer on distance

Based on the assumption that electron transfer between globular proteins occurs by a collective excitation of polaron type, the dependence of the rate of this process on the distance between the donor and acceptor centers with regard to their detailed electron structure was calculated. The electron structure of the heme was calculated by the quantum-chemical MNDO-PM3 method. The results were compared with experimental data on interprotein and intraglobular electron transfer. It is shown that, in the framework of this model, the electron transfer is not exponential and does not require a particular transfer pathway since the whole protein macromolecule is involved in the formation of the electron excited state.     

Photo-induced cyclic electron transfer involving cytochrome bc1 complex and reaction center in the obligate aerobic phototroph Roseobacter denitrificans.

Flash-induced redox changes of b-type and c-type cytochromes have been studied in chromatophores from the aerobic photosynthetic bacterium Roseobacter denitrificans under redox-controlled conditions. The flash-oxidized primary donor P+ of the reaction center (RC) is rapidly re-reduced by heme H1 (Em,7 = 290 mV), heme H2 (Em,7 = 240 mV) or low-potential hemes L1/L2 (Em,7 = 90 mV) of the RC-bound tetraheme, depending on their redox state before photoexcitation. By titrating the extent of flash-induced low-potential heme oxidation, a midpoint potential equal to -50 mV has been determined for the primary quinone acceptor QA. Only the photo-oxidized heme H2 is re-reduced in tens of milliseconds, in a reaction sensitive to inhibitors of the bc1 complex, leading to the concomitant oxidation of a cytochrome c spectrally distinct from the RC-bound hemes. This reaction involves cytochrome c551 in a diffusional process. Participation of the bc1 complex in a cyclic electron transfer chain has been demonstrated by detection of flash-induced reduction of cytochrome b561, stimulated by antimycin and inhibited by myxothiazol. Cytochrome b561, reduced upon flash excitation, is re-oxidized slowly even in the absence of antimycin. The rate of reduction of cytochrome b561 in the presence of antimycin increases upon lowering the ambient redox potential, most likely reflecting the progressive prereduction of the ubiquinone pool. Chromatophores contain approximately 20 ubiquinone-10 molecules per RC. At the optimal redox poise, approximately 0.3 cytochrome b molecules per RC are reduced following flash excitation. Cytochrome b reduction titrates out at Eh < 100 mV, when low-potential heme(s) rapidly re-reduce P+ preventing cyclic electron transfer. Results can be rationalized in the framework of a Q-cycle-type model.     

Spectral components of the α-band of cytochrome oxidase

Oxidative redox titrations of the mitochondrial cytochromes were performed in near-anoxic RAW 264.7 cells by inhibiting complex I. Cytochrome oxidation changes were measured with multi-wavelength spectroscopy and the ambient redox potential was calculated from the oxidation state of endogenous cytochrome c. Two spectral components were separated in the α-band range of cytochrome oxidase and they were identified as the difference spectrum of heme a when it has a high (a(H)) or low (a(L)) midpoint potential (E(m)) by comparing their occupancy during redox titrations carried out when the membrane potential (ΔΨ) was dissipated with a protonophore to that predicted by the neoclassical model of redox cooperativity. The difference spectrum of a(L) has a maximum at 605nm whereas the spectrum of a(H) has a maximum at 602nm. The ΔΨ-dependent shift in the E(m) of a(H) was too great to be accounted for by electron transfer from cytochrome c to heme a against ΔΨ but was consistent with a model in which a(H) is formed after proton uptake against ΔΨ suggesting that the spectral changes are the result of protonation. A stochastic simulation was implemented to model oxidation states, proton uptake and E(m) changes during redox titrations. The redox anti-cooperativity between heme a and heme a(3), and proton binding, could be simulated with a model where the pump proton interacted with heme a and the substrate proton interacted with heme a(3) with anti-cooperativity between proton binding sites, but not with a single proton binding site coupled to both hemes.