Suppression of CTC1 inhibits hepatocellular carcinoma cell growth and enhances RHPS4 cytotoxicity.

BACKGROUND: Although DNA repair mechanisms function to maintain genomic integrity, in cancer cells these mechanisms may negatively affect treatment efficiency. The strategy of targeting cancer cells via inhibiting DNA damage repair has been successfully used in breast and ovarian cancer using PARP inhibitors. Unfortunately, such strategies have not yet yielded results in liver cancer. Hepatocellular carcinoma (HCC), the most common type of liver cancer, is a treatment-resistant malignancy. Despite the development of guided therapies, treatment regimens for advanced HCC patients still fall short of the current need and significant problems such as cancer relapse with resistance still exist. In this paper, we targeted telomeric replication protein CTC1, which is responsible for telomere maintenance. METHODS: CTC expression was analyzed using tumor and matched-tissue RNA-sequencing data from TCGA and GTEx. In HCC cell lines, q-RT-PCR and Western blotting were used to detect CTC1 expression. The knock-down of CTC1 was achieved using lentiviral plasmids. The effects of CTC1 silencing on HCC cells were analyzed by flow cytometry, MTT, spheroid and colony formation assays. RESULTS: CTC1 is significantly downregulated in HCC tumor samples. However, CTC1 protein levels were higher in sorafenib-resistant cell lines compared to the parental groups. CTC1 inhibition reduced cell proliferation in sorafenib-resistant HCC cell lines and diminished their spheroid and colony forming capacities. Moreover, these cells were more sensitive to single and combined drug treatment with G4 stabilizer RHPS4 and sorafenib. CONCLUSION: Our results suggest that targeting CTC1 might be a viable option for combinational therapies designed for sorafenib resistant HCC patients.

Identification of PPA1 inhibitor candidates for potential repurposing in cancer medicine.

Inorganic pyrophosphatase 1 (PPA1) is pivotal to cellular metabolism as it facilitates the hydrolysis of PPi-a by-product of various metabolic processes that influence cell growth and differentiation. Overexpression of PPA1 enzyme has been linked to diminished patient survival and was shown to influence tumor cell dynamics, thereby positioning it as a potential therapy target for a variety of cancers including colorectal cancer, diffuse large B-cell lymphoma, and lung adenocarcinoma. Despite this therapeutic promise, there are no known inhibitors of PPA1 as of today. In this study, we searched for potential PPA1 inhibitors using a molecular docking screen of 30 470 compounds with a history of clinical trials and/or US Food and Drug Administration approval. We specifically targeted the active pocket that coincides with the established catalytic domain. Our screen identified promising hits, which we further subjected to ADMET (absorption, distribution, metabolism, excretion, and toxicity) filtering. Subsequent molecular dynamics (MD) analyses were conducted on devazepide, quinotolast, and tarazepide-the three substances that successfully navigated all filters. MD analyses reinforced the stability of the protein-ligand complexes and confirmed ligand binding, as substantiated by our root mean square deviation, radius of gyration and secondary structures of proteins analyses. Furthermore, Molecular Mechanics Poisson-Boltzmann Surface Area calculations post-MD identified devazepide and quinotolast as showing higher binding affinities; being supported by principal component analysis, free energy landscape, and dynamic cross-correlation matrix results. Overall, our study reveals devazepide and quinotolast as potential candidates for PPA1 inhibition which could be considered for repurposing studies that need further experimental validation. These results not only reveal a potential for clinical repurposing for PPA1 inhibition but they also offer valuable insights into the development of future compounds for targeting the crucial PPA1 enzyme.

A systematic molecular and pharmacologic evaluation of AKT inhibitors reveals new insight into their biological activity.

