Dynamics of water molecules in the active-site cavity of human cytochromes P450

We have studied the dynamics of water molecules in six crystal structures of four human cytochromes P450, 2A6, 2C8, 2C9, and 3A4, with molecular dynamics simulations. In the crystal structures, only a few water molecules are seen and the reported sizes of the active-site cavity vary a lot. In the simulations, the cavities are completely filled with water molecules, although with approximately 20% lower density than in bulk water. The 2A6 protein differs from the other three in that it has a very small cavity with only two water molecules and no exchange with the surroundings. The other three proteins have quite big cavities, with 41 water molecules on average in 2C8 and 54-58 in 2C9 and 3A4, giving a water volume of 1500-2100 A3. The two crystal structures of 2C9 differ quite appreciably, whereas those of 3A4 are quite similar. The active-site cavity is connected to the surroundings by three to six channels, through which there is a quite frequent exchange of water molecules (one molecule is exchanged every 30-200 ps), except in 2A6. Most of the channels are observed also in the crystal structures, but two to three channels in each protein open only during the simulations. There are no water molecules close to the heme iron ion in these simulations of the high-spin ferric state (the average distance to the closest water molecule is 3.3-5 A), and there are few ordered water molecules in the active sites, none of which is conserved in all proteins.