Intermolecular potential models for water and alkanes describe pure component properties fairly well, but fail to reproduce properties of water-alkane mixtures. Understanding interactions between water and non-polar molecules like alkanes is important not only for the hydrocarbon industry but has implications to biological processes as well. Although non-polar solutes in water have been widely studied, much less work has focused on water in non-polar solvents. In this study we calculate the solubility of water in different alkanes (methane to dodecane) at ambient conditions where the water content in alkanes is very low so that the non-polar water-alkane interactions determine solubility. Only the alkane-rich phase is simulated since the fugacity of water in the water rich phase is calculated from an accurate equation of state. Using the SPC/E model for water and TraPPE model for alkanes along with Lorentz-Berthelot mixing rules for the cross parameters produces a water solubility that is an order of magnitude lower than the experimental value. It is found that an effective water Lennard-Jones energy εW/k = 220 K is required to match the experimental water solubility in TraPPE alkanes. This number is much higher than used in most simulation water models (SPC/E—εW/k = 78.2 K). It is surprising that the interaction energy obtained here is also higher than the water-alkane interaction energy predicted by studies on solubility of alkanes in water. The reason for this high water-alkane interaction energy is not completely understood. Some factors that might contribute to the large interaction energy, such as polarizability of alkanes, octupole moment of methane, and clustering of water at low concentrations in alkanes, are examined. It is found that, though important, these factors do not completely explain the anomalously strong attraction between alkanes and water observed experimentally.