A computer-aided drug design approach to discover tumour suppressor p53 protein activators for colorectal cancer therapy.

Colorectal cancer (CRC) is the third most detected cancer and the second foremost cause of cancer deaths in the world. Intervention targeting p53 provides potential therapeutic strategies, but thus far no p53-based therapy has been successfully translated into clinical cancer treatment. Here we developed a Quantitative Structure-Activity Relationships (QSAR) classification models using empirical molecular descriptors and fingerprints to predict the activity against the p53 protein, using the potency value with the active or inactive label, were developed. These models were built using in total 10,505 molecules that were extracted from the ChEMBL, ZINC and Reaxys® databases, and recent literature. Three machine learning (ML) techniques e.g., Random Forest, Support Vector Machine, Convolutional Neural Network were explored to build models for p53 inhibitor prediction. The performances of the models were successfully evaluated by internal and external validation. Moreover, based on the best in silico p53 model, a virtual screening campaign was carried out using 1443 FDA-approved drugs that were extracted from the ZINC database. A list of virtual screening hits was assented on base of some limits established in this approach, such as: (1) probability of being active against p53; (2) applicability domain; (3) prediction of the affinity between the p53, and ligands, through molecular docking. The most promising according to the limits established above was dihydroergocristine. This compound revealed cytotoxic activity against a p53-expressing CRC cell line with an IC50 of 56.8 µM. This study demonstrated that the computer-aided drug design approach can be used to identify previously unknown molecules for targeting p53 protein with anti-cancer activity and thus pave the way for the study of a therapeutic solution for CRC.

The impact of p53 on aristolochic acid I-induced nephrotoxicity and DNA damage in vivo and in vitro.

Exposure to aristolochic acid (AA) is associated with human nephropathy and urothelial cancer. The tumour suppressor TP53 is a critical gene in carcinogenesis and frequently mutated in AA-induced urothelial tumours. We investigated the impact of p53 on AAI-induced nephrotoxicity and DNA damage in vivo by treating Trp53(+/+), Trp53(+/-) and Trp53(-/-) mice with 3.5 mg/kg body weight (bw) AAI daily for 2 or 6 days. Renal histopathology showed a gradient of intensity in proximal tubular injury from Trp53(+/+) to Trp53(-/-) mice, especially after 6 days. The observed renal injury was supported by nuclear magnetic resonance (NMR)-based metabonomic measurements, where a consistent Trp53 genotype-dependent trend was observed for urinary metabolites that indicate aminoaciduria (i.e. alanine), lactic aciduria (i.e. lactate) and glycosuria (i.e. glucose). However, Trp53 genotype had no impact on AAI-DNA adduct levels, as measured by 32P-postlabelling, in either target (kidney and bladder) or non-target (liver) tissues, indicating that the underlying mechanisms of p53-related AAI-induced nephrotoxicity cannot be explained by differences in AAI genotoxicity. Performing gas chromatography-mass spectrometry (GC-MS) on kidney tissues showed metabolic pathways affected by AAI treatment, but again Trp53 status did not clearly impact on such metabolic profiles. We also cultured primary mouse embryonic fibroblasts (MEFs) derived from Trp53(+/+), Trp53(+/-) and Trp53(-/-) mice and exposed them to AAI in vitro (50 µM for up to 48 h). We found that Trp53 genotype impacted on the expression of NAD(P)H:quinone oxidoreductase (Nqo1), a key enzyme involved in AAI bioactivation. Nqo1 induction was highest in Trp53(+/+) MEFs and lowest in Trp53(-/-) MEFs; and it correlated with AAI-DNA adduct formation, with lowest adduct levels being observed in AAI-exposed Trp53(-/-) MEFs. Overall, our results clearly demonstrate that p53 status impacts on AAI-induced renal injury, but the underlying mechanism(s) involved remain to be further explored. Despite the impact of p53 on AAI bioactivation and DNA damage in vitro, such effects were not observed in vivo.

The Impact of p53 on Aristolochic Acid I-Induced Gene Expression In Vivo.

Exposure to aristolochic acid (AA) is linked to kidney disease and urothelial cancer in humans. The major carcinogenic component of the AA plant extract is aristolochic acid I (AAI). The tumour suppressor p53 is frequently mutated in AA-induced tumours. We previously showed that p53 protects from AAI-induced renal proximal tubular injury, but the underlying mechanism(s) involved remain to be further explored. In the present study, we investigated the impact of p53 on AAI-induced gene expression by treating Trp53(+/+), Trp53(+/-), and Trp53(-/-) mice with 3.5 mg/kg body weight (bw) AAI daily for six days. The Clariom™ S Assay microarray was used to elucidate gene expression profiles in mouse kidneys after AAI treatment. Analyses in Qlucore Omics Explorer showed that gene expression in AAI-exposed kidneys is treatment-dependent. However, gene expression profiles did not segregate in a clear-cut manner according to Trp53 genotype, hence further investigations were performed by pathway analysis with MetaCore™. Several pathways were significantly altered to varying degrees for AAI-exposed kidneys. Apoptotic pathways were modulated in Trp53(+/+) kidneys; whereas oncogenic and pro-survival pathways were significantly altered for Trp53(+/-) and Trp53(-/-) kidneys, respectively. Alterations of biological processes by AAI in mouse kidneys could explain the mechanisms by which p53 protects from or p53 loss drives AAI-induced renal injury in vivo.

