Simulating flow induced migration in vascular remodelling

Author summary Shear stress is known to play a leading role in endothelial cell (EC) migration and hence, vascular remodelling. Vascular remodelling is, however, more complicated than only EC migration. To achieve a better understanding of this process, we developed a computational model in which, shear stress mediated EC migration has the leading role and other factors, such as avascular regions and EC elongation, are also accounted for. We have tested this model for different vessel shapes during remodelling and could study the role that each of these factors play in remodelling. This model gives us the possibility of addition of other factors such as biochemical signals and angiogenesis which will help us in the study of vascular remodelling in both development and disease. Shear stress induces directed endothelial cell (EC) migration in blood vessels leading to vessel diameter increase and induction of vascular maturation. Other factors, such as EC elongation and interaction between ECs and non-vascular areas are also important. Computational models have previously been used to study collective cell migration. These models can be used to predict EC migration and its effect on vascular remodelling during embryogenesis. We combined live time-lapse imaging of the remodelling vasculature of the quail embryo yolk sac with flow quantification using a combination of micro-Particle Image Velocimetry and computational fluid dynamics. We then used the flow and remodelling data to inform a model of EC migration during remodelling. To obtain the relation between shear stress and velocityin vitrofor EC cells, we developed a flow chamber to assess how confluent sheets of ECs migrate in response to shear stress. Using these data as an input, we developed a multiphase, self-propelled particles (SPP) model where individual agents are driven to migrate based on the level of shear stress while maintaining appropriate spatial relationship to nearby agents. These agents elongate, interact with each other, and with avascular agents at each time-step of the model. We compared predicted vascular shape to real vascular shape after 4 hours from our time-lapse movies and performed sensitivity analysis on the various model parameters. Our model shows that shear stress has the largest effect on the remodelling process. Importantly, however, elongation played an especially important part in remodelling. This model provides a powerful tool to study the input of different biological processes on remodelling.