Binding Motif Research Articles

A network of transcription factors in complex with a regulating cell cycle cyclin orchestrates fungal oxidative stress responses

BackgroundResponse to oxidative stress is universal in almost all organisms and the mitochondrial membrane protein, BbOhmm, negatively affects oxidative stress responses and virulence in the insect fungal pathogen, Beauveria bassiana. Nothing further, however, is known concerning how BbOhmm and this phenomenon is regulated.ResultsThree oxidative stress response regulating Zn2Cys6 transcription factors (BbOsrR1, 2, and 3) were identified and verified via chromatin immunoprecipitation (ChIP)-qPCR analysis as binding to the BbOhmm promoter region, with BbOsrR2 showing the strongest binding. Targeted gene knockout of BbOsrR1 or BbOsrR3 led to decreased BbOhmm expression and consequently increased tolerances to free radical generating compounds (H2O2 and menadione), whereas the ΔBbOsrR2 strain showed increased BbOhmm expression with concomitant decreased tolerances to these compounds. RNA and ChIP sequencing analysis revealed that BbOsrR1 directly regulated a wide range of antioxidation and transcription-associated genes, negatively affecting the expression of the BbClp1 cyclin and BbOsrR2. BbClp1 was shown to localize to the cell nucleus and negatively mediate oxidative stress responses. BbOsrR2 and BbOsrR3 were shown to feed into the Fus3-MAPK pathway in addition to regulating antioxidation and detoxification genes. Binding motifs for the three transcription factors were found to partially overlap in the promoter region of BbOhmm and other target genes. Whereas BbOsrR1 appeared to function independently, co-immunoprecipitation revealed complex formation between BbClp1, BbOsrR2, and BbOsrR3, with BbClp1 partially regulating BbOsrR2 phosphorylation.ConclusionsThese findings reveal a regulatory network mediated by BbOsrR1 and the formation of a BbClp1-BbOsrR2-BbOsrR3 complex that orchestrates fungal oxidative stress responses.

Nanostructure and dynamics of N-truncated copper amyloid-β peptides from advanced X-ray absorption fine structure.

An X-ray absorption spectroscopy (XAS) electrochemical cell was used to collect high-quality XAS measurements of N-truncated Cu:amyloid-β (Cu:Aβ) samples under near-physiological conditions. N-truncated Cu:Aβ peptide complexes contribute to oxidative stress and neurotoxicity in Alzheimer's patients' brains. However, the redox properties of copper in different Aβ peptide sequences are inconsistent. Therefore, the geometry of binding sites for the copper binding in Aβ4-8/12/16 was determined using novel advanced extended X-ray absorption fine structure (EXAFS) analysis. This enables these peptides to perform redox cycles in a manner that might produce toxicity in human brains. Fluorescence XAS measurements were corrected for systematic errors including defective-pixel data, monochromator glitches and dispersion of pixel spectra. Experimental uncertainties at each data point were measured explicitly from the point-wise variance of corrected pixel measurements. The copper-binding environments of Aβ4-8/12/16 were precisely determined by fitting XAS measurements with propagated experimental uncertainties, advanced analysis and hypothesis testing, providing a mechanism to pursue many similarly complex questions in bioscience. The low-temperature XAS measurements here determine that CuII is bound to the first amino acids in the high-affinity amino-terminal copper and nickel (ATCUN) binding motif with an oxygen in a tetragonal pyramid geometry in the Aβ4-8/12/16 peptides. Room-temperature XAS electrochemical-cell measurements observe metal reduction in the Aβ4-16 peptide. Robust investigations of XAS provide structural details of CuII binding with a very different bis-His motif and a water oxygen in a quasi-tetrahedral geometry. Oxidized XAS measurements of Aβ4-12/16 imply that both CuII and CuIII are accommodated in an ATCUN-like binding site. Hypotheses for these CuI, CuII and CuIII geometries were proven and disproven using the novel data and statistical analysis including F tests. Structural parameters were determined with an accuracy some tenfold better than literature claims of past work. A new protocol was also developed using EXAFS data analysis for monitoring radiation damage. This gives a template for advanced analysis of complex biosystems.

