The unusual redox properties of C-type oxidases

The structure of the diheme cytochrome c4 from Neisseria gonorrhoeae reveals multiple contributors to tuning reduction potentials

Cytochrome c4 (c4) is a diheme protein implicated as an electron donor to cbb3 oxidases in multiple pathogenic bacteria. Despite its prevalence, understanding of how specific structural features of c4 optimize its function is lacking. The human pathogen Neisseria gonorrhoeae (Ng) thrives in low oxygen environments owing to the activity of its cbb3 oxidase. Herein, we report characterization of Ng c4. Spectroelectrochemistry experiments of the wild-type (WT) protein have shown that the two Met/His-ligated hemes differ in potentials by ∼100 mV, and studies of the two His/His-ligated variants provided unambiguous assignment of heme A from the N-terminal domain of the protein as the high-potential heme. The crystal structure of the WT protein at 2.45 Å resolution has revealed that the two hemes differ in their solvent accessibility. In particular, interactions made by residues His57 and Ser59 in Loop1 near the axial ligand Met63 contribute to the tight enclosure of heme A, working together with the surface charge, to raise the reduction potential of the heme iron in this domain. The structure reveals a prominent positively-charged patch, which encompasses surfaces of both domains. In contrast to prior findings with c4 from Pseudomonas stutzeri, the interdomain interface of Ng c4 contributes minimally to the values of the heme iron potentials in the two domains. Analyses of the heme solvent accessibility, interface properties, and surface charges offer insights into the interplay of these structural elements in tuning redox properties of c4 and other multiheme proteins.

Mutating the environment of heme b595 of E. coli cytochrome bd-I oxidase shifts its redox potential by 200 mV without inactivating the enzyme

Cytochrome bd-I catalyzes the reduction of oxygen to water with the aid of hemes b₅₅₈, b₅₉₅ and d. Here, effects of a mutation of E445, a ligand of heme b₅₉₅ and of R448, hydrogen bonded to E445 are studied electrochemically in the E. coli enzyme. The equilibrium potential of the three hemes are shifted by up to 200 mV in these mutants. Strikingly the E445D and the R448N mutants show a turnover of 41 ± 2 % and 20 ± 4 %, respectively. Electrocatalytic studies confirm that the mutants react with oxygen and bind and release NO. These results point towards the ability of cytochrome bd to react even if the electron transfer is less favorable.

Enhanced Phenazine-1-Carboxamide Production in Pseudomonas chlororaphis H5△fleQ△relA through Fermentation Optimization

Phenazine-1-carboxamide (PCN) is effective to control many plant pathogens, and improving PCN production would be of great significance in promoting its development as a biopesticide. This study was conducted to improve the PCN production of Pseudomonas chlororaphis H5△fleQ△relA through fermentation optimization in both shake flask and bioreactor. The PCN production of H5△fleQ△relA was improved from 2.75 ± 0.23 g/L to 5.51 ± 0.17 g/L by medium optimization in shake flask using Plackett-Burman design, the path of steepest ascent experiment and central composite design. Then, PCN production reached 8.58 ± 0.25 g/L through optimizing pH in 1 L bioreactor. After pH optimization, the transcriptional levels of ccoO_2 and ccoQ_2 genes related to microbial aerobic respiration were significantly upregulated, and the relative abundance of 3-oxo-C14-HSL was significantly enhanced 15-fold, and these changes were vital for cell activity and metabolites production. Furthermore, the PCN production reached 9.58 ± 0.57 g/L after optimization of the fed-batch fermentation strategy in 1 L bioreactor. Finally, the fermentation scale-up of the optimal medium and optimal feeding strategy were conducted in 30 L bioreactor at the optimal pH, and their PCN production reached 9.17 g/L and 9.62 g/L respectively, which were comparable to that in 1 L bioreactor. In this study, the high PCN production was achieved from the shake-flask fermentation to 30 L bioreactor, and the optimal feeding strategy improved PCN production in bioreactor without increasing total glycerol compared with in shake flask. It provides promising pathways for the optimization of processes for the production of other phenazines.

Molecular Details on Multiple Cofactor Containing Redox Metalloproteins Revealed by Infrared and Resonance Raman Spectroscopies