BACKGROUND: AKT, a critical effector of the phosphoinositide 3-kinase (PI3K) signalling cascade, is an intensely pursued therapeutic target in oncology. Two distinct classes of AKT inhibitors have been in clinical development, ATP-competitive and allosteric. Class-specific differences in drug activity are likely the result of differential structural and conformational requirements governing efficient target binding, which ultimately determine isoform-specific potency, selectivity profiles and activity against clinically relevant AKT mutant variants. METHODS: We have carried out a systematic evaluation of clinical AKT inhibitors using in vitro pharmacology, molecular profiling and biochemical assays together with structural modelling to better understand the context of drug-specific and drug-class-specific cell-killing activity. RESULTS: Our data demonstrate clear differences between ATP-competitive and allosteric AKT inhibitors, including differential effects on non-catalytic activity as measured by a novel functional readout. Surprisingly, we found that some mutations can cause drug resistance in an isoform-selective manner despite high structural conservation across AKT isoforms. Finally, we have derived drug-class-specific phosphoproteomic signatures and used them to identify effective drug combinations. CONCLUSIONS: These findings illustrate the utility of individual AKT inhibitors, both as drugs and as chemical probes, and the benefit of AKT inhibitor pharmacological diversity in providing a repertoire of context-specific therapeutic options.

Expression dynamics of repetitive DNA in early human embryonic development.

BACKGROUND: The last decade witnessed a number of genome-wide studies on human pre-implantation, which mostly focused on genes and provided only limited information on repeats, excluding the satellites. Considering the fact that repeats constitute a large portion of our genome with reported links to human physiology and disease, a thorough understanding of their spatiotemporal regulation during human embryogenesis will give invaluable clues on chromatin dynamics across time and space. Therefore, we performed a detailed expression analysis of all repetitive DNA elements including the satellites across stages of human pre-implantation and embryonic stem cells. RESULTS: We uncovered stage-specific expressions of more than a thousand repeat elements whose expressions fluctuated with a mild global decrease at the blastocyst stage. Most satellites were highly expressed at the 4-cell level and expressions of ACRO1 and D20S16 specifically peaked at this point. Whereas all members of the SVA elements were highly upregulated at 8-cell and morula stages, other transposons and small RNA repeats exhibited a high level of variation among their specific subtypes. Our repeat enrichment analysis in gene promoters coupled with expression correlations highlighted potential links between repeat expressions and nearby genes, emphasising mostly 8-cell and morula specific genes together with SVA_D, LTR5_Hs and LTR70 transposons. The DNA methylation analysis further complemented the understanding on the mechanistic aspects of the repeatome's regulation per se and revealed critical stages where DNA methylation levels are negatively correlating with repeat expression. CONCLUSIONS: Taken together, our study shows that specific expression patterns are not exclusive to genes and long non-coding RNAs but the repeatome also exhibits an intriguingly dynamic pattern at the global scale. Repeats identified in this study; particularly satellites, which were historically associated with heterochromatin, and those with potential links to nearby gene expression provide valuable insights into the understanding of key events in genomic regulation and warrant further research in epigenetics, genomics and developmental biology.

Heterochromatinization induced by GAA-repeat hyperexpansion in Friedreich's ataxia can be reduced upon HDAC inhibition by vitamin B3.

Large intronic expansions of the triplet-repeat sequence (GAA.TTC) cause transcriptional repression of the Frataxin gene (FXN) leading to Friedreich's ataxia (FRDA). We previously found that GAA-triplet expansions stimulate heterochromatinization in vivo in transgenic mice. We report here using chromosome conformation capture (3C) coupled with high-throughput sequencing that the GAA-repeat expansion in FRDA cells stimulates a higher-order structure as a fragment containing the GAA-repeat expansion showed an increased interaction frequency with genomic regions along the FXN locus. This is consistent with a more compacted chromatin and coincided with an increase in both constitutive H3K9me3 and facultative H3K27me3 heterochromatic marks in FRDA. Consistent with this, DNase I accessibility in regions flanking the GAA repeats in patients was decreased compared with healthy controls. Strikingly, this effect could be antagonized with the class III histone deactylase (HDAC) inhibitor vitamin B3 (nicotinamide) which activated the silenced FXN gene in several FRDA models. Examination of the FXN locus revealed a reduction of H3K9me3 and H3K27me3, an increased accessibility to DNase I and an induction of euchromatic H3 and H4 histone acetylations upon nicotinamide treatment. In addition, transcriptomic analysis of nicotinamide treated and untreated FRDA primary lymphocytes revealed that the expression of 67% of genes known to be dysregulated in FRDA was ameliorated by the treatment. These findings show that nictotinamide can up-regulate the FXN gene and reveal a potential mechanism of action for nicotinamide in reactivating the epigenetically silenced FXN gene and therefore support the further assessment of HDAC inhibitors (HDACi's) in FRDA and diseases caused by a similar mechanism.

