9
views
0
recommends
+1 Recommend
0 collections
    0
    shares
      • Record: found
      • Abstract: found
      • Article: found
      Is Open Access

      Upregulation of p53 through induction of MDM2 degradation: improved potency through the introduction of an alkylketone sidechain on the anthraquinone core

      research-article
      , , , ,
      Journal of Enzyme Inhibition and Medicinal Chemistry
      Taylor & Francis
      p53, MDM2, anthraquinone, cancer, structure–activity relationship (SAR)

      Read this article at

      Bookmark
          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

          Abstract

          Overexpression of ubiquitin ligase MDM2 causes depletion of the p53 tumour-suppressor and thus leads to cancer progression. In recent years, anthraquinone analogs have received significant attention due to their ability to downregulate MDM2, thereby promoting p53-induced apoptosis. Previously, we have developed potent anthraquinone compounds having the ability to upregulate p53 via inhibition of MDM2 in both cell culture and animal models of acute lymphocytic leukaemia. Earlier work was focussed on mechanistic work, pharmacological validation of this class of compounds in animal models, and mapping out structural space that allows for further modification and optimisation. Herein, we describe our work in optimising the substituents on the two phenol hydroxyl groups. It was found that the introduction of an alkylketone moiety led to a potent series of analogs with BW-AQ-350 being the most potent compound yet (IC 50 = 0.19 ± 0.01 µM) which exerts cytotoxicity by inducing MDM2 degradation and p53 upregulation.

          Related collections

          Most cited references63

          • Record: found
          • Abstract: found
          • Article: not found

          Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings.

          Experimental and computational approaches to estimate solubility and permeability in discovery and development settings are described. In the discovery setting 'the rule of 5' predicts that poor absorption or permeation is more likely when there are more than 5 H-bond donors, 10 H-bond acceptors, the molecular weight (MWT) is greater than 500 and the calculated Log P (CLogP) is greater than 5 (or MlogP > 4.15). Computational methodology for the rule-based Moriguchi Log P (MLogP) calculation is described. Turbidimetric solubility measurement is described and applied to known drugs. High throughput screening (HTS) leads tend to have higher MWT and Log P and lower turbidimetric solubility than leads in the pre-HTS era. In the development setting, solubility calculations focus on exact value prediction and are difficult because of polymorphism. Recent work on linear free energy relationships and Log P approaches are critically reviewed. Useful predictions are possible in closely related analog series when coupled with experimental thermodynamic solubility measurements.
            Bookmark
            • Record: found
            • Abstract: found
            • Article: not found

            The resurgence of covalent drugs.

            Covalent drugs have proved to be successful therapies for various indications, but largely owing to safety concerns, they are rarely considered when initiating a target-directed drug discovery project. There is a need to reassess this important class of drugs, and to reconcile the discordance between the historic success of covalent drugs and the reluctance of most drug discovery teams to include them in their armamentarium. This review surveys the prevalence and pharmacological advantages of covalent drugs, discusses how potential risks and challenges may be addressed through innovative design, and presents the broad opportunities provided by targeted covalent inhibitors.
              Bookmark
              • Record: found
              • Abstract: found
              • Article: found
              Is Open Access

              Genome dynamics of the human embryonic kidney 293 lineage in response to cell biology manipulations

              The human embryonic kidney (HEK) 293 cell line and its derivatives are used in experiments ranging from signal transduction and protein interaction studies over viral packaging to rapid small-scale protein expression and biopharmaceutical production. The original 293 cells1 2 3 were derived in 1973 from the kidney of an aborted human embryo of unknown parenthood by transformation with sheared Adenovirus 5 DNA. The human embryonic kidney cells at first seemed recalcitrant to transformation. After many attempts, cell growth took off only several months after the isolation of a single transformed clone. This cell line is known as HEK293 or 293 cells (ATCC accession number CRL-1573). A 4-kbp adenoviral genome fragment is known to have integrated in chromosome 19 (ref. 4) and encodes for the E1A/E1B proteins, which interfere with the cell cycle control pathways and counteract apoptosis5 6. Cytogenetic analysis established that the 293 line is pseudotriploid7. Given the broad use of 293 cells for biomedical research and virus/protein production, we decided to perform a comprehensive genomic characterization of the 293 cell line and the most commonly used derived lines (Fig. 1a) to better understand the dynamics of the 293 genome under the procedures commonly used in biotechnological engineering of mammalian cell lines. First among these derived lines, we analysed 293T, which expresses a temperature-sensitive allele of the SV40 T antigen8 9. This enables the amplification of vectors containing the SV40 ori and thus considerably increases the expression levels obtained with transient transfection. SV40 T forms a complex with and inhibits p53, possibly further compromising genome integrity10. The original 293 line was suspension growth-adapted through serial passaging in Joklik’s modified minimal Eagle’s medium11. Full adaptation took about 7 months, and the first passages were so difficult that the few cells that grew through are likely to have been almost clonal (Dr Bruce Stillman, personal communication). The fully adapted cell line is known as 293S and is also analysed here. Subsequently, this line was mutagenized with ethylmethanesulfonate (EMS) and a Ricin toxin-resistant clone was selected out. The line lacked N-acetylglucosaminyltransferase I activity (encoded by the MGAT1 gene) and accordingly predominantly modifies glycoproteins with the Man5GlcNAc2 N-glycan. Then, a stable tetR repressor–expressing clone of this glyco-engineered cell line was derived to enable tetracyclin-inducible protein expression12. This cell line is widely used for the production of homogenously N-glycosylated proteins and will be referred to as 293SG. Apart from these four cell lines in common use, we also analysed the genome of two 293-derived lines used in our laboratory for protein–protein interaction screening (293FTM) and glyco-engineering (293SGGD; details in Supplementary Information). In our study, following genomic studies of other human cell lines13 14 15, we aim to provide a full-genome resource for these cell biology ‘workhorse’ cell lines while developing the necessary tools to make such resources easily available. This enables all researchers using the 293 cell lines to make fully informed analyses of genomic regions of interest to their studies, without expert bioinformatics skills. We also map the genomic changes accumulating after standard laboratory cell culturing (passaging and freezing), providing a way to assess genomic stability of each line. Furthermore, we present a workflow for determining the insertion sites of viral sequences and plasmids based on the genome sequencing data. The extreme chromosome structure diversity/plasticity in the 293 cell line underlies a novel application: selection of 293 clones surviving stringent selective conditions (in our case: ricin toxin), followed by whole-genome analysis of copy number alterations, can effectively pinpoint the genomic region(s) that contain the gene(s) that is required for adaptation to those selective conditions. Results 293 cell lineage genome, karyotype and transcriptome For genome resequencing, we used complete genomics (CG) high-coverage genome sequencing technology16 (Supplementary Methods; data set summary in Supplementary Tables 1 and 2, and sequencing quality overview in Supplementary Fig. 1). 293 cells are of female provenance, as we find no trace of Y-chromosome-derived sequence in our data sets. The mitochondrial sequence belongs to the oldest European haplogroup U5a1 (refs 17, 18). Furthermore, we applied multiplex fluorescence in situ hybridization analysis to our 293 lines (Supplementary Data 1). A wide diversity of karyotypes was found, also within each clone, with some chromosomal alterations relative to the human reference genome present in almost all cells, and others in only a small proportion of cells. Overall, the pseudotriploidy of the 293 lineage was confirmed both by CG sequencing and karyotyping. To further define the 293 cell lineage and to enable the future development of cell line authentication genotyping assays, we analysed which single-nucleotide polymorphisms (SNPs) in protein-coding regions were common to the six sequenced 293 cell lines (Supplementary Data 2) and we manually curated the functional annotation of all novel (that is, not present in dbSNP) 293-defining SNPs (Supplementary Data 2). The genome-wide 2-kb-resolution sequencing coverage depth analysis provides a 2-kb-window copy number that is relative to the genome-averaged copy number in that particular genome. To obtain the absolute copy number, an independent data source is required. For this purpose, we used the Illumina SNP-array-determined genome-averaged ploidy number. The resulting calibrated 2-kb-resolution copy number shows very good consistency with the lower-resolution Illumina SNP-array copy number variant (CNV) results (Spearman rho=0.67–0.80, depending on the cell line; P 300,000 Broad Institute mouse/human genome-wide shRNAs mapped uniquely to the human RefSeq gene collection, visualized these in an IGV annotation track (Fig. 6b) and investigated which of these targets are mutated in our 293 cell lines. Depending on the cell line, this was the case for 9,608–11,534 (~6% of the ones that aligned) of these shRNAs, which may render these nonfunctional in gene silencing. The 293 line was also one of the many cell lines selected for analysis by the ENCODE project37. Several data sets that are highly complementary to ours and deal, for example, with epigenomics are becoming available in this way. We will be updating our web interfaces for the 293 genome with these and other generated data sets on an ongoing basis. Discussion Cell lines are instrumental for our growing understanding of mammalian biology and for biopharmaceutical production. 293 cells are second only to HeLa cells in the frequency of their use in cell biology (a search in PubMed for this cell line and its most popular derivatives yields ~20,000 hits). They are second only to CHO cells for their use in biopharmaceutical production (and take the prime spot for use in small-scale protein production and in viral vector propagation). However, 293 cells were at some point derived from an individual human embryo with a genome different from the reference. Moreover, the establishment of the cell line and its continuous growth in vitro impose selective conditions on the cells, which are often adapted to through mutation. Thus, the human reference genome sequence provides only a partial understanding of the genome of human cell lines. As genome-wide short interfering RNA resources are now available for human cells38 39, and as sequence-specific genome-engineering tools are rapidly becoming standard tools for mammalian cell genetic engineering40 41 42, a sequence and average copy number level knowledge of the entire genomes of the cell lines under study is of great advantage. Furthermore, the cell-line-specific genome sequences reported here will also be beneficial in the interpretation of RNA-seq and proteomics experiments that make use of these cells. 293 cells have been cultivated for decades in different laboratories, which most likely has led to different progressive genome structure alterations. This may underlie the sometimes different conclusions drawn from experimentation with 293 cell lines (and many other cell lines). All cell lines sequenced here are available to the research community. Up to the level of sensitivity afforded by our sequencing approach (single copy plasmid insertions were easily detected), these cell lines have no inadvertent virus insertions, which should help to put at rest some of the concerns towards the use of the 293 cells for biopharmaceutical production. The analytical tools we provide here for integrated plasmids and viral sequences will be very valuable in fully characterizing cell lines used for the production of biopharmaceuticals, both towards the copy number and stability of the inserted plasmids and the validation that such cell lines are free of inadvertent viral sequence contamination. We have shown that comparative sequencing of several 293 lines of the same descent reveal genomic copy number alterations that explain diverse phenotypes of the lineage and its subclones. Extensive further experimentation is now required to validate the role of these CNVs in cellular transformation, suspension growth adaptation and metabolism. We hope that such studies will contribute to the design of new generations of 293 cells that are even better adapted to experimental and pharmaceutical production requirements, and the knowledge gained may be instructive in how to directly engineer other human cell lines. Furthermore, it is clear from our data that the standard practice of generating a stable clone through transfection and selection will result in the isolation of one geno/karyotype present in the parental cell line. Thus, any phenotype of the resulting stable transfectant may be because of the integrated transgene, or may be because of a genomic difference between the new line and its parental line. Consequently, such experiments should be interpreted with great caution and these data argue for the use of efficient transient transfection or propagation of a polyclonal pool of stable transfectants (in which case a more representative population of the parental cells is analysed) in, for example, quantitative signal transduction studies that use 293 cells (as is used in many drug screening and ‘omics’ experiments). However, the other side of the medal is that there is promise in a potential forward genetics approach offered by analysing phenotype-causative focal copy number variations (in particular full deletions) in 293-derived clones selected for adaptation to new growth conditions (such as high-cell density cultivation while producing biopharmaceuticals, virus infection, activation of particular signal transduction pathways and so on). This approach is made possible by the apparent property of 293 cells to have lost control over chromosomal structure to a great extent. Consequently, a culture of 293 cells should be considered as an entire 'population' of individual cells with different chromosomal structure makeup. Copy number variations are easy to identify at high resolution using high-coverage resequencing. Further experimentation will reveal whether phenotype-selected copy number variations can always be distinguished from such variations that occur randomly. In this perspective, genomic diversity of the 293 cell line might prove to be an experimental opportunity and might further enhance its role as a provider of knowledge on human cell biology. Methods Cell cultivation for DNA and RNA preparation All cell lines were cultured from frozen stocks at 37 °C in Dulbecco’s Modified Eagle Medium (DMEM; Invitrogen) supplemented with 10% (v/v) fetal calf serum, 2 mM L-glutamine, 100 U ml−1 penicillin G, 110 mg l−1 sodium pyruvate and 100 μg ml−1 streptomycin. All lines were routinely split twice a week, when ~80% confluency was reached. Depending on the cell line, the dilution was between 1:3 (293A) and 1:20 (293T). To prepare genomic DNA, ~30 million cells were harvested for each line. The genomic DNA was extracted and purified using the Gentra Puregene Cell kit (Qiagen GmbH, Hilden, Germany) with RNAse treatment of the samples, according to the manufacturer’s instructions. DNA concentrations were determined fluorimetrically with the Quant-iT PicoGreen dsDNA Reagent (Molecular Probes, Life Technologies Ltd., Paisley, UK). For RNA preparation, the cell lines were cultured in 75-cm2 filter cap flasks in a humidified, 8% CO2 atmosphere incubator in DMEM/Ham’s F12 (DMEM/F12; Invitrogen) supplemented with 10% (v/v) fetal calf serum, 2 mM L-glutamine, 100 U ml−1 penicillin G and 100 μg ml−1 streptomycin. Flask positions in the incubator were randomized daily to correct for potential temperature biases. Total RNA was extracted from three replicates of each cell line using Qiagen’s RNeasy Midi kit according to the manufacturer’s instructions, including an on-column DNase-I digest. Concentrations were determined with a NanoDrop ND-1000 spectrophotometer (Thermo Scientific), and RNA quality was assessed on a 2100 Bioanalyzer using RNA 6000 Pico chips (Agilent Technologies). All samples had an RNA integrity number of 9.5 or better. For the RT–qPCR validation of miRNA expression levels, procedures were identical except that the small RNAs were isolated using the miRCURY RNA isolation kit Cell and Plant (Exiqon), again according to the manufacturer’s instructions. Exon arrays After spiking total RNA from each cell line with bacterial poly-A RNA-positive controls (Affymetrix), every sample was reverse-transcribed, converted to double-stranded cDNA, in vitro-transcribed and amplified using the Ambion WT Expression Kit. The obtained single-stranded cDNA was biotinylated after fragmentation with the Affymetrix WT Terminal Labeling kit as outlined in the manufacturer’s instructions. The resulting samples were mixed with hybridization controls (Affymetrix) and hybridized on GeneChip Human Exon 1.0 ST Arrays (Affymetrix). The arrays were stained and washed in a GeneChip Fluidics Station 450 (Affymetrix) and scanned for raw probe signal intensities with the GeneChip Scanner 3000 (Affymetrix). For the processing of the data, see extended experimental procedures. Exon-array data analysis We used a combination of the R Statistical Software Package (www.r-project.org) and Affymetrix Power Tools (APT; Affymetrix) for the quality control and differential expression analysis of the exon-array data, partly as described earlier43. The full R code and APT commands are available as in Supplementary Data 9 and 10). Briefly, exon- and gene-level intensity estimates were generated by background correction, normalization and probe summarization using the robust multi-array average algorithm with APT. At the gene level, after quality control of the raw data in R, genes of which the expression was undetected in all six lines were removed from further analysis, as were the genes of which expression was below the estimated noise level in all lines. This noise level threshold was set at the signal intensity level that eliminated ‘detection’ of expression of more than 95% of the genes on the Y-chromosome, which is absent from the HEK293 lineage (which was derived from a female embryo) and thus serves as an appropriate internal negative control. Differential gene expression analysis was performed for the relevant cell line pairs using a linear model fit implemented in the R Bioconductor package Limma44, considering only core probe sets. The Benjamini–Hochberg (BH) method was applied to correct for multiple testing. Lists of significantly up- and downregulated genes (BH-adjusted P values<0.01) with a minimal twofold change in expression were subjected to functional enrichment analysis using DAVID45 and IPA (Ingenuity Systems, www.ingenuity.com), transcription factor regulation prediction using DiRE46 and manual inspection. Those lists are available as Supplementary Materials. For integration in the IGV genome browser36, we chose to display all genes found to be differentially expressed (BH-adjusted P value<0.01) in the pairwise comparison of interest, irrespective of their log2-fold change, which is displayed as a function of the bar height. The ‘web link to gene expression data’ track links every gene of which expression was detected to a table with the statistical details. The mean exon expression values in the IGV ‘mean probe set expression’ tracks represent the log2 signal values of the filtered extended exon probe sets, that is, after removal of undetected, cross-hybridizing and noisy probes. CG sequencing and analysis Anticipating the pseudotriploidy of the HEK293 genome, genomic DNA from each cell line was submitted to CG’s sequencing service16 (detailed in Supplementary Information) with the request to maximize the sequencing machine’s output to achieve the highest coverage possible, yielding 158~287 Gb of mapped reads of which 122~190 Gb of reads mapped with an expected paired distance (Supplementary Tables 1 and 2). The raw data were analysed with version 1.11 of the company’s analysis software and processed with CGAtools v1.5 (http://cgatools.sourceforge.net/). This pipeline entails read mapping followed by local reassembly of reads that map to a region in which deviation from the reference sequence is suspected from the mapping results. This is then used as the input for SNP and small indel calling. A second analysis focuses on copy number variation (see Supplementary Note 1) and uses the genome-normalized average sequence coverage as input, together with the genome-normalized sequence coverage of 46 normal diploid human genome-resequencing data sets (baseline genome) for the area under analysis. These latter data are used to correct the coverage for sequence-specific biases in the sequencing workflow. The output of this analysis is 2-kbp-resolution copy number expressed as a factor relative to a copy number of 2. As described in the main text, we derived true copy number from these data through calibration with genome-weighted average ploidy as derived from Illumina SNP-array data (Supplementary Table 6). A third analysis uses the paired-end reads of which the mate pairs do not map to a continuous stretch of the human reference genome sequence, and which thus provide evidence for chromosomal rearrangements. These reads are de novo assembled into ‘junction sequence contigs’ that contain the information about the breakpoints involved in such chromosomal rearrangements. The CG raw data and initial analysis results were processed by CGAtools v1.5 (http://cgatools.sourceforge.net/) with scripts from the CG user community tool repository and our in-house scripts (see Supplementary Note 2). To enable independent analysis of the data, we mapped the sequencing reads to the human reference genome, build hg18, using RTG Investigator from Real Time Genomics (http://www.realtimegenomics.com/) with default setting (maximum mate-pair insert size: 1,000, minimum insert size 0 and report the maximum best five matches). Upon mapping, SNP and small indel calling were also performed using the RTG software Investigator. Only SNP/indels passing the quality filter (called in more than half of the reads and covered by less than 200 × coverage to avoid variant calling in highly repetitive regions) were kept for further analysis. The lists of SNPs and indels called either by CG or RTG were merged by vcftools47. To remove platform-specific artifacts from the CG sequencing, the extended variant list was filtered using ANNOVAR48, to remove variants located in a region where less than 30% of the CG69 data sets had sequencing information. We then functionally annotated this filtered extended variant list by ANNOVAR. We used GenomeComb (http://genomecomb.sourceforge.net/) to reformat the SNV calling results from CG for the six cell lines49. In order to increase the number of concordants between cell lines and reduce the false-positive SNV calling rate, we used the obligatory filtering strategy: remove uncertain calls and filtered based on the variant score reported from CG in each cell line. Variant scores lower than the reported average variant score were removed. The SVs detected from CG analysis were first filtered with cgatools against the publicly available Yoruban (NA19238) CG genome data set, to remove frequently occurring SVs. SVs in the 293-derived cell lines were further filtered against the 293 line and we only retained those with low frequency (<10%) in the CG69 population for further manual inspection. SNP-array procedures Genomic DNA (same sample as used for genome sequencing) of each cell line was analysed using the Illumina HumanCytoSNP-12 v2.1 SNP-array, entirely according to the manufacturer’s instructions. For analysis, we used the ASCAT algorithm, which accurately determines allele-specific copy numbers in tumours and aneuploid cell lines by estimating and adjusting for overall ploidy and effective tumour fraction in the sample50. ASCAT uses the raw BAF and logR data of the Illumina HumanCytoSNP-12 v2.1. Author contributions Y.-C.L. designed experiments and analysed the CG data, plasmid insertion site detection and data integration, under the scientific supervision of Y.V.d.P. M.B. carried out exon-array experiments and data analysis, qPCR validation of array data and general data mining. L.M. conducted mitochondrial haplotype study, PCR validation of plasmid insertion sites and general data mining. I.L. performed general data mining, under the scientific supervision of J.T. N.V.R. carried out multiplex fluorescence in situ hybridization data generation and analysis, under the scientific supervision of F.S. A.S.: 293 Variant Viewer website construction. J.R. assisted in GenomeComb analysis of CG data. M.M. carried out SNP arrays, under the scientific supervision of D.L. S.P. helped with general bioinformatics assistance. R.D. and J.C. performed CG data acquisition. N.C. carried out project initiation and design and scientific supervision. N.C., M.B., Y.C.-L. and L.M. co-wrote the manuscript. Additional information Accession codes: Complete Genomics sequencing data have been deposited in the European Nucleotide Archive (ENA) under the accession code PRJEB3209. The Affymetrix exon-array data have been deposited in the ArrayExpress Archive under the accession code E-MEXP-3516. How to cite this article: Lin, Y.-C. et al. Genome dynamics of the human embryonic kidney 293 lineage in response to cell biology manipulations. Nat. Commun. 5:4767 doi: 10.1038/ncomms5767 (2014). Supplementary Material Supplementary Information Supplementary Figures 1-11, Supplementary Tables 1-7, Supplementary Notes 1-5, Supplementary Methods and Supplementary References Supplementary Data 1 M-FISH (Multicolor-FISH) analysis Supplementary Data 2 Manually annotated single-nucleotide polymorphisms (SNPs) and structural variants (SVs) Supplementary Data 3 Copy number variation per chromosome Supplementary Data 4 Approximate copy number of citric acid cycle-related genes in the 6 sequenced cell lines Supplementary Data 5 Gene-level differential expression Supplementary Data 6 Copy number variation, minor allele frequency, ploidy level determined by CG and SNP array on multiple passages Supplementary Data 7 Overview of plasmid insertion sites Supplementary Data 8 Sanger sequencing validation of plasmid insertion sites Supplementary Data 9 R code used for gene level differential expression analysis Supplementary Data 10 R code used for exon level differential expression analysis
                Bookmark

                Author and article information

                Journal
                J Enzyme Inhib Med Chem
                J Enzyme Inhib Med Chem
                Journal of Enzyme Inhibition and Medicinal Chemistry
                Taylor & Francis
                1475-6366
                1475-6374
                31 August 2022
                2022
                31 August 2022
                : 37
                : 1
                : 2370-2381
                Affiliations
                Department of Chemistry, Center for Diagnostics and Therapeutics, Georgia State University , Atlanta, GA, USA
                Author notes

                Supplemental data for this article is available online at https://doi.org/10.1080/14756366.2022.2116699

                CONTACT Binghe Wang wang@ 123456gsu.edu Department of Chemistry, Center for Diagnostics and Therapeutics, Georgia State University , Atlanta, GA, 30303, USA
                Xiaoxiao Yang xyang20@ 123456gsu.edu Department of Chemistry, Center for Diagnostics and Therapeutics, Georgia State University , Atlanta, GA, 30303, USA
                Author information
                https://orcid.org/0000-0002-2200-5270
                Article
                2116699
                10.1080/14756366.2022.2116699
                9448394
                36043494
                dc146bcd-8253-421d-bb6d-fcefbdbec9e5
                © 2022 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.

                This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

                History
                Page count
                Figures: 5, Tables: 3, Pages: 12, Words: 9694
                Categories
                Research Article
                Research Paper

                Pharmaceutical chemistry
                p53,mdm2,anthraquinone,cancer,structure–activity relationship (sar)
                Pharmaceutical chemistry
                p53, mdm2, anthraquinone, cancer, structure–activity relationship (sar)

                Comments

                Comment on this article