Construction and analysis of competitive endogenous RNA networks and prognostic models associated with ovarian cancer based on the exoRBase database.

To construct a competitive endogenous RNA (ceRNA) regulatory network in blood exosomes of patients with ovarian cancer (OC) using bioinformatics and explore its pathogenesis. The exoRbase2.0 database was used to download blood exosome gene sequencing data from patients OC and normal controls and the expression profiles of exosomal mRNA, long non-coding RNA (lncRNA), and circular RNA (circRNA) were detected independently using R language for differential expression analysis. TargetScan and miRanda databases were combined for the prediction and differential expression of mRNA-binding microRNAs (miRNA). The miRcode and starBase databases were used to predict miRNAs that bind to differentially expressed lncRNAs and circRNAs repectively. The relevant mRNA, circRNA, lncRNA and their corresponding miRNA prediction data were imported into Cytoscape software for visualization of the ceRNA network. The R language and KEGG Orthology-based Annotation System (KOBAS) were used to execute and illustrate the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Hub genes were identified using The CytoHubba plugin. Thirty-one differentially expressed mRNAs, 17 differentially expressed lncRNAs, and 24 differentially expressed circRNAs were screened. Cytoscape software was used to construct the ceRNA network with nine mRNA nodes, two lncRNA nodes, eight circRNA nodes, and 51 miRNA nodes. Both GO and KEGG were focused on the Spliceosome pathway, indicating that spliceosomes are closely linked with the development of OC, while heterogenous nuclear ribonucleoprotein K and RNA binding motif protein X-linked genes were the top 10 score Hub genes screened by Cytoscape software, including two lncRNAs, four mRNAs, and four circRNAs. In patients with OC, the expression of eukaryotic translation initiation factor 4 gamma 2 (EIF4G2), SERPINE 1 mRNA binding protein 1 (SERBP1), ribosomal protein L15 (RPL15) and human leukocyte antigen complex P5 (HCP5) was significantly higher whereas that of testis expressed transcript, Y-linked 15 and DEAD-box helicase 3 Y-linked genes was lower compared to normal controls Immunocorrelation scores revealed that SERBP1 was significantly and negatively correlated with endothelial cells and CD4+ T cells and positively correlated with natural killer (NK) cells and macrophages, respectively; RPL15 was significantly positively correlated with macrophages and endothelial cells and negatively correlated with CD8+ T cells and uncharacterized cells, respectively. EIF4G2 was significantly and negatively correlated with endothelial cells and CD4+ T cells, and positively correlated with uncharacterized cells, respectively. Based on the survival data and the significant correlation characteristics derived from the multifactorial Cox analysis (P < 0.05), the survival prediction curves demonstrated that the prognostic factors associated with 3-year survival in patients with OC were The prognostic factors associated with survival were Macrophage, RPL15. This study successfully constructs a ceRNA regulatory network in blood exosomes of OV patients, which provides the specific targets for diagnosis and treatment of OC.

Functional Integrity of Radical SAM Enzyme Dph1•Dph2 Requires Non-Canonical Cofactor Motifs with Tandem Cysteines.

