Crystallographic studies of Xe and Kr binding within the large internal cavity of cytochrome ba3 from Thermus thermophilus: structural analysis and role of oxygen transport channels in the heme-Cu oxidases

Investigation of the Mechanism of Membrane Potential Generation by Heme-Copper Respiratory Oxidases in a Real Time Mode

Heme-copper respiratory oxidases are highly efficient molecular machines. These membrane enzymes catalyze the final step of cellular respiration in eukaryotes and many prokaryotes: the transfer of electrons from cytochromes or quinols to molecular oxygen and oxygen reduction to water. The free energy released in this redox reaction is converted by heme-copper respiratory oxidases into the transmembrane gradient of the electrochemical potential of hydrogen ions H+). Heme-copper respiratory oxidases have a unique mechanism for generating H+, namely, a redox-coupled proton pump. A combination of direct electrometric method for measuring the kinetics of membrane potential generation with the methods of prestationary kinetics and site-directed mutagenesis in the studies of heme-copper oxidases allows to obtain a unique information on the translocation of protons inside the proteins in real time. The review summarizes the data of studies employing time-resolved electrometry to decipher the mechanisms of functioning of these important bioenergetic enzymes.

Structural basis for safe and efficient energy conversion in a respiratory supercomplex

Proton-translocating respiratory complexes assemble into supercomplexes that are proposed to increase the efficiency of energy conversion and limit the production of harmful reactive oxygen species during aerobic cellular respiration. Cytochrome bc complexes and cytochrome aa₃ oxidases are major drivers of the proton motive force that fuels ATP generation via respiration, but how wasteful electron- and proton transfer is controlled to enhance safety and efficiency in the context of supercomplexes is not known. Here, we address this question with the 2.8 Å resolution cryo-EM structure of the cytochrome bcc-aa₃ (III₂-IV₂) supercomplex from the actinobacterium Corynebacterium glutamicum. Menaquinone, substrate mimics, lycopene, an unexpected Q_c site, dioxygen, proton transfer routes, and conformational states of key protonable residues are resolved. Our results show how safe and efficient energy conversion is achieved in a respiratory supercomplex through controlled electron and proton transfer. The structure may guide the rational design of drugs against actinobacteria that cause diphtheria and tuberculosis.

Aerobic energy metabolism is driven by proton-pumping respiratory supercomplexes. The study reports the structural basis for energy conversion in such supercomplex. It may aid metabolic engineering and drug design against diphtheria and tuberculosis.

Comparative study of the effects of high hydrostatic pressure per se and high argon pressure on urate oxidase ligand stabilization

The stability of the tetrameric enzyme urate oxidase in complex with excess of 8-azaxanthine was investigated either under high hydrostatic pressure per se or under a high pressure of argon. The active site is located at the interface of two subunits, and the catalytic activity is directly related to the integrity of the tetramer. This study demonstrates that applying pressure to a protein–ligand complex drives the thermodynamic equilibrium towards ligand saturation of the complex, revealing a new binding site. A transient dimeric intermediate that occurs during the pressure-induced dissociation process was characterized under argon pressure and excited substates of the enzyme that occur during the catalytic cycle can be trapped by pressure. Comparison of the different structures under pressure infers an allosteric role of the internal hydrophobic cavity in which argon is bound, since this cavity provides the necessary flexibility for the active site to function.

Cryo-EM structures of intermediates suggest an alternative catalytic reaction cycle for cytochrome c oxidase

Cytochrome c oxidases are among the most important and fundamental enzymes of life. Integrated into membranes they use four electrons from cytochrome c molecules to reduce molecular oxygen (dioxygen) to water. Their catalytic cycle has been considered to start with the oxidized form. Subsequent electron transfers lead to the E-state, the R-state (which binds oxygen), the P-state (with an already split dioxygen bond), the F-state and the O-state again. Here, we determined structures of up to 1.9 Å resolution of these intermediates by single particle cryo-EM. Our results suggest that in the O-state the active site contains a peroxide dianion and in the P-state possibly an intact dioxygen molecule, the F-state may contain a superoxide anion. Thus, the enzyme's catalytic cycle may have to be turned by 180 degrees.

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.

Coupled transport of electrons and protons in a bacterial cytochrome c oxidase-DFT calculated properties compared to structures and spectroscopies

After a general introduction to the features and mechanisms of cytochrome c oxidases (CcOs) in mitochondria and aerobic bacteria, we present DFT calculated physical and spectroscopic properties for the catalytic reaction cycle compared with experimental observations in bacterial ba3 type CcO, also with comparisons/contrasts to aa3 type CcOs. The Dinuclear Complex (DNC) is the active catalytic reaction center, containing a heme a3 Fe center and a near lying Cu center (called CuB) where by successive reduction and protonation, molecular O2 is transformed to two H2O molecules, and protons are pumped from an inner region across the membrane to an outer region by transit through the CcO integral membrane protein. Structures, energies and vibrational frequencies for Fe-O and O-O modes are calculated by DFT over the catalytic cycle. The calculated DFT frequencies in the DNC of CcO are compared with measured frequencies from Resonance Raman spectroscopy to clarify the composition, geometry, and electronic structures of different intermediates through the reaction cycle, and to trace reaction pathways. X-ray structures of the resting oxidized state are analyzed with reference to the known experimental reaction chemistry and using DFT calculated structures in fitting observed electron density maps. Our calculations lead to a new proposed reaction pathway for coupling the PR → F → OH (ferryl-oxo → ferric-hydroxo) pathway to proton pumping by a water shift mechanism. Through this arc of the catalytic cycle, major shifts in pKa's of the special tyrosine and a histidine near the upper water pool activate proton transfer. Additional mechanisms for proton pumping are explored, and the role of the CuB+ (cuprous state) in controlling access to the dinuclear reaction site is proposed.

Paramagnetic Shifts and Guest Exchange Kinetics in ConFe4-n Metal-Organic Capsules

