Two critical positions in zinc finger domains are heavily mutated in three human cancer types

Intracellular pH differentially regulates transcription of metabolic and signaling pathways in normal epithelial cells

Intracellular pH (pHi) dynamics regulate normal cell function, and dysregulated pHi dynamics is an emerging hallmark of cancer (constitutively increased pHi) and neurodegeneration (constitutively decreased pHi). However, the molecular mechanisms by which pHi dynamics regulate cell biology are poorly understood. Here, we discovered that altering pHi in normal human breast epithelial cells triggers global transcriptional changes. We identified 176 genes differentially regulated by pHi, with pHi-dependent genes clustering in signaling and glycolytic pathways. Using various normal epithelial cell models, we showed pH-dependent Notch homolog 1 protein expression, with increased protein abundance at high pHi. This resulted in pH-dependent downstream signaling, with increased Notch homolog 1 signaling at high pHi. We also found that high pHi increased the expression of glycolytic enzymes and regulators of pyruvate fate, including lactate dehydrogenase and pyruvate dehydrogenase kinase. These transcriptional changes were sufficient to alter lactate production, with high pHi shifting these normal epithelial cells toward a glycolytic metabolism and increasing lactate production. Thus, pHi dynamics transcriptionally regulate signaling and metabolic pathways in normal epithelial cells. Our data reveal new molecular regulators of pHi-dependent biology and a role for increased pHi in driving the acquisition of cancer-associated signaling and metabolic changes in normal human epithelial cells.

From driver genes to gene families: A computational analysis of oncogenic mutations and ubiquitination anomalies in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a widespread primary liver cancer with a high fatality rate. Despite several genes with oncogenic effects in HCC have been identified, many remain undiscovered. In this study, we conducted a comprehensive computational analysis to explore the involvement of genes within the same families as known driver genes in HCC. Specifically, we expanded the concept beyond single-gene mutations to encompass gene families sharing homologous structures, integrating various omics data to comprehensively understand gene abnormalities in cancer. Our analysis identified 74 domains with an enriched mutation burden, 404 domain mutation hotspots, and 233 dysregulated driver genes. We observed that specific low-frequency somatic mutations may contribute to HCC occurrence, potentially overlooked by single-gene algorithms. Furthermore, we systematically analyzed how abnormalities in the ubiquitinated proteasome system (UPS) impact HCC, finding that abnormal genes in E3, E2, DUB families, and Degron genes often result in HCC by affecting the stability of oncogenic or tumor suppressor proteins. In conclusion, expanding the exploration of driver genes to include gene families with homologous structures emerges as a promising strategy for uncovering additional oncogenic alterations in HCC.

Formerly degenerate seventh zinc finger domain from transcription factor ZNF711 rehabilitated by experimental NMR structure

Domain Z7 of nuclear transcription factor ZNF711 has the consensus last metal-ligand H23 found in odd-numbered zinc fingers of this protein replaced by a phenylalanine. Ever since the discovery of ZNF711, it has been thought that Z7 is probably non-functional because of the H23F substitution. The presence of H26 three positions downstream prompted us to examine if this histidine could substitute as the last metal-ligand. The Z7 domain adopts a stable tertiary structure upon metal-binding. The NMR structure of Zn2+-bound Z7 shows the classical ββα-fold of CCHH zinc fingers. Mutagenesis and pH titration experiments indicate that H26 is not involved in metal binding and that Z7 has a tridentate metal-binding site comprised of only residues C3, C6, and H19. By contrast, an F23H mutation that introduces a histidine in the consensus position forms a tetradentate ligand. The structure of the WT Z7 is stable causing restricted ring-flipping of phenylalanines 10 and 23. Dynamics are increased with either the H26A or F23H substitutions and aromatic ring rotation is no longer hindered in the two mutants. The mutations have only small effects on the Kd values for Zn2+ and Co2+ and retain the high thermal stability of the WT domain above 80°C. Like two previously reported designed zinc fingers with the last ligand replaced by water, the WT Z7 domain is catalytically active, hydrolyzing 4-nitrophenyl acetate. We discuss the implications of naturally occurring tridentate zinc fingers for cancer mutations and drug targeting of notoriously undruggable transcription factors.

Review Article on Molecular Basis of Zinc and Copper Interactions in Cancer Physiology

Various clinical manifestations associated with measurable abnormalities of Zn and Cu in serum and tissue were determined in Cancer-Patients (CP), and therefore, these two metals are drawing more and more attention presently than ever before. Cancer is a disease of uncontrolled-abnormal-cell-division with invasion-potential which was exhibited to occur due to dys-regulation/dys-homeostasis of fundamental-biological-pathways (FBP) including antioxidant-enzyme-defense-system, anti-inflammatory and immune-systems, and DNA-damage-repair-system in the human-body resulting in generation of chronic-oxidative-stress induced DNA-damage and gene-mutations, inflammation and compromised immune-system, tumor-induced increased angiogenesis, and inhibition of apoptosis processes. Zn and Cu were recognized to be the most crucial components of FBP and imbalance in Zn/Cu ratios in CP asserted to generate chronic toxicity in human body through various mechanisms including increased chronic oxidative stress linked compromised DNA integrity and gene mutations due to malfunctioning of DNA damage repair enzymes; increased angiogenesis process due to Zn- and Cu-binding proteins metallothionein and ceruloplasmin-induced enhanced expression of tumor growth factors; and elevation in inflammatory response which was further shown to down/upregulate gene expression of multiple Zn transporter proteins leading to dys-homeostasis of intracellular Zn concentrations, and it was determined to disturb the equilibrium between cell growth and division, proliferation, differentiation, and apoptosis processes which lead to cancer progression. Moreover, Zn was reported to affect matrix metalloproteinase activity and influence immune system cells to respond differently to different cytokines and enhance immune-suppressive effects accelerating the angiogenesis, invasion, and metastasis potential in cancer. Further, the most significant use of serum Cu/Zn ratio was recommended in clinical diagnosis, prognosis, tumor stage, patient survival, and cancer follow-up studies which need further investigations to elucidate and explore their roles in cancer physiology for clinical perspective.