The Dph1•Dph2 heterodimer from yeast is a radical SAM (RS) enzyme that generates the 3-amino-3-carboxy-propyl (ACP) precursor for diphthamide, a clinically relevant modification on eukaryotic elongation factor 2 (eEF2). ACP formation requires SAM cleavage and atypical Cys-bound Fe-S clusters in each Dph1 and Dph2 subunit. Intriguingly, the first Cys residue in each motif is found next to another ill-defined cysteine that we show is conserved across eukaryotes. As judged from structural modeling, the orientation of these tandem cysteine motifs (TCMs) suggests a candidate Fe-S cluster ligand role. Hence, we generated, by site-directed DPH1 and DPH2 mutagenesis, Dph1•Dph2 variants with cysteines from each TCM replaced individually or in combination by serines. Assays diagnostic for diphthamide formation in vivo reveal that while single substitutions in the TCM of Dph2 cause mild defects, double mutations almost entirely inactivate the RS enzyme. Based on enhanced Dph1 and Dph2 subunit instability in response to cycloheximide chases, the variants with Cys substitutions in their cofactor motifs are particularly prone to protein degradation. In sum, we identify a fourth functionally cooperative Cys residue within the Fe-S motif of Dph2 and show that the Cys-based cofactor binding motifs in Dph1 and Dph2 are critical for the structural integrity of the dimeric RS enzyme in vivo.

Kinesin Regulation in the Proximal Axon is Essential for Dendrite-selective Transport.

Neurons are polarized and typically extend multiple dendrites and one axon. To maintain polarity, vesicles carrying dendritic proteins are arrested upon entering the axon. To determine whether kinesin regulation is required for terminating anterograde axonal transport, we overexpressed the dendrite-selective kinesin KIF13A. This caused mistargeting of dendrite-selective vesicles to the axon and a loss of dendritic polarity. Polarity was not disrupted if the kinase MARK2/Par1b was coexpressed. MARK2/Par1b is concentrated in the proximal axon, where it maintains dendritic polarity-likely by phosphorylating S1371 of KIF13A, which lies in a canonical 14-3-3 binding motif. We probed for interactions of KIF13A with 14-3-3 isoforms and found that 14-3-3β and 14-3-3ζ bound KIF13A. Disruption of MARK2 or 14-3-3 activity by small molecule inhibitors caused a loss of dendritic polarity. These data show that kinesin regulation is integral for dendrite-selective transport. We propose a new model in which KIF13A that moves dendrite-selective vesicles in the proximal axon is phosphorylated by MARK2. Phosphorylated KIF13A is then recognized by 14-3-3, which causes dissociation of KIF13A from the vesicle and termination of transport. These findings define a new paradigm for the regulation of vesicle transport by localized kinesin tail phosphorylation, to restrict dendrite-selective vesicles from entering the axon.

Type I interferons induce an epigenetically distinct memory B cell subset in chronic viral infection

Memory B cells (MBCs) are key providers of long-lived immunity against infectious disease, yet in chronic viral infection, they do not produce effective protection. How chronic viral infection disrupts MBC development and whether such changes are reversible remain unknown. Through single-cell (sc)ATAC-seq and scRNA-seq during acute versus chronic lymphocytic choriomeningitis viral infection, we identified a memory subset enriched for interferon (IFN)-stimulated genes (ISGs) during chronic infection that was distinct from the T-bet+ subset normally associated with chronic infection. Blockade of IFNAR-1 early in infection transformed the chromatin landscape of chronic MBCs, decreasing accessibility at ISG-inducing transcription factor binding motifs and inducing phenotypic changes in the dominating MBC subset, with a decrease in the ISG subset and an increase in CD11c+CD80+ cells. However, timing was critical, with MBCs resistant to intervention at 4weeks post-infection. Together, our research identifies a key mechanism to instruct MBC identity during viral infection.

Self‐Adaptive Aromatic Cation‐π Driven Dimensional Polymorphism in Supramolecular Polymers for the Photocatalytic Oxidation and Separation of Aromatic/Cyclic Aliphatic Compounds