We investigate the magnetic resonance properties and exchange kinetics of guest molecules in a series of hetero-bimetallic capsules, [Co_nFe_4-nL₆]^4- (n = 1-3), where L^2- = 4,4'-bis[(2-pyridinylmethylene)amino]-[1,1'-biphenyl]-2,2'-disulfonate. H bond networks between capsule sulfonates and guanidinium cations promote the crystallization of [Co_nFe_4-nL₆]^4-. The following four isostructural crystals are reported: two guest-free forms, (C(NH₂)₃)₄[Co_1.8Fe_2.2L₆]·69H₂O (1) and (C(NH₂)₃)₄[Co_2.7Fe_1.3L₆]·73H₂O (2), and two Xe- and CFCl₃-encapsulated forms, (C(NH₂)₃)₄[(Xe)_0.8Co_1.8Fe_2.2L₆]·69H₂O (3) and (C(NH₂)₃)₄[(CFCl₃)Co_2.0Fe_2.0L₆]·73H₂O (4), respectively. Structural analyses reveal that Xe induces negligible structural changes in 3, while the angles between neighboring phenyl groups expand by ca. 3° to accommodate the much larger guest, CFCl₃, in 4. These guest-encapsulated [Co_nFe_4-nL₆]^4- molecules reveal ¹²⁹Xe and ¹⁹F chemical shift changes of ca. -22 and -10 ppm at 298 K, respectively, per substitution of low-spin Fe^II by high-spin Co^II. Likewise, the temperature dependence of the ¹²⁹Xe and ¹⁹F NMR resonances increases by 0.1 and 0.06 ppm/K, respectively, with each additional paramagnetic Co^II center. The optimal temperature for hyperpolarized (hp) ¹²⁹Xe chemical exchange saturation transfer (hyper-CEST) with [Co_nFe_4-nL₆]^4- capsules was found to be inversely proportional to the number of Co^II centers, n, which is consistent with the Xe chemical exchange accelerating as the portals expand. The systematic study was facilitated by the tunability of the [M₄L₆]^4- capsules, further highlighting these metal-organic systems for developing responsive sensors with highly shifted ¹²⁹Xe resonances.

Krypton-derivatization highlights O2-channeling in a four-electron reducing oxidase

Kr-derivatization and X-ray structures indicated O₂-channel and gating-loop that prevent side-reaction in reduction of O₂ to water in F₄₂₀H₂ oxidase.

Dimeric structures of quinol-dependent nitric oxide reductases (qNORs) revealed by cryo-electron microscopy

Quinol-dependent nitric oxide reductases (qNORs) are membrane-integrated, iron-containing enzymes of the denitrification pathway, which catalyze the reduction of nitric oxide (NO) to the major ozone destroying gas nitrous oxide (N₂O). Cryo-electron microscopy structures of active qNOR from Alcaligenes xylosoxidans and an activity-enhancing mutant have been determined to be at local resolutions of 3.7 and 3.2 Å, respectively. They unexpectedly reveal a dimeric conformation (also confirmed for qNOR from Neisseria meningitidis) and define the active-site configuration, with a clear water channel from the cytoplasm. Structure-based mutagenesis has identified key residues involved in proton transport and substrate delivery to the active site of qNORs. The proton supply direction differs from cytochrome c-dependent NOR (cNOR), where water molecules from the cytoplasm serve as a proton source similar to those from cytochrome c oxidase.

Investigations of Accessibility of T2/T3 Copper Center of Two-Domain Laccase from Streptomyces griseoflavus Ac-993

Laccases (EC 1.10.3.2) are multicopper oxidoreductases acting on diphenols and related substances. Laccases are highly important for biotechnology and environmental remediation. These enzymes contain mononuclear one T2 copper ion and two T3 copper ions (Cu3_α and Cu3_β), which form the so-called trinuclear center (TNC). Along with the typical three-domain laccases Bacteria produce two-domain (2D) enzymes, which are active at neutral and basic pH, thermostable, and resistant to inhibitors. In this work we present the comparative analysis of crystal structures and catalytic properties of recombinant 2D laccase from Streptomyces griseoflavus Ac-993 (SgfSL) and its four mutant forms with replacements of two amino acid residues, located at the narrowing of the presumable T3-solvent tunnels. We obtained inactive enzymes with substitutions of His165, with Phe, and Ile170 with Ala or Phe. His165Ala variant was more active than the wild type. We suggest that His165 is a “gateway” at the O₂-tunnel leading from solvent to the Cu3_β of the enzyme. The side chain of Ile170 could be indirectly involved in the coordination of copper ions at the T3 center by maintaining the position of the imidazole ring of His157 that belongs to the first coordination sphere of Cu3_α.

Oxygen Activation Switch in the Copper Amine Oxidase of Escherichia coli

Copper amine oxidases (CuAOs) are metalloenzymes that reduce molecular oxygen to hydrogen peroxide during catalytic turnover of primary amines. In addition to Cu²⁺ in the active site, two peripheral calcium sites, ∼32 Å from the active site, have roles in Escherichia coli amine oxidase (ECAO). The buried Ca²⁺ (Asp533, Leu534, Asp535, Asp678, and Ala679) is essential for full-length protein production, while the surface Ca²⁺ (Glu573, Tyr667, Asp670, and Glu672) modulates biogenesis of the 2,4,5-trihydroxyphenylalanine quinone (TPQ) cofactor. The E573Q mutation at the surface site prevents calcium binding and TPQ biogenesis. However, TPQ biogenesis can be restored by a suppressor mutation (I342F) in the proposed oxygen delivery channel to the active site. While supporting TPQ biogenesis (∼60% WTECAO TPQ), I342F/E573Q has almost no amine oxidase activity (∼4.6% WTECAO activity). To understand how these long-range mutations have major effects on TPQ biogenesis and catalysis, we employed ultraviolet-visible spectroscopy, steady-state kinetics, inhibition assays, and X-ray crystallography. We show that the surface metal site controls the equilibrium (disproportionation) of the Cu²⁺-substrate reduced TPQ (TPQ_AMQ) Cu⁺-TPQ semiquinone (TPQ_SQ) couple. Removal of the calcium ion from this site by chelation or mutagenesis shifts the equilibrium to Cu²⁺-TPQ_AMQ or destabilizes Cu⁺-TPQ_SQ. Crystal structure analysis shows that TPQ biogenesis is stalled at deprotonation in the Cu²⁺-tyrosinate state. Our findings support WTECAO using the inner sphere electron transfer mechanism for oxygen reduction during catalysis, and while a Cu⁺-tyrosyl radical intermediate is not essential for TPQ biogenesis, it is required for efficient biogenesis.

Bacterial denitrifying nitric oxide reductases and aerobic respiratory terminal oxidases use similar delivery pathways for their molecular substrates

