No abstract

Hybrid density functional theory using the B3LYP functional has been applied to study the energetics of proton translocation in cytochrome oxidase. Redox potentials of the metal centers and the tyrosyl radical have been computed, as well as pK a values of important groups along the translocation path. Models with more than 100 atoms have been used with geometries fully optimized in the presence of a dielectric continuum to represent the surrounding enzyme. The models are built from the bovine X-ray structure (2OCC), which does not contain water molecules. The main result obtained is that for these models the pK a values of propionate A of heme a3 are actually larger than those of the oxo and hydroxyl groups of the binuclear center, even after an electron has been transferred to heme a3 from heme a. The consequence of this rather surprising finding is that a proton coming from the inside of the membrane will actually thermodynamically prefer to go to the propionate for further pumping rather than go to the binuclear center for consumption in the oxygen chemistry. Full energetic cycles are constructed for both the case without and the case with a proton gradient across the membrane.

The mechanism of the nitric oxide reduction in a bacterial nitric oxide reductase (NOR) has been investigated in two model systems of the heme-b(3)-Fe(B) active site using density functional theory (B3LYP). A model with an octahedral coordination of the non-heme Fe(B) consisting of three histidines, one glutamate and one water molecule gave an energetically feasible reaction mechanism. A tetrahedral coordination of the non-heme iron, corresponding to the one of Cu(B) in cytochrome oxidase, gave several very high barriers which makes this type of coordination unlikely. The first nitric oxide coordinates to heme b(3) and is partly reduced to a more nitroxyl anion character, which activates it toward an attack from the second NO. The product in this reaction step is a hyponitrite dianion coordinating in between the two irons. Cleaving an NO bond in this intermediate forms an Fe(B) (IV)O and nitrous oxide, and this is the rate determining step in the reaction mechanism. In the model with an octahedral coordination of Fe(B) the intrinsic barrier of this step is 16.3 kcal/mol, which is in good agreement with the experimental value of 15.9 kcal/mol. However, the total barrier is 21.3 kcal/mol, mainly due to the endergonic reduction of heme b(3) taken from experimental reduction potentials. After nitrous oxide has left the active site the ferrylic Fe(B) will form a mu-oxo bridge to heme b(3) in a reaction step exergonic by 45.3 kcal/mol. The formation of a quite stable mu-oxo bridge between heme b(3) and Fe(B) is in agreement with this intermediate being the experimentally observed resting state in oxidized NOR. The formation of a ferrylic non-heme Fe(B) in the proposed reaction mechanism could be one reason for having an iron as the non-heme metal ion in NOR instead of a Cu as in cytochrome oxidase.

The mechanism for the reaction between nitric oxide (NO) and O(2) bound to the heme iron of myoglobin (Mb), including the following isomerization to nitrate, has been investigated using hybrid density functional theory (B3LYP). Myoglobin working as a NO scavenger could be of importance, since NO reversibly inhibits the terminal enzyme in the respiration chain, cytochrome c oxidase. The concentration of NO in the cell will thus affect the respiration and thereby the synthesis of ATP. The calculations show that the reaction between NO and the heme-bound O(2) gives a peroxynitrite intermediate whose O-O bond undergoes a homolytic cleavage, forming a NO(2) radical and myoglobin in the oxo-ferryl state. The NO(2) radical then recombines with the oxo-ferryl, forming heme-bound nitrate. Nine different models have been used in the present study to examine the effect on the reaction both by the presence and the protonation state of the distal His64, and by the surroundings of the proximal His93. The barriers going from the oxy-Mb and nitric oxide reactant to the peroxynitrite intermediate and further to the oxo-ferryl and NO(2) radical are around 10 and 7 kcal/mol, respectively. Forming the product, nitrate bound to the heme iron has a barrier of less than approximately 7 kcal/mol. The overall reaction going from a free nitric oxide and oxy-Mb to the heme bound nitrate is exergonic by more than 30 kcal/mol.

The mechanism for the reduction of nitric oxide to nitrous oxide and water in an A-type flavoprotein (FprA) in Moorella thermoacetica, which has been proposed to be a scavenging type of nitric oxide reductase, has been investigated using density functional theory (B3LYP). A dinitrosyl complex, [{FeNO}(7)](2), has previously been proposed to be a key intermediate in the NO reduction catalyzed by FprA. The electrons and protons involved in the reduction were suggested to "super-reduce" the dinitrosyl intermediate to [{FeNO}(8)](2) or the corresponding diprotonated form, [{FeNO(H)}(8)](2). In this type of mechanism the electron and/or proton transfers will be a part of the rate-determining step. In the present study, on the other hand, a reaction mechanism is suggested in which N(2)O can be formed before the protons and electrons enter the catalytic cycle. One of the irons in the diiron center is used to stabilize the formation of a hyponitrite dianion, instead of binding a second NO. Cleaving the N-O bond in the hyponitrite dianion intermediate is the rate-determining step in the proposed reaction mechanism. The barrier of 16.5 kcal mol(-1) is in good agreement with the barrier height of the experimental rate-determining step of 14.8 kcal mol(-1). The energetics of some intermediates in the "super-reduction" mechanism and the mechanism proceeding via a hyponitrite dianion are compared, favoring the latter. It is also discussed how to experimentally discriminate between the two mechanisms.