Artificial cytochrome b: computer modeling and evaluation of redox potentials

rBet v 1a-BanLecwt induce upregulation of IL-10 and IFN-γ gene expression in Caco-2/THP-1 co-culture and secretion of IL-10 and IFN-γ/IL-4 levels in PBMCs of birch pollen allergic donors

Novel allergen immunotherapy (AIT) approaches necessitate the use of more effective and safe therapeutics, which can be accomplished by employing novel adjuvants for improved innate immune cell activation, as well as hypoallergenic allergen forms. In this study, we investigate the immunomodulatory effects of a chimera rBet v 1a-BanLecwt (rBv1a-BLwt; Cwt) composed of the major birch pollen allergen Bet v 1a and banana lectin (BanLecwt; BLwt) and two novel chimeras, rBv1l-BLH84T (rBet v 1l-BanLecH84T; C1) and rBLH84T-Bv1l (rBanLecH84T-Bet v 1l; C2), both composed of BLH84T and hypoallergenic birch pollen allergen Bv1l in the co-culture model Caco-2/THP-1, and PBMCs from donors with birch pollen allergy. The chimeric molecules rBv1l-BLH84T (C1) and rBLH84T-Bv1l (C2) were created in silico and then produced in E. coli using recombinant DNA technology. Real-time PCR analysis of gene expression following compound treatment in the co-culture model revealed that all three chimeras have the potential to induce the anti-inflammatory cytokine IL-10 gene expression in Caco-2 cells and IFN-γ gene expression in THP-1 cells. Sandwich ELISA revealed that Cwt increased IL-10 secretion and IFN-/IL-4 levels in PBMCs from birch pollen allergic donors, whereas C1 and C2 were less effective. The findings suggest that Cwt should be analyzed further due to its potential benefit in AIT.

Accelerated molecular dynamics study to compare the thermostability of Bacillus licheniformis and Aspergillus niger α-amylase

The thermostability of enzymes plays a significant role in the starch hydrolysis process in the industry. The structural difference between thermostable Bacillus licheniformis α-amylase (BLA) and thermolabile Aspergillus niger α-amylase (ANA) is interesting to be explored. This work aimed to study the thermostability-determining factor of BLA as compared to a non-thermostable enzyme, ANA, using molecular dynamics (MD) simulation at a high temperature. A 100 ns of classical MD, which was followed by 200 ns accelerated MD was conducted to explore the conformational changes of the enzyme. It is revealed that the intramolecular interactions through salt bridge interactions and the presence of calcium ions dominates the stability effect of BLA, despite the absence of a disulfide bond in the structure. These results should be useful in designing a thermostable enzyme that can be used in industrial processes.Communicated by Ramaswamy H. Sarma.

The evolution of the human mitochondrial bc1 complex- adaptation for reduced rate of superoxide production?

The mitochondrial bc1 complex is a major source of mitochondrial superoxide. While bc1-generated superoxide plays a beneficial signaling role, excess production of superoxide lead to aging and degenerative diseases. The catalytic core of bc1 comprises three peptides -cytochrome b, Fe-S protein, and cytochrome c1. All three core peptides exhibit accelerated evolution in anthropoid primates. It has been suggested that the evolution of cytochrome b in anthropoids was driven by a pressure to reduce the production of superoxide. In humans, the bc1 core peptides exhibit anthropoid-specific substitutions that are clustered near functionally critical sites that may affect the production of superoxide. Here we compare the high-resolution structures of bovine, mouse, sheep and human bc1 to identify structural changes that are associated with human-specific substitutions. Several cytochrome b substitutions in humans alter its interactions with other subunits. Most significantly, there is a cluster of seven substitutions, in cytochrome b, the Fe-S protein, and cytochrome c1 that affect the interactions between these proteins at the tether arm of the Fe-S protein and may alter the rate of ubiquinone oxidation and the rate of superoxide production. Another cluster of substitutions near heme bH and the ubiquinone reduction site, Qi, may affect the rate of ubiquinone reduction and thus alter the rate of superoxide production. These results are compatible with the hypothesis that cytochrome b in humans (and other anthropoid primates) evolve to reduce the rate of production of superoxide thus enabling the exceptional longevity and exceptional cognitive ability of humans.

Large-scale comparison between the diffraction-component precision indexes favors Cruickshank’s Rfree function

This study aims to provide a first large-scale comparison between the various diffraction-component precision index (DPI) equations, assess the applicability of the parameter, and make recommendations on DPI computation. The DPI estimates the average accuracy of the atomic coordinates obtained by the structural refinement of protein diffraction data, with application in crystallography and cheminformatics. Although, Cruickshank and Blow proposed DPI equations based on R and Rfree in order to calculate DPI values, which remain scarcely employed in the quality assessment of the Protein Data Base (PDB) files, due to the unclear data extraction protocols (to assign variables), the complex equations, the lack of extensive applicability studies and the limited access to automated computations. In order to address these shortcomings, the entire RCSB PDB database was evaluated using Cruickshank?s and Blow?s R and Rfree DPI variations. Computations of 143070 X-ray structures indicate that Rfree-based DPI equations apply to 30 % more protein structures compared to R-based DPI equations, with Cruickshank Rfree-based DPI (CRF) exceeding the number of successful Blow?s Rfree-based DPI (BRF) computations. Although our results indicate that, in general, the resolutions < 2 ? assure consistency among the various DPIs computations (differences <0.05 ?), we recommend the use of CRF DPI because of its wider applicability.