The superfamily of heme‑copper oxidoreductases (HCOs) include both NO and O₂ reductases. Nitric oxide reductases (NORs) are bacterial membrane enzymes that catalyze an intermediate step of denitrification by reducing nitric oxide (NO) to nitrous oxide (N₂O). They are structurally similar to heme‑copper oxygen reductases (HCOs), which reduce O₂ to water. The experimentally observed apparent bimolecular rate constant of NO delivery to the deeply buried catalytic site of NORs was previously reported to approach the diffusion-controlled limit (10⁸-10⁹ M^-1 s^-1). Using the crystal structure of cytochrome-c dependent NOR (cNOR) from Pseudomonas aeruginosa, we employed several protocols of molecular dynamics (MD) simulation, which include flooding simulations of NO molecules, implicit ligand sampling and umbrella sampling simulations, to elucidate how NO in solution accesses the catalytic site of this cNOR. The results show that NO partitions into the membrane, enters the enzyme from the lipid bilayer and diffuses to the catalytic site via a hydrophobic tunnel that is resolved in the crystal structures. This is similar to what has been found for O₂ diffusion through the closely related O₂ reductases. The apparent second order rate constant approximated using the simulation data is ~5 × 10⁸ M^-1 s^-1, which is optimized by the dynamics of the amino acid side chains lining in the tunnel. It is concluded that both NO and O₂ reductases utilize well defined hydrophobic tunnels to assure that substrate diffusion to the buried catalytic sites is not rate limiting under physiological conditions.

Xenon-inhibition of the MscL mechano-sensitive channel and the CopB copper ATPase under different conditions suggests direct effects on these proteins

Xenon is frequently used as a general anesthetic in humans, but the mechanism remains an issue of debate. While for some membrane proteins, a direct interaction of xenon with the protein has been shown to be the inhibitory mechanism, other membrane protein functions could be affected by changes of membrane properties due to partitioning of the gas into the lipid bilayer. Here, the effect of xenon on a mechanosensitive ion channel and a copper ion-translocating ATPase was compared under different conditions. Xenon inhibited spontaneous gating of the Escherichia coli mechano-sensitive mutant channel MscL-G22E, as shown by patch-clamp recording techniques. Under high hydrostatic pressure, MscL-inhibition was reversed. Similarly, the activity of the Enterococcus hirae CopB copper ATPase, reconstituted into proteoliposomes, was inhibited by xenon. However, the CopB ATPase activity was also inhibited by xenon when CopB was in a solubilized state. These findings suggest that xenon acts by directly interacting with these proteins, rather than via indirect effects by altering membrane properties. Also, inhibition of copper transport may be a novel effect of xenon that contributes to anesthesia.

Cytochrome aa3 Oxygen Reductase Utilizes the Tunnel Observed in the Crystal Structures To Deliver O2 for Catalysis

Cytochrome aa₃ is the terminal respiratory enzyme of all eukaryotes and many bacteria and archaea, reducing O₂ to water and harnessing the free energy from the reaction to generate the transmembrane electrochemical potential. The diffusion of O₂ to the heme-copper catalytic site, which is buried deep inside the enzyme, is the initiation step of the reaction chemistry. Our previous molecular dynamics (MD) study with cytochrome ba₃, a homologous enzyme of cytochrome aa₃ in Thermus thermophilus, demonstrated that O₂ diffuses from the lipid bilayer to its reduction site through a 25 Å long tunnel inferred by Xe binding sites detected by X-ray crystallography [Mahinthichaichan, P., Gennis, R., and Tajkhorshid, E. (2016) Biochemistry 55, 1265-1278]. Although a similar tunnel is observed in cytochrome aa₃, this putative pathway appears partially occluded between the entrances and the reduction site. Also, the experimentally determined second-order rate constant for O₂ delivery in cytochrome aa₃ (∼10⁸ M^-1 s^-1) is 10 times slower than that in cytochrome ba₃ (∼10⁹ M^-1 s^-1). A question to be addressed is whether cytochrome aa₃ utilizes this X-ray-inferred tunnel as the primary pathway for O₂ delivery. Using complementary computational methods, including multiple independent flooding MD simulations and implicit ligand sampling calculations, we probe the O₂ delivery pathways in cytochrome aa₃ of Rhodobacter sphaeroides. All of the O₂ molecules that arrived in the reduction site during the simulations were found to diffuse through the X-ray-observed tunnel, despite its apparent constriction, supporting its role as the main O₂ delivery pathway in cytochrome aa₃. The rate constant for O₂ delivery in cytochrome aa₃, approximated using the simulation results, is 10 times slower than in cytochrome ba₃, in agreement with the experimentally determined rate constants.

Oxygen Activation and Energy Conservation by Cytochrome c Oxidase

This review focuses on the type A cytochrome c oxidases (CcO), which are found in all mitochondria and also in several aerobic bacteria. CcO catalyzes the respiratory reduction of dioxygen (O₂) to water by an intriguing mechanism, the details of which are fairly well understood today as a result of research for over four decades. Perhaps even more intriguingly, the membrane-bound CcO couples the O₂ reduction chemistry to translocation of protons across the membrane, thus contributing to generation of the electrochemical proton gradient that is used to drive the synthesis of ATP as catalyzed by the rotary ATP synthase in the same membrane. After reviewing the structure of the core subunits of CcO, the active site, and the transfer paths of electrons, protons, oxygen, and water, we describe the states of the catalytic cycle and point out the few remaining uncertainties. Finally, we discuss the mechanism of proton translocation and the controversies in that area that still prevail.

Inert Gas Deactivates Protein Activity by Aggregation

Biologically inert gases play important roles in the biological functionality of proteins. However, researchers lack a full understanding of the effects of these gases since they are very chemically stable only weakly absorbed by biological tissues. By combining X-ray fluorescence, particle sizing and molecular dynamics (MD) simulations, this work shows that the aggregation of these inert gases near the hydrophobic active cavity of pepsin should lead to protein deactivation. Micro X-ray fluorescence spectra show that a pepsin solution can contain a high concentration of Xe or Kr after gassing, and that the gas concentrations decrease quickly with degassing time. Biological activity experiments indicate a reversible deactivation of the protein during this gassing and degassing. Meanwhile, the nanoparticle size measurements reveal a higher number of “nanoparticles” in gas-containing pepsin solution, also supporting the possible interaction between inert gases and the protein. Further, MD simulations indicate that gas molecules can aggregate into a tiny bubble shape near the hydrophobic active cavity of pepsin, suggesting a mechanism for reducing their biological function.

Dynamics of nitric oxide controlled by protein complex in bacterial system

Nitric oxide (NO) plays diverse and significant roles in biological processes despite its cytotoxicity, raising the question of how biological systems control the action of NO to minimize its cytotoxicity in cells. As a great example of such a system, we found a possibility that NO-generating nitrite reductase (NiR) forms a complex with NO-decomposing membrane-integrated NO reductase (NOR) to efficiently capture NO immediately after its production by NiR in anaerobic nitrate respiration called denitrification. The 3.2-Å resolution structure of the complex of one NiR functional homodimer and two NOR molecules provides an idea of how these enzymes interact in cells, while the structure may not reflect the one in cells due to the membrane topology. Subsequent all-atom molecular dynamics (MD) simulations of the enzyme complex model in a membrane and structure-guided mutagenesis suggested that a few interenzyme salt bridges and coulombic interactions of NiR with the membrane could stabilize the complex of one NiR homodimer and one NOR molecule and contribute to rapid NO decomposition in cells. The MD trajectories of the NO diffusion in the NiR:NOR complex with the membrane showed that, as a plausible NO transfer mechanism, NO released from NiR rapidly migrates into the membrane, then binds to NOR. These results help us understand the mechanism of the cellular control of the action of cytotoxic NO.

