The stress distribution along the trajectories of passive particles released in turbulent flow were computed with the use of Lagrangian methods and direct numerical simulations. The flow fields selected were transitional Poiseuille-Couette flow situations found in ventricular assist devices and turbulent flows at conditions found in blood pumps. The passive particle properties were selected to represent molecules of the von Willebrand factor (vWF) protein. Damage to the vWF molecule can cause disease, most often related to hemostasis. The hydrodynamic shear stresses along the trajectories of the particles were calculated and the changes in the distribution of stresses were determined for proteins released in different locations in the flow field and as a function of exposure time. The stress distributions indicated that even when the average applied stress was within a safe operating regime, the proteins spent part of their trajectories in flow areas of damaging stress. Further examination showed that the history of the distribution of stresses applied on the vWF molecules, rather than the average, should be used to evaluate hydrodynamically-induced damage.