The present work is dedicated to the first theoretical study of the rotationally inelastic collisions of Ne with H2O and its isotopologue D2O in an attempt to analyze the effect on the dynamics of H substitution by deuterium. To this aim two new potential energy surfaces are developed. Their quality is tested by computing the bound states of the complexes and comparing them with those most recently reported by other teams. System-specific collisional propensity rules are inferred for these two systems by analyzing the computed state-to-state cross sections at low and higher collision energy. The application of the Alexander parity index propensity rule is also discussed, and the present results are compared with those obtained for the collisions with other noble gases.