Quantification of Local Electric Field Changes at the Active Site of Cytochrome c Oxidase by Fourier Transform Infrared Spectroelectrochemical Titrations

Cytochrome c oxidase (CcO) is a transmembrane protein complex that reduces molecular oxygen to water while translocating protons across the mitochondrial membrane. Changes in the redox states of its cofactors trigger both O₂ reduction and vectorial proton transfer, which includes a proton-loading site, yet unidentified. In this work, we exploited carbon monoxide (CO) as a vibrational Stark effect (VSE) probe at the binuclear center of CcO from Rhodobacter sphaeroides. The CO stretching frequency was monitored as a function of the electrical potential, using Fourier transform infrared (FTIR) absorption spectroelectrochemistry. We observed three different redox states (R₄CO, R₂CO, and O), determined their midpoint potential, and compared the resulting electric field to electrostatic calculations. A change in the local electric field strength of +2.9 MV/cm was derived, which was induced by the redox transition from R₄CO to R₂CO. We performed potential jump experiments to accumulate the R₂CO and R₄CO species and studied the FTIR difference spectra in the protein fingerprint region. The comparison of the experimental and computational results reveals that the key glutamic acid residue E286 is protonated in the observed states, and that its hydrogen-bonding environment is disturbed upon the redox transition of heme a₃. Our experiments also suggest propionate A of heme a₃ changing its protonation state in concert with the redox state of a second cofactor, heme a. This supports the role of propionic acid side chains as part of the proton-loading site.

Role of AP-endonuclease (Ape1) active site residues in stabilization of the reactant enzyme-DNA complex

Apurinic/apyrimidinic endonuclease 1 (Ape1) is an important metal-dependent enzyme in the base excision repair mechanism, responsible for the backbone cleavage of abasic DNA through a phosphate hydrolysis reaction. Molecular dynamics simulations of Ape1 complexed to its substrate DNA performed for models containing 1 or 2 Mg²⁺ -ions as cofactor located at different positions show a complex with 1 metal ion bound on the leaving group site of the scissile phosphate to be the most likely reaction-competent conformation. Active-site residue His309 is found to be protonated based on pKa calculations and the higher conformational stability of the Ape1-DNA substrate complex compared to scenarios with neutral His309. Simulations of the D210N mutant further support the prevalence of protonated His309 and strongly suggest Asp210 as the general base for proton acceptance by a nucleophilic water molecule.

Current advances in research of cytochrome c oxidase

The function of cytochrome c oxidase as a biomolecular nanomachine that transforms energy of redox reaction into protonmotive force across a biological membrane has been subject of intense research, debate, and controversy. The structure of the enzyme has been solved for several organisms; however details of its molecular mechanism of proton pumping still remain elusive. Particularly, the identity of the proton pumping site, the key element of the mechanism, is still open to dispute. The pumping mechanism has been for a long time one of the key unsolved issues of bioenergetics and biochemistry, but with the accelerating progress in this field many important details and principles have emerged. Current advances in cytochrome oxidase research are reviewed here, along with a brief discussion of the most complete proton pumping mechanism proposed to date, and a molecular basis for control of its efficiency.

Coupled electron and proton transfer reactions during the O→E transition in bovine cytochrome c oxidase

A combined DFT/electrostatic approach is employed to study the coupling of proton and electron transfer reactions in cytochrome c oxidase (CcO) and its proton pumping mechanism. The coupling of the chemical proton to the internal electron transfer within the binuclear center is examined for the O→E transition. The novel features of the His291 pumping model are proposed, which involve timely well-synchronized sequence of the proton-coupled electron transfer reactions. The obtained pK(a)s and E(m)s of the key ionizable and redox-active groups at the different stages of the O→E transition are consistent with available experimental data. The PT step from E242 to H291 is examined in detail for various redox states of the hemes and various conformations of E242 side-chain. Redox potential calculations of the successive steps in the reaction cycle during the O→E transition are able to explain a cascade of equilibria between the different intermediate states and electron redistribution between the metal centers during the course of the catalytic activity. All four electrometric phases are discussed in the light of the obtained results, providing a robust support for the His291 model of proton pumping in CcO.

Multivalent binding of formin-binding protein 21 (FBP21)-tandem-WW domains fosters protein recognition in the pre-spliceosome

The high abundance of repetitive but nonidentical proline-rich sequences in spliceosomal proteins raises the question of how these known interaction motifs recruit their interacting protein domains. Whereas complex formation of these adaptors with individual motifs has been studied in great detail, little is known about the binding mode of domains arranged in tandem repeats and long proline-rich sequences including multiple motifs. Here we studied the interaction of the two adjacent WW domains of spliceosomal protein FBP21 with several ligands of different lengths and composition to elucidate the hallmarks of multivalent binding for this class of recognition domains. First, we show that many of the proteins that define the cellular proteome interacting with FBP21-WW1-WW2 contain multiple proline-rich motifs. Among these is the newly identified binding partner SF3B4. Fluorescence resonance energy transfer (FRET) analysis reveals the tandem-WW domains of FBP21 to interact with splicing factor 3B4 (SF3B4) in nuclear speckles where splicing takes place. Isothermal titration calorimetry and NMR shows that the tandem arrangement of WW domains and the multivalency of the proline-rich ligands both contribute to affinity enhancement. However, ligand exchange remains fast compared with the NMR time scale. Surprisingly, a N-terminal spin label attached to a bivalent ligand induces NMR line broadening of signals corresponding to both WW domains of the FBP21-WW1-WW2 protein. This suggests that distinct orientations of the ligand contribute to a delocalized and semispecific binding mode that should facilitate search processes within the spliceosome.

