Interaction of cytochrome c with cytochrome c oxidase: an NMR study on two soluble fragments derived from Paracoccus denitrificans

H-bond refinement for electron transfer membrane-bound protein-protein complexes: cytochrome c oxidase and cytochrome c552

In this study we propose a protocol to evaluate membrane-bound cytochrome c oxidase-cytochrome c552 docking candidates. An initial rigid docking algorithm generates docking poses of the cytochrome c oxidase-cytochrome c552, candidates are then aggregated into a 512-DPPC membrane model and solvated in explicit solvent. Molecular dynamic simulations are performed to induce conformational changes to membrane-bound protein complexes. Lastly each protein-protein complex is optimized in terms of its hydrogen bond network, evaluated energetically and ranked. The protocol is directly applicable to other membrane-protein complexes, such as protein-ligand systems.

Structural insights into electron transfer in caa3-type cytochrome oxidase

Cytochrome c oxidase is a member of the heme copper oxidase superfamily (HCO)¹. HCOs function as the terminal enzymes in the respiratory chain of mitochondria and aerobic prokaryotes, coupling molecular oxygen reduction to transmembrane proton pumping. Integral to the enzyme’s function is the transfer of electrons from cytochrome c to the oxidase via a transient association of the two proteins. Electron entry and exit are proposed to occur from the same site on cytochrome c^2–4. Here we report the crystal structure of the caa₃-type cytochrome oxidase from Thermus thermophilus, which has a covalently tethered cytochrome c domain. Crystals were grown in a bicontinuous mesophase using a synthetic short-chain monoacylglycerol as the hosting lipid. From the electron density map, at 2.36 Å resolution, a novel integral membrane subunit and a native glycoglycerophospholipid embedded in the complex were identified. Contrary to previous electron transfer mechanisms observed for soluble cytochrome c, the structure reveals the architecture of the electron transfer complex for the fused cupredoxin/cytochrome c domain which implicates different sites on cytochrome c for electron entry and exit. Support for an alternative to the classical proton gate characteristic of this HCO class is presented.

Electrochemical and infrared spectroscopic analysis of the interaction of the Cu(A) domain and cytochrome c(552) from Thermus thermophilus

The hydrophobically guided complex formation between the Cu(A) fragment from Thermus thermophilus ba(3) terminal oxidase and its electron transfer substrate, cytochrome c(552), was investigated electrochemically. In the presence of the purified Cu(A) fragment, a clear downshift of the c(552) redox potential from 171 to 111mV±10mV vs SHE' was found. Interestingly, this potential change fully matches complex formation with this electron acceptor site in other oxidases guided by electrostatic or covalent interactions. Redox induced FTIR difference spectra revealed conformational changes associated with complex formation and indicated the involvement of heme propionates. This article is part of a Special Issue entitled: 17th European Bioenergetics Conference (EBEC 2012).

Advances in nuclear magnetic resonance for drug discovery

Nuclear Magnetic Resonance (NMR) techniques are widely used in the drug discovery process. The primary feature exploited in these investigations is the large difference in mass between drugs and receptors (usually proteins) and the effect this has on the rotational or translational correlation times for drugs bound to their targets. Many NMR parameters, such as the diffusion coefficient, spin diffusion, nuclear Overhauser enhancement, and transverse and longitudinal relaxation times, are strong functions of either the overall tumbling or translation of molecules in solution. This has led to the development of a wide variety of NMR techniques applicable to the elucidation of protein and nucleic acid structure in solution, the screening of drug candidates for binding to a target of choice, and the study of the conformational changes which occur in a target upon drug binding. High-throughput screening by NMR methods has recently received a boost from the introduction of sophisticated computational techniques for reducing the time needed for the acquisition of the primary NMR data for multidimensional studies.

NMR basis for interprotein electron transfer gating between cytochrome c and cytochrome c oxidase

