The influence of biophysical parameters in a biomechanical model of cortical folding patterns

Abnormal cortical folding patterns, such as lissencephaly, pachygyria and polymicrogyria malformations, may be related to neurodevelopmental disorders. In this context, computational modeling is a powerful tool to provide a better understanding of the early brain folding process. Recent studies based on biomechanical modeling have shown that mechanical forces play a crucial role in the formation of cortical convolutions. However, the effect of biophysical parameters in these models remain unclear. In this paper, we investigate the effect of the cortical growth, the initial geometry and the initial cortical thickness on folding patterns. In addition, we not only use several descriptors of the folds such as the dimensionless mean curvature, the surface-based three-dimensional gyrification index and the sulcal depth, but also propose a new metric to quantify the folds orientation. The results demonstrate that the cortical growth mode does almost not affect the complexity degree of surface morphology; the variation in the initial geometry changes the folds orientation and depth, and in particular, the slenderer the shape is, the more folds along its longest axis could be seen and the deeper the sulci become. Moreover, the thinner the initial cortical thickness is, the higher the spatial frequency of the folds is, but the shallower the sulci become, which is in agreement with the previously reported effects of cortical thickness.