Similarity of cytochrome c oxidases in different organisms

Most of biological oxygen reduction is catalyzed by the heme-copper oxygen reductases. These enzymes are redox-driven proton pumps that take part in generating the proton gradient in both prokaryotes and mitochondria that drives synthesis of ATP. The enzymes have been divided into three evolutionarily-related groups: the A-, B-, and C-families. Recent comparative studies suggest that all oxygen reductases perform the same chemistry for oxygen reduction and comprise the same essential elements of the proton pumping mechanism, such as the proton loading and kinetic gating sites, which, however, appear to be different in different families. All species of the A-family, however, demonstrate remarkable similarity of the central processing unit of the enzyme, as revealed by their recent crystal structures. Here we demonstrate that cytochrome c oxidases (CcO) of such diverse organisms as a mammal (bovine heart mitochondrial CcO), photosynthetic bacteria (Rhodobacter sphaeroides CcO), and soil bacteria (Paracoccus denitrificans CcO) are not only structurally similar, but almost identical in microscopic electrostatics and thermodynamics properties of their key amino-acids. By using pK(a) calculations of some of the key residues of the catalytic site, D- and K- proton input, and putative proton output channels of these three different enzymes, we demonstrate that the microscopic properties of key residues are almost identical, which strongly suggests the same mechanism in these species. The quantitative precision with which the microscopic physical properties of these enzymes have remained constant despite different evolutionary routes undertaken is striking.

Understanding properties of cofactors in proteins: redox potentials of synthetic cytochromes b

Haehnel et al. synthesized 399 different artificial cytochrome b (aCb) models. They consist of a template-assisted four-helix bundle with one embedded heme group. Their redox potentials were measured and cover the range from -148 to -89 mV. No crystal structures of these aCb are available. Therefore, we use the chemical composition and general structural principles to generate atomic coordinates of 31 of these aCb mutants, which are chosen to cover the whole interval of redox potentials. We start by modeling the coordinates of one aCb from scratch. Its structure remains stable after energy minimization and during molecular dynamics simulation over 2 ns. Based on this structure, coordinates of the other 30 aCb mutants are modeled. The calculated redox potentials for these 31 aCb agree within 10 mV with the experimental values in terms of root mean square deviation. Analysis of the dependence of heme redox potential on protein environment shows that the shifts in redox potentials relative to the model systems in water are due to the low-dielectric medium of the protein and the protonation states of the heme propionic acid groups, which are influenced by the surrounding amino acids. Alternatively, we perform a blind prediction of the same redox potentials using an empirical approach based on a linear scoring function and reach a similar accuracy. Both methods are useful to understand and predict heme redox potentials. Based on the modeled structure we can understand the detailed structural differences between aCb mutants that give rise to shifts in heme redox potential. On the other hand, one can explore the correlation between sequence variations and aCb redox potentials more directly and on much larger scale using the empirical prediction scheme, which--thanks to its simplicity--is much faster.

DFT/electrostatic calculations of pK(a) values in cytochrome c oxidase

Using classical electrostatic calculations, earlier we examined the dependence of the protonation state of bovine cytochrome c oxidase (CcO) on its redox state. Based on these calculations, we have proposed a model of CcO proton pumping that involves His291, one of the Cu(B) histidine ligands, which was found to respond to redox changes of the enzyme Fe(a)(3)-Cu(B) catalytic center. In this work, we employ combined density functional and continuum electrostatic calculations to evaluate the pK(a)() values of His291 and Glu242, two key residues of the model. The pK(a) values are calculated for different redox states of the enzyme, and the influence of different factors on the pK(a)'s is analyzed in detail. The calculated pK(a)() values of Glu242 are between 9.4 and 12.0, depending on the redox state of the protein, which is in excellent agreement with recent experimental measurements. Assuming the reduced state of heme a(3), His291 of the oxidized Cu(B) center possesses a pK(a)() between 2.1 and 4.0, while His291 of the reduced Cu(B) center has a pK(a) above 17. The obtained results support the proposal that the His291 ligand of the Cu(B) center in CcO is a proton pump element.

Detection of heme oxygenase activity in a library of four-helix bundle proteins: towards the de novo synthesis of functional heme proteins

Design and chemical synthesis of de novo heme proteins with enzymatic activity on cellulose membranes is described. 352 antiparallel four-helix bundle proteins with a single histidine for heme ligation were assembled from three different sets of short amphipathic helices on membrane-bound peptide templates. The templates were coupled by linkers to cellulose membranes of microplate format, which could be cleaved for control of intermediate and final products. The incorporation of heme and the heme oxygenase activity of the 352 proteins were monitored by measuring UV-visible spectra directly on the cellulose. The kinetics of the heme oxygenase reaction was studied by monitoring the decrease of the Soret band and the transient absorbance of verdoheme being an intermediate product in the formation of biliverdin. Four of the proteins covering a broad range of the enzymatic rate constants were selected and synthesized in solution for further characterization. Detailed studies by redox potentiometry, analytical ultracentrifugation, and electron paramagnetic resonance spectroscopy yielded information about the aggregation state of the proteins, the spin state and the putative coordination environment of the iron. The amount of five-coordinated high-spin iron and a positive reduction potential were found to promote the oxygenase activity of the proteins.

Combined DFT and electrostatics study of the proton pumping mechanism in cytochrome c oxidase

