Structure-specific nucleases in genome dynamics and strategies for targeting cancers.

Nucleases are a super family of enzymes that hydrolyze phosphodiester bonds present in genomes. They widely vary in substrates, causing differentiation in cleavage patterns and having a diversified role in maintaining genetic material. Through cellular evolution of prokaryotic to eukaryotic, nucleases become structure-specific in recognizing its own or foreign genomic DNA/RNA configurations as its substrates, including flaps, bubbles, and Holliday junctions. These special structural configurations are commonly found as intermediates in processes like DNA replication, repair, and recombination. The structure-specific nature and diversified functions make them essential to maintaining genome integrity and evolution in normal and cancer cells. In this article, we review their roles in various pathways, including Okazaki fragment maturation during DNA replication, end resection in homology-directed recombination repair of DNA double-strand breaks, DNA excision repair and apoptosis DNA fragmentation in response to exogeneous DNA damage, and HIV life cycle. As the nucleases serve as key points for the DNA dynamics, cellular apoptosis, and cancer cell survival pathways, we discuss the efforts in the field in developing the therapeutic regimens, taking advantage of recently available knowledge of their diversified structures and functions.

Targeting PARG induces tumor cell growth inhibition and antitumor immune response by reducing phosphorylated STAT3 in ovarian cancer.

BACKGROUND: Ovarian cancer is the most lethal gynecological malignancy, with limited treatment options after failure of standard therapies. Despite the potential of poly(ADP-ribose) polymerase inhibitors in treating DNA damage response (DDR)-deficient ovarian cancer, the development of resistance and immunosuppression limit their efficacy, necessitating alternative therapeutic strategies. Inhibitors of poly(ADP-ribose) glycohydrolase (PARG) represent a novel class of inhibitors that are currently being assessed in preclinical and clinical studies for cancer treatment. METHODS: By using a PARG small-molecule inhibitor, COH34, and a cell-penetrating antibody targeting the PARG's catalytic domain, we investigated the effects of PARG inhibition on signal transducer and activator of transcription 3 (STAT3) in OVCAR8, PEO1, and Brca1-null ID8 ovarian cancer cell lines, as well as in immune cells. We examined PARG inhibition-induced effects on STAT3 phosphorylation, nuclear localization, target gene expression, and antitumor immune responses in vitro, in patient-derived tumor organoids, and in an immunocompetent Brca1-null ID8 ovarian mouse tumor model that mirrors DDR-deficient human high-grade serous ovarian cancer. We also tested the effects of overexpressing a constitutively activated STAT3 mutant on COH34-induced tumor cell growth inhibition. RESULTS: Our findings show that PARG inhibition downregulates STAT3 activity through dephosphorylation in ovarian cancer cells. Importantly, overexpression of a constitutively activated STAT3 mutant in tumor cells attenuates PARG inhibitor-induced growth inhibition. Additionally, PARG inhibition reduces STAT3 phosphorylation in immune cells, leading to the activation of antitumor immune responses, shown in immune cells cocultured with ovarian cancer patient tumor-derived organoids and in immune-competent mice-bearing mouse ovarian tumors. CONCLUSIONS: We have identified a novel antitumor mechanism underlying PARG inhibition beyond its primary antitumor effects through blocking DDR in ovarian cancer. Furthermore, targeting PARG activates antitumor immune responses, thereby potentially increasing response rates to immunotherapy in patients with ovarian cancer.

Effect of genetic distances of different genotypes of maize on the authenticity of single seeds detected by NIR spectroscopy.

Introduction: NIR spectroscopy combined with chemometric algorithms has been widely used for seed authenticity detection. However, the study of seed genetic distance, an internal feature that affects the discriminative performance of classification models, has rarely been reported. Methods: Therefore, maize seed samples of different genotypes were selected to investigate the effect of genetic distance on the authenticity of single seeds detected by NIR spectroscopy. Firstly, the Support vector machine (SVM) model was established using spectral information combined with a preprocessing algorithm, and then the DNA of the samples was extracted to study the correlation between genetic and relative spectral distances, the model identification performance, and finally to compare the similarities and differences between the results of genetic clustering and relative spectral clustering. Results: The results were as follows: the average accuracy of the models was 93.6% (inbred lines) and 93.7% (hybrids), respectively; Genetic distance and correlation spectral distance exhibited positive correlation significantly (inbred lines: r=0.177, p<0.05; hybrids: r=0.238, p<0.05), likewise genetic distance and model accuracy also showed positive correlation (inbred lines: r=0.611, p<0.01; hybrids: r=0.6158, p<0.01); Genetic clustering and spectral clustering results were essentially uniform for 94.3% (inbred lines) and 93.9% (hybrids), respectively. Discussion: This study reveals the relationship between the genetic and relative spectral distances of seeds and the accuracy of the model, which provides theoretical basis for the establishment of the standardized system for detecting the authenticity of seeds by NIR spectroscopic techniques.

