Hashimoto thyroiditis (HT) is a pathology that often causes a gradual thyroid insufficiency in affected patients due to the autoimmune destruction of this gland. The cellular immune response mediated by T helper lymphocytes TH1 and TH17 can induce the HT disease. In this pathologic condition, there is an imbalance between the TH17 and Treg lymphocytes as well as a gut microbiota dysfunction. The objective of this work was to describe the interactions of the cell subpopulations that participate in HT. To achieve this goal, we generated a mathematical model that allowed the simulation of different scenarios for the dynamic interaction between thyroid cells, the immune system, and the gut microbiota. We used a hypothetical-deductive design of mathematical modeling based on a system of ordinary differential equations, where the state variables are the TH1, TH17, and Treg lymphocytes, the thyrocytes, and the bacteria from gut microbiota. This work generated a compartmental model of the cellular immune response occurring in the thyroid gland. It was observed that TH1 and TH17 lymphocytes could increase the immune cells’ activity, as well as activate effector cells directly and trigger the apoptosis and inflammation processes of healthy thyrocytes indirectly. Likewise, the model showed that a reduction in Treg lymphocytes could increase the activity of TH17 lymphocytes when an imbalance of the gut microbiota composition occurred. The numerical results highlight the TH1, TH17, and bacterial balance of the gut microbiota activities as important factors for the development of HT disease.