In this paper, we study the system size expansion of a stochastic model for radiation-induced DNA damage kinetics and repair. In particular, we characterize both the macroscopic deterministic limit and the fluctuation around it. We further show that such fluctuations are Gaussian-distributed. In deriving such results, we provide further insights into the relationship between stochastic and deterministic mathematical models for radiation-induced DNA damage repair. Specifically, we demonstrate how the governing deterministic equations commonly employed in the field arise naturally within the stochastic framework as a macroscopic limit. Additionally, by examining the fluctuations around this macroscopic limit, we uncover deviations from a Poissonian behavior driven by interactions and clustering among DNA damages. Although such behaviors have been empirically observed, our derived results represent the first rigorous derivation that incorporates these deviations from a Poissonian distribution within a mathematical model, eliminating the need for specific ad hoc corrections.