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%.

Using electrostatic calculations, we have examined the dependence of the protonation state of cytochrome c oxidase from bovine heart on its redox state. Based on these calculations, we propose a possible scheme of redox-linked proton pumping. The scheme involves His291 - one of the ligands of the Cu(B) redox center - which plays the role of the proton loading site (PLS) of the pump. The mechanism of pumping is based on ET reaction between two hemes of the enzyme, which is coupled to a transfer of two protons. Upon ET, the first proton (fast reaction) is transferred to the PLS (His291), while subsequent transfer of the second "chemical" proton to the binuclear center (slow reaction) is accompanied by the ejection of the first (pumped) proton. Within the proposed model, we discuss the catalytic cycle of the enzyme.

Cation-pi interactions between aromatic residues and cationic amino groups in side chains and have been recognized as noncovalent bonding interactions relevant for molecular recognition and for stabilization and definition of the native structure of proteins. We propose a novel type of cation-pi interaction in metalloproteins; namely interaction between ligands coordinated to a metal cation--which gain positive charge from the metal--and aromatic groups in amino acid side chains. Investigation of crystal structures of metalloproteins in the Protein Data Bank (PDB) has revealed that there exist quite a number of metalloproteins in which aromatic rings of phenylalanine, tyrosine, and tryptophan are situated close to a metal center interacting with coordinated ligands. Among these ligands are amino acids such as asparagine, aspartate, glutamate, histidine, and threonine, but also water and substrates like ethanol. These interactions play a role in the stability and conformation of metalloproteins, and in some cases may also be directly involved in the mechanism of enzymatic reactions, which occur at the metal center. For the enzyme superoxide dismutase, we used quantum chemical computation to calculate that Trp163 has an interaction energy of 10.09 kcal mol(-1) with the ligands coordinated to iron.

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.

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.

Factors determining conformations of imidazole axially coordinated to heme in heme proteins were investigated by analyzing 693 hemes in 432 different crystal structures of heme proteins from the Protein Data Bank (PDB), where at least one histidine is ligated to heme. The results from a search of the PDB for protein structures were interpreted with molecular force field computations. Analysis of data from these crystal structures indicated that there are two main factors that determine the orientations of imidazole ligated to heme. These are the interactions of imidazole with the propionic acid side chains of heme and with the histidine backbone. From the analysis of the crystal structures of heme proteins, it turned out that the hydrogen bonding pattern is often not decisive, though it is probably used by nature to fine-tune the orientation of imidazole axially ligated to heme. We found that in many heme proteins the NdeltaH group of imidazole ligated to heme can assume a number of different hydrogen bonds and that in mutant structures the orientation of the ligated imidazole often does not change significantly, although the mutant altered the hydrogen bonding scheme involving the imidazole. Data from crystal structures of heme proteins show that there are preferred orientations of imidazoles with respect to heme. Generally, the NdeltaH group of imidazole is oriented toward the propionic acid groups of the heme. In some cases, the NdeltaH group of imidazole is close to only one of the propionic acid groups, but it is practically never oriented in the opposite direction. The imidazole also adopts a preferred orientation with respect to its histidine backbone such that the plane of the imidazole ring is practically never parallel to the Calpha-Cbeta bond of its histidine backbone. For a given conformation of histidine backbone with respect to heme, as well as imidazole with respect to histidine backbone, the orientation of the imidazole with respect to heme is uniquely determined, since the three orientations depend on each other. Hence, the interaction of the imidazole with the backbone also influences the orientation of the imidazole with respect to the heme. Force field computations are in agreement with experimental data. With this method, we showed that there is an energy minimum when the NdeltaH group of the imidazole is oriented toward the propionic acid groups and that there are minima of energy for orientations where the imidazole ring is orthogonal to the plane defined by the Calpha-Cbeta and Cbeta-Cgamma bonds of the histidine. The computations also demonstrated that these interactions are mainly of electrostatic origin. By taking into account these two major factors, we were able to understand the orientations of axially coordinated imidazoles for all groups of heme proteins, except for the group of cytochrome c peroxidase. In this group, the orientation of the imidazole is determined by a strong hydrogen bond of the NdeltaH group with Asp235.

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.

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.