Developing a model of DNA replication to be used for Monte Carlo calculations that predict the sizes and shapes of molecules resulting from DNA double-strand breaks induced by X irradiation during DNA synthesis.

A Monte Carlo computer program was written to introduce double-strand breaks (DSBs) randomly into cellular DNA that is configured according to different models of DNA replication. Then, from a review of the literature using DNA fiber autoradiography and other studies relating to rates of replication of DNA that is organized in approximately 3-Mbp regions or bands, a particular model for DNA replication was developed. Using this model, Monte Carlo calculations were made to predict the types and sizes of molecules that would result from introducing DSBs into DNA when synchronous cells are irradiated in the middle of S phase. Then results of the Monte Carlo calculations were compared with migration profiles obtained by pulsed-field gel electrophoresis (PFGE) for molecular size distributions of linear DNA molecules. For these comparisons, CHO cells irradiated in S phase also were pulse-labeled at the time of irradiation with [3H]dThd for 15 min to compare the migration patterns of 3H-labeled replicating DNA with those of the mass of S-phase DNA, measured by imaging with a CCD camera. For the Monte Carlo calculations, we assumed from the reports in the literature that molecules containing replication bubbles with and without forks would be trapped in the PFGE plug. We also assumed that those molecules that are < or = 8 Mbp, both linear and with replication forks, would be released into the lane. However, approximately 75% of the 3H-labeled DNA that is released from the plug migrated much more slowly than linear molecules, which we attributed to the slow migration of 3H-labeled molecules having replication forks not attached to bubbles. The percentages of both mass of S-phase DNA and 3H-labeled replicating DNA released from the plug, as determined by PFGE, were compared with comparable values determined from Monte Carlo calculations. A DNA replication model that provides good agreement between the PFGE results and Monte Carlo calculations is described. Furthermore, Monte Carlo methodology is presented that can be used for comparing data obtained with PFGE with results of Monte Carlo calculations that are based on different models of DNA replication and different assumptions for the migration of various types of replicating molecules.