Machine learning and multi-omics data reveal driver gene-based molecular subtypes in hepatocellular carcinoma for precision treatment

The heterogeneity of Hepatocellular Carcinoma (HCC) poses a barrier to effective treatment. Stratifying highly heterogeneous HCC into molecular subtypes with similar features is crucial for personalized anti-tumor therapies. Although driver genes play pivotal roles in cancer progression, their potential in HCC subtyping has been largely overlooked. This study aims to utilize driver genes to construct HCC subtype models and unravel their molecular mechanisms. Utilizing a novel computational framework, we expanded the initially identified 96 driver genes to 1192 based on mutational aspects and an additional 233 considering driver dysregulation. These genes were subsequently employed as stratification markers for further analyses. A novel multi-omics subtype classification algorithm was developed, leveraging mutation and expression data of the identified stratification genes. This algorithm successfully categorized HCC into two distinct subtypes, CLASS A and CLASS B, demonstrating significant differences in survival outcomes. Integrating multi-omics and single-cell data unveiled substantial distinctions between these subtypes regarding transcriptomics, mutations, copy number variations, and epigenomics. Moreover, our prognostic model exhibited excellent predictive performance in training and external validation cohorts. Finally, a 10-gene classification model for these subtypes identified TTK as a promising therapeutic target with robust classification capabilities. This comprehensive study provides a novel perspective on HCC stratification, offering crucial insights for a deeper understanding of its pathogenesis and the development of promising treatment strategies.

Formerly degenerate seventh zinc finger domain from transcription factor ZNF711 rehabilitated by experimental NMR structure

Domain Z7 of nuclear transcription factor ZNF711 has the consensus last metal-ligand H23 found in odd-numbered zinc-fingers of this protein replaced by a phenylalanine. Ever since the discovery of ZNF711 it has been thought that Z7 is probably non-functional because of the H23F substitution. The presence of H26 three positions downstream prompted us to examine if this histidine could substitute as the last metal ligand. The Z7 domain adopts a stable tertiary structure upon metal binding. The NMR structure of Zn2+-bound Z7 shows the classical ββα-fold of CCHH zinc fingers. Mutagenesis and pH titration experiments indicate that H26 is not involved in metal binding and that Z7 has a tridentate metal-binding site comprised of only residues C3, C6, and H19. By contrast, an F23H mutation that introduces a histidine in the consensus position forms a tetradentate ligand. The structure of the WT Z7 is stable causing restricted ring-flipping of phenyalanines 10 and 23. Dynamics are increased with either the H26A or F23H substitutions and aromatic ring rotation is no longer hindered in the two mutants. The mutations have only small effects on the Kd values for Zn2+ and Co2+ and retain the high thermal stability of the WT domain above 80 °C. Like two previously reported designed zinc fingers with the last ligand replaced by water, the WT Z7 domain is catalytically active, hydrolyzing 4-nitophenyl acetate. We discuss the implications of naturally occurring tridentate zinc fingers for cancer mutations and drug targeting of notoriously undruggable transcription factors. Our findings that Z7 can fold with only a subset of three metal ligands suggests the recent view that most everything about protein structure can be predicted through homology modeling might be premature for at least the resilient and versatile zinc-finger motif.

A comprehensive meta-analysis of transcriptome data to identify signature genes associated with pancreatic ductal adenocarcinoma

PURPOSE

Pancreatic ductal adenocarcinoma (PDAC) has a five-year survival rate of less than 5%. Absence of symptoms at primary tumor stages, as well as high aggressiveness of the tumor can lead to high mortality in cancer patients. Most patients are recognized at the advanced or metastatic stage without surgical symptom, because of the lack of reliable early diagnostic biomarkers. The objective of this work was to identify potential cancer biomarkers by integrating transcriptome data.

METHODS

Several transcriptomic datasets comprising of 11 microarrays were retrieved from the GEO database. After pre-processing, a meta-analysis was applied to identify differentially expressed genes (DEGs) between tumor and nontumor samples for datasets. Next, co-expression analysis, functional enrichment and survival analyses were used to determine the functional properties of DEGs and identify potential prognostic biomarkers. In addition, some regulatory factors involved in PDAC including transcription factors (TFs), protein kinases (PKs), and miRNAs were identified.

RESULTS

