The plasmas in space, cosmic and astrophysical environments are long known to consist of numerous massive ionic components contributing to various wave instability fluctuation phenomena. Indeed, the ion-inertial effects need to be incorporated into realistic analyses, rather than treating the gravitating ionic species traditionally as a Boltzmann distributed fluid. Herein, we present an atypical theoretical model setup to study gravito-electrostatic mode-fluctuations in self-gravitating inhomogeneous interstellar dust molecular clouds (DMCs) on the astrophysical fluid scales of space and time. The main goal is focused on investigating the influence of self-consistent dynamic ion-inertial effects on the stability. Methodological application of standard multiple scaling techniques reduces the basic plasma structure equations into a unique pair of decoupled Korteweg–de Vries (KdV) equations for the weak fluctuations. In contrast, the fully nonlinear counterparts are shown to evolve as a new gravito-electrostatically coupled pair of the Sagdeev energy-integral equations. In both the perturbation regimes, excitation of two distinct eigenmode classes – electrostatic compressive solitons and self-gravitational rarefactive solitons with unusual and unique parametric features – is demonstrated and portrayed. The graphical shape analysis reflects new plasma conditions for such eigenspectral patterns to coevolve in realistic interstellar parameter windows hitherto remaining unexplored. It is seen that the inertial ions play a destabilizing influential role leading to enhanced fluctuations toward establishing a reorganized gravito-electrostatic equilibrium structure. Finally, we discuss the consistency of our results in the framework of existing inertialess ion theories, experimental findings and multiple space satellite-based observations, together with new implications.