Cytochrome c oxidase is a redox-driven proton pump which converts atmospheric oxygen to water and couples the oxygen reduction reaction to the creation of a membrane proton gradient. The structure of the enzyme has been solved; however, the mechanism of proton pumping is still poorly understood. Recent calculations from this group indicate that one of the histidine ligands of enzyme's CuB center, His291, may play the role of the pumping element. In this paper, we report on the results of calculations that combined first principles DFT and continuum electrostatics to evaluate the energetics of the key energy generating step of the model-the transfer of the chemical proton to the binuclear center of the enzyme, where the hydroxyl group is converted to water, and the concerted expulsion of the proton from delta-nitrogen of His291 ligand of CuB center. We show that the energy generated in this step is sufficient to push a proton against an electrochemical membrane gradient of about 200 mV. We have also re-calculated the pKa of His291 for an extended model in which the whole Fe(a3)-CuB center with their ligands is treated by DFT. Two different DFT functionals (B3LYP and PBE0), and various dielectric models of the protein have been used in an attempt to estimate potential errors of the calculations. Although current methods of calculations do not allow unambiguous predictions of energetics in proteins within few pKa units, as required in this case, the present calculation provides further support for the proposed His291 model of CcO pump and makes a specific prediction that could be targeted in the experimental test.

Function of redox-active tyrosine in photosystem II

Water oxidation at photosystem II Mn-cluster is mediated by the redox-active tyrosine Y(Z). We calculated the redox potential (E(m)) of Y(Z) and its symmetrical counterpart Y(D), by solving the linearized Poisson-Boltzmann equation. The calculated E(m)(Y( )/Y(-)) were +926 mV/+694 mV for Y(Z)/Y(D) with the Mn-cluster in S2 state. Together with the asymmetric position of the Mn-cluster relative to Y(Z/D), differences in H-bond network between Y(Z) (Y(Z)/D1-His(190)/D1-Asn(298)) and Y(D) (Y(D)/D2-His(189)/D2-Arg(294)/CP47-Glu(364)) are crucial for E(m)(Y(Z/D)). When D1-His(190) is protonated, corresponding to a thermally activated state, the calculated E(m)(Y(Z)) was +1216 mV, which is as high as the E(m) for P(D1/D2). We observed deprotonation at CP43-Arg(357) upon S-state transition, which may suggest its involvement in the proton exit pathway. E(m)(Y(D)) was affected by formation of P(D2)(+) (but not P(D1)(+)) and sensitive to the protonation state of D2-Arg(180). This points to an electrostatic link between Y(D) and P(D2).

Two conformational states of Glu242 and pKas in bovine cytochrome c oxidase

Cytochrome c oxidase (CcO) is the terminal enzyme in the respiratory electron transport chain of aerobic organisms. It catalyses the reduction of atmospheric oxygen to water, and couples this reaction to proton pumping across the membrane; this process generates the electrochemical gradient that subsequently drives the synthesis of ATP. The molecular details of the mechanism by which electron transfer is coupled to proton pumping in CcO is poorly understood. Recent calculations from our group indicate that His291, a ligand of the Cu(B) center of the enzyme, may play the role of the pumping element. In this paper we describe calculations in which a DFT/continuum electrostatic method is used to explore the coupling of the conformational changes of Glu242 residue, the main proton donor of both chemical and pump protons, to its pKa, and the pKa of His291, a putative proton loading site of our pumping model. The computations are done for several redox states of metal centers, different protonation states of Glu242 and His291, and two well-defined conformations of the Glu242 side chain. Thus, in addition to equilibrium redox/protonation states of the catalytic cycle, we also examine the transient and intermediate states. Different dielectric models are employed to investigate the robustness of the results, and their viability in the light of the proposed proton pumping mechanism of CcO. The main results are in agreement with the experimental measurements and support the proposed pumping mechanism. Additionally, the present calculations indicate a possibility of gating through conformational changes of Glu242; namely, in the pumping step, we find that Glu242 needs to be reprotonated before His291 can eject a proton to the P-site of membrane. As a result, the reprotonation of Glu can control proton release from the proton loading site.

Conformational transitions and redox potential shifts of cytochrome P450 induced by immobilization

Cytochrome P450 (P450) from Pseudomonas putida was immobilized on Ag electrodes coated with self-assembled monolayers (SAMs) via electrostatic and hydrophobic interactions as well as by covalent cross-linking. The redox and conformational equilibria of the immobilized protein were studied by potential-dependent surface-enhanced resonance Raman spectroscopy. All immobilization conditions lead to the formation of the cytochrome P420 (P420) form of the enzyme. The redox potential of the electrostatically adsorbed P420 is significantly more positive than in solution and shows a steady downshift upon shortening of the length of the carboxyl-terminated SAMs, i.e., upon increasing the strength of the local electric field. Thus, two opposing effects modulate the redox potential of the adsorbed enzyme. First, the increased hydrophobicity of the heme environment brought about by immobilization on the SAM tends to upshift the redox potential by stabilizing the formally neutral ferrous form. Second, increasing electric fields tend to stabilize the positively charged ferric form, producing the opposite effect. The results provide insight into the parameters that control the structure and redox properties of heme proteins and contribute to the understanding of the apparently anomalous behavior of P450 enzymes in bioelectronic devices.

Heterologous expression and site-directed mutagenesis of an ascorbate-reducible cytochrome b561

