Modelling and measuring effects of shear stress in extrusion bioprinting of endothelial- epithelial cell co-cultures

Extrusion-based bioprinting enables the development of tissue-like constructs; however, the impact of printing-associated shear stress on cell viability and function remains a critical consideration. To address this, we developed a comprehensive workflow combining rheological characterization, computational fluid dynamics (CFD) modelling, and experimental validation to predict and assess shear stress effects during bioprinting. The rheological properties of gelatin methacryloyl (GelMA) at 5 % (w/v, 20 °C) and 10 % (30 °C) concentrations were modelled, comparing various non- Newtonian regression models. CFD simulations were validated using micro-particle image velocimetry, showing agreement between predicted and measured velocities. The impact of bioprinting-associated shear stress on cell viability was assessed using a co-culture of human umbilical vein endothelial cells and breast epithelial cells. Immediate post-printing analysis revealed increased apoptosis in GelMA 5 % (w/v, 20 °C), although 10 % (w/v) GelMA demonstrated higher shear stress levels compared to 5 % GelMA. After 1 day of culture in crosslinked hydrogels, apoptosis increased in extrusion pressure, demonstrating the impact of low levels of acute shear stress. This workflow provides a robust methodology for predicting acute shear stress impacts during bioprinting, laying the foundation for future optimization studies.