Targeting PRMT9-mediated arginine methylation suppresses cancer stem cell maintenance and elicits cGAS-mediated anticancer immunity.

Current anticancer therapies cannot eliminate all cancer cells, which hijack normal arginine methylation as a means to promote their maintenance via unknown mechanisms. Here we show that targeting protein arginine N-methyltransferase 9 (PRMT9), whose activities are elevated in blasts and leukemia stem cells (LSCs) from patients with acute myeloid leukemia (AML), eliminates disease via cancer-intrinsic mechanisms and cancer-extrinsic type I interferon (IFN)-associated immunity. PRMT9 ablation in AML cells decreased the arginine methylation of regulators of RNA translation and the DNA damage response, suppressing cell survival. Notably, PRMT9 inhibition promoted DNA damage and activated cyclic GMP-AMP synthase, which underlies the type I IFN response. Genetically activating cyclic GMP-AMP synthase in AML cells blocked leukemogenesis. We also report synergy of a PRMT9 inhibitor with anti-programmed cell death protein 1 in eradicating AML. Overall, we conclude that PRMT9 functions in survival and immune evasion of both LSCs and non-LSCs; targeting PRMT9 may represent a potential anticancer strategy.

KDM3B Single-Nucleotide Polymorphisms Impact Radiation Therapy Toxicity Through Circular RNA-Mediated KDM3B Expression and Inflammatory Responses.

PURPOSE: Genome-wide association studies have identified single-nucleotide polymorphisms (SNPs) associated with radiation therapy (RT) toxicities in patients with prostate cancer. SNP rs17599026 in intron 21 of KDM3B is significantly associated with the development of late urinary toxicity, specifically in the increase in urinary frequency 2 years after RT compared with pretreatment conditions. The present study aimed to provide mechanistic insights for this association. METHODS AND MATERIALS: Using human tissues and cell lines, we examined the protein expression of KDM3B and molecular mechanisms underlying the SNP modulation by variants of KDM3B SNP alleles. In animals with normal and heterozygous expressions of Kdm3b, we examined the relationship between Kdm3b expression and radiation toxicity. RESULTS: KDM3B rs17599026 lies in a motif important for circular RNA expression that is responsible for sponging miRNAs to regulate KDM3B expression. Using a murine model with heterozygous deletion of the Kdm3b gene, we found that lower Kdm3b expression is associated with altered pattern of urination after bladder irradiation, which is related to differential degrees of tissue inflammation as measured by analyses of gene expression, lymphocyte infiltration, and noninvasive ultrasound imaging. CONCLUSIONS: KDM3B SNPs can impact its expression through regulating noncoding RNA expression. Differential KDM3B expression underlies radiation toxicity through tissue inflammation at the molecular and physiological level. Our study outcome offers a foundation for mechanism-based mitigation for radiation toxicity for prostate cancer survivors.

POLD1 DEDD Motif Mutation Confers Hypermutation in Endometrial Cancer and Durable Response to Pembrolizumab.

BACKGROUND: Mutations in the DNA polymerase delta 1 (POLD1) exonuclease domain cause DNA proofreading defects, hypermutation, hereditary colorectal and endometrial cancer, and are predictive of immunotherapy response. Exonuclease activity is carried out by two magnesium cations, bound to four highly conserved, negatively charged amino acids (AA) consisting of aspartic acid at amino acid position 316 (p.D316), glutamic acid at position 318 (p.E318), p.D402, and p.D515 (termed DEDD motif). Germline polymorphisms resulting in charge-discordant AA substitutions in the DEDD motif are classified as variants of uncertain significance (VUSs) by laboratories and thus would be considered clinically inactionable. We hypothesize this mutation class is clinically pathogenic. METHODS: A review of clinical presentation was performed in our index patient with a POLD1(p.D402N) heterozygous proband with endometrial cancer. Implications of this mutation class were evaluated by a Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA)-guided systematic review, in silico analysis with orthogonal biochemical confirmation, and whole-exome and RNA sequencing analysis of the patient's tumor and engineered cell lines. RESULTS: Our systematic review favored a Mendelian disease mutation class associated with endometrial and colorectal cancers. In silico analysis predicted defective protein function, confirmed by biochemical assay demonstrating loss of nuclease activity. A POLD1-specific mutational signature was found in both the patient's tumor and POLD1(p.D402N) overexpressing cell. Furthermore, paired whole-exome/transcriptome analysis of endometrial tumor demonstrated hypermutation and T cell-inflamed gene expression profile (GEP), which are joint predictive biomarkers for pembrolizumab. Our patient showed a deep, durable response to immune checkpoint inhibitor (ICI). CONCLUSION: Charge-discordant AA substitution in the DEDD motif of POLD1 is detrimental to DNA proofreading and should be reclassified as likely pathogenic and possibly predictive of ICI sensitivity.