Epigenetic and neurological effects and safety of high-dose nicotinamide in patients with Friedreich's ataxia: an exploratory, open-label, dose-escalation study.

BACKGROUND: Friedreich's ataxia is a progressive degenerative disorder caused by deficiency of the frataxin protein. Expanded GAA repeats within intron 1 of the frataxin (FXN) gene lead to its heterochromatinisation and transcriptional silencing. Preclinical studies have shown that the histone deacetylase inhibitor nicotinamide (vitamin B3) can remodel the pathological heterochromatin and upregulate expression of FXN. We aimed to assess the epigenetic and neurological effects and safety of high-dose nicotinamide in patients with Friedreich's ataxia. METHODS: In this exploratory, open-label, dose-escalation study in the UK, male and female patients (aged 18 years or older) with Friedreich's ataxia were given single doses (phase 1) and repeated daily doses of 2-8 g oral nicotinamide for 5 days (phase 2) and 8 weeks (phase 3). Doses were gradually escalated during phases 1 and 2, with individual maximum tolerated doses used in phase 3. The primary outcome was the upregulation of frataxin expression. We also assessed the safety and tolerability of nicotinamide, used chromatin immunoprecipitation to investigate changes in chromatin structure at the FXN gene locus, and assessed the effect of nicotinamide treatment on clinical scales for ataxia. This study is registered with ClinicalTrials.gov, number NCT01589809. FINDINGS: Nicotinamide was generally well tolerated; the main adverse event was nausea, which in most cases was mild, dose-related, and resolved spontaneously or after dose reduction, use of antinausea drugs, or both. Phase 1 showed a dose-response relation for proportional change in frataxin protein concentration from baseline to 8 h post-dose, which increased with increasing dose (p=0·0004). Bayesian analysis predicted that 3·8 g would result in a 1·5-times increase and 7·5 g in a doubling of frataxin protein concentration. Phases 2 and 3 showed that daily dosing at 3·5-6 g resulted in a sustained and significant (p<0·0001) upregulation of frataxin expression, which was accompanied by a reduction in heterochromatin modifications at the FXN locus. Clinical measures showed no significant changes. INTERPRETATION: Nicotinamide was associated with a sustained improvement in frataxin concentrations towards those seen in asymptomatic carriers during 8 weeks of daily dosing. Further investigation of the long-term clinical benefits of nicotinamide and its ability to ameliorate frataxin deficiency in Friedreich's ataxia is warranted. FUNDING: Ataxia UK, Ataxia Ireland, Association Suisse de l'Ataxie de Friedreich, Associazione Italiana per le Sindromi Atassiche, UK National Institute for Health Research, European Friedreich's Ataxia Consortium for Translational Studies, and Imperial Biomedical Research Centre.

Thyroid hormone T3 augments the cytotoxicity of sorafenib in Huh7 hepatocellular carcinoma cells by suppressing AKT expression.

BACKGROUND AND OBJECTIVES: Hepatocellular carcinoma (HCC) is a primary cancer that poorly responds to treatment. Molecular cancer studies led to the development of kinase inhibitors, among which sorafenib stands out as a multi-kinase inhibitor approved by FDA for first line use in HCC patients. However, the efficiency of sorafenib was shown to be counteracted by numerous subcellular pathways involving the effector kinase AKT, causing resistance and limiting its survival benefit. On the way of breaking such resistance mechanisms and increase the efficiency of sorafenib, deeper understanding of hepatocellular physiology is essential. Thyroid hormones were shown to be metabolized in liver and inevitably affect the molecular behaviour of hepatocytes. Interestingly, thyroid hormone T3 was also demonstrated to be potentially influential in liver regeneration and treatment with this hormone reportedly led to a decrease in HCC tumor growths. In this study, we aimed to uncover the impact of T3 hormone on the cytotoxic response to sorafenib in HCC in vitro. MATERIALS AND METHODS: We pre-treated the HCC cell line Huh-7 with T3 prior to sorafenib exposure both in 2D and 3D culture. We checked cell viability with MTT assay in 2D culture and measured the sizes of 3D spheroids with bright-field microscopy followed by a surface analysis with ImageJ. We also performed scratch assay to measure cell migration as well as western blot and qPCR to uncover affected pathways. RESULTS: We observed an additive effect to sorafenib's cytotoxicity both in 2D and 3D culture. Cell migration assay also confirmed our finding and pointed out a benefit of T3 hormone in HCC cell migration. Western blot experiments showed that T3 exerts its additive effect by suppressing AKT expression upon sorafenib treatment both at protein and gene expression levels. CONCLUSION: Our results open a promising new avenue in increasing sorafenib's cytotoxicity where thyroid hormone T3 is utilized to modulate AKT expression to combat resistance, and warrant further studies in the field.

