Methods and systems for determining orthodontic treatment based on 3D digital model of subject's teeth. The method includes determining, for a given trajectory segment of the orthodontic treatment of a given tooth in which the given tooth is moved from a start position to an end position, whether the force required to move the given tooth from the start position to the end position causes an induced stress associated with at least one other of the plurality of subject's teeth. In response to the induced stress being outside of a predetermined threshold level, determining a counter force to be applied to the at least one other of the plurality of subject's teeth. The orthodontic treatment can be thus determined as including the force to be applied to the given tooth and the determined counter force to be applied to the at least one other of the plurality if subject's teeth.