S-phase Duration Research Articles

Determination of S-Phase Duration Using 5-Ethynyl-2'-deoxyuridine Incorporation in Saccharomyces cerevisiae.

Eukaryotic DNA replication is a highly regulated process that ensures that the genetic blueprint of a cell is correctly duplicated prior to chromosome segregation. As DNA synthesis defects underlie chromosome rearrangements, monitoring DNA replication has become essential to understand the basis of genome instability. Saccharomyces cerevisiae is a classical model to study cell cycle regulation, but key DNA replication parameters, such as the fraction of cells in the S phase or the S-phase duration, are still difficult to determine. This protocol uses short and non-toxic pulses of 5-ethynyl-2'-deoxyuridine (EdU), a thymidine analog, in engineered TK-hENT1 yeast cells, followed by its detection by Click reaction to allow the visualization and quantification of DNA replication with high spatial and temporal resolution at both the single-cell and population levels by microscopy and flow cytometry. This method may identify previously overlooked defects in the S phase and cell cycle progression of yeast mutants, thereby allowing the characterization of new players essential for ensuring genome stability.

Measuring S-Phase Duration from Asynchronous Cells Using Dual EdU-BrdU Pulse-Chase Labeling Flow Cytometry

Eukaryotes duplicate their chromosomes during the cell cycle S phase using thousands of initiation sites, tunable fork speed and megabase-long spatio-temporal replication programs. The duration of S phase is fairly constant within a given cell type, but remarkably plastic during development, cell differentiation or various stresses. Characterizing the dynamics of S phase is important as replication defects are associated with genome instability, cancer and ageing. Methods to measure S-phase duration are so far indirect, and rely on mathematical modelling or require cell synchronization. We describe here a simple and robust method to measure S-phase duration in cell cultures using a dual EdU-BrdU pulse-labeling regimen with incremental thymidine chases, and quantification by flow cytometry of cells entering and exiting S phase. Importantly, the method requires neither cell synchronization nor genome engineering, thus avoiding possible artifacts. It measures the duration of unperturbed S phases, but also the effect of drugs or mutations on it. We show that this method can be used for both adherent and suspension cells, cell lines and primary cells of different types from human, mouse and Drosophila. Interestingly, the method revealed that several commonly-used cancer cell lines have a longer S phase compared to untransformed cells.

Increased expression of Polδ does not alter the canonical replication program in vivo.

Background: In vitro experiments utilising the reconstituted Saccharomyces cerevisiae eukaryotic replisome indicated that the efficiency of the leading strand replication is impaired by a moderate increase in Polδ concentration. It was hypothesised that the slower rate of the leading strand synthesis characteristic for reactions containing two-fold and four-fold increased concentration of Polδ represented a consequence of a relatively rare event, during which Polδ stochastically outcompeted Polε and, in an inefficient manner, temporarily facilitated extension of the leading strand. Inspired by this observation, we aimed to determine whether similarly increased Polδ levels influence replication dynamics in vivo using the fission yeast Schizosaccharomyces pombe as a model system. Methods: To generate S. pombe strains over-expressing Polδ, we utilised Cre-Lox mediated cassette exchange and integrated one or three extra genomic copies of all four Polδ genes. To estimate expression of respective Polδ genes in Polδ-overexpressing mutants, we measured relative transcript levels of cdc1 + , cdc6 + (or cdc6 L591G ), cdc27 + and cdm1 + by reverse transcription followed by quantitative PCR (RT-qPCR). To assess the impact of Polδ over-expression on cell physiology and replication dynamics, we used standard cell biology techniques and polymerase usage sequencing. Results: We provide an evidence that two-fold and four-fold over-production of Polδ does not significantly alter growth rate, cellular morphology and S-phase duration. Polymerase usage sequencing analysis further indicates that increased Polδ expression does not change activities of Polδ, Polε and Polα at replication initiation sites and across replication termination zones. Additionally, we show that mutants over-expressing Polδ preserve WT-like distribution of replication origin efficiencies. Conclusions: Our experiments do not disprove the existence of opportunistic polymerase switches; however, the data indicate that, if stochastic replacement of Polε for Polδ does occur i n vivo, it represents a rare phenomenon that does not significantly influence canonical replication program.