Gene regulation and epigenetics in Friedreich's ataxia.

This is an exciting time in the study of Friedreich's ataxia. Over the last 10 years much progress has been made in uncovering the mechanisms, whereby the Frataxin gene is silenced by (GAA)n repeat expansions and several of the findings are now ripe for testing in the clinic. The discovery that the Frataxin gene is heterochromatinised and that this can be antagonised in vivo has led to the tantalizing possibility that the disease might be amenable to a more radical therapeutic approach involving epigenetic modifiers. Here, we set out to review progress in the understanding of the fundamental mechanisms whereby genes are regulated at this level and how these findings have been applied to achieve a deeper understanding of the dysregulation that occurs as the primary genetic lesion in Friedreich's ataxia.

In silico drug screen reveals potential competitive MTHFR inhibitors for clinical repurposing.

MTHFR (Methylenetetrahydrofolate reductase) is a pivotal enzyme involved in one-carbon metabolism, which is critical for the proliferation of cancer cells. In line with this, published literature showed that MTHFR knockdown caused impaired growth of multiple types of cancer cells. Moreover, higher MTHFR expression levels were linked to shorter overall survival in hepatocellular carcinoma, adrenocortical carcinoma, and low-grade glioma, bringing the need to design MTHFR inhibitors as a possible treatment option. No competitive inhibitors of MTHFR have been reported as of today. This study aimed to identify potential competitive MTHFR inhibitor candidates using an in silico drug screen. A total of 30470 molecules containing biogenic compounds, FDA-approved drugs, and those in clinical trials were screened against the catalytic pocket of MTHFR in the presence and absence of cofactors. Binding energy and ADMET analysis revealed that Vilanterol (ß2-adrenergic agonist), Selexipag (prostacyclin receptor agonist), and Ramipril Diketopiperazine (ACE inhibitor) are potential competitive inhibitors of MTHFR. Molecular dynamics analyses and MM-PBSA calculations with these compounds particularly revealed the amino acids between 285-290 for ligand binding and highlighted Vilanterol as the strongest candidate for MTHFR inhibition. Our results could guide the development of novel MTHFR inhibitor compounds, which could be inspired by the drugs brought into the spotlight here. More importantly, these potential candidates could be quhickly tested as a repurposing strategy in pre-clinical and clinical studies of the cancers mentioned above.Communicated by Ramaswamy H. Sarma.

Repeat expression is linked to patient survival and exhibits single nucleotide variation in pancreatic cancer revealing LTR70:r.879A>G.

Despite an overwhelming number of cancer literature reporting the links between patient survival and the expression levels of genes or mutations/single nucleotide variations (SNVs) on them, there is only limited information on repeat elements, which make at least half the human genome. Here, we analysed RNA-seq data obtained from primary pancreatic cancer tissues of 51 patients and revealed that two transposons, HERVI-int and X6A_LINE, showed an upregulation trend in the patients who lived shorter, along with 56 other potential repeats which were linked to survival. We also detected expressed single nucleotide variations (SNVs) on repeats, among which LTR70:r.879A>G stands out with the effect of its presence on this particular repeat's expression levels and a significant link to overall patient survival. Interestingly, the expression of LTR70:r.879A>G correlated with different cancer genes in comparison to its reference version highlighting the involvement of BRAF and Fumerate Hydratase with this expressed SNV. This is one of the first studies revealing possible links between repeat expression and survival in cancer and it warrants further research in this avenue.

Comparative molecular dynamics analyses on PIK3CA hotspot mutations with PI3Kα specific inhibitors and ATP.

