Abstract
A precise 100-ns molecular dynamics simulation in aquo was performed for the heterotetrameric sarcosine oxidase bound with a substrate analogue, dimethylglycine. The spatial region including the protein was divided into small rectangular cells. The average number of the water molecules locating within each cell was calculated based on the simulation trajectory. The clusters of the cells filled with water molecules were used to determine the water channels. The narrowness of the channels, the average hydropathy indices of the residues of the channels, and the number of migration events of water molecules through the channels were consistent with the selective transport hypothesis whereby tunnel T3 is the pathway for the exit of the iminium intermediate of the enzyme reaction.