Increased expression of Polδ does not alter the canonical replication program in vivo

Background: In vitro experiments utilising the reconstituted Saccharomyces cerevisiae eukaryotic replisome indicated that the efficiency of the leading strand replication is impaired by a moderate increase in Polδ concentration. It was hypothesised that the slower rate of the leading strand synthesis characteristic for reactions containing two-fold and four-fold increased concentration of Polδ represented a consequence of a relatively rare event, during which Polδ stochastically outcompeted Polε and, in an inefficient manner, temporarily facilitated extension of the leading strand. Inspired by this observation, we aimed to determine whether similarly increased Polδ levels influence replication dynamics in vivo using the fission yeast Schizosaccharomyces pombe as a model system. Methods: To generate S. pombe strains over-expressing Polδ, we utilised Cre-Lox mediated cassette exchange and integrated one or three extra genomic copies of all four Polδ genes. To estimate expression of respective Polδ genes in Polδ-overexpressing mutants, we measured relative transcript levels of cdc1 +, cdc6 + (or cdc6 L591G), cdc27 + and cdm1 + by reverse transcription followed by quantitative PCR (RT-qPCR). To assess the impact of Polδ over-expression on cell physiology and replication dynamics, we used standard cell biology techniques and polymerase usage sequencing. Results: We provide an evidence that two-fold and four-fold over-production of Polδ does not significantly alter growth rate, cellular morphology and S-phase duration. Polymerase usage sequencing analysis further indicates that increased Polδ expression does not change activities of Polδ, Polε and Polα at replication initiation sites and across replication termination zones. Additionally, we show that mutants over-expressing Polδ preserve WT-like distribution of replication origin efficiencies. Conclusions: Our experiments do not disprove the existence of opportunistic polymerase switches; however, the data indicate that, if stochastic replacement of Polε for Polδ does occur i n vivo, it represents a rare phenomenon that does not significantly influence canonical replication program.

Increased expression of Polδ does not alter the canonical replication program in vivo.

Background: In vitro experiments utilising the reconstituted Saccharomyces cerevisiae eukaryotic replisome indicated that the efficiency of the leading strand replication is impaired by a moderate increase in Polδ concentration. It was hypothesised that the slower rate of the leading strand synthesis characteristic for reactions containing two-fold and four-fold increased concentration of Polδ represented a consequence of a relatively rare event, during which Polδ stochastically outcompeted Polε and, in an inefficient manner, temporarily facilitated extension of the leading strand. Inspired by this observation, we aimed to determine whether similarly increased Polδ levels influence replication dynamics in vivo using the fission yeast Schizosaccharomyces pombe as a model system. Methods: To generate S. pombe strains over-expressing Polδ, we utilised Cre-Lox mediated cassette exchange and integrated one or three extra genomic copies of all four Polδ genes. To estimate expression of respective Polδ genes in Polδ-overexpressing mutants, we measured relative transcript levels of cdc1 + , cdc6 + (or cdc6 L591G ), cdc27 + and cdm1 + by reverse transcription followed by quantitative PCR (RT-qPCR). To assess the impact of Polδ over-expression on cell physiology and replication dynamics, we used standard cell biology techniques and polymerase usage sequencing. Results: We provide an evidence that two-fold and four-fold over-production of Polδ does not significantly alter growth rate, cellular morphology and S-phase duration. Polymerase usage sequencing analysis further indicates that increased Polδ expression does not change activities of Polδ, Polε and Polα at replication initiation sites and across replication termination zones. Additionally, we show that mutants over-expressing Polδ preserve WT-like distribution of replication origin efficiencies. Conclusions: Our experiments do not disprove the existence of opportunistic polymerase switches; however, the data indicate that, if stochastic replacement of Polε for Polδ does occur i n vivo, it represents a rare phenomenon that does not significantly influence canonical replication program.

PRIM1 deficiency causes a distinctive primordial dwarfism syndrome.

