A structural and functional perspective on the evolution of the heme–copper oxidases

The heme–copper oxidases (HCOs) catalyze the reduction of O₂ to water, and couple the free energy to proton pumping across the membrane. HCOs are divided into three sub-classes, A, B and C, whose order of emergence in evolution has been controversial. Here we have analyzed recent structural and functional data on HCOs and their homologues, the nitric oxide reductases (NORs). We suggest that the C-type oxidases are ancient enzymes that emerged from the NORs. In contrast, the A-type oxidases are the most advanced from both structural and functional viewpoints, which we interpret as evidence for having evolved later.     

Functional proton transfer pathways in the heme–copper oxidase superfamily

Heme–copper oxidases (HCuOs) terminate the respiratory chain in mitochondria and most bacteria. They are transmembrane proteins that catalyse the reduction of oxygen and use the liberated free energy to maintain a proton-motive force across the membrane. The HCuO superfamily has been divided into the oxygen-reducing A-, B- and C-type oxidases as well as the bacterial NO reductases (NOR), catalysing the reduction of NO in the denitrification process. Proton transfer to the catalytic site in the mitochondrial-like A family occurs through two well-defined pathways termed the D- and K-pathways. The B, C, and NOR families differ in the pathways as well as the mechanisms for proton transfer to the active site and across the membrane. Recent structural and functional investigations, focussing on proton transfer in the B, C and NOR families will be discussed in this review. This article is part of a Special Issue entitled: Respiratory Oxidases.     

Roles for heme–copper oxidases in extreme high-light and oxidative stress response in the cyanobacterium Synechococcus sp. PCC 7002

The ctaCIDIEI and ctaCIIDIIEII gene clusters that encode heme–copper cytochrome oxidases have been characterized in the marine cyanobacterium Synechococcus sp. PCC 7002 and the inactivation of ctaDI was shown to affect high-light adaptation. In this study, Synechococcus sp. PCC 7002 wild-type, ctaDI, ctaDII, and ctaDI–ctaDII double mutants were grown under extreme high-light and oxidative stress to further assess the roles of cytochrome oxidases in cyanobacteria. Cells of the ctaDI mutant strain barely grew under extreme high-light illumination of 4.5 mE m⁻² s⁻¹, suggesting that CtaDI is required for high-light acclimation in Synechococcus sp. PCC 7002. The ctaDI–ctaDII double mutant cells unexpectedly tolerated extreme high-light intensity, indicating that the disruption of ctaDII gene suppresses the high-light sensitivity phenotype of the ctaDI single mutant. The ctaDII mutant cells also exhibited higher tolerance to the oxidative stress compound, methyl viologen, in the growth media. The ctaDII mutant and the ctaDI–ctaDII double mutant cells had approximately twofold higher levels of superoxide dismutase (SOD) activity, indicating that the disruption of ctaDII gene increased the capacity to decompose active oxygen species. These results suggest that the CtaII cytochrome oxidase may be involved with the oxidative stress response, including the control of SOD expression.     

A realistic in silico model for structure/function studies of molybdenum–copper CO dehydrogenase

CO dehydrogenase (CODH) is an environmentally crucial bacterial enzyme that oxidizes CO to CO₂ at a Mo–Cu active site. Despite the close to atomic resolution structure (1.1 Å), significant uncertainties have remained with regard to the protonation state of the water-derived equatorial ligand coordinated at the Mo-center, as well as the nature of intermediates formed during the catalytic cycle. To address the protonation state of the equatorial ligand, we have developed a realistic in silico QM model (~179 atoms) containing structurally essential residues surrounding the active site. Using our QM model, we examined each plausible combination of redox states (Mo^VI–Cu^I, Mo^V–Cu^II, Mo^V–Cu^I, and Mo^IV–Cu^I) and Mo-coordinated equatorial ligands (O^2?, OH^?, H₂O), as well as the effects of second-sphere residues surrounding the active site. Herein, we present a refined computational model for the Mo(VI) state in which Glu763 acts as an active site base, leading to a MoO₂-like core and a protonated Glu763. Calculated structural and spectroscopic data (hyperfine couplings) are in support of a MoO₂-like core in agreement with XRD data. The calculated two-electron reduction potential (E = ?467 mV vs. SHE) is in reasonable agreement with the experimental value (E = ?558 mV vs. SHE) for the redox couple comprising an equatorial oxo ligand and protonated Glu763 in the Mo^VI–Cu^I state and an equatorial water in the Mo^IV–Cu^I state. We also suggest a potential role of second-sphere residues (e.g., Glu763, Phe390) based on geometric changes observed upon exclusion of these residues in the most plausible oxidized states.     

