The objective of this numerical study was to evaluate the stress induced by orthodontic loading in anchorage implants and surrounding tissues. Orthodontic implants were included in this study. 3D geometrical models were constructed and material characteristics were taken from the literature. Finite element models were created based on the geometry and material characteristics of the screws. Orthodontic horizontal loads of 2 N were applied, and the biomechanical parameters were evaluated by colored scales. The highest von Mises values were recorded around the implant neck area and at the bone-implant interface.