The internal structure constants k_j_ and the radius of gyration are useful tools for investigating the apsidal motion and tidal evolution of close binaries and planetary systems. These parameters are available for various evolutionary phases but they are scarce for the late stages of stellar evolution. To cover this gap, we present here the calculations of the apsidal-motion constants, the fractional radius of gyration, and the gravitational potential energy for two grids of cooling evolutionary sequences of white dwarfs and for neutron star models. The cooling sequences of white dwarfs were computed with LPCODE. An additional alternative to the white dwarf models was also adopted with the MESA code which allows non-stop calculations from the pre main-sequence (PMS) to the white dwarf cooling sequences. Neutron star models were acquired from the NSCool/TOV subroutines. The apsidal-motion constants, the moment of inertia and the gravitational potential energy were computed with a fourth-order Runge-Kutta method. The parameters are made available for four cooling sequences of white dwarfs (DA and DB types): 0.52, 0.57, 0.837 and 1.0M_{sun}_ and for neutron star models covering a mass range from 1.0 up to 2.183M_{sun}_, in 0.1 mass step. We show that, contrary to previously established opinion, the product of the form-factors {beta} and {alpha}, which are related to the moment of inertia, and gravitational potential energy, is not constant during some evolutionary phases. Regardless of the final products of stellar evolution (white dwarfs, neutron stars and perhaps black holes), we found that they recover the initial value of product {alpha}{beta} at the pre main-sequence phase (~0.4). These results may have important consequences for the investigation of the Jacobi virial equation.