FANCD2 and RAD51 recombinase directly inhibit DNA2 nuclease at stalled replication forks and FANCD2 acts as a novel RAD51 mediator in strand exchange to promote genome stability.

FANCD2 protein, a key coordinator and effector of the interstrand crosslink repair pathway, is also required to prevent excessive nascent strand degradation at hydroxyurea-induced stalled forks. The RAD51 recombinase has also been implicated in regulation of resection at stalled replication forks. The mechanistic contributions of these proteins to fork protection are not well understood. Here, we used purified FANCD2 and RAD51 to study how each protein regulates DNA resection at stalled forks. We characterized three mechanisms of FANCD2-mediated fork protection: (1) The N-terminal domain of FANCD2 inhibits the essential DNA2 nuclease activity by directly binding to DNA2 accounting for over-resection in FANCD2 defective cells. (2) Independent of dimerization with FANCI, FANCD2 itself stabilizes RAD51 filaments to inhibit multiple nucleases, including DNA2, MRE11 and EXO1. (3) Unexpectedly, we uncovered a new FANCD2 function: by stabilizing RAD51 filaments, FANCD2 acts to stimulate the strand exchange activity of RAD51. Our work biochemically explains non-canonical mechanisms by which FANCD2 and RAD51 protect stalled forks. We propose a model in which the strand exchange activity of FANCD2 provides a simple molecular explanation for genetic interactions between FANCD2 and BRCA2 in the FA/BRCA fork protection pathway.

Okazaki fragment maturation: DNA flap dynamics for cell proliferation and survival.

Unsuccessful processing of Okazaki fragments leads to the accumulation of DNA breaks which are associated with many human diseases including cancer and neurodegenerative disorders. Recently, Okazaki fragment maturation (OFM) has received renewed attention regarding how unprocessed Okazaki fragments are sensed and repaired, and how inappropriate OFM impacts on genome stability and cell viability, especially in cancer cells. We provide an overview of the highly efficient and faithful canonical OFM pathways and their regulation of genomic integrity and cell survival. We also discuss how cells induce alternative error-prone OFM processes to promote cell survival in response to environmental stresses. Such stress-induced OFM processes may be important mechanisms driving mutagenesis, cellular evolution, and resistance to radio/chemotherapy and targeted therapeutics in human cancers.

RNA N6 -methyladenosine reader YTHDC1 is essential for TGF-beta-mediated metastasis of triple negative breast cancer.

RNA N6 -methyladenosine (m6A) modification and its regulators fine tune gene expression and contribute to tumorigenesis. This study aims to uncover the essential role and the underlying molecular mechanism(s) of the m6A reader YTHDC1 in promoting triple negative breast cancer (TNBC) metastasis. METHODS: In vitro and in vivo models were employed to determine the pathological function of YTHDC1 in TNBC metastasis. To identify bona fide YTHDC1 target RNAs, we conducted RNA-seq, m6A-seq, and RIP-seq, followed by integrative data analysis and validation assays. RESULTS: By analyzing The Cancer Genome Atlas (TCGA) dataset, we found that elevated expression of YTHDC1 is positively correlated with poor prognosis in breast cancer patients. Using a mammary fat pad mouse model of TNBC, YTHDC1 significantly promoted lung metastasis of TNBC cells. Through multiple transcriptome-wide sequencing and integrative data analysis, we revealed dysregulation of metastasis-related pathways following YTHDC1 depletion and identified SMAD3 as a bona fide YTHDC1 target RNA. Depletion of YTHDC1 caused nuclear retention of SMAD3 mRNA, leading to lower SMAD3 protein levels. Loss of YTHDC1 led to impaired TGF-ß-induced gene expression, leading to inhibition of epithelial-mesenchymal transition (EMT) and suppressed TNBC cell migration and invasion. SMAD3 overexpression was able to restore the response to TGF-ß in YTHDC1 depleted TNBC cells. Furthermore, we demonstrated that the oncogenic role of YTHDC1 is mediated through its recognition of m6A as m6A-binding defective mutants of YTHDC1 were unable to rescue the impaired cell migration and invasion of YTHDC1 knockout TNBC cells. CONCLUSIONS: We show that YTHDC1 plays a critical oncogenic role in TNBC metastasis through promoting the nuclear export and expression of SMAD3 to augment the TGF-ß signaling cascade. Overall, our study demonstrates that YTHDC1 is vital for TNBC progression by enhancing TNBC cell survival and TGF-ß-mediated EMT via SMAD3 to enable the formation of distant metastasis and highlights the therapeutic potential of targeting the YTHDC1/m6A/SMAD3 axis for TNBC treatment.