AbstractThe phenomenon of polymorphism is ubiquitous in nature, the controlled manipulation of which not only increases our ontological understanding of nature but also facilitates the conceptualization and realization of novel functional materials. However, achieving targeted polymorphism in supramolecular assemblies (SAs) remains a formidable challenge, largely because of the constraints inherent in controlling the specific binding motifs of noncovalent interactions. Herein, we propose self‐adaptive aromatic cation‐π binding motifs to construct polymorphic SAs in both the solid and solution states. Using distinct discrete cation‐π‐cation and long‐range cation‐π binding motifs enables control of the self‐assembly directionality of a C2h‐symmetric bifunctional monomer, resulting in the successful formation of both two‐dimensional and three‐dimensional crystalline SAs (2D‐CSA and 3D‐CSA). The differences in the molecular packing of 3D‐CSA compared with that of 2D‐CSA significantly improve the charge separation and carrier mobility, leading to enhanced photocatalytic activity for the aerobic oxidation of thioanisole to methyl phenyl sulfoxide (yield of 99 % vs 57 %). 2D‐CSA, which has a vertical extended structure with favorable stronger interaction with toluene though face‐to‐face cation‐π interactions than methylcyclohexane, shows higher toluene/methylcyclohexane separation efficiency than 3D‐CSA (96.9 % for 2D‐CSA vs 56.3 % for 3D‐CSA).

Abstract LB293: O-GlcNAc transferase drives the proliferative state in adaptive drug tolerance