Splitting of the O-O bond at the heme-copper catalytic site of respiratory oxidases

Heme-copper oxidases catalyze the four-electron reduction of O₂ to H₂O at a catalytic site that is composed of a heme group, a copper ion (Cu_B), and a tyrosine residue. Results from earlier experimental studies have shown that the O-O bond is cleaved simultaneously with electron transfer from a low-spin heme (heme a/b), forming a ferryl state (P_R ; Fe⁴⁺=O^2-, Cu_B²⁺-OH^-). We show that with the Thermus thermophilus ba₃ oxidase, at low temperature (10°C, pH 7), electron transfer from the low-spin heme b to the catalytic site is faster by a factor of ~10 (τ ≅ 11 μs) than the formation of the P_R ferryl (τ ≅110 μs), which indicates that O₂ is reduced before the splitting of the O-O bond. Application of density functional theory indicates that the electron acceptor at the catalytic site is a high-energy peroxy state [Fe³⁺-O^--O^-(H⁺)], which is formed before the P_R ferryl. The rates of heme b oxidation and P_R ferryl formation were more similar at pH 10, indicating that the formation of the high-energy peroxy state involves proton transfer within the catalytic site, consistent with theory. The combined experimental and theoretical data suggest a general mechanism for O₂ reduction by heme-copper oxidases.

Migration of small ligands in globins: Xe diffusion in truncated hemoglobin N

In heme proteins, the efficient transport of ligands such as NO or O₂ to the binding site is achieved via ligand migration networks. A quantitative assessment of ligand diffusion in these networks is thus essential for a better understanding of the function of these proteins. For this, Xe migration in truncated hemoglobin N (trHbN) of Mycobacterium Tuberculosis was studied using molecular dynamics simulations. Transitions between pockets of the migration network and intra-pocket relaxation occur on similar time scales (10 ps and 20 ps), consistent with low free energy barriers (1-2 kcal/mol). Depending on the pocket from where Xe enters a particular transition, the conformation of the side chains lining the transition region differs which highlights the coupling between ligand and protein degrees of freedom. Furthermore, comparison of transition probabilities shows that Xe migration in trHbN is a non-Markovian process. Memory effects arise due to protein rearrangements and coupled dynamics as Xe moves through it.

Binding and transport of ligands in proteins is essential, in particular in globular proteins which often exhibit internal cavities. In truncated Hemoglobin N (trHbN) these cavities are arranged as a network with particular connectivities. The present work supports the notion that ligand diffusion in trHbN is an active process and coupled to the protein dynamics. Furthermore, transition probabilities between neighboring pockets depend on the location from where the ligand entered the transition, which is typical for non-Markovian processes. Hence, ligand migration in trHbN exhibits memory effects due to dynamical coupling between the protein and ligand motion.

The cellular membrane as a mediator for small molecule interaction with membrane proteins

The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Analyzing the electrogenicity of cytochrome c oxidase

Measurements of voltage changes in response to charge separation within membrane proteins can offer fundamental information on spectroscopically "invisible" steps. For example, results from studies of voltage changes associated with electron and proton transfer in cytochrome c oxidase could, in principle, be used to discriminate between different theoretical models describing the molecular mechanism of proton pumping. Earlier analyses of data from these measurements have been based on macroscopic considerations that may not allow for exploring the actual molecular mechanisms. Here, we have used a coarse-grained model describing the relation between observed voltage changes and specific charge-transfer reactions, which includes an explicit description of the membrane, the electrolytes, and the electrodes. The results from these calculations offer mechanistic insights at the molecular level. Our main conclusion is that previously assumed mechanistic evidence that was based on electrogenic measurements is not unique. However, the ability of our calculations to obtain reliable voltage changes means that we have a tool that can be used to describe a wide range of electrogenic charge transfers in channels and transporters, by combining voltage measurements with other experiments and simulations to analyze new mechanistic proposals.

Structure of a bd oxidase indicates similar mechanisms for membrane-integrated oxygen reductases

The cytochrome bd oxidases are terminal oxidases that are present in bacteria and archaea. They reduce molecular oxygen (dioxygen) to water, avoiding the production of reactive oxygen species. In addition to their contribution to the proton motive force, they mediate viability under oxygen-related stress conditions and confer tolerance to nitric oxide, thus contributing to the virulence of pathogenic bacteria. Here we present the atomic structure of the bd oxidase from Geobacillus thermodenitrificans, revealing a pseudosymmetrical subunit fold. The arrangement and order of the heme cofactors support the conclusions from spectroscopic measurements that the cleavage of the dioxygen bond may be mechanistically similar to that in the heme-copper-containing oxidases, even though the structures are completely different.

Control of the Position of Oxygen Delivery in Soybean Lipoxygenase-1 by Amino Acid Side Chains within a Gas Migration Channel

Understanding gas migration pathways is critical to unraveling structure-function relationships in enzymes that process gaseous substrates such as O2, H2, and N2 This work investigates the role of a defined pathway for O2 in regulating the peroxidation of linoleic acid by soybean lipoxygenase 1. Computational and mutagenesis studies provide strong support for a dominant delivery channel that shuttles molecular oxygen to a specific region of the active site, thereby ensuring the regio- and stereospecificity of product. Analysis of reaction kinetics and product distribution in channel mutants also reveals a plasticity to the gas migration pathway. The findings show that a single site mutation (I553W) limits oxygen accessibility to the active site, greatly increasing the fraction of substrate that reacts with oxygen free in solution. They also show how a neighboring site mutation (L496W) can result in a redirection of oxygen toward an alternate position of the substrate, changing the regio- and stereospecificity of peroxidation. The present data indicate that modest changes in a protein scaffold may modulate the access of small gaseous molecules to enzyme-bound substrates.

The CO Photodissociation and Recombination Dynamics of the W172Y/F282T Ligand Channel Mutant of Rhodobacter sphaeroides aa3 Cytochrome c Oxidase

