Discovery of nonsteroidal anti-inflammatory drug and anticancer drug enhancing reprogramming and induced pluripotent stem cell generation.

Recent breakthroughs in creating induced pluripotent stem cells (iPSCs) provide alternative means to obtain embryonic stem-like cells without destroying embryos by introducing four reprogramming factors (Oct3/4, Sox2, and Klf4/c-Myc or Nanog/Lin28) into somatic cells. iPSCs are versatile tools for investigating early developmental processes and could become sources of tissues or cells for regenerative therapies. Here, for the first time, we describe a strategy to analyze genomics datasets of mouse embryonic fibroblasts (MEFs) and embryonic stem cells to identify genes constituting barriers to iPSC reprogramming. We further show that computational chemical biology combined with genomics analysis can be used to identify small molecules regulating reprogramming. Specific downregulation by small interfering RNAs (siRNAs) of several key MEF-specific genes encoding proteins with catalytic or regulatory functions, including WISP1, PRRX1, HMGA2, NFIX, PRKG2, COX2, and TGFß3, greatly increased reprogramming efficiency. Based on this rationale, we screened only 17 small molecules in reprogramming assays and discovered that the nonsteroidal anti-inflammatory drug Nabumetone and the anticancer drug 4-hydroxytamoxifen can generate iPSCs without Sox2. Nabumetone could also produce iPSCs in the absence of c-Myc or Sox2 without compromising self-renewal and pluripotency of derived iPSCs. In summary, we report a new concept of combining genomics and computational chemical biology to identify new drugs useful for iPSC generation. This hypothesis-driven approach provides an alternative to shot-gun screening and accelerates understanding of molecular mechanisms underlying iPSC induction.

Cyclin E and histone H3 levels are regulated by 5-fluorouracil in a DNA mismatch repair-dependent manner.

Several studies indicate that the DNA mismatch repair (MMR) system may trigger cytotoxicity upon 5-fluorouracil (5-FU) recognition, but signaling pathways regulated by MMR in response to 5-FU are unknown. We hypothesize that recognition of 5-FU in DNA by MMR proteins trigger specific signaling cascades that results in slowing of the cell cycle and cell death. Whole human genome cDNA microarrays were used to examine relative signaling responses induced in MMR-proficient cells after 5-FU (5 µM) treatment for 24 hours. Analysis revealed 43 pathways differentially affected by 5-FU compared to control (P 1.4-fold) and downregulated cdc25C, cyclins B1 and B2, histone H2A, H2B, and H3 (< -1.4-fold) over control. Cell cycle analysis revealed a G1/S arrest by 5-FU that was congruent with increased cyclin E and decreased cdc25C protein expression. Importantly, with knockdown of hMLH1 and hMSH2, we observed that decreased histone H3 expression by 5-FU was dependent on hMLH1. Additionally, 5-FU treatment dramatically decreased levels of several histone H3 modifications. Our data suggest that 5-FU induces a G1/S arrest by regulating cyclin E and cdc25C expression, and MMR recognition of 5-FU in DNA may modulate cyclin E to affect the cell cycle. Furthermore, MMR recognition of 5-FU reduces histone H3 levels that could be related to DNA access by proteins and/or cell death during the G1/S phase of the cell cycle.

Both microsatellite length and sequence context determine frameshift mutation rates in defective DNA mismatch repair.

It is generally accepted that longer microsatellites mutate more frequently in defective DNA mismatch repair (MMR) than shorter microsatellites. Indeed, we have previously observed that the A10 microsatellite of transforming growth factor beta type II receptor (TGFBR2) frameshifts -1 bp at a faster rate than the A8 microsatellite of activin type II receptor (ACVR2), although both genes become frameshift-mutated in >80% of MMR-defective colorectal cancers. To experimentally determine the effect of microsatellite length upon frameshift mutation in gene-specific sequence contexts, we altered the microsatellite length within TGFBR2 exon 3 and ACVR2 exon 10, generating A7, A10 and A13 constructs. These constructs were cloned 1 bp out of frame of EGFP, allowing a -1 bp frameshift to drive EGFP expression, and stably transfected into MMR-deficient cells. Subsequent non-fluorescent cells were sorted, cultured for 7-35 days and harvested for EGFP analysis and DNA sequencing. Longer microsatellites within TGFBR2 and ACVR2 showed significantly higher mutation rates than shorter ones, with TGFBR2 A13, A10 and A7 frameshifts measured at 22.38x10(-4), 2.17x10(-4) and 0.13x10(-4), respectively. Surprisingly, shorter ACVR2 constructs showed three times higher mutation rates at A7 and A10 lengths than identical length TGFBR2 constructs but comparably lower at the A13 length, suggesting influences from both microsatellite length as well as the sequence context. Furthermore, the TGFBR2 A13 construct mutated into 33% A11 sequences (-2 bp) in addition to expected A12 (-1 bp), indicating that this construct undergoes continual subsequent frameshift mutation. These data demonstrate experimentally that both the length of a mononucleotide microsatellite and its sequence context influence mutation rate in defective DNA MMR.

Mutation rates of TGFBR2 and ACVR2 coding microsatellites in human cells with defective DNA mismatch repair.

Microsatellite instability promotes colonic tumorigenesis through generating frameshift mutations at coding microsatellites of tumor suppressor genes, such as TGFBR2 and ACVR2. As a consequence, signaling through these TGFbeta family receptors is abrogated in DNA Mismatch repair (MMR)-deficient tumors. How these mutations occur in real time and mutational rates of these human coding sequences have not previously been studied. We utilized cell lines with different MMR deficiencies (hMLH1-/-, hMSH6-/-, hMSH3-/-, and MMR-proficient) to determine mutation rates. Plasmids were constructed in which exon 3 of TGFBR2 and exon 10 of ACVR2 were cloned +1 bp out of frame, immediately after the translation initiation codon of an enhanced GFP (EGFP) gene, allowing a -1 bp frameshift mutation to drive EGFP expression. Mutation-resistant plasmids were constructed by interrupting the coding microsatellite sequences, preventing frameshift mutation. Stable cell lines were established containing portions of TGFBR2 and ACVR2, and nonfluorescent cells were sorted, cultured for 7-35 days, and harvested for flow cytometric mutation detection and DNA sequencing at specific time points. DNA sequencing revealed a -1 bp frameshift mutation (A9 in TGFBR2 and A7 in ACVR2) in the fluorescent cells. Two distinct fluorescent populations, M1 (dim, representing heteroduplexes) and M2 (bright, representing full mutants) were identified, with the M2 fraction accumulating over time. hMLH1 deficiency revealed 11 (5.91 x 10(-4)) and 15 (2.18 x 10(-4)) times higher mutation rates for the TGFBR2 and ACVR2 microsatellites compared to hMSH6 deficiency, respectively. The mutation rate of the TGFBR2 microsatellite was approximately 3 times higher in both hMLH1 and hMSH6 deficiencies than the ACVR2 microsatellite. The -1 bp frameshift mutation rates of TGFBR2 and ACVR2 microsatellite sequences are dependent upon the human MMR background.