DNA replication is fundamental for cell proliferation in all organisms. Nonetheless, components of the replisome have been implicated in human disease, and here we report PRIM1 encoding the catalytic subunit of DNA primase as a novel disease gene. Using a variant classification agnostic approach, biallelic mutations in PRIM1 were identified in five individuals. PRIM1 protein levels were markedly reduced in patient cells, accompanied by replication fork asymmetry, increased interorigin distances, replication stress, and prolonged S-phase duration. Consequently, cell proliferation was markedly impaired, explaining the patients' extreme growth failure. Notably, phenotypic features distinct from those previously reported with DNA polymerase genes were evident, highlighting differing developmental requirements for this core replisome component that warrant future investigation.

Sonic hedgehog accelerates DNA replication to cause replication stress promoting cancer initiation in medulloblastoma.

The mechanisms generating cancer-initiating mutations are not well understood. Sonic hedgehog (SHH) pathway activation is frequent in medulloblastoma (MB), with PTCH1 mutations being a common initiating event. Here we investigated the role of the developmental mitogen SHH in initiating carcinogenesis in the cells of origin: granule cell progenitors (GCPs). We delineate a molecular mechanism for tumor initiation in MB. Exposure of GCPs to Shh causes a distinct form of DNA replication stress, increasing both origin firing and fork velocity. Shh promotes DNA helicase loading and activation, with increased Cdc7-dependent origin firing. The S-phase duration is reduced and hyper-recombination occurs, causing copy number neutral loss of heterozygosity-a frequent event at the PTCH1/ptch1 locus. Moreover, Cdc7 inhibition to attenuate origin firing reduces recombination and preneoplastic tumor formation in mice. Therefore, tissue-specific replication stress induced by Shh promotes loss of heterozygosity, which in tumor-prone Ptch1+/- GCPs results in loss of this tumor suppressor-an early cancer-initiating event.

Estimation of the Minimum Number of Replication Origins Per Chromosome in any Organism.

Eukaryote nuclear genomes predominantly replicate through multiple replication origins. The number of replication origins activated per chromosome during the S-phase duration may vary according to many factors, but the predominant one is replication stress. Several studies have applied different approaches to estimate the number and map the positions of the replication origins in various organisms. However, without a parameter to restrict the minimum of necessary origins, less sensitive techniques may suggest conflicting results. The estimation of the minimum number of replication origins (MO) per chromosome is an innovative method that allows the establishment of a threshold, which serves as a parameter for genomic approaches that map origins. For this, the MO can be easily obtained through a formula that requires as parameters: chromosome size, S-phase duration, and replication rate. The chromosome size for any organism can be acquired in genomic databanks (such as NCBI), the S-phase duration can be estimated by monitoring DNA replication, and the replication rate is obtained through the DNA combing approach. The estimation of MO is a simple, quick, and easy method that provides a new methodological framework to assist studies of mapping replication origins in any organism.

Transcription activity contributes to the firing of non-constitutive origins in African trypanosomes helping to maintain robustness in S-phase duration

The co-synthesis of DNA and RNA potentially generates conflicts between replication and transcription, which can lead to genomic instability. In trypanosomatids, eukaryotic parasites that perform polycistronic transcription, this phenomenon and its consequences are still little studied. Here, we showed that the number of constitutive origins mapped in the Trypanosoma brucei genome is less than the minimum required to complete replication within S-phase duration. By the development of a mechanistic model of DNA replication considering replication-transcription conflicts and using immunofluorescence assays and DNA combing approaches, we demonstrated that the activation of non-constitutive (backup) origins are indispensable for replication to be completed within S-phase period. Together, our findings suggest that transcription activity during S phase generates R-loops, which contributes to the emergence of DNA lesions, leading to the firing of backup origins that help maintain robustness in S-phase duration. The usage of this increased pool of origins, contributing to the maintenance of DNA replication, seems to be of paramount importance for the survival of this parasite that affects million people around the world.

Replication timing in Drosophila and its peculiarities in polytene chromosomes

