科研之友
微信
新浪微博
Facebook
分享链接
摘要
Myoglobin has important biological functions in storing and transporting small diatomic molecules in human body. Two possible orientations of carbon monoxide (CO) in the heme distal pocket (named as B-1 and B-2 states) of myoglobin have been experimentally indicated. In this study, ab initio quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulation of CO in myoglobin was carried out to investigate the two possible B states. Our results demonstrate that the B-1 and B-2 states correspond to Fe center dot center dot center dot CO (with carbon atom closer to iron center of heme) and Fe center dot center dot center dot OC (with oxygen atom closer to Fe), by comparing with the experimental infrared spectrum. QM electrostatic polarization effect on CO brought from the protein and solvent environment is the main driving force, which anchors CO in two distinctive orientations and hinders its rotation. The calculated vibrational frequency shift between the state B-1 and B-2 is 13.1 cm(-1), which is in good agreement with experimental value of 11.5 cm(-1). This study also shows that the electric field produced by the solvent plays an important role in assisting protein functions by exerting directional electric field at the active site of the protein. From residue-based electric field decomposition, several residues were found to have most contributions to the total electric field at the CO center, including a few charged residues and three adjacent uncharged polar residues (namely, HIS64, ILE107, and PHE43). This study provides new physical insights on rational design of enzyme with higher electric field at the active site.