In the ligand channel of the cytochrome c oxidase from Rhodobacter sphaeroides (Rs aa3 ) W172 and F282 have been proposed to generate a constriction that may slow ligand access to and from the active site. To explore this issue, the tryptophan and phenylalanine residues in Rs aa3 were mutated to the less bulky tyrosine and threonine residues, respectively, which occupy these sites in Thermus thermophilus (Tt) ba3 cytochrome oxidase. The CO photolysis and recombination dynamics of the reduced wild-type Rs aa3 and the W172Y/F282T mutant were investigated using time-resolved optical absorption spectroscopy. The spectral changes associated with the multiple processes are attributed to different conformers. The major CO recombination process (44 μs) in the W172Y/F282T mutant is ~500 times faster than the predominant CO recombination process in the wild-type enzyme (~23 ms). Classical dynamic simulations of the wild-type enzyme and double mutant showed significant structural changes at the active site in the mutant, including movement of the heme a3 ring-D propionate toward CuB and reduced binuclear center cavity volume. These structural changes effectively close the ligand exit pathway from the binuclear center, providing a basis for the faster CO recombination in the double mutant.

A broken-symmetry density functional study of structures, energies, and protonation states along the catalytic O-O bond cleavage pathway in ba3 cytochrome c oxidase from Thermus thermophilus

Broken-symmetry density functional calculations have been performed on the [Fea3, CuB] dinuclear center (DNC) of ba3 cytochrome c oxidase from Thermus thermophilus in the states of [Fea3(3+)-(HO2)(-)-CuB(2+), Tyr237(-)] and [Fea3(4+)[double bond, length as m-dash]O(2-), OH(-)-CuB(2+), Tyr237˙], using both PW91-D3 and OLYP-D3 functionals. Tyr237 is a special tyrosine cross-linked to His233, a ligand of CuB. The calculations have shown that the DNC in these states strongly favors the protonation of His376, which is above propionate-A, but not of the carboxylate group of propionate-A. The energies of the structures obtained by constrained geometry optimizations along the O-O bond cleavage pathway between [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] and [Fea3(4+)[double bond, length as m-dash]O(2-)HO(-)-CuB(2+), Tyr237˙] have also been calculated. The transition of [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] → [Fea3(4+)[double bond, length as m-dash]O(2-)HO(-)-CuB(2+), Tyr237˙] shows a very small barrier, which is less than 3.0/2.0 kcal mol(-1) in PW91-D3/OLYP-D3 calculations. The protonation state of His376 does not affect this O-O cleavage barrier. The rate limiting step of the transition from state A (in which O2 binds to Fea3(2+)) to state PM ([Fea3(4+)[double bond, length as m-dash]O(2-), OH(-)-CuB(2+), Tyr237˙], where the O-O bond is cleaved) in the catalytic cycle is, therefore, the proton transfer originating from Tyr237 to O-O to form the hydroperoxo [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] state. The importance of His376 in proton uptake and the function of propionate-A/neutral-Asp372 as a gate to prevent the proton from back-flowing to the DNC are also shown.

All the O2 Consumed by Thermus thermophilus Cytochrome ba3 Is Delivered to the Active Site through a Long, Open Hydrophobic Tunnel with Entrances within the Lipid Bilayer

Cytochrome ba3 is a proton-pumping heme-copper oxygen reductase from the extreme thermophile Thermus thermophilus. Despite the fact that the enzyme's active site is buried deep within the protein, the apparent second order rate constant for the initial binding of O2 to the active-site heme has been experimentally found to be 10(9) M(-1) s(-1) at 298 K, at or near the diffusion limit, and 2 orders of magnitude faster than for O2 binding to myoglobin. To provide quantitative and microscopic descriptions of the O2 delivery pathway and mechanism in cytochrome ba3, extensive molecular dynamics simulations of the enzyme in its membrane-embedded form have been performed, including different protocols of explicit ligand sampling (flooding) simulations with O2, implicit ligand sampling analysis, and in silico mutagenesis. The results show that O2 diffuses to the active site exclusively via a Y-shaped hydrophobic tunnel with two 25-Å long membrane-accessible branches that coincide with the pathway previously suggested by the crystallographically identified xenon binding sites. The two entrances of the bifurcated tunnel of cytochrome ba3 are located within the lipid bilayer, where O2 is preferentially partitioned from the aqueous phase. The largest barrier to O2 migration within the tunnel is estimated to be only 1.5 kcal/mol, allowing O2 to reach the enzyme active site virtually impeded by one-dimensional diffusion once it reaches a tunnel entrance at the protein surface. Unlike other O2-utilizing proteins, the tunnel is "open" with no transient barriers observed due to protein dynamics. This unique low-barrier passage through the protein ensures that O2 transit through the protein is never rate-limiting.

Manganese Compounds as Water-Oxidizing Catalysts: From the Natural Water-Oxidizing Complex to Nanosized Manganese Oxide Structures

All cyanobacteria, algae, and plants use a similar water-oxidizing catalyst for water oxidation. This catalyst is housed in Photosystem II, a membrane-protein complex that functions as a light-driven water oxidase in oxygenic photosynthesis. Water oxidation is also an important reaction in artificial photosynthesis because it has the potential to provide cheap electrons from water for hydrogen production or for the reduction of carbon dioxide on an industrial scale. The water-oxidizing complex of Photosystem II is a Mn-Ca cluster that oxidizes water with a low overpotential and high turnover frequency number of up to 25-90 molecules of O2 released per second. In this Review, we discuss the atomic structure of the Mn-Ca cluster of the Photosystem II water-oxidizing complex from the viewpoint that the underlying mechanism can be informative when designing artificial water-oxidizing catalysts. This is followed by consideration of functional Mn-based model complexes for water oxidation and the issue of Mn complexes decomposing to Mn oxide. We then provide a detailed assessment of the chemistry of Mn oxides by considering how their bulk and nanoscale properties contribute to their effectiveness as water-oxidizing catalysts.

PrinCCes: Continuity-based geometric decomposition and systematic visualization of the void repertoire of proteins

Grooves and pockets on the surface, channels through the protein, the chambers or cavities, and the tunnels connecting the internal points to each other or to the external fluid environment are fundamental determinants of a wide range of biological functions. PrinCCes (Protein internal Channel & Cavity estimation) is a computer program supporting the visualization of voids. It includes a novel algorithm for the decomposition of the entire void volume of the protein or protein complex to individual entities. The decomposition is based on continuity. An individual void is defined by uninterrupted extension in space: a spherical probe can freely move between any two internal locations of a continuous void. Continuous voids are detected irrespective of their topological complexity, they may contain any number of holes and bifurcations. The voids of a protein can be visualized one by one or in combinations as triangulated surfaces. The output is automatically exported to free VMD (Visual Molecular Dynamics) or Chimera software, allowing the 3D rotation of the surfaces and the production of publication quality images. PrinCCes with graphic user interface and command line versions are available for MS Windows and Linux. The source code and executable can be downloaded at any of the following links: http://scholar.semmelweis.hu/czirjakgabor/s/princces/#t1 https://github.com/CzirjakGabor/PrinCCes http://1drv.ms/1bP9iJ3.

