Background: Thanks to the development of the mathematical/statistical reverse engineering and the high-throughput measuring biotechnology, lots of biologically meaningful gene-gene interaction networks have been revealed. Steady-state analysis of these systems provides an important clue to understand and to predict the systematic behaviours of the biological system. However, modeling such a complex and large-scale system is one of the challenging difficulties in systems biology. Results: We introduce a new stochastic modeling approach that can describe gene regulatory mechanisms by dividing two (DNA and protein) layers. Simple queuing system is employed to explain the DNA layer and the protein layer is modeled using G-networks which enable us to account for the post-translational protein interactions. Our method is applied to a transcription repression system and an active protein degradation system. The steadystate results suggest that the active protein degradation system is more sensitive but the transcription repression system might be more reliable than the transcription repression system. Conclusions: Our two layer stochastic model successfully describes the long-run behaviour of gene regulatory networks which consist of various mRNA/protein processes. The analytic solution of the G-networks enables us to extend our model to a large-scale system. A more reliable modeling approach could be achieved by cooperating with a real experimental study in synthetic biology.