The final interprotein electron transfer (ET) in the mammalian respiratory chain, from cytochrome c (Cyt c) to cytochrome c oxidase (CcO) is investigated by (1)H-(15)N heteronuclear single quantum coherence spectral analysis. The chemical shift perturbation in isotope-labeled Cyt c induced by addition of unlabeled CcO indicates that the hydrophobic heme periphery and adjacent hydrophobic amino acid residues of Cyt c dominantly contribute to the complex formation, whereas charged residues near the hydrophobic core refine the orientation of Cyt c to provide well controlled ET. Upon oxidation of Cyt c, the specific line broadening of N-H signals disappeared and high field (1)H chemical shifts of the N-terminal helix were observed, suggesting that the interactions of the N-terminal helix with CcO are reduced by steric constraint in oxidized Cyt c, while the chemical shift perturbations in the C-terminal helix indicate notable interactions of oxidized Cyt c with CcO. These results suggest that the overall affinity of oxidized Cyt c for CcO is significantly, but not very much weaker than that of reduced Cyt c. Thus, electron transfer is gated by dissociation of oxidized Cyt c from CcO, the rate of which is controlled by the affinity of oxidized Cyt c to CcO for providing an appropriate electron transfer rate for the most effective energy coupling. The conformational changes in Lys13 upon CcO binding to oxidized Cyt c, shown by (1)H- and (1)H, (15)N-chemical shifts, are also expected to gate intraprotein ET by a polarity control of heme c environment.

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(3) 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(1) and M18(1). The rubredoxin binding surface in the complex with cytochrome c(3) was identified in a heteronuclear 2D NMR titration. The identified heme methyls on cytochrome c(3) 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 (1)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.

Association of cytochrome c with membrane-bound cytochrome c oxidase proceeds parallel to the membrane rather than in bulk solution

Electron transfer between the water-soluble cytochrome c and the integral membrane protein cytochrome c oxidase (COX) is the terminal reaction in the respiratory chain. The first step in this reaction is the diffusional association of cytochrome c toward COX, and it is still not completely clear whether cytochrome c diffuses in the bulk solution while encountering COX, or whether it prefers to diffuse laterally on the membrane surface. This is a rather crucial question, since in the latter case the association would be strongly dependent on the lipid composition and the presence of additional membrane proteins. We applied Brownian dynamics simulations to investigate the effect of an atomistically modeled dipalmitoyl phosphatidylcholine membrane on the association behavior of cytochrome c toward COX from Paracoccus denitrificans. We studied the negatively charged, physiological electron-transfer partner of COX, cytochrome c(552), and the positively charged horse-heart cytochrome c. As expected, both cytochrome c species prefer diffusion in bulk solution while associating toward COX embedded in a membrane, where the partial charges of the lipids were switched off, and the corresponding optimal association pathways largely overlap with the association toward fully solvated COX. Remarkably, after switching on the lipid partial charges, both cytochrome c species were strongly attracted by the inhomogeneous charge distribution caused by the zwitterionic lipid headgroups. This effect is particularly enhanced for horse-heart cytochrome c and is stronger at lower ionic strength. We therefore conclude that in the presence of a polar or even a charged membrane, cytochrome c diffuses laterally rather than in three dimensions.

Probing the Paracoccus denitrificans cytochrome c(1)-cytochrome c(552) interaction by mutagenesis and fast kinetics

