Kinetic studies of the reactions of O(2) and NO with reduced Thermus thermophilus ba(3) and bovine aa(3) using photolabile carriers

Characterization and evaluation of photolabile (µ-peroxo)(µ-hydroxo)bis[bis(bipyridyl)cobalt caged oxygen compounds to facilitate time-resolved crystallographic studies of cytochrome c oxidase

Photolabile (µ-peroxo)(µ-hydroxo)bis[bis(bipyridyl)-cobalt-based caged oxygen compounds have been synthesized and characterized by optical absorbance spectroscopy, X-ray crystallography. and the quantum yield and redox stability were investigated. Furthermore, conditions were established where redox incompatibilities encountered between caged oxygen compounds and oxygen-dependant cytochrome c oxidase (CcO) could be circumvented. Herein, we demonstrate that millimolar concentrations of molecular oxygen can be released from a caged oxygen compound with spatio-temporal control upon laser excitation, triggering enzymatic turnover in cytochrome c oxidase. Spectroscopic evidence confirms the attainment of a homogeneous reaction initiation at concentrations and conditions relevant for further crystallography studies. This was demonstrated by the oxidizing microcrystals of reduced CcO by liberation of millimolar concentrations of molecular oxygen from a caged oxygen compound. We believe this will expand the scope of available techniques for the detailed investigation of oxygen-dependant enzymes with its native substrate and facilitate further time-resolved X-ray based studies such as wide/small angle X-ray scattering and serial femtosecond crystallography.

Time-resolved serial crystallography to track the dynamics of carbon monoxide in the active site of cytochrome c oxidase

Cytochrome c oxidase (CcO) is part of the respiratory chain and contributes to the electrochemical membrane gradient in mitochondria as well as in many bacteria, as it uses the energy released in the reduction of oxygen to pump protons across an energy-transducing biological membrane. Here, we use time-resolved serial femtosecond crystallography to study the structural response of the active site upon flash photolysis of carbon monoxide (CO) from the reduced heme a3 of ba3-type CcO. In contrast with the aa3-type enzyme, our data show how CO is stabilized on CuB through interactions with a transiently ordered water molecule. These results offer a structural explanation for the extended lifetime of the CuB-CO complex in ba3-type CcO and, by extension, the extremely high oxygen affinity of the enzyme.

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.

The Reactions of O2 and NO with Mixed-Valence ba3 Cytochrome c Oxidase from Thermus thermophilus

Earlier CO flow-flash experiments on the fully reduced Thermus thermophilus ba₃ (Tt ba₃) cytochrome oxidase revealed that O₂ binding was slowed down by a factor of 10 in the presence of CO (Szundi et al., 2010, PNAS 107, 21010-21015). The goal of the current study is to explore whether the long apparent lifetime (∼50 ms) of the Cu_B⁺-CO complex generated upon photolysis of the CO-bound mixed-valence Tt ba₃ (Koutsoupakis et al., 2019, Acc. Chem. Res. 52, 1380-1390) affects O₂ and NO binding and the ability of Cu_B to act as an electron donor during O-O bond splitting. The CO recombination, NO binding, and the reaction of mixed-valence Tt ba₃ with O₂ were investigated by time-resolved optical absorption spectroscopy using the CO flow-flash approach and photolabile O₂ and NO carriers. No electron backflow was detected after photolysis of the mixed-valence CO-bound Tt ba₃. The rate of O₂ and NO binding was two times slower than in the fully reduced enzyme in the presence of CO and 20 times slower than in the absence of CO. The purported long-lived Cu_B⁺-CO complex did not prevent O-O bond splitting and the resulting P_M formation, which was significantly faster (5-10 times) than in the bovine heart enzyme. We propose that O₂ binding to heme a₃ in Tt ba₃ causes CO to dissociate from Cu_B⁺ in a concerted manner through steric and/or electronic effects, thus allowing Cu_B⁺ to act as an electron donor in the mixed-valence enzyme. The significantly faster O₂ binding and O-O bond cleavage in Tt ba₃ compared to analogous steps in the aa₃ oxidases could reflect evolutionary adaptation of the enzyme to the microaerobic conditions of the T. thermophilus HB8 species.

The proton pumping bo oxidase from Vitreoscilla

The cytochrome bo3 quinol oxidase from Vitreoscilla (vbo3) catalyses oxidation of ubiquinol and reduction of O2 to H2O. Data from earlier studies suggested that the free energy released in this reaction is used to pump sodium ions instead of protons across a membrane. Here, we have studied the functional properties of heterologously expressed vbo3 with a variety of methods. (i) Following oxygen consumption with a Clark-type electrode, we did not observe a measurable effect of Na+ on the oxidase activity of purified vbo3 solubilized in detergent or reconstituted in liposomes. (ii) Using fluorescent dyes, we find that vbo3 does not pump Na+ ions, but H+ across the membrane, and that H+-pumping is not influenced by the presence of Na+. (iii) Using an oxygen pulse method, it was found that 2 H+/e- are ejected from proteoliposomes, in agreement with the values found for the H+-pumping bo3 oxidase of Escherichia coli (ecbo3). This coincides with the interpretation that 1 H+/e- is pumped across the membrane and 1 H+/e- is released during quinol oxidation. (iv) When the electron transfer kinetics of vbo3 upon reaction with oxygen were followed in single turnover experiments, a similar sequence of reaction steps was observed as reported for the E. coli enzyme and none of these reactions was notably affected by the presence of Na+. Overall the data show that vbo3 is a proton pumping terminal oxidase, behaving similarly to the Escherichia coli bo3 quinol oxidase.

Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function

Heme-copper oxidases (HCOs) are terminal enzymes on the mitochondrial or bacterial respiratory electron transport chain, which utilize a unique heterobinuclear active site to catalyze the 4H+/4e- reduction of dioxygen to water. This process involves a proton-coupled electron transfer (PCET) from a tyrosine (phenolic) residue and additional redox events coupled to transmembrane proton pumping and ATP synthesis. Given that HCOs are large, complex, membrane-bound enzymes, bioinspired synthetic model chemistry is a promising approach to better understand heme-Cu-mediated dioxygen reduction, including the details of proton and electron movements. This review encompasses important aspects of heme-O2 and copper-O2 (bio)chemistries as they relate to the design and interpretation of small molecule model systems and provides perspectives from fundamental coordination chemistry, which can be applied to the understanding of HCO activity. We focus on recent advancements from studies of heme-Cu models, evaluating experimental and computational results, which highlight important fundamental structure-function relationships. Finally, we provide an outlook for future potential contributions from synthetic inorganic chemistry and discuss their implications with relevance to biological O2-reduction.

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.

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.

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.

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.

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.

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.

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.

Excitation wavelength dependent O2 release from copper(II)-superoxide compounds: laser flash-photolysis experiments and theoretical studies

Irradiation of the copper(II)-superoxide synthetic complexes [(TMG3tren)Cu(II)(O2)](+) (1) and [(PV-TMPA)Cu(II)(O2)](+) (2) with visible light resulted in direct photogeneration of O2 gas at low temperature (from -40 °C to -70 °C for 1 and from -125 to -135 °C for 2) in 2-methyltetrahydrofuran (MeTHF) solvent. The yield of O2 release was wavelength dependent: λexc = 436 nm, ϕ = 0.29 (for 1), ϕ = 0.11 (for 2), and λexc = 683 nm, ϕ = 0.035 (for 1), ϕ = 0.078 (for 2), which was followed by fast O2-recombination with [(TMG3tren)Cu(I)](+) (3) and [(PV-TMPA)Cu(I)](+) (4). Enthalpic barriers for O2 rebinding to the copper(I) center (∼10 kJ mol(-1)) and for O2 dissociation from the superoxide compound 1 (45 kJ mol(-1)) were determined. TD-DFT studies, carried out for 1, support the experimental results confirming the dissociative character of the excited states formed upon blue- or red-light laser excitation.

An engineered heme-copper center in myoglobin: CO migration and binding

We have investigated CO migration and binding in CuBMb, a copper-binding myoglobin double mutant (L29H-F43H), by using Fourier transform infrared spectroscopy and flash photolysis over a wide temperature range. This mutant was originally engineered with the aim to mimic the catalytic site of heme-copper oxidases. Comparison of the wild-type protein Mb and CuBMb shows that the copper ion in the distal pocket gives rise to significant effects on ligand binding to the heme iron. In Mb and copper-free CuBMb, primary and secondary ligand docking sites are accessible upon photodissociation. In copper-bound CuBMb, ligands do not migrate to secondary docking sites but rather coordinate to the copper ion. Ligands entering the heme pocket from the outside normally would not be captured efficiently by the tight distal pocket housing the two additional large imidazole rings. Binding at the Cu ion, however, ensures efficient trapping in CuBMb. The Cu ion also restricts the motions of the His64 side chain, which is the entry/exit door for ligand movement into the active site, and this restriction results in enhanced geminate and slow bimolecular CO rebinding. These results support current mechanistic views of ligand binding in hemoglobins and the role of the CuB in the active of heme-copper oxidases. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

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.

Spectral identification of intermediates generated during the reaction of dioxygen with the wild-type and EQ(I-286) mutant of Rhodobacter sphaeroides cytochrome c oxidase

Cytochrome c oxidase from Rhodobacter sphaeroides is frequently used to model the more complex mitochondrial enzyme. The O(2) reduction in both enzymes is generally described by a unidirectional mechanism involving the sequential formation of the ferrous-oxy complex (compound A), the P(R) state, the oxyferryl F form, and the oxidized state. In this study we investigated the reaction of dioxygen with the wild-type reduced R. sphaeroides cytochrome oxidase and the EQ(I-286) mutant using the CO flow-flash technique. Singular value decomposition and multiexponential fitting of the time-resolved optical absorption difference spectra showed that three apparent lifetimes, 18 μs, 53 μs, and 1.3 ms, are sufficient to fit the kinetics of the O(2) reaction of the wild-type enzyme. A comparison of the experimental intermediate spectra with the corresponding intermediate spectra of the bovine enzyme revealed that P(R) is not present in the reaction mechanism of the wild-type R. sphaeroides aa(3). Transient absorbance changes at 440 and 610 nm support this conclusion. For the EQ(I-286) mutant, in which a key glutamic residue in the D proton pathway is replaced by glutamine, two lifetimes, 16 and 108 μs, were observed. A spectral analysis of the intermediates shows that the O(2) reaction in the EQ(I-286) mutant terminates at the P(R) state, with 70% of heme a becoming oxidized. These results indicate significant differences in the kinetics of O(2) reduction between the bovine and wild-type R. sphaeroides aa(3) oxidases, which may arise from differences in the relative rates of internal electron and proton movements in the two enzymes.