Effects of protein-protein interactions on electron transfer: docking and electron transfer calculations for complexes between flavodoxin and c-type cytochromes

Human Cytochrome P450 3A4 as a Biocatalyst: Effects of the Engineered Linker in Modulation of Coupling Efficiency in 3A4-BMR Chimeras

Human liver cytochrome P450 3A4 is the main enzyme involved in drug metabolism. This makes it an attractive target for biocatalytic applications, such as the synthesis of pharmaceuticals and drug metabolites. However, its poor solubility, stability and low coupling have limited its application in the biotechnological context. We previously demonstrated that the solubility of P450 3A4 can be increased by creating fusion proteins between the reductase from Bacillus megaterium BM3 (BMR) and the N-terminally modified P450 3A4 (3A4-BMR). In this work, we aim at increasing stability and coupling efficiency by varying the length of the loop connecting the two domains to allow higher inter-domain flexibility, optimizing the interaction between the domains. Starting from the construct 3A4-BMR containing the short linker Pro-Ser-Arg, two constructs were generated by introducing a 3 and 5 glycine hinge (3A4-3GLY-BMR and 3A4-5GLY-BMR). The three fusion proteins show the typical absorbance at 450 nm of the reduced heme-CO adduct as well as the correct incorporation of the FAD and FMN cofactors. Each of the three chimeric proteins were more stable than P450 3A4 alone. Moreover, the 3A4-BMR-3-GLY enzyme showed the highest NADPH oxidation rate in line with the most positive reduction potential. On the other hand, the 3A4-BMR-5-GLY fusion protein showed a Vmax increased by 2-fold as well as a higher coupling efficiency when compared to 3A4-BMR in the hydroxylation of the marker substrate testosterone. This protein also showed the highest rate value of cytochrome c reduction when this external electron acceptor is used to intercept electrons from BMR to P450. The data suggest that the flexibility and the interaction between domains in the chimeric proteins is a key parameter to improve turnover and coupling efficiency. These findings provide important guidelines in engineering catalytically self-sufficient human P450 for applications in biocatalysis.

Predicting Protein-Protein Interactions Using BiGGER: Case Studies

The importance of understanding interactomes makes preeminent the study of protein interactions and protein complexes. Traditionally, protein interactions have been elucidated by experimental methods or, with lower impact, by simulation with protein docking algorithms. This article describes features and applications of the BiGGER docking algorithm, which stands at the interface of these two approaches. BiGGER is a user-friendly docking algorithm that was specifically designed to incorporate experimental data at different stages of the simulation, to either guide the search for correct structures or help evaluate the results, in order to combine the reliability of hard data with the convenience of simulations. Herein, the applications of BiGGER are described by illustrative applications divided in three Case Studies: (Case Study A) in which no specific contact data is available; (Case Study B) when different experimental data (e.g., site-directed mutagenesis, properties of the complex, NMR chemical shift perturbation mapping, electron tunneling) on one of the partners is available; and (Case Study C) when experimental data are available for both interacting surfaces, which are used during the search and/or evaluation stage of the docking. This algorithm has been extensively used, evidencing its usefulness in a wide range of different biological research fields.

Effects of Film Morphology and Surface Chemistry on the Direct Electrochemistry of Cytochrome c at Boron-Doped Diamond Electrodes

The effects of film morphology and surface termination on the direct electron transfer of horse heart cytochrome c on boron-doped ultrananocrystalline (B-UNCD) and microcrystalline (B-MCD) diamond thin-film electrodes were investigated. Quasi-reversible, diffusion-controlled cyclic voltammetric responses were observed on oxygen-terminated (atomic O/C ~0.015), but not hydrogen-terminated (atomic O/C ~0.02) diamond thin films. The effect of the surface termination was the same for both the nanostructured B-UNCD film with sp2-bonded carbon atoms in the grain boundaries and the well faceted B-MCD film with micron-sized grains and largely devoid of sp2 carbon. Stable cyclic voltammetric i-E curves were recorded with cycling for both oxygen-terminated films indicating the absence of protein denaturation and electrode fouling. The peak currents increased linearly with the square root of the scan rate and the protein concentration; both indicative of a reaction rate limited by semi-infinite linear diffusion of the protein. Similar heterogeneous electron-transfer rate constants were observed for oxygen-terminated B-UNCD (3.48 (± 1.25) × 10-3 cm/s) and B-MCD films (2.38 (± 0.72) × 10-3 cm/s). The results clearly reveal that the oxygen-terminated surface is more active for electron-transfer with this soluble redox protein than is the hydrogen-terminated surface. The film morphology does not influence the diffusion-controlled response of the redox protein.

