In the frame of an inertial observer at some position, which we can take to be $\vec r=0$, the Newtonian gravitational field at position $\vec r$ is
$$\vec g\equiv \ddot{\vec r} = -\frac{4\pi G}{3}\rho \,\vec r,\tag{1}$$
where $\rho$ is the density of the universe. This follows straightforwardly from the shell theorem. Of course, it's also just the second Friedmann equation for the case of a universe containing only nonrelativistic matter.
That seems too easy! What about the homogeneity of the universe, which equation (1) does not seem to respect? Well, there are two important points to be made.
Different parts of the universe are accelerating with respect to each other.
This is due to their mutual gravitational attraction and is precisely in accordance with equation (1). A static universe is inconsistent with Newtonian gravity.
Indeed, while it's often suggested that cosmic expansion is a general relativistic effect, it's also predicted by Newtonian gravity.
The gravitational field is frame dependent.
Specifically, it depends on the reference frame's acceleration. Since different parts of the universe are accelerating with respect to each other, their inertial frames (that is, frames only accelerated by gravity) measure different gravitational fields. In particular, the local gravitational field is always zero for an inertial observer. As a function of position, it follows equation (1) for every observer (if they define their own position to be $\vec r = 0$).
This is how equation (1) respects the symmetry of a uniform mass distribution.