Electron transfer (ET) between Paracoccus denitrificans cytochrome (cyt) c(1) and cytochrome c(552) was studied using the soluble redox fragments cyt c(1CF) and cyt c(552F). A new ruthenium cyt c(552F) derivative labeled at C23 (Ru(z)-23-c(552F)) was designed to measure rapid electron transfer with cyt c(1CF) in the physiological direction using flash photolysis. The bimolecular rate constant k(12) decreased rapidly with ionic strength above 40 mM, consistent with a diffusional process guided by long-range electrostatic interactions between the two proteins. However, a new kinetic phase was detected at an ionic strength of <35 mM with the ruthenium photoexcitation technique in which k(12) became very rapid (3 x 10(9) M(-1) s(-1)) and nearly independent of ionic strength, suggesting that the reaction became so fast that it was controlled by short-range diffusion along the protein surfaces guided by hydrophobic interactions. These results are consistent with a two-step model for formation of the final encounter complex. No intracomplex electron transfer between Ru(z)-23-c(552F) and c(1CF) was observed even at the lowest ionic strength, indicating that the dissociation constant of the complex was >30 microM. On the other hand, the ruthenium-labeled yeast cytochrome c derivative Ru(z)-39-Cc formed a tight 1:1 complex with cyt c(1CF) at ionic strengths of <60 mM with an intracomplex electron transfer rate constant of 50000 s(-1). A group of cyt c(1CF) variants in the presumed docking site were generated on the basis of information from the yeast cyt bc(1)-cyt c cocrystal structure. Kinetic analysis of cyt c(1CF) mutants located near the heme crevice provided preliminary identification of the interaction site for cyt c(552F) and suggested that formation of the encounter complex is guided primarily by the overall electrostatic surface potential rather than by defined ions.

Electron transfer and energy transduction in the terminal part of the respiratory chain - lessons from bacterial model systems

This review focuses on the terminal part of the respiratory chain where, macroscopically speaking, electron transfer (ET) switches from the two-electron donor, ubiquinol, to the single-electron carrier, cytochrome c, to finally reduce the four-electron acceptor dioxygen. With 3-D structures of prominent representatives of such multi-subunit membrane complexes known for some time, this section of the ET chain still leaves a number of key questions unanswered. The two relevant enzymes, ubiquinol:cytochrome c oxidoreductase and cytochrome c oxidase, appear as rather diverse modules, differing largely in their design for substrate interaction, internal ET, and moreover, in their mechanisms of energy transduction. While the canonical mitochondrial complexes have been investigated for almost five decades, the corresponding bacterial enzymes have been established only recently as attractive model systems to address basic reactions in ET and energy transduction. Lacking the intricate coding background and mitochondrial assembly pathways, bacterial respiratory enzymes typically offer a much simpler subunit composition, while maintaining all fundamental functions established for their complex "relatives". Moreover, related issues ranging from primary steps in cofactor insertion to supramolecular architecture of ET complexes, can also be favourably addressed in prokaryotic systems to hone our views on prototypic structures and mechanisms common to all family members.

Electron transfer kinetics between soluble modules of Paracoccus denitrificans cytochrome c1 and its physiological redox partners

The transient electron transfer (ET) interactions between cytochrome c1 of the bc1-complex from Paracoccus denitrificans and its physiological redox partners cytochrome c552 and cytochrome c550 have been characterized functionally by stopped-flow spectroscopy. Two different soluble fragments of cytochrome c1 were generated and used together with a soluble cytochrome c552 module as a model system for interprotein ET reactions. Both c1 fragments lack the membrane anchor; the c1 core fragment (c1CF) consists of only the hydrophilic heme-carrying domain, whereas the c1 acidic fragment (c1AF) additionally contains the acidic domain unique to P. denitrificans. In order to determine the ionic strength dependencies of the ET rate constants, an optimized stopped-flow protocol was developed to overcome problems of spectral overlap, heme autoxidation and the prevalent non-pseudo first order conditions. Cytochrome c1 reveals fast bimolecular rate constants (10(7) to 10(8) M(-1) s(-1)) for the ET reaction with its physiological substrates c552 and c550, thus approaching the limit of a diffusion-controlled process, with 2 to 3 effective charges of opposite sign contributing to these interactions. No direct involvement of the N-terminal acidic c1-domain in electrostatically attracting its substrates could be detected. However, a slight preference for cytochrome c550 over c552 reacting with cyochrome c1 was found and attributed to the different functions of both cytochromes in the respiratory chain of P. denitrificans.

Protein−Protein Interactions Studied by EPR Relaxation Measurements: Cytochromecand CytochromecOxidase