Abstract Drivers of the drug tolerant proliferative persister (DTPP) state have not been well investigated. Histone H3 lysine-4 trimethylation (H3K4me3), an active histone mark, might enable slow cycling drug tolerant persisters (DTP) to regain proliferative capacity. This study aimed to determine H3K4me3 transcriptionally active sites identifying a key regulator of DTPPs in melanoma. Deploying a model of adaptive cancer drug tolerance H3K4me3 ChIP-seq data of DTPPs guided identification of top transcription factor binding motifs. These suggested involvement of O-linked N-acetylglucosamine transferase (OGT), which was confirmed by metabolomics analysis and biochemical assays. OGT impact on DTPPs and adaptive resistance was explored in vitro and in vivo. H3K4me3 remodeling was widespread in CPG island regions and DNA binding motifs associated with O-GlcNAc marked chromatin. Accordingly, we observed an upregulation of OGT, O-GlcNAc and its binding partner TET1 in chronically treated cancer cells. Inhibition of OGT led to loss of H3K4me3 and downregulation of genes contributing to drug resistance. Genetic ablation of OGT prevented acquired drug resistance in vivo. Our findings uncover a fundamental mechanism of adaptive drug resistance that governs cancer cell reprogramming towards acquired drug resistance, a process that can be exploited to improve response duration and patient outcomes. Citation Format: Dinoop Ravindran Menon, Heinz Hammerlindl, Gregory Gimenez, Elmar Zuegner, Christoph Magens, Michael Eccles, Mayumi Fujita, Helmut Schaider. O-GlcNAc transferase drives the proliferative state in adaptive drug tolerance [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB293.

Self-Adaptive Aromatic Cation-π Driven Dimensional Polymorphism in Supramolecular Polymers for the Photocatalytic Oxidation and Separation of Aromatic/Cyclic Aliphatic Compounds.

The phenomenon of polymorphism is ubiquitous in nature, the controlled manipulation of which not only increases our ontological understanding of nature but also facilitates the conceptualization and realization of novel functional materials. However, achieving targeted polymorphism in supramolecular assemblies (SAs) remains a formidable challenge, largely because of the constraints inherent in controlling the specific binding motifs of noncovalent interactions. Herein, we propose self-adaptive aromatic cation-π binding motifs to construct polymorphic SAs in both the solid and solution states. Using distinct discrete cation-π-cation and long-range cation-π binding motifs enables control of the self-assembly directionality of a C2h-symmetric bifunctional monomer, resulting in the successful formation of both two-dimensional and three-dimensional crystalline SAs (2D-CSA and 3D-CSA). The differences in the molecular packing of 3D-CSA compared with that of 2D-CSA significantly improve the charge separation and carrier mobility, leading to enhanced photocatalytic activity for the aerobic oxidation of thioanisole to methyl phenyl sulfoxide (yield of 99 % vs 57 %). 2D-CSA, which has a vertical extended structure with favorable stronger interaction with toluene though face-to-face cation-π interactions than methylcyclohexane, shows higher toluene/methylcyclohexane separation efficiency than 3D-CSA (96.9 % for 2D-CSA vs 56.3 % for 3D-CSA).

Abstract NG08: Dissecting and quantifying pancreatic cancer plasticity using single-cell multiomics, lineage tracing and functional genomics reveals novel mediators of therapy resistance

Abstract NG08: Dissecting and quantifying pancreatic cancer plasticity using single-cell multiomics, lineage tracing and functional genomics reveals novel mediators of therapy resistance

A comprehensive analysis of pneumococcal two-component system regulatory networks.

Two-component systems are key signal-transduction systems that enable bacteria to respond to a wide variety of environmental stimuli. The human pathogen, Streptococcus pneumoniae (pneumococcus) encodes 13 two-component systems and a single orphan response regulator, most of which are significant for pneumococcal pathogenicity. Mapping the regulatory networks governed by these systems is key to understand pneumococcal host adaptation. Here we employ a novel bioinformatic approach to predict the regulons of each two-component system based on publicly available whole-genome sequencing data. By employing pangenome-wide association studies (panGWAS) to predict genotype-genotype associations for each two-component system, we predicted regulon genes of 11 of the pneumococcal two-component systems. Through validation via next-generation RNA-sequencing on response regulator overexpression mutants, several top candidate genes predicted by the panGWAS analysis were confirmed as regulon genes. The present study presents novel details on multiple pneumococcal two-component systems, including an expansion of regulons, identification of candidate response regulator binding motifs, and identification of candidate response regulator-regulated small non-coding RNAs. We also demonstrate a use for panGWAS as a complementary tool in target gene identification via identification of genotype-to-genotype links. Expanding our knowledge on two-component systems in pathogens is crucial to understanding how these bacteria sense and respond to their host environment, which could prove useful in future drug development.

The Putative GATA Transcription Factor SbGATA22 as a Novel Regulator of Dhurrin Biosynthesis.

Cyanogenic glucosides are specialized metabolites produced by over 3000 species of higher plants from more than 130 families. The deployment of cyanogenic glucosides is influenced by biotic and abiotic factors in addition to being developmentally regulated, consistent with their roles in plant defense and stress mitigation. Despite their ubiquity, very little is known regarding the molecular mechanisms that regulate their biosynthesis. The biosynthetic pathway of dhurrin, the cyanogenic glucoside found in the important cereal crop sorghum (Sorghum bicolor (L.) Moench), was described over 20 years ago, and yet no direct regulator of the biosynthetic genes has been identified. To isolate regulatory proteins that bind to the promoter region of the key dhurrin biosynthetic gene of sorghum, SbCYP79A1, yeast one-hybrid screens were performed. A bait fragment containing 1204 base pairs of the SbCYP79A1 5' regulatory region was cloned upstream of a reporter gene and introduced into Saccharomyces cerevisiae. Subsequently, the yeast was transformed with library cDNA representing RNA from two different sorghum developmental stages. From these screens, we identified SbGATA22, an LLM domain B-GATA transcription factor that binds to the putative GATA transcription factor binding motifs in the SbCYP79A1 promoter region. Transient assays in Nicotiana benthamiana show that SbGATA22 localizes to the nucleus. The expression of SbGATA22, in comparison with SbCYP79A1 expression and dhurrin concentration, was analyzed over 14 days of sorghum development and in response to nitrogen application, as these conditions are known to affect dhurrin levels. Collectively, these findings suggest that SbGATA22 may act as a negative regulator of SbCYP79A1 expression and provide a preliminary insight into the molecular regulation of dhurrin biosynthesis in sorghum.

WNK kinase is a vasoactive chloride sensor in endothelial cells

Endothelial cells (ECs) line the wall of blood vessels and regulate arterial contractility to tune regional organ blood flow and systemic pressure. Chloride (Cl-) is the most abundant anion in ECs and the Cl- sensitive With-No-Lysine (WNK) kinase is expressed in this cell type. Whether intracellular Cl- signaling and WNK kinase regulate EC function to alter arterial contractility is unclear. Here, we tested the hypothesis that intracellular Cl- signaling in ECs regulates arterial contractility and examined the signaling mechanisms involved, including the participation of WNK kinase. Our data obtained using two-photon microscopy and cell-specific inducible knockout mice indicated that acetylcholine, a prototypical vasodilator, stimulated a rapid reduction in intracellular Cl- concentration ([Cl-]i) due to the activation of TMEM16A, a Cl- channel, in ECs of resistance-size arteries. TMEM16A channel-mediated Cl- signaling activated WNK kinase, which phosphorylated its substrate proteins SPAK and OSR1 in ECs. OSR1 potentiated transient receptor potential vanilloid 4 (TRPV4) currents in a kinase-dependent manner and required a conserved binding motif located in the channel C terminus. Intracellular Ca2+ signaling was measured in four dimensions in ECs using a high-speed lightsheet microscope. WNK kinase-dependent activation of TRPV4 channels increased local intracellular Ca2+ signaling in ECs and produced vasodilation. In summary, we show that TMEM16A channel activation reduces [Cl-]i, which activates WNK kinase in ECs. WNK kinase phosphorylates OSR1 which then stimulates TRPV4 channels to produce vasodilation. Thus, TMEM16A channels regulate intracellular Cl- signaling and WNK kinase activity in ECs to control arterial contractility.

Extensive Loss and Gain of Conserved Noncoding Elements During Early Teleost Evolution.

Conserved noncoding elements in vertebrates are enriched around transcription factor loci associated with development. However, loss and rapid divergence of conserved noncoding elements has been reported in teleost fish, albeit taking only few genomes into consideration. Taking advantage of the recent increase in high-quality teleost genomes, we focus on studying the evolution of teleost conserved noncoding elements, carrying out targeted genomic alignments and comparisons within the teleost phylogeny to detect conserved noncoding elements and reconstruct the ancestral teleost conserved noncoding elements repertoire. This teleost-centric approach confirms previous observations of extensive vertebrate conserved noncoding elements loss early in teleost evolution, but also reveals massive conserved noncoding elements gain in the teleost stem-group over 300 million years ago. Using synteny-based association to link conserved noncoding elements to their putatively regulated target genes, we show the most teleost gained conserved noncoding elements are found in the vicinity of orthologous loci involved in transcriptional regulation and embryonic development that are also associated with conserved noncoding elements in other vertebrates. Moreover, teleost and vertebrate conserved noncoding elements share a highly similar motif and transcription factor binding site vocabulary. We suggest that early teleost conserved noncoding element gains reflect a restructuring of the ancestral conserved noncoding element repertoire through both extreme divergence and de novo emergence. Finally, we support newly identified pan-teleost conserved noncoding elements have potential for accurate resolution of teleost phylogenetic placements in par with coding sequences, unlike ancestral only elements shared with spotted gar. This work provides new insight into conserved noncoding element evolution with great value for follow-up work on phylogenomics, comparative genomics, and the study of gene regulation evolution in teleosts.

14-3-3 binding motif phosphorylation disrupts Hdac4-organized condensates to stimulate cardiac reprogramming

Cell fate conversion is associated with extensive post-translational modifications (PTMs) and architectural changes of sub-organelles, yet how these events are interconnected remains unknown. We report here the identification of a phosphorylation code in 14-3-3 binding motifs (PC14-3-3) that greatly stimulates induced cardiomyocyte (iCM) formation from fibroblasts. PC14-3-3 is identified in pivotal functional proteins for iCM reprogramming, including transcription factors and chromatin modifiers. Akt1 kinase and protein phosphatase 2A are the key writer and key eraser of the PC14-3-3 code, respectively. PC14-3-3 activation induces iCM formation with the presence of only Tbx5. In contrast, PC14-3-3 inhibition by mutagenesis or inhibitor-mediated code removal abolishes reprogramming. We discover that key PC14-3-3-embedded factors, such as histone deacetylase 4 (Hdac4), Mef2c, and Foxo1, form Hdac4-organized inhibitory nuclear condensates. PC14-3-3 activation disrupts Hdac4 condensates to promote cardiac gene expression. Our study suggests that sub-organelle dynamics regulated by a PTM code could be a general mechanism for stimulating cell reprogramming.

The electrostatic landscape of MHC-peptide binding revealed using inception networks

Predictive modeling of macromolecular recognition and protein-protein complementarity represents one of the cornerstones of biophysical sciences. However, such models are often hindered by the combinatorial complexity of interactions at the molecular interfaces. Exemplary of this problem is peptide presentation by the highly polymorphic major histocompatibility complex class I (MHC-I) molecule, a principal component of immune recognition. We developed human leukocyte antigen (HLA)-Inception, a deep biophysical convolutional neural network, which integrates molecular electrostatics to capture non-bonded interactions for predicting peptide binding motifs across 5,821 MHC-I alleles. These predictions of generated motifs correlate strongly with experimental peptide binding and presentation data. Beyond molecular interactions, the study demonstrates the application of predicted motifs in analyzing MHC-I allele associations with HIV disease progression and patient response to immune checkpoint inhibitors. A record of this paper's transparent peer review process is included in the supplemental information.

Wemics: A Single-Base Resolution Methylation Quantification Method for Enhanced Prediction of Epigenetic Regulation.

DNA methylation, an epigenetic mechanism that alters gene expression without changing DNA sequence, is essential for organism development and key biological processes like genomic imprinting and X-chromosome inactivation. Despite tremendous efforts in DNA methylation research, accurate quantification of cytosine methylation remains a challenge. Here, a single-base methylation quantification approach is introduced by weighting methylation of consecutive CpG sites (Wemics) in genomic regions. Wemics quantification of DNA methylation better predicts its regulatory impact on gene transcription and identifies differentially methylated regions (DMRs) with more biological relevance. Most Wemics-quantified DMRs in lung cancer are epigenetically conserved and recurrently occurred in other primary cancers from The Cancer Genome Atlas (TCGA), and their aberrant alterations can serve as promising pan-cancer diagnostic markers. It is further revealed that these detected DMRs are enriched in transcription factor (TF) binding motifs, and methylation of these TF binding motifs and TF expression synergistically regulate target gene expression. Using Wemics on epigenomic-transcriptomic data from the large lung cancer cohort, a dozen novel genes with oncogenic potential are discovered that are upregulated by hypomethylation but overlooked by other quantification methods. These findings increase the understanding of the epigenetic mechanism by which DNA methylation regulates gene expression.

TripHLApan: predicting HLA molecules binding peptides based on triple coding matrix and transfer learning.

Human leukocyte antigen (HLA) recognizes foreign threats and triggers immune responses by presenting peptides to T cells. Computationally modeling the binding patterns between peptide and HLA is very important for the development of tumor vaccines. However, it is still a big challenge to accurately predict HLA molecules binding peptides. In this paper, we develop a new model TripHLApan for predicting HLA molecules binding peptides by integrating triple coding matrix, BiGRU + Attention models, and transfer learning strategy. We have found the main interaction site regions between HLA molecules and peptides, as well as the correlation between HLA encoding and binding motifs. Based on the discovery, we make the preprocessing and coding closer to the natural biological process. Besides, due to the input being based on multiple types of features and the attention module focused on the BiGRU hidden layer, TripHLApan has learned more sequence level binding information. The application of transfer learning strategies ensures the accuracy of prediction results under special lengths (peptides in length 8) and model scalability with the data explosion. Compared with the current optimal models, TripHLApan exhibits strong predictive performance in various prediction environments with different positive and negative sample ratios. In addition, we validate the superiority and scalability of TripHLApan's predictive performance using additional latest data sets, ablation experiments and binding reconstitution ability in the samples of a melanoma patient. The results show that TripHLApan is a powerful tool for predicting the binding of HLA-I and HLA-II molecular peptides for the synthesis of tumor vaccines. TripHLApan is publicly available at https://github.com/CSUBioGroup/TripHLApan.git.

Single-cell multi-omics analysis identifies context-specific gene regulatory gates and mechanisms.

There is a growing interest in inferring context specific gene regulatory networks from single-cell RNA sequencing (scRNA-seq) data. This involves identifying the regulatory relationships between transcription factors (TFs) and genes in individual cells, and then characterizing these relationships at the level of specific cell types or cell states. In this study, we introduce scGATE (single-cell gene regulatory gate) as a novel computational tool for inferring TF-gene interaction networks and reconstructing Boolean logic gates involving regulatory TFs using scRNA-seq data. In contrast to current Boolean models, scGATE eliminates the need for individual formulations and likelihood calculations for each Boolean rule (e.g. AND, OR, XOR). By employing a Bayesian framework, scGATE infers the Boolean rule after fitting the model to the data, resulting in significant reductions in time-complexities for logic-based studies. We have applied assay for transposase-accessible chromatin with sequencing (scATAC-seq) data and TF DNA binding motifs to filter out non-relevant TFs in gene regulations. By integrating single-cell clustering with these external cues, scGATE is able to infer context specific networks. The performance of scGATE is evaluated using synthetic and real single-cell multi-omics data from mouse tissues and human blood, demonstrating its superiority over existing tools for reconstructing TF-gene networks. Additionally, scGATE provides a flexible framework for understanding the complex combinatorial and cooperative relationships among TFs regulating target genes by inferring Boolean logic gates among them.

Cardiac GR Mediates the Diurnal Rhythm in Ventricular Arrhythmia Susceptibility.

Ventricular arrhythmias (VAs) demonstrate a prominent day-night rhythm, commonly presenting in the morning. Transcriptional rhythms in cardiac ion channels accompany this phenomenon, but their role in the morning vulnerability to VAs and the underlying mechanisms are not understood. We investigated the recruitment of transcription factors that underpins transcriptional rhythms in ion channels and assessed whether this mechanism was pertinent to the heart's intrinsic diurnal susceptibility to VA. Assay for transposase-accessible chromatin with sequencing performed in mouse ventricular myocyte nuclei at the beginning of the animals' inactive (ZT0) and active (ZT12) periods revealed differentially accessible chromatin sites annotating to rhythmically transcribed ion channels and distinct transcription factor binding motifs in these regions. Notably, motif enrichment for the glucocorticoid receptor (GR; transcriptional effector of corticosteroid signaling) in open chromatin profiles at ZT12 was observed, in line with the well-recognized ZT12 peak in circulating corticosteroids. Molecular, electrophysiological, and in silico biophysically-detailed modeling approaches demonstrated GR-mediated transcriptional control of ion channels (including Scn5a underlying the cardiac Na+ current, Kcnh2 underlying the rapid delayed rectifier K+ current, and Gja1 responsible for electrical coupling) and their contribution to the day-night rhythm in the vulnerability to VA. Strikingly, both pharmacological block of GR and cardiomyocyte-specific genetic knockout of GR blunted or abolished ion channel expression rhythms and abolished the ZT12 susceptibility to pacing-induced VA in isolated hearts. Our study registers a day-night rhythm in chromatin accessibility that accompanies diurnal cycles in ventricular myocytes. Our approaches directly implicate the cardiac GR in the myocyte excitability rhythm and mechanistically link the ZT12 surge in glucocorticoids to intrinsic VA propensity at this time.