Docking and migration of carbon monoxide in nitrogenase: the case for gated pockets from infrared spectroscopy and molecular dynamics

Evidence of a CO docking site near the FeMo cofactor in nitrogenase has been obtained by Fourier transform infrared spectroscopy-monitored low-temperature photolysis. We investigated the possible migration paths for CO from this docking site using molecular dynamics calculations. The simulations support the notion of a gas channel with multiple internal pockets from the active site to the protein exterior. Travel between pockets is gated by the motion of protein residues. Implications for the mechanism of nitrogenase reactions with CO and N2 are discussed.

Substrate pathways in the nitrogenase MoFe protein by experimental identification of small molecule binding sites

In the nitrogenase molybdenum-iron (MoFe) protein, we have identified five potential substrate access pathways from the protein surface to the FeMo-cofactor (the active site) or the P-cluster using experimental structures of Xe pressurized into MoFe protein crystals from Azotobacter vinelandii and Clostridium pasteurianum. Additionally, all published structures of the MoFe protein, including those from Klebsiella pneumoniae, were analyzed for the presence of nonwater, small molecules bound to the protein interior. Each pathway is based on identification of plausible routes from buried small molecule binding sites to both the protein surface and a metallocluster. Of these five pathways, two have been previously suggested as substrate access pathways. While the small molecule binding sites are not conserved among the three species of MoFe protein, residues lining the pathways are generally conserved, indicating that the proposed pathways may be accessible in all three species. These observations imply that there is unlikely a unique pathway utilized for substrate access from the protein surface to the active site; however, there may be preferred pathways such as those described here.

Mutation of a single residue in the ba3 oxidase specifically impairs protonation of the pump site

The ba3-type cytochrome c oxidase from Thermus thermophilus is a membrane-bound protein complex that couples electron transfer to O2 to proton translocation across the membrane. To elucidate the mechanism of the redox-driven proton pumping, we investigated the kinetics of electron and proton transfer in a structural variant of the ba3 oxidase where a putative "pump site" was modified by replacement of Asp372 by Ile. In this structural variant, proton pumping was uncoupled from internal electron transfer and O2 reduction. The results from our studies show that proton uptake to the pump site (time constant ∼65 μs in the wild-type cytochrome c oxidase) was impaired in the Asp372Ile variant. Furthermore, a reaction step that in the wild-type cytochrome c oxidase is linked to simultaneous proton uptake and release with a time constant of ∼1.2 ms was slowed to ∼8.4 ms, and in Asp372Ile was only associated with proton uptake to the catalytic site. These data identify reaction steps that are associated with protonation and deprotonation of the pump site, and point to the area around Asp372 as the location of this site in the ba3 cytochrome c oxidase.

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.

Conserved glycine 232 in the ligand channel of ba3 cytochrome oxidase from Thermus thermophilus

Knowing how the protein environment modulates ligand pathways and redox centers in the respiratory heme-copper oxidases is fundamental for understanding the relationship between the structure and function of these enzymes. In this study, we investigated the reactions of O2 and NO with the fully reduced G232V mutant of ba3 cytochrome c oxidase from Thermus thermophilus (Tt ba3) in which a conserved glycine residue in the O2 channel of the enzyme was replaced with a bulkier valine residue. Previous studies of the homologous mutant of Rhodobacter sphaeroides aa3 cytochrome c oxidase suggested that the valine completely blocked the access of O2 to the active site [Salomonsson, L., et al. (2004) Proc. Natl. Acad. Sci. U.S.A. 101, 11617-11621]. Using photolabile O2 and NO carriers, we find by using time-resolved optical absorption spectroscopy that the rates of O2 and NO binding are not significantly affected in the Tt ba3 G232V mutant. Classical molecular dynamics simulations of diffusion of O2 to the active site in the wild-type enzyme and G232V mutant show that the insertion of the larger valine residue in place of the glycine appears to open up other O2 and NO exit/entrance pathways that allow these ligands unhindered access to the active site, thus compensating for the larger valine residue.

Linking chemical electron-proton transfer to proton pumping in cytochrome c oxidase: broken-symmetry DFT exploration of intermediates along the catalytic reaction pathway of the iron-copper dinuclear complex

After a summary of the problem of coupling electron and proton transfer to proton pumping in cytochrome c oxidase, we present the results of our earlier and recent density functional theory calculations for the dinuclear Fe-a3-CuB reaction center in this enzyme. A specific catalytic reaction wheel diagram is constructed from the calculations, based on the structures and relative energies of the intermediate states of the reaction cycle. A larger family of tautomers/protonation states is generated compared to our earlier work, and a new lowest-energy pathway is proposed. The entire reaction cycle is calculated for the new smaller model (about 185-190 atoms), and two selected arcs of the wheel are chosen for calculations using a larger model (about 205 atoms). We compare the structural and redox energetics and protonation calculations with available experimental data. The reaction cycle map that we have built is positioned for further improvement and testing against experiment.

The pathway of O₂to the active site in heme-copper oxidases

The route of O₂to and from the high-spin heme in heme-copper oxidases has generally been believed to emulate that of carbon monoxide (CO). Time-resolved and stationary infrared experiments in our laboratories of the fully reduced CO-bound enzymes, as well as transient optical absorption saturation kinetics studies as a function of CO pressure, have provided strong support for CO binding to CuB⁺ on the pathway to and from the high-spin heme. The presence of CO on CuB⁺ suggests that O₂binding may be compromised in CO flow-flash experiments. Time-resolved optical absorption studies show that the rate of O₂and NO binding in the bovine enzyme (1 × 10⁸M⁻¹s⁻¹) is unaffected by the presence of CO, which is consistent with the rapid dissociation (t½ = 1.5μs) of CO from CuB⁺. In contrast, in Thermus thermophilus (Tt) cytochrome ba3 the O₂and NO binding to heme a3 slows by an order of magnitude in the presence of CO (from 1 × 10⁹ to 1 × 10⁸M⁻¹s⁻¹), but is still considerably faster (~10μs at 1atm O₂) than the CO off-rate from CuB in the absence of O₂(milliseconds). These results show that traditional CO flow-flash experiments do not give accurate results for the physiological binding of O₂and NO in Tt ba3, namely, in the absence of CO. They also raise the question whether in CO flow-flash experiments on Tt ba3 the presence of CO on CuB⁺ impedes the binding of O₂to CuB⁺ or, if O₂does not bind to CuB⁺ prior to heme a3, whether the CuB⁺-CO complex sterically restricts access of O₂to the heme. Both possibilities are discussed, and we argue that O₂binds directly to heme a3 in Tt ba3, causing CO to dissociate from CuB⁺ in a concerted manner through steric and/or electronic effects. This would allow CuB⁺ to function as an electron donor during the fast (5μs) breaking of the OO bond. These results suggest that the binding of CO to CuB⁺ on the path to and from heme a3 may not be applicable to O₂and NO in all heme-copper oxidases. This article is part of a Special Issue entitled: Vibrational spectroscopies and bioenergetic systems.

