Physical properties of the Sun (orientation of rotation axis, oblateness coefficient J2⊙, and change rate of the gravitational parameter μ˙⊙) are determined using a dynamical model describing the motion of the Sun, planets, the Moon, asteroids, and Trans-Neptunian objects (TNOs). Among the many kinds of observations used to determine the orbits and physical properties of the bodies, the most important for our study are precise interplanetary ranging data: Earth–Mercury ranges from MESSENGER spacecraft and Earth–Mars ranges from Odyssey and MRO. The findings allow us to improve the model of the Sun in modern planetary ephemerides. First, the dynamically determined direction of the Sun’s pole is ≈2° off the visible axis of rotation of the Sun’s surface, which is corroborated by present knowledge of the Sun’s interior. Second, the change rate of the Sun’s gravitational parameter is found to be smaller (in absolute value) than the nominal value derived from the estimate of mass loss through radiation and solar wind. Possible interpretations are discussed.