Drosophila melanogaster is one of the popular model organisms in DNA replication studies. Since the 1960s, DNA replication of polytene chromosomes has been extensively studied by cytological methods. In the recent two decades, the progress in our understanding of DNA replication was associated with new techniques. Use of fluorescent dyes increased the resolution of cytological methods significantly. High-throughput methods allowed analysis of DNA replication on a genome scale, as well as its correlation with chromatin structure and gene activi ty. Precise mapping of the cytological structures of polytene chromosomes to the genome assembly allowed comparison of replication between polytene chromosomes and chromosomes of diploid cells. New features of replication characteristic for D. melanogaster were described for both diploid and polytene chromosomes. Comparison of genomic replication profiles revealed a significant similarity between Drosophila and other well-studi ed eukaryotic species, such as human. Early replication is often confined to intensely transcribed gene-dense regions characterized by multiple replication initiation sites. Features of DNA replication in Drosophila might be explained by a compact genome. The organization of replication in polytene chromosomes has much in common with the organization of replication in chromosomes in diploid cells. The most important feature of replication in polytene chromosomes is its low rate and the dependence of S-phase duration on many factors: external and internal, local and global. The speed of replication forks in D. melanogaster polytene chromosomes is affected by SUUR and Rif1 proteins. It is not known yet how universal the mechanisms associated with these factors are, but their study is very promising.

Cell cycle and differentiation of Sca-1+ and Sca-1− hematopoietic stem and progenitor cells

ABSTRACTHematopoietic stem and progenitor cells (HSPCs) are crucial for lifelong blood cell production. We analyzed the cell cycle and cell production rate in HSPCs in murine hematopoiesis. The labeling of DNA-synthesizing cells by two thymidine analogues, optimized for in-vivo use, enabled determination of the cell cycle flow rate into G2-phase, the duration of S-phase and the average cell cycle time in Sca-1+ and Sca-1− HSPCs. Determination of cells with 2n DNA content labeled in preceding S-phase was then used to establish the cell flow rates in G1-phase. Our measurements revealed a significant difference in how Sca-1+ and Sca-1− myeloid progenitors self-renew and differentiate. Division of the Sca-1+ progenitors led to loss of the Sca-1 marker in about half of newly produced cells, corresponding to asymmetric cell division. Sca-1− cells arising from cell division entered a new round of the cell cycle, corresponding to symmetric self-renewing cell division. The novel data also enabled the estimation of the cell production rates in Sca-1+ and in three subtypes of Sca-1− HSPCs and revealed Sca-1 negative cells as the major amplification stage in the blood cell development.

Hierarchy of Sca-1positive and Sca-1negative hematopoietic stem and progenitor cells and their contribution to blood cell production

The hematopoiesis steadily generates a large number of blood cells by intensive cell proliferation. A hierarchy of hematopoietic stem and progenitor cells (HSPCs) has been well characterized in mouse. The aim of our study was to find out to which extent these distinct populations contributed to blood cell production. We estimated the production rates of various types of HSPCs according to the percentage of S phase cells, their number in bone marrow, and the S-phase duration. Flow cytometry was used for immunophenotyping of HSPCs according to Lin markers, c-Kit, Sca-1, CD48, CD150, CD34, CD16/32, CD71, and IL7R expression. The dual BrdU/EdU labelling technique was applied to determine the S-phase duration and the percentage of DNA-synthesizing (S-phase) HSPCs. The HSPC apoptotic rate was detected by Ghost Dye™ Red 780 staining. The S-phase fraction significantly differed in various types of HSPCs. The values ranged from ≈ 1 % to ≈ 90 % of cells in S-phase, and were a characterizing feature of particular types of HSPCs. The CD71 (transferrin receptor) expression level closely correlated with the proliferation rate in all types of HSPCs. The highest S-phase fraction was found in the Lin c Kit+Sca-1- progenitors committed to the megakaryocyte erythroid development, followed by the progenitors committed to the granulocyte macrophage development. The HSPCs lacking the Sca-1 antigen generated approximately 22 times more cells than the Sca 1 expressing ones. Research was supported by the Grant Agency of the Czech Republic (GACR 17-01897S).

Ensa controls S-phase length by modulating Treslin levels

