A quantitative model for oxygen uptake and release in a family of hemeproteins

MOTIVATION Hemeproteins have many diverse functions that largely depend on the rate at which they uptake or release small ligands, like oxygen. These proteins have been extensively studied using either simulations or experiments, albeit only qualitatively and one or two proteins at a time. RESULTS We present a physical-chemical model, which uses data obtained exclusively from computer simulations, to describe the uptake and release of oxygen in a family of hemeproteins, called truncated hemoglobins (trHbs). Through a rigorous statistical analysis we demonstrate that our model successfully recaptures all the reported experimental oxygen association and dissociation kinetic rate constants, thus allowing us to establish the key factors that determine the rates at which these hemeproteins uptake and release oxygen. We found that internal tunnels as well as the distal site water molecules control ligand uptake, whereas oxygen stabilization by distal site residues controls ligand release. Because these rates largely determine the functions of these hemeproteins, these approaches will also be important tools in characterizing the trHbs members with unknown functions. CONTACT lboechi@ic.fcen.uba.ar SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

[1]  Juan P. Bustamante,et al.  Ligand uptake modulation by internal water molecules and hydrophobic cavities in hemoglobins. , 2014, The journal of physical chemistry. B.

[2]  K. Schulten,et al.  Single-Molecule Experiments in Vitro and in Silico , 2007, Science.

[3]  Markus Meuwly,et al.  Nitric oxide dynamics in truncated hemoglobin: docking sites, migration pathways, and vibrational spectroscopy from molecular dynamics simulations. , 2009, Biophysical journal.

[4]  F. J. Luque,et al.  Role of Pre-A Motif in Nitric Oxide Scavenging by Truncated Hemoglobin, HbN, of Mycobacterium tuberculosis* , 2009, Journal of Biological Chemistry.

[5]  Kasper P Kepp,et al.  O2 binding to heme is strongly facilitated by near-degeneracy of electronic states. , 2013, Chemphyschem : a European journal of chemical physics and physical chemistry.

[6]  R. Gnanadesikan,et al.  Probability plotting methods for the analysis of data. , 1968, Biometrika.

[7]  F Javier Luque,et al.  Structural determinants of ligand migration in Mycobacterium tuberculosis truncated hemoglobin O , 2008, Proteins.

[8]  Victor Guallar,et al.  Ligand Migration in the Truncated Hemoglobin-II from Mycobacterium tuberculosis , 2009, Journal of Biological Chemistry.

[9]  J. Harvey,et al.  Spin-forbidden ligand binding to the ferrous-heme group: ab initio and DFT studies. , 2007, The journal of physical chemistry. B.

[10]  Suman Kundu,et al.  The Crystal Structure of Synechocystis Hemoglobin with a Covalent Heme Linkage* , 2004, Journal of Biological Chemistry.

[11]  Martino Bolognesi,et al.  The roles of Tyr(CD1) and Trp(G8) in Mycobacterium tuberculosis truncated hemoglobin O in ligand binding and on the heme distal site architecture. , 2007, Biochemistry.

[12]  P. Hunt,et al.  A hemoglobin from plants homologous to truncated hemoglobins of microorganisms , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[13]  Paul D. W. Kirk,et al.  Model Selection in Systems Biology Depends on Experimental Design , 2014, PLoS Comput. Biol..

[14]  Leonardo Boechi,et al.  Dioxygen affinity in heme proteins investigated by computer simulation. , 2006, Journal of inorganic biochemistry.

[15]  Rafael C. Bernardi,et al.  Molecular dynamics simulations of large macromolecular complexes. , 2015, Current opinion in structural biology.

[16]  Markus Meuwly,et al.  Atomistic simulation of NO dioxygenation in group I truncated hemoglobin. , 2010, Journal of the American Chemical Society.

[17]  Richard Daigle,et al.  Ligand Binding to Truncated Hemoglobin N from Mycobacterium tuberculosis Is Strongly Modulated by the Interplay between the Distal Heme Pocket Residues and Internal Water* , 2008, Journal of Biological Chemistry.

[18]  F. Yang,et al.  Crystal structures of CO-, deoxy- and met-myoglobins at various pH values. , 1996, Journal of molecular biology.