Vibrational spectroscopy and in particular, resonance Raman (RR) spectroscopy, can provide molecular details on metalloproteins containing multiple cofactors, which are often challenging for other spectroscopies. Due to distinct spectroscopic fingerprints, RR spectroscopy has a unique capacity to monitor simultaneously and independently different metal cofactors that can have particular roles in metalloproteins. These include e.g., (i) different types of hemes, for instance hemes c, a and a3 in caa3-type oxygen reductases, (ii) distinct spin populations, such as electron transfer (ET) low-spin (LS) and catalytic high-spin (HS) hemes in nitrite reductases, (iii) different types of Fe-S clusters, such as 3Fe-4S and 4Fe-4S centers in di-cluster ferredoxins, and (iv) bi-metallic center and ET Fe-S clusters in hydrogenases. IR spectroscopy can provide unmatched molecular details on specific enzymes like hydrogenases that possess catalytic centers coordinated by CO and CN- ligands, which exhibit spectrally well separated IR bands. This article reviews the work on metalloproteins for which vibrational spectroscopy has ensured advances in understanding structural and mechanistic properties, including multiple heme-containing proteins, such as nitrite reductases that house a notable total of 28 hemes in a functional unit, respiratory chain complexes, and hydrogenases that carry out the most fundamental functions in cells.

Activation of O2 and NO in heme-copper oxidases - mechanistic insights from computational modelling

Heme-copper oxidases are transmembrane enzymes involved in aerobic and anaerobic respiration. The largest subgroup contains the cytochrome c oxidases (CcO), which reduce molecular oxygen to water. A significant part of the free energy released in this exergonic process is conserved as an electrochemical gradient across the membrane, via two processes, electrogenic chemistry and proton pumping. A deviant subgroup is the cytochrome c dependent NO reductases (cNOR), which reduce nitric oxide to nitrous oxide and water. This is also an exergonic reaction, but in this case none of the released free energy is conserved. Computational studies applying hybrid density functional theory to cluster models of the bimetallic active sites in the heme-copper oxidases are reviewed. To obtain a reliable description of the reaction mechanisms, energy profiles of the entire catalytic cycles, including the reduction steps have to be constructed. This requires a careful combination of computational results with certain experimental data. Computational studies have elucidated mechanistic details of the chemical parts of the reactions, involving cleavage and formation of covalent bonds, which have not been obtainable from pure experimental investigations. Important insights regarding the mechanisms of energy conservation have also been gained. The computational studies show that the reduction potentials of the active site cofactors in the CcOs are large enough to afford electrogenic chemistry and proton pumping, i.e. efficient energy conservation. These results solve a conflict between different types of experimental data. A mechanism for the proton pumping, involving a specific and crucial role for the active site tyrosine, conserved in all CcOs, is suggested. For the cNORs, the calculations show that the low reduction potentials of the active site cofactors are optimized for fast elimination of the toxic NO molecules. At the same time, the low reduction potentials lead to endergonic reduction steps with high barriers. To prevent even higher barriers, which would lead to a too slow reaction, when the electrochemical gradient across the membrane is present, the chemistry must occur in a non-electrogenic manner. This explains why there is no energy conservation in cNOR.

Electrocatalytic evidence of the diversity of the oxygen reaction in the bacterial bd oxidase from different organisms

Cytochrome bd oxidase is a bacterial terminal oxygen reductase that was suggested to enable adaptation to different environments and to confer resistance to stress conditions. An electrocatalytic study of the cyt bd oxidases from Escherichia coli, Corynebacterium glutamicum and Geobacillus thermodenitrificans gives evidence for a different reactivity towards oxygen. An inversion of the redox potential values of the three hemes is found when comparing the enzymes from different bacteria. This inversion can be correlated with different protonated glutamic acids as evidenced by reaction induced FTIR spectroscopy. The influence of the microenvironment of the hemes on the reactivity towards oxygen is discussed.

Spectroscopic Evidence of Hyponitrite Radical Intermediate in NO Disproportionation at a MOF-Supported Mononuclear Copper Site

Dianionic hyponitrite (N₂ O₂ ^2- ) is often proposed, based on model complexes, as the key intermediate in reductive coupling of nitric oxide to nitrous oxide at the bimetallic active sites of heme-copper oxidases and nitric oxide reductases. In this work, we examine the gas-solid reaction of nitric oxide with the metal-organic framework Cu^I -ZrTpmC* with a suite of in situ spectroscopies and density functional theory simulations, and identify an unusual chelating N₂ O₂ ^.- intermediate. These results highlight the advantage provided by site-isolation in metal-organic frameworks (MOFs) for studying important reaction intermediates, and provide a mechanistic scenario compatible with the proposed one-electron couple in these enzymes.

