Gamma-ray binaries are binary systems that show non-thermal broadband emission from radio to gamma rays. If the system comprises a massive star and a young non-accreting pulsar, their winds collide producing non-thermal emission, most likely from the shocked pulsar wind. Thermal X-rays are expected from the shocked stellar wind, with a spectrum akin to the one observed in massive star binaries. The goal of this work is, through the study of the thermal X-ray emission from the shocked stellar wind in pulsar gammaray binaries, constrain the pulsar spin-down luminosity and the stellar wind properties. A semi-analytic model is developed to compute the thermal X-ray emission from the shocked stellar wind in pulsar gamma-ray binaries. The model results are compared with XMM-Newton observations of LS 5039, a candidate pulsar gamma-ray binary with a strong stellar wind. Exploring the range of possible values for the stellar mass-loss rate and orbital inclination, we obtain an upper limit on the pulsar spin-down luminosity of 6 × 10 36 erg s −1 . We conclude that, to explain the non-thermal luminosity of LS 5039 in the pulsar wind scenario, a non-thermal to spin-down luminosity ratio very close to unity may be required.