[19]  Michael H. Kutner Applied Linear Statistical Models , 1974 .

[20]  Pier Remigio Salvi,et al.  Following Ligand Migration Pathways from Picoseconds to Milliseconds in Type II Truncated Hemoglobin from Thermobifida fusca , 2012, PloS one.

[21]  M. Brunori,et al.  The role of cavities in protein dynamics: crystal structure of a photolytic intermediate of a mutant myoglobin. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[22]  Massimo Coletta,et al.  Ligand‐ and proton‐linked conformational changes of the ferrous 2/2 hemoglobin of Pseudoalteromonas haloplanktis TAC125 , 2011, IUBMB life.

[23]  David L Mobley,et al.  Alchemical free energy methods for drug discovery: progress and challenges. , 2011, Current opinion in structural biology.

[24]  David Dantsker,et al.  Reactions of Mycobacterium tuberculosis truncated hemoglobin O with ligands reveal a novel ligand-inclusive hydrogen bond network. , 2003, Biochemistry.

[25]  Amy E. Keating,et al.  Data-Driven Prediction and Design of bZIP Coiled-Coil Interactions , 2015, PLoS Comput. Biol..

[26]  J. S. Olson,et al.  Myoglobin discriminates between O2, NO, and CO by electrostatic interactions with the bound ligand , 1997, JBIC Journal of Biological Inorganic Chemistry.

[27]  Veronica Morea,et al.  A novel thermostable hemoglobin from the actinobacterium Thermobifida fusca , 2005, The FEBS journal.

[28]  Markus Meuwly,et al.  Oxygen migration pathways in NO-bound truncated hemoglobin. , 2012, Chemphyschem : a European journal of chemical physics and physical chemistry.

[29]  Nicola Marzari,et al.  Simulation of heme using DFT + U: a step toward accurate spin-state energetics. , 2007, The journal of physical chemistry. B.

[30]  Juan P. Bustamante,et al.  Ligand uptake in Mycobacterium tuberculosis truncated hemoglobins is controlled by both internal tunnels and active site water molecules , 2015, F1000Research.

[31]  S. Franzen,et al.  Spin-dependent mechanism for diatomic ligand binding to heme , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[32]  B. Wittenberg,et al.  Reaction of Mycobacterium tuberculosis Truncated Hemoglobin O with Hydrogen Peroxide , 2007, Journal of Biological Chemistry.

[33]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[34]  Paolo Ascenzi,et al.  A TyrCD1/TrpG8 hydrogen bond network and a TyrB10—TyrCD1 covalent link shape the heme distal site of Mycobacterium tuberculosis hemoglobin O , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[35]  M. Perutz,et al.  An x-ray study of azide methaemoglobin. , 1966, Journal of molecular biology.

[36]  Alessandra Pesce,et al.  Structural characterization of a group II 2/2 hemoglobin from the plant pathogen Agrobacterium tumefaciens. , 2011, Biochimica et biophysica acta.

[37]  Juan P. Bustamante,et al.  Structural flexibility of the heme cavity in the cold‐adapted truncated hemoglobin from the Antarctic marine bacterium Pseudoalteromonas haloplanktis TAC125 , 2015, The FEBS journal.

[38]  Martino Bolognesi,et al.  Truncated Hemoglobins: A New Family of Hemoglobins Widely Distributed in Bacteria, Unicellular Eukaryotes, and Plants* 210 , 2002, The Journal of Biological Chemistry.

[39]  Noam Agmon,et al.  Structure and energetics of the hydronium hydration shells. , 2007, The journal of physical chemistry. A.

[40]  Klaus Schulten,et al.  Finding gas migration pathways in proteins using implicit ligand sampling. , 2008, Methods in enzymology.

[41]  M Bolognesi,et al.  Mycobacterium tuberculosis hemoglobin N displays a protein tunnel suited for O2 diffusion to the heme , 2001, The EMBO journal.

[42]  António M. Baptista,et al.  Exploring O2 Diffusion in A-Type Cytochrome c Oxidases: Molecular Dynamics Simulations Uncover Two Alternative Channels towards the Binuclear Site , 2014, PLoS Comput. Biol..

