Computer simulation of tumour resection-induced brain deformation by a meshless approach

Spatiotemporal single-cell tracking analysis in 3D tissues to reveal heterogeneous cellular response to mechanical stimuli

Mechanical stimuli have been recognized as important for tissue maturation, homeostasis and constructing engineered three-dimensional (3D) tissues. However, we know little about the cellular mechanical response in tissues that could be considerably heterogeneous and spatiotemporally dynamic due to the complex structure of tissues. Here, we report a spatiotemporal single-cell tracking analysis of in vitro 3D tissues under mechanical stretch, to reveal the heterogeneous cellular behavior by using a developed stretch and optical live imaging system. The system could affect the cellular orientation and directly measure the distance of cells in in vitro 3D myoblast tissues (3DMTs) at the single-cell level. Moreover, we observed the spatiotemporal heterogeneous cellular locomotion and shape changes under mechanical stretch in 3DMTs. This single-cell tracking analysis can become a principal method to investigate the heterogeneous cellular response in tissues and provide insights that conventional analyses have not yet offered.

SlicerCBM: automatic framework for biomechanical analysis of the brain

PURPOSE

Brain shift that occurs during neurosurgery disturbs the brain's anatomy. Prediction of the brain shift is essential for accurate localisation of the surgical target. Biomechanical models have been envisaged as a possible tool for such predictions. In this study, we created a framework to automate the workflow for predicting intra-operative brain deformations.

METHODS

We created our framework by uniquely combining our meshless total Lagrangian explicit dynamics (MTLED) algorithm for computing soft tissue deformations, open-source software libraries and built-in functions within 3D Slicer, an open-source software package widely used for medical research. Our framework generates the biomechanical brain model from the pre-operative MRI, computes brain deformation using MTLED and outputs results in the form of predicted warped intra-operative MRI.

RESULTS

Our framework is used to solve three different neurosurgical brain shift scenarios: craniotomy, tumour resection and electrode placement. We evaluated our framework using nine patients. The average time to construct a patient-specific brain biomechanical model was 3 min, and that to compute deformations ranged from 13 to 23 min. We performed a qualitative evaluation by comparing our predicted intra-operative MRI with the actual intra-operative MRI. For quantitative evaluation, we computed Hausdorff distances between predicted and actual intra-operative ventricle surfaces. For patients with craniotomy and tumour resection, approximately 95% of the nodes on the ventricle surfaces are within two times the original in-plane resolution of the actual surface determined from the intra-operative MRI.

CONCLUSION

Our framework provides a broader application of existing solution methods not only in research but also in clinics. We successfully demonstrated the application of our framework by predicting intra-operative deformations in nine patients undergoing neurosurgical procedures.

A position- and time-dependent pressure profile to model viscoelastic mechanical behavior of the brain tissue due to tumor growth

This study proposed a computational framework to calculate the resultant position- and time-dependent pressure profile on the brain tissue due to tumor growth. A finite element (FE) patch of the brain tissue was constructed and an inverse dynamic FE-optimization algorithm was used to calculate its viscoelastic mechanical properties under compressive uniaxial loading. Two patient-specific post-tumor resection FE models were input to the FE-optimization algorithm to calculate the optimized 3rd-order position-dependent and normal distribution time-dependent pressure profile parameters. The optimized viscoelastic material properties, the most suitable simulation time, and the optimized 3rd-order position- and -time-dependent pressure profiles were calculated.