Our model adopts the Gilliespie algorithm to randomly simulate from the microscopic perspective of molecules, so as to reflect the actual situation in the most real way. The specific practice is divided into the following steps: 1) For all reactions that change the amount and state of the effective substance in the system, such as the production of each protein, the interaction of each protein, and the modification of DNA by protein or small molecules. Each chemical reaction is analyzed separately. Conditions, rates, changes, etc. under different conditions are set for each reaction in combination with literature and experimental data. 2) Determine the amount of each substance at the initial moment of the system, the parameters of each reaction obtained from part I, and calculate each time interval of the system through conditional probability formula and multiplication rule. Also, calculate the time of the next occurrence of the system at each reaction time, so as to obtain the amount of different time points and different effective substances.