[43]  Andrea Ilari,et al.  Crystal structure and ligand binding properties of the truncated hemoglobin from Geobacillus stearothermophilus. , 2007, Archives of biochemistry and biophysics.

[44]  Alessandra Pesce,et al.  Structural Determinants in the Group III Truncated Hemoglobin from Campylobacter jejuni* , 2006, Journal of Biological Chemistry.

[45]  Marianne Rooman,et al.  Stability Curve Prediction of Homologous Proteins Using Temperature-Dependent Statistical Potentials , 2014, PLoS Comput. Biol..

[46]  Veronica Morea,et al.  The Truncated Oxygen-avid Hemoglobin from Bacillus subtilis , 2005, Journal of Biological Chemistry.

[47]  F. J. Luque,et al.  Ligand‐induced dynamical regulation of NO conversion in Mycobacterium tuberculosis truncated hemoglobin‐N , 2006, Proteins.

[48]  V. Barnett,et al.  Applied Linear Statistical Models , 1975 .

[49]  Syun-Ru Yeh,et al.  Structural and Functional Properties of a Truncated Hemoglobin from a Food-borne Pathogen Campylobacter jejuni* , 2007, Journal of Biological Chemistry.

[50]  K. Schulten,et al.  Imaging the migration pathways for O2, CO, NO, and Xe inside myoglobin. , 2006, Biophysical journal.

[51]  F. J. Luque,et al.  Ligand migration in Methanosarcina acetivorans protoglobin: effects of ligand binding and dimeric assembly. , 2011, The journal of physical chemistry. B.

[52]  Martino Bolognesi,et al.  Ligand interactions in the distal heme pocket of Mycobacterium tuberculosis truncated hemoglobin N: roles of TyrB10 and GlnE11 residues. , 2006, Biochemistry.

[53]  M Brunori,et al.  Structural dynamics of myoglobin. , 2000, Biophysical chemistry.

[54]  R. A. Goldbeck,et al.  Water and ligand entry in myoglobin: Assessing the speed and extent of heme pocket hydration after CO photodissociation , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[55]  Leonardo Boechi,et al.  Small ligand-globin interactions: reviewing lessons derived from computer simulation. , 2013, Biochimica et biophysica acta.

[56]  D. Rousseau,et al.  Ligand binding in the ferric and ferrous states of Paramecium hemoglobin. , 2000, Biochemistry.

[57]  Lorna J. Smith,et al.  Understanding protein folding via free-energy surfaces from theory and experiment. , 2000, Trends in biochemical sciences.

[58]  Q. Gibson,et al.  Mapping the Pathways for O2 Entry Into and Exit from Myoglobin* , 2001, The Journal of Biological Chemistry.

[59]  Joel Friedman,et al.  Heme-Ligand Tunneling in Group I Truncated Hemoglobins*[boxs] , 2004, Journal of Biological Chemistry.

[60]  D. Rousseau,et al.  A cooperative oxygen-binding hemoglobin from Mycobacterium tuberculosis. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[61]  M Bolognesi,et al.  A novel two‐over‐two α‐helical sandwich fold is characteristic of the truncated hemoglobin family , 2000, The EMBO journal.

[62]  L. E. Laverman,et al.  Mechanistic studies of nitric oxide reactions with water soluble iron(II), cobalt(II), and iron(III) porphyrin complexes in aqueous solutions: implications for biological activity. , 2001, Journal of the American Chemical Society.

[63]  D. Rousseau,et al.  Structural investigations of the hemoglobin of the cyanobacterium Synechocystis PCC6803 reveal a unique distal heme pocket. , 2000, European journal of biochemistry.

[64]  Peter C. Ford,et al.  A Dissociative Mechanism for Reactions of Nitric Oxide with Water Soluble Iron(III) Porphyrins , 1997 .

[65]  R. Elber Ligand diffusion in globins: simulations versus experiment. , 2010, Current opinion in structural biology.

[66]  Dario A. Estrin,et al.  Evolutionary and Functional Relationships in the Truncated Hemoglobin Family , 2016, PLoS Comput. Biol..