We perform calculations of structure functions for elastic and inelastic spin-dependent scattering of weakly interacting massive particles (WIMPs) off 125Te, 129Xe, and 131Xe. The nuclear structure calculations are performed in the microscopic interacting boson-fermion model (IBFM-2). In our calculations we employ one-body and leading long-range two-body WIMP-nucleus currents derived from chiral effective field theory. We demonstrate that the relevant matrix elements can be reliably computed in the IBFM-2, which will allow investigation of heavy deformed nuclei previously inaccessible to theoretical calculations.