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
|
Yang X, Liu S, Yin Z, Chen M, Song J, Li P, Yang L. New insights into the proton pumping mechanism of ba 3 cytochrome c oxidase: the functions of key residues and water. Phys Chem Chem Phys 2023; 25:25105-25115. [PMID: 37461851 DOI: 10.1039/d3cp01334k] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/28/2023]
Abstract
As the terminal oxidase of cell respiration in mitochondria and aerobic bacteria, the proton pumping mechanism of ba3-type cytochrome c oxidase (CcO) of Thermus thermophiles is still not fully understood. Especially, the functions of key residues which were considered as the possible proton loading sites (PLSs) above the catalytic center, as well as water located above and within the catalytic center, remain unclear. In this work, molecular dynamic simulations were performed on a set of designed mutants of key residues (Asp287, Asp372, His376, and Glu126II). The results showed that Asp287 may not be a PLS, but it could modulate the ability of the proton transfer pathway to transfer protons through its salt bridge with Arg225. Maintaining the closed state of the water pool above the catalytic center is necessary for the participation of inside water molecules in proton transfer. Water molecules inside the water pool can form hydrogen bond chains with PLS to facilitate proton transfer. Additional quantum cluster models of the Fe-Cu metal catalytic center are established, indicating that when the proton is transferred from Tyr237, it is more likely to reach the OCu atom directly through only one water molecule. This work provides a more profound understanding of the functions of important residues and specific water molecules in the proton pumping mechanism of CcO.
Collapse
Affiliation(s)
- Xiaoyue Yang
- School of Pharmaceutical Sciences & Key Laboratory of Advanced Drug Preparation Technologies, Zhengzhou University, Henan 450001, China.
| | - Shaohui Liu
- School of Pharmaceutical Sciences & Key Laboratory of Advanced Drug Preparation Technologies, Zhengzhou University, Henan 450001, China.
| | - Zhili Yin
- School of Pharmaceutical Sciences & Key Laboratory of Advanced Drug Preparation Technologies, Zhengzhou University, Henan 450001, China.
| | - Mengguo Chen
- School of Pharmaceutical Sciences & Key Laboratory of Advanced Drug Preparation Technologies, Zhengzhou University, Henan 450001, China.
| | - Jinshuai Song
- Green Catalysis Center, and College of Chemistry, Zhengzhou University, Henan 450001, China
| | - Pengfei Li
- Department of Chemistry and Biochemistry, Loyola University Chicago, Illinois 60660, USA
| | - Longhua Yang
- School of Pharmaceutical Sciences & Key Laboratory of Advanced Drug Preparation Technologies, Zhengzhou University, Henan 450001, China.
| |
Collapse
|
3
|
Fickenscher Z, Hey-Hawkins E. Added Complexity!-Mechanistic Aspects of Heterobimetallic Complexes for Application in Homogeneous Catalysis. Molecules 2023; 28:molecules28104233. [PMID: 37241974 DOI: 10.3390/molecules28104233] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2023] [Revised: 05/16/2023] [Accepted: 05/17/2023] [Indexed: 05/28/2023] Open
Abstract
Inspired by multimetallic assemblies and their role in enzyme catalysis, chemists have developed a plethora of heterobimetallic complexes for application in homogeneous catalysis. Starting with small heterobimetallic complexes with σ-donating and π-accepting ligands, such as N-heterocyclic carbene and carbonyl ligands, more and more complex systems have been developed over the past two decades. These systems can show a significant increase in catalytic activity compared with their monometallic counterparts. This increase can be attributed to new reaction pathways enabled by the presence of a second metal center in the active catalyst. This review focuses on mechanistic aspects of heterobimetallic complexes in homogeneous catalysis. Depending on the type of interaction of the second metal with the substrates, heterobimetallic complexes can be subdivided into four classes. Each of these classes is illustrated with multiple examples, showcasing the versatility of both, the types of interactions possible, and the reactions accessible.
Collapse
Affiliation(s)
- Zeno Fickenscher
- Institute of Inorganic Chemistry, Universität Leipzig, Johannisallee 29, D-04103 Leipzig, Germany
| | - Evamarie Hey-Hawkins
- Institute of Inorganic Chemistry, Universität Leipzig, Johannisallee 29, D-04103 Leipzig, Germany
| |
Collapse
|
4
|
Ghatak A, Samanta S, Nayek A, Mukherjee S, Dey SG, Dey A. Second-Sphere Hydrogen-Bond Donors and Acceptors Affect the Rate and Selectivity of Electrochemical Oxygen Reduction by Iron Porphyrins Differently. Inorg Chem 2022; 61:12931-12947. [PMID: 35939766 DOI: 10.1021/acs.inorgchem.2c02170] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The factors that control the rate and selectivity of 4e-/4H+ O2 reduction are important for efficient energy transformation as well as for understanding the terminal step of respiration in aerobic organisms. Inspired by the design of naturally occurring enzymes which are efficient catalysts for O2 and H2O2 reduction, several artificial systems have been generated where different second-sphere residues have been installed to enhance the rate and efficiency of the 4e-/4H+ O2 reduction. These include hydrogen-bonding residues like amines, carboxylates, ethers, amides, phenols, etc. In some cases, improvements in the catalysis were recorded, whereas in some cases improvements were marginal or nonexistent. In this work, we use an iron porphyrin complex with pendant 1,10-phenanthroline residues which show a pH-dependent variation of the rate of the electrochemical O2 reduction reaction (ORR) over 2 orders of magnitude. In-situ surface-enhanced resonance Raman spectroscopy reveals the presence of different intermediates at different pH's reflecting different rate-determining steps at different pH's. These data in conjunction with density functional theory calculations reveal that when the distal 1,10-phenanthroline is neutral it acts as a hydrogen-bond acceptor which stabilizes H2O (product) binding to the active FeII state and retards the reaction. However, when the 1,10-phenanthroline is protonated, it acts as a hydrogen-bond donor which enhances O2 reduction by stabilizing FeIII-O2.- and FeIII-OOH intermediates and activating the O-O bond for cleavage. On the basis of these data, general guidelines for controlling the different possible rate-determining steps in the complex multistep 4e-/4H+ ORR are developed and a bioinspired principle-based design of an efficient electrochemical ORR is presented.
Collapse
Affiliation(s)
- Arnab Ghatak
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A Raja SC Mullick Road, Kolkata, West Bengal 700032, India
| | - Soumya Samanta
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A Raja SC Mullick Road, Kolkata, West Bengal 700032, India
| | - Abhijit Nayek
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A Raja SC Mullick Road, Kolkata, West Bengal 700032, India
| | - Sudipta Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A Raja SC Mullick Road, Kolkata, West Bengal 700032, India
| | - Somdatta Ghosh Dey
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A Raja SC Mullick Road, Kolkata, West Bengal 700032, India
| | - Abhishek Dey
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A Raja SC Mullick Road, Kolkata, West Bengal 700032, India
| |
Collapse
|
5
|
Bhunia S, Ghatak A, Dey A. Second Sphere Effects on Oxygen Reduction and Peroxide Activation by Mononuclear Iron Porphyrins and Related Systems. Chem Rev 2022; 122:12370-12426. [PMID: 35404575 DOI: 10.1021/acs.chemrev.1c01021] [Citation(s) in RCA: 30] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Activation and reduction of O2 and H2O2 by synthetic and biosynthetic iron porphyrin models have proved to be a versatile platform for evaluating second-sphere effects deemed important in naturally occurring heme active sites. Advances in synthetic techniques have made it possible to install different functional groups around the porphyrin ligand, recreating artificial analogues of the proximal and distal sites encountered in the heme proteins. Using judicious choices of these substituents, several of the elegant second-sphere effects that are proposed to be important in the reactivity of key heme proteins have been evaluated under controlled environments, adding fundamental insight into the roles played by these weak interactions in nature. This review presents a detailed description of these efforts and how these have not only demystified these second-sphere effects but also how the knowledge obtained resulted in functional mimics of these heme enzymes.
Collapse
Affiliation(s)
- Sarmistha Bhunia
- School of Chemical Science, Indian Association for the Cultivation of Science, 2A Raja SC Mullick Road, Kolkata 700032, India
| | - Arnab Ghatak
- School of Chemical Science, Indian Association for the Cultivation of Science, 2A Raja SC Mullick Road, Kolkata 700032, India
| | - Abhishek Dey
- School of Chemical Science, Indian Association for the Cultivation of Science, 2A Raja SC Mullick Road, Kolkata 700032, India
| |
Collapse
|
6
|
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
|
7
|
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
|
8
|
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]
|
9
|
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
|
10
|
Han Du WG, McRee D, Götz AW, Noodleman L. A Water Molecule Residing in the Fe a33+···Cu B2+ Dinuclear Center of the Resting Oxidized as-Isolated Cytochrome c Oxidase: A Density Functional Study. Inorg Chem 2020; 59:8906-8915. [PMID: 32525689 PMCID: PMC8114904 DOI: 10.1021/acs.inorgchem.0c00724] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2020] [Indexed: 11/30/2022]
Abstract
Although the dinuclear center (DNC) of the resting oxidized "as-isolated" cytochrome c oxidase (CcO) is not a catalytically active state, its detailed structure, especially the nature of the bridging species between the Fea33+ and CuB2+ metal sites, is still both relevant and unsolved. Recent crystallographic work has shown an extended electron density for a peroxide type dioxygen species (O1-O2) bridging the Fea3 and CuB centers. In this paper, our density functional theory (DFT) calculations show that the observed peroxide type electron density between the two metal centers is most likely a mistaken analysis due to overlap of the electron density of a water molecule located at different positions between apparent O1 and O2 sites in DNCs of different CcO molecules with almost the same energy. Because the diffraction pattern and the resulting electron density map represent the effective long-range order averaged over many molecules and unit cells in the X-ray structure, this averaging can lead to an apparent observed superposition of different water positions between the Fea33+ and CuB2+ metal sites.
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
| | - Duncan McRee
- 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
|
11
|
Bím D, Chalupský J, Culka M, Solomon EI, Rulíšek L, Srnec M. Proton-Electron Transfer to the Active Site Is Essential for the Reaction Mechanism of Soluble Δ 9-Desaturase. J Am Chem Soc 2020; 142:10412-10423. [PMID: 32406236 DOI: 10.1021/jacs.0c01786] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
A full understanding of the catalytic action of non-heme iron (NHFe) and non-heme diiron (NHFe2) enzymes is still beyond the grasp of contemporary computational and experimental techniques. Many of these enzymes exhibit fascinating chemo-, regio-, and stereoselectivity, in spite of employing highly reactive intermediates which are necessary for activations of most stable chemical bonds. Herein, we study in detail one intriguing representative of the NHFe2 family of enzymes: soluble Δ9 desaturase (Δ9D), which desaturates rather than performing the thermodynamically favorable hydroxylation of substrate. Its catalytic mechanism has been explored in great detail by using QM(DFT)/MM and multireference wave function methods. Starting from the spectroscopically observed 1,2-μ-peroxo diferric P intermediate, the proton-electron uptake by P is the favored mechanism for catalytic activation, since it allows a significant reduction of the barrier of the initial (and rate-determining) H-atom abstraction from the stearoyl substrate as compared to the "proton-only activated" pathway. Also, we ruled out that a Q-like intermediate (high-valent diamond-core bis-μ-oxo-[FeIV]2 unit) is involved in the reaction mechanism. Our mechanistic picture is consistent with the experimental data available for Δ9D and satisfies fairly stringent conditions required by Nature: the chemo-, stereo-, and regioselectivity of the desaturation of stearic acid. Finally, the mechanisms evaluated are placed into a broader context of NHFe2 chemistry, provided by an amino acid sequence analysis through the families of the NHFe2 enzymes. Our study thus represents an important contribution toward understanding the catalytic action of the NHFe2 enzymes and may inspire further work in NHFe(2) biomimetic chemistry.
Collapse
Affiliation(s)
- Daniel Bím
- J. Heyrovský Institute of Physical Chemistry, The Czech Academy of Sciences, Dolejškova 3, Prague 8 182 23, Czech Republic.,Institute of Organic Chemistry and Biochemistry, The Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 166 10, Czech Republic
| | - Jakub Chalupský
- Institute of Organic Chemistry and Biochemistry, The Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 166 10, Czech Republic
| | - Martin Culka
- Institute of Organic Chemistry and Biochemistry, The Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 166 10, Czech Republic
| | - Edward I Solomon
- Department of Chemistry, Stanford University, 333 Campus Drive, Stanford, California 94305-5080, United States
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry, The Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 166 10, Czech Republic
| | - Martin Srnec
- J. Heyrovský Institute of Physical Chemistry, The Czech Academy of Sciences, Dolejškova 3, Prague 8 182 23, Czech Republic
| |
Collapse
|
12
|
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
|
13
|
Importance of the iron-sulfur component and of the siroheme modification in the resting state of sulfite reductase. J Inorg Biochem 2019; 203:110928. [PMID: 31756559 DOI: 10.1016/j.jinorgbio.2019.110928] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2019] [Revised: 11/12/2019] [Accepted: 11/13/2019] [Indexed: 11/24/2022]
Abstract
The active site of sulfite reductase (SiR) consists of an unusual siroheme-Fe4S4 assembly coupled via a cysteinate sulfur, and serves for multi-electron reduction reactions. Clear explanations have not been demonstrated for the reasons behind the choice of siroheme (vs. other types of heme) or for the single-atom coupling to an Fe4S4 center (as opposed to simple adjacency or to coupling via chains consisting of more than one atom). Possible explanations for these choices have previously been invoked, relating to the control of the spin state of the substrate-binding (siro)heme iron, modulation of the trans effect of the (Fe4S4-bound) cysteinate, or modulation of the redox potential. Reported here is a density functional theory (DFT) investigation of the structural interplay (in terms of geometry, molecular orbitals and magnetic interactions) between the siroheme and the Fe4S4 center as well as the importance of the covalent modifications within siroheme compared to the more common heme b, aiming to verify the role of the siroheme modification and of the Fe4S4 cluster at the SiR active site, with focus on previously-formulated hypotheses (geometrical/sterics, spin state, redox and electron-transfer control). A calibration of various DFT methods/variants for the correct description of ground state spin multiplicity is performed using a set of problematic cases of bioinorganic Fe centers; out of 11 functionals tested, M06-L and B3LYP offer the best results - though none of them correctly predict the spin state for all test cases. Upon examination of the relative energies of spin states, reduction potentials, energy decomposition (electrostatic, exchange-repulsion, orbital relaxation, correlation and dispersion interactions) and Mayer bond indices in SiR models, the following main roles of the siroheme and cubane are identified: (1) the cubane cofactor decreases the reduction potential of the siroheme and stabilizes the siroheme-cysteine bond interaction, and (2) the siroheme removes the quasi-degeneracy between the intermediate and high-spin states found in ferrous systems by preserving the latter as ground state; the higher-spin preference and the increased accessibility of multiple spin states are likely to be important in selective binding of the substrate and of the subsequent reaction intermediates, and in efficient changes in redox states throughout the catalytic cycle.
Collapse
|
14
|
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
|
15
|
Ehudin MA, Schaefer AW, Adam SM, Quist DA, Diaz DE, Tang JA, Solomon EI, Karlin KD. Influence of intramolecular secondary sphere hydrogen-bonding interactions on cytochrome c oxidase inspired low-spin heme-peroxo-copper complexes. Chem Sci 2019; 10:2893-2905. [PMID: 30996867 PMCID: PMC6431958 DOI: 10.1039/c8sc05165h] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2018] [Accepted: 01/03/2019] [Indexed: 11/21/2022] Open
Abstract
Dioxygen reduction by heme-copper oxidases is a critical biochemical process, wherein hydrogen bonding is hypothesized to participate in the critical step involving the active-site reductive cleavage of the O-O bond. Sixteen novel synthetic heme-(μ-O2 2-)-Cu(XTMPA) complexes, whose design is inspired by the cytochrome c oxidase active site structure, were generated in an attempt to form the first intramolecular H-bonded complexes. Derivatives of the "parent" ligand (XTMPA, TMPA = (tris((2-pyridyl)methyl)amine)) possessing one or two amine pendants preferentially form an H-bond with the copper-bound O-atom of the peroxide bridge. This is evidenced by a characteristic blue shift in the ligand-to-metal charge transfer (LMCT) bands observed in UV-vis spectroscopy (consistent with lowering of the peroxo π* relative to the iron orbitals) and a weakening of the O-O bond determined by resonance Raman spectroscopy (rR), with support from Density Functional Theory (DFT) calculations. Remarkably, with the TMPA-based infrastructure (versus similar heme-peroxo-copper complexes with different copper ligands), the typically undetected Cu-O stretch for these complexes was observed via rR, affording critical insights into the nature of the O-O peroxo core for the complexes studied. While amido functionalities have been shown to have greater H-bonding capabilities than their amino counterparts, in these heme-peroxo-copper complexes amido substituents distort the local geometry such that H-bonding with the peroxo core only imparts a weak electronic effect; optimal H-bonding interactions are observed by employing two amino groups on the copper ligand. The amino-substituted systems presented in this work reveal a key orientational anisotropy in H-bonding to the peroxo core for activating the O-O bond, offering critical insights into effective O-O cleavage chemistry. These findings indirectly support computational and protein structural studies suggesting the presence of an interstitial H-bonding water molecule in the CcO active site, which is critical for the desired reactivity. The results are evaluated with appropriate controls and discussed with respect to potential O2-reduction capabilities.
Collapse
Affiliation(s)
- Melanie A Ehudin
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , USA .
| | - Andrew W Schaefer
- Department of Chemistry , Stanford University , Stanford , California 94305 , USA .
| | - Suzanne M Adam
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , USA .
| | - David A Quist
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , USA .
| | - Daniel E Diaz
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , USA .
| | - Joel A Tang
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , USA .
| | - Edward I Solomon
- Department of Chemistry , Stanford University , Stanford , California 94305 , USA .
| | - Kenneth D Karlin
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , USA .
| |
Collapse
|
16
|
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
|
17
|
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
|
18
|
Hung MU, Yang ST, Ramanathan M, Liu ST. Hydrogenation of nitroarenes catalyzed by a dipalladium complex. Appl Organomet Chem 2017. [DOI: 10.1002/aoc.3976] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Ming-Uei Hung
- Department of Chemistry; National Taiwan University; Taipei 106 Taiwan
| | - Shu-Ting Yang
- Department of Chemistry; National Taiwan University; Taipei 106 Taiwan
| | - Mani Ramanathan
- Department of Chemistry; National Taiwan University; Taipei 106 Taiwan
| | - Shiuh-Tzung Liu
- Department of Chemistry; National Taiwan University; Taipei 106 Taiwan
| |
Collapse
|
19
|
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
|
20
|
Cui SX, Wang HY, Xu J, Zhang JP. Theoretical studies on magnetic properties of a binuclear paddle wheel Cu(II) complex {Cu2(μ2-O2CCH3)4}(OCNH2CH3)2. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2017. [DOI: 10.1134/s0036024417060085] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
21
|
Bhagi-Damodaran A, Michael MA, Zhu Q, Reed J, Sandoval BA, Mirts EN, Chakraborty S, Moënne-Loccoz P, Zhang Y, Lu Y. Why copper is preferred over iron for oxygen activation and reduction in haem-copper oxidases. Nat Chem 2017; 9:257-263. [PMID: 28221360 PMCID: PMC5321616 DOI: 10.1038/nchem.2643] [Citation(s) in RCA: 91] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2016] [Accepted: 09/09/2016] [Indexed: 01/13/2023]
Abstract
Haem-copper oxidase (HCO) catalyses the natural reduction of oxygen to water using a haem-copper centre. Despite decades of research on HCOs, the role of non-haem metal and the reason for nature's choice of copper over other metals such as iron remains unclear. Here, we use a biosynthetic model of HCO in myoglobin that selectively binds different non-haem metals to demonstrate 30-fold and 11-fold enhancements in the oxidase activity of Cu- and Fe-bound HCO mimics, respectively, as compared with Zn-bound mimics. Detailed electrochemical, kinetic and vibrational spectroscopic studies, in tandem with theoretical density functional theory calculations, demonstrate that the non-haem metal not only donates electrons to oxygen but also activates it for efficient O-O bond cleavage. Furthermore, the higher redox potential of copper and the enhanced weakening of the O-O bond from the higher electron density in the d orbital of copper are central to its higher oxidase activity over iron. This work resolves a long-standing question in bioenergetics, and renders a chemical-biological basis for the design of future oxygen-reduction catalysts.
Collapse
Affiliation(s)
| | - Matthew A. Michael
- Department of Biomedical Engineering, Chemistry, and
Biological Sciences, Stevens Institute of Technology, Hoboken, NJ, USA
| | - Qianhong Zhu
- Division of Environmental & Biomolecular Systems, Institute
of Environmental Health, Oregon Health & Science University, 3181 SW Sam Jackson Park
Road, Portland, OR, USA
| | - Julian Reed
- Department of Biochemistry, University of Illinois at
Urbana-Champaign, Urbana, IL, USA
| | - Braddock A. Sandoval
- Department of Chemistry, University of Illinois at
Urbana-Champaign, Urbana, IL, USA
| | - Evan N. Mirts
- Center for Biophysics and Quantitative Biology, University
of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Saumen Chakraborty
- Department of Chemistry, University of Illinois at
Urbana-Champaign, Urbana, IL, USA
| | - Pierre Moënne-Loccoz
- Division of Environmental & Biomolecular Systems, Institute
of Environmental Health, Oregon Health & Science University, 3181 SW Sam Jackson Park
Road, Portland, OR, USA
| | - Yong Zhang
- Department of Biomedical Engineering, Chemistry, and
Biological Sciences, Stevens Institute of Technology, Hoboken, NJ, USA
| | - Yi Lu
- Department of Chemistry, University of Illinois at
Urbana-Champaign, Urbana, IL, USA
- Department of Biochemistry, University of Illinois at
Urbana-Champaign, Urbana, IL, USA
- Center for Biophysics and Quantitative Biology, University
of Illinois at Urbana-Champaign, Urbana, IL, USA
| |
Collapse
|
22
|
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
|
23
|
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
|
24
|
Ohta T, Nagaraju P, Liu JG, Ogura T, Naruta Y. The secondary coordination sphere and axial ligand effects on oxygen reduction reaction by iron porphyrins: a DFT computational study. J Biol Inorg Chem 2016; 21:745-55. [PMID: 27501847 DOI: 10.1007/s00775-016-1380-9] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2016] [Accepted: 07/26/2016] [Indexed: 11/25/2022]
Abstract
Oxygen reduction reaction (ORR) catalyzed by a bio-inspired iron porphyrin bearing a hanging carboxylic acid group over the porphyrin ring, and a tethered axial imidazole ligand was studied by DFT calculations. BP86 free energy calculations of the redox potentials and pK a's of reaction components involved in the proton coupled electron transfer (PCET) reactions of the ferric-hydroxo and -superoxo complexes were performed based on Born-Haber thermodynamic cycle in conjunction with a continuum solvation model. The comparison was made with iron porphyrins that lack either in the hanging acid group or axial ligand, suggesting that H-bond interaction between the carboxylic acid and iron-bound hydroxo, aquo, superoxo, and peroxo ligands (de)stabilizes the Fe-O bonding, resulting in the increase in the reduction potential of the ferric complexes. The axial ligand interaction with the imidazole raises the affinity of the iron-bound superoxo and peroxo ligands for proton. In addition, a low-spin end-on ferric-hydroperoxo intermediate, a key precursor for O-O cleavage, can be stabilized in the presence of axial ligation. Thus, selective and efficient ORR of iron porphyrin can be achieved with the aid of the secondary coordination sphere and axial ligand interactions.
Collapse
Affiliation(s)
- Takehiro Ohta
- Picobiology Institute, Graduate School of Life Science, University of Hyogo, RSC-UH LP Center, Hyogo, 679-5148, Japan. .,Institute for Materials Chemistry and Engineering, Kyushu University, Fukuoka, 812-8581, Japan.
| | - Perumandla Nagaraju
- Institute of Science and Technology Research, Chubu University, Kasugai, Aichi, 487-8501, Japan.,Institute for Materials Chemistry and Engineering, Kyushu University, Fukuoka, 812-8581, Japan
| | - Jin-Gang Liu
- Department of Chemistry, East China University of Science and Technology, Shanghai, 200237, China
| | - Takashi Ogura
- Picobiology Institute, Graduate School of Life Science, University of Hyogo, RSC-UH LP Center, Hyogo, 679-5148, Japan
| | - Yoshinori Naruta
- Institute of Science and Technology Research, Chubu University, Kasugai, Aichi, 487-8501, Japan.,Institute for Materials Chemistry and Engineering, Kyushu University, Fukuoka, 812-8581, Japan
| |
Collapse
|
25
|
Liao RZ, Chen SL, Siegbahn PEM. Unraveling the Mechanism and Regioselectivity of the B12-Dependent Reductive Dehalogenase PceA. Chemistry 2016; 22:12391-9. [DOI: 10.1002/chem.201601575] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2016] [Indexed: 01/09/2023]
Affiliation(s)
- 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 P. R. China
| | - Shi-Lu Chen
- School of Chemistry; Beijing Institute of Technology; Beijing 100081 P. R. China
| | - Per E. M. Siegbahn
- Department of Organic Chemistry; Arrhenius Laboratory; Stockholm University; 10691 Stockholm Sweden
| |
Collapse
|
26
|
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
|
27
|
Huang DW, Liu YH, Peng SM, Liu ST. Dinickel complexes with anthyridine-based ligands. Dalton Trans 2016; 45:8265-71. [PMID: 27103235 DOI: 10.1039/c6dt00567e] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Two new dinickel complexes with anthyridine-based ligands, 5-phenyl-2,8-bis(2-pyridinyl)-1,9,10-anthyridine () and 5-phenyl-2,8-bis(6'-bipyridinyl)-1,9,10-anthyridine (), are reported. Complexation of Ni(OAc)2 with and in trifluoroacetic acid provided the corresponding dinickel complexes [{Ni2()(H2O)6(CF3COO)2}(CF3COO)2] () and [Ni2()(CF3COO)4(H2O)] (), respectively. Both complexes were characterized by spectroscopic methods and further confirmed by X-ray crystallography. Structural analyses of and revealed the NiNi distances in the complexes to be 5.4086(6) and 5.0138(7) Å, respectively. The catalytic activities of complexes and toward the reduction of carboxylic acids were evaluated. It appears that complex shows a good catalytic activity toward the reduction of carboxylates into the corresponding alcohols by diphenylsilane.
Collapse
Affiliation(s)
- Da-Wei Huang
- Department of Chemistry, National Taiwan University, Taipei 106, Taiwan.
| | - Yi-Hung Liu
- Department of Chemistry, National Taiwan University, Taipei 106, Taiwan.
| | - Shie-Ming Peng
- Department of Chemistry, National Taiwan University, Taipei 106, Taiwan.
| | - Shiuh-Tzung Liu
- Department of Chemistry, National Taiwan University, Taipei 106, Taiwan.
| |
Collapse
|
28
|
Han Du WG, Götz AW, Yang L, Walker RC, Noodleman L. 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. Phys Chem Chem Phys 2016; 18:21162-71. [PMID: 27094074 DOI: 10.1039/c6cp00349d] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
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.
Collapse
Affiliation(s)
- 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.
| | | | | | | | | |
Collapse
|
29
|
Huang DW, Liu YH, Peng SM, Liu ST. Dicopper Complexes with Anthyridine-Based Ligands: Coordination and Catalytic Activity. Organometallics 2016. [DOI: 10.1021/acs.organomet.5b00895] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Da-wei Huang
- Department of Chemistry, National Taiwan University, Taipei 106, Taiwan
| | - Yi-Hung Liu
- Department of Chemistry, National Taiwan University, Taipei 106, Taiwan
| | - Shie-Ming Peng
- Department of Chemistry, National Taiwan University, Taipei 106, Taiwan
| | - Shiuh-Tzung Liu
- Department of Chemistry, National Taiwan University, Taipei 106, Taiwan
| |
Collapse
|
30
|
Han Du WG, Noodleman L. Broken Symmetry DFT Calculations/Analysis for Oxidized and Reduced Dinuclear Center in Cytochrome c Oxidase: Relating Structures, Protonation States, Energies, and Mössbauer Properties in ba3 Thermus thermophilus. Inorg Chem 2015; 54:7272-90. [PMID: 26192749 PMCID: PMC4525772 DOI: 10.1021/acs.inorgchem.5b00700] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2015] [Indexed: 12/22/2022]
Abstract
The Fea3(3+)···CuB(2+) dinuclear center (DNC) structure of the as-isolated oxidized ba3 cytochrome c oxidase (CcO) from Thermus thermophilus (Tt) is still not fully understood. When the proteins are initially crystallized in the oxidized state, they typically become radiolyticly reduced through X-ray irradiation. Several X-ray crystal structures of reduced ba3 CcO from Tt are available. However, depending on whether the crystals were prepared in a lipidic cubic phase environment or in detergent micelles, and whether the CcO's were chemically or radiolyticly reduced, the X-ray diffraction analysis of the crystals showed different Fea3(2+)···CuB(+) DNC structures. On the other hand, Mössbauer spectroscopic experiments on reduced and oxidized ba3 CcOs from Tt (Zimmermann et al., Proc. Natl. Acad. Sci. USA 1988, 85, 5779-5783) revealed multiple (57)Fea3(2+) and (57)Fea3(3+) components. Moreover, one of the (57)Fea3(3+) components observed at 4.2 K transformed from a proposed "low-spin" state to a different high-spin species when the temperature was increased above 190 K, whereas the other high-spin (57)Fea3(3+) component remained unchanged. In the current Article, in order to understand the heterogeneities of the DNC in both Mössbauer spectra and X-ray crystal structures, the spin crossover of one of the (57)Fea3(3+) components, and how the coordination and spin states of the Fea3(3+/2+) and Cu(2+/1+) sites relate to the heterogeneity of the DNC structures, we have applied density functional OLYP calculations to the DNC clusters established based on the different X-ray crystal structures of ba3 CcO from Tt. As a result, specific oxidized and reduced DNC structures related to the observed Mössbauer spectra and to spectral changes with temperature have been proposed. Our calculations also show that, in certain intermediate states, the His233 and His283 ligand side chains may dissociate from the CuB(+) site, and they may become potential proton loading sites during the catalytic cycle.
Collapse
Affiliation(s)
- Wen-Ge Han Du
- Department of Integrative Structural and Computational
Biology, CB213, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| | - Louis Noodleman
- Department of Integrative Structural and Computational
Biology, CB213, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| |
Collapse
|
31
|
Proton transfer in the K-channel analog of B-type Cytochrome c oxidase from Thermus thermophilus. Biophys J 2015; 107:2177-84. [PMID: 25418102 DOI: 10.1016/j.bpj.2014.09.010] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2014] [Revised: 09/03/2014] [Accepted: 09/11/2014] [Indexed: 01/08/2023] Open
Abstract
A key enzyme in aerobic metabolism is cytochrome c oxidase (CcO), which catalyzes the reduction of molecular oxygen to water in the mitochondrial and bacterial membranes. Substrate electrons and protons are taken up from different sides of the membrane and protons are pumped across the membrane, thereby generating an electrochemical gradient. The well-studied A-type CcO uses two different entry channels for protons: the D-channel for all pumped and two consumed protons, and the K-channel for the other two consumed protons. In contrast, the B-type CcO uses only a single proton input channel for all consumed and pumped protons. It has the same location as the A-type K-channel (and thus is named the K-channel analog) without sharing any significant sequence homology. In this study, we performed molecular-dynamics simulations and electrostatic calculations to characterize the K-channel analog in terms of its energetic requirements and functionalities. The function of Glu-15B as a proton sink at the channel entrance is demonstrated by its rotational movement out of the channel when it is deprotonated and by its high pKA value when it points inside the channel. Tyr-244 in the middle of the channel is identified as the valve that ensures unidirectional proton transfer, as it moves inside the hydrogen-bond gap of the K-channel analog only while being deprotonated. The electrostatic energy landscape was calculated for all proton-transfer steps in the K-channel analog, which functions via proton-hole transfer. Overall, the K-channel analog has a very stable geometry without large energy barriers.
Collapse
|
32
|
Blachly PG, Sandala GM, Giammona DA, Bashford D, McCammon JA, Noodleman L. Broken-Symmetry DFT Computations for the Reaction Pathway of IspH, an Iron-Sulfur Enzyme in Pathogenic Bacteria. Inorg Chem 2015; 54:6439-61. [PMID: 26098647 PMCID: PMC4568833 DOI: 10.1021/acs.inorgchem.5b00751] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The recently discovered methylerythritol phosphate (MEP) pathway provides new targets for the development of antibacterial and antimalarial drugs. In the final step of the MEP pathway, the [4Fe-4S] IspH protein catalyzes the 2e(-)/2H(+) reductive dehydroxylation of (E)-4-hydroxy-3-methyl-but-2-enyl diphosphate (HMBPP) to afford the isoprenoid precursors isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP). Recent experiments have attempted to elucidate the IspH catalytic mechanism to drive inhibitor development. Two competing mechanisms have recently emerged, differentiated by their proposed HMBPP binding modes upon 1e(-) reduction of the [4Fe-4S] cluster: (1) a Birch reduction mechanism, in which HMBPP remains bound to the [4Fe-4S] cluster through its terminal C4-OH group (ROH-bound) until the -OH is cleaved as water; and (2) an organometallic mechanism, in which the C4-OH group rotates away from the [4Fe-4S] cluster, allowing the HMBPP olefin group to form a metallacycle complex with the apical iron (η(2)-bound). We perform broken-symmetry density functional theory computations to assess the energies and reduction potentials associated with the ROH- and η(2)-bound states implicated by these competing mechanisms. Reduction potentials obtained for ROH-bound states are more negative (-1.4 to -1.0 V) than what is typically expected of [4Fe-4S] ferredoxin proteins. Instead, we find that η(2)-bound states are lower in energy than ROH-bound states when the [4Fe-4S] cluster is 1e(-) reduced. Furthermore, η(2)-bound states can already be generated in the oxidized state, yielding reduction potentials of ca. -700 mV when electron addition occurs after rotation of the HMBPP C4-OH group. We demonstrate that such η(2)-bound states are kinetically accessible both when the IspH [4Fe-4S] cluster is oxidized and 1e(-) reduced. The energetically preferred pathway gives 1e(-) reduction of the cluster after substrate conformational change, generating the 1e(-) reduced intermediate proposed in the organometallic mechanism.
Collapse
Affiliation(s)
| | - Gregory M Sandala
- ‡Department of Chemistry and Biochemistry, Mount Allison University, 63C York Street, Sackville, New Brunswick E4L 1G8, Canada
| | - Debra Ann Giammona
- §Department of Structural Biology, St. Jude Children's Research Hospital, 262 Danny Thomas Place, Memphis, Tennessee 38105, United States
| | - Donald Bashford
- §Department of Structural Biology, St. Jude Children's Research Hospital, 262 Danny Thomas Place, Memphis, Tennessee 38105, United States
| | | | - Louis Noodleman
- #Department of Integrative Structural and Computational Biology, CB213, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| |
Collapse
|
33
|
Blachly PG, Sandala GM, Giammona D, Liu T, Bashford D, McCammon JA, Noodleman L. Use of Broken-Symmetry Density Functional Theory To Characterize the IspH Oxidized State: Implications for IspH Mechanism and Inhibition. J Chem Theory Comput 2014; 10:3871-3884. [PMID: 25221444 PMCID: PMC4159220 DOI: 10.1021/ct5005214] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2014] [Indexed: 12/31/2022]
Abstract
With current therapies becoming less efficacious due to increased drug resistance, new inhibitors of both bacterial and malarial targets are desperately needed. The recently discovered methylerythritol phosphate (MEP) pathway for isoprenoid synthesis provides novel targets for the development of such drugs. Particular attention has focused on the IspH protein, the final enzyme in the MEP pathway, which uses its [4Fe-4S] cluster to catalyze the formation of the isoprenoid precursors IPP and DMAPP from HMBPP. IspH catalysis is achieved via a 2e-/2H+ reductive dehydroxylation of HMBPP; the mechanism by which catalysis is achieved, however, is highly controversial. The work presented herein provides the first step in assessing different routes to catalysis by using computational methods. By performing broken-symmetry density functional theory (BS-DFT) calculations that employ both the conductor-like screening solvation model (DFT/COSMO) and a finite-difference Poisson-Boltzmann self-consistent reaction field methodology (DFT/SCRF), we evaluate geometries, energies, and Mössbauer signatures of the different protonation states that may exist in the oxidized state of the IspH catalytic cycle. From DFT/SCRF computations performed on the oxidized state, we find a state where the substrate, HMBPP, coordinates the apical iron in the [4Fe-4S] cluster as an alcohol group (ROH) to be one of two, isoenergetic, lowest-energy states. In this state, the HMBPP pyrophosphate moiety and an adjacent glutamate residue (E126) are both fully deprotonated, making the active site highly anionic. Our findings that this low-energy state also matches the experimental geometry of the active site and that its computed isomer shifts agree with experiment validate the use of the DFT/SCRF method to assess relative energies along the IspH reaction pathway. Additional studies of IspH catalytic intermediates are currently being pursued.
Collapse
Affiliation(s)
- Patrick G. Blachly
- Department
of Chemistry and Biochemistry, University
of California San Diego, 9500 Gilman Drive, Mail Code 0365, La Jolla, California 92093-0365, United States
| | - Gregory M. Sandala
- Department
of Chemistry and Biochemistry, Mount Allison
University, 63C York
Street, Sackville, New Brunswick E4L 1G8, Canada
| | - Debra
Ann Giammona
- Department
of Structural Biology, St. Jude Children’s
Research Hospital, 262
Danny Thomas Place, Memphis, Tennessee 38105, United States
| | - Tiqing Liu
- Skaggs School of Pharmacy and Pharmaceutical
Sciences, Howard Hughes Medical
Institute, and Department of Pharmacology, University
of California San Diego, La Jolla, California 92093-0365, United States
| | - Donald Bashford
- Department
of Structural Biology, St. Jude Children’s
Research Hospital, 262
Danny Thomas Place, Memphis, Tennessee 38105, United States
| | - J. Andrew McCammon
- Department
of Chemistry and Biochemistry, University
of California San Diego, 9500 Gilman Drive, Mail Code 0365, La Jolla, California 92093-0365, United States
- Skaggs School of Pharmacy and Pharmaceutical
Sciences, Howard Hughes Medical
Institute, and Department of Pharmacology, University
of California San Diego, La Jolla, California 92093-0365, United States
| | - Louis Noodleman
- Department
of Integrative Structural and Computational Biology, The Scripps Research Institute, TPC15, 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| |
Collapse
|
34
|
Gagliardi L, Solomon EI. Preface for the Forum on Insights into Spectroscopy and Reactivity from Electronic Structure Theory. Inorg Chem 2014; 53:6357-60. [DOI: 10.1021/ic5013654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Laura Gagliardi
- Department
of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Edward I. Solomon
- Department
of Chemistry, Stanford University, Stanford, California 94305-4401, United States
| |
Collapse
|