Molecular understanding of heteronuclear active sites in heme-copper oxidases, nitric oxide reductases, and sulfite reductases through biomimetic modelling

Heme-copper oxidases (HCO), nitric oxide reductases (NOR), and sulfite reductases (SiR) catalyze the multi-electron and multi-proton reductions of O2, NO, and SO32-, respectively. Each of these reactions is important to drive cellular energy production through respiratory metabolism and HCO, NOR, and SiR evolved to contain heteronuclear active sites containing heme/copper, heme/nonheme iron, and heme-[4Fe-4S] centers, respectively. The complexity of the structures and reactions of these native enzymes, along with their large sizes and/or membrane associations, make it challenging to fully understand the crucial structural features responsible for the catalytic properties of these active sites. In this review, we summarize progress that has been made to better understand these heteronuclear metalloenzymes at the molecular level though study of the native enzymes along with insights gained from biomimetic models comprising either small molecules or proteins. Further understanding the reaction selectivity of these enzymes is discussed through comparisons of their similar heteronuclear active sites, and we offer outlook for further investigations.

Nitric Oxide Does Not Inhibit but Is Metabolized by the Cytochrome bcc-aa3 Supercomplex

Nitric oxide (NO) is a well-known active site ligand and inhibitor of respiratory terminal oxidases. Here, we investigated the interaction of NO with a purified chimeric bcc-aa₃ supercomplex composed of Mycobacterium tuberculosis cytochrome bcc and Mycobacterium smegmatisaa₃-type terminal oxidase. Strikingly, we found that the enzyme in turnover with O₂ and reductants is resistant to inhibition by the ligand, being able to metabolize NO at 25 °C with an apparent turnover number as high as ≈303 mol NO (mol enzyme)⁻¹ min⁻¹ at 30 µM NO. The rate of NO consumption proved to be proportional to that of O₂ consumption, with 2.65 ± 0.19 molecules of NO being consumed per O₂ molecule by the mycobacterial bcc-aa₃. The enzyme was found to metabolize the ligand even under anaerobic reducing conditions with a turnover number of 2.8 ± 0.5 mol NO (mol enzyme)⁻¹ min⁻¹ at 25 °C and 8.4 µM NO. These results suggest a protective role of mycobacterial bcc-aa₃ supercomplexes against NO stress.

Cytochrome c Reductase is a Key Enzyme Involved in the Extracellular Electron Transfer Pathway towards Transition Metal Complexes in Pseudomonas Putida