Density functional study for the bridged dinuclear center based on a high-resolution X-ray crystal structure of ba3 cytochrome c oxidase from Thermus thermophilus

Strong electron density for a peroxide type dioxygen species bridging the Fea3 and CuB dinuclear center (DNC) was observed in the high-resolution (1.8 Å) X-ray crystal structures (PDB entries 3S8G and 3S8F) of ba3 cytochrome c oxidase (CcO) from Thermus thermophilus. The crystals represent the as-isolated X-ray photoreduced CcO structures. The bridging peroxide was proposed to arise from the recombination of two radiation-produced HO(•) radicals formed either very near to or even in the space between the two metals of the DNC. It is unclear whether this peroxide species is in the O2(2-), O2(•)(-), HO2(-), or the H2O2 form and what is the detailed electronic structure and binding geometry including the DNC. In order to answer what form of this dioxygen species was observed in the DNC of the 1.8 Å X-ray CcO crystal structure (3S8G), we have applied broken-symmetry density functional theory (BS-DFT) geometric and energetic calculations (using OLYP potential) on large DNC cluster models with different Fea3-CuB oxidation and spin states and with O2(2-), O2(•)(-), HO2(-), or H2O2 in the bridging position. By comparing the DFT optimized geometries with the X-ray crystal structure (3S8G), we propose that the bridging peroxide is HO2(-). The X-ray crystal structure is likely to represent the superposition of the Fea3(2+)-(HO2(-))-CuB(+) DNC's in different states (Fe(2+) in low spin (LS), intermediate spin (IS), or high spin (HS)) with the majority species having the proton of the HO2(-) residing on the oxygen atom (O1) which is closer to the Fea3(2+) site in the Fea3(2+)-(HO-O)(-)-CuB(+) conformation. Our calculations show that the side chain of Tyr237 is likely trapped in the deprotonated Tyr237(-) anion form in the 3S8G X-ray crystal structure.

Computational prediction and in vitro analysis of potential physiological ligands of the bile acid binding site in cytochrome c oxidase

A conserved bile acid site has been crystallographically defined in the membrane domain of mammalian and Rhodobacter sphaeroides cytochrome c oxidase (RsCcO). Diverse amphipathic ligands were shown previously to bind to this site and affect the electron transfer equilibrium between heme a and a3 cofactors by blocking the K proton uptake path. Current studies identify physiologically relevant ligands for the bile acid site using a novel three-pronged computational approach: ROCS comparison of ligand shape and electrostatics, SimSite3D comparison of ligand binding site features, and SLIDE screening of potential ligands by docking. Identified candidate ligands include steroids, nicotinamides, flavins, nucleotides, retinoic acid, and thyroid hormones, which are predicted to make key protein contacts with the residues involved in bile acid binding. In vitro oxygen consumption and ligand competition assays on RsCcO wildtype and its Glu101Ala mutant support regulatory activity and specificity of some of these ligands. An ATP analog and GDP inhibit RsCcO under low substrate conditions, while fusidic acid, cholesteryl hemisuccinate, retinoic acid, and T3 thyroid hormone are more potent inhibitors under both high and low substrate conditions. The sigmoidal kinetics of RsCcO inhibition in the presence of certain nucleotides is reminiscent of previously reported ATP inhibition of mammalian CcO, suggesting regulation involving the conserved core subunits of both mammalian and bacterial oxidases. Ligand binding to the bile acid site is noncompetitive with respect to cytochrome c and appears to arrest CcO in a semioxidized state with some resemblance to the "resting" state of the enzyme.

Single mutations that redirect internal proton transfer in the ba3 oxidase from Thermus thermophilus

The ba3-type cytochrome c oxidase from Thermus thermophilus is a membrane-bound proton pump. Results from earlier studies have shown that with the aa3-type oxidases proton uptake to the catalytic site and "pump site" occurs simultaneously. However, with ba3 oxidase the pump site is loaded before proton transfer to the catalytic site because the proton transfer to the latter is slower than that with the aa3 oxidases. In addition, the timing of formation and decay of catalytic intermediates is different in the two types of oxidases. In the present study, we have investigated two mutant ba3 CytcOs in which residues of the proton pathway leading to the catalytic site as well as the pump site were exchanged, Thr312Val and Tyr244Phe. Even though ba3 CytcO uses only a single proton pathway for transfer of the substrate and "pumped" protons, the amino-acid residue substitutions had distinctly different effects on the kinetics of proton transfer to the catalytic site and the pump site. The results indicate that the rates of these reactions can be modified independently by replacement of single residues within the proton pathway. Furthermore, the data suggest that the Thr312Val and Tyr244Phe mutations interfere with a structural rearrangement in the proton pathway that is rate limiting for proton transfer to the catalytic site.

Current advances in research of cytochrome c oxidase

The function of cytochrome c oxidase as a biomolecular nanomachine that transforms energy of redox reaction into protonmotive force across a biological membrane has been subject of intense research, debate, and controversy. The structure of the enzyme has been solved for several organisms; however details of its molecular mechanism of proton pumping still remain elusive. Particularly, the identity of the proton pumping site, the key element of the mechanism, is still open to dispute. The pumping mechanism has been for a long time one of the key unsolved issues of bioenergetics and biochemistry, but with the accelerating progress in this field many important details and principles have emerged. Current advances in cytochrome oxidase research are reviewed here, along with a brief discussion of the most complete proton pumping mechanism proposed to date, and a molecular basis for control of its efficiency.

Ligand access to the active site in Thermus thermophilus ba(3) and bovine heart aa(3) cytochrome oxidases

