Extended molecular dynamics simulation of the carbon monoxide migration in sperm whale myoglobin.

We report the results of an extended molecular dynamics simulation on the migration of photodissociated carbon monoxide in wild-type sperm whale myoglobin. Our results allow following one possible ligand migration dynamics from the distal pocket to the Xe1 cavity via a path involving the other xenon binding cavities and momentarily two additional packing defects along the pathway. Comparison with recent time resolved structural data obtained by Laue crystallography with subnanosecond to millisecond resolution shows a more than satisfactory agreement. In fact, according to time resolved crystallography, CO, after photolysis, can occupy the Xe1 and Xe4 cavities. However, no information on the trajectory of the ligand from the distal pocket to the Xe1 is available. Our results clearly show one possible path within the protein. In addition, although our data refer to a single trajectory, the local dynamics of the ligand in each cavity is sufficiently equilibrated to obtain local structural and thermodynamic information not accessible to crystallography. In particular, we show that the CO motion and the protein fluctuations are strictly correlated: free energy calculations of the migration between adjacent cavities show that the migration is not a simple diffusion but is kinetically or thermodynamically driven by the collective motions of the protein; conversely, the protein fluctuations are influenced by the ligand in such a way that the opening/closure of the passage between adjacent cavities is strictly correlated to the presence of CO in its proximity. The compatibility between time resolved crystallographic experiments and molecular dynamics simulations paves the way to a deeper understanding of the role of internal dynamics and packing defects in the control of ligand binding in heme proteins.

[1]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[2]  I. Schlichting,et al.  Crystal structure of photolysed carbonmonoxy-myoglobin , 1994, Nature.

[3]  M. Brunori,et al.  Cavities and packing defects in the structural dynamics of myoglobin , 2001, EMBO reports.

[4]  References , 1971 .

[5]  J. Straub,et al.  MOLECULAR DYNAMICS SIMULATION STUDY OF THE B-STATES OF SOLVATED CARBON MONOXYMYOGLOBIN , 1997 .

[6]  G. Kachalova,et al.  A steric mechanism for inhibition of CO binding to heme proteins. , 1999, Science.

[7]  P. Weiner,et al.  Computer Simulation of Biomolecular Systems , 1997 .

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

[9]  Andrea Amadei,et al.  Molecular dynamics simulations with constrained roto-translational motions: Theoretical basis and statistical mechanical consistency , 2000 .

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

[11]  I. Schlichting,et al.  Trapping intermediates in the crystal: ligand binding to myoglobin. , 2000, Current opinion in structural biology.

[12]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[13]  Aleksandr V. Smirnov,et al.  Watching a Protein as it Functions with 150-ps Time-Resolved X-ray Crystallography , 2003, Science.

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

[15]  D Bourgeois,et al.  Photolysis of the Carbon Monoxide Complex of Myoglobin: Nanosecond Time-Resolved Crystallography , 1996, Science.

[16]  R. Laskowski.  SURFNET: a program for visualizing molecular surfaces, cavities, and intermolecular interactions. , 1995, Journal of molecular graphics.

[17]  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 .

[18]  H Frauenfelder,et al.  Dynamics of ligand binding to myoglobin. , 1975, Biochemistry.

[19]  Alan E. Mark,et al.  The GROMOS96 Manual and User Guide , 1996 .

[20]  D Bourgeois,et al.  Protein conformational relaxation and ligand migration in myoglobin: a nanosecond to millisecond molecular movie from time-resolved Laue X-ray diffraction. , 2001, Biochemistry.

[21]  I. Kuntz,et al.  Protein-ligand dynamics. A 96 picosecond simulation of a myoglobin-xenon complex. , 1988, Journal of molecular biology.

[22]  M. Brunori,et al.  Nitric oxide moves myoglobin centre stage. , 2001, Trends in biochemical sciences.

[23]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[24]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[25]  M. Karplus,et al.  Substrate conformational transitions in the active site of chorismate mutase: Their role in the catalytic mechanism , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[26]  H. Frauenfelder,et al.  Conformational substates in proteins. , 1988, Annual review of biophysics and biophysical chemistry.

[27]  Q H Gibson,et al.  Ligand migration in sperm whale myoglobin. , 1997, Biochemistry.

[28]  F M Richards,et al.  Areas, volumes, packing and protein structure. , 1977, Annual review of biophysics and bioengineering.

[29]  M. Karplus,et al.  Molecular dynamics study of the photodissociation of carbon monoxide from myoglobin : ligand dynamics in the first 10 ps , 1991 .

[30]  R. Elber,et al.  Computer simulations of carbon monoxide photodissociation in myoglobin: structural interpretation of the B states. , 1998, Biophysical journal.

[31]  K. Moffat,et al.  Photolysis-induced structural changes in single crystals of carbonmonoxy myoglobin at 40 K , 1994, Nature Structural Biology.

[32]  G. Nienhaus,et al.  X-ray structure determination of a metastable state of carbonmonoxy myoglobin after photodissociation. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[33]  A. Miele,et al.  Complex landscape of protein structural dynamics unveiled by nanosecond Laue crystallography , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[34]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[35]  W. Eaton,et al.  Nonexponential structural relaxations in proteins , 1996 .

[36]  Berk Hess,et al.  Improving Efficiency of Large Time-Scale Molecular Dynamics Simulations of Hydrogen-Rich Systems , 1999 .

[37]  I. Kuntz,et al.  Cavities in proteins: structure of a metmyoglobin-xenon complex solved to 1.9 A. , 1984, Biochemistry.

[38]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[39]  P. Wolynes,et al.  The energy landscapes and motions of proteins. , 1991, Science.

[40]  M. Lim,et al.  Mid-infrared vibrational spectrum of CO after photodissociation from heme: Evidence for a ligand docking site in the heme pocket of hemoglobin and myoglobin , 1995 .