Structure-specific nucleases: role in Okazaki fragment maturation.

Proper function of structure-specific nucleases is key for faithful Okazaki fragment maturation (OFM) process completion. Deregulation of such nucleases leads to aberrant OFM and causes a spectrum of mutations, some of which may confer survival outcomes under specific stresses and serve as attractive targets for therapeutic intervention in human cancers.

Error-prone, stress-induced 3' flap-based Okazaki fragment maturation supports cell survival.

How cells with DNA replication defects acquire mutations that allow them to escape apoptosis under environmental stress is a long-standing question. Here, we report that an error-prone Okazaki fragment maturation (OFM) pathway is activated at restrictive temperatures in rad27Δ yeast cells. Restrictive temperature stress activated Dun1, facilitating transformation of unprocessed 5' flaps into 3' flaps, which were removed by 3' nucleases, including DNA polymerase δ (Polδ). However, at certain regions, 3' flaps formed secondary structures that facilitated 3' end extension rather than degradation, producing alternative duplications with short spacer sequences, such as pol3 internal tandem duplications. Consequently, little 5' flap was formed, suppressing rad27Δ-induced lethality at restrictive temperatures. We define a stress-induced, error-prone OFM pathway that generates mutations that counteract replication defects and drive cellular evolution and survival.

Role of 53BP1 in end protection and DNA synthesis at DNA breaks.

Double-strand break (DSB) repair choice is greatly influenced by the initial processing of DNA ends. 53BP1 limits the formation of recombinogenic single-strand DNA (ssDNA) in BRCA1-deficient cells, leading to defects in homologous recombination (HR). However, the exact mechanisms by which 53BP1 inhibits DSB resection remain unclear. Previous studies have identified two potential pathways: protection against DNA2/EXO1 exonucleases presumably through the Shieldin (SHLD) complex binding to ssDNA, and localized DNA synthesis through the CTC1-STN1-TEN1 (CST) and DNA polymerase α (Polα) to counteract resection. Using a combinatorial approach of END-seq, SAR-seq, and RPA ChIP-seq, we directly assessed the extent of resection, DNA synthesis, and ssDNA, respectively, at restriction enzyme-induced DSBs. We show that, in the presence of 53BP1, Polα-dependent DNA synthesis reduces the fraction of resected DSBs and the resection lengths in G0/G1, supporting a previous model that fill-in synthesis can limit the extent of resection. However, in the absence of 53BP1, Polα activity is sustained on ssDNA yet does not substantially counter resection. In contrast, EXO1 nuclease activity is essential for hyperresection in the absence of 53BP1. Thus, Polα-mediated fill-in partially limits resection in the presence of 53BP1 but cannot counter extensive hyperresection due to the loss of 53BP1 exonuclease blockade. These data provide the first nucleotide mapping of DNA synthesis at resected DSBs and provide insight into the relationship between fill-in polymerases and resection exonucleases.

Role of the microbiota in ileitis of a mouse model of inflammatory bowel disease-Glutathione peroxide isoenzymes 1 and 2-double knockout mice on a C57BL background.

