To test a new interaction potential, molecular dynamics simulations of zeolite natrolite were performed for the structures under ambient conditions hydrated by perdeuterated water and at high pressure (1.87 GPa) in the superhydrated phase, which were recently studied by neutron diffraction. The experimental structures were reproduced with reasonable accuracy, and the hydrogen bond features are discussed. As in ordinary natrolite, a flip motion of water molecules around the HOH bisector is found, which, together with translational oscillations, gives rise to transient hydrogen bonds between water molecules, which do not appear from experimental equilibrium coordinates. The dynamics of water molecules can explain some problems encountered in refining the experimental structure. Vibrational spectra of natrolite containing perdeuterated water, which are not yet measured, were simulated, and their qualitative trend is discussed.