Investigation into the mechanism of λmax shifts and their dependence on pH for the 2-hydrazinopyridine derivatives of two copper amine oxidases

Copper amine oxidases (CuAO), from Escherichia coli (ECAO) and pea seedling (PSAO) were reacted with an excess of the hydrazine derivative 2-hydrazinopyridine (2HP) to form an initial, strongly absorbing adduct, (adduct 1; λ_max 420–430 nm) formed by the covalent binding of 2HP with the active site cofactor 2,4,5-trihydroxyphenylalanine quinone (TPQ). Thermal incubation of buffered solutions of adduct 1 (pH 5.65–10.7) or addition of KOH solution (giving a final pH of 13–15) led isosbestically to a dramatic λ_max shift yielding adduct 2 (λ_max 520–530 nm). For both ECAO and PSAO, an increase in pH resulted in increased formation of adduct 2 with concomitant loss of adduct 1. Maximum adduct 2 formation occurred at pH 9.84 in ECAO and at pH 10.7 in PSAO. Beyond these pH levels, adduct 2 formation occurred to a much lesser extent which was independent of pH, suggesting enzyme denaturation. It is proposed that the conversion of adduct 1 to adduct 2 occurs as a result of hydrazone to azo conversion mediated by loss of a single proton, possibly to the active site base. It is further postulated that adduct formation and subsequent deprotonation can be likened to the substrate and product Schiff base complexes in the reductive half cycle of copper/TPQ containing amine oxidases. As part of this study an extinction coefficient at 280 nm was determined for ECAO by gravimetric analysis. This yielded a value of 2.1×10⁵ M⁻¹ cm⁻¹ giving rise to the need of a correction factor when estimating the protein concentration from an absorbance reading at 280 nm. Using the estimated molecular mass of 160 kDa for the homodimeric ECAO, a correction factor of 0.76 must be applied.     

A minimal model for the mitochondrial rapid mode of Ca²+ uptake mechanism

Mitochondria possess a remarkable ability to rapidly accumulate and sequester Ca²⁺. One of the mechanisms responsible for this ability is believed to be the rapid mode (RaM) of Ca²⁺ uptake. Despite the existence of many models of mitochondrial Ca²⁺ dynamics, very few consider RaM as a potential mechanism that regulates mitochondrial Ca²⁺ dynamics. To fill this gap, a novel mathematical model of the RaM mechanism is developed herein. The model is able to simulate the available experimental data of rapid Ca²⁺ uptake in isolated mitochondria from both chicken heart and rat liver tissues with good fidelity. The mechanism is based on Ca²⁺ binding to an external trigger site(s) and initiating a brief transient of high Ca²⁺ conductivity. It then quickly switches to an inhibited, zero-conductive state until the external Ca²⁺ level is dropped below a critical value (∼100–150 nM). RaM''s Ca²⁺- and time-dependent properties make it a unique Ca²⁺ transporter that may be an important means by which mitochondria take up Ca²⁺ in situ and help enable mitochondria to decode cytosolic Ca²⁺ signals. Integrating the developed RaM model into existing models of mitochondrial Ca²⁺ dynamics will help elucidate the physiological role that this unique mechanism plays in mitochondrial Ca²⁺-homeostasis and bioenergetics.     

A model for the rate of enzymatic hydrolysis of cellulose in heterogeneous solid–liquid systems

