A Dynamic Stochastic Model for DNA Replication Initiation in Early Embryos

Background: Eukaryotic cells seem unable to monitor replication completion during normal S phase, yet must ensure a reliable replication completion time. This is an acute problem in early Xenopus embryos since DNA replication origins are located and activated stochastically, leading to the random completion problem. DNA combing, kinetic modelling and other studies using Xenopus egg extracts have suggested that potential origins are much more abundant than actual initiation events and that the time-dependent rate of initiation, I(t), markedly increases through S phase to ensure the rapid completion of unreplicated gaps and a narrow distribution of completion times. However, the molecular mechanism that underlies this increase has remained obscure.Methodology/Principal Findings: Using both previous and novel DNA combing data we have confirmed that I(t) increases through S phase but have also established that it progressively decreases before the end of S phase. To explore plausible biochemical scenarios that might explain these features, we have performed comparisons between numerical simulations and DNA combing data. Several simple models were tested: i) recycling of a limiting replication fork component from completed replicons; ii) time-dependent increase in origin efficiency; iii) time-dependent increase in availability of an initially limiting factor, e. g. by nuclear import. None of these potential mechanisms could on its own account for the data. We propose a model that combines time-dependent changes in availability of a replication factor and a fork-density dependent affinity of this factor for potential origins. This novel model quantitatively and robustly accounted for the observed changes in initiation rate and fork density.Conclusions/Significance: This work provides a refined temporal profile of replication initiation rates and a robust, dynamic model that quantitatively explains replication origin usage during early embryonic S phase. These results have significant implications for the organisation of replication origins in higher eukaryotes

Modeling Inhomogeneous DNA Replication Kinetics

In eukaryotic organisms, DNA replication is initiated at a series of chromosomal locations called origins, where replication forks are assembled proceeding bidirectionally to replicate the genome. The distribution and firing rate of these origins, in conjunction with the velocity at which forks progress, dictate the program of the replication process. Previous attempts at modeling DNA replication in eukaryotes have focused on cases where the firing rate and the velocity of replication forks are homogeneous, or uniform, across the genome. However, it is now known that there are large variations in origin activity along the genome and variations in fork velocities can also take place. Here, we generalize previous approaches to modeling replication, to allow for arbitrary spatial variation of initiation rates and fork velocities. We derive rate equations for left- and right-moving forks and for replication probability over time that can be solved numerically to obtain the mean-field replication program. This method accurately reproduces the results of DNA replication simulation. We also successfully adapted our approach to the inverse problem of fitting measurements of DNA replication performed on single DNA molecules. Since such measurements are performed on specified portion of the genome, the examined DNA molecules may be replicated by forks that originate either within the studied molecule or outside of it. This problem was solved by using an effective flux of incoming replication forks at the model boundaries to represent the origin activity outside the studied region. Using this approach, we show that reliable inferences can be made about the replication of specific portions of the genome even if the amount of data that can be obtained from single-molecule experiments is generally limited

DNA replication timing is deterministic at the level of chromosomal domains but stochastic at the level of replicons in Xenopus egg extracts

Replication origins in Xenopus egg extracts are located at apparently random sequences but are activated in clusters that fire at different times during S phase under the control of ATR/ATM kinases. We investigated whether chromosomal domains and single sequences replicate at distinct times during S phase in egg extracts. Replication foci were found to progressively appear during early S phase and foci labelled early in one S phase colocalized with those labelled early in the next S phase. However, the distribution of these two early labels did not coincide between single origins or origin clusters on single DNA fibres. The 4 Mb Xenopus rDNA repeat domain was found to replicate later than the rest of the genome and to have a more nuclease-resistant chromatin structure. Replication initiated more frequently in the transcription unit than in the intergenic spacer. These results suggest for the first time that in this embryonic system, where transcription does not occur, replication timing is deterministic at the scale of large chromatin domains (1–5 Mb) but stochastic at the scale of replicons (10 kb) and replicon clusters (50–100 kb)

Computational Methods to Study Kinetics of DNA Replication

