The aim of metabolic modeling was to understand the cause effect interaction and reliance linked with the complex interactions of biological networks and molecular systems. Drugs that therapeutically modulate the biological processes of disease are often developed with limited knowledge of the underlying complexity of their specific targets. The robustness for systemic modulating behavior of thyroid hormone secretion during the course of different time unit simulation is explained in this study. In this work, a computational model has been developed which mimics the in vivo simulation. The model was constructed with the help of cell designer 4.1 used for analyzing the effect of perturbed amount of drugs at 0.1, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 amounts as a unit, targeting Thyroid Peroxidase (TPO), Thyroid Stimulating Hormone Receptor (TSHR), and Sodium Iodide Symporter (NIS). The rate kinetic equations were defined with each reaction to simulate the molecular species dynamic behavior. The modulating behavior of thyroid hormone secretion was analyzed by the process of activation and inactivation states of TPO, TSHR, and NIS at various amounts of drugs.Obtained results explain suitably the entire observable fact of the drug effects and are capable to proceed in response to the perturbations of the natural cell. TSHR was found as the most potent molecular therapeutic target in this study.