It's well-known/documented that you can construct the vertices of a regular dodecahedron from the following Cartesian coordinates (with appropriate rotation/scaling):
$(±1, ±1, ±1)\\ (0, ±\varphi, ±\frac{1}{\varphi})\\ (±\frac{1}{\varphi}, 0, ±\varphi)\\ (±\varphi, ±\frac{1}{\varphi}, 0)$
It's also well-known that you can construct the regular icosahedron in spherical coordinates with one point at each of the latitudes $±\pi/2$, and 5 points at each of the latitudes $±\arctan(1/2)$, rotationally symmetric except with the top and bottom hemispheres offset by $\pi/5$ radians.
Since the icosahedron and dodecahedron are duals of each other, this approach clearly should extend to the dodecahedron. For two different latitudes $±x$ and $±y$, each will have 5 vertices evenly spaced around them, for a total of the 20 vertices. The upper hemisphere will be offset in longitude $\pi/5$ from the lower hemisphere. The question is, what are $x$ and $y$?
The orientation of the dodecahedron in Cartesian coordinates you cite is of course not the same as the orientation that would be produced by taking the dual of the icosahedron in spherical coordinates.
The desired latitude angles are, $$\pm \tan^{-1} \frac{3-\sqrt{5}}{4}, \quad \pm \tan^{-1} \frac{3+\sqrt{5}}{4}.$$ Equivalently, these are $$\pm \cot^{-1} (3 + \sqrt{5}), \quad \pm \cot^{-1} (3 - \sqrt{5}).$$ The computation is rather tedious to derive from the spherical coordinates of the icosahedron, but the way to do it is to observe that the centroid of an icosahedral face may be computed as the average of its three vertices in Cartesian coordinates, and after converting this to spherical coordinates, ignore the radius and longitudinal angle.
So instead of doing such a calculation, it may be better to use the Cartesian coordinates for the dodecahedron that you cited. Specifically, we want to calculate the angle at the origin between the centroid of a dodecahedral face to any of its vertices. If we call this angle $\alpha$, then $\pi/2 - \alpha$ will be the more extreme latitude of the two pairs of five vertices. Then if we calculate the angle at the origin between any two adjacent vertices, and call this $\beta$, then $\pi/2 - (\alpha + \beta)$ will be the less extreme latitude of the two pairs. Since for your given coordinates the dodecahedral circumradius is $R = \sqrt{3}$, in my opinion these calculations are much more tractable than the first approach.