After applying meta-analysis, 1074 DEGs including 539 down- and 535 up-regulated genes were identified. Pathway enrichment analyzes using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed that DEGs were significantly enriched in the HIF-1 signaling pathway and focal adhesion. The results also showed that some of the DEGs were assigned to TFs that belonged to 23 conserved families. Sixty-four PKs were identified among the DEGs that showed the CAMK family was the most abundant group. Moreover, investigation of corresponding upstream regions of DEGs identified 11 conserved sequence motifs. Furthermore, weighted gene co-expression network analysis (WGCNA) identified 8 modules, more of them were significantly enriched in Ras signaling, p53 signaling, MAPK signaling pathways. In addition, several hubs in modules were identified, including EMP1, EVL, ELP5, DEF8, MTERF4, GLUP1, CAPN1, IGF1R, HSD17B14, TOM1L2 and RAB11FIP3. According to survival analysis, it was identified that the expression levels of two genes, EMP1 and RAB11FIP3 are related to prognosis.

CONCLUSION

We identified several genes critical for PDAC based on meta-analysis and system biology approach. These genes may serve as potential targets for the treatment and prognosis of PDAC.

Comprehensive characterization of coding and non-coding single nucleotide polymorphisms of the Myoneurin (MYNN) gene using molecular dynamics simulation and docking approaches

Genome-wide association studies (GWAS) identified a coding single nucleotide polymorphism, MYNN rs10936599, at chromosome 3q. MYNN gene encodes myoneurin protein, which has been associated with several cancer pathogenesis and disease development processes. However, there needed to be a more detailed characterization of this polymorphism's (and other coding and non-coding polymorphisms) structural, functional, and molecular impact. The current study addressed this gap and analyzed different properties of rs10936599 and non-coding SNPs of MYNN via a thorough computational method. The variant, rs10936599, was predicted functionally deleterious by nine functionality prediction approaches, like SIFT, PolyPhen-2, and REVEL, etc. Following that, structural modifications were estimated through the HOPE server and Mutation3D. Moreover, the mutation was found in a conserved and active residue, according to ConSurf and CPORT. Further, the secondary structures were predicted, followed by tertiary structures, and there was a significant deviation between the native and variant models. Similarly, molecular simulation also showed considerable differences in the dynamic pattern of the wildtype and mutant structures. Molecular docking revealed that the variant binds with better docking scores with ligand NOTCH2. In addition to that, non-coding SNPs located at the MYNN locus were retrieved from the ENSEMBL database. These were found to disrupt the transcription factor binding regulatory regions; nonetheless, only two affect miRNA target sites. Again, eight non-coding variants were detected in the testes with normalized expression, whereas HaploReg v4.1 unveiled annotations for non-coding variants. In summary, in silico comprehensive characterization of coding and non-coding single nucleotide polymorphisms of MYNN gene will assist researchers to work on MYNN gene and establish their association with certain types of cancers.

Z-flipon variants reveal the many roles of Z-DNA and Z-RNA in health and disease

Identifying roles for Z-DNA remains challenging given their dynamic nature. Here, we perform genome-wide interrogation with the DNABERT transformer algorithm trained on experimentally identified Z-DNA forming sequences (Z-flipons). The algorithm yields large performance enhancements (F1 = 0.83) over existing approaches and implements computational mutagenesis to assess the effects of base substitution on Z-DNA formation. We show Z-flipons are enriched in promoters and telomeres, overlapping quantitative trait loci for RNA expression, RNA editing, splicing, and disease-associated variants. We cross-validate across a number of orthogonal databases and define BZ junction motifs. Surprisingly, many effects we delineate are likely mediated through Z-RNA formation. A shared Z-RNA motif is identified in SCARF2, SMAD1, and CACNA1 transcripts, whereas other motifs are present in noncoding RNAs. We provide evidence for a Z-RNA fold that promotes adaptive immunity through alternative splicing of KRAB domain zinc finger proteins. An analysis of OMIM and presumptive gnomAD loss-of-function datasets reveals an overlap of Z-flipons with disease-causing variants in 8.6% and 2.9% of Mendelian disease genes, respectively, greatly extending the range of phenotypes mapped to Z-flipons.

The Roles of Zinc Finger Proteins in Colorectal Cancer

Despite colorectal cancer remaining a leading worldwide cause of cancer-related death, there remains a paucity of effective treatments for advanced disease. The molecular mechanisms underlying the development of colorectal cancer include altered cell signaling and cell cycle regulation that may result from epigenetic modifications of gene expression and function. Acting as important transcriptional regulators of normal biological processes, zinc finger proteins also play key roles in regulating the cellular mechanisms underlying colorectal neoplasia. These actions impact cell differentiation and proliferation, epithelial-mesenchymal transition, apoptosis, homeostasis, senescence, and maintenance of stemness. With the goal of highlighting promising points of therapeutic intervention, we review the oncogenic and tumor suppressor roles of zinc finger proteins with respect to colorectal cancer tumorigenesis and progression.

Characterization of Phosphorylation-Dependent Antibody Binding to Cancer-Mutated Linkers of C2H2 Zinc Finger Proteins by Intact Transition Epitope Mapping-Thermodynamic Weak-Force Order Analysis