The complex formed between cytochrome c oxidase from Paracoccus denitrificans and its electron-transfer partner cytochrome c has been studied by multi-frequency pulse electron paramagnetic resonance spectroscopy. The dipolar relaxation of a fast-relaxing paramagnetic center induced on a more slowly relaxing center can be used to measure their distance in the range of 1-4 nm. This method has been used here for the first time to study transient protein-protein complex formation, employing soluble fragments for both interacting species. We observed significantly enhanced transversal relaxation of the CuA center in cytochrome c oxidase due to the fast-relaxing iron of cytochrome c upon complex formation. The possibility to measure cytochrome c oxidase in the presence and absence of cytochrome c permitted us to separate the dipolar relaxation from other relaxation contributions. This allowed a quantitative simulation and interpretation of the relaxation data. The specific temperature dependence of the dipolar relaxation together with the high orientational selectivity achieved at high magnetic field values may provide detailed information on distance and relative orientation of the two proteins with respect to each other in the complex. Our experimental results cannot be explained by any single well-defined structure of the complex of cytochrome c oxidase with cytochrome c, but rather suggest that a broad distribution in distances and relative orientations between the two proteins exist within this complex.

Green tea catechins inhibit bacterial DNA gyrase by interaction with its ATP binding site

Catechins are the main ingredients of green tea extracts and have been shown to possess versatile biological activities, including antimicrobial. We determined that the catechins inhibit bacterial DNA gyrase by binding to the ATP binding site of the gyrase B subunit. In the group of four tested catechins, epigallocatechin gallate (EGCG) had the highest activity, followed by epicatechin gallate (ECG) and epigallocatechin (EGC). Specific binding to the N-terminal 24 kDa fragment of gyrase B was determined by fluorescence spectroscopy and confirmed using heteronuclear two-dimensional NMR spectroscopy of the EGCG-15N-labeled gyrase B fragment complex. Protein residues affected by binding to EGCG were identified through chemical shift perturbation. Molecular docking calculations suggest that the benzopyran ring of EGCG penetrates deeply into the active site while the galloyl moiety anchors it to the cleft through interactions with its hydroxyl groups, which explains the higher activity of EGCG and ECG.

Conformation of thec552:aa3electron transfer complex inParacoccus denitrificansstudied by EPR on oriented samples

The EPR spectral parameters of aa(3) oxidase and cyt c(552) from Paracoccus denitrificans were studied in purified oxidase and enriched cyt c(552). The orientation of the g-tensors of hemes a and c(552) were determined on partially ordered membranes, enriched cyt c(552) and a c(552):aa(3) subcomplex. The known correlation of g-tensor to molecular axes in histidine/methionine ligated hemes permits us to position cyt c(552) with respect to the parent membrane. Taken together with previous data on the interaction surface between aa(3) oxidase and cyt c(552), these results allow us to arrive at a single conformation for the c(552):aa(3) electron transfer complex.

An NMR-based identification of a peptide fragment from the beta-subunit of a G-protein showing specific interactions with the GBB domain of the Ste20 kinase in budding yeast

In mitogen-activated protein kinase (MAPK) cascades of budding yeast, pheromone-induced mating signal is transmitted by interactions between the beta-subunit of a G-protein (G-beta) and the G-beta binding (GBB) domain of Ste20 kinase. Previously, mutational analyses of the beta-subunit of G-protein had identified two critical mutations which abrogate binding of the GBB domain of Ste20. In this work, we have identified, by use of NMR spectroscopy, a peptide fragment from the G-beta that shows specific interactions with the isolated GBB domain of Ste20. A model structure of the Ste20/G-beta complex reveals that the interface of the hetero-complex may be sustained by parallel orientation of two potentially interacting helical segments that are further stabilized by ionic, hydrogen bond, and helix macro-dipole interactions.

The Electron Transfer Complex between Cytochrome c552 and the CuA Domain of the Thermus thermophilus ba3 Oxidase