Cytochromes b561 (Cyts b561) are ubiquitous membrane proteins catalyzing ascorbate-mediated trans-membrane electron transfer. A heterologous expression system in Saccharomyces cerevisiae was developed to study their structure-function relationship. Recombinant mouse chromaffin granule Cyt b561 (CGCytb) shows spectral characteristics, ascorbate reducibility, and redox potentials identical to that of the native bovine protein. Moreover, the reconstituted recombinant protein mediated trans-membrane electron transport with kinetic characteristics similar to that of bovine CGCytb. Site-directed mutant analysis supports the presence of two hemes coordinated by the highly conserved His pairs H52/H120 and H86/H159. Reduction of CGCytb by ascorbate showed biphasic kinetics (Kd1: 0.016 +/- 0.005 mM, Kd2: 1.24 +/- 0.19 mM). Mutation of a well-conserved Arg residue (R72) abolished high affinity CGCytb reduction by ascorbate, indicating that this residue may be critical for substrate binding. On the other hand, mutation of a Lys previously suggested to play a role in ascorbate binding (K83), did not affect the ascorbate-mediated reduction of the protein.

Computational studies on imidazole heme conformations

Density functional theory computations of heme with ionized propionic acid groups, axially coordinated with two imidazoles, were performed for different mutual orientations of the imidazole planes. Environmental influences from water or protein were considered with a continuum dielectric medium by solving the Poisson equation. In vacuum, optimized geometries yielded imidazole-heme conformations where the NH groups of imidazoles are oriented toward the heme propionic groups in agreement with data from crystal structures of heme proteins. Conformational free-energy dependencies of the mutual orientation of axially ligated imidazoles calculated in protein (epsilon=10) and water (epsilon=80) environments confirmed the vacuum results, albeit the energy difference between the preferred and the 180 degrees opposite orientations of the imidazole ligand decreased from 3.84 kcal/mol in vacuum to 2.35 and 2.40 kcal/mol in protein and water, respectively. Two main factors determine the imidazole orientation: (1) the direct intramolecular electrostatic interactions of propionic groups with the polar NH groups of imidazole and (2) the electrostatic interaction of the total dipole moment of the imidazole-heme complex with the reaction field. In vacuum, only the first type of interaction is present, while in a dielectric medium the latter effect becomes competitive at high dielectric constant, resulting in a decrease of the orientational preference. Interestingly, the orientational preference of the imidazole axially ligated to heme becomes even more pronounced, if the negatively charged propionates are neutralized by counter charges that mimic salt bridges or protonation of the propionates.

Redox equilibria in hydroxylamine oxidoreductase. Electrostatic control of electron redistribution in multielectron oxidative processes

We report results of continuum electrostatics calculations of the cofactor redox potentials, and of the titratable group pK(a) values, in hydroxylamine oxidoreductase (HAO). A picture of a sophisticated multicomponent control of electron flow in the protein emerged from the studies. First, we found that neighboring heme cofactors strongly interact electrostatically, with energies of 50-100 mV. Thus, cofactor redox potentials depend on the oxidation state of other cofactors, and cofactor redox potentials in the active (partially oxidized) enzyme differ substantially from the values obtained in electrochemical redox titration experiments. We found that, together, solvent-exposed heme 1 (having a large negative redox potential) and heme 2 (having a large positive redox potential) form a lock for electrons generated during the oxidation reaction The attachment of HAO's physiological electron transfer partner cytochrome c(554) results in a positive shift in the redox potential of heme 1, and "opens the electron gate". Electrons generated as a result of hydroxylamine oxidation travel to heme 3 and heme 8, which have redox potentials close to 0 mV versus NHE (this result is in partial disagreement with an existing experimental redox potential assignment). The closeness of hemes 3 and 8 from different enzyme subunits allows redistribution of the four electrons generated as a result of hydroxylamine oxidation, among the three enzyme subunits. For the multielectron oxidation process to be maximally efficient, the redox potentials of the electron-accepting cofactors should be roughly equal, and electrostatic interactions between extra electrons on these cofactors should be minimal. The redox potential assignments presented in the paper satisfy this general rule.

Calculation of the Redox Potential of the Protein Azurin and Some Mutants

Azurin from Pseudomonas aeruginosa is a small 128-residue, copper-containing protein. Its redox potential can be modified by mutating the protein. Free-energy calculations based on classical molecular-dynamics simulations of the protein and from mutants in aqueous solution at different pH values were used to compute relative redox potentials. The precision of the free-energy calculations with the lambda coupling-parameter approach is evaluated as function of the number and sequence of lambda values, the sampling time and initial conditions. It is found that the precision is critically dependent on the relaxation of hydrogen-bonding networks when changing the atomic-charge distribution due to a change of redox state or pH value. The errors in the free energies range from 1 to 10 k(B)T, depending on the type of process. Only qualitative estimates of the change in redox potential by protein mutation can be obtained.

Proton Exit Channels in Bovine Cytochrome c Oxidase