With Intact Transition Epitope Mapping-Thermodynamic Weak-force Order (ITEM-TWO) analysis in combination with molecular modeling, the phosphorylation-dependent molecular recognition motif of the anti-HpTGEKP antibody has been investigated with binary and ternary component mixtures consisting of antibody and (phospho-) peptides. Amino acid sequences have been selected to match either the antibody's recognition motif or the cancer-related zinc finger protein mutations and phosphorylations of the respective amino acid residues. Upon electrospraying of all the components of the mixtures, that is, hexapeptides, antibody, and intact immune complexes, the produced ions were subjected to mass spectrometric mass filtering. The antibody ions as well as the immune complex ions traversed into the mass spectrometer's collision chamber, whereas paths of unbound peptide ions were blocked prior to entering the collision cell. After dissociation of the multiply charged immune complexes in the gas phase, the complex-released peptide ions were recorded after having traversed the second mass filter. Complex-released peptides were unambiguously identified by their masses using mass analysis with isotope resolution. From the results of our studies with seven (phospho-) peptides with distinct amino acid sequences, which resembled either the antibody's binding motif or mutations, we conclude the following: (i) A negatively charged phospho group, located near the peptide's N-terminus is mandatory for antibody binding when placed on the peptide surface at a precise distance to the C-terminally located positively charged ε-amino group of a lysinyl residue. (ii) A bulky amino acid residue, such as the tyrosinyl residue at the N-terminal position of the (phospho-) threoninyl residue, abolishes antibody binding. (iii) Two closely spaced phospho groups negatively interfere with the surface polarity pattern and abolish antibody binding as well. (iv) Non-phosphorylated peptides are not binding partners of the anti-HpTGEKP antibody.

Epitope Fine Mapping by Mass Spectrometry: Investigations of Immune Complexes Consisting of Monoclonal Anti-HpTGEKP Antibody and Zinc Finger Protein Linker Phospho-Hexapeptides

Accurate formation of antibody-antigen complexes has been relied on in both, multitudes of scientific projects and ample therapeutic and diagnostic applications. Mass spectrometrically determined dissociation behavior of immune complexes with the anti-HpTGEKP antibody revealed that the ten most frequently occurring phospho-hexapeptide linker sequences from C2H2 zinc finger proteins could be divided into two classes: orthodox binders, where strong noncovalent interactions developed as anticipated, and unorthodox binders with deviating structures and weaker binding. Phosphorylation of threonine was compulsory for antibody binding in an orthodox manner. Gas phase dissociation energy determinations of seven C2H2 zinc finger protein linker phospho-hexapeptides with orthodox binding properties revealed a bipolar binding motif of the antibody paratope. Epitope peptides, which in addition to the negatively charged phospho-threonine residue were C-terminally flanked by positively charged residues provided stronger binding, i. e. dissociation was endothermic, than peptides with acidic amino acid residues at these positions, for which dissociation was exothermic.

Familial Predisposition to Leiomyomata: Searching for Protective Genetic Factors

In order to determine genetic loci associated with decreasing risk of uterine leiomyomata (UL), a genome-wide association study (GWAS) was performed. We analyzed a group of patients with a family history of UL and a control group consisting of patients without uterine fibroids and a family predisposition to this pathology. Six significant single nucleotide polymorphisms were selected for PCR-genotyping of a large data set of patients with UL. All investigated loci (rs3020434, rs11742635, rs124577644, rs12637801, rs2861221, and rs17677069) demonstrated the lower frequency of minor alleles within a group of women with UL, especially in a subgroup consisting of patients with UL and a familial history of leiomyomata. We also found that the minor allele frequencies of these SNPs in our control group were higher than those across the Caucasian population in all. Based on the obtained data, an evaluation of the common risk of UL was performed. Further work will pave the way to create a specific SNP-panel and allow us to estimate a genotype-based leiomyoma incidence risk. Subsequent studies of genetic variability in a group of patients with a familial predisposition to UL will allow us to make the prediction of the development and course of the disease more individualized, as well as to give our patients personalized recommendations about individual reproductive strategies.

Zinc finger protein 852 is essential for the proliferation, drug sensitivity, and self-renewal of gastric cancer cells

Exploring cellular and molecular mechanisms responsible for gastric cancer growth, survival, self-renewal, and metastasis helps develop efficacious therapeutic strategies. In this study, the expression and function of Zinc finger protein 852 (ZNF852) in human gastric cancer cell lines were characterized. ZNF852 was up-regulated in gastric cancer cell lines relative to normal gastric epithelial cell line GES-1. When the ZNF852 gene was ablated in gastric cancer cell line MGC-803 using the CRISPR/Cas9-encoding lentivirus, the proliferation of MGC-803 was suppressed. ZNF852 deficiency also resulted in the inhibition of MGC-803 sphere formation, along with decreases in SRY-box 2 (SOX2), Octamer-binding transcription factor 4 (OCT4), and Nanog homeobox (NANOG), suggesting that ZNF852 sustains self-renewal of MGC-803 cells. Furthermore, ZNF852 deficiency increased oxaliplatin-induced MGC-803 cell death, implying the role of ZNF852 in drug sensitivity. Subcutaneous infusion of MGC-803 cells into nude mice illustrated the same effects of ZNF852 on the proliferation and self-renewal of gastric cancer cells. Similar effects of ANF852 were also seen in gastric cancer cell line SNU-1. Interestingly, ZNF852 deficiency caused down-regulation of epidermal growth factor receptor (EGFR) on gastric cancer cells. In summary, this study uncovers the positive regulatory role of ZNF852 in gastric cancer growth and maintenance. ZNF852 could be a potential therapeutic target for inhibiting gastric cancer initiation or progression. This article is protected by copyright. All rights reserved.

Intratumor heterogeneity and clonal evolution revealed in castration-resistant prostate cancer by longitudinal genomic analysis