Knowledge of the structure and dynamics of the ligand channel(s) in heme-copper oxidases is critical for understanding how the protein environment modulates the functions of these enzymes. Using photolabile NO and O(2) carriers, we recently found that NO and O(2) binding in Thermus thermophilus (Tt) ba(3) is ~10 times faster than in the bovine enzyme, indicating that inherent structural differences affect ligand access in these enzymes. Using X-ray crystallography, time-resolved optical absorption measurements, and theoretical calculations, we investigated ligand access in wild-type Tt ba(3) and the mutants, Y133W, T231F, and Y133W/T231F, in which tyrosine and threonine in the O(2) channel of Tt ba(3) are replaced by the corresponding bulkier tryptophan and phenylalanine, respectively, present in the aa(3) enzymes. NO binding in Y133W and Y133W/T231F was found to be 5 times slower than in wild-type ba(3) and the T231F mutant. The results show that the Tt ba(3) Y133W mutation and the bovine W126 residue physically impede NO access to the binuclear center. In the bovine enzyme, there is a hydrophobic "way station", which may further slow ligand access to the active site. Classical simulations of diffusion of Xe to the active sites in ba(3) and bovine aa(3) show conformational freedom of the bovine F238 and the F231 side chain of the Tt ba(3) Y133W/T231F mutant, with both residues rotating out of the ligand channel, resulting in no effect on ligand access in either enzyme.

From static structure to living protein: computational analysis of cytochrome c oxidase main-chain flexibility

Crystallographic structure and deuterium accessibility comparisons of CcO in different redox states have suggested conformational changes of mechanistic significance. To predict the intrinsic flexibility and low energy motions in CcO, this work has analyzed available high-resolution crystallographic structures with ProFlex and elNémo computational methods. The results identify flexible regions and potential conformational changes in CcO that correlate well with published structural and biochemical data and provide mechanistic insights. CcO is predicted to undergo rotational motions on the interior and exterior of the membrane, driven by transmembrane helical tilting and bending, coupled with rocking of the β-sheet domain. Consequently, the proton K-pathway becomes sufficiently flexible for internal water molecules to alternately occupy upper and lower parts of the pathway, associated with conserved Thr-359 and Lys-362 residues. The D-pathway helices are found to be relatively rigid, with a highly flexible entrance region involving the subunit I C-terminus, potentially regulating the uptake of protons. Constriction and dilation of hydrophobic channels in RsCcO suggest regulation of the oxygen supply to the binuclear center. This analysis points to coupled conformational changes in CcO and their potential to influence both proton and oxygen access.

Observation of ligand transfer in ba3 oxidase from Thermus thermophilus: simultaneous FTIR detection of photolabile heme a3(2+)-CN and transient Cu(B)(2+)-CN complexes

FTIR and light-minus-dark FTIR spectroscopy have been employed to investigate the reaction of oxidized and fully reduced ba(3) oxidase with cyanide. The characterization of the structures of the bound CN(-) in the binuclear heme Fe-Cu(B) center is essential, given that a central issue in the function of ba(3) oxidase is the extent to which the partially reduced substrates interact with the two metals. In the reaction of oxidized ba(3) oxidase with cyanide the initially formed heme a(3)(3+)-C≡N-Cu(B)(2+) species with ν(CN) frequency at 2152 cm(-1) was replaced by a photolabile complex with a frequency at 2075 cm(-1) characteristic of heme a(3)(2+)-CN(-). Photolysis of the heme a(3)(2+)-CN(-) adduct produced a band at 2146 cm(-1) attributed to the formation of a transient Cu(B)(2+)-CN(-) complex. All forms are pH independent between pH 5.5-9.5 and at pD 7.5 indicating the absence of ionizable groups that influence the properties of the cyanide complexes. In contrast to previous reports, our results show that CN(-) does not bind simultaneously to both heme a(3)(2+) and Cu(B)(2+) to form the mixed valence a(3)(2+)-CN·Cu(B)(2+)CN species. The photolysis products of the heme a(3)(2+)-CN(-)/Cu(B)(2+) and heme a(3)(2+)-CN(-)/Cu(B)(1+) species are different suggesting that relaxation dynamics in the binuclear center following ligand photodissociation are dependent on the oxidation state of Cu(B).

Mobility of Xe atoms within the oxygen diffusion channel of cytochrome ba(3) oxidase

We use a form of "freeze-trap, kinetic crystallography" to explore the migration of Xe atoms away from the dinuclear heme a(3)/Cu(B) center in Thermus thermophilus cytochrome ba(3) oxidase. This enzyme is a member of the heme-copper oxidase superfamily and is thus crucial for dioxygen-dependent life. The mechanisms involved in the migration of oxygen, water, electrons, and protons into and/or out of the specialized channels of the heme-copper oxidases are generally not well understood. Pressurization of crystals with Xe gas previously revealed a O(2) diffusion channel in cytochrome ba(3) oxidase that is continuous, Y-shaped, 18-20 Å in length and comprised of hydrophobic residues, connecting the protein surface within the bilayer to the a(3)-Cu(B) center in the active site. To understand movement of gas molecules within the O(2) channel, we performed crystallographic analysis of 19 Xe laden crystals freeze-trapped in liquid nitrogen at selected times between 0 and 480 s while undergoing outgassing at room temperature. Variation in Xe crystallographic occupancy at five discrete sites as a function of time leads to a kinetic model revealing relative degrees of mobility of Xe atoms within the channel. Xe egress occurs primarily through the channel formed by the Xe1 → Xe5 → Xe3 → Xe4 sites, suggesting that ingress of O(2) is likely to occur by the reverse of this process. The channel itself appears not to undergo significant structural changes during Xe migration, thereby indicating a passive role in this important physiological function.

Bioenergetics at extreme temperature: Thermus thermophilus ba(3)- and caa(3)-type cytochrome c oxidases

Seven years into the completion of the genome sequencing projects of the thermophilic bacterium Thermus thermophilus strains HB8 and HB27, many questions remain on its bioenergetic mechanisms. A key fact that is occasionally overlooked is that oxygen has a very limited solubility in water at high temperatures. The HB8 strain is a facultative anaerobe whereas its relative HB27 is strictly aerobic. This has been attributed to the absence of nitrate respiration genes from the HB27 genome that are carried on a mobilizable but highly-unstable plasmid. In T. thermophilus, the nitrate respiration complements the primary aerobic respiration. It is widely known that many organisms encode multiple biochemically-redundant components of the respiratory complexes. In this minireview, the presence of the two cytochrome c oxidases (CcO) in T. thermophilus, the ba(3)- and caa(3)-types, is outlined along with functional considerations. We argue for the distinct evolutionary histories of these two CcO including their respective genetic and molecular organizations, with the caa(3)-oxidase subunits having been initially 'fused'. Coupled with sequence analysis, the ba(3)-oxidase crystal structure has provided evolutionary and functional information; for example, its subunit I is more closely related to archaeal sequences than bacterial and the substrate-enzyme interaction is hydrophobic as the elevated growth temperature weakens the electrostatic interactions common in mesophiles. Discussion on the role of cofactors in intra- and intermolecular electron transfer and proton pumping mechanism is also included.