New technologies such as DNA combing have led to the availability of large quanti-ties of data that describe the state of DNA while undergoing replication in S phase. In this chapter, we describe methods used to extract various parameters of replica-tion — fork velocity, origin initiation rate, fork density, numbers of potential and utilized origins — from such data. We first present a version of the technique that applies to “ideal ” data. We then show how to deal with a number of real-world complications, such as the asynchrony of starting times of a population of cells, the finite length of fragments used in the analysis, and the finite amount of DNA in a chromosome. Key words: DNA replication, replication fork velocity, origin initiation

Plx1 is required for chromosomal DNA replication under stressful conditions

Polo-like kinase (Plk)1 is required for mitosis progression. However, although Plk1 is expressed throughout the cell cycle, its function during S-phase is unknown. Using Xenopus laevis egg extracts, we demonstrate that Plx1, the Xenopus orthologue of Plk1, is required for DNA replication in the presence of stalled replication forks induced by aphidicolin, etoposide or reduced levels of DNA-bound Mcm complexes. Plx1 binds to chromatin and suppresses the ATM/ATR-dependent intra-S-phase checkpoint that inhibits origin firing. This allows Cdc45 loading and derepression of DNA replication initiation. Checkpoint activation increases Plx1 binding to the Mcm complex through its Polo box domain. Plx1 recruitment to chromatin is independent of checkpoint mediators Tipin and Claspin. Instead, ATR-dependent phosphorylation of serine 92 of Mcm2 is required for the recruitment of Plx1 to chromatin and for the recovery of DNA replication under stress. Depletion of Plx1 leads to accumulation of chromosomal breakage that is prevented by the addition of recombinant Plx1. These data suggest that Plx1 promotes genome stability by regulating DNA replication under stressful conditions

Evidence for Sequential and Increasing Activation of Replication Origins along Replication Timing Gradients in the Human Genome

Genome-wide replication timing studies have suggested that mammalian chromosomes consist of megabase-scale domains of coordinated origin firing separated by large originless transition regions. Here, we report a quantitative genome-wide analysis of DNA replication kinetics in several human cell types that contradicts this view. DNA combing in HeLa cells sorted into four temporal compartments of S phase shows that replication origins are spaced at 40 kb intervals and fire as small clusters whose synchrony increases during S phase and that replication fork velocity (mean 0.7 kb/min, maximum 2.0 kb/min) remains constant and narrowly distributed through S phase. However, multi-scale analysis of a genome-wide replication timing profile shows a broad distribution of replication timing gradients with practically no regions larger than 100 kb replicating at less than 2 kb/min. Therefore, HeLa cells lack large regions of unidirectional fork progression. Temporal transition regions are replicated by sequential activation of origins at a rate that increases during S phase and replication timing gradients are set by the delay and the spacing between successive origin firings rather than by the velocity of single forks. Activation of internal origins in a specific temporal transition region is directly demonstrated by DNA combing of the IGH locus in HeLa cells. Analysis of published origin maps in HeLa cells and published replication timing and DNA combing data in several other cell types corroborate these findings, with the interesting exception of embryonic stem cells where regions of unidirectional fork progression seem more abundant. These results can be explained if origins fire independently of each other but under the control of long-range chromatin structure, or if replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA. These findings shed a new light on the replication timing program of mammalian genomes and provide a general model for their replication kinetics

High level production and localization of bovine rod cGMP-gated cation channel subunit 1 in baculovirus-infected insect cells and Saccharomyces cerevisiae

The cGMP-gated cation channel is responsible for the last step in the vertebrate phototransduction cascade which couples light activation of rhodopsin to a change in membrane permeability. Two different expression systems, baculovirus-infected Spodoptera frugiperda (Sf9) insect cells and Saccharomyces cerevisiae, were used for the overproduction of the cGMP-gated cation channel subunit 1 of bovine rod cells. Presence of recombinant channel protein was monitored by SDS-PAGE and Western-blot analysis. Through immunogold labeling, the heterologously expressed cation channel was found to be localized predominantly in the inner compartments of the infected insect cells and in the vacuole of recombinant yeast cells