Bone remodeling process has been used in orthodontics to treat malposition of teeth in patients by applying stimuli outside of usual everyday loads to promote tooth movement by affecting equilibrium state of the surrounding bone tissue. Accurate modeling of long term orthodontic tooth movement (OTM) is crucial in the field of dental biomechanical research since it allows to predict the behavior and interaction of bone-tooth environment in a non-destructive way, and helps to gain more insight on how exactly tooth motion progresses over time. Existence of such predictive tools might help to avoid the adverse effects of OTM on teeth and the surrounding tissues during this clinical procedure. In this study a new numerical approach to simulating long-term OTM is proposed, that involves external bone adaptation with strain energy density of the bone taken as the stimulus parameter and bone adaptation modeled by nodal movements at the bone-tooth interface using Abaqus UMESHMOTION subroutine. Contrary to conventional re-meshing algorithms, where the mesh of resorbed-apposed bone region is constantly updated and element deletion/creation is performed for each increment, the proposed method only moves nodes without changing the initial mesh topology. For this study, a 3D model of right central maxillary incisor tooth and its surrounding maxillary bone was used for the modeling of OTM for a duration of 1 week. Two test cases were performed and the results from induced tooth motion were investigated. Results indicate tooth movement values that were quite close to clinical values provided in the literature and this method is easily applicable to validate various postulates of OTM via adapting the stimulus-adaption rate relation and patient-specific planning of orthodontic patients as well.