Intratumor heterogeneity is a key driver for local relapse and treatment failure. Thus, using multifocal prostate cancer as a model to investigate tumor inter-clonal relationships and tumor evolution could aid in our understanding of drug resistance. Previous studies discovered genomic alterations by comparing hormone-sensitive prostate cancer (HSPC) with castration-resistant prostate cancer (CRPC) in large cohorts. However, most studies did not sequentially sample tumors from the same patient. In our study, we performed whole-exome sequencing (WES) on 14 specimens from five locally relapsed patients before and after androgen-deprivation therapy. We described the landscape of genomic alterations before and after treatment and identified critical driver events that could have contributed to the evolution of CRPC. In addition to confirming known cancer genes such as TP53 and CDK12, we also identified new candidate genes that may play a role in the progression of prostate cancer, including MYO15A, CHD6 and LZTR1. At copy number alteration (CNA) level, gain of 8q24.13-8q24.3 was observed in 60% of patients and was the most commonly altered locus in both HSPC and CRPC tumors. Finally, utilizing phylogenetic reconstruction, we explored the clonal progression pattern from HSPC to CRPC in each patient. Our findings highlight the complex and heterogeneous mechanisms underlying the development of drug resistance, and underscore the potential value of monitoring tumor clonal architectures during disease progression in a clinical setting.

Histidine phosphorylation in metalloprotein binding sites

Post-translational modifications (PTMs) are invaluable regulatory tools for the control of catalytic functionality, protein-protein interactions, and signaling pathways. Historically, the study of phosphorylation as a PTM has been focused on serine, threonine, and tyrosine residues. In contrast, the significance of mammalian histidine phosphorylation remains largely unexplored. This gap in knowledge regarding the molecular basis for histidine phosphorylation as a regulatory agent exists in part because of the relative instability of phosphorylated histidine as compared with phosphorylated serine, threonine and tyrosine. However, the unique metal binding abilities of histidine make it one of the most common metal coordinating ligands in nature, and it is interesting to consider how phosphorylation would change the metal coordinating ability of histidine, and consequently, the properties of the phosphorylated metalloprotein. In this review, we examine eleven metalloproteins that have been shown to undergo reversible histidine phosphorylation at or near their metal binding sites. These proteins are described with respect to their biological activity and structure, with a particular emphasis on how phosphohistidine may tune the primary coordination sphere and protein conformation. Furthermore, several common methods, challenges, and limitations of studying sensitive, high affinity metalloproteins are discussed.

The Identification of Zinc-Finger Protein 433 as a Possible Prognostic Biomarker for Clear-Cell Renal Cell Carcinoma

Clear-cell renal cell carcinoma (ccRCC) is the most common and aggressive form of all urological cancers, with poor prognosis and high mortality. At late stages, ccRCC is known to be mainly resistant to chemotherapy and radiotherapy. Therefore, it is urgent and necessary to identify biomarkers that can facilitate the early detection of ccRCC in patients. In this study, the levels of transcripts of ccRCC from The Cancer Genome Atlas (TCGA) dataset were used to identify prognostic biomarkers in this disease. Analyzing the data obtained indicated that the KRAB-ZNF protein is significantly suppressed in clear-cell carcinomas. Furthermore, ZNF433 is differentially expressed in ccRCC in a stage- and histological-grade-specific manner. In addition, ZNF433 expression was correlated with metastasis, with greater node involvement associated with lower ZNF433 expression (p < 0.01) and with a more unsatisfactory overall survival outcome (HR, 0.45; 95% CI, 0.33–0.6; p = 8.5 × 10⁻⁸). Since ccRCC is characterized by mutations in proteins that alter epigenetic modifications and /or chromatin remodeling, we examined the expression of ZNF433 transcripts in ccRCC with wildtype and mutated forms of BAP1, KDMC5, MTOR, PBRM1, SETD2, and VHL. Analysis revealed that ZNF433 expression was significantly reduced in ccRCC with mutations in the BAP1, SETD2, and KDM5C genes (p < 0.05). In addition, the ZNF433 promoter region was highly methylated, and hypermethylation was significantly associated with mRNA suppression (p < 2.2 × 10⁻¹⁶). In silico analysis of potential ZNF target genes found that the largest group of target genes are involved in cellular metabolic processes, which incidentally are particularly impaired in ccRCC. It was concluded from this study that gene expression of ZNF433 is associated with cancer progression and poorer prognosis, and that ZNF433 behaves in a manner that suggests that it is a prognostic marker and a possible tumor-suppressor gene in clear-cell renal cell carcinoma.

Signatures of Discriminative Copy Number Aberrations in 31 Cancer Subtypes

Copy number aberrations (CNA) are one of the most important classes of genomic mutations related to oncogenetic effects. In the past three decades, a vast amount of CNA data has been generated by molecular-cytogenetic and genome sequencing based methods. While this data has been instrumental in the identification of cancer-related genes and promoted research into the relation between CNA and histo-pathologically defined cancer types, the heterogeneity of source data and derived CNV profiles pose great challenges for data integration and comparative analysis. Furthermore, a majority of existing studies have been focused on the association of CNA to pre-selected "driver" genes with limited application to rare drivers and other genomic elements. In this study, we developed a bioinformatics pipeline to integrate a collection of 44,988 high-quality CNA profiles of high diversity. Using a hybrid model of neural networks and attention algorithm, we generated the CNA signatures of 31 cancer subtypes, depicting the uniqueness of their respective CNA landscapes. Finally, we constructed a multi-label classifier to identify the cancer type and the organ of origin from copy number profiling data. The investigation of the signatures suggested common patterns, not only of physiologically related cancer types but also of clinico-pathologically distant cancer types such as different cancers originating from the neural crest. Further experiments of classification models confirmed the effectiveness of the signatures in distinguishing different cancer types and demonstrated their potential in tumor classification.

