The interactions with water of the diol epoxides (DEs) of both a planar and a nonplanar PAH have been examined using molecular dynamics. To determine probable water locations around the DE for later use in the study of DE protonation, molecular dynamics simulations using the OPLS force field were carried out on diol epoxides surrounded by a 22 A box of explicit water molecules. Results for 30 ps simulations indicate that 10-60% of the time, depending strongly on the conformation and type of the DE, there is a water molecule forming a hydrogen bond with the epoxide oxygen. The patterns seen in the frequency at which a DE binds a water molecule reflect patterns seen in the relationship between the type of PAH DE and amount of DNA adduct formation. Examination of the orientations and arrangements of the water and DEs during the simulations showed that the bound waters existed in several preferred configurations which are also dependent upon the PAH DE geometry.