Carcinogenic polycyclic aromatic hydrocarbons induce CYP1A1 in human cells via a p53-dependent mechanism.

The tumour suppressor gene TP53 is mutated in more than 50 % of human tumours, making it one of the most important cancer genes. We have investigated the role of TP53 in cytochrome P450 (CYP)-mediated metabolic activation of three polycyclic aromatic hydrocarbons (PAHs) in a panel of isogenic colorectal HCT116 cells with differing TP53 status. Cells that were TP53(+/+), TP53(+/-), TP53(-/-), TP53(R248W/+) or TP53(R248W/-) were treated with benzo[a]pyrene (BaP), dibenz[a,h]anthracene and dibenzo[a,l]pyrene, and the formation of DNA adducts was measured by (32)P-postlabelling analysis. Each PAH formed significantly higher DNA adduct levels in TP53(+/+) cells than in the other cell lines. There were also significantly lower levels of PAH metabolites in the culture media of these other cell lines. Bypass of the need for metabolic activation by treating cells with the corresponding reactive PAH-diol-epoxide metabolites resulted in similar adduct levels in all cell lines, which confirms that the influence of p53 is on the metabolism of the parent PAHs. Western blotting showed that CYP1A1 protein expression was induced to much greater extent in TP53(+/+) cells than in the other cell lines. CYP1A1 is inducible via the aryl hydrocarbon receptor (AHR), but we did not find that expression of AHR was dependent on p53; rather, we found that BaP-induced CYP1A1 expression was regulated through p53 binding to a p53 response element in the CYP1A1 promoter region, thereby enhancing its transcription. This study demonstrates a new pathway for CYP1A1 induction by environmental PAHs and reveals an emerging role for p53 in xenobiotic metabolism.

The impact of p53 on DNA damage and metabolic activation of the environmental carcinogen benzo[a]pyrene: effects in Trp53(+/+), Trp53(+/-) and Trp53(-/-) mice.

The tumour suppressor p53 is one of the most important cancer genes. Previous findings have shown that p53 expression can influence DNA adduct formation of the environmental carcinogen benzo[a]pyrene (BaP) in human cells, indicating a role for p53 in the cytochrome P450 (CYP) 1A1-mediated biotransformation of BaP in vitro. We investigated the potential role of p53 in xenobiotic metabolism in vivo by treating Trp53(+/+), Trp53(+/-) and Trp53(-/-) mice with BaP. BaP-DNA adduct levels, as measured by (32)P-postlabelling analysis, were significantly higher in liver and kidney of Trp53(-/-) mice than of Trp53(+/+) mice. Complementarily, significantly higher amounts of BaP metabolites were also formed ex vivo in hepatic microsomes from BaP-pretreated Trp53(-/-) mice. Bypass of the need for metabolic activation by treating mice with BaP-7,8-dihydrodiol-9,10-epoxide resulted in similar adduct levels in liver and kidney in all mouse lines, confirming that the influence of p53 is on the biotransformation of the parent compound. Higher BaP-DNA adduct levels in the livers of Trp53(-/-) mice correlated with higher CYP1A protein levels and increased CYP1A enzyme activity in these animals. Our study demonstrates a role for p53 in the metabolism of BaP in vivo, confirming previous in vitro results on a novel role for p53 in CYP1A1-mediated BaP metabolism. However, our results also suggest that the mechanisms involved in the altered expression and activity of the CYP1A1 enzyme by p53 in vitro and in vivo are different.

Characterising Mutational Spectra of Carcinogens in the Tumour Suppressor Gene TP53 Using Human TP53 Knock-in (Hupki) Mouse Embryo Fibroblasts.

DNA in dividing cells is prone to mutagenesis, with mutations making key contributions to human disease including cancer. The tumour suppressor gene TP53 is the most frequently mutated gene in human tumours. Here, we present a robust protocol for studying TP53 mutagenesis utilising human TP53 knock-in (Hupki) mouse embryonic fibroblasts (HUFs). In the HUF immortalisation assay (HIMA), primary HUFs are treated with known or suspected carcinogens at 3% oxygen and then transferred to 20% atmospheric oxygen to induce senescence. Cells containing mutations (e.g., in TP53) that allow bypassing of senescence eventually emerge as immortalised clonal cell lines after 2-3 months of serial passaging. As not all immortalised HUF cells contain TP53 mutations, we developed a Nutlin-3a counter-screen to select for TP53-mutated clones prior to sequencing. TP53 mutation spectra generated can be compared with those of human tumours recorded in the International Agency for Research on Cancer TP53 mutation database. Environmental mutagens that have demonstrated and validated the utility of the HIMA include ultraviolet radiation, aristolochic acid, and benzo[a]pyrene. The TP53 mutation patterns induced by these mutagens in the HIMA corresponded to those found in human tumours from patients exposed to these mutagens. The approach presented helps to deepen our understanding of human cancer aetiology.

