1
|
Sohraby F, Nunes-Alves A. Characterization of the Bottlenecks and Pathways for Inhibitor Dissociation from [NiFe] Hydrogenase. J Chem Inf Model 2024; 64:4193-4203. [PMID: 38728115 PMCID: PMC11134402 DOI: 10.1021/acs.jcim.4c00187] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Revised: 04/24/2024] [Accepted: 04/30/2024] [Indexed: 05/12/2024]
Abstract
[NiFe] hydrogenases can act as efficient catalysts for hydrogen oxidation and biofuel production. However, some [NiFe] hydrogenases are inhibited by gas molecules present in the environment, such as O2 and CO. One strategy to engineer [NiFe] hydrogenases and achieve O2- and CO-tolerant enzymes is by introducing point mutations to block the access of inhibitors to the catalytic site. In this work, we characterized the unbinding pathways of CO in the complex with the wild-type and 10 different mutants of [NiFe] hydrogenase from Desulfovibrio fructosovorans using τ-random accelerated molecular dynamics (τRAMD) to enhance the sampling of unbinding events. The ranking provided by the relative residence times computed with τRAMD is in agreement with experiments. Extensive data analysis of the simulations revealed that from the two bottlenecks proposed in previous studies for the transit of gas molecules (residues 74 and 122 and residues 74 and 476), only one of them (residues 74 and 122) effectively modulates diffusion and residence times for CO. We also computed pathway probabilities for the unbinding of CO, O2, and H2 from the wild-type [NiFe] hydrogenase, and we observed that while the most probable pathways are the same, the secondary pathways are different. We propose that introducing mutations to block the most probable paths, in combination with mutations to open the main secondary path used by H2, can be a feasible strategy to achieve CO and O2 resistance in the [NiFe] hydrogenase from Desulfovibrio fructosovorans.
Collapse
Affiliation(s)
- Farzin Sohraby
- Institute of Chemistry, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Ariane Nunes-Alves
- Institute of Chemistry, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| |
Collapse
|
2
|
Gomez A, Erb TJ, Grubmüller H, Vöhringer-Martinez E. Conformational Dynamics of the Most Efficient Carboxylase Contributes to Efficient CO 2 Fixation. J Chem Inf Model 2023; 63:7807-7815. [PMID: 38049384 DOI: 10.1021/acs.jcim.3c01447] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/06/2023]
Abstract
Crotonyl-CoA carboxylase/reductase (Ccr) is one of the fastest CO2 fixing enzymes and has become part of efficient artificial CO2-fixation pathways in vitro, paving the way for future applications. The underlying mechanism of its efficiency, however, is not yet completely understood. X-ray structures of different intermediates in the catalytic cycle reveal tetramers in a dimer of dimers configuration with two open and two closed active sites. Upon binding a substrate, this active site changes its conformation from the open state to the closed state. It is challenging to predict how these coupled conformational changes will alter the CO2 binding affinity to the reaction's active site. To determine whether the open or closed conformations of Ccr affect binding of CO2 to the active site, we performed all-atom molecular simulations of the various conformations of Ccr. The open conformation without a substrate showed the highest binding affinity. The CO2 binding sites are located near the catalytic relevant Asn81 and His365 residues and in an optimal position for CO2 fixation. Furthermore, they are unaffected by substrate binding, and CO2 molecules stay in these binding sites for a longer time. Longer times at these reactive binding sites facilitate CO2 fixation through the nucleophilic attack of the reactive enolate in the closed conformation. We previously demonstrated that the Asn81Leu variant cannot fix CO2. Simulations of the Asn81Leu variant explain the loss of activity through the removal of the Asn81 and His365 binding sites. Overall, our findings show that the conformational dynamics of the enzyme controls CO2 binding. Conformational changes in Ccr increase the level of CO2 in the open subunit before the substrate is bound, the active site closes, and the reaction starts. The full catalytic Ccr cycle alternates among CO2 addition, conformational change, and chemical reaction in the four subunits of the tetramer coordinated by communication between the two dimers.
Collapse
Affiliation(s)
- Aharon Gomez
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepión 4030000, Chile
| | - Tobias J Erb
- Department of Biochemistry and Synthetic Metabolism, Max Planck Institute for Terrestrial Microbiology, Karl-von-Frisch-Straße 10, Marburg D-35043, Germany
- LOEWE Center for Synthetic Microbiology (SYNMIKRO), Marburg 35032, Germany
| | - Helmut Grubmüller
- Department of Theoretical and Computational Biophysics, Max-Planck-Institute for Multidisciplinary Sciences, Am Fassberg 11, Göttingen D-37073, Germany
| | - Esteban Vöhringer-Martinez
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepión 4030000, Chile
| |
Collapse
|
3
|
Stripp ST, Duffus BR, Fourmond V, Léger C, Leimkühler S, Hirota S, Hu Y, Jasniewski A, Ogata H, Ribbe MW. Second and Outer Coordination Sphere Effects in Nitrogenase, Hydrogenase, Formate Dehydrogenase, and CO Dehydrogenase. Chem Rev 2022; 122:11900-11973. [PMID: 35849738 PMCID: PMC9549741 DOI: 10.1021/acs.chemrev.1c00914] [Citation(s) in RCA: 55] [Impact Index Per Article: 27.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Gases like H2, N2, CO2, and CO are increasingly recognized as critical feedstock in "green" energy conversion and as sources of nitrogen and carbon for the agricultural and chemical sectors. However, the industrial transformation of N2, CO2, and CO and the production of H2 require significant energy input, which renders processes like steam reforming and the Haber-Bosch reaction economically and environmentally unviable. Nature, on the other hand, performs similar tasks efficiently at ambient temperature and pressure, exploiting gas-processing metalloenzymes (GPMs) that bind low-valent metal cofactors based on iron, nickel, molybdenum, tungsten, and sulfur. Such systems are studied to understand the biocatalytic principles of gas conversion including N2 fixation by nitrogenase and H2 production by hydrogenase as well as CO2 and CO conversion by formate dehydrogenase, carbon monoxide dehydrogenase, and nitrogenase. In this review, we emphasize the importance of the cofactor/protein interface, discussing how second and outer coordination sphere effects determine, modulate, and optimize the catalytic activity of GPMs. These may comprise ionic interactions in the second coordination sphere that shape the electron density distribution across the cofactor, hydrogen bonding changes, and allosteric effects. In the outer coordination sphere, proton transfer and electron transfer are discussed, alongside the role of hydrophobic substrate channels and protein structural changes. Combining the information gained from structural biology, enzyme kinetics, and various spectroscopic techniques, we aim toward a comprehensive understanding of catalysis beyond the first coordination sphere.
Collapse
Affiliation(s)
- Sven T Stripp
- Freie Universität Berlin, Experimental Molecular Biophysics, Berlin 14195, Germany
| | | | - Vincent Fourmond
- Laboratoire de Bioénergétique et Ingénierie des Protéines, Institut de Microbiologie de la Méditerranée, Institut Microbiologie, Bioénergies et Biotechnologie, CNRS, Aix Marseille Université, Marseille 13402, France
| | - Christophe Léger
- Laboratoire de Bioénergétique et Ingénierie des Protéines, Institut de Microbiologie de la Méditerranée, Institut Microbiologie, Bioénergies et Biotechnologie, CNRS, Aix Marseille Université, Marseille 13402, France
| | - Silke Leimkühler
- University of Potsdam, Molecular Enzymology, Potsdam 14476, Germany
| | - Shun Hirota
- Nara Institute of Science and Technology, Division of Materials Science, Graduate School of Science and Technology, Nara 630-0192, Japan
| | - Yilin Hu
- Department of Molecular Biology & Biochemistry, University of California, Irvine, California 92697-3900, United States
| | - Andrew Jasniewski
- Department of Molecular Biology & Biochemistry, University of California, Irvine, California 92697-3900, United States
| | - Hideaki Ogata
- Nara Institute of Science and Technology, Division of Materials Science, Graduate School of Science and Technology, Nara 630-0192, Japan
- Hokkaido University, Institute of Low Temperature Science, Sapporo 060-0819, Japan
- Graduate School of Science, University of Hyogo, Hyogo 678-1297, Japan
| | - Markus W Ribbe
- Department of Molecular Biology & Biochemistry, University of California, Irvine, California 92697-3900, United States
- Department of Chemistry, University of California, Irvine, California 92697-2025, United States
| |
Collapse
|
4
|
Banerjee R, Lipscomb JD. Small-Molecule Tunnels in Metalloenzymes Viewed as Extensions of the Active Site. Acc Chem Res 2021; 54:2185-2195. [PMID: 33886257 DOI: 10.1021/acs.accounts.1c00058] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Rigorous substrate selectivity is a hallmark of enzyme catalysis. This selectivity is generally ascribed to a thermodynamically favorable process of substrate binding to the enzyme active site based upon complementary physiochemical characteristics, which allows both acquisition and orientation. However, this chemical selectivity is more difficult to rationalize for diminutive molecules that possess too narrow a range of physical characteristics to allow either precise positioning or discrimination between a substrate and an inhibitor. Foremost among these small molecules are dissolved gases such as H2, N2, O2, CO, CO2, NO, N2O, NH3, and CH4 so often encountered in metalloenzyme catalysis. Nevertheless, metalloenzymes have evolved to metabolize these small-molecule substrates with high selectivity and efficiency.The soluble methane monooxygenase enzyme (sMMO) acts upon two of these small molecules, O2 and CH4, to generate methanol as part of the C1 metabolic pathway of methanotrophic organisms. sMMO is capable of oxidizing many alternative hydrocarbon substrates. Remarkably, however, it will preferentially oxidize methane, the substrate with the fewest discriminating physical characteristics and the strongest C-H bond. Early studies led us to broadly attribute this specificity to the formation of a "molecular sieve" in which a methane- and oxygen-sized tunnel provides a size-selective route from bulk solvent to the completely buried sMMO active site. Indeed, recent cryogenic and serial femtosecond ambient temperature crystallographic studies have revealed such a route in sMMO. A detailed study of the sMMO tunnel considered here in the context of small-molecule tunnels identified in other metalloenzymes reveals three discrete characteristics that contribute to substrate selectivity and positioning beyond that which can be provided by the active site itself. Moreover, the dynamic nature of many tunnels allows an exquisite coordination of substrate binding and reaction phases of the catalytic cycle. Here we differentiate between the highly selective molecular tunnel, which allows only the one-dimensional transit of small molecules, and the larger, less-selective channels found in typical enzymes. Methods are described to identify and characterize tunnels as well as to differentiate them from channels. In metalloenzymes which metabolize dissolved gases, we posit that the contribution of tunnels is so great that they should be considered to be extensions of the active site itself. A full understanding of catalysis by these enzymes requires an appreciation of the roles played by tunnels. Such an understanding will also facilitate the use of the enzymes or their synthetic mimics in industrial or pharmaceutical applications.
Collapse
Affiliation(s)
- Rahul Banerjee
- Department of Biochemistry, Molecular Biology and Biophysics, University of Minnesota, Minneapolis, Minnesota 55391, United States
| | - John D. Lipscomb
- Department of Biochemistry, Molecular Biology and Biophysics, University of Minnesota, Minneapolis, Minnesota 55391, United States
| |
Collapse
|
5
|
Rovaletti A, Greco C, Ryde U. QM/MM study of the binding of H 2 to MoCu CO dehydrogenase: development and applications of improved H 2 van der Waals parameters. J Mol Model 2021; 27:68. [PMID: 33538901 PMCID: PMC7862525 DOI: 10.1007/s00894-020-04655-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Accepted: 12/20/2020] [Indexed: 11/28/2022]
Abstract
The MoCu CO dehydrogenase enzyme not only transforms CO into CO2 but it can also oxidise H2. Even if its hydrogenase activity has been known for decades, a debate is ongoing on the most plausible mode for the binding of H2 to the enzyme active site and the hydrogen oxidation mechanism. In the present work, we provide a new perspective on the MoCu-CODH hydrogenase activity by improving the in silico description of the enzyme. Energy refinement—by means of the BigQM approach—was performed on the intermediates involved in the dihydrogen oxidation catalysis reported in our previously published work (Rovaletti, et al. “Theoretical Insights into the Aerobic Hydrogenase Activity of Molybdenum–Copper CO Dehydrogenase.” Inorganics 7 (2019) 135). A suboptimal description of the H2–HN(backbone) interaction was observed when the van der Waals parameters described in previous literature for H2 were employed. Therefore, a new set of van der Waals parameters is developed here in order to better describe the hydrogen–backbone interaction. They give rise to improved binding modes of H2 in the active site of MoCu CO dehydrogenase. Implications of the resulting outcomes for a better understanding of hydrogen oxidation catalysis mechanisms are proposed and discussed.
Collapse
Affiliation(s)
- Anna Rovaletti
- Department of Earth and Environmental Sciences, Milano-Bicocca University, Piazza della Scienza 1, 20126, Milan, Italy
| | - Claudio Greco
- Department of Earth and Environmental Sciences, Milano-Bicocca University, Piazza della Scienza 1, 20126, Milan, Italy.
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00, Lund, Sweden.
| |
Collapse
|
6
|
Orio M, Pantazis DA. Successes, challenges, and opportunities for quantum chemistry in understanding metalloenzymes for solar fuels research. Chem Commun (Camb) 2021; 57:3952-3974. [DOI: 10.1039/d1cc00705j] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Overview of the rich and diverse contributions of quantum chemistry to understanding the structure and function of the biological archetypes for solar fuel research, photosystem II and hydrogenases.
Collapse
Affiliation(s)
- Maylis Orio
- Aix-Marseille Université
- CNRS
- iSm2
- Marseille
- France
| | - Dimitrios A. Pantazis
- Max-Planck-Institut für Kohlenforschung
- Kaiser-Wilhelm-Platz 1
- 45470 Mülheim an der Ruhr
- Germany
| |
Collapse
|
7
|
Zacarias S, Temporão A, Carpentier P, van der Linden P, Pereira IAC, Matias PM. Exploring the gas access routes in a [NiFeSe] hydrogenase using crystals pressurized with krypton and oxygen. J Biol Inorg Chem 2020; 25:863-874. [PMID: 32865640 DOI: 10.1007/s00775-020-01814-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2020] [Accepted: 08/20/2020] [Indexed: 11/25/2022]
Abstract
Hydrogenases are metalloenzymes that catalyse both H2 evolution and uptake. They are gas-processing enzymes with deeply buried active sites, so the gases diffuse through channels that connect the active site to the protein surface. The [NiFeSe] hydrogenases are a special class of hydrogenases containing a selenocysteine as a nickel ligand; they are more catalytically active and less O2-sensitive than standard [NiFe] hydrogenases. Characterisation of the channel system of hydrogenases is important to understand how the inhibitor oxygen reaches the active site to cause oxidative damage. To this end, crystals of Desulfovibrio vulgaris Hildenborough [NiFeSe] hydrogenase were pressurized with krypton and oxygen, and a method for tracking labile O2 molecules was developed, for mapping a hydrophobic channel system similar to that of the [NiFe] enzymes as the major route for gas diffusion.
Collapse
Affiliation(s)
- Sónia Zacarias
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157, Oeiras, Portugal
| | - Adriana Temporão
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157, Oeiras, Portugal
| | - Philippe Carpentier
- European Synchrotron Radiation Facility, Grenoble, France
- Institut de Recherche Interdisciplinaire de Grenoble (IRIG), Laboratoire Chimie et Biologie des Métaux (LCBM), Université Grenoble Alpes, CNRS, CEA, Grenoble, France
| | - Peter van der Linden
- Partnership for Soft Condensed Matter, European Synchrotron Radiation Facility, 71 Avenue des Martyrs, CS 40220, 38043, Grenoble, France
| | - Inês A C Pereira
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157, Oeiras, Portugal
| | - Pedro M Matias
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157, Oeiras, Portugal.
- iBET, Instituto de Biologia Experimental e Tecnológica, Apartado 12, 2780-901, Oeiras, Portugal.
| |
Collapse
|
8
|
Javanainen M, Vattulainen I, Monticelli L. Correction to "On Atomistic Models for Molecular Oxygen". J Phys Chem B 2020; 124:6943-6946. [PMID: 32786927 DOI: 10.1021/acs.jpcb.0c06376] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
9
|
Barbosa TM, Baltazar CSA, Cruz DR, Lousa D, Soares CM. Studying O 2 pathways in [NiFe]- and [NiFeSe]-hydrogenases. Sci Rep 2020; 10:10540. [PMID: 32601316 PMCID: PMC7324405 DOI: 10.1038/s41598-020-67494-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Accepted: 06/02/2020] [Indexed: 12/25/2022] Open
Abstract
Hydrogenases are efficient biocatalysts for H2 production and oxidation with various potential biotechnological applications.[NiFe]-class hydrogenases are highly active in both production and oxidation processes—albeit primarily biased to the latter—but suffer from being sensitive to O2.[NiFeSe] hydrogenases are a subclass of [NiFe] hydrogenases with, usually, an increased insensitivity to aerobic environments. In this study we aim to understand the structural causes of the low sensitivity of a [NiFeSe]-hydrogenase, when compared with a [NiFe] class enzyme, by studying the diffusion of O2. To unravel the differences between the two enzymes, we used computational methods comprising Molecular Dynamics simulations with explicit O2 and Implicit Ligand Sampling methodologies. With the latter, we were able to map the free energy landscapes for O2 permeation in both enzymes. We derived pathways from these energy landscapes and selected the kinetically more relevant ones with reactive flux analysis using transition path theory. These studies evidence the existence of quite different pathways in both enzymes and predict a lower permeation efficiency for O2 in the case of the [NiFeSe]-hydrogenase when compared with the [NiFe] enzyme. These differences can explain the experimentally observed lower inhibition by O2 on [NiFeSe]-hydrogenases, when compared with [NiFe]-hydrogenases. A comprehensive map of the residues lining the most important O2 pathways in both enzymes is also presented.
Collapse
Affiliation(s)
- Tiago M Barbosa
- ITQB NOVA, Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157, Oeiras, Portugal
| | - Carla S A Baltazar
- ITQB NOVA, Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157, Oeiras, Portugal
| | - Davide R Cruz
- ITQB NOVA, Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157, Oeiras, Portugal
| | - Diana Lousa
- ITQB NOVA, Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157, Oeiras, Portugal
| | - Cláudio M Soares
- ITQB NOVA, Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157, Oeiras, Portugal.
| |
Collapse
|
10
|
Rozza AM, Menyhárd DK, Oláh J. Gas Sensing by Bacterial H-NOX Proteins: An MD Study. Molecules 2020; 25:molecules25122882. [PMID: 32585836 PMCID: PMC7356049 DOI: 10.3390/molecules25122882] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2020] [Revised: 06/17/2020] [Accepted: 06/17/2020] [Indexed: 01/08/2023] Open
Abstract
Gas sensing is crucial for both prokaryotes and eukaryotes and is primarily performed by heme-based sensors, including H-NOX domains. These systems may provide a new, alternative mode for transporting gaseous molecules in higher organisms, but for the development of such systems, a detailed understanding of the ligand-binding properties is required. Here, we focused on ligand migration within the protein matrix: we performed molecular dynamics simulations on three bacterial (Ka, Ns and Cs) H-NOX proteins and studied the kinetics of CO, NO and O2 diffusion. We compared the response of the protein structure to the presence of ligands, diffusion rate constants, tunnel systems and storage pockets. We found that the rate constant for diffusion decreases in the O2 > NO > CO order in all proteins, and in the Ns > Ks > Cs order if single-gas is considered. Competition between gases seems to seriously influence the residential time of ligands spent in the distal pocket. The channel system is profoundly determined by the overall fold, but the sidechain pattern has a significant role in blocking certain channels by hydrophobic interactions between bulky groups, cation-π interactions or hydrogen bonding triads. The majority of storage pockets are determined by local sidechain composition, although certain functional cavities, such as the distal and proximal pockets are found in all systems. A major guideline for the design of gas transport systems is the need to chemically bind the gas molecule to the protein, possibly joining several proteins with several heme groups together.
Collapse
Affiliation(s)
- Ahmed M. Rozza
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology, Budapest Szent Gellért tér 4, H-1111 Budapest, Hungary;
- Department of Biotechnology, Faculty of Agriculture, Al-Azhar University, Cairo 11651, Egypt
| | - Dóra K. Menyhárd
- Laboratory of Structural Chemistry and Biology & MTA-ELTE Protein Modelling Research Group, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary
- Correspondence: (D.K.M.); (J.O.)
| | - Julianna Oláh
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology, Budapest Szent Gellért tér 4, H-1111 Budapest, Hungary;
- Correspondence: (D.K.M.); (J.O.)
| |
Collapse
|
11
|
Tai H, Hirota S. Mechanism and Application of the Catalytic Reaction of [NiFe] Hydrogenase: Recent Developments. Chembiochem 2020; 21:1573-1581. [DOI: 10.1002/cbic.202000058] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2020] [Revised: 02/25/2020] [Indexed: 01/28/2023]
Affiliation(s)
- Hulin Tai
- MOE Key Laboratory of Natural Resources of the Changbai Mountain and Functional MoleculesDepartment of ChemistryYanbian University Park Road 977 Yanji 133002 Jilin China
| | - Shun Hirota
- Division of Materials ScienceGraduate School of Science and TechnologyNara Institute of Science and Technology 8916-5 Takayama Ikoma Nara 630-0192 Japan
| |
Collapse
|
12
|
Domene C, Jorgensen C, Schofield CJ. Mechanism of Molecular Oxygen Diffusion in a Hypoxia-Sensing Prolyl Hydroxylase Using Multiscale Simulation. J Am Chem Soc 2020; 142:2253-2263. [DOI: 10.1021/jacs.9b09236] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Carmen Domene
- Chemistry Research Laboratory, Mansfield Road, University of Oxford, Oxford OX1 3TA, United Kingdom
- Department of Chemistry, Britannia House, King’s College London, 7 Trinity Street, London SE1 1DB, United Kingdom
- Department of Chemistry, University of Bath, Claverton Down Bath BA2 7AY, United Kingdom
| | - Christian Jorgensen
- Department of Chemistry, Britannia House, King’s College London, 7 Trinity Street, London SE1 1DB, United Kingdom
| | - Christopher J. Schofield
- Chemistry Research Laboratory, Mansfield Road, University of Oxford, Oxford OX1 3TA, United Kingdom
| |
Collapse
|
13
|
Carof A, Giannini S, Blumberger J. How to calculate charge mobility in molecular materials from surface hopping non-adiabatic molecular dynamics - beyond the hopping/band paradigm. Phys Chem Chem Phys 2019; 21:26368-26386. [PMID: 31793569 DOI: 10.1039/c9cp04770k] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Charge transport in high mobility organic semiconductors is in an intermediate regime between small polaron hopping and band transport limits. We have recently shown that surface hopping non-adiabatic molecular dynamics is a powerful method for prediction of charge transport mechanisms in organic materials and for near-quantitative prediction of charge mobilities at room temperature where the effects of nuclear zero-point motion and tunneling are still relatively small [S. Giannini et al., Nat. Commun., 2019, 10, 3843]. Here we assess and critically discuss the extensions to Tully's original method that have led to this success: (i) correction for missing electronic decoherence, (ii) detection of trivial crossings and (iii) removal of decoherence correction-induced spurious charge transfer. If any one of these corrections is not included, the charge mobility diverges with system size, each for different physical reasons. Yet if they are included, convergence with system size, detailed balance and good internal consistency are achieved.
Collapse
Affiliation(s)
- Antoine Carof
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK.
| | - Samuele Giannini
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK.
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK. and Institute for Advanced Study, Technische Universität München, Lichtenbergstrasse 2 a, D-85748 Garching, Germany
| |
Collapse
|
14
|
Bringas M, Lombardi LE, Luque FJ, Estrin DA, Capece L. Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations. Chemphyschem 2019; 20:2451-2460. [PMID: 31365183 DOI: 10.1002/cphc.201900589] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2019] [Revised: 07/30/2019] [Indexed: 11/07/2022]
Abstract
Computer simulation studies of the molecular basis for ligand migration in proteins allow the description of key events such as the transition between docking sites, displacement of existing ligands and solvent molecules, and open/closure of specific "gates", among others. In heme proteins, ligand migration from the solvent to the active site preludes the binding to the heme iron and triggers different functions. In this work, molecular dynamics simulations, a Markov State Model of migration and empirical kinetic equations are combined to study the migration of O2 and NO in two truncated hemoglobins of Mycobacterium tuberculosis (Mt-TrHbN and Mt-TrHbO). For Mt-TrHbN, we show that the difference in the association constant in the oxy and deoxy states relies mainly in the displacement of water molecules anchored in the distal cavity in the deoxy form. The results here provide a valuable approach to study ligand migration in globins.
Collapse
Affiliation(s)
- Mauro Bringas
- Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, C1428EGA, Buenos Aires, Argentina.,Instituto de Química Física de los Materiales, Medio Ambiente y Energía (INQUIMAE-CONICET), C1428EGA, Buenos Aires, Argentina
| | - Leandro E Lombardi
- Instituto de Cálculo, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires - CONICET, C1428EGA, Buenos Aires, Argentina
| | - F Javier Luque
- Department of Nutrition, Food Sciences and Gastronomy, Faculty of Pharmacy and Food Sciences, University of Barcelona, Campus Torribera, 08921, Santa Coloma de Gramenet, Spain.,Institute of Biomedicine (IBUB) and Institute of Theoretical and Computational Chemistry (IQTCUB), University of Barcelona, 08028, Barcelona, Spain
| | - Darío A Estrin
- Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, C1428EGA, Buenos Aires, Argentina.,Instituto de Química Física de los Materiales, Medio Ambiente y Energía (INQUIMAE-CONICET), C1428EGA, Buenos Aires, Argentina
| | - Luciana Capece
- Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, C1428EGA, Buenos Aires, Argentina.,Instituto de Química Física de los Materiales, Medio Ambiente y Energía (INQUIMAE-CONICET), C1428EGA, Buenos Aires, Argentina
| |
Collapse
|
15
|
Lu Y, Koo J. O 2 sensitivity and H 2 production activity of hydrogenases-A review. Biotechnol Bioeng 2019; 116:3124-3135. [PMID: 31403182 DOI: 10.1002/bit.27136] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2019] [Revised: 07/23/2019] [Accepted: 08/05/2019] [Indexed: 01/24/2023]
Abstract
Hydrogenases are metalloproteins capable of catalyzing the interconversion between molecular hydrogen and protons and electrons. The iron-sulfur clusters within the enzyme enable rapid relay of electrons which are either consumed or generated at the active site. Their unparalleled catalytic efficiency has attracted attention, especially for potential use in H2 production and/or fuel cell technologies. However, there are limitations to using hydrogenases, especially due to their high O2 sensitivity. The subclass, called [FeFe] hydrogenases, are particularly more vulnerable to O2 but proficient in H2 production. In this review, we provide an overview of mechanistic and protein engineering studies focused on understanding and enhancing O2 tolerance of the enzyme. The emphasis is on ongoing studies that attempt to overcome O2 sensitivity of the enzyme while it catalyzes H2 production in an aerobic environment. We also discuss pioneering attempts to utilize the enzyme in biological H2 production and other industrial processes, as well as our own perspective on future applications.
Collapse
Affiliation(s)
- Yuan Lu
- Department of Chemical Engineering, Tsinghua University, Beijing, China
| | - Jamin Koo
- Department of Chemical Engineering, Hongik University, Seoul, Republic of Korea
| |
Collapse
|
16
|
del Barrio M, Guendon C, Kpebe A, Baffert C, Fourmond V, Brugna M, Léger C. Valine-to-Cysteine Mutation Further Increases the Oxygen Tolerance of Escherichia coli NiFe Hydrogenase Hyd-1. ACS Catal 2019. [DOI: 10.1021/acscatal.9b00543] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Melisa del Barrio
- CNRS, Bioénergétique et Ingénierie des Protéines, Aix Marseille University, 13009 Marseille, France
| | - Chloé Guendon
- CNRS, Bioénergétique et Ingénierie des Protéines, Aix Marseille University, 13009 Marseille, France
| | - Arlette Kpebe
- CNRS, Bioénergétique et Ingénierie des Protéines, Aix Marseille University, 13009 Marseille, France
| | - Carole Baffert
- CNRS, Bioénergétique et Ingénierie des Protéines, Aix Marseille University, 13009 Marseille, France
| | - Vincent Fourmond
- CNRS, Bioénergétique et Ingénierie des Protéines, Aix Marseille University, 13009 Marseille, France
| | - Myriam Brugna
- CNRS, Bioénergétique et Ingénierie des Protéines, Aix Marseille University, 13009 Marseille, France
| | - Christophe Léger
- CNRS, Bioénergétique et Ingénierie des Protéines, Aix Marseille University, 13009 Marseille, France
| |
Collapse
|
17
|
Nunes-Alves A, Zuckerman DM, Arantes GM. Escape of a Small Molecule from Inside T4 Lysozyme by Multiple Pathways. Biophys J 2019. [PMID: 29539393 DOI: 10.1016/j.bpj.2018.01.014] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
The T4 lysozyme L99A mutant is often used as a model system to study small-molecule binding to proteins, but pathways for ligand entry and exit from the buried binding site and the associated protein conformational changes have not been fully resolved. Here, molecular dynamics simulations were employed to model benzene exit from its binding cavity using the weighted ensemble (WE) approach to enhance sampling of low-probability unbinding trajectories. Independent WE simulations revealed four pathways for benzene exit, which correspond to transient tunnels spontaneously formed in previous simulations of apo T4 lysozyme. Thus, benzene unbinding occurs through multiple pathways partially created by intrinsic protein structural fluctuations. Motions of several α-helices and side chains were involved in ligand escape from metastable microstates. WE simulations also provided preliminary estimates of rate constants for each exit pathway. These results complement previous works and provide a semiquantitative characterization of pathway heterogeneity for binding of small molecules to proteins.
Collapse
Affiliation(s)
- Ariane Nunes-Alves
- Department of Biochemistry, Instituto de Química, Universidade de São Paulo, São Paulo, Brazil
| | - Daniel M Zuckerman
- Department of Biomedical Engineering, School of Medicine, Oregon Health & Science University, Portland, Oregon.
| | | |
Collapse
|
18
|
Chen G, Xu W, Lu D, Wu J, Liu Z. Markov-state model for CO 2 binding with carbonic anhydrase under confinement. J Chem Phys 2018; 148:035101. [PMID: 29352785 DOI: 10.1063/1.5003298] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
Enzyme immobilization with a nanostructure material can enhance its stability and facilitate reusability. However, the apparent activity is often compromised due to additional diffusion barriers and complex interactions with the substrates and solvent molecules. The present study elucidates the effects of the surface hydrophobicity of nano-confinement on CO2 diffusion to the active site of human carbonic anhydrase II (CA), an enzyme that is able to catalyze CO2 hydration at extremely high turnover rates. Using the Markov-state model in combination with coarse-grained molecular dynamics simulations, we demonstrate that a hydrophobic cage increases CO2 local density but hinders its diffusion towards the active site of CA under confinement. By contrast, a hydrophilic cage hinders CO2 adsorption but promotes its binding with CA. An optimal surface hydrophobicity can be identified to maximize both the CO2 occupation probability and the diffusion rate. The simulation results offer insight into understanding enzyme performance under nano-confinement and help us to advance broader applications of CA for CO2 absorption and recovery.
Collapse
Affiliation(s)
- Gong Chen
- Ministry of Education Key Laboratory of Industrial Biocatalysis, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Weina Xu
- Ministry of Education Key Laboratory of Industrial Biocatalysis, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Diannan Lu
- Ministry of Education Key Laboratory of Industrial Biocatalysis, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Jianzhong Wu
- Department of Chemical and Environmental Engineering, University of California, Riverside, California 92521, USA
| | - Zheng Liu
- Ministry of Education Key Laboratory of Industrial Biocatalysis, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| |
Collapse
|
19
|
Chen G, Kong X, Lu D, Wu J, Liu Z. Kinetics of CO 2 diffusion in human carbonic anhydrase: a study using molecular dynamics simulations and the Markov-state model. Phys Chem Chem Phys 2018; 19:11690-11697. [PMID: 28435964 DOI: 10.1039/c7cp00887b] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Molecular dynamics (MD) simulations, in combination with the Markov-state model (MSM), were applied to probe CO2 diffusion from an aqueous solution into the active site of human carbonic anhydrase II (hCA-II), an enzyme useful for enhanced CO2 capture and utilization. The diffusion process in the hydrophobic pocket of hCA-II was illustrated in terms of a two-dimensional free-energy landscape. We found that CO2 diffusion in hCA-II is a rate-limiting step in the CO2 diffusion-binding-reaction process. The equilibrium distribution of CO2 shows its preferential accumulation within a hydrophobic domain in the protein core region. An analysis of the committors and reactive fluxes indicates that the main pathway for CO2 diffusion into the active site of hCA-II is through a binding pocket where residue Gln136 contributes to the maximal flux. The simulation results offer a new perspective on the CO2 hydration kinetics and useful insights toward the development of novel biochemical processes for more efficient CO2 sequestration and utilization.
Collapse
Affiliation(s)
- Gong Chen
- Key Laboratory of Industrial Biocatalysis, Ministry of Education, China.
| | | | | | | | | |
Collapse
|
20
|
Abstract
Redox enzymes, which catalyze reactions involving electron transfers in living organisms, are very promising components of biotechnological devices, and can be envisioned for sensing applications as well as for energy conversion. In this context, one of the most significant challenges is to achieve efficient direct electron transfer by tunneling between enzymes and conductive surfaces. Based on various examples of bioelectrochemical studies described in the recent literature, this review discusses the issue of enzyme immobilization at planar electrode interfaces. The fundamental importance of controlling enzyme orientation, how to obtain such orientation, and how it can be verified experimentally or by modeling are the three main directions explored. Since redox enzymes are sizable proteins with anisotropic properties, achieving their functional immobilization requires a specific and controlled orientation on the electrode surface. All the factors influenced by this orientation are described, ranging from electronic conductivity to efficiency of substrate supply. The specificities of the enzymatic molecule, surface properties, and dipole moment, which in turn influence the orientation, are introduced. Various ways of ensuring functional immobilization through tuning of both the enzyme and the electrode surface are then described. Finally, the review deals with analytical techniques that have enabled characterization and quantification of successful achievement of the desired orientation. The rich contributions of electrochemistry, spectroscopy (especially infrared spectroscopy), modeling, and microscopy are featured, along with their limitations.
Collapse
|
21
|
Kalms J, Schmidt A, Frielingsdorf S, Utesch T, Gotthard G, von Stetten D, van der Linden P, Royant A, Mroginski MA, Carpentier P, Lenz O, Scheerer P. Tracking the route of molecular oxygen in O 2-tolerant membrane-bound [NiFe] hydrogenase. Proc Natl Acad Sci U S A 2018; 115:E2229-E2237. [PMID: 29463722 PMCID: PMC5877991 DOI: 10.1073/pnas.1712267115] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
[NiFe] hydrogenases catalyze the reversible splitting of H2 into protons and electrons at a deeply buried active site. The catalytic center can be accessed by gas molecules through a hydrophobic tunnel network. While most [NiFe] hydrogenases are inactivated by O2, a small subgroup, including the membrane-bound [NiFe] hydrogenase (MBH) of Ralstonia eutropha, is able to overcome aerobic inactivation by catalytic reduction of O2 to water. This O2 tolerance relies on a special [4Fe3S] cluster that is capable of releasing two electrons upon O2 attack. Here, the O2 accessibility of the MBH gas tunnel network has been probed experimentally using a "soak-and-freeze" derivatization method, accompanied by protein X-ray crystallography and computational studies. This combined approach revealed several sites of O2 molecules within a hydrophobic tunnel network leading, via two tunnel entrances, to the catalytic center of MBH. The corresponding site occupancies were related to the O2 concentrations used for MBH crystal derivatization. The examination of the O2-derivatized data furthermore uncovered two unexpected structural alterations at the [4Fe3S] cluster, which might be related to the O2 tolerance of the enzyme.
Collapse
Affiliation(s)
- Jacqueline Kalms
- Charité - Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin, and Berlin Institute of Health, Institute of Medical Physics and Biophysics, Group Protein X-ray Crystallography and Signal Transduction, D-10117 Berlin, Germany
| | - Andrea Schmidt
- Charité - Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin, and Berlin Institute of Health, Institute of Medical Physics and Biophysics, Group Protein X-ray Crystallography and Signal Transduction, D-10117 Berlin, Germany
| | | | - Tillmann Utesch
- Institut für Chemie, Technische Universität Berlin, D-10623 Berlin, Germany
| | | | | | - Peter van der Linden
- European Synchrotron Radiation Facility, F-38043 Grenoble, France
- Partnership for Soft Condensed Matter (PSCM), F-38043 Grenoble, France
| | - Antoine Royant
- European Synchrotron Radiation Facility, F-38043 Grenoble, France
- Univ. Grenoble Alpes, CNRS, CEA, Institut de Biologie Structurale (IBS), F-38000 Grenoble, France
| | | | - Philippe Carpentier
- European Synchrotron Radiation Facility, F-38043 Grenoble, France
- Univ. Grenoble Alpes, CNRS, CEA, Institut de Biosciences et Biotechnologies de Grenoble (BIG)-Laboratoire Chimie et Biologie des Métaux (LCBM), F-38000 Grenoble, France
| | - Oliver Lenz
- Institut für Chemie, Technische Universität Berlin, D-10623 Berlin, Germany
| | - Patrick Scheerer
- Charité - Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin, and Berlin Institute of Health, Institute of Medical Physics and Biophysics, Group Protein X-ray Crystallography and Signal Transduction, D-10117 Berlin, Germany;
| |
Collapse
|
22
|
Lábas A, Menyhárd DK, Harvey JN, Oláh J. First Principles Calculation of the Reaction Rates for Ligand Binding to Myoglobin: The Cases of NO and CO. Chemistry 2018; 24:5350-5358. [PMID: 29285802 DOI: 10.1002/chem.201704867] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2017] [Indexed: 12/12/2022]
Abstract
Ligand binding by proteins is among the most fundamental processes in nature. Among these processes the binding of small gas molecules, such as O2 , CO and NO to heme proteins has traditionally received vivid interest, which was further boosted by their recently recognized significant role in gas sensing in the body. At the heart of the binding of these ligands to the heme group is the spinforbidden reaction between high-spin iron(II) and the ligand yielding a low-spin adduct. We use computational means to address the complete mechanism of CO and NO binding by myoglobin. Considering that it involves several steps occurring on different time scales, molecular dynamics simulations were performed to address the diffusion of the ligand through the enzyme, and DFT calculations in combination with statistical rate calculation to investigate the spin-forbidden reaction. The calculations yielded rate constants in qualitative agreement with experiments and revealed that the bottleneck of NO and CO binding is different; for NO, diffusion was found to be rate-limiting, whereas for CO, the spin-forbidden step is the slowest.
Collapse
Affiliation(s)
- Anikó Lábas
- Department of Inorganic Chemistry, Budapest University of Technology and Economics, H-1111, Budapest, Szent Gellért tér 4., Hungary
| | - Dóra K Menyhárd
- MTA-ELTE Protein Modelling Research Group, H-1117, Budapest, Pázmány Péter st. 1/A, Hungary
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven, B-3001, Leuven Celestijnenlaan 200F- box 2404, Belgium
| | - Julianna Oláh
- Department of Inorganic Chemistry, Budapest University of Technology and Economics, H-1111, Budapest, Szent Gellért tér 4., Hungary
| |
Collapse
|
23
|
Quantum chemical approaches to [NiFe] hydrogenase. Essays Biochem 2017; 61:293-303. [PMID: 28487405 DOI: 10.1042/ebc20160079] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2016] [Revised: 02/22/2017] [Accepted: 03/01/2017] [Indexed: 11/17/2022]
Abstract
The mechanism by which [NiFe] hydrogenase catalyses the oxidation of molecular hydrogen is a significant yet challenging topic in bioinorganic chemistry. With far-reaching applications in renewable energy and carbon mitigation, significant effort has been invested in the study of these complexes. In particular, computational approaches offer a unique perspective on how this enzyme functions at an electronic and atomistic level. In this article, we discuss state-of-the art quantum chemical methods and how they have helped deepen our comprehension of [NiFe] hydrogenase. We outline the key strategies that can be used to compute the (i) geometry, (ii) electronic structure, (iii) thermodynamics and (iv) kinetic properties associated with the enzymatic activity of [NiFe] hydrogenase and other bioinorganic complexes.
Collapse
|
24
|
Wang PH, Yu I, Feig M, Sugita Y. Influence of protein crowder size on hydration structure and dynamics in macromolecular crowding. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.01.012] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
|
25
|
Javanainen M, Vattulainen I, Monticelli L. On Atomistic Models for Molecular Oxygen. J Phys Chem B 2017; 121:518-528. [PMID: 28004930 DOI: 10.1021/acs.jpcb.6b11183] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Molecular oxygen (O2) is key to all life on earth, as it is constantly cycled via photosynthesis and cellular respiration. Substantial scientific effort has been devoted to understanding every part of this cycle. Classical molecular dynamics (MD) simulations have been used to study some of the key processes involved in cellular respiration: O2 permeation through alveolar monolayers and cellular membranes, its binding to hemoglobin during transport in the bloodstream, as well as its transport along optimal pathways toward its reduction sites in proteins. Moreover, MD simulations can help interpret the results of several imaging techniques in which O2 is used because of its paramagnetic nature. However, despite the widespread use of computational models for the O2 molecule, their performances have never been systematically evaluated. In this paper, we assess the performances of 14 different models of O2 available in the literature by calculating four thermodynamic properties: density, heat of vaporization, free energy of hydration, and free energy of solvation in hexadecane. For each property, reliable experimental data are available. Most models perform reasonably well in predicting the correct trends, but they fail to reproduce the experimental data quantitatively. We then develop new models for O2, with and without a quadrupole moment, and compare their behavior with the behavior of previously published models. The new models show significant improvement in terms of density, heat of vaporization, and free energy of hydration. However, quantitative agreement with water-oil partitioning is not reached due to discrepancies between the calculated and measured free energies of solvation in hexadecane. We suggest that classical pairwise-additive models may be inadequate to properly describe the thermodynamics of solvation of apolar species, such as O2, in apolar solvents.
Collapse
Affiliation(s)
- Matti Javanainen
- Department of Physics, Tampere University of Technology , 33720 Tampere, Finland.,Department of Physics, University of Helsinki , 00100 Helsinki, Finland
| | - Ilpo Vattulainen
- Department of Physics, Tampere University of Technology , 33720 Tampere, Finland.,Department of Physics, University of Helsinki , 00100 Helsinki, Finland.,MEMPHYS - Centre for Biomembrane Physics, University of Southern Denmark , 5230 Odense, Denmark
| | - Luca Monticelli
- University of Lyon, CNRS, UMR 5086 MMSB , 69367 Lyon, France
| |
Collapse
|
26
|
Abstract
Markov state models (MSMs) are an extremely useful tool for understanding the conformational dynamics of macromolecules and for analyzing MD simulations in a quantitative fashion. They have been extensively used for peptide and protein folding, for small molecule binding, and for the study of native ensemble dynamics. Here, we adapt the MSM methodology to gain insight into the dynamics of misfolded states. To overcome possible flaws in root-mean-square deviation (RMSD)-based metrics, we introduce a novel discretization approach, based on coarse-grained contact maps. In addition, we extend the MSM methodology to include "sink" states in order to account for the irreversibility (on simulation time scales) of processes like protein misfolding. We apply this method to analyze the mechanism of misfolding of tandem repeats of titin domains, and how it is influenced by confinement in a chaperonin-like cavity.
Collapse
Affiliation(s)
- Anshul Sirur
- Department of Chemistry, University of Cambridge, Lensfield Road, CB2 1EW Cambridge, United Kingdom
| | - David De Sancho
- Department of Chemistry, University of Cambridge, Lensfield Road, CB2 1EW Cambridge, United Kingdom
| | - Robert B Best
- Department of Chemistry, University of Cambridge, Lensfield Road, CB2 1EW Cambridge, United Kingdom
| |
Collapse
|
27
|
Mayne CG, Arcario MJ, Mahinthichaichan P, Baylon JL, Vermaas JV, Navidpour L, Wen PC, Thangapandian S, Tajkhorshid E. The cellular membrane as a mediator for small molecule interaction with membrane proteins. BIOCHIMICA ET BIOPHYSICA ACTA 2016; 1858:2290-2304. [PMID: 27163493 PMCID: PMC4983535 DOI: 10.1016/j.bbamem.2016.04.016] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/18/2016] [Revised: 04/26/2016] [Accepted: 04/27/2016] [Indexed: 01/05/2023]
Abstract
The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.
Collapse
Affiliation(s)
- Christopher G Mayne
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States.
| | - Mark J Arcario
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, United States; College of Medicine, University of Illinois at Urbana-Champaign, United States.
| | - Paween Mahinthichaichan
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Department of Biochemistry, University of Illinois at Urbana-Champaign, United States.
| | - Javier L Baylon
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, United States.
| | - Josh V Vermaas
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, United States.
| | - Latifeh Navidpour
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States.
| | - Po-Chao Wen
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States.
| | - Sundarapandian Thangapandian
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Department of Biochemistry, University of Illinois at Urbana-Champaign, United States.
| | - Emad Tajkhorshid
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Department of Biochemistry, University of Illinois at Urbana-Champaign, United States; Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, United States; College of Medicine, University of Illinois at Urbana-Champaign, United States.
| |
Collapse
|
28
|
Kubas A, Orain C, De Sancho D, Saujet L, Sensi M, Gauquelin C, Meynial-Salles I, Soucaille P, Bottin H, Baffert C, Fourmond V, Best RB, Blumberger J, Léger C. Mechanism of O 2 diffusion and reduction in FeFe hydrogenases. Nat Chem 2016; 9:88-95. [PMID: 27995927 DOI: 10.1038/nchem.2592] [Citation(s) in RCA: 84] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2015] [Accepted: 07/08/2016] [Indexed: 11/09/2022]
Abstract
FeFe hydrogenases are the most efficient H2-producing enzymes. However, inactivation by O2 remains an obstacle that prevents them being used in many biotechnological devices. Here, we combine electrochemistry, site-directed mutagenesis, molecular dynamics and quantum chemical calculations to uncover the molecular mechanism of O2 diffusion within the enzyme and its reactions at the active site. We propose that the partial reversibility of the reaction with O2 results from the four-electron reduction of O2 to water. The third electron/proton transfer step is the bottleneck for water production, competing with formation of a highly reactive OH radical and hydroxylated cysteine. The rapid delivery of electrons and protons to the active site is therefore crucial to prevent the accumulation of these aggressive species during prolonged O2 exposure. These findings should provide important clues for the design of hydrogenase mutants with increased resistance to oxidative damage.
Collapse
Affiliation(s)
- Adam Kubas
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK.,Institute of Physical Chemistry, Polish Academy of Science, ul. Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Christophe Orain
- Aix Marseille Univ, CNRS, Laboratoire de Bioénergétique et Ingénierie des Protéines, Institut de Microbiologie de la Méditerranée, Marseille, France
| | - David De Sancho
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.,CIC nanoGUNE, Tolosa Hiribidea 76, 20018 Donostia-San Sebastián, Spain.,IKERBASQUE; Basque Foundation for Science, María Díaz de Haro 3, 48013 Bilbao, Spain
| | - Laure Saujet
- Institut de Biologie et de Technologies de Saclay IBITECS, SB2SM, F-91191 Gif sur Yvette, France.,Institut de Biologie Intégrative de la Cellule I2BC, UMR 9198, CEA, CNRS, Université Paris Sud, F-91191 Gif sur Yvette, France
| | - Matteo Sensi
- Aix Marseille Univ, CNRS, Laboratoire de Bioénergétique et Ingénierie des Protéines, Institut de Microbiologie de la Méditerranée, Marseille, France
| | - Charles Gauquelin
- Université de Toulouse, INSA, UPS, INP, LISBP, INRA:UMR792, CNRS:UMR 5504, 135 avenue de Rangueil, Toulouse 31077 Cedex 04, France
| | - Isabelle Meynial-Salles
- Université de Toulouse, INSA, UPS, INP, LISBP, INRA:UMR792, CNRS:UMR 5504, 135 avenue de Rangueil, Toulouse 31077 Cedex 04, France
| | - Philippe Soucaille
- Université de Toulouse, INSA, UPS, INP, LISBP, INRA:UMR792, CNRS:UMR 5504, 135 avenue de Rangueil, Toulouse 31077 Cedex 04, France
| | - Hervé Bottin
- Institut de Biologie et de Technologies de Saclay IBITECS, SB2SM, F-91191 Gif sur Yvette, France.,Institut de Biologie Intégrative de la Cellule I2BC, UMR 9198, CEA, CNRS, Université Paris Sud, F-91191 Gif sur Yvette, France
| | - Carole Baffert
- Aix Marseille Univ, CNRS, Laboratoire de Bioénergétique et Ingénierie des Protéines, Institut de Microbiologie de la Méditerranée, Marseille, France
| | - Vincent Fourmond
- Aix Marseille Univ, CNRS, Laboratoire de Bioénergétique et Ingénierie des Protéines, Institut de Microbiologie de la Méditerranée, Marseille, France
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892-0520, USA
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK
| | - Christophe Léger
- Aix Marseille Univ, CNRS, Laboratoire de Bioénergétique et Ingénierie des Protéines, Institut de Microbiologie de la Méditerranée, Marseille, France
| |
Collapse
|
29
|
Ogata H, Lubitz W, Higuchi Y. Structure and function of [NiFe] hydrogenases. J Biochem 2016; 160:251-258. [PMID: 27493211 DOI: 10.1093/jb/mvw048] [Citation(s) in RCA: 68] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2016] [Accepted: 07/06/2016] [Indexed: 12/22/2022] Open
Abstract
Hydrogenases catalyze the reversible conversion of molecular hydrogen to protons and electrons via a heterolytic splitting mechanism. The active sites of [NiFe] hydrogenases comprise a dinuclear Ni-Fe center carrying CO and CN- ligands. The catalytic activity of the standard (O2-sensitive) [NiFe] hydrogenases vanishes under aerobic conditions. The O2-tolerant [NiFe] hydrogenases can sustain H2 oxidation activity under atmospheric conditions. These hydrogenases have very similar active site structures that change the ligand sphere during the activation/catalytic process. An important structural difference between these hydrogenases has been found for the proximal iron-sulphur cluster located in the vicinity of the active site. This unprecedented [4Fe-3S]-6Cys cluster can supply two electrons, which lead to rapid recovery of the O2 inactivation, to the [NiFe] active site.
Collapse
Affiliation(s)
- Hideaki Ogata
- Max Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, Mülheim an der Ruhr 45470, Germany
| | - Wolfgang Lubitz
- Max Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, Mülheim an der Ruhr 45470, Germany
| | - Yoshiki Higuchi
- Department of Life Science, Graduate School of Life Science, University of Hyogo, 3-2-1 Koto, Kamigori-cho, Ako-gun, Hyogo 678-1297, Japan .,RIKEN SPring-8 Center, 1-1-1 Koto, Sayo-cho, Sayo-gun, Hyogo 679-5148, Japan.,CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
30
|
Vermaas JV, Trebesch N, Mayne CG, Thangapandian S, Shekhar M, Mahinthichaichan P, Baylon JL, Jiang T, Wang Y, Muller MP, Shinn E, Zhao Z, Wen PC, Tajkhorshid E. Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation. Methods Enzymol 2016; 578:373-428. [PMID: 27497175 PMCID: PMC6404235 DOI: 10.1016/bs.mie.2016.05.042] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.
Collapse
Affiliation(s)
- J V Vermaas
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - N Trebesch
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - C G Mayne
- University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - S Thangapandian
- University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - M Shekhar
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - P Mahinthichaichan
- University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - J L Baylon
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - T Jiang
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - Y Wang
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - M P Muller
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; College of Medicine, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - E Shinn
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - Z Zhao
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - P-C Wen
- University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - E Tajkhorshid
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; College of Medicine, University of Illinois at Urbana-Champaign, Urbana, IL, United States.
| |
Collapse
|
31
|
Mahinthichaichan P, Gennis RB, Tajkhorshid E. All the O2 Consumed by Thermus thermophilus Cytochrome ba3 Is Delivered to the Active Site through a Long, Open Hydrophobic Tunnel with Entrances within the Lipid Bilayer. Biochemistry 2016; 55:1265-78. [PMID: 26845082 DOI: 10.1021/acs.biochem.5b01255] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Cytochrome ba3 is a proton-pumping heme-copper oxygen reductase from the extreme thermophile Thermus thermophilus. Despite the fact that the enzyme's active site is buried deep within the protein, the apparent second order rate constant for the initial binding of O2 to the active-site heme has been experimentally found to be 10(9) M(-1) s(-1) at 298 K, at or near the diffusion limit, and 2 orders of magnitude faster than for O2 binding to myoglobin. To provide quantitative and microscopic descriptions of the O2 delivery pathway and mechanism in cytochrome ba3, extensive molecular dynamics simulations of the enzyme in its membrane-embedded form have been performed, including different protocols of explicit ligand sampling (flooding) simulations with O2, implicit ligand sampling analysis, and in silico mutagenesis. The results show that O2 diffuses to the active site exclusively via a Y-shaped hydrophobic tunnel with two 25-Å long membrane-accessible branches that coincide with the pathway previously suggested by the crystallographically identified xenon binding sites. The two entrances of the bifurcated tunnel of cytochrome ba3 are located within the lipid bilayer, where O2 is preferentially partitioned from the aqueous phase. The largest barrier to O2 migration within the tunnel is estimated to be only 1.5 kcal/mol, allowing O2 to reach the enzyme active site virtually impeded by one-dimensional diffusion once it reaches a tunnel entrance at the protein surface. Unlike other O2-utilizing proteins, the tunnel is "open" with no transient barriers observed due to protein dynamics. This unique low-barrier passage through the protein ensures that O2 transit through the protein is never rate-limiting.
Collapse
Affiliation(s)
- Paween Mahinthichaichan
- Department of Biochemistry, and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign , Urbana, Illinois 61801, United States
| | - Robert B Gennis
- Department of Biochemistry, and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign , Urbana, Illinois 61801, United States
| | - Emad Tajkhorshid
- Department of Biochemistry, and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign , Urbana, Illinois 61801, United States
| |
Collapse
|
32
|
Computation of Rate Constants for Diffusion of Small Ligands to and from Buried Protein Active Sites. Methods Enzymol 2016; 578:299-326. [DOI: 10.1016/bs.mie.2016.05.039] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
|
33
|
Pesce A, Bustamante JP, Bidon-Chanal A, Boechi L, Estrin DA, Luque FJ, Sebilo A, Guertin M, Bolognesi M, Ascenzi P, Nardini M. The N-terminal pre-A region of Mycobacterium tuberculosis 2/2HbN promotes NO-dioxygenase activity. FEBS J 2015; 283:305-22. [PMID: 26499089 DOI: 10.1111/febs.13571] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2015] [Revised: 08/09/2015] [Accepted: 08/16/2015] [Indexed: 01/24/2023]
Abstract
UNLABELLED A unique defense mechanisms by which Mycobacterium tuberculosis protects itself from nitrosative stress is based on the O2 -dependent NO-dioxygenase (NOD) activity of truncated hemoglobin 2/2HbN (Mt2/2HbN). The NOD activity largely depends on the efficiency of ligand migration to the heme cavity through a two-tunnel (long and short) system; recently, it was also correlated with the presence at the Mt2/2HbN N-terminus of a short pre-A region, not conserved in most 2/2HbNs, whose deletion results in a drastic reduction of NO scavenging. In the present study, we report the crystal structure of Mt2/2HbN-ΔpreA, lacking the pre-A region, at a resolution of 1.53 Å. We show that removal of the pre-A region results in long range effects on the protein C-terminus, promoting the assembly of a stable dimer, both in the crystals and in solution. In the Mt2/2HbN-ΔpreA dimer, access of heme ligands to the short tunnel is hindered. Molecular dynamics simulations show that the long tunnel branch is the only accessible pathway for O2 -ligand migration to/from the heme, and that the gating residue Phe(62)E15 partly restricts the diameter of the tunnel. Accordingly, kinetic measurements indicate that the kon value for peroxynitrite isomerization by Mt2/2HbN-ΔpreA-Fe(III) is four-fold lower relative to the full-length protein, and that NO scavenging by Mt2/2HbN-ΔpreA-Fe(II)-O2 is reduced by 35-fold. Therefore, we speculate that Mt2/2HbN evolved to host the pre-A region as a mechanism for preventing dimerization, thus reinforcing the survival of the microorganism against the reactive nitrosative stress in macrophages. DATABASE Coordinates and structure factors have been deposited in the Protein Data Bank under accession number 5AB8.
Collapse
Affiliation(s)
| | - Juan P Bustamante
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, University of Buenos Aires, Argentina
| | - Axel Bidon-Chanal
- Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, University of Barcelona, Santa Coloma de Gramenet, Spain
| | - Leonardo Boechi
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, University of Buenos Aires, Argentina
| | - Darío A Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, University of Buenos Aires, Argentina
| | - Francisco Javier Luque
- Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, University of Barcelona, Santa Coloma de Gramenet, Spain
| | - Anne Sebilo
- Department of Biochemistry, Microbiology and Bioinformatics, Laval University, Quebec, Canada
| | - Michel Guertin
- Department of Biochemistry, Microbiology and Bioinformatics, Laval University, Quebec, Canada
| | - Martino Bolognesi
- Department of Biosciences, University of Milan, Italy.,CNR-IBF and CIMAINA, University of Milan, Italy
| | - Paolo Ascenzi
- Interdepartmental Laboratory of Electron Microscopy, Roma Tre University, Rome, Italy.,National Institute of Biostructures and Biosystems, Rome, Italy
| | - Marco Nardini
- Department of Biosciences, University of Milan, Italy
| |
Collapse
|
34
|
Blumberger J. Recent Advances in the Theory and Molecular Simulation of Biological Electron Transfer Reactions. Chem Rev 2015; 115:11191-238. [DOI: 10.1021/acs.chemrev.5b00298] [Citation(s) in RCA: 238] [Impact Index Per Article: 26.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Jochen Blumberger
- Department of Physics and
Astronomy, University College London, Gower Street, London WC1E 6BT, U.K
| |
Collapse
|
35
|
De Sancho D, Kubas A, Wang PH, Blumberger J, Best RB. Identification of Mutational Hot Spots for Substrate Diffusion: Application to Myoglobin. J Chem Theory Comput 2015; 11:1919-27. [PMID: 26574395 PMCID: PMC6132223 DOI: 10.1021/ct5011455] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The pathways by which small molecules (substrates or inhibitors) access active sites are a key aspect of the function of enzymes and other proteins. A key problem in designing or altering such proteins is to identify sites for mutation that will have the desired effect on the substrate transport properties. While specific access channels have been invoked in the past, molecular simulations suggest that multiple routes are possible, complicating the analysis. This complexity, however, can be captured by a Markov State Model (MSM) of the ligand diffusion process. We have developed a sensitivity analysis of the resulting rate matrix, which identifies the locations where mutations should have the largest effect on the diffusive on rate. We apply this method to myoglobin, which is the best characterized example both from experiment and simulation. We validate the approach by translating the sensitivity parameter obtained from this method into the CO binding rates in myoglobin upon mutation, resulting in a semi-quantitative correlation with experiments. The model is further validated against an explicit simulation for one of the experimental mutants.
Collapse
Affiliation(s)
- David De Sancho
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, United Kingdom
- CIC nanoGUNE , Tolosa Hiribidea 76, 20018 Donostia-San Sebastian, Spain
- IKERBASQUE, Basque Foundation for Science , Maria Diaz de Haro 3, 48013 Bilbao, Spain
| | - Adam Kubas
- Department of Physics and Astronomy, University College London , Gower Street, London WC1E 6BT, United Kingdom
| | - Po-Hung Wang
- Department of Physics and Astronomy, University College London , Gower Street, London WC1E 6BT, United Kingdom
- Theoretical Molecular Science Laboratory , 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London , Gower Street, London WC1E 6BT, United Kingdom
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health , Bethesda, Maryland 20892-0520, United States
| |
Collapse
|
36
|
Oteri F, Baaden M, Lojou E, Sacquin-Mora S. Multiscale Simulations Give Insight into the Hydrogen In and Out Pathways of [NiFe]-Hydrogenases from Aquifex aeolicus and Desulfovibrio fructosovorans. J Phys Chem B 2014; 118:13800-11. [DOI: 10.1021/jp5089965] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Francesco Oteri
- Laboratoire
de Biochimie Théorique, CNRS UPR9080, Institut de Biologie Physico-Chimique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Marc Baaden
- Laboratoire
de Biochimie Théorique, CNRS UPR9080, Institut de Biologie Physico-Chimique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Elisabeth Lojou
- Bioénergétique
et Ingénierie des Protéines, Institut de Microbiologie
de la Méditerranée, CNRS, Aix Marseille University, 31 Chemin Joseph Aiguier, 13402 Marseille Cedex, France
| | - Sophie Sacquin-Mora
- Laboratoire
de Biochimie Théorique, CNRS UPR9080, Institut de Biologie Physico-Chimique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
37
|
Finkelmann AR, Stiebritz MT, Reiher M. Activation Barriers of Oxygen Transformation at the Active Site of [FeFe] Hydrogenases. Inorg Chem 2014; 53:11890-902. [DOI: 10.1021/ic501049z] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Arndt R. Finkelmann
- Laboratorium
für Physikalische
Chemie, ETH Zürich, Valdimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Martin T. Stiebritz
- Laboratorium
für Physikalische
Chemie, ETH Zürich, Valdimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratorium
für Physikalische
Chemie, ETH Zürich, Valdimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
38
|
de Poulpiquet A, Ranava D, Monsalve K, Giudici-Orticoni MT, Lojou E. Biohydrogen for a New Generation of H2/O2Biofuel Cells: A Sustainable Energy Perspective. ChemElectroChem 2014. [DOI: 10.1002/celc.201402249] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
|
39
|
Damas JM, Baptista AM, Soares CM. The Pathway for O2 Diffusion inside CotA Laccase and Possible Implications on the Multicopper Oxidases Family. J Chem Theory Comput 2014; 10:3525-31. [DOI: 10.1021/ct500196e] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- João M. Damas
- Instituto de Tecnologia Química
Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157 Oeiras, Portugal
| | - António M. Baptista
- Instituto de Tecnologia Química
Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157 Oeiras, Portugal
| | - Cláudio M. Soares
- Instituto de Tecnologia Química
Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780-157 Oeiras, Portugal
| |
Collapse
|
40
|
Vedha SA, Solomon RV, Venuvanalingam P. Atomic partitioning of M-H2 bonds in [NiFe] hydrogenase--a test case of concurrent binding. Phys Chem Chem Phys 2014; 16:10698-707. [PMID: 24756140 DOI: 10.1039/c4cp00526k] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The possibility of simultaneous addition of η(2)-H2 to both the metals (Ni and Fe) in the active site of the as isolated state of the enzyme (Ni-SI) is examined here by an atom-by-atom electronic energy partitioning based on the QTAIM method. Results show that the 4LS state prefers H2 removal than addition. Destabilization of the atomic basins of the thiolate bridges and decrease of the electrophilicity of the Fe and Ni, resulting in poor back donation to the CO ligand, are the bottlenecks that hamper dihydrogen activation simultaneously. The study helps to understand why such states are seldom accessed in the activation of dihydrogen. Moreover, Ni has been found to be the natural choice for the dihydrogen binding.
Collapse
Affiliation(s)
- Swaminathan Angeline Vedha
- Theoretical & Computational Chemistry Laboratory, School of Chemistry, Bharathidasan University, Tiruchirappalli 24, India.
| | | | | |
Collapse
|
41
|
Kubas A, De Sancho D, Best RB, Blumberger J. Aerobic damage to [FeFe]-hydrogenases: activation barriers for the chemical attachment of O2. Angew Chem Int Ed Engl 2014; 53:4081-4. [PMID: 24615978 PMCID: PMC4143129 DOI: 10.1002/anie.201400534] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2014] [Indexed: 11/11/2022]
Abstract
[FeFe]-hydrogenases are the best natural hydrogen-producing enzymes but their biotechnological exploitation is hampered by their extreme oxygen sensitivity. The free energy profile for the chemical attachment of O2 to the enzyme active site was investigated by using a range-separated density functional re-parametrized to reproduce high-level ab initio data. An activation free-energy barrier of 13 kcal mol(-1) was obtained for chemical bond formation between the di-iron active site and O2, a value in good agreement with experimental inactivation rates. The oxygen binding can be viewed as an inner-sphere electron-transfer process that is strongly influenced by Coulombic interactions with the proximal cubane cluster and the protein environment. The implications of these results for future mutation studies with the aim of increasing the oxygen tolerance of this enzyme are discussed.
Collapse
Affiliation(s)
- Adam Kubas
- Department of Physics and Astronomy, University College LondonGower Street, London WC1E 6BT (UK)
| | - David De Sancho
- Department of Chemistry, Cambridge UniversityLensfield Road, Cambridge CB2 1EW (UK)
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of HealthBethesda, MD 20892-0520 (USA)
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College LondonGower Street, London WC1E 6BT (UK)
| |
Collapse
|
42
|
Kubas A, De Sancho D, Best RB, Blumberger J. Aerobic Damage to [FeFe]-Hydrogenases: Activation Barriers for the Chemical Attachment of O2. Angew Chem Int Ed Engl 2014. [DOI: 10.1002/ange.201400534] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
|
43
|
Topin J, Diharce J, Fiorucci S, Antonczak S, Golebiowski J. O2 Migration Rates in [NiFe] Hydrogenases. A Joint Approach Combining Free-Energy Calculations and Kinetic Modeling. J Phys Chem B 2014; 118:676-81. [DOI: 10.1021/jp4093964] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Jérémie Topin
- Institut de Chimie de Nice, UMR 7272, Université de Nice-Sophia Antipolis, CNRS, Parc Valrose, 06108 Nice Cedex 2, France
| | - Julien Diharce
- Institut de Chimie de Nice, UMR 7272, Université de Nice-Sophia Antipolis, CNRS, Parc Valrose, 06108 Nice Cedex 2, France
| | - Sébastien Fiorucci
- Institut de Chimie de Nice, UMR 7272, Université de Nice-Sophia Antipolis, CNRS, Parc Valrose, 06108 Nice Cedex 2, France
| | - Serge Antonczak
- Institut de Chimie de Nice, UMR 7272, Université de Nice-Sophia Antipolis, CNRS, Parc Valrose, 06108 Nice Cedex 2, France
| | - Jérôme Golebiowski
- Institut de Chimie de Nice, UMR 7272, Université de Nice-Sophia Antipolis, CNRS, Parc Valrose, 06108 Nice Cedex 2, France
| |
Collapse
|
44
|
Wang PH, Bruschi M, De Gioia L, Blumberger J. Uncovering a dynamically formed substrate access tunnel in carbon monoxide dehydrogenase/acetyl-CoA synthase. J Am Chem Soc 2013; 135:9493-502. [PMID: 23713976 DOI: 10.1021/ja403110s] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The transport of small ligands to active sites of proteins is the basis of vital processes in biology such as enzymatic catalysis and cell signaling, but also of more destructive ones including enzyme inhibition and oxidative damage. Here, we show how a diffusion-reaction model solved by means of molecular dynamics and density functional theory calculations provides novel insight into the transport of small ligands in proteins. In particular, we unravel the existence of an elusive, dynamically formed gas channel, which CO2 takes to diffuse from the solvent to the active site (C-cluster) of the bifunctional multisubunit enzyme complex carbon monoxide dehydrogenase/acetyl-CoA synthase (CODH/ACS). Two cavities forming this channel are temporarily created by protein fluctuations and are not apparent in the X-ray structures. The ligand transport is controlled by two residues at the end of this tunnel, His113 and His116, and occurs on the same time scale on which chemical binding to the active site takes place (0.1-1 ms), resulting in an overall binding rate on the second time scale. We find that upon reduction of CO2 to CO, the newly formed Fe-hydroxy ligand greatly strengthens the hydrogen-bond network, preventing CO from exiting the protein through the same way that CO2 takes to enter the protein. This is the basis for directional transport of CO from the production site (C-cluster of CODH subunit) to the utilization site (A-cluster of ACS subunit). In view of these results, a general picture emerges of how large proteins guide small ligands toward their active sites.
Collapse
Affiliation(s)
- Po-hung Wang
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | | | | | | |
Collapse
|
45
|
Abbruzzetti S, Spyrakis F, Bidon-Chanal A, Luque FJ, Viappiani C. Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations. Phys Chem Chem Phys 2013; 15:10686-701. [PMID: 23733145 DOI: 10.1039/c3cp51149a] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.
Collapse
Affiliation(s)
- Stefania Abbruzzetti
- Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Parma, viale delle Scienze 7A, 43124, Parma, Italy
| | | | | | | | | |
Collapse
|
46
|
Winkler M, Esselborn J, Happe T. Molecular basis of [FeFe]-hydrogenase function: an insight into the complex interplay between protein and catalytic cofactor. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2013; 1827:974-85. [PMID: 23507618 DOI: 10.1016/j.bbabio.2013.03.004] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2012] [Revised: 02/21/2013] [Accepted: 03/08/2013] [Indexed: 12/20/2022]
Abstract
The precise electrochemical features of metal cofactors that convey the functions of redox enzymes are essentially determined by the specific interaction pattern between cofactor and enclosing protein environment. However, while biophysical techniques allow a detailed understanding of the features characterizing the cofactor itself, knowledge about the contribution of the protein part is much harder to obtain. [FeFe]-hydrogenases are an interesting class of enzymes that catalyze both, H2 oxidation and the reduction of protons to molecular hydrogen with significant efficiency. The active site of these proteins consists of an unusual prosthetic group (H-cluster) with six iron and six sulfur atoms. While H-cluster architecture and catalytic states during the different steps of H2 turnover have been thoroughly investigated during the last 20 years, possible functional contributions from the polypeptide framework were only assumed according to the level of conservancy and X-ray structure analyses. Due to the recent development of simpler and more efficient expression systems the role of single amino acids can now be experimentally investigated. This article summarizes, compares and categorizes the results of recent investigations based on site directed and random mutagenesis according to their informative value about structure function relationships in [FeFe]-hydrogenases. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.
Collapse
Affiliation(s)
- Martin Winkler
- Ruhr-Universität Bochum, Fakultät für Biologie und Biotechnologie, Lehrstuhl für Biochemie der Pflanzen, AG Photobiotechnologie, Bochum, Germany
| | | | | |
Collapse
|
47
|
Structural features of [NiFeSe] and [NiFe] hydrogenases determining their different properties: a computational approach. J Biol Inorg Chem 2012; 17:543-55. [PMID: 22286956 DOI: 10.1007/s00775-012-0875-2] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2011] [Accepted: 01/12/2012] [Indexed: 10/14/2022]
Abstract
Hydrogenases are metalloenzymes that catalyze the reversible reaction H(2)<->2H(+) + 2e(-), being potentially useful in H(2) production or oxidation. [NiFeSe] hydrogenases are a particularly interesting subgroup of the [NiFe] class that exhibit tolerance to O(2) inhibition and produce more H(2) than standard [NiFe] hydrogenases. However, the molecular determinants responsible for these properties remain unknown. Hydrophobic pathways for H(2) diffusion have been identified in [NiFe] hydrogenases, as have proton transfer pathways, but they have never been studied in [NiFeSe] hydrogenases. Our aim was, for the first time, to characterize the H(2) and proton pathways in a [NiFeSe] hydrogenase and compare them with those in a standard [NiFe] hydrogenase. We performed molecular dynamics simulations of H(2) diffusion in the [NiFeSe] hydrogenase from Desulfomicrobium baculatum and extended previous simulations of the [NiFe] hydrogenase from Desulfovibrio gigas (Teixeira et al. in Biophys J 91:2035-2045, 2006). The comparison showed that H(2) density near the active site is much higher in [NiFeSe] hydrogenase, which appears to have an alternative route for the access of H(2) to the active site. We have also determined a possible proton transfer pathway in the [NiFeSe] hydrogenase from D. baculatum using continuum electrostatics and Monte Carlo simulation and compared it with the proton pathway we found in the [NiFe] hydrogenase from D. gigas (Teixeira et al. in Proteins 70:1010-1022, 2008). The residues constituting both proton transfer pathways are considerably different, although in the same region of the protein. These results support the hypothesis that some of the special properties of [NiFeSe] hydrogenases could be related to differences in the H(2) and proton pathways.
Collapse
|
48
|
Smith DMA, Xiong Y, Straatsma TP, Rosso KM, Squier TC. Force-Field Development and Molecular Dynamics of [NiFe] Hydrogenase. J Chem Theory Comput 2012; 8:2103-14. [PMID: 26593842 DOI: 10.1021/ct300185u] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Classical molecular force-field parameters describing the structure and motion of metal clusters in [NiFe] hydrogenase enzymes can be used to compare the dynamics and thermodynamics of [NiFe] under different oxidation, protonation, and ligation circumstances. Using density functional theory (DFT) calculations of small model clusters representative of the active site and the proximal, medial, and distal Fe/S metal centers and their attached protein side chains, we have calculated classical force-field parameters for [NiFe] in reduced and oxidized states, including internal coordinates, force constants, and atom-centered charges. Derived force constants revealed that cysteinate ligands bound to the metal ions are more flexible in the Ni-B active site, which has a bridging hydroxide ligand, than in the Ni-C active site, which has a bridging hydride. Ten nanosecond all-atom, explicit-solvent MD simulations of [NiFe] hydrogenase in oxidized and reduced catalytic states established the stability of the derived force-field parameters in terms of Cα and metal cluster fluctuations. Average active site structures from the protein MD simulations are consistent with [NiFe] structures from the Protein Data Bank, suggesting that the derived force-field parameters are transferrable to other hydrogenases beyond the structure used for testing. A comparison of experimental H2-production rates demonstrated a relationship between cysteinate side chain rotation and activity, justifying the use of a fully dynamic model of [NiFe] metal cluster motion.
Collapse
Affiliation(s)
- Dayle M A Smith
- Pacific Northwest National Laboratory, P.O. Box 999, MSIN J4-33, Richland, Washington 99352, United States
| | - Yijia Xiong
- Pacific Northwest National Laboratory, P.O. Box 999, MSIN J4-33, Richland, Washington 99352, United States
| | - T P Straatsma
- Pacific Northwest National Laboratory, P.O. Box 999, MSIN J4-33, Richland, Washington 99352, United States
| | - Kevin M Rosso
- Pacific Northwest National Laboratory, P.O. Box 999, MSIN J4-33, Richland, Washington 99352, United States
| | - Thomas C Squier
- Pacific Northwest National Laboratory, P.O. Box 999, MSIN J4-33, Richland, Washington 99352, United States
| |
Collapse
|
49
|
Abou Hamdan A, Dementin S, Liebgott PP, Gutierrez-Sanz O, Richaud P, De Lacey AL, Rousset M, Bertrand P, Cournac L, Léger C. Understanding and tuning the catalytic bias of hydrogenase. J Am Chem Soc 2012; 134:8368-71. [PMID: 22540997 DOI: 10.1021/ja301802r] [Citation(s) in RCA: 92] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
When enzymes are optimized for biotechnological purposes, the goal often is to increase stability or catalytic efficiency. However, many enzymes reversibly convert their substrate and product, and if one is interested in catalysis in only one direction, it may be necessary to prevent the reverse reaction. In other cases, reversibility may be advantageous because only an enzyme that can operate in both directions can turnover at a high rate even under conditions of low thermodynamic driving force. Therefore, understanding the basic mechanisms of reversibility in complex enzymes should help the rational engineering of these proteins. Here, we focus on NiFe hydrogenase, an enzyme that catalyzes H(2) oxidation and production, and we elucidate the mechanism that governs the catalytic bias (the ratio of maximal rates in the two directions). Unexpectedly, we found that this bias is not mainly determined by redox properties of the active site, but rather by steps which occur on sites of the proteins that are remote from the active site. We evidence a novel strategy for tuning the catalytic bias of an oxidoreductase, which consists in modulating the rate of a step that is limiting only in one direction of the reaction, without modifying the properties of the active site.
Collapse
Affiliation(s)
- Abbas Abou Hamdan
- Laboratoire de Bioénergétique et Ingénierie des Protéines, Institut de Microbiologie de la Méditerranée, CNRS and Aix-Marseille Université , 31 chemin Joseph Aiguier, 13402 Marseille Cedex 20, France
| | | | | | | | | | | | | | | | | | | |
Collapse
|
50
|
Mechanistic insight into the blocking of CO diffusion in [NiFe]-hydrogenase mutants through multiscale simulation. Proc Natl Acad Sci U S A 2012; 109:6399-404. [PMID: 22493222 DOI: 10.1073/pnas.1121176109] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
[NiFe]-hydrogenases are fascinating biological catalysts with potential application in biofuel cells. However, a severe problem in practical application is the strong sensitivity of hydrogenase to gaseous inhibitor molecules such as CO and O(2). Recently, a number of successful protein engineering studies have been reported that aimed at lowering the access of diatomic inhibitors to the active site pocket, but the molecular mechanism conferring increased resistance remained unclear. Here we use a multiscale simulation approach combining molecular dynamics with a master equation formalism to explain the steady drop in CO diffusion rate observed for the mutants V74M L122A, V74M L122M, and V74M of Desulfovibrio fructosovorans [NiFe]-hydrogenase. We find that diffusion in these variants is controlled by two gates, one between residues 74 and 476 and the other between residues 74 and 122. The existence of two control points in different locations explains why the reduction in the experimental diffusion rate does not simply correlate with the width of the main gas channel. We also find that in the more effective mutation (V74M) CO molecules are still able to reach the active site through transitions that are gated by the microsecond dihedral motions of the side chain of R476 and the thermal fluctuations of the width of the gas channel defined by M74 and L122. Reflecting on the molecular information gained from simulation, we discuss future mutation experiments that could further lower the diffusion rates of small ligands inhibiting [NiFe]-hydrogenase.
Collapse
|