The effect of ionic strength on the electron-transfer rate of surface immobilized cytochrome C

Horse heart cytochrome c was immobilized on four different self-assembled monolayer (SAM) films. The electron tunneling kinetics were studied in the different assemblies as a function of the ionic strength of the buffer solution using cyclic voltammetry. When cytochrome c is electrostatically immobilized, the standard electron exchange rate constant k0 decreases with the increase of the solution's ionic strength. In contrast, the protein covalently attached or ligated has a rate constant independent of the ionic strength. The inhomogeneity of electrostatically immobilized cytochrome c increases with the increase of the solution's ionic strength whereas that of the covalently attached protein is independent of the ionic strength. A comparison of these different electron-transfer behaviors suggests that the thermodynamically stable geometry of cytochrome c in the electrostatic assemblies is also an electron transfer favorable one. It suggests that the surface charges of cytochrome c are capable of guiding it into geometries in which its front surface faces the electron-transfer partner. The inhomogeneity observed in this study indicates that a distribution of cytochrome c orientations and thus a distribution of electron transfer rate constants exists.

Protein electron transfer (mechanism and reproductive toxicity): iminium, hydrogen bonding, homoconjugation, amino acid side chains (redox and charged), and cell signaling

This contribution presents novel biochemical perspectives of protein electron transfer (ET) with focus on the iminium nature of the peptide link, along with relationships to reproductive toxicity. The favorable influence of hydrogen bonding on protein ET has been widely documented. Hydrogen bonding of the zwitterionic peptide enhances iminium character. A wide array of such bonding agents is available in vivo, with many reports on the peptide link itself. ET proceeds along the backbone, due in part, to homoconjugation. Redox amino acids (AAs), mainly tyrosine (Tyr), tryptophan (Typ), histidine (His), cysteine (Cys), disulfide, and methionine (Met), are involved in the competing processes for radical formation: direct hydrogen atom abstraction versus electron and proton loss. It appears that the radical or radical cation generated during the redox process is capable of interacting with n-electrons of the backbone. Beneficial effects of cationic AAs impact the conduction process. A relationship apparently exists involving cell signaling, protein conduction, and radicals or electrons. In addition, the link between protein ET and reproductive toxicity is examined. A key element is the role of reactive oxygen species (ROS) generated by protein ET. There is extensive evidence for involvement of ROS in generation of birth defects. The radical species arise in protein mainly by ET transformations by enzymes, as illustrated in the case of alcoholism.

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.

Sulphate respiration from hydrogen in Desulfovibrio bacteria: a structural biology overview

Sulphate-reducing organisms are widespread in anaerobic enviroments, including the gastrointestinal tract of man and other animals. The study of these bacteria has attracted much attention over the years, due also to the fact that they can have important implications in industry (in biocorrosion and souring of oil and gas deposits), health (in inflamatory bowel diseases) and the environment (bioremediation). The characterization of the various components of the electron transport chain associated with the hydrogen metabolism in Desulfovibrio has generated a large and comprehensive list of studies. This review summarizes the more relevant aspects of the current information available on the structural data of various molecules associated with hydrogen metabolism, namely hydrogenases and cytochromes. The transmembrane redox complexes known to date are also described and discussed. Redox-Bohr and cooperativity effects, observed in a few cytochromes, and believed to be important for their functional role, are discussed. Kinetic studies performed with these redox proteins, showing clues to their functional inter-relationship, are also addressed. These provide the groundwork for the application of a variety of molecular modelling approaches to understanding electron transfer and protein interactions among redox partners, leading to the characterization of several transient periplasmic complexes. In contrast to the detailed understanding of the periplasmic hydrogen oxidation process, very little is known about the cytoplasmic side of the respiratory electron transfer chain, in terms of molecular components (with exception of the terminal reductases), their structure and the protein-protein interactions involved in sulphate reduction. Therefore, a thorough understanding of the sulphate respiratory chain in Desulfovibrio remains a challenging task.

