Theoretical identification of proton channels in the quinol oxidase aa3 from Acidianus ambivalens

Coupling between protonation and conformation in cytochrome c oxidase: Insights from constant-pH MD simulations

Cytochrome c oxidases (CcOs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes reduce dioxygen (O(2)) to water and, simultaneously, generate a transmembrane electrochemical proton gradient. Despite their importance in the aerobic metabolism and the large amount of structural and biochemical data available for the A1-type CcO family, there is still no consensually accepted description of the molecular mechanisms operating in this protein. A substantial number of questions about the CcO's working mechanism remain to be answered, including how the protonation behavior of some key residues is modulated during a reduction cycle and how is the conformation of the protein affected by protonation. The main objective of this work was to study the protonation-conformation coupling in CcOs and identify the molecular factors that control the protonation state of some key residues. In order to directly capture the interplay between protonation and conformational effects, we have performed constant-pH MD simulations of an A1-type CcO inserted into a lipid bilayer in two redox states (oxidized and reduced) at physiological pH. From the simulations, we were able to identify several groups with unusual titration behavior that are highly dependent on the protein redox state, including the A-propionate from heme a and the D-propionate from heme a3, two key groups possibly involved in proton pumping. The protonation state of these two groups is heavily influenced by subtle conformational changes in the protein (notably of R481(I) and R482(I)) and by small changes in the hydrogen bond network.

Exploring O2 diffusion in A-type cytochrome c oxidases: molecular dynamics simulations uncover two alternative channels towards the binuclear site

Cytochrome c oxidases (Ccoxs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes couple dioxygen (O₂) reduction to the generation of a transmembrane electrochemical proton gradient. Despite decades of research and the availability of a large amount of structural and biochemical data available for the A-type Ccox family, little is known about the channel(s) used by O₂ to travel from the solvent/membrane to the heme a₃-Cu_B binuclear center (BNC). Moreover, the identification of all possible O₂ channels as well as the atomic details of O₂ diffusion is essential for the understanding of the working mechanisms of the A-type Ccox. In this work, we determined the O₂ distribution within Ccox from Rhodobacter sphaeroides, in the fully reduced state, in order to identify and characterize all the putative O₂ channels leading towards the BNC. For that, we use an integrated strategy combining atomistic molecular dynamics (MD) simulations (with and without explicit O₂ molecules) and implicit ligand sampling (ILS) calculations. Based on the 3D free energy map for O₂ inside Ccox, three channels were identified, all starting in the membrane hydrophobic region and connecting the surface of the protein to the BNC. One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the BNC. Both alternative O₂ channels start in the membrane spanning region and terminate close to Y288_I. These channels are a combination of multiple transiently interconnected hydrophobic cavities, whose opening and closure is regulated by the thermal fluctuations of the lining residues. Furthermore, our results show that, in this Ccox, the most likely (energetically preferred) routes for O₂ to reach the BNC are the alternative channels, rather than the X-ray inferred pathway.

Cytochrome c oxidases (Ccoxs), the terminal enzymes of the respiratory electron transport chain in eukaryotes and many prokaryotes, are key enzymes in aerobic respiration. These proteins couple the reduction of molecular dioxygen to water with the creation of a transmembrane electrochemical proton gradient. Over the last decades, most of the Ccoxs research focused on the mechanisms and energetics of reduction and/or proton pumping, and little emphasis has been given to the pathways used by dioxygen to reach the binuclear center, where dioxygen reduction takes place. In particular, the existence and the characteristics of the channel(s) used by O₂ to travel from the solvent/membrane to the binuclear site are still unclear. In this work, we combine all-atom molecular dynamics simulations and implicit ligand sampling calculations in order to identify and characterize the O₂ delivery channels in the Ccox from Rhodobacter sphaeroides. Altogether, our results suggest that, in this Ccox, O₂ can diffuse via three well-defined channels that start in membrane region (where O₂ solubility is higher than in the water). One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the binuclear center.

The superfamily of heme-copper oxygen reductases: types and evolutionary considerations

Heme-copper oxygen reductases (HCO) reduce O(2) to water being the last enzymatic complexes of most aerobic respiratory chains. These enzymes promote energy conservation coupling the catalytic reaction to charge separation and charge translocation across the prokaryotic cytoplasmatic or mitochondrial membrane. In this way they contribute to the establishment and maintenance of the transmembrane difference of electrochemical potential, which is vital for solute/nutrient cell import, synthesis of ATP and motility. The HCO enzymes most probably share with the nitric oxide reductases, NORs, a common ancestor. We have proposed the classification of HCOs into three different types, A, B and C; based on the constituents of their proton channels (Pereira, Santana and Teixeira (2001) Biochim Biophys Acta, 1505, 185-208). This classification was recently challenged by the suggestion of other different types of HCOs. Using an enlarged sampling we performed an exhaustive bioinformatic reanalysis of HCOs family. Our results strengthened our previously proposed classification and showed no need for the existence of more divisions. Now, we analyze the taxonomic distribution of HCOs and NORs and the congruence of their sequence trees with the 16S rRNA tree. We observed that HCOs are widely distributed in the two prokaryotic domains and that the different types of enzymes are not confined to a specific taxonomic group or environmental niche.

Dioxygen and nitric oxide pathways and affinity to the catalytic site of rubredoxin:oxygen oxidoreductase from Desulfovibrio gigas

Rubredoxin:oxygen oxidoreductase (ROO) is the terminal oxidase of a soluble electron transfer chain found in Desulfovibrio gigas. This protein belongs to the flavodiiron family and was initially described as an oxygen reductase, converting this substrate to water and acting as an oxygen-detoxifying system. However, more recent studies evidenced also the ability for this protein to act as a nitric oxide reductase, suggesting an alternative physiological role. To clarify the apparent bifunctional nature of this protein, we performed molecular dynamics simulations of the protein, in different redox states, together with O(2) and NO molecules in aqueous solution. The two small molecules were parameterized using free-energy calculations of the hydration process. With these simulations we were able to identify specific protein paths that allow the diffusion of both these molecules through the protein towards the catalytic centers. Also, we have tried to characterize the preference of ROO towards the presence of O(2) and/or NO at the active site. By using free-energy simulations, we did not find any significant preference for ROO to accommodate both O(2) and NO. Also, from our molecular dynamics simulations we were able to identify similar diffusion profiles for both O(2) and NO molecules. These two conclusions are in good agreement with previous experimental works stating that ROO is able to catalyze both O(2) and NO.

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.

Diversity of the heme-copper superfamily in archaea: insights from genomics and structural modeling

Recent advances in DNA sequencing technologies have provided unprecedented access into the diversity of the microbial world. Herein we use the comparative genomic analysis of microbial genomes and environmental metagenomes coupled with structural modelling to explore the diversity of aerobic respiration in Archaea. We focus on the heme-copper oxidoreductase superfamily which is responsible for catalyzing the terminal reaction in aerobic respiration-the reduction of molecular oxygen to water. Sequence analyses demonstrate that there are at least eight heme-copper oxygen reductase families: A-, B-, C-, D-, E-, F-, G-, and H-families. Interestingly, five of these oxygen reductase families (D-, E-, F-, G-, and H-families) are currently found exclusively in Archaea. We review the structural properties of all eight families focusing on the members found within Archaea. Structural modelling coupled with sequence analysis suggests that many of the oxygen reductases identified from thermophilic Archaea have modified proton channel properties compared to the currently studied mesophilic bacterial oxygen reductases. These structural differences may be due to adaptation to the specific environments in which these enzymes function. We conclude with a brief analysis of the phylogenetic distribution and evolution of Archaeal heme-copper oxygen reductases.

Pathways of H2 toward the active site of [NiFe]-hydrogenase

Hydrogenases catalyze the reversible oxidation of molecular hydrogen (H(2)), but little is known about the diffusion of H(2) toward the active site. Here we analyze pathways for H(2) permeation using molecular dynamics (MD) simulations in explicit solvent. Various MD simulation replicates were done, to improve the sampling of the system states. H(2) easily permeates hydrogenase in every simulation and it moves preferentially in channels. All H(2) molecules that reach the active site made their approach from the side of the Ni ion. H(2) is able to reach distances of <4 A from the active site, although after 6 A permeation is difficult. In this region we mutated Val-67 into alanine and perform new MD simulations. These simulations show an increase of H(2) inside the protein and at lower distances from the active site. This valine can be a control point in the H(2) access to the active center.

Reorganization and conformational changes in the reduction of tetraheme cytochromes

Molecular dynamics simulation (MD) constitutes an alternative to time-consuming experiments for studying conformational changes. We apply MD on a redox system where experimental information exists for the fully oxidized and fully reduced states: tetraheme cytochrome c3. Instead of doing one simulation for each state, we apply 10 4-ns replicas for both states, which provides robust statistics to characterize the redox changes. Besides these long simulations, we perform 120 short ones (50 ps), where an equilibrated oxidized state is perturbed to a reduced state. This allows the application of a nonequilibrium method, the subtraction technique, which makes it possible to characterize the different timescales of conformational changes. Reduction induces conformational changes in the N-terminus and on the loops spanning residues 36-42 and 88-93, which correlate very well with experiments, demonstrating the applicability of this methodology. We also analyze the effect of reduction on hydrogen bonds, solvent accessible surface and bound water, the changes being found to involve the hemes and propionate groups. Redox-induced protonation is also investigated, by protonating the propionates D from hemes I and IV. Although this change in the former does not have major conformational consequences, it induces in the latter conformational changes beyond the ones obtained with reduction.