The contribution of essentially quantum internal molecular motions to the second virial coefficient B2 of water vapor is analyzed in the framework of the path integral approach. A general purpose ab initio polarizable force field QMPFF2 or a nonpolarizable three-site water model are used with oscillator and Morse valence potentials. It is demonstrated that the contribution may be significant but depends strongly on the form of the intramolecular potential. In the case of the more realistic stretching Morse potential, inclusion of quantum molecular flexibility into the simulation reduces the virial coefficient by 20%-40%. Also, the internal modes make a contribution to the difference in the virial coefficient for light and heavy water, which is opposite to that of the intermolecular motions, so that the net effect can even change the sign at higher temperatures.