A Computational Template for Three-Dimensional Modeling of the Vascular Scaffold of the Human Thyroid Gland

We recently designed an innovative scaffold-bioreactor unit for the bioengineering of a three-dimensional (3D) bioartificial human thyroid gland or its miniaturized replica as a part of a microfluidic chip test system. This device is based on the evidence that the 3D geometry of the intraglandular stromal/vascular scaffold (SVS; i.e., the fibrous and vascular matrix) of mammalian viscera plays a key role in guiding growth and differentiation of in vitro seeded cells. Therefore, we initiated a research program focused on computer-aided reconstruction of the 2nd to 4th order intralobar arterial network (IAN) of the human thyroid gland as a reliable surrogate for its 3D SVS, to be used as an input for rapid prototyping of a biomaterial replica. To this end, we developed a computational template that works within the Mathematica environment, giving rise to a quasi-fractal growth of the IAN distribution, constrained within an approximation of the thyroid lobe shape as a closed surface. Starting from edge detection of planar images of real human thyroid lobes acquired by in vivo real-time ultrasonography, we performed data approximation of the lobar profiles based on splines and Bezier curves, providing 3D lobar shapes as geometric boundaries for vessel growth by a diffusion-limited aggregation model. Our numerical procedures allowed for a robust connection between development of lobar arterial trees and thyroid lobe shape, led to a vascular self-similarity consistent with that of a cadaveric lobar arterial cast, and reproduced arterial vessels in a proportion not statistically different from that described for the real human thyroid gland. We conclude that our algorithmic template offers a reliable reproduction of the extremely complex IAN of the adult human thyroid lobe, potentially useful as a computational guidance for bioprinting of thyroid lobe matrix replicas. In addition, due to the simplicity and limited number of morphometrical parameters required by our system, we predict its application to the design of a number of patient-tailored human bioartificial organs and organs-on-chip, including parenchymal viscera and bones. Impact statement The study introduces the computer simulation of the three-dimensional (3D) intrinsic vascular matrix of the human thyroid gland, offering a general concept applicable to a number of other human viscera. Indeed, it provides a flexible software tool for reproduction of a 3D surrogate of the organ's 3D stromal matrix, suitable for eventual 3D bioprinting with biomaterials, and recellularization with organ-specific stem cells/progenitors. The final expectation is the design of patient-tailored 3D organ's matrices upon clinical request.