Hydrolysis of cellulose by cellulase enzymes has been studied in a stirred batch reactor at 50°C. A kinetic model has been devised by which the behaviour of such a reaction could be described. The model has been developed on the basis of shrinking particle theory and Langmuir isotherm concept. The applicability of the model has been tested by comparing the experimental results for diverse reaction systems, obtained in the present study or taken from the literature, and those predicted from the model. The degree of agreement was within ±2–11%.     

Structure and coordination of CuB in the Acidianus ambivalens aa 3 quinol oxidase heme–copper center

The coordination environment of the Cu_B center of the quinol oxidase from Acidianus ambivalens, a type B heme–copper oxygen reductase, was investigated by Fourier transform (FT) IR and extended X-ray absorption fine structure (EXAFS) spectroscopy. The comparative structural chemistry of dinuclear Fe–Cu sites of the different types of oxygen reductases is of great interest. Fully reduced A. ambivalens quinol oxidase binds CO at the heme a ₃ center, with ν(CO)=1,973 cm⁻¹. On photolysis, the CO migrated to the Cu_B center, forming a Cu_B^I–CO complex with ν(CO)=2,047 cm⁻¹. Raising the temperature of the samples to 25°C did not result in a total loss of signal in the FTIR difference spectrum although the intensity of these signals was reduced sevenfold. This observation is consistent with a large energy barrier against the geminate rebinding of CO to the heme iron from Cu_B, a restricted limited access at the active-site pocket for a second binding, and a kinetically stable Cu_B–CO complex in A. ambivalens aa ₃. The Cu_B center was probed in a number of different states using EXAFS spectroscopy. The oxidized state was best simulated by three histidines and a solvent O scatterer. On reduction, the site became three-coordinate, but in contrast to the bo ₃ enzyme, there was no evidence for heterogeneity of binding of the coordinated histidines. The Cu_B centers in both the oxidized and the reduced enzymes also appeared to contain substoichiometric amounts (0.2 mol equiv) of nonlabile chloride ion. EXAFS data of the reduced carbonylated enzyme showed no difference between dark and photolyzed forms. The spectra could be well fit by 2.5 imidazoles, 0.5 Cl⁻ and 0.5 CO ligands. This arrangement of scatterers would be consistent with about half the sites remaining as unligated Cu(his)₃ and half being converted to Cu(his)₂Cl⁻CO, a 50/50 ratio of Cu(his)₂Cl⁻ and Cu(his)₃CO, or some combination of these formulations. Electronic Supplementary Material Supplementary material is available for this article at .     

Protection of the pancreatic β-cell glucoreceptor mechanism for insulin secretion during culture in chemically defined medium

High concentrations of glucose have a protective effect on the glucoreceptor mechanism for insulin secretion during culture of pancreatic islets in chemically defined media. To study at what level glucose exerts this effect, insulin secretion from beta-cell-rich mouse pancreatic islets was measured before and after culture for 1 week in the presence of different substances. Before culture, glucose and inosine were potent stimulators, mannose and fructose were less potent and xylitol had no effect on secretion. Culture in 3mm-glucose resulted in a 10-fold decrease in the insulin response to glucose stimulation. A less marked decrease was noted after culture in 20mm- or 30mm-glucose. Inosine-stimulated secretion was much decreased after culture in high concentrations of glucose, whereas the responses to mannose or fructose were unchanged. After culture in 30mm-mannose, glucose-stimulated secretion was similar to that observed after culture in high concentrations of glucose, whereas the response to mannose had much decreased. There were no secretory responses to glucose or fructose after culture in 30mm-fructose, or to glucose or xylitol after culture in 30mm-xylitol. Culture in 10mm-inosine did not preserve any significant response to glucose or inosine. The insulin contents of islets and culture media were higher after culture in high concentrations of glucose, mannose or inosine than after culture in fructose, xylitol or low concentrations of glucose. It is suggested that glucose, and to some extent mannose, preserves the glucoreceptor mechanism for insulin secretion by influencing an early stage in glucose metabolism, presumably glucokinase activity.     