The Greatwall/Ensa/PP2A-B55 pathway is essential for controlling mitotic substrate phosphorylation and mitotic entry. Here, we investigate the effect of the knockdown of the Gwl substrate, Ensa, in human cells. Unexpectedly, Ensa knockdown promotes a dramatic extension of S phase associated with a lowered density of replication forks. Notably, Ensa depletion results in a decrease of Treslin levels, a pivotal protein for the firing of replication origins. Accordingly, the extended S phase in Ensa-depleted cells is completely rescued by the overexpression of Treslin. Our data herein reveal a new mechanism by which normal cells regulate S-phase duration by controlling the ubiquitin-proteasome degradation of Treslin in a Gwl/Ensa-dependent pathway.

Cell-to-cell variability and robustness in S-phase duration from genome replication kinetics.

Genome replication, a key process for a cell, relies on stochastic initiation by replication origins, causing a variability of replication timing from cell to cell. While stochastic models of eukaryotic replication are widely available, the link between the key parameters and overall replication timing has not been addressed systematically. We use a combined analytical and computational approach to calculate how positions and strength of many origins lead to a given cell-to-cell variability of total duration of the replication of a large region, a chromosome or the entire genome. Specifically, the total replication timing can be framed as an extreme-value problem, since it is due to the last region that replicates in each cell. Our calculations identify two regimes based on the spread between characteristic completion times of all inter-origin regions of a genome. For widely different completion times, timing is set by the single specific region that is typically the last to replicate in all cells. Conversely, when the completion time of all regions are comparable, an extreme-value estimate shows that the cell-to-cell variability of genome replication timing has universal properties. Comparison with available data shows that the replication program of three yeast species falls in this extreme-value regime.

A flow cytometric method for estimating S-phase duration in plants

The duration of the DNA synthesis stage (S phase) of the cell cycle is fundamental in our understanding of cell cycle kinetics, cell proliferation, and DNA replication timing programs. Most S-phase duration estimates that exist for plants are based on indirect measurements. We present a method for directly estimating S-phase duration by pulse-labeling root tips or actively dividing suspension cells with the halogenated thymidine analog 5-ethynl-2'-deoxyuridine (EdU) and analyzing the time course of replication with bivariate flow cytometry. The transition between G1 and G2 DNA contents can be followed by measuring the mean DNA content of EdU-labeled S-phase nuclei as a function of time after the labeling pulse. We applied this technique to intact root tips of maize (Zea mays L.), rice (Oryza sativa L.), barley (Hordeum vulgare L.), and wheat (Triticum aestivum L.), and to actively dividing cell cultures of Arabidopsis (Arabidopsis thaliana (L.) Heynh.) and rice. Estimates of S-phase duration in root tips were remarkably consistent, varying only by ~3-fold, although the genome sizes of the species analyzed varied >40-fold.

Orderly progression through S-phase requires dynamic ubiquitylation and deubiquitylation of PCNA.

Proliferating-cell nuclear antigen (PCNA) is a DNA sliding clamp with an essential function in DNA replication and a key role in tolerance to DNA damage by ensuring the bypass of lesions. In eukaryotes, DNA damage tolerance is regulated by ubiquitylation of lysine 164 of PCNA through a well-known control mechanism; however, the regulation of PCNA deubiquitylation remains poorly understood. Our work is a systematic and functional study on PCNA deubiquitylating enzymes (DUBs) in Schizosaccharomyces pombe. Our study reveals that the deubiquitylation of PCNA in fission yeast cells is a complex process that requires several ubiquitin proteases dedicated to the deubiquitylation of a specific subnuclear fraction of mono- and di-ubiquitylated PCNA or a particular type of poly-ubiquitylated PCNA and that there is little redundancy among these enzymes. To understand how DUB activity regulates the oscillatory pattern of ubiquitylated PCNA in fission yeast, we assembled multiple DUB mutants and found that a quadruple mutation of ubp2+, ubp12+, ubp15+, and ubp16+ leads to the stable accumulation of mono-, di-, and poly-ubiquitylated forms of PCNA, increases S-phase duration, and sensitizes cells to DNA damage. Our data suggest that the dynamic ubiquitylation and deubiquitylation of PCNA occurs during S-phase to ensure processive DNA replication.