C57Bl6 (B6) mice devoid of glutathione peroxidases 1 and 2 (Gpx1/2-DKO) develop ileitis after weaning. We previously showed germ-free Gpx1/2-DKO mice of mixed B6.129 background did not develop ileocolitis. Here, we examine the composition of the ileitis provoking microbiota in B6 Gpx1/2-DKO mice. DNA was isolated from the ileum fecal stream and subjected to high-throughput sequencing of the V3 and V4 regions of the 16S rRNA gene to determine the abundance of operational taxonomic units (OTUs). We analyzed the role of bacteria by comparing the microbiomes of the DKO and pathology-free non-DKO mice. Mice were treated with metronidazole, streptomycin, and vancomycin to alter pathology and correlate the OTU abundances with pathology levels. Principal component analysis based on Jaccard distance of abundance showed 3 distinct outcomes relative to the source Gpx1/2-DKO microbiome. Association analyses of pathology and abundance of OTUs served to rule out 7-11 of 24 OTUs for involvement in the ileitis. Collections of OTUs were identified that appeared to be linked to ileitis in this animal model and would be classified as commensals. In Gpx1/2-DKO mice, host oxidant generation from NOX1 and DUOX2 in response to commensals may compromise the ileum epithelial barrier, a role generally ascribed to oxidants generated from mitochondria, NOX2 and endoplasmic reticulum stress in response to presumptive pathogens in IBD. Elevated oxidant levels may contribute to epithelial cell shedding, which is strongly associated with progress toward inflammation in Gpx1/2-DKO mice and predictive of relapse in IBD by allowing leakage of microbial components into the submucosa.

Arginine methylation of APE1 promotes its mitochondrial translocation to protect cells from oxidative damage.

Apurinic/apyrimidinic endonuclease 1 (APE1) is an essential multifunctional protein in mammals that plays critical roles in DNA repair and redox signaling within the cell. Impaired APE1 function or dysregulation is associated with disease susceptibility and poor cancer prognosis. Orchestrated regulatory mechanisms are crucial to ensure its function in a specific subcellular location at specific time. Here, we report arginine methylation as a post-translational modification (PTM) that regulates APE1 translocation to mitochondria in HeLa and HEK-293 cells. Protein arginine methyl-transferase 1 (PRMT1) was shown to methylate APE1 in vitro. Site-directed mutagenesis identified R301 as the major methylation site. We confirmed that APE1 is methylated in cells and that the R301K mutation significantly reduces its methylation. Baseline mitochondrial APE1 levels were low under standard culture conditions, but they could be induced by oxidative agents. Methylation-deficient APE1 showed reduced mitochondrial translocation. Methylation affected the interaction of APE1 with Tom20, translocase of the outer mitochondrial membrane. Methylation-deficient APE1 resulted in increased mitochondrial DNA damage and increased cytochrome c release after stimuli. These data suggest that methylation of APE1 promotes its mitochondrial translocation and protects cells from oxidative damage. This work describes a novel PTM regulation model of APE1 subcellular distribution through arginine methylation.

Ubiquitinated-PCNA protects replication forks from DNA2-mediated degradation by regulating Okazaki fragment maturation and chromatin assembly.

Upon genotoxic stress, PCNA ubiquitination allows for replication of damaged DNA by recruiting lesion-bypass DNA polymerases. However, PCNA is also ubiquitinated during normal S-phase progression. By employing 293T and RPE1 cells deficient in PCNA ubiquitination, generated through CRISPR/Cas9 gene editing, here, we show that this modification promotes cellular proliferation and suppression of genomic instability under normal growth conditions. Loss of PCNA-ubiquitination results in DNA2-dependent but MRE11-independent nucleolytic degradation of nascent DNA at stalled replication forks. This degradation is linked to defective gap-filling in the wake of the replication fork and incomplete Okazaki fragment maturation, which interferes with efficient PCNA unloading by ATAD5 and subsequent nucleosome deposition by CAF-1. Moreover, concomitant loss of PCNA-ubiquitination and the BRCA pathway results in increased nascent DNA degradation and PARP inhibitor sensitivity. In conclusion, we show that by ensuring efficient Okazaki fragment maturation, PCNA-ubiquitination protects fork integrity and promotes the resistance of BRCA-deficient cells to PARP-inhibitors.

Multiple roles of DNA2 nuclease/helicase in DNA metabolism, genome stability and human diseases.