The structural analysis of the redox complex between the soluble cytochrome c552 and the membrane-integral cytochrome ba3 oxidase of Thermus thermophilus is complicated by the transient nature of this protein-protein interaction. Using NMR-based chemical shift perturbation mapping, however, we identified the contact regions between cytochrome c552 and the CuA domain, the fully functional water-soluble fragment of subunit II of the ba3 oxidase. First we determined the complete backbone resonance assignments of both proteins for each redox state. Subsequently, two-dimensional [15N,1H]TROSY spectra recorded for each redox partner both in free and complexed state indicated those surface residues affected by complex formation between the two proteins. This chemical shift analysis performed for both redox states provided a topological description of the contact surface on each partner molecule. Remarkably, very pronounced indirect effects, which were observed on the back side of the heme cleft only in the reduced state, suggested that alterations of the electron distribution in the porphyrin ring due to formation of the protein-protein complex are apparently sensed even beyond the heme propionate groups. The contact residues of each redox partner, as derived from the chemical shift perturbation mapping, were employed for a protein-protein docking calculation that provided a structure ensemble of 10 closely related conformers representing the complex between cytochrome c552 and the CuA domain. Based on these structures, the electron transfer pathway from the heme of cytochrome c552 to the CuA center of the ba3 oxidase has been predicted.

The study of bimolecular reactions under non-pseudo-first order conditions

In this work a new equation which describes the time evolution of bimolecular reactions is derived and tested by experiment. The equation is general and the results show that second-order reactions of any simple type may be accurately described by a quotient of exponential functions. The model and reagent concentration dependent observed rate constants show a complex non-linear behaviour when experimental conditions deviate from pseudo-first order nevertheless reducing to the well-known linear dependence when pseudo-first order conditions are met.

NMR study of the electron transfer complex of plant ferredoxin and sulfite reductase: mapping the interaction sites of ferredoxin

Plant ferredoxin serves as the physiological electron donor for sulfite reductase, which catalyzes the reduction of sulfite to sulfide. Ferredoxin and sulfite reductase form an electrostatically stabilized 1:1 complex for the intermolecular electron transfer. The protein-protein interaction between these proteins from maize leaves was analyzed by nuclear magnetic resonance spectroscopy. Chemical shift perturbation and cross-saturation experiments successfully mapped the location of two major interaction sites of ferredoxin: region 1 including Glu-29, Glu-30, and Asp-34 and region 2 including Glu-92, Glu-93, and Glu-94. The importance of these two acidic patches for interaction with sulfite reductase was confirmed by site-specific mutation of acidic ferredoxin residues in regions 1 and 2, separately and in combination, by which the ability of mutant ferredoxins to transfer electrons and bind to sulfite reductase was additively lowered. Taken together, this study gives a clear illustration of the molecular interaction between ferredoxin and sulfite reductase. We also present data showing that this interaction surface of ferredoxin significantly differs from that when ferredoxin-NADP(+) reductase is the interaction partner.

A structural model for the adduct between cytochrome c and cytochrome c oxidase

An ensemble of structural models of the adduct between cytochrome c and cytochrome c oxidase from Paracoccus denitrificans has been calculated based on the experimental data from site-directed mutagenesis and NMR experiments that have accumulated over the last years of research on this system. The residues from each protein that are at the protein-protein interface have been identified by the above experimental work, and this information has been converted in a series of restraints explicitly used in calculations. It is found that a single static structural model cannot satisfy all experimental data simultaneously. Therefore, it is proposed that the adduct exists as a dynamic ensemble of different orientations in equilibrium, and may be represented by a combination or average of the various limiting conformations calculated here. The equilibrium involves both conformations that are competent for electron transfer and conformations that are not. Long-range recognition of the partners is driven by non-specific electrostatic interactions, while at shorter distances hydrophobic contacts tune the reciprocal orientation. Electron transfer from cytochrome bc (1) to cytochrome c oxidase is mediated through cytochrome c experiencing multiple encounters with both of its partners, only part of which are productive. The number of encounters, and thus the electron transfer rate, may be increased by the formation of a cytochrome bc (1)-cytochrome c oxidase supercomplex and/or (in human) by increasing the concentration of the two enzymes in the membrane space.

Transient complexes of redox proteins: structural and dynamic details from NMR studies

