CO escape from myoglobin with metadynamics simulations

The relatively small size of myoglobin makes it suitable for the investigation of the ligand escape process in respiratory proteins and, in general, an ideal model system for the study of the more general structure-function paradigm. In this work, we use Molecular Dynamics simulations combined with an accelerated algorithm, the metadynamics, to probe the escape of CO from myoglobin. Our approach permits to quantitatively describe the escape process via the reconstruction of the associated free energy surface. Additionally, hints on the involvement of a larger numbers of residues than hitherto assumed in the gating process are extracted from our data.