Reductive activation of the heme iron–nitrosyl intermediate in the reaction mechanism of cytochrome c nitrite reductase: a theoretical study

Cytochrome c nitrite reductase catalyzes the six-electron, seven-proton reduction of nitrite to ammonia without release of any detectable reaction intermediate. This implies a unique flexibility of the active site combined with a finely tuned proton and electron delivery system. In the present work, we employed density functional theory to study the recharging of the active site with protons and electrons through the series of reaction intermediates based on nitrogen monoxide [Fe(II)-NO(+), Fe(II)-NO·, Fe(II)-NO(-), and Fe(II)-HNO]. The activation barriers for the various proton and electron transfer steps were estimated in the framework of Marcus theory. Using the barriers obtained, we simulated the kinetics of the reduction process. We found that the complex recharging process can be accomplished in two possible ways: either through two consecutive proton-coupled electron transfers (PCETs) or in the form of three consecutive elementary steps involving reduction, PCET, and protonation. Kinetic simulations revealed the recharging through two PCETs to be a means of overcoming the predicted deep energetic minimum that is calculated to occur at the stage of the Fe(II)-NO· intermediate. The radical transfer role for the active-site Tyr(218), as proposed in the literature, cannot be confirmed on the basis of our calculations. The role of the highly conserved calcium located in the direct proximity of the active site in proton delivery has also been studied. It was found to play an important role in the substrate conversion through the facilitation of the proton transfer steps.     

A model for the evolution of crèching behaviour in gulls

A crèche is an aggregation of chicks outside nesting territories, within chicks continue to be fed only by their own parents. Several adaptive functions of crèching have been proposed, the most frequent being a reduction in predator pressure. Using an evolutionary stable strategy approach based on the computation of individuals' fecundity, we examined which regime of aerial and terrestrial predation is likely to favour the evolution and stability of the crèching strategy (CS) in gulls. Our results confirm the hypothesis that habitat instability associated with high levels of terrestrial predation favours the evolution and maintenance of crèching behaviour. Moreover, our results suggest that a low aggressiveness against predators may be a pre-adaptation to a CS. In contrast, the high synchronisation often observed in crèching species does not favour the evolution of a crèching behaviour and is thus probably under selection pressures different from those modelled here.     

A possible rôle for microsomal oxidases in metamorphosis and reproduction in the housefly

Housefly larvae and adults of the Isolan-B strain, insecticide resistant because of high microsomal oxidase activity, and the susceptible WHO Standard Reference (SR) strain, were reared on diets containing phenobarbital or piperonyl butoxide. Periodically throughout the treatments groups of insects were assayed for microsomal aldrin epoxidase activity. This activity was compared with such growth and development parameters as pupation, adult emergence, and reproduction in similar groups of insects. At diet levels of 0·01 to 0·50 per cent both chemicals caused large increases, as much as elevenfold in microsomal oxidase activity in third instar larvae and 15 to 100 per cent inhibition of pupation and emergence. In the adult diet at 1 per cent, both compounds caused at least 50 per cent decrease in egg production. Phenobarbital enhanced the enzyme system in both larval and adult stages of both strains but piperonyl butoxide, while enhancing enzyme activity in larvae of the Isolan-B strain, was an inhibitor in the adult stage. The results are interpreted as an indication of a direct connexion between microsomal oxidase activity and the action of hormones in the housefly, probably through regulation of hormone titre by these enzymes.     

Evolution of quinol oxidation within the heme‑copper oxidoreductase superfamily

