1
|
Noodleman L, Götz AW, Han Du WG, Hunsicker-Wang L. Reaction pathways, proton transfer, and proton pumping in ba3 class cytochrome c oxidase: perspectives from DFT quantum chemistry and molecular dynamics. Front Chem 2023; 11:1186022. [PMID: 38188931 PMCID: PMC10766771 DOI: 10.3389/fchem.2023.1186022] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Accepted: 11/27/2023] [Indexed: 01/09/2024] Open
Abstract
After drawing comparisons between the reaction pathways of cytochrome c oxidase (CcO, Complex 4) and the preceding complex cytochrome bc1 (Complex 3), both being proton pumping complexes along the electron transport chain, we provide an analysis of the reaction pathways in bacterial ba3 class CcO, comparing spectroscopic results and kinetics observations with results from DFT calculations. For an important arc of the catalytic cycle in CcO, we can trace the energy pathways for the chemical protons and show how these pathways drive proton pumping of the vectorial protons. We then explore the proton loading network above the Fe heme a3-CuB catalytic center, showing how protons are loaded in and then released by combining DFT-based reaction energies with molecular dynamics simulations over states of that cycle. We also propose some additional reaction pathways for the chemical and vector protons based on our recent work with spectroscopic support.
Collapse
Affiliation(s)
- Louis Noodleman
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, CA, United States
| | - Andreas W. Götz
- San Diego Supercomputer Center, University of California San Diego, La Jolla, CA, United States
| | - Wen-Ge Han Du
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, CA, United States
| | | |
Collapse
|
2
|
Cao YC, Liao RZ. QM Calculations Revealed that Outer-Sphere Electron Transfer Boosted O-O Bond Cleavage in the Multiheme-Dependent Cytochrome bd Oxygen Reductase. Inorg Chem 2023; 62:4066-4075. [PMID: 36857027 DOI: 10.1021/acs.inorgchem.2c03742] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2023]
Abstract
The cytochrome bd oxygen reductase catalyzes the four-electron reduction of dioxygen to two water molecules. The structure of this enzyme reveals three heme molecules in the active site, which differs from that of heme-copper cytochrome c oxidase. The quantum chemical cluster approach was used to uncover the reaction mechanism of this intriguing metalloenzyme. The calculations suggested that a proton-coupled electron transfer reduction occurs first to generate a ferrous heme b595. This is followed by the dioxygen binding at the heme d center coupled with an outer-sphere electron transfer from the ferrous heme b595 to the dioxygen moiety, affording a ferric ion superoxide intermediate. A second proton-coupled electron transfer produces a heme d ferric hydroperoxide, which undergoes efficient O-O bond cleavage facilitated by an outer-sphere electron transfer from the ferrous heme b595 to the O-O σ* orbital and an inner-sphere proton transfer from the heme d hydroxyl group to the leaving hydroxide. The synergistic benefits of the two types of hemes rationalize the highly efficient oxygen reduction repertoire for the multi-heme-dependent cytochrome bd oxygen reductase family.
Collapse
Affiliation(s)
- Yu-Chen Cao
- Key Laboratory of Material Chemistry for Energy Conversion and Storage, Ministry of Education, Hubei Key Laboratory of Bioinorganic Chemistry and Materia Medica, Hubei Key Laboratory of Materials Chemistry and Service Failure, School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, Wuhan 430074, People's Republic of China
| | - Rong-Zhen Liao
- Key Laboratory of Material Chemistry for Energy Conversion and Storage, Ministry of Education, Hubei Key Laboratory of Bioinorganic Chemistry and Materia Medica, Hubei Key Laboratory of Materials Chemistry and Service Failure, School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, Wuhan 430074, People's Republic of China
| |
Collapse
|
3
|
Du WGH, Götz AW, Noodleman L. Mössbauer Property Calculations on Fea33+∙∙∙H2O∙∙∙CuB2+ Dinuclear Center Models of the Resting Oxidized as-Isolated Cytochrome c Oxidase. Chemphyschem 2022; 23:e202100831. [PMID: 35142420 PMCID: PMC9054037 DOI: 10.1002/cphc.202100831] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Revised: 02/03/2022] [Indexed: 11/24/2022]
Abstract
Mössbauer isomer shift and quadrupole splitting properties have been calculated using the OLYP‐D3(BJ) density functional method on previously obtained (W.‐G. Han Du, et al., Inorg Chem. 2020, 59, 8906–8915) geometry optimized Fea33+−H2O−CuB2+ dinuclear center (DNC) clusters of the resting oxidized (O state) “as‐isolated” cytochrome c oxidase (CcO). The calculated results are highly consistent with the available experimental observations. The calculations have also shown that the structural heterogeneities of the O state DNCs implicated by the Mössbauer experiments are likely consequences of various factors, particularly the variable positions of the central H2O molecule between the Fea33+ and CuB2+ sites in different DNCs, whether or not this central H2O molecule has H‐bonding interaction with another H2O molecule, the different spin states having similar energies for the Fea33+ sites, and whether the Fea33+ and CuB2+ sites are ferromagnetically or antiferromagnetically spin‐coupled.
Collapse
Affiliation(s)
- Wen-Ge Han Du
- The Scripps Research Institute, Integrative Structural and Computational Biology, UNITED STATES
| | | | - Louis Noodleman
- The Scripps Research Institute, Department of Integrative Structural and Computational Biology, Hz112, 10550 North Torrey Pines Road, 92037, La Jolla, UNITED STATES
| |
Collapse
|
4
|
Siletsky SA, Borisov VB. Proton Pumping and Non-Pumping Terminal Respiratory Oxidases: Active Sites Intermediates of These Molecular Machines and Their Derivatives. Int J Mol Sci 2021; 22:10852. [PMID: 34639193 PMCID: PMC8509429 DOI: 10.3390/ijms221910852] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Revised: 10/04/2021] [Accepted: 10/05/2021] [Indexed: 11/16/2022] Open
Abstract
Terminal respiratory oxidases are highly efficient molecular machines. These most important bioenergetic membrane enzymes transform the energy of chemical bonds released during the transfer of electrons along the respiratory chains of eukaryotes and prokaryotes from cytochromes or quinols to molecular oxygen into a transmembrane proton gradient. They participate in regulatory cascades and physiological anti-stress reactions in multicellular organisms. They also allow microorganisms to adapt to low-oxygen conditions, survive in chemically aggressive environments and acquire antibiotic resistance. To date, three-dimensional structures with atomic resolution of members of all major groups of terminal respiratory oxidases, heme-copper oxidases, and bd-type cytochromes, have been obtained. These groups of enzymes have different origins and a wide range of functional significance in cells. At the same time, all of them are united by a catalytic reaction of four-electron reduction in oxygen into water which proceeds without the formation and release of potentially dangerous ROS from active sites. The review analyzes recent structural and functional studies of oxygen reduction intermediates in the active sites of terminal respiratory oxidases, the features of catalytic cycles, and the properties of the active sites of these enzymes.
Collapse
Affiliation(s)
- Sergey A. Siletsky
- Department of Bioenergetics, Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Leninskie Gory, Moscow 119991, Russia
| | - Vitaliy B. Borisov
- Department of Molecular Energetics of Microorganisms, Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Leninskie Gory, Moscow 119991, Russia;
| |
Collapse
|
5
|
Blomberg MRA. Activation of O 2 and NO in heme-copper oxidases - mechanistic insights from computational modelling. Chem Soc Rev 2021; 49:7301-7330. [PMID: 33006348 DOI: 10.1039/d0cs00877j] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
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.
Collapse
Affiliation(s)
- Margareta R A Blomberg
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91, Stockholm, Sweden.
| |
Collapse
|
6
|
Reed CJ, Lam QN, Mirts EN, Lu Y. Molecular understanding of heteronuclear active sites in heme-copper oxidases, nitric oxide reductases, and sulfite reductases through biomimetic modelling. Chem Soc Rev 2021; 50:2486-2539. [PMID: 33475096 PMCID: PMC7920998 DOI: 10.1039/d0cs01297a] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
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.
Collapse
Affiliation(s)
- Christopher J Reed
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urban, IL 61801, USA.
| | - Quan N Lam
- Department of Biochemistry, University of Illinois at Urbana-Champaign, Urban, IL 61801, USA
| | - Evan N Mirts
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA and Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Yi Lu
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urban, IL 61801, USA. and Department of Biochemistry, University of Illinois at Urbana-Champaign, Urban, IL 61801, USA and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA and Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
7
|
Ghosh AC, Duboc C, Gennari M. Synergy between metals for small molecule activation: Enzymes and bio-inspired complexes. Coord Chem Rev 2021. [DOI: 10.1016/j.ccr.2020.213606] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
8
|
Noodleman L, Han Du WG, McRee D, Chen Y, Goh T, Götz AW. Coupled transport of electrons and protons in a bacterial cytochrome c oxidase-DFT calculated properties compared to structures and spectroscopies. Phys Chem Chem Phys 2021; 22:26652-26668. [PMID: 33231596 DOI: 10.1039/d0cp04848h] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Louis Noodleman
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA.
| | | | | | | | | | | |
Collapse
|
9
|
Electrocatalytic Oxygen Reduction at Multinuclear Metal Active Sites Inspired by Metalloenzymes. E-JOURNAL OF SURFACE SCIENCE AND NANOTECHNOLOGY 2020. [DOI: 10.1380/ejssnt.2020.81] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
|
10
|
Kim H, Sharma SK, Schaefer AW, Solomon EI, Karlin KD. Heme-Cu Binucleating Ligand Supports Heme/O 2 and Fe II-Cu I/O 2 Reactivity Providing High- and Low-Spin Fe III-Peroxo-Cu II Complexes. Inorg Chem 2019; 58:15423-15432. [PMID: 31657921 DOI: 10.1021/acs.inorgchem.9b02521] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The focus of this study is in the description of synthetic heme/copper/O2 chemistry employing a heme-containing binucleating ligand which provides a tridentate chelate for copper ion binding. The addition of O2 (-80 °C, tetrahydrofuran (THF) solvent) to the reduced heme compound (PImH)FeII (1), gives the oxy-heme adduct, formally a heme-superoxide complex FeIII-(O2•-) (2) (resonance Raman spectroscopy (rR): νO-O, 1171 cm-1 (Δ18O2, -61 cm-1); νFe-O, 575 cm-1 (Δ18O2, -24 cm-1)). Simple warming of 2 to room temperature regenerates reduced complex 1; this reaction is reversible, as followed by UV-vis spectroscopy. Complex 2 is electron paramagnetic resonance (EPR)-silent and exhibits upfield-shifted pyrrole resonances (δ 9.12 ppm) in 2H NMR spectroscopy, indicative of a six-coordinate low-spin heme. The coordination of the tethered imidazolyl arm to the heme-superoxide complex as an axial base ligand is suggested. We also report the new fully reduced heme-copper complex [(PImH)FeIICuI]+ (3), where the copper ion is bound to the tethered tridentate portion of PImH. This reacts with O2 to give a distinctive low-temperature-stable, high-spin (S = 2, overall) peroxo-bridged complex [(PImH)FeIII-(O22-)-CuII]+ (3a): λmax, 420 (Soret), 545, 565 nm; δpyrr, 93 ppm; νO-O, 799 cm-1 (Δ18O2, -48 cm-1); νFe-O, 524 cm-1 (Δ18O2, -23 cm-1). To 3a, the addition of dicyclohexylimidazole (DCHIm), which serves as a heme axial base, leads to low-spin (S = 0 overall) species complex [(DCHIm)(PImH)FeIII-(O22-)-CuII]+ (3b): λmax, 425 (Soret), 538 nm; δpyrr, 10.2 ppm; νO-O, 817 cm-1 (Δ18O2, -55 cm-1); νFe-O, 610 cm-1 (Δ18O2, -26 cm-1). These investigations into the characterization of the O2-adducts from (PImH)FeII (1) with/without additional copper chelation advance our understanding of the dioxygen reactivity of heme-only and heme/Cu-ligand heterobinuclear system, thus potentially relevant to O2 reduction in heme-copper oxidases or fuel-cell chemistry.
Collapse
Affiliation(s)
- Hyun Kim
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , United States
| | - Savita K Sharma
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , United States
| | - Andrew W Schaefer
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | - Edward I Solomon
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | - Kenneth D Karlin
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , United States
| |
Collapse
|
11
|
Han Du WG, Götz AW, Noodleman L. DFT Fe a3-O/O-O Vibrational Frequency Calculations over Catalytic Reaction Cycle States in the Dinuclear Center of Cytochrome c Oxidase. Inorg Chem 2019; 58:13933-13944. [PMID: 31566371 PMCID: PMC6839913 DOI: 10.1021/acs.inorgchem.9b01840] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Density functional vibrational frequency calculations have been performed on eight geometry optimized cytochrome c oxidase (CcO) dinuclear center (DNC) reaction cycle intermediates and on the oxymyoglobin (oxyMb) active site. The calculated Fe-O and O-O stretching modes and their frequency shifts along the reaction cycle have been compared with the available resonance Raman (rR) measurements. The calculations support the proposal that in state A[Fea33+-O2-•···CuB+] of CcO, O2 binds with Fea32+ in a similar bent end-on geometry to that in oxyMb. The calculations show that the observed 20 cm-1 shift of the Fea3-O stretching mode from the PR to F state is caused by the protonation of the OH- ligand on CuB2+ (PR[Fea34+═O2-···HO--CuB2+] → F[Fea34+═O2-···H2O-CuB2+]), and that the H2O ligand is still on the CuB2+ site in the rR identified F[Fea34+═O2-···H2O-CuB2+] state. Further, the observed rR band at 356 cm-1 between states PR and F is likely an O-Fea3-porphyrin bending mode. The observed 450 cm-1 low Fea3-O frequency mode for the OH active oxidized state has been reproduced by our calculations on a nearly symmetrically bridged Fea33+-OH-CuB2+ structure with a relatively long Fea3-O distance near 2 Å. Based on Badger's rule, the calculated Fea3-O distances correlate well with the calculated νFe-O-2/3 (νFe-O is the Fea3-O stretching frequency) with correlation coefficient R = 0.973.
Collapse
Affiliation(s)
- Wen-Ge Han Du
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037
| | - Andreas W. Götz
- San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive MC0505, La Jolla, CA 92093
| | - Louis Noodleman
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037
| |
Collapse
|
12
|
Geometric and Electronic Structure Contributions to O-O Cleavage and the Resultant Intermediate Generated in Heme-Copper Oxidases. J Am Chem Soc 2019; 141:10068-10081. [PMID: 31146528 DOI: 10.1021/jacs.9b04271] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This study investigates the mechanism of O-O bond cleavage in heme-copper oxidase (HCO) enzymes, combining experimental and computational insights from enzyme intermediates and synthetic models. It is determined that HCOs undergo a proton-initiated O-O cleavage mechanism where a single water molecule in the active site enables proton transfer (PT) from the cross-linked tyrosine to a peroxo ligand bridging the heme FeIII and CuII, and multiple H-bonding interactions lower the tyrosine p Ka. Due to sterics within the active site, the proton must either transfer initially to the O(Fe) (a high-energy intermediate), or from another residue over a ∼10 Å distance to reach the O(Cu) atom directly. While the distance between the H+ donor (Tyr) and acceptor (O(Cu)) results in a barrier to PT, this separation is critical for the low barrier to O-O cleavage as it enhances backbonding from Fe into the O22- σ* orbital. Thus, PT from Tyr precedes O-O elongation and is rate-limiting, consistent with available kinetic data. The electron transfers from tyrosinate after the barrier via a superexchange pathway provided by the cross-link, generating intermediate PM. PM is evaluated using available experimental data. The geometric structure contains an FeIV═O that is H-bonded to the CuII-OH. The electronic structure is a singlet, where the FeIV and CuII are antiferromagnetically coupled through the H-bond between the oxo(Fe) and hydroxo(Cu) ligands, while the CuII and Tyr• are ferromagnetically coupled due their delocalization into orthogonal magnetic orbitals on the cross-linked His residue. These findings provide critical insights into the mechanism of efficient O2 reduction in HCOs, and the nature of the PM intermediate that couples this reaction to proton pumping.
Collapse
|
13
|
Adam SM, Wijeratne GB, Rogler PJ, Diaz DE, Quist DA, Liu JJ, Karlin KD. Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function. Chem Rev 2018; 118:10840-11022. [PMID: 30372042 PMCID: PMC6360144 DOI: 10.1021/acs.chemrev.8b00074] [Citation(s) in RCA: 132] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
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.
Collapse
Affiliation(s)
- Suzanne M. Adam
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Gayan B. Wijeratne
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Patrick J. Rogler
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Daniel E. Diaz
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - David A. Quist
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Jeffrey J. Liu
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Kenneth D. Karlin
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
14
|
Han Du WG, Götz AW, Noodleman L. A Water Dimer Shift Activates a Proton Pumping Pathway in the P R → F Transition of ba 3 Cytochrome c Oxidase. Inorg Chem 2018; 57:1048-1059. [PMID: 29308889 DOI: 10.1021/acs.inorgchem.7b02461] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Broken-symmetry density functional calculations have been performed on the [Fea34+,CuB2+] state of the dinuclear center (DNC) for the PR → F part of the catalytic cycle of ba3 cytochrome c oxidase (CcO) from Thermus thermophilus (Tt), using the OLYP-D3-BJ functional. The calculations show that the movement of the H2O molecules in the DNC affects the pKa values of the residue side chains of Tyr237 and His376+, which are crucial for proton transfer/pumping in ba3 CcO from Tt. The calculated lowest energy structure of the DNC in the [Fea34+,CuB2+] state (state F) is of the form Fea34+═O2-···CuB2+, in which the H2O ligand that resulted from protonation of the OH- ligand in the PR state is dissociated from the CuB2+ site. The calculated Fea34+═O2- distance in F (1.68 Å) is 0.03 Å longer than that in PR (1.65 Å), which can explain the different Fea34+═O2- stretching modes in P (804 cm-1) and F (785 cm-1) identified by resonance Raman experiments. In this F state, the CuB2+···O2- (ferryl-oxygen) distance is only around 2.4 Å. Hence, the subsequent OH state [Fea33+-OH--CuB2+] with a μ-hydroxo bridge can be easily formed, as shown by our calculations.
Collapse
Affiliation(s)
- Wen-Ge Han Du
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute , 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego , 9500 Gilman Drive MC0505, La Jolla, California 92093, United States
| | - Louis Noodleman
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute , 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| |
Collapse
|
15
|
Schaefer AW, Kieber-Emmons MT, Adam SM, Karlin KD, Solomon EI. Phenol-Induced O-O Bond Cleavage in a Low-Spin Heme-Peroxo-Copper Complex: Implications for O 2 Reduction in Heme-Copper Oxidases. J Am Chem Soc 2017; 139:7958-7973. [PMID: 28521498 PMCID: PMC5605297 DOI: 10.1021/jacs.7b03292] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
This study evaluates the reaction of a biomimetic heme-peroxo-copper complex, {[(DCHIm)(F8)FeIII]-(O22-)-[CuII(AN)]}+ (1), with a phenolic substrate, involving a net H-atom abstraction to cleave the bridging peroxo O-O bond that produces FeIV═O, CuII-OH, and phenoxyl radical moieties, analogous to the chemistry carried out in heme-copper oxidases (HCOs). A 3D potential energy surface generated for this reaction reveals two possible reaction pathways: one involves nearly complete proton transfer (PT) from the phenol to the peroxo ligand before the barrier; the other involves O-O homolysis, where the phenol remains H-bonding to the peroxo OCu in the transition state (TS) and transfers the H+ after the barrier. In both mechanisms, electron transfer (ET) from phenol occurs after the PT (and after the barrier); therefore, only the interaction with the H+ is involved in lowering the O-O cleavage barrier. The relative barriers depend on covalency (which governs ET from Fe), and therefore vary with DFT functional. However, as these mechanisms differ by the amount of PT at the TS, kinetic isotope experiments were conducted to determine which mechanism is active. It is found that the phenolic proton exhibits a secondary kinetic isotope effect, consistent with the calculations for the H-bonded O-O homolysis mechanism. The consequences of these findings are discussed in relation to O-O cleavage in HCOs, supporting a model in which a peroxo intermediate serves as the active H+ acceptor, and both the H+ and e- required for O-O cleavage derive from the cross-linked Tyr residue present at the active site.
Collapse
Affiliation(s)
- Andrew W Schaefer
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| | - Matthew T Kieber-Emmons
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
- Department of Chemistry, The University of Utah , Salt Lake City, Utah 84112, United States
| | - Suzanne M Adam
- Department of Chemistry, The Johns Hopkins University , Baltimore, Maryland 21218, United States
| | - Kenneth D Karlin
- Department of Chemistry, The Johns Hopkins University , Baltimore, Maryland 21218, United States
| | - Edward I Solomon
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| |
Collapse
|
16
|
Carvalheda CA, Pisliakov AV. Insights into proton translocation in cbb 3 oxidase from MD simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2017; 1858:396-406. [PMID: 28259641 DOI: 10.1016/j.bbabio.2017.02.013] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Received: 10/05/2016] [Revised: 02/03/2017] [Accepted: 02/27/2017] [Indexed: 01/18/2023]
Abstract
Heme-copper oxidases are membrane protein complexes that catalyse the final step of the aerobic respiration, namely the reduction of oxygen to water. The energy released during catalysis is coupled to the active translocation of protons across the membrane, which contributes to the establishment of an electrochemical gradient that is used for ATP synthesis. The distinctive C-type (or cbb3) cytochrome c oxidases, which are mostly present in proteobacteria, exhibit a number of unique structural and functional features, including high catalytic activity at low oxygen concentrations. At the moment, the functioning mechanism of C-type oxidases, in particular the proton transfer/pumping mechanism presumably via a single proton channel, is still poorly understood. In this work we used all-atom molecular dynamics simulations and continuum electrostatics calculations to obtain atomic-level insights into the hydration and dynamics of a cbb3 oxidase. We provide the details of the water dynamics and proton transfer pathways for both the "chemical" and "pumped" protons, and show that formation of protonic connections is strongly affected by the protonation state of key residues, namely H243, E323 and H337.
Collapse
Affiliation(s)
- Catarina A Carvalheda
- Computational Biology, School of Life Sciences, University of Dundee, Dow Street, Dundee, DD1 5EH, United Kingdom; Physics, School of Sciences and Engineering, University of Dundee, Nethergate, Dundee, DD1 4HN, United Kingdom.
| | - Andrei V Pisliakov
- Computational Biology, School of Life Sciences, University of Dundee, Dow Street, Dundee, DD1 5EH, United Kingdom; Physics, School of Sciences and Engineering, University of Dundee, Nethergate, Dundee, DD1 4HN, United Kingdom.
| |
Collapse
|
17
|
Adam SM, Garcia-Bosch I, Schaefer AW, Sharma SK, Siegler MA, Solomon EI, Karlin KD. Critical Aspects of Heme-Peroxo-Cu Complex Structure and Nature of Proton Source Dictate Metal-O(peroxo) Breakage versus Reductive O-O Cleavage Chemistry. J Am Chem Soc 2017; 139:472-481. [PMID: 28029788 PMCID: PMC5274545 DOI: 10.1021/jacs.6b11322] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The 4H+/4e- reduction of O2 to water, a key fuel-cell reaction also carried out in biology by oxidase enzymes, includes the critical O-O bond reductive cleavage step. Mechanistic investigations on active-site model compounds, which are synthesized by rational design to incorporate systematic variations, can focus on and resolve answers to fundamental questions, including protonation and/or H-bonding aspects, which accompany electron transfer. Here, we describe the nature and comparative reactivity of two low-spin heme-peroxo-Cu complexes, LS-4DCHIm, [(DCHIm)F8FeIII-(O22-)-CuII(DCHIm)4]+, and LS-3DCHIm, [(DCHIm)F8FeIII-(O22-)-CuII(DCHIm)3]+ (F8 = tetrakis(2,6-difluorophenyl)-porphyrinate; DCHIm = 1,5-dicyclohexylimidazole), toward different proton (4-nitrophenol and [DMF·H+](CF3SO3-)) (DMF = dimethyl-formamide) or electron (decamethylferrocene (Fc*)) sources. Spectroscopic reactivity studies show that differences in structure and electronic properties of LS-3DCHIm and LS-4DCHIm lead to significant differences in behavior. LS-3DCHIm is resistant to reduction, is unreactive toward weakly acidic 4-NO2-phenol, and stronger acids cleave the metal-O bonds, releasing H2O2. By contrast, LS-4DCHIm forms an adduct with 4-NO2-phenol, which includes an H-bond to the peroxo O-atom distal to Fe (resonance Raman (rR) spectroscopy and DFT). With addition of Fc* (2 equiv overall required), O-O reductive cleavage occurs, giving water, Fe(III), and Cu(II) products; however, a kinetic study reveals a one-electron rate-determining process, ket = 1.6 M-1 s-1 (-90 °C). The intermediacy of a high-valent [(DCHIm)F8FeIV═O] species is thus implied, and separate experiments show that one-electron reduction-protonation of [(DCHIm)F8FeIV═O] occurs faster (ket2 = 5.0 M-1 s-1), consistent with the overall postulated mechanism. The importance of the H-bonding interaction as a prerequisite for reductive cleavage is highlighted.
Collapse
Affiliation(s)
- Suzanne M. Adam
- Johns Hopkins University, Baltimore, Maryland 21218, United States
| | | | - Andrew W. Schaefer
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Savita K. Sharma
- Johns Hopkins University, Baltimore, Maryland 21218, United States
| | | | - Edward I. Solomon
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | | |
Collapse
|
18
|
Yang L, Skjevik ÅA, Han Du WG, Noodleman L, Walker RC, Götz AW. Water exit pathways and proton pumping mechanism in B-type cytochrome c oxidase from molecular dynamics simulations. BIOCHIMICA ET BIOPHYSICA ACTA 2016; 1857:1594-1606. [PMID: 27317965 PMCID: PMC4995112 DOI: 10.1016/j.bbabio.2016.06.005] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2016] [Revised: 05/21/2016] [Accepted: 06/14/2016] [Indexed: 01/22/2023]
Abstract
Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. While proton uptake channels as well as water exit channels have been identified for A-type CcOs, the means by which water and protons exit B-type CcOs remain unclear. In this work, we investigate potential mechanisms for proton transport above the dinuclear center (DNC) in ba3-type CcO of Thermus thermophilus. Using long-time scale, all-atom molecular dynamics (MD) simulations for several relevant protonation states, we identify a potential mechanism for proton transport that involves propionate A of the active site heme a3 and residues Asp372, His376 and Glu126(II), with residue His376 acting as the proton-loading site. The proposed proton transport process involves a rotation of residue His376 and is in line with experimental findings. We also demonstrate how the strength of the salt bridge between residues Arg225 and Asp287 depends on the protonation state and that this salt bridge is unlikely to act as a simple electrostatic gate that prevents proton backflow. We identify two water exit pathways that connect the water pool above the DNC to the outer P-side of the membrane, which can potentially also act as proton exit transport pathways. Importantly, these water exit pathways can be blocked by narrowing the entrance channel between residues Gln151(II) and Arg449/Arg450 or by obstructing the entrance through a conformational change of residue Tyr136, respectively, both of which seem to be affected by protonation of residue His376.
Collapse
Affiliation(s)
- Longhua Yang
- Department of Chemistry, Nanchang University, 999 Xuefudadao, Nanchang, Jiangxi 330031, China; San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive MC0505, La Jolla, CA 92093, USA
| | - Åge A Skjevik
- San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive MC0505, La Jolla, CA 92093, USA; Department of Biomedicine, University of Bergen, N-5009 Bergen, Norway
| | - Wen-Ge Han Du
- Department of Integrative Structural and Computational Biology, GAC1118, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA
| | - Louis Noodleman
- Department of Integrative Structural and Computational Biology, GAC1118, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA
| | - Ross C Walker
- San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive MC0505, La Jolla, CA 92093, USA; Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA.
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive MC0505, La Jolla, CA 92093, USA.
| |
Collapse
|
19
|
Yang L, Skjevik ÅA, Han Du WG, Noodleman L, Walker RC, Götz AW. Data for molecular dynamics simulations of B-type cytochrome c oxidase with the Amber force field. Data Brief 2016; 8:1209-14. [PMID: 27547799 PMCID: PMC4979044 DOI: 10.1016/j.dib.2016.07.043] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2016] [Revised: 07/01/2016] [Accepted: 07/19/2016] [Indexed: 11/25/2022] Open
Abstract
Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. This article presents parameters for the cofactors of ba3-type CcO that are compatible with the all-atom Amber ff12SB and ff14SB force fields. Specifically, parameters were developed for the CuA pair, heme b, and the dinuclear center that consists of heme a3 and CuB bridged by a hydroperoxo group. The data includes geometries in XYZ coordinate format for cluster models that were employed to compute proton transfer energies and derive bond parameters and point charges for the force field using density functional theory. Also included are the final parameter files that can be employed with the Amber leap program to generate input files for molecular dynamics simulations with the Amber software package. Based on the high resolution (1.8 Å) X-ray crystal structure of the ba3-type CcO from Thermus thermophilus (Protein Data Bank ID number PDB: 3S8F), we built a model that is embedded in a POPC lipid bilayer membrane and solvated with TIP3P water molecules and counterions. We provide PDB data files of the initial model and the equilibrated model that can be used for further studies.
Collapse
Affiliation(s)
- Longhua Yang
- Department of Chemistry, Nanchang University, 999 Xuefudadao, Nanchang, Jiangxi 330031, China; San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA
| | - Åge A Skjevik
- San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA; Department of Biomedicine, University of Bergen, N-5009 Bergen, Norway
| | - Wen-Ge Han Du
- Department of Integrative Structural and Computational Biology, GAC1118, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA, USA
| | - Louis Noodleman
- Department of Integrative Structural and Computational Biology, GAC1118, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA, USA
| | - Ross C Walker
- San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA; Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA
| |
Collapse
|