A nonlinear propagation of modified ion-acoustic (mIA) shock waves in a relativistic degenerate plasma (containing inertial viscous positive and negative ion fluids, relativistic electron fluids, and negatively charged immobile heavy ions) has been investigated theoretically. The modified Burgers (mB) and further modified Burgers (FmB) equations have been derived by adopting reductive perturbation technique. The solutions of both mB and FmB equations have been numerically analyzed to characterize the basic features of mIA shock waves. The basic properties (speed, amplitude, width, etc.) of these electrostatic shock waves are found to be significantly modified by the effects of negatively charged static heavy ions and the plasma particle number densities. It is found that the properties of these shock waves obtained from this analysis are significantly different from those obtained from the analysis of standard Burgers equation. The implications of our results in space and interstellar compact objects like non-rotating white dwarfs, neutron stars, etc. are briefly discussed.