Cytochrome c oxidase (CcO) is the terminal transmembrane enzyme of the respiratory electron transport chain in aerobic cells. It catalyzes the reduction of oxygen to water and utilizes the free energy of the reduction reaction for proton pumping, a process which results in a membrane electrochemical proton gradient. Although the structure of the enzyme has been solved for several organisms, the molecular mechanism of proton pumping and proton exit pathways remain unknown. In our previous work, the continuum electrostatic calculations were employed to evaluate the electrostatic potential, energies, and protonation state of bovine cytochrome c oxidase for different redox states of the enzyme. A possible mechanism of oxygen reduction and proton pumping via His291 was proposed. In this paper, using electrostatic calculations, we examine the proton exit pathways in the enzyme. By monitoring the changes of the protonation states, proton affinities, and energies of electrostatic interactions between the titratable groups in different redox states of CcO, we identified the clusters of strongly interacting residues. Using these data, we detected four possible proton exit points on the periplasmic side of the membrane (Lys171B/Asp173B, His24B/Asp25B, Asp51, and Asp300). We then were able to trace the proton exit pathways and to evaluate the energy profiles along the paths. On the basis of energetic considerations and the conservation of the residues in a protein sequence, the most likely exit pathway is one via the Lys171B/Asp173B site. The obtained results are fully consistent with our His291 model of proton pumping, and provide a rationale for the absence of proton leaking in CcO between the pumping strokes.

Redox and redox-coupled processes of heme proteins and enzymes at electrochemical interfaces

Modern bioelectrochemical methods rely upon the immobilisation of redox proteins and enzymes on electrodes coated with biocompatible materials to prevent denaturation. However, even when protein denaturation is effectively avoided, heterogeneous protein electron transfer is often coupled to non-Faradaic processes like reorientation, conformational transitions or acid-base equilibria. Disentangling these processes requires methods capable of probing simultaneously the structure and reaction dynamics of the adsorbed species. Here we provide an overview of the recent developments in Raman and infrared surface-enhanced spectroelectrochemical techniques applied to the study of soluble and membrane bound redox heme proteins and enzymes. Possible biological implications of the findings are critically discussed.

A Kunitz type protease inhibitor related protein is synthesized in Drosophila prepupal salivary glands and released into the moulting fluid during pupation

From the Drosophila virilis late puff region 31C, we microcloned two neighbouring genes, Kil-1 and Kil-2, that encode putative Kunitz serine protease inhibitor like proteins. The Kil-1 gene is expressed exclusively in prepupal salivary glands. Using a size mutant of the KIL-1 protein and MALDI-TOF analysis, we demonstrate that during pupation this protein is released from the prepupal salivary glands into the pupation fluid covering the surface of the pupa. 3-D-structure predictions are consistent with the known crystal structure of the human Kunitz type protease inhibitor 2KNT. This is the first experimental proof for the extracorporal presence of a distinct Drosophila prepupal salivary gland protein. Possible functions of KIL-1 in the context of the control of proteolytic activities in the pupation fluid are discussed.

Electrostatic study of the proton pumping mechanism in bovine heart cytochrome C oxidase

Cytochrome c oxidase (CcO) is the terminal enzyme of the cell respiratory chain in mitochondria and aerobic bacteria. It catalyzes the reduction of oxygen to water and utilizes the free energy of the reduction reaction for proton pumping across the inner-mitochondrial membrane, a process that results in a membrane electrochemical proton gradient. Although the structure of the enzyme has been solved for several organisms, the molecular mechanism of proton pumping remains unknown. In the present paper, continuum electrostatic calculations were employed to evaluate the electrostatic potential, energies, and protonation state of bovine heart cytochrome c oxidase for different redox states of the enzyme along its catalytic cycle. Three different computational models of the enzyme were employed to test the stability of the results. The energetics and pH dependence of the P-->F, F-->O, and O-->E steps of the cycle have been investigated. On the basis of electrostatic calculations, two possible schemes of redox-linked proton pumping are discussed. The first scheme involves His291 as a pump element, whereas the second scheme involves a group linked to propionate D of heme a(3). In both schemes, loading of the pump site is coupled to ET between the two hemes of the enzyme, while transfer of a chemical proton is accompanied by ejection of the pumped H(+). The two models, as well as the energetics results are compared with recent experimental kinetic data. The proton pumping across the membrane is an endergonic process, which requires a sufficient amount of energy to be provided by the chemical reaction in the active site. In our calculations, the conversion of OH(-) to H(2)O provides 520 meV of energy to displace pump protons from a loading site and overall about 635 meV for each electron passing through the system. Assuming that the two charges are translocated per electron against the membrane potential of 200 meV, the model predicts an overall efficiency of 63%.

Redox potential of quinones in both electron transfer branches of photosystem I

The redox potentials of the two electron transfer (ET) active quinones in the central part of photosystem I (PSI) were determined by evaluating the electrostatic energies from the solution of the Poisson-Boltzmann equation based on the crystal structure. The calculated redox potentials are -531 mV for A1A and -686 mV for A1B. From these results we conclude the following. (i) Both branches are active with a much faster ET in the B-branch than in the A-branch. (ii) The measured lifetime of 200-290 ns of reduced quinones agrees with the estimate for the A-branch and corroborates with an uphill ET from this quinone to the iron-sulfur cluster as observed in recent kinetic measurements. (iii) The electron paramagnetic resonance spectroscopic data refer to the A-branch quinone where the corresponding ET is uphill in energy. The negative redox potential of A1 in PSI is primarily because of the influence from the negatively charged FX, in contrast to the positive shift on the quinone redox potential in bacterial reaction center and PSII that is attributed to the positively charged non-heme iron atom. The conserved residue Asp-B575 changes its protonation state after quinone reduction. The difference of 155 mV in the quinone redox potentials of the two branches were attributed to the conformation of the backbone with a large contribution from Ser-A692 and Ser-B672 and to the side chain of Asp-B575, whose protonation state couples differently with the formation of the quinone radicals.

Tuning heme redox potentials in the cytochrome C subunit of photosynthetic reaction centers

The photosynthetic reaction center (RC) from Rhodopseudomonas viridis contains four cytochrome c hemes. They establish the initial part of the electron transfer (ET) chain through the RC. Despite their chemical identity, their midpoint potentials cover an interval of 440 mV. The individual heme midpoint potentials determine the ET kinetics and are therefore tuned by specific interactions with the protein environment. Here, we use an electrostatic approach based on the solution of the linearized Poisson-Boltzmann equation to evaluate the determinants of individual heme redox potentials. Our calculated redox potentials agree within 25 meV with the experimentally measured values. The heme redox potentials are mainly governed by solvent accessibility of the hemes and propionic acids, by neutralization of the negative charges at the propionates through either protonation or formation of salt bridges, by interactions with other hemes, and to a lesser extent, with other titratable protein side chains. In contrast to earlier computations on this system, we used quantum chemically derived atomic charges, considered an equilibrium-distributed protonation pattern, and accounted for interdependencies of site-site interactions. We provide values for the working potentials of all hemes as a function of the solution redox potential, which are crucial for calculations of ET rates. We identify residues whose site-directed mutation might significantly influence ET processes in the cytochrome c part of the RC. Redox potentials measured on a previously generated mutant could be reproduced by calculations based on a model structure of the mutant generated from the wild type RC.

Investigation of protonatable residues in Rhodothermus marinus caa3 haem-copper oxygen reductase: comparison with Paracoccus denitrificans aa3 haem-copper oxygen reductase

The Rhodothermus marinus caa(3 )haem-copper oxygen reductase contains all the residues of the so-called D- and K-proton channels, with the notable exception of the helix VI glutamate residue (Glu278(I) in Paracoccus denitrificans aa(3)), being nevertheless a true oxygen reductase reducing O(2) to water, and an efficient proton pump. Instead, in the same helix, but one turn below, it has a tyrosine residue (Tyr256(I), R. marinus caa(3) numbering), whose hydroxyl group occupies the same spatial position as the carboxylate group of Glu278(I), as deduced by comparative modelling techniques. Therefore, we proposed previously that this tyrosine residue could play an important role in the proton pathway. In this article we further study this hypothesis, by investigating the equilibrium thermodynamics of protonation in R. marinus caa(3), using theoretical methodologies based on the structural model previously obtained. Control calculations are also performed for the P. denitrificans aa(3) oxygen reductase. In both oxygen reductases we find several residues that are proton active (i.e., that display partial protonation) at physiological pH, some of them being redox sensitive (i.e., sensitive to the protein redox state). However, the caa(3 )Tyr256(I) is not proton active at physiological pH, in contrast to the aa(3) Glu278(I) which is both proton active at physiological pH and shows a high redox sensitivity. In R. marinus caa(3) we do not find any other residues in the same protein zone that can have this property. Therefore, there are no putative D-channel residues that are proton active in this oxidase. The protonatable residues of the K-channel are much more functionally conserved in both oxygen reductases than the same type of residues in the D-channel. Two (Tyr262(I) and Lys336(I), caa(3) numbering) out of three protonatable K-channel residues are proton active and redox sensitive in both proteins.

Molecular basis for redox-Bohr and cooperative effects in cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774: Crystallographic and modeling studies of oxidized and reduced high-resolution structures at pH 7.6

The tetraheme cytochrome c3 is a small metalloprotein with ca. 13,000 Da found in sulfate-reducing bacteria, which is believed to act as a partner of hydrogenase. The three-dimensional structure of the oxidized and reduced forms of cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774 at pH 7.6 were determined using high-resolution X-ray crystallography and were compared with the previously determined oxidized form at pH 4.0. Theoretical calculations were performed with both structures, using continuum electrostatic calculations and Monte Carlo sampling of protonation and redox states, in order to understand the molecular basis of the redox-Bohr and cooperativity effects related to the coupled transfer of electrons and protons. We were able to identify groups that showed redox-linked conformational changes. In particular, Glu61, His76, and propionate D of heme II showed important contributions to the redox-cooperativity, whereas His76, propionate A of heme I, and propionate D of heme IV were the key residues for the redox-Bohr effect. Upon reduction, an important movement of the backbone region surrounding hemes I and II was also identified, that, together with a few redox-linked conformational changes in side-chain residues, results in a significant decrease in the solvent accessibility of hemes I and II.

Identification of residues crucially involved in the binding of the heme moiety of soluble guanylate cyclase

Soluble guanylate cyclase (sGC), a heterodimeric hemeprotein, is the only receptor for the biological messenger nitric oxide (NO) identified to date and is intimately involved in various signal transduction pathways. By using the recently discovered NO- and heme-independent sGC activator BAY 58-2667 and a novel cGMP reporter cell, we could distinguish between heme-containing and heme-free sGC in an intact cellular system. Using these novel tools, we identified the invariant amino acids tyrosine 135 and arginine 139 of the beta(1)-subunit as crucially important for both the binding of the heme moiety and the activation of sGC by BAY 58-2667. The heme is displaced by BAY 58-2667 due to a competition between the carboxylic groups of this compound and the heme propionic acids for the identified residues tyrosine 135 and arginine 139. This displacement results in the release of the axial heme ligand histidine 105 and to the observed activation of sGC. Based on these findings we postulate a signal transmission triad composed of histidine 105, tyrosine 135, and arginine 139 responsible for the enzyme activation by this compound and probably also for transducing changes in heme status and porphyrin geometry upon NO binding into alterations of sGC catalytic activity.

Redox-Bohr and other cooperativity effects in the nine-heme cytochrome C from Desulfovibrio desulfuricans ATCC 27774: crystallographic and modeling studies

The nine-heme cytochrome c is a monomeric multiheme cytochrome found in Desulfovibrio desulfuricans ATCC 27774. The polypeptide chain comprises 296 residues and wraps around nine hemes of type c. It is believed to take part in the periplasmic assembly of proteins involved in the mechanism of hydrogen cycling, receiving electrons from the tetraheme cytochrome c3. With the purpose of understanding the molecular basis of electron transfer processes in this cytochrome, we have determined the crystal structures of its oxidized and reduced forms at pH 7.5 and performed theoretical calculations of the binding equilibrium of protons and electrons in these structures. This integrated study allowed us to observe that the reduction process induced relevant conformational changes in several residues, as well as protonation changes in some protonatable residues. In particular, the surroundings of hemes I and IV constitute two areas of special interest. In addition, we were able to ascertain the groups involved in the redox-Bohr effect present in this cytochrome and the conformational changes that may underlie the redox-cooperativity effects on different hemes. Furthermore, the thermodynamic simulations provide evidence that the N- and C-terminal domains function in an independent manner, with the hemes belonging to the N-terminal domain showing, in general, a lower redox potential than those found in the C-terminal domain. In this way, electrons captured by the N-terminal domain could easily flow to the C-terminal domain, allowing the former to capture more electrons. A notable exception is heme IX, which has low redox potential and could serve as the exit path for electrons toward other proteins in the electron transfer pathway.

Redox potential of quinones in photosynthetic reaction centers from Rhodobacter sphaeroides: dependence on protonation of Glu-L212 and Asp-L213

The absolute values of the one-electron redox potentials of the two quinones (Q(A) and Q(B)) in bacterial photosynthetic reaction centers from Rhodobacter sphaeroides were calculated by evaluating the electrostatic energies from the solution of the linearized Poisson-Boltzmann equation at pH 7.0. The redox potential for Q(A) was calculated to be between -173 and -160 mV, which is close to the lowest measured values that are assumed to refer to nonequilibrated protonation patterns in the redox state Q(A)(-). The redox potential of quinone Q(B) is found to be about 160-220 mV larger for the light-exposed than for the dark-adapted structure. These values of the redox potentials are obtained if Asp-L213 is nearly protonated (probability 0.75-1.0) before and after electron transfer from Q(A) to Q(B), while Glu-L212 is partially protonated (probability 0.6) in the initial state Q(A)(-)Q(B)(0) and fully protonated in the final state Q(A)(0)Q(B)(-). Conversely, if the charge state of the quinones is varied from Q(A)(-)Q(B)(0) to Q(A)(0)Q(B)(-) corresponding to the electron transfer from Q(A) to Q(B), Asp-L213 remains protonated, while Glu-L212 changes its protonation state from 0.15 H(+) to fully protonated. In agreement with results from FTIR spectra, there is proton uptake at Glu-L212 going along with the electron transfer, whereas Asp-L213 does not change its protonation state. However, in our simulations Asp-L213 is considered to be protonated rather than ionized as deduced from FTIR spectra. The calculated redox potential of Q(A) shows little dependence on the charge state of Asp-L213, which is due to a strong coupling with the protonation state of Asp-M17 but increases by 50 mV if Glu-L212 changes from the ionized to the protonated charge state. Both are in agreement with fluorescence measurements observing the decay of SP(+)Q(A)(-) in a wide pH regime. The computed difference in redox potential of Q(B) in the light-exposed and dark-adapted structure was traced back to the hydrogen bond of Q(B) with His-L190 that is lost in the dark-adapted structure and the charge of the non-heme iron atom, which is closer to Q(B) in the light-exposed than in the dark-adapted structure.

Energetics of radical transfer in DNA photolyase

Charge separation and radical transfer in DNA photolyase from Escherichia coli is investigated by computing electrostatic free energies from a solution of the Poisson-Boltzmann equation. For the initial charge separation 450 meV are available. According to recent experiments [Aubert et al. Nature 2000, 405, 586-590] the flavin receives an electron from the proximal tryptophan W382, which consequently forms a cationic radical WH(*)(+)382. The radical state is subsequently transferred along the triad W382-W359-W306 of conserved tryptophans. The radical transfer to the intermediate tryptophan W359 is nearly isoenergetic (58 meV uphill); the radical transfer from the intermediate W359 to the distal W306 is 200 meV downhill in energy, funneling and stabilizing the radical state at W306. The resulting cationic radical WH(*)(+)306 is further stabilized by deprotonation, yielding the neutral radical W(*)306, which is 214 meV below WH(*)(+)306. The time scale of the charge recombination process yielding back the resting enzyme with FADH(*) is governed by reprotonation of W306, with a calculated lifetime of 1.2 ms that correlates well with the measured lifetime of 17 ms. In photolyase from Anacystis nidulans the radical state is partially transferred to a tyrosine [Aubert et al. Proc. Natl. Acad. Sci. U.S.A. 1999, 96, 5423-5427]. In photolyase from Escherichia coli, there is a tyrosine (Y464) close to the distal tryptophan W306 that could play this role. We show that this tyrosine cannot be involved in radical transfer, because the electron transfer from tyrosine to W306 is much too endergonic (750 meV) and a direct hydrogen transfer is likely too slow. Coupling of specific charge states of the tryptophan triad with protonation patterns of titratable residues of photolyase is small.