DNA2 nuclease/helicase is a structure-specific nuclease, 5'-to-3' helicase, and DNA-dependent ATPase. It is involved in multiple DNA metabolic pathways, including Okazaki fragment maturation, replication of 'difficult-to-replicate' DNA regions, end resection, stalled replication fork processing, and mitochondrial genome maintenance. The participation of DNA2 in these different pathways is regulated by its interactions with distinct groups of DNA replication and repair proteins and by post-translational modifications. These regulatory mechanisms induce its recruitment to specific DNA replication or repair complexes, such as DNA replication and end resection machinery, and stimulate its efficient cleavage of various structures, for example, to remove RNA primers or to produce 3' overhangs at telomeres or double-strand breaks. Through these versatile activities at replication forks and DNA damage sites, DNA2 functions as both a tumor suppressor and promoter. In normal cells, it suppresses tumorigenesis by maintaining the genomic integrity. Thus, DNA2 mutations or functional deficiency may lead to cancer initiation. However, DNA2 may also function as a tumor promoter, supporting cancer cell survival by counteracting replication stress. Therefore, it may serve as an ideal target to sensitize advanced DNA2-overexpressing cancers to current chemo- and radiotherapy regimens.

Functional deficiency of DNA repair gene EXO5 results in androgen-induced genomic instability and prostate tumorigenesis.

Germline mutations of DNA double-strand break (DSB) response and repair genes that drive tumorigenesis could be a major cause of prostate cancer (PCa) heritability. In this study, we demonstrated the role of novel exonuclease 5 (EXO5) gene in androgen-induced double strand breaks repair via homology-directed repair pathway and prostate tumorigenesis. Using whole-exome sequencing of samples from 20 PCa families, with three or more siblings diagnosed with metastatic PCa, we identified mutations in 31 genes involved in DSB response and repair. Among them, the L151P mutation in the exonuclease 5 (EXO5) gene was present in all affected siblings in three PCa families. We found two other EXO5 SNPs significantly associated with risk of PCa in cases-controls study from databases of genotype and phenotype (dbGaP), which are in linkage disequilibrium (D' = 1) with Exo5 L151P found in PCa family. The L151 residue is conserved across diverse species and its mutation is deleterious for protein functions, as demonstrated by our bioinformatics analyses. The L151P mutation impairs the DNA repair function of EXO5 due to loss of nuclease activity, as well as failure of nuclear localization. CRISPR elimination of EXO5 in a PCa cell line impaired homology-directed recombination repair (HDR) and caused androgen-induced genomic instability, as indicated by frequent occurrence of the oncogenic fusion transcript TMPRSS2-ERG. Genetic and functional validation of the EXO5 mutations indicated that EXO5 is a risk gene for prostate tumorigenesis, likely due to its functions in HDR.

Dexamethasone and Tofacitinib suppress NADPH oxidase expression and alleviate very-early-onset ileocolitis in mice deficient in GSH peroxidase 1 and 2.

C57BL6/J (B6) mice lacking Se-dependent GSH peroxidase 1 and 2 (GPx1/2-DKO) develop mild to moderate ileocolitis around weaning. These DKO mice have a disease resembling human very-early-onset inflammatory bowel disease (VEOIBD), which is associated with mutations in NADPH oxidase genes. Drugs including dexamethasone (Dex), Tofacitinib (Tofa; a Janus kinase/JAK inhibitor) and anti-TNF antibody are effective to treat adult, but not pediatric IBD. AIMS: To test the efficacy of hydrophobic Dex and hydrophilic Dex phosphate (Dex phos), Tofa, anti-Tnf Ab, Noxa1ds-TAT and gp91ds-TAT peptides (inhibiting NOX1 and NOX2 assembly respectively), antioxidant MJ33 and ML090, and pifithrin-α (p53 inhibitor) on alleviation of gut inflammation in DKO weanlings. MAIN METHODS: All treatments began on 22-day-old GPx1/2-DKO mice. The mouse intestine pathology was compared between the drug- and vehicle-treated groups after six or thirteen days of treatment. KEY FINDINGS: Among all drugs tested, Dex, Dex phos and Tofa were the strongest to suppress ileocolitis in the DKO weanlings. Dex, Dex phos and Tofa inhibited crypt apoptosis and increased crypt density. Dex or Dex phos alone also inhibited cell proliferation, exfoliation and crypt abscess in the ileum. Dex, but not Tofa, retarded mouse growth. Both Dex and Tofa inhibited ileum Nox1, Nox4 and Duox2, but not Nox2 gene expression. Noxa1ds-TAT and gp91ds-TAT peptides as well as MJ33 had subtle effect on suppressing pathology, while others had negligible effect. SIGNIFICANCE: These findings suggest that NADPH oxidases can be novel drug targets for pediatric IBD therapy, and Tofa may be considered for treating VEOIBD.