In this paper, we continue exploring the consequences of the general equation of motion (EOM) governing all Lagrangian perturbations ξ about a time-dependent, ideal magnetohydrodynamic (MHD) configuration, which includes self-gravity, external gravity, pressure gradients, compressibility, inertial effects, and anisotropic Lorentz force. We here address the specific case of MHD stability for 3D stationary equilibria, where the perturbed EOM features a symmetric operator F and an antisymmetric Doppler-Coriolis operator v·∇. For this case, we state and prove the general properties for the solutions ξ of the governing dynamical system. For axisymmetric perturbations about axisymmetric equilibria with purely toroidal, or purely poloidal magnetic fields, specific stability theorems can be formulated. We derive a useful integral expression for the quadratic quantity given by the inner product ⟨ξ,F[ξ]⟩. For deriving stability statements on MHD states where self-gravity is involved as well, we provide an upper bound on the perturbed self-gravitational energy associated with the displacement ξ. The resulting expression elucidates the role of potentially stabilizing versus destabilizing contributions and shows the role of gravity, entropy gradients, velocity shear, currents, Lorentz forces, inertia, and pressure gradients in offering many routes to unstable behavior in flowing gases and plasmas. These have historically mostly been studied for static v=0 configurations, looking at stability of exactly force-balanced states, or by assuming stationarity similar to our approach here (i.e., ∂t≡0 for the state we perturb), but typically in combination with some reduced dimensionality on the configuration of interest (translational or axisymmetry). We show that in these limits, we find and generalize expressions well-known from, e.g., the study of ideal MHD stability of tokamak plasmas or from Schwarzschild's criteria controlling convection in hydrodynamic, (external) gravitating systems. When applied to stationary, axisymmetric configurations in motion, we can use our upper bound to derive a sufficient stability criterion for perturbations of arbitrary azimuthal mode number m used in ξ(r)=η(r,z) exp(imφ).