This Article Full Text Full Text (PDF) Supplement All Versions of this Article: biophysj.106.085746v1 91/5/1844    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Cohen, J. Articles by Schulten, K. Search for Related Content PubMed PubMed Citation Articles by Cohen, J. Articles by Schulten, K.

Imaging the Migration Pathways for O₂, CO, NO, and Xe Inside Myoglobin

Department of Physics and Beckman Institute, University of Illinois, Urbana, Illinois

Correspondence: Address reprint requests to Klaus Schulten, Tel.: 217-244-1604; Fax: 217-244-6078; E-mail: kschulte{at}ks.uiuc.edu.

Myoglobin (Mb) is perhaps the most studied protein, experimentally and theoretically. Despite the wealth of known details regarding the gas migration processes inside Mb, there exists no fully conclusive picture of these pathways. We address this deficiency by presenting a complete map of all the gas migration pathways inside Mb for small gas ligands (O₂, NO, CO, and Xe). To accomplish this, we introduce a computational approach for studying gas migration, which we call implicit ligand sampling. Rather than simulating actual gas migration events, we infer the location of gas migration pathways based on a free-energy perturbation approach applied to simulations of Mb's dynamical fluctuations at equilibrium in the absence of ligand. The method provides complete three-dimensional maps of the potential of mean force of gas ligand placement anywhere inside a protein-solvent system. From such free-energy maps we identify each gas docking site, the pathways between these sites, to the heme and to the external solution. Our maps match previously known features of these pathways in Mb, but also point to the existence of additional exits from the protein matrix in regions that are not easily probed by experiment. We also compare the pathway maps of Mb for different gas ligands and for different animal species.

This article has been cited by other articles:

				J. Z. Ruscio, D. Kumar, M. Shukla, M. G. Prisant, T. M. Murali, and A. V. Onufriev Atomic level computational identification of ligand migration pathways between solvent and binding site in myoglobin PNAS, July 8, 2008; 105(27): 9204 - 9209. [Abstract] [Full Text] [PDF]

				M. Anselmi, A. Di Nola, and A. Amadei The Kinetics of Ligand Migration in Crystallized Myoglobin as Revealed by Molecular Dynamics Simulations Biophys. J., June 1, 2008; 94(11): 4277 - 4281. [Abstract] [Full Text] [PDF]

				K. Kubiak and W. Nowak Molecular Dynamics Simulations of the Photoactive Protein Nitrile Hydratase Biophys. J., May 15, 2008; 94(10): 3824 - 3838. [Abstract] [Full Text] [PDF]

				R. E. Alcantara, C. Xu, T. G. Spiro, and V. Guallar A quantum-chemical picture of hemoglobin affinity PNAS, November 20, 2007; 104(47): 18451 - 18455. [Abstract] [Full Text] [PDF]

				J. Cohen and K. Schulten O2 Migration Pathways Are Not Conserved across Proteins of a Similar Fold Biophys. J., November 15, 2007; 93(10): 3591 - 3600. [Abstract] [Full Text] [PDF]