The heme?copper oxidoreductase (HCO) superfamily is a large superfamily of terminal respiratory enzymes that are widely distributed across the three domains of life. The superfamily includes biochemically diverse oxygen reductases and nitric oxide reductases that are pivotal in the pathways of aerobic respiration and denitrification. The adaptation of HCOs to use quinol as the electron donor instead of cytochrome c has significant implication for the respiratory flexibility and energetic efficiency of the respiratory chains that include them. In this work, we explore the adaptation of this scaffold to two different electron donors, cytochromes c and quinols, with extensive sequence analysis of these enzymes from publicly available datasets. Our work shows that quinol oxidation evolved independently within the HCO superfamily at least seven times. Enzymes from only two of these independently evolved clades have been biochemically well-characterized. Combining structural modeling with sequence analysis, we identify putative quinol binding sites in each of the novel quinol oxidases. Our analysis of experimental and modeling data suggests that the quinol binding site appears to have evolved at the same structural position within the scaffold more than once.     

A phenomenological cohesive model for the macroscopic simulation of cell–matrix adhesions

Cell adhesion is crucial for cells to not only physically interact with each other but also sense their microenvironment and respond accordingly. In fact, adherent cells can generate physical forces that are transmitted to the surrounding matrix, regulating the formation of cell–matrix adhesions. The main purpose of this work is to develop a computational model to simulate the dynamics of cell–matrix adhesions through a cohesive formulation within the framework of the finite element method and based on the principles of continuum damage mechanics. This model enables the simulation of the mechanical adhesion between cell and extracellular matrix (ECM) as regulated by local multidirectional forces and thus predicts the onset and growth of the adhesion. In addition, this numerical approach allows the simulation of the cell as a whole, as it models the complete mechanical interaction between cell and ECM. As a result, we can investigate and quantify how different mechanical conditions in the cell (e.g., contractile forces, actin cytoskeletal properties) or in the ECM (e.g., stiffness, external forces) can regulate the dynamics of cell–matrix adhesions.     

A bioreaction–diffusion model for growth of marine sponge explants in bioreactors

Marine sponges are sources of high-value bioactives. Engineering aspects of in vitro culture of sponges from cuttings (explants) are poorly understood. This work develops a diffusion-controlled growth model for sponge explants. The model assumes that the explant growth is controlled by diffusive transport of at least some nutrients from the surrounding medium into the explant that generally has a poorly developed aquiferous system for internal irrigation during early stages of growth. Growth is assumed to obey Monod-type kinetics. The model is shown to satisfactorily explain the measured growth behavior of the marine sponge Crambe crambe in two different growth media. In addition, the model is generally consistent with published data for growth of explants of the sponges Disidea avara and Hemimycale columella. The model predicted that nutrient concentration profiles for nutrients, such as dissolved oxygen within the explant, are consistent with data published by independent researchers. In view of the proposed model’s ability to explain available data for growth of several species of sponge explants, diffusive transport does play a controlling role in explant growth at least until a fully developed aquiferous system has become established. According to the model and experimental observations, the instantaneous growth rate depends on the size of the explant and all those factors that influence the diffusion of critical nutrients within the explant. Growth follows a hyperbolic profile that is consistent with the Monod kinetics.     

Looking for the minimum common denominator in haem–copper oxygen reductases: Towards a unified catalytic mechanism

Haem–copper oxygen reductases are transmembrane protein complexes that reduce dioxygen to water and pump protons across the mitochondrial or periplasmatic membrane, contributing to the transmembrane difference of electrochemical potential. Seven years ago we proposed a classification of these enzymes into three different families (A, B and C), based on the amino acid residues of their proton channels and amino acid sequence comparison, later supported by the so far identified characteristics of the catalytic centre of members from each family. The three families have in common the same general structural fold of the catalytic subunit, which contains the same or analogous prosthetic groups, and proton channels. These observations raise the hypothesis that the mechanisms for dioxygen reduction, proton pumping and the coupling of the two processes may be the same for all these enzymes. Under this hypothesis, they should be performed and controlled by the same or equivalent elements/events, and the identification of retained elements in all families will reveal their importance and may prompt the definition of the enzyme operating mode. Thus, we believe that the search for a minimum common denominator has a crucial importance, and in this article we highlight what is already established for the haem–copper oxygen reductases and emphasize the main questions still unanswered in a comprehensive basis.