The gravity fields of celestial bodies that possess an atmosphere are periodically perturbed by the redistribution of fluid mass associated with atmospheric dynamics. A component of this perturbation is due to the gravitational response of the body to the deformation of its surface induced by the atmospheric pressure loading. The magnitude of this effect depends on the relation between the loading and the response in terms of geopotential variations measured by the load Love numbers. In this work, we simulate and analyze the gravity field generated by the atmospheres of Venus and Mars by accounting for different models of their internal structure. By precisely characterizing the phenomena that drive the mass transportation in the atmosphere through general circulation models, we determine the effect of the interior structure on the response to the atmospheric loading. An accurate estimation of the time-varying gravity field, which measures the atmospheric contribution, may provide significant constraints on the interior structure through the measurement of the load Love numbers. A combined determination of tidal and load Love numbers would enhance our knowledge of the interior of planetary bodies, providing further geophysical constraints in the inversion of internal structure models.