With a computational model combining large eddy simulation and a discrete element model, detailed quantification of the bed load effects on turbulent open‐channel flows is presented. The objective is the revelation of bed load particle impact on the mean flow properties and coherent structures. Two comparative numerical experiments with mobile and immobile beds are conducted. Mean properties (e.g., velocity and Reynolds stress profiles) show good agreement with experimental data. Comparing the mobile and immobile cases, the effective bed position is nearly the same, whereas the equivalent sand roughness is changed. The flow experiences higher bottom shear stress over immobile bed. To quantify impact on turbulent structures, a revised quadrant analysis is performed to calculate four key parameters of ejection and sweep events (duration, maximum shear stress, transported momentum, and period). Results show that the ejection and sweep events have comparable importance in the outer region. However, sweep becomes dominant in the near‐wall region. The motion of particles enhances the sweep dominance by breaking up the ejection structures and decreasing their occurrence ratio. The results also suggest that the ejection events are easier to be influenced by the particle motions because they originate from the near‐wall region. The duration, maximum shear stress, and transported momentum decrease close to the bed. The period remains relatively constant in the outer region but decreases near the bed. Visualization of the coherent structure reveals that the instantaneous particle motion has strong correlation with the bursting cycle events.