| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Department of Chemistry, University of Basel, CH-4056 Basel, Switzerland
Correspondence: Address reprint requests to Markus Meuwly, Klingelbergstrasse 80, CH-4056 Basel, Switzerland. Tel.: 41-61-267-3821; E-mail: m.meuwly{at}unibas.ch.
Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm-1 (compared to the experimental value of 10 cm-1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol-1 and agrees well with the experimentally estimated value of 0.29 kcal mol-1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.
This article has been cited by other articles:
 |
|
 |
|
  N. Plattner and M. Meuwly The Role of Higher CO-Multipole Moments in Understanding the Dynamics of Photodissociated Carbonmonoxide in Myoglobin Biophys. J., April 1, 2008; 94(7): 2505 - 2515. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. Anselmi, M. Aschi, A. Di Nola, and A. Amadei Theoretical Characterization of Carbon Monoxide Vibrational Spectrum in Sperm Whale Myoglobin Distal Pocket Biophys. J., May 15, 2007; 92(10): 3442 - 3447. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 E. Pinakoulaki, H. Yoshimura, V. Daskalakis, S. Yoshioka, S. Aono, and C. Varotsis Two ligand-binding sites in the O2-sensing signal transducer HemAT: Implications for ligand recognition/discrimination and signaling PNAS, October 3, 2006; 103(40): 14796 - 14801. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
D. R. Nutt and M. Meuwly Studying Reactive Processes with Classical Dynamics: Rebinding Dynamics in MbNO Biophys. J., February 15, 2006; 90(4): 1191 - 1201. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
G. Hummer, F. Schotte, and P. A. Anfinrud Unveiling functional protein motions with picosecond x-ray crystallography and molecular dynamics simulations PNAS, October 26, 2004; 101(43): 15330 - 15334. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
D. R. Nutt and M. Meuwly CO migration in native and mutant myoglobin: Atomistic simulations for the understanding of protein function PNAS, April 20, 2004; 101(16): 5998 - 6002. [Abstract] [Full Text] [PDF]
|
|
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
