We establish nonlinear L1∩H3 → Lp orbital stability, 2 p ∞, with sharp rates of decay, of large-amplitude Lax-type shock profiles for a class of symmetric hyperbolic-parabolic systems including compressible gas dynamics and magnetohydrodynamics (MHD) under the necessary conditions of strong spectral stability, i.e., a stable point spectrum of the linearized operator about the wave, transversality of the profile, and hyperbolic stability of the associated ideal shock. This yields in particular, together with the spectral stability results of [50], the nonlinear stability of arbitrarily large-amplitude shock profiles of isentropic Navier-Stokes equations for a gamma-law gas as γ → 1: the first complete large-amplitude stability result for a shock profile of a system with real (i.e., partial) viscosity. A corresponding small-amplitude result was established in [53, 54] for general systems of Kawashima class by a combination of "Kawashima-type" energy estimates and pointwise Green function bounds, where the small-amplitude assumption was used only to close the energy estimates. Here, under the mild additional assumption that hyperbolic characteristic speeds (relative to the shock) are not only nonzero but of a common sign, we close the estimates instead by use of a Goodman-type weighted norm [25, 26] designed to control estimates in the crucial hyperbolic modes.