PI3K pathway is heavily emphasized in cancer where PIK3CA, which encodes for the p110α subunit of PI3Kα, presents itself as the second most common mutated gene. A lot of effort has been put in developing PI3K inhibitors, opening promising avenues for the treatment of cancer. Among these, PI3Kα specific inhibitor alpelisib was approved by FDA for breast cancer and other α-isoform specific inhibitors such as inavolisib and serabelisib reached clinical trials. However, the mode of action of these inhibitors on mutated PI3Kα and how they interact with mutant structures has not been fully elucidated yet. In this study, we are revealing the calculated interactions and binding affinities of these inhibitors within the context of PIK3CA hotspot mutations (E542K, E545K and H1047R) by employing molecular dynamics (MD) simulations. We performed principal component analysis to understand the motions of the protein complex during our simulations and also checked the correlated motions of all amino acids. Binding affinity calculations with MM-PBSA confirmed the consistent binding of alpelisib across mutations and revealed relatively higher affinities for inavolisib towards wild-type and H1047R mutant structures in comparison to other inhibitors. On the other hand, E542K mutation significantly impaired the interaction of inavolisib and serabelisib with PI3Kα. We also investigated the structural relationship of the natural ligand ATP with PI3Kα, and interestingly realized a significant reduction in binding affinity for the mutants, with potentially unexpected implications on the mechanisms that render these mutations oncogenic. Moreover, correlated motions of all residues were generally higher for ATP except the H1047R mutation which exhibited a distinguishable reduction. The results presented here could be guiding for pre-clinical and clinical studies of personalized medicine where individual mutations are a strong consideration point.

Identification of differentially expressed genomic repeats in primary hepatocellular carcinoma and their potential links to biological processes and survival.

Hepatocellular carcinoma (HCC) is one of the deadliest cancers. Research on HCC so far primarily focused on genes and provided limited information on genomic repeats, which constitute more than half of the human genome and contribute to genomic stability. In line with this, repeat dysregulation was significantly shown to be pathological in various cancers and other diseases. In this study, we aimed to determine the full repeat expression profile of HCC for the first time. We utilised two independent RNA-seq datasets obtained from primary HCC tumours with matched normal tissues of 20 and 17 HCC patients, respectively. We quantified repeat expressions and analysed their differential expression. We also identified repeats that are cooperatively expressed with genes by constructing a gene coexpression network. Our results indicated that HCC tumours in both datasets harbour 24 differentially expressed repeats and even more elements were coexpressed with genes involved in various metabolic pathways. We discovered that two L1 elements (L1M3b, L1M3de) were downregulated and a handful of HERV subfamily repeats (HERV-Fc1-int, HERV3-int, HERVE_a-int, HERVK11D-int, HERVK14C-int, HERVL18-int) were upregulated with the exception of HERV1_LTRc, which was downregulated. Various LTR elements (LTR32, LTR9, LTR4, LTR52-int, LTR70) and MER elements (MER11C, MER11D, MER57C1, MER9a1, MER74C) were implicated along with few other subtypes including Charlie12, MLT2A2, Tigger15a, Tigger 17b. The only satellite repeat differentially expressed in both datasets was GSATII, whose expression was upregulated in 33 (>90%) out of 37 patients. Notably, GSATII expression correlated with HCC survival genes. Elements discovered here promise future studies to be considered for biomarker and HCC therapy research. The coexpression pattern of the GSATII satellite with HCC survival genes and the fact that it has been upregulated in the vast majority of patients make this repeat particularly stand out for HCC.

Signature changes in the expressions of protein-coding genes, lncRNAs, and repeat elements in early and late cellular senescence.

Replicative cellular senescence is the main cause of aging. It is important to note that early senescence is linked to tissue regeneration, whereas late senescence is known to trigger a chronically inflammatory phenotype. Despite the presence of various genome-wide studies, there is a lack of information on distinguishing early and late senescent phenotypes at the transcriptome level. Particularly, the changes in the noncoding RNA portion of the aging cell have not been fully elucidated. By utilising RNA sequencing data of fibroblasts, hereby, we are not only reporting changes in gene expression profiles and relevant biological processes in the early and late senescent phenotypes but also presenting significant differences in the expressions of many unravelled long noncoding RNAs (lncRNAs) and transcripts arisen from repetitive DNA. Our results indicate that, in addition to previously reported L1 elements, various LTR and DNA transposons, as well as members of the classical satellites including HSAT5 and α-satellites (ALR/Alpha), are expressed at higher levels in late senescence. Moreover, we revealed finer links between the expression levels of repeats with the genes located near them and known to be involved in cell cycle and senescence. Noncoding elements reported here provide a new perspective to be explored in further experimental studies.

