Abstract Putative natural massive satellites (exomoons) have gained increasing attention when they orbit Jupiter-like planets within the habitable zone of their host main-sequence star S. An exomoon s is expected to move within the equatorial plane of its host planet p, with its spin S s aligned with its orbital angular momentum L , which, in turn, is parallel to the planetary spin S p. If, in particular, the common tilt ε of such angular momenta to the plane of the satellite–planet motion about the star, assumed fixed, has certain values, the stellar latitudinal irradiation experienced on the exomoon may allow it to sustain life as we know it, at least for certain orbital configurations. An Earth analog (similar in mass, radius, oblateness, and obliquity) is considered, which orbits within 5–10 planetary radii R p from its Jupiter-like host planet. The de Sitter and Lense–Thirring spin precessions due to the general relativistic post-Newtonian (pN) field of the host planet have an impact on an exomoon’s habitability for a variety of different initial spin–orbit configurations. Here I show it by identifying long-term variations in the satellite’s obliquity ε s, where variations can be ≲10°–100°, depending on the initial spin–orbit configuration, with a timescale of ≃0.1–1 million years. Also the satellite’s quadrupole mass moment J 2 s induces obliquity variations that are faster than the pN ones but do not cancel them. Tidal dissipations, which may potentially have a relevant impact on the outlined pattern, are not included in the present analysis.