Electric-field-induced redox potential shifts of tetraheme cytochromes c3 immobilized on self-assembled monolayers: surface-enhanced resonance Raman spectroscopy and simulation studies

The tetraheme protein cytochrome c(3) (Cyt-c(3)) from Desulfovibrio gigas, immobilized on a self-assembled monolayer (SAM) of 11-mercaptoundecanoic acid, is studied by theoretical and spectroscopic methods. Molecular dynamics simulations indicate that the protein docks to the negatively charged SAM via its lysine-rich domain around the exposed heme IV. Complex formation is associated with only little protein structural perturbations. This finding is in line with the resonance Raman and surface-enhanced resonance Raman (SERR) spectroscopic results that indicate essentially the same heme pocket structures for the protein in solution and adsorbed on SAM-coated Ag electrodes. Electron- and proton-binding equilibrium calculations reveal substantial negative shifts of the redox potentials compared to the protein in solution. The magnitude of these shifts decreases in the order heme IV (-161 mV) > heme III (-73 mV) > heme II (-57 mV) > heme I (-26 mV), resulting in a change of the order of reduction. These shifts originate from the distance-dependent electrostatic interactions between the SAM headgroups and the individual hemes, leading to a stabilization of the oxidized forms. The results of the potential-dependent SERR spectroscopic analyses are consistent with the theoretical predictions and afford redox potential shifts of -160 mV (heme IV), -90 mV (heme III), -70 mV (heme II), and +20 mV (heme I) relative to the experimental redox potentials for Cyt-c(3) in solution. SERR spectroscopic experiments reveal electric-field-induced changes of the redox potentials also for the structurally very similar Cyt-c(3) from Desulfovibrio vulgaris, although the shifts are somewhat smaller compared to Cyt-c(3) from D. gigas. This study suggests that electric-field-induced redox potential shifts may also occur upon binding to biomembranes or partner proteins and thus may affect biological electron transfer processes.

Modeling electron transfer thermodynamics in protein complexes: interaction between two cytochromes c(3)

Redox protein complexes between type I and type II tetraheme cytochromes c(3) from Desulfovibrio vulgaris Hildenborough are here analyzed using theoretical methodologies. Various complexes were generated using rigid-body docking techniques, and the two lowest energy complexes (1 and 2) were relaxed using molecular dynamics simulations with explicit solvent and subjected to further characterization. Complex 1 corresponds to an interaction between hemes I from both cytochromes c(3). Complex 2 corresponds to an interaction between the heme IV from type I and the heme I from type II cytochrome c(3). Binding free energy calculations using molecular mechanics, Poisson-Boltzmann, and surface accessibility methods show that complex 2 is more stable than complex 1. Thermodynamic calculations on complex 2 show that complex formation induces changes in the reduction potential of both cytochromes c(3), but the changes are larger in the type I cytochrome c(3) (the largest one occurring on heme IV, of approximately 80 mV). These changes are sufficient to invert the global titration curves of both cytochromes, generating directionally in electron transfer from type I to type II cytochrome c(3), a phenomenon of obvious thermodynamic origin and consequences, but also with kinetic implications. The existence of processes like this occurring at complex formation may constitute a natural design of efficient redox chains.

Docking and electron transfer studies between rubredoxin and rubredoxin:oxygen oxidoreductase

