Oxygen diffusion in minihemoglobin from Cerebratulus lacteus: a locally enhanced sampling study

Heme proteins serve as a source of oxygen in nervous tissue during anoxia. The functional routes of a dioxygen (O2) diffusion in a novel structure of a minihemoglobin (CerHb) molecule present in worm Cerebratulus lacteus are not known. In this paper, the results of 1ns molecular dynamics simulations of this process are presented. The locally enhanced sampling method (LES) and CHARMM force field were used for simulations of CerHb with 1–15 copies of O2. It was found, that several alternative routes are possible. The dominant path consists of two steps. Firstly, ligands move from the heme pocket to a different cavity through the barrier defined by the residues Phe10 and Tyr48. Secondly, ligands leave the protein passing through the more complex barrier situated between the E/F loop and the H helix. We note that the number of paths observed depends on a number of LES copies of (O2).

[1]  T. Burmester,et al.  Cytoglobin: a novel globin type ubiquitously expressed in vertebrate tissues. , 2002, Molecular biology and evolution.

[2]  Gerhard Hummer,et al.  Unveiling functional protein motions with picosecond x-ray crystallography and molecular dynamics simulations. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[3]  M. Brunori,et al.  Molecular dynamics simulation of sperm whale myoglobin: effects of mutations and trapped CO on the structure and dynamics of cavities. , 2005, Biophysical journal.

[4]  Andrea Amadei,et al.  Extended molecular dynamics simulation of the carbon monoxide migration in sperm whale myoglobin. , 2004, Biophysical journal.

[5]  R. Elber,et al.  Modeling side chains in peptides and proteins: Application of the locally enhanced sampling and the simulated annealing methods to find minimum energy conformations , 1991 .

[6]  M. Teeter,et al.  Myoglobin cavities provide interior ligand pathway , 2004, Protein science : a publication of the Protein Society.

[7]  M. Karplus,et al.  Enhanced sampling in molecular dynamics: use of the time-dependent Hartree approximation for a simulation of carbon monoxide diffusion through myoglobin , 1990 .

[8]  M. Nardini,et al.  Thr-E11 Regulates O2 Affinity in Cerebratulus lacteus Mini-hemoglobin* , 2004, Journal of Biological Chemistry.

[9]  Wieslaw Nowak,et al.  Locally enhanced sampling in free energy calculations: Application of mean field approximation to accurate calculation of free energy differences , 1992 .

[10]  M. Karplus,et al.  A comparison between molecular dynamics and X-ray results for dissociated CO in myoglobin , 1997, Nature Structural Biology.

[11]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[12]  R. G. Hart,et al.  Structure of Myoglobin: A Three-Dimensional Fourier Synthesis at 2 Å. Resolution , 1960, Nature.

[13]  M. Perutz,et al.  Regulation of oxygen affinity of hemoglobin: influence of structure of the globin on the heme iron. , 1979, Annual review of biochemistry.

[14]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[15]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[16]  Thomas Hankeln,et al.  A vertebrate globin expressed in the brain , 2000, Nature.

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

[18]  M. Meuwly,et al.  CO migration in native and mutant myoglobin: atomistic simulations for the understanding of protein function. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[19]  M. Nardini,et al.  The 109 residue nerve tissue minihemoglobin from Cerebratulus lacteus highlights striking structural plasticity of the alpha-helical globin fold. , 2002, Structure.