Mediator-based extracellular electron transfer (EET) pathways can balance the redox metabolism of microbes. However, such electro-biosynthesis processes are constrained by the unknown underlying EET mechanisms. In this paper, Pseudomonas putida was studied to systematically investigate its EET pathway to transition metal complexes (i. e., [Fe(CN)6 ]3-/4- and [Co(bpy)3 ]3+/2+ ; bpy=2,2'-bipyridyl) under anaerobic conditions. Comparative proteomics showed the aerobic respiratory components were upregulated in a bioelectrochemical system without oxygen, suggesting their potential contribution to EET. Further tests found inhibiting cytochrome c oxidase activity by NaN3 and NADH dehydrogenase by rotenone did not significantly change the current output. However, the EET pathway was completely blocked, while cytochrome c reductase activity was inhibited by antimycin A. Although it cannot be excluded that cytochrome c and the periplasmic subunit of cytochrome c oxidase donate electrons to the transition metal complexes, these results strongly demonstrate that cytochrome c reductase is a key complex for the EET pathway.

Role of the Two Metals in the Active Sites of Heme Copper Oxidases-A Study of NO Reduction in cbb3 Cytochrome c Oxidase

The superfamily of heme copper oxidases reduces molecular oxygen or nitric oxide, and the active sites comprise a high-spin heme group (a₃ or b₃) and a non-heme metal (Cu_B or Fe_B). The cbb₃ C family of cytochrome c oxidases, with the high-spin heme b₃ and Cu_B in the active site, is a subfamily of the heme copper oxidases that can reduce both molecular oxygen, which is the main substrate, and nitric oxide. The mechanism for NO reduction in cbb₃ oxidase is studied here using hybrid density functional theory and compared to other cytochrome c oxidases (A and B families), with a high-spin heme a₃ and Cu_B in the active site, and to cytochrome c dependent NO reductase, with a high-spin heme b₃ and a non-heme Fe_B in the active site. It is found that the reaction mechanism and the detailed reaction energetics of the cbb₃ oxidases are not similar to those of cytochrome c dependent NO reductase, which has the same type of high-spin heme group but a different non-heme metal. This is in contrast to earlier expectations. Instead, the NO reduction mechanism in cbb₃ oxidases is very similar to that in the other cytochrome c oxidases, with the same non-heme metal, Cu_B, and is independent of the type of high-spin heme group. The conclusion is that the type of non-heme metal (Cu_B or Fe_B) in the active site of the heme copper oxidases is more important for the reaction mechanisms than the type of high-spin heme, at least for the NO reduction reaction. The reason is that the proton-coupled reduction potentials of the active site cofactors determine the energetics for the NO reduction reaction, and they depend to a larger extent on the non-heme metal. Observed differences in NO reduction reactivity among the various cytochrome c oxidases may be explained by differences outside the BNC, affecting the rate of proton transfer, rather than in the BNC itself.

In situ absorbance measurements: a new means to study respiratory electron transfer in chemolithotrophic microorganisms

Absorbance measurements on intact chemolithotrophic microorganisms that respire aerobically on soluble iron are described that used a novel integrating cavity absorption meter to eliminate the effects of light scattering on the experimental results. Steady state kinetic measurements on ferric iron production by intact cells revealed that the Michaelis Menten equation described the initial rates of product formation for at least 8 different chemolithotrophic microorganisms in 6 phyla distributed equally among the archaea and the Gram negative and Gram positive eubacteria. Cell-monitored turnover measurements during aerobic respiration on soluble iron by the same 12 intact microorganisms revealed six different patterns of iron-dependent absorbance changes, suggesting that there may be at least six different sets of prosthetic groups and biomolecules that can accomplish aerobic respiration on soluble iron. Detailed kinetic studies revealed that the 3-component iron respiratory chain of Acidithiobacillus ferrooxidans functioned as an ensemble with a single macroscopic rate constant when the iron-reduced proteins were oxidized in the presence of excess molecular oxygen. The principal member of this 3-component system was a cupredoxin called rusticyanin that was present in the periplasm of At. ferrooxidans at an approximate concentration of 350 mg/mL, an observation that provides new insights into the crowded environments in the periplasms of Gram negative eubacteria that conduct electrons across their periplasm. The ability to conduct direct spectrophotometric measurements under noninvasive physiological conditions represents a new and powerful approach to examine the rates and extents of biological events in situ without disrupting the complexity of the live cellular environment.

The mechanism for oxygen reduction in the C family cbb3 cytochrome c oxidases - Implications for the proton pumping stoichiometry

Cytochrome c oxidases (CcOs) couple the exergonic reduction of molecular oxygen to proton pumping across the membrane in which they are embedded, thereby conserving a significant part of the free energy. The A family CcOs are known to pump four protons per oxygen molecule, while there is no consensus regarding the proton pumping stoichiometry for the C family cbb₃ oxidases. Hybrid density functional theory is used here to investigate the catalytic mechanism for oxygen reduction in cbb₃ oxidases. A surprising result is that the barrier for O O bond cleavage at the mixed valence reduction level seems to be too high compared to the overall reaction rate of the enzyme. It is therefore suggested that the O O bond is cleaved only after the first proton coupled reduction step, and that this reduction step most likely is not coupled to proton pumping. Furthermore, since the cbb₃ oxidases have only one proton channel leading to the active site, it is proposed that the activated E_H intermediate, suggested to be responsible for proton pumping in one of the reduction steps in the A family, cannot be involved in the catalytic cycle for cbb₃, which results in the lack of proton pumping also in the E to R reduction step. In summary, the calculations indicate that only two protons are pumped per oxygen molecule in cbb₃ oxidases. However, more experimental information on this divergent enzyme is needed, e.g. whether the flow of electrons resembles that in the other more well-studied CcO families.

Bacterial Respiratory Chain Diversity Reveals a Cytochrome c Oxidase Reducing O2 at Low Overpotentials

Cytochrome c oxidases (CcOs) are the terminal enzymes in energy-converting chains of microorganisms, where they reduce oxygen into water. Their affinity for O₂ makes them attractive biocatalysts for technological devices in which O₂ concentration is limited, but the high overpotentials they display on electrodes severely limit their applicative use. Here, the CcO of the acidophilic bacterium Acidithiobacillus ferrooxidans is studied on various carbon materials by direct protein electrochemistry and mediated one with redox mediators either diffusing or co-immobilized at the electrode surface. The entrapment of the CcO in a network of hydrophobic carbon nanofibers permits a direct electrochemical communication between the enzyme and the electrode. We demonstrate that the CcO displays a μM affinity for O₂ and reduces O₂ at exceptionally high electrode potentials in the range of +700 to +540 mV vs NHE over a pH range of 4-6. The kinetics of interactions between the enzyme and its physiological partners are fully quantified. Based on these results, an electron transfer pathway allowing O₂ reduction in the acidic metabolic chain is proposed.

A low-potential terminal oxidase associated with the iron-only nitrogenase from the nitrogen-fixing bacterium Azotobacter vinelandii

The biological route for nitrogen gas entering the biosphere is reduction to ammonia by the nitrogenase enzyme, which is inactivated by oxygen. Three types of nitrogenase exist, the least-studied of which is the iron-only nitrogenase. The Anf3 protein in the bacterium Rhodobacter capsulatus is essential for diazotrophic (i.e. nitrogen-fixing) growth with the iron-only nitrogenase, but its enzymatic activity and function are unknown. Here, we biochemically and structurally characterize Anf3 from the model diazotrophic bacterium Azotobacter vinelandii Determining the Anf3 crystal structure to atomic resolution, we observed that it is a dimeric flavocytochrome with an unusually close interaction between the heme and the FAD cofactors. Measuring the reduction potentials by spectroelectrochemical redox titration, we observed values of -420 ± 10 and -330 ± 10 mV for the two FAD potentials and -340 ± 1 mV for the heme. We further show that Anf3 accepts electrons from spinach ferredoxin and that Anf3 consumes oxygen without generating superoxide or hydrogen peroxide. We predict that Anf3 protects the iron-only nitrogenase from oxygen inactivation by functioning as an oxidase in respiratory protection, with flavodoxin or ferredoxin as the physiological electron donors.

Active Site Midpoint Potentials in Different Cytochrome c Oxidase Families: A Computational Comparison

Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain, reducing molecular oxygen to water. The binuclear active site in C cO comprises a high-spin heme associated with a Cu_B complex and a redox active tyrosine. The electron transport in the respiratory chain is driven by increasing midpoint potentials of the involved cofactors, resulting in a release of free energy, which is stored by coupling the electron transfer to proton translocation across a membrane, building up an electrochemical gradient. In this context, the midpoint potentials of the active site cofactors in the C cOs are of special interest, since they determine the driving forces for the individual oxygen reduction steps and thereby affect the efficiency of the proton pumping. It has been difficult to obtain useful information on some of these midpoint potentials from experiments. However, since each of the reduction steps in the catalytic cycle of oxygen reduction to water corresponds to the formation of an O-H bond, they can be calculated with a reasonably high accuracy using quantum chemical methods. From the calculated O-H bond strengths, the proton-coupled midpoint potentials of the active site cofactors can be estimated. Using models representing the different families of C cO's (A, B, and C), the calculations give midpoint potentials that should be relevant during catalytic turnover. The calculations also suggest possible explanations for why some experimentally measured potentials deviate significantly from the calculated ones, i.e., for Cu_B in all oxidase families, and for heme b₃ in the C family.

The regulation of oxidative phosphorylation pathway on Vibrio alginolyticus adhesion under adversities

Vibrio alginolyticus is one of the most important pathogens in mariculture and leading to heavy losses. After treatment with Cu²⁺, Pb²⁺, and low pH, the expression of oxidative phosphorylation pathway genes, including coxA, coxB, coxC, ccoN, ccoO, and ccoQ, was found commonly downregulated by RNA‐seq as well as quantitative real‐time PCR. RNAi significantly reduced the expression of coxA, coxB, coxC, ccoN, ccoO, and ccoQ in V. alginolyticus. Compared with the wild‐type strain, the adhesion abilities of RNAi strains of V. alginolyticus were significantly impaired, as well as their cytochrome C oxidase activity. ccoQ appeared to be more important in the regulation of bacterial adhesion in these target genes, while ccoO was relatively weak in the regulation of the adhesion. Meanwhile, the changes of temperature, salinity, pH, and starvation affected coxA, coxB, coxC, ccoN, ccoO, and ccoQ expression remarkably. These findings indicated that: the oxidative phosphorylation pathway is a critical regulator of adhesion in V. alginolyticus; coxA, coxB, coxC, ccoN, ccoO, and ccoQ regulate the bacterial adhesion in response to environmental changes such as temperature, salinity, pH, and starvation.

Heme redox potentials hold the key to reactivity differences between nitric oxide reductase and heme-copper oxidase

Despite high structural homology between NO reductases (NORs) and heme-copper oxidases (HCOs), factors governing their reaction specificity remain to be understood. Using a myoglobin-based model of NOR (Fe_BMb) and tuning its heme redox potentials (E°') to cover the native NOR range, through manipulating hydrogen bonding to the proximal histidine ligand and replacing heme b with monoformyl (MF-) or diformyl (DF-) hemes, we herein demonstrate that the E°' holds the key to reactivity differences between NOR and HCO. Detailed electrochemical, kinetic, and vibrational spectroscopic studies, in tandem with density functional theory calculations, demonstrate a strong influence of heme E°' on NO reduction. Decreasing E°' from +148 to -130 mV significantly impacts electronic properties of the NOR mimics, resulting in 180- and 633-fold enhancements in NO association and heme-nitrosyl decay rates, respectively. Our results indicate that NORs exhibit finely tuned E°' that maximizes their enzymatic efficiency and helps achieve a balance between opposite factors: fast NO binding and decay of dinitrosyl species facilitated by low E°' and fast electron transfer facilitated by high E°'. Only when E°' is optimally tuned in Fe_BMb(MF-heme) for NO binding, heme-nitrosyl decay, and electron transfer does the protein achieve multiple (>35) turnovers, previously not achieved by synthetic or enzyme-based NOR models. This also explains a long-standing question in bioenergetics of selective cross-reactivity in HCOs. Only HCOs with heme E°' in a similar range as NORs (between -59 and 200 mV) exhibit NOR reactivity. Thus, our work demonstrates efficient tuning of E°' in various metalloproteins for their optimal functionality.

The electron distribution in the "activated" state of cytochrome c oxidase

Cytochrome c oxidase catalyzes reduction of O₂ to H₂O at a catalytic site that is composed of a copper ion and heme group. The reaction is linked to translocation of four protons across the membrane for each O₂ reduced to water. The free energy associated with electron transfer to the catalytic site is unequal for the four electron-transfer events. Most notably, the free energy associated with reduction of the catalytic site in the oxidized cytochrome c oxidase (state O) is not sufficient for proton pumping across the energized membrane. Yet, this electron transfer is mechanistically linked to proton pumping. To resolve this apparent discrepancy, a high-energy oxidized state (denoted O_H) was postulated and suggested to be populated only during catalytic turnover. The difference between states O and O_H was suggested to be manifested in an elevated midpoint potential of Cu_B in the latter. This proposal predicts that one-electron reduction of cytochrome c oxidase after its oxidation would yield re-reduction of essentially only Cu_B. Here, we investigated this process and found ~5% and ~6% reduction of heme a₃ and Cu_B, respectively, i.e. the apparent redox potentials for heme a₃ and Cu_B are lower than that of heme a.

On the role of subunit M in cytochrome cbb3 oxidase

Cytochrome cbb₃ (or C-type) oxidases are a highly divergent group and the least studied members of the heme-copper oxidases (HCOs) superfamily. HCOs couple the reduction of oxygen at the end of the respiratory chain to the active proton translocation across the membrane, contributing to establishment of an electrochemical gradient essential for ATP synthesis. Cbb₃ oxidases exhibit unique structural and functional features and have an essential role in the metabolism of many clinically relevant human pathogens. Such characteristics make them a promising therapeutic target. Three subunits, N, O and P, comprise the core cbb₃ complex, with N, the catalytic subunit, being highly conserved among all members of the HCO superfamily, including the A-type (aa₃, mitochondrial-like) oxidases. An additional fourth subunit containing a single transmembrane (TM) helix was present in the first crystal structure of cbb₃. This TM segment was recently proposed to be part of a novel protein CcoM, which was shown to have a putative role in the complex stability and assembly. In this work, we performed large-scale all-atom molecular dynamics simulations of the CcoNOPM complex to further characterize the interactions between subunit M and the core subunits and to determine whether the presence of the fourth subunit influences the water/proton channels previously described for the core complex. The previously proposed putative CcoNOPH complex is also assessed, and the potential functional redundancy of CcoM and CcoQ is discussed.