Recent experimental observations at JET show evidence of reduced ion temperature profile stiffness. An extensive set of nonlinear gyrokinetic simulations are performed based on the experimental discharges, investigating the physical mechanism behind the observations. The impact on the ion heat flux of various parameters that differ within the data-set are explored. These parameters include the safety factor, magnetic shear, toroidal flow shear, effect of rotation on the magnetohydrodynamic equilibrium, R/L_n, β_e, Z_eff, T_e/T_i, and the fast-particle content. While previously hypothesized to be an important factor in the stiffness reduction, the combined effect of toroidal flow shear and low magnetic shear is not predicted by the simulations to lead to a significant reduction in ion heat flux, due both to an insufficient magnitude of flow shear and significant parallel velocity gradient destabilization. It is however found that nonlinear electromagnetic effects due to both thermal and fast-particle pressure gradients, even at low β_e, can significantly reduce the ion heat flux, and is a key factor in explaining the experimental observations. A total of four discharges are examined, at both inner and outer radii. For all cases studied, the simulated and experimental ion heat flux values agree within reasonable variations of input parameters around the experimental uncertainties.