[en] Orthodontic tooth movement is achieved by the process of repeated alveolar bone resorption on the pressure side and new bone formation on the tension side. In order to optimize orthodontic treatment, it is important to identify and study the biological processes involved. This article presents a mechanobiological model using partial differential equations to describe cell densities, growth factor concentrations, and matrix densities occurring during orthodontic tooth movement. We hypothesize that such a model can predict tooth movement based on the mechanobiological activity of cells in the PDL. The developed model consists of nine coupled non-linear partial differential equations, and two distinct signaling pathways were modeled: the RANKL-RANK-OPG pathway regulating the communication between osteoblasts and osteoclasts and the TGF-beta pathway mediating the differentiation of mesenchymal stem cells into osteoblasts. The predicted concentrations and densities were qualitatively validated by comparing the results to experiments reported in the literature. In the current form, the model supports our hypothesis, as it is capable of conceptually simulating important features of the biological interactions in the alveolar bone-PDL complex during orthodontic tooth movement.