Recurrent ZNF83-E293V Mutation Promotes Bladder Cancer Progression through the NF-κB Pathway via Transcriptional Dysregulation of S100A8

Urothelial carcinoma (UC) is the predominant form of bladder cancer. Significant molecular heterogeneity caused by diverse molecular alterations brings about large variations in the response to treatment in UC. An improved understanding of the genetic mechanisms underlying the development and progression of UC is essential. Through deep analysis of next-generation sequencing data of 99 UC patients, we found that 18% of cases had recurrent somatic mutations in zinc finger protein gene zinc finger protein 83 (ZNF83). ZNF83 mutations were correlated with poor prognosis of UC. We also found a hotspot mutation, p.E293V, in the evolutionarily well-conserved region of ZNF83. ZNF83-E293V increased tumor growth and reduced the apoptosis of UC cells compared to wild-type ZNF83 both in vitro and in mice xenografted tumors. ZNF83-E293V activated nuclear factor κB (NF-κB) more potently than did the wild-type protein owing to its decreased transcriptional repression for S100A8. The NF-κB inhibitors could pharmacologically block the tumor growth in mice engrafted with ZNF83-E293V-transfected UC cells. These findings provide a mechanistic insight and a potential therapeutic strategy for UC, which established a foundation for using the ZNF83-E293V mutation as a predictive biomarker of therapeutic response from NF-κB inhibitors.

The Role of Zinc and Copper in Gynecological Malignancies

Zinc (Zn) and copper (Cu) are essential microelements, which take part in cellular metabolism, feature in enzymatic systems, and regulate enzyme activity. Homeostasis of these micronutrients is tightly regulated by multiple compensatory mechanisms that balance their concentrations including transporters, importers, and metallothioneins. An altered intake of only one of these trace elements may cause an imbalance in their levels and result in their competition for absorption. Relatively low levels of zinc and increased levels of copper may result in an increased level of oxidative stress and impair the antioxidant properties of multiple enzymes. Altered levels of trace elements were discovered in various pathologies including immunological, degenerative, and inflammatory diseases. Moreover, due to the role of Zn and Cu in oxidative stress and chronic inflammation, they were found to influence cancerogenesis. We review the roles of zinc and copper and their mechanisms in tumor growth, metastasis potential, microenvironment remodeling, and drug resistance. We highlight their role as potential biomarkers for cancer diagnosis, treatment, and prognosis, concentrating on their impact on gynecological malignancies.

Cancer and pH Dynamics: Transcriptional Regulation, Proteostasis, and the Need for New Molecular Tools

An emerging hallmark of cancer cells is dysregulated pH dynamics. Recent work has suggested that dysregulated intracellular pH (pHi) dynamics enable diverse cancer cellular behaviors at the population level, including cell proliferation, cell migration and metastasis, evasion of apoptosis, and metabolic adaptation. However, the molecular mechanisms driving pH-dependent cancer-associated cell behaviors are largely unknown. In this review article, we explore recent literature suggesting pHi dynamics may play a causative role in regulating or reinforcing tumorigenic transcriptional and proteostatic changes at the molecular level, and discuss outcomes on tumorigenesis and tumor heterogeneity. Most of the data we discuss are population-level analyses; lack of single-cell data is driven by a lack of tools to experimentally change pHi with spatiotemporal control. Data is also sparse on how pHi dynamics play out in complex in vivo microenvironments. To address this need, at the end of this review, we cover recent advances for live-cell pHi measurement at single-cell resolution. We also discuss the essential role for tool development in revealing mechanisms by which pHi dynamics drive tumor initiation, progression, and metastasis.

PertInInt: An Integrative, Analytical Approach to Rapidly Uncover Cancer Driver Genes with Perturbed Interactions and Functionalities

A major challenge in cancer genomics is to identify genes with functional roles in cancer and uncover their mechanisms of action. We introduce an integrative framework that identifies cancer-relevant genes by pinpointing those whose interaction or other functional sites are enriched in somatic mutations across tumors. We derive analytical calculations that enable us to avoid time-prohibitive permutation-based significance tests, making it computationally feasible to simultaneously consider multiple measures of protein site functionality. Our accompanying software, PertInInt, combines knowledge about sites participating in interactions with DNA, RNA, peptides, ions, or small molecules with domain, evolutionary conservation, and gene-level mutation data. When applied to 10,037 tumor samples, PertInInt uncovers both known and newly predicted cancer genes, while additionally revealing what types of interactions or other functionalities are disrupted. PertInInt’s analysis demonstrates that somatic mutations are frequently enriched in interaction sites and domains and implicates interaction perturbation as a pervasive cancer-driving event.

A fast, analytical framework called PertInInt enables efficient integration of multiple measures of protein site functionality—including interaction, domain, and evolutionary conservation—with gene-level mutation data in order to rapidly detect cancer driver genes along with their disrupted functionalities.