PCB-126 effects on aryl hydrocarbon receptor, ubiquitin and p53 expression levels in a fish product (Sparus aurata L.).

The aim of the present study is to determine if Ahr ligands as PCB-126, a dioxin-like, might contribute to inhibition of the tumour suppressor p53 by promoting its degradation through proteasome-ubiquitin system (UPS). The findings show, in the presence of PCB-126, a significant increase in p53 immunoreactivity in fish compared to the control. Subsequently, there is a decrease in p53 immunoreactivity at 24 h which is maintained even at 72 h. There is also a slight decrease in ubiquitin immunoreactivity to 12 h compared to the control and a marked decrease to 24 and 72 h. It's very important to underline as in this study we demonstrate a marked decrease in ubiquitin and p53 immunoreactivity at 24 and 72 h. Our result emphasise the need to deeply the role of this receptor in UPS regulation as potential therapeutic target for cancer treatment.

The impact of chemotherapeutic drugs on the CYP1A1-catalysed metabolism of the environmental carcinogen benzo[a]pyrene: Effects in human colorectal HCT116 TP53(+/+), TP53(+/-) and TP53(-/-) cells.

Polycyclic aromatic hydrocarbons such as benzo[a]pyrene (BaP) can induce cytochrome P450 1A1 (CYP1A1) via a p53-dependent mechanism. The effect of different p53-activating chemotherapeutic drugs on CYP1A1 expression, and the resultant effect on BaP metabolism, was investigated in a panel of isogenic human colorectal HCT116 cells with differing TP53 status. Cells that were TP53(+/+), TP53(+/-) or TP53(-/-) were treated for up to 48 h with 60 µM cisplatin, 50 µM etoposide or 5 µM ellipticine, each of which caused high p53 induction at moderate cytotoxicity (60-80% cell viability). We found that etoposide and ellipticine induced CYP1A1 in TP53(+/+) cells but not in TP53(-/-) cells, demonstrating that the mechanism of CYP1A1 induction is p53-dependent; cisplatin had no such effect. Co-incubation experiments with the drugs and 2.5 µM BaP showed that: (i) etoposide increased CYP1A1 expression in TP53(+/+) cells, and to a lesser extent in TP53(-/-) cells, compared to cells treated with BaP alone; (ii) ellipticine decreased CYP1A1 expression in TP53(+/+) cells in BaP co-incubations; and (iii) cisplatin did not affect BaP-mediated CYP1A1 expression. Further, whereas cisplatin and etoposide had virtually no influence on CYP1A1-catalysed BaP metabolism, ellipticine treatment strongly inhibited BaP bioactivation. Our results indicate that the underlying mechanisms whereby etoposide and ellipticine regulate CYP1A1 expression must be different and may not be linked to p53 activation alone. These results could be relevant for smokers, who are exposed to increased levels of BaP, when prescribing chemotherapeutic drugs. Beside gene-environment interactions, more considerations should be given to potential drug-environment interactions during chemotherapy.

Quantitative analysis of gene expression changes in response to genotoxic compounds.

Techniques that quantify molecular endpoints sufficiently sensitive to identify and classify potentially toxic compounds have wide potential for high-throughput in vitro screening. Expression of three genes, RAD51C, TP53 and cystatin A (CSTA), in HEPG2 cells was measured by Q-PCR amplification. In parallel, we developed alternative assays for the same 3 gene signature based on an acridinium-ester chemiluminescent reporter molecule. HEPG2 cells were challenged with eighteen different compounds (n=18) chosen to represent compounds that are genotoxic (n=8), non-genotoxic non-carcinogenic (n=2) or have a less well defined mechanism of action with respect to genotoxicity (n=8). At least one of the three genes displayed dysregulated expression in the majority of compounds tested by Q-PCR and ten compounds changed the CSTA expression significantly. Acridinium-ester labelled probes for the three genes were synthesised and tested. Analytical sensitivity was characterised and suggested a limit of detection generally better than 0.1fmol but often 10-50 attomol. A linear amplification step was optimised and this quantitative method detected statistically significant increases in RAD51C and CSTA expression in agreement with the Q-PCR results, demonstrating the potential of this technology. The broad agreement of the amplified chemiluminescent method and Q-PCR in measuring gene expression suggests wider potential application for this chemiluminescent technology.