In this study, we introduce a semi-classical finite-difference time-domain (FDTD) framework for modeling molecule-plasmon interactions in three dimensions. By integrating molecular rate equations with electromagnetic field calculations on the Yee grid, the approach provides a self-consistent analysis of the dynamic coupling between molecular processes and plasmonic fields.