Multiparametric Cell Cycle Analysis Using the Operetta High-Content Imager and Harmony Software with PhenoLOGIC.

High-content imaging is a powerful tool for determining cell phenotypes at the single cell level. Characterising the effect of small molecules on cell cycle distribution is important for understanding their mechanism of action especially in oncology drug discovery but also for understanding potential toxicology liabilities. Here, a high-throughput phenotypic assay utilising the PerkinElmer Operetta high-content imager and Harmony software to determine cell cycle distribution is described. PhenoLOGIC, a machine learning algorithm within Harmony software was employed to robustly separate single cells from cell clumps. DNA content, EdU incorporation and pHH3 (S10) expression levels were subsequently utilised to separate cells into the various phases of the cell cycle. The assay is amenable to multiplexing with an additional pharmacodynamic marker to assess cell cycle changes within a specific cellular sub-population. Using this approach, the cell cycle distribution of γH2AX positive nuclei was determined following treatment with DNA damaging agents. Likewise, the assay can be multiplexed with Ki67 to determine the fraction of quiescent cells and with BrdU dual labelling to determine S-phase duration. This methodology therefore provides a relatively cheap, quick and high-throughput phenotypic method for determining accurate cell cycle distribution for small molecule mechanism of action and drug toxicity studies.

S-phase duration is the main target of cell cycle regulation in neural progenitors of developing ferret neocortex.

ABSTRACTThe evolutionary expansion of the neocortex primarily reflects increases in abundance and proliferative capacity of cortical progenitors and in the length of the neurogenic period during development. Cell cycle parameters of neocortical progenitors are an important determinant of cortical development. The ferret (Mustela putorius furo), a gyrencephalic mammal, has gained increasing importance as a model for studying corticogenesis. Here, we have studied the abundance, proliferation, and cell cycle parameters of different neural progenitor types, defined by their differential expression of the transcription factors Pax6 and Tbr2, in the various germinal zones of developing ferret neocortex. We focused our analyses on postnatal day 1, a late stage of cortical neurogenesis when upper‐layer neurons are produced. Based on cumulative 5‐ethynyl‐2′‐deoxyuridine (EdU) labeling as well as Ki67 and proliferating cell nuclear antigen (PCNA) immunofluorescence, we determined the duration of the various cell cycle phases of the different neocortical progenitor subpopulations. Ferret neocortical progenitors were found to exhibit longer cell cycles than those of rodents and little variation in the duration of G1 among distinct progenitor types, also in contrast to rodents. Remarkably, the main difference in cell cycle parameters among the various progenitor types was the duration of S‐phase, which became shorter as progenitors progressively changed transcription factor expression from patterns characteristic of self‐renewal to those of neuron production. Hence, S‐phase duration emerges as major target of cell cycle regulation in cortical progenitors of this gyrencephalic mammal. J. Comp. Neurol. 524:456–470, 2016. © 2015 The Authors The Journal of Comparative Neurology Published by Wiley Periodicals, Inc.

Deficiency of the Arabidopsis helicase RTEL1 triggers a SOG1-dependent replication checkpoint in response to DNA cross-links.

To maintain genome integrity, DNA replication is executed and regulated by a complex molecular network of numerous proteins, including helicases and cell cycle checkpoint regulators. Through a systematic screening for putative replication mutants, we identified an Arabidopsis thaliana homolog of human Regulator of Telomere Length 1 (RTEL1), which functions in DNA replication, DNA repair, and recombination. RTEL1 deficiency retards plant growth, a phenotype including a prolonged S-phase duration and decreased cell proliferation. Genetic analysis revealed that rtel1 mutant plants show activated cell cycle checkpoints, specific sensitivity to DNA cross-linking agents, and increased homologous recombination, but a lack of progressive shortening of telomeres, indicating that RTEL1 functions have only been partially conserved between mammals and plants. Surprisingly, RTEL1 deficiency induces tolerance to the deoxynucleotide-depleting drug hydroxyurea, which could be mimicked by DNA cross-linking agents. This resistance does not rely on the essential replication checkpoint regulator WEE1 but could be blocked by a mutation in the SOG1 transcription factor. Taken together, our data indicate that RTEL1 is required for DNA replication and that its deficiency activates a SOG1-dependent replication checkpoint.