Redox proteins participate in many metabolic routes, in particular those related to energy conversion. Protein-protein complexes of redox proteins are characterized by a weak affinity and a short lifetime. Two-dimensional NMR spectroscopy has been applied to many redox protein complexes, providing a wealth of information about the process of complex formation, the nature of the interface and the dynamic properties of the complex. These studies have shown that some complexes are non-specific and exist as a dynamic ensemble of orientations while in other complexes the proteins assume a single orientation. The binding interface in these complexes consists of a small hydrophobic patch for specificity, surrounded by polar, uncharged residues that may enhance dissociation, and, in most complexes, a ring or patch of charged residues that enhances the association by electrostatic interactions. The entry and exit port of the electrons is located within the hydrophobic interaction site, ensuring rapid electron transfer from one redox centre to the next.

A Brownian dynamics study: the effect of a membrane environment on an electron transfer system

During the past few years, three-dimensional crystal structures of many of the important integral membrane proteins responsible for the bioenergetic processes of photosynthesis and respiration have been determined. Moreover, a few crystal structures of protein-protein complexes have become available that characterize the interaction between those membrane proteins and the electron carrier protein cytochrome c. Here, we address the association kinetics for binding of cytochrome c to cytochrome c oxidase (COX) from Paracoccus denitrificans by Brownian dynamics simulations. The effects of ionic strength and protein mutations were studied for two different cytochrome c species: the positively charged, dipolar horse heart cytochrome c and the negatively charged physiological electron transfer partner cytochrome c(552). We studied association toward "naked" COX and toward membrane-embedded COX where the membrane is represented as an uncharged DPPC bilayer modeled in atomistic detail. For the nonnatural association toward "naked" COX, the association rates are >100 times larger for horse heart cytochrome c than for cytochrome c(552). Interestingly, the presence of the lipid bilayer leads to a dramatic decrease of the association rate of horse heart cytochrome c, but slightly enhances association of cytochrome c(552), leading to very similar association rates of both proteins to membrane-embedded COX. This finding from computational modeling studies may reflect the optimization of surface patches and of the total net charge on electron transfer pairs in nature.

The architecture of the binding site in redox protein complexes: Implications for fast dissociation

Interprotein electron transfer is characterized by protein interactions on the millisecond time scale. Such transient encounters are ensured by extremely high rates of complex dissociation. Computational analysis of the available crystal structures of redox protein complexes reveals features of the binding site that favor fast dissociation. In particular, the complex interface is shown to have low geometric complementarity and poor packing. These features are consistent with the necessity for fast dissociation since the absence of close packing facilitates solvation of the interface and disruption of the complex.

Different interaction modes of two cytochrome-c oxidase soluble CuA fragments with their substrates

Cytochrome-c oxidase is the terminal enzyme in the respiratory chains of mitochondria and many bacteria and catalyzes the formation of water by reduction of dioxygen. The first step in the cytochrome oxidase reaction is the bimolecular electron transfer from cytochrome c to the homobinuclear mixed-valence CuA center of subunit II. In Thermus thermophilus a soluble cytochrome c552 acts as the electron donor to ba3 cytochrome-c oxidase, an interaction believed to be mainly hydrophobic. In Paracoccus denitrificans, electrostatic interactions appear to play a major role in the electron transfer process from the membrane-spanning cytochrome c552. In the present study, soluble fragments of the CuA domains and their respective cytochrome c electron donors were analyzed by stopped-flow spectroscopy to further characterize the interaction modes. The forward and the reverse electron transfer reactions were studied as a function of ionic strength and temperature, in all cases yielding monoexponential time-dependent reaction profiles in either direction. From the apparent second-order rate constants, equilibrium constants were calculated, with values of 4.8 and of 0.19, for the T. thermophilus and P. denitrificans c552 and CuA couples, respectively. Ionic strength strongly affects the electron transfer reaction in P. denitrificans indicating that about five charges on the protein interfaces control the interaction, when analyzed according to the Brønsted equation, whereas in the T. thermophilus only 0.5 charges are involved. Overall the results indicate that the soluble CuA domains are excellent models for the initial electron transfer processes in cytochrome-c oxidases.