Modeling and analysis of site-specific mutations in cancer identifies known plus putative novel hotspots and bias due to contextual sequences

In cancer, recurrently mutated sites in DNA and proteins, called hotspots, are thought to be raised by positive selection and therefore important due to its potential functional impact. Although recent evidence for APOBEC enzymatic activity have shown that specific types of sequences are likely to be false, the identification of putative hotspots is important to confirm either its functional role or its mechanistic bias. In this work, an algorithm and a statistical model is presented to detect hotspots. The model consists of a beta-binomial component plus fixed effects that efficiently fits the distribution of mutated sites. The algorithm employs an optimal stepwise approach to find the model parameters. Simulations show that the proposed algorithmic model is highly accurate for common hotspots. The approach has been applied to TCGA mutational data from 33 cancer types. The results show that well-known cancer hotspots are easily detected. Besides, novel hotspots are also detected. An analysis of the sequence context of detected hotspots show a preference for TCG sites that may be related to APOBEC or other unknown mechanistic biases. The detected hotspots are available online in http://bioinformatica.mty.itesm.mx/HotSpotsAnnotations.

The geometric influence on the Cys2His2 zinc finger domain and functional plasticity

The Cys₂His₂ zinc finger is the most common DNA-binding domain expanding in metazoans since the fungi human split. A proposed catalyst for this expansion is an arms race to silence transposable elements yet it remains poorly understood how this domain is able to evolve the required specificities. Likewise, models of its DNA binding specificity remain error prone due to a lack of understanding of how adjacent fingers influence each other's binding specificity. Here, we use a synthetic approach to exhaustively investigate binding geometry, one of the dominant influences on adjacent finger function. By screening over 28 billion protein–DNA interactions in various geometric contexts we find the plasticity of the most common natural geometry enables more functional amino acid combinations across all targets. Further, residues that define this geometry are enriched in genomes where zinc fingers are prevalent and specificity transitions would be limited in alternative geometries. Finally, these results demonstrate an exhaustive synthetic screen can produce an accurate model of domain function while providing mechanistic insight that may have assisted in the domains expansion.

Evaluation of the Anti-Tumor Activity of the Humanized Monoclonal Antibody NEO-201 in Preclinical Models of Ovarian Cancer

Purpose: Despite high initial response rates with cytoreductive surgery, conventional chemotherapy and the incorporation of biologic agents, ovarian cancer patients often relapse and die from their disease. New approaches are needed to improve patient outcomes. This study was designed to evaluate the antitumor activity of NEO-201 monoclonal antibody (mAb) in preclinical models of ovarian cancer where the NEO-201 target is highly expressed.

Experimental Design: Functional analysis of NEO-201 against tumor cell lines was performed by antibody-dependent cellular cytotoxicity (ADCC) assays. Binding of NEO-201 to tumor tissues and cell lines were determined by immunohistochemistry (IHC) and flow cytometry, respectively. Further characterization of the antigen recognized by NEO-201 was performed by mass spectrometry. Ovarian cancer models were used to evaluate the anti-tumor activity of NEO-201 in vivo. NEO-201 at a concentration of 250 g/mouse was injected intraperitoneally (IP) on days 1, 4, and 8. Human PBMCs were injected IP simultaneously as effector cells.

Results: Both IHC and flow cytometry revealed that NEO-201 binds prominently to the colon, pancreatic, and mucinous ovarian cancer tissues and cell lines. Immunoprecipitation of the antigen recognized by NEO-201 was performed in human ovarian, colon, and pancreatic cancer cell lines. From these screening, carcinoembryonic antigen-related cell adhesion molecule 5 (CEACAM5) and CEACAM6 were identified as the most likely targets of NEO-201. Our results confirmed that NEO-201 binds different types of cancers; the binding is highly selective for the tumor cells without cross reactivity with the surrounding healthy tissue. Functional analysis revealed that NEO-201 mediates ADCC killing against human ovarian and colorectal carcinoma cell lines in vitro. In addition, NEO-201 inhibited tumor growth in the presence of activated human PBMCs in orthotopic mouse models of both primary and metastatic ovarian cancer. Importantly, NEO-201 prolonged survival of tumor-bearing mice.

Conclusions: These data suggested that NEO-201 has an antitumor activity against tumor cells expressing its antigen. Targeting an antigen expressed in tumors, but not in normal tissues, allows patient selection for optimal treatment. These findings strongly indicate that NEO-201 warrants clinical testing as both a novel therapeutic and diagnostic agent for treatment of ovarian carcinomas. A first in human clinical trial evaluating NEO-201 in adults with chemo-resistant solid tumors is ongoing at the NIH clinical Center.

Arsenite Exposure Displaces Zinc from ZRANB2 Leading to Altered Splicing