Dysregulated expression of repetitive DNA in ER+/HER2- breast cancer.

Limited studies on breast cancer indicated pathogenic changes in the expressions of some repeat elements. A global analysis was much needed within this context to distinguish the most significant repeats from more than a thousand repeat motifs. Utilising a previously presented RNA-seq dataset, we studied expression changes of all repeats in ER+/HER2- human breast tumour samples obtained from 22 patients in comparison to matched normal tissues. Fifty six (56) repeat subtypes including satellites and transposons were found to be differentially expressed and most of them were novel for breast cancer. HERVKC4-int and HERV1_LTRc, whose expressions correlated well with that of the estrogen receptor gene ESR1, were upregulated at the highest level. REP522 and D20S16 satellites were also significantly upregulated along with insignificant increases in the expressions of other satellites including HSATI and BSR/beta. Interestingly, expressions of REP522 and D20S16 correlated with many key breast cancer pathway (e.g. BRCA1, BRCA2, AKT1, MTOR, KRAS) and survival genes; possibly highlighting their importance in the carcinogenesis of breast. Additional differentially expressed elements such as L1P and various MER transposons also exhibited a similar pattern. Finally, our repeat enrichment analysis on the promoters of differentially expressed genes revealed further links between additional repeats and nearby genes.

TEffectR: an R package for studying the potential effects of transposable elements on gene expression with linear regression model.

INTRODUCTION: Recent studies highlight the crucial regulatory roles of transposable elements (TEs) on proximal gene expression in distinct biological contexts such as disease and development. However, computational tools extracting potential TE -proximal gene expression associations from RNA-sequencing data are still missing. IMPLEMENTATION: Herein, we developed a novel R package, using a linear regression model, for studying the potential influence of TE species on proximal gene expression from a given RNA-sequencing data set. Our R package, namely TEffectR, makes use of publicly available RepeatMasker TE and Ensembl gene annotations as well as several functions of other R-packages. It calculates total read counts of TEs from sorted and indexed genome aligned BAM files provided by the user, and determines statistically significant relations between TE expression and the transcription of nearby genes under diverse biological conditions. AVAILABILITY: TEffectR is freely available at https://github.com/karakulahg/TEffectR along with a handy tutorial as exemplified by the analysis of RNA-sequencing data including normal and tumour tissue specimens obtained from breast cancer patients.

Transcriptional Activation of Pericentromeric Satellite Repeats and Disruption of Centromeric Clustering upon Proteasome Inhibition.

Heterochromatinisation of pericentromeres, which in mice consist of arrays of major satellite repeats, are important for centromere formation and maintenance of genome stability. The dysregulation of this process has been linked to genomic stress and various cancers. Here we show in mice that the proteasome binds to major satellite repeats and proteasome inhibition by MG132 results in their transcriptional de-repression; this de-repression is independent of cell-cycle perturbation. The transcriptional activation of major satellite repeats upon proteasome inhibition is accompanied by delocalisation of heterochromatin protein 1 alpha (HP1α) from chromocentres, without detectable change in the levels of histone H3K9me3, H3K4me3, H3K36me3 and H3 acetylation on the major satellite repeats. Moreover, inhibition of the proteasome was found to increase the number of chromocentres per cell, reflecting destabilisation of the chromocentre structures. Our findings suggest that the proteasome plays a role in maintaining heterochromatin integrity of pericentromeres.

Sexual dimorphism in mammalian autosomal gene regulation is determined not only by Sry but by sex chromosome complement as well.

Differences between males and females are normally attributed to developmental and hormonal differences between the sexes. Here, we demonstrate differences between males and females in gene silencing using a heterochromatin-sensitive reporter gene. Using "sex-reversal" mouse models with varying sex chromosome complements, we found that this differential gene silencing was determined by X chromosome complement, rather than sex. Genome-wide transcription profiling showed that the expression of hundreds of autosomal genes was also sensitive to sex chromosome complement. These genome-wide analyses also uncovered a role for Sry in modulating autosomal gene expression in a sex chromosome complement-specific manner. The identification of this additional layer in the establishment of sexual dimorphisms has implications for understanding sexual dimorphisms in physiology and disease.