Collisionless plasma shocks are a common feature of many space and astrophysical systems. They are sources of high-energy particles and nonthermal emission, channeling as much as 20% of the shock’s energy into nonthermal particles. The generation and acceleration of these nonthermal particles have been previously studied and shown to affect shock hydrodynamics to the zeroth order. In this work, we use self-consistent hybrid particle-in-cell simulations to examine the effect of self-generated nonthermal ion populations on the nature of collisionless, quasi-parallel shocks. Accelerated nonthermal particles downstream of the shock diffuse into the upstream region, taking energy away from the shock, which increases the compression ratio, slows the shock down, and flattens the nonthermal population’s spectral index for lower-Mach-number shocks. We show that this enhances shock compressibility when the heat flux is included in the Rankine–Hugoniot jump conditions, results that are roughly consistent with previous theories of “cosmic-ray-modified shocks.” Additionally, the simulation data show that heat flux and enthalpy flux cancels out in the upstream region, yielding a relatively simple, alternative closure for the jump conditions which accurately predict for the shock speed and compression ratio. The results have the potential to explain discrepancies between predictions and observations in a wide range of systems, such as inaccuracies in predictions of the arrival times of coronal mass ejections and the conflicting radio and X-ray observations of intracluster shocks. These effects will likely need to be included in fluid modeling to predict shock evolution accurately.