Exposure to arsenic, a class I carcinogen, affects 200 million people globally. Skin is the major target organ, but the molecular etiology of arsenic-induced skin carcinogenesis remains unclear. Arsenite (As³⁺)-induced disruption of alternative splicing could be involved, but the mechanism is unknown. Zinc finger proteins play key roles in alternative splicing. As³⁺ can displace zinc (Zn²⁺) from C3H1 and C4 zinc finger motifs (zfm's), affecting protein function. ZRANB2, an alternative splicing regulator with two C4 zfm's integral to its structure and splicing function, was chosen as a candidate for this study. We hypothesized that As³⁺ could displace Zn²⁺ from ZRANB2, altering its structure, expression, and splicing function. As³⁺/Zn²⁺ binding and mutual displacement experiments were performed with synthetic apo-peptides corresponding to each ZRANB2 zfm, employing a combination of intrinsic fluorescence, ultraviolet spectrophotometry, zinc colorimetric assay, and liquid chromatography-tandem mass spectrometry. ZRANB2 expression in HaCaT cells acutely exposed to As³⁺ (0 or 5 μM, 0-72 h; or 0-5 μM, 6 h) was examined by RT-qPCR and immunoblotting. ZRANB2-dependent splicing of TRA2B mRNA, a known ZRANB2 target, was monitored by reverse transcription-polymerase chain reaction. As³⁺ bound to, as well as displaced Zn²⁺ from, each zfm. Also, Zn²⁺ displaced As³⁺ from As³⁺-bound zfm's acutely, albeit transiently. As³⁺ exposure induced ZRANB2 protein expression between 3 and 24 h and at all exposures tested but not ZRANB2 mRNA expression. ZRANB2-directed TRA2B splicing was impaired between 3 and 24 h post-exposure. Furthermore, ZRANB2 splicing function was also compromised at all As³⁺ exposures, starting at 100 nm. We conclude that As³⁺ exposure displaces Zn²⁺ from ZRANB2 zfm's, changing its structure and compromising splicing of its targets, and increases ZRANB2 protein expression as a homeostatic response both at environmental/toxicological exposures and therapeutically relevant doses.

Epigenetic mechanisms are behind the regulation of the key genes associated with the osteoblastic differentiation of the mesenchymal stem cells: The role of zoledronic acid on tuning the epigenetic changes

Mesenchymal stem cells (MSCs) are multipotent stem cells and show distinct features such as capability for self-renewal and differentiation into several lineages of cells including osteoblasts, chondrocytes, and adipocytes. In this study, the methylation status of the promoter region of zinc finger and BTB domain containing 16 (ZBTB16), twist-related protein 1(Twist1), de novo DNA methyltransferases 3A (DNMT3A), SRY-box 9 (Sox9), osteocalcin (OCN), and peroxisome proliferator-activated receptor γ2 (PPARγ2) genes and their messenger RNA (mRNA) expression levels were evaluated during the osteoblastic differentiation of MSCs (ODMSCs). We planned two experimental groups including zoledronic acid (ZA)-treated and nontreated cells (negative control) which both were differentiated into the osteoblasts. Methylation level of DNA in the promoter regions was assayed by methylation-specific-quantitative polymerase chain reaction (MS-qPCR), and mRNA levels of the target inhibitory/stimulatory genes during osteoblastic differentiation of MSCs were measured using real-time PCR. During the experimental induction of ODMSCs, the mRNA expression of the OCN gene was upregulated and methylation level of its promoter region was decreased. Moreover, Sox9 and PPARγ2 mRNA levels were attenuated and their promoter regions methylation levels were significantly augmented. However, the mRNA expression of the DNMT3A was not affected during the ODMSCs though its methylation rate was increased. In addition, ZA could enhance the expression of the ZBTB16 and decrease its promoter regions methylation and on the opposite side, it diminished mRNA expression of Sox9, Twist1, and PPARγ2 genes and increased their methylation rates. Intriguingly, ZA did not show a significant impact on gene expression and methylation levels the OCN and DNMT3A. We found that methylation of the promoter regions of Sox9, OCN, and PPARγ2 genes might be one of the main mechanisms adjusting the genes expression during the ODMSCs. Furthermore, we noticed that ZA can accelerate the MSCs differentiation to the osteoblast cells via two regulatory processes; suppression of osteoblastic differentiation inhibitor genes including Sox9, Twist1, and PPARγ2, and through promotion of the ZBTB16 expression.

Precision medicine review: rare driver mutations and their biophysical classification

How can biophysical principles help precision medicine identify rare driver mutations? A major tenet of pragmatic approaches to precision oncology and pharmacology is that driver mutations are very frequent. However, frequency is a statistical attribute, not a mechanistic one. Rare mutations can also act through the same mechanism, and as we discuss below, “latent driver” mutations may also follow the same route, with “helper” mutations. Here, we review how biophysics provides mechanistic guidelines that extend precision medicine. We outline principles and strategies, especially focusing on mutations that drive cancer. Biophysics has contributed profoundly to deciphering biological processes. However, driven by data science, precision medicine has skirted some of its major tenets. Data science embodies genomics, tissue- and cell-specific expression levels, making it capable of defining genome- and systems-wide molecular disease signatures. It classifies cancer driver genes/mutations and affected pathways, and its associated protein structural data guide drug discovery. Biophysics complements data science. It considers structures and their heterogeneous ensembles, explains how mutational variants can signal through distinct pathways, and how allo-network drugs can be harnessed. Biophysics clarifies how one mutation—frequent or rare—can affect multiple phenotypic traits by populating conformations that favor interactions with other network modules. It also suggests how to identify such mutations and their signaling consequences. Biophysics offers principles and strategies that can help precision medicine push the boundaries to transform our insight into biological processes and the practice of personalized medicine. By contrast, “phenotypic drug discovery,” which capitalizes on physiological cellular conditions and first-in-class drug discovery, may not capture the proper molecular variant. This is because variants of the same protein can express more than one phenotype, and a phenotype can be encoded by several variants.