The interaction and electron transfer (ET) between rubredoxin (Rd) and rubredoxin:oxygen oxidoreductase (ROO) from Desulfovibrio gigas is studied by molecular modelling techniques. Experimental kinetic assays using recombinant proteins show that the Rd reoxidation by ROO displays a bell-shaped dependence on ionic strength, suggesting a non-trivial electrostatic dependence of the interaction between these two proteins. Rigid docking studies reveal a prevalence for Rd to interact, in a very specific way, with the surface of the ROO dimer near its FMN cofactors. The optimization of the lowest energy complexes, using molecular dynamics simulation, shows a very tight interaction between the surface of the two proteins, with a high probability for Rd residues (but not the iron centre directly) to be in direct contact with the FMN cofactors of ROO. Both electrostatics and van der Waals interactions contribute to the final energy of the complex. In these complexes, the major contributions for complex formation are polar interactions between acidic residues of Rd and basic residues of ROO, plus substantial non-polar interactions between different groups. Important residues for this process are identified. ET estimates (using the Pathways model), in the optimized lowest energy complexes, suggest that these configurations are efficient for transferring electrons. The experimental bell-shaped dependence of kinetics on ionic strength is analysed in view of the molecular modelling results, and hypotheses for the molecular basis of this phenomenon are discussed.

Protein-protein interactions: mechanisms and modification by drugs

Protein-protein interactions form the proteinaceous network, which plays a central role in numerous processes in the cell. This review highlights the main structures, properties of contact surfaces, and forces involved in protein-protein interactions. The properties of protein contact surfaces depend on their functions. The characteristics of contact surfaces of short-lived protein complexes share some similarities with the active sites of enzymes. The contact surfaces of permanent complexes resemble domain contacts or the protein core. It is reasonable to consider protein-protein complex formation as a continuation of protein folding. The contact surfaces of the protein complexes have unique structure and properties, so they represent prospective targets for a new generation of drugs. During the last decade, numerous investigations have been undertaken to find or design small molecules that block protein dimerization or protein(peptide)-receptor interaction, or on the other hand, induce protein dimerization.

Gene cluster of Rhodothermus marinus high-potential iron-sulfur Protein: oxygen oxidoreductase, a caa(3)-type oxidase belonging to the superfamily of heme-copper oxidases

The respiratory chain of the thermohalophilic bacterium Rhodothermus marinus contains an oxygen reductase, which uses HiPIP (high potential iron-sulfur protein) as an electron donor. The structural genes encoding the four subunits of this HiPIP:oxygen oxidoreductase were cloned and sequenced. The genes for subunits II, I, III, and IV (named rcoxA to rcoxD) are found in this order and seemed to be organized in an operon of at least five genes with a terminator structure a few nucleotides downstream of rcoxD. Examination of the amino acid sequence of the Rcox subunits shows that the subunits of the R. marinus enzyme have homology to the corresponding subunits of oxidases belonging to the superfamily of heme-copper oxidases. RcoxB has the conserved histidines involved in binding the binuclear center and the low-spin heme. All of the residues proposed to be involved in proton transfer channels are conserved, with the exception of the key glutamate residue of the D-channel (E(278), Paracoccus denitrificans numbering). Analysis of the homology-derived structural model of subunit I shows that the phenol group of a tyrosine (Y) residue and the hydroxyl group of the following serine (S) may functionally substitute the glutamate carboxyl in proton transfer. RcoxA has an additional sequence for heme C binding, after the Cu(A) domain, that is characteristic of caa(3) oxidases belonging to the superfamily. Homology modeling of the structure of this cytochrome domain of subunit II shows no marked electrostatic character, especially around the heme edge region, suggesting that the interaction with a redox partner is not of an electrostatic nature. This observation is analyzed in relation to the electron donor for this caa(3) oxidase, the HiPIP. In conclusion, it is shown that an oxidase, which uses an iron-sulfur protein as an electron donor, is structurally related to the caa(3) class of heme-copper cytochrome c oxidases. The data are discussed in the framework of the evolution of oxidases within the superfamily of heme-copper oxidases.