Optimized Patient- and Leukemia-Specific Chemotherapy Protocols For The Treatment Of Acute Myeloid Leukemia

▪Treatment for Acute Myeloid Leukaemia (AML) with chemotherapy may result in acute and long-term life-threatening complications due to drug toxicity. Only relatively few patient-and leukemia-specific factors are taken into consideration in current protocols and choice of treatment often depends on the treating physician's experience. With the advent of novel treatments and large amounts of patient- and leukemia-specific genomic data, there is a clear need for a systematic approach to the design and execution of chemotherapy regimens. Mathematical modeling is a tool that can be used for the automation of chemotherapy treatment due to its advantages in systematically exploring extensive datasets in order to capture a system's dynamics and subsequently provide better insight for process enhancement. An AML-specific model that combines the leukemia-specific actions of the cell cycle (i.e. drug target) with patient-specific pharmacology of the drugs (pharmacokinetics) was developed. Specifically, it simulates the response of patients with AML undergoing treatment with two standard chemotherapy protocols, one intensive and the other non-intensive: (a) Daunarubicin (DNR) and Cytosine Arabinoside (Ara-C) used in standard intravenous (iv) doses (DA – 3+7 or 3+10) and (b) low dose Ara-C (LDAC) administered subcutaneously (sc). Sensitivity analysis of the model identifies cell cycle times as the critical parameters that control treatment outcome. For model analysis, clinical data from 6 patients who underwent chemotherapy are used for the estimation of cell cycle time distribution. The patient data are comprised of disease characteristics (tumor burden, cell cycle times, normal cell population) as well as patient-specific characteristics (gender, age, weight and height). The estimated mean S-phase duration (Ts) is 15 hrs (range: 9-21 hrs) and mean whole cell cycle duration (Tc) is 47.5 hrs (range: 33-68 hrs). The estimated data reveal a clear relationship of cell cycle times to treatment outcomes. Specifically, low Ts duration combined with high Tc duration indicates worse treatment outcomes, whereas, the reverse combination is indicative of a good response to treatment. In order to improve effectiveness of AML therapy and reduction of toxicity, chemotherapy treatment is presented as an optimal control problem with the main aim of obtaining a treatment schedule which could maximize leukemic cell kill yet minimize death of the normal cell population in the bone marrow. By the end of treatment, the leukemic population should be reduced to a level of approximately 109 cells at which point BM hypoplasia is achieved. Out of the 6 patients studied, 2 patients had a successful treatment with leukemic hypoplasia achieved, 2 had a reduction of leukemic cells without achieving the hypoplasia target and 2 had disease progression on chemotherapy. The optimization algorithm is formulated and solved for all patients for both intensive and non-intensive treatment protocols with maximal and minimal thresholds set for efficacy and toxicity, respectively. For iv Ara-C, total drug administration is set between 50mg – 4000mg with infusion duration between 1 min to 24 hours. The window for DNR dose optimization is stricter due to potential toxic effects and the only independent variable is dose with 30mg – 90mg per infusion. For sc Ara-C, the maximum dose per day is 40mg and doses are permitted up to four times daily for a maximum period of 20 days. Optimization results obtained for the 6 patients indicate that continuous infusions are more effective for leukaemia cell kill than rapid infusions. For non-intensive chemotherapy, 40mg of Ara-C in continuous infusion over 10 days is better than daily divided doses with leukemic hypoplasia achieved for all patient case studies. For the intensive protocol, dose increase of DNR to 90 mg/m2 combined with Ara-C daily infusion is the optimal chemotherapy regimen. Ara-C doses differ between patients and the optimal dose range is between 200 to 250 mg/m2. In summary, this work presents the potential for improved treatment design in AML therapy, dependent on disease and patient characteristics, defined on a case-by-case basis. This design would provide the opportunity to personalize treatment protocols for gold standard intensive and non-intensive therapies as well as for novel drugs through the use of optimization methods. Disclosures:No relevant conflicts of interest to declare.