We study a general relativistic gravitomagnetic 3-body effect induced by the spin angular momentum S X of a rotating mass M X orbited at distance r X by a local gravitationally bound restricted two-body system S of size r ≪ r X consisting of a test particle revolving around a massive body M. At the lowest post-Newtonian order, we analytically work out the doubly averaged rates of change of the Keplerian orbital elements of the test particle by finding non-vanishing long-term effects for the inclination I, the node Ω and the pericenter ω. Such theoretical results are confirmed by a numerical integration of the equations of motion for a fictitious 3-body system. We numerically calculate the magnitudes of the post-Newtonian gravitomagnetic 3-body precessions for some astronomical scenarios in our solar system. For putative man-made orbiters of the natural moons Enceladus and Europa in the external fields of Saturn and Jupiter, the relativistic precessions due to the angular momenta of the gaseous giant planets can be as large as ≃ 10 − 50 milliarcseconds per year mas yr −1 . A preliminary numerical simulation shows that, for certain orbital configurations of a hypothetical Europa orbiter, its range-rate signal ∆ρ can become larger than the current Doppler accuracy of the existing spacecraft Juno at Jupiter, i.e. σρ = 0.015 mm s −1 , after 1 d. The effects induced by the Sun's angular momentum on artificial probes of Mercury and the Earth are at the level of ≃ 1 − 0.1 microarcseconds per year µas yr −1 .