Structural insights into the interaction between adenovirus C5 hexon and human lactoferrin.

Adenovirus (AdV) infection of the respiratory epithelium is common but poorly understood. Human AdV species C types, such as HAdV-C5, utilize the Coxsackie-adenovirus receptor (CAR) for attachment and subsequently integrins for entry. CAR and integrins are however located deep within the tight junctions in the mucosa where they would not be easily accessible. Recently, a model for CAR-independent AdV entry was proposed. In this model, human lactoferrin (hLF), an innate immune protein, aids the viral uptake into epithelial cells by mediating interactions between the major capsid protein, hexon, and yet unknown host cellular receptor(s). However, a detailed understanding of the molecular interactions driving this mechanism is lacking. Here, we present a new cryo-EM structure of HAdV-5C hexon at high resolution alongside a hybrid structure of HAdV-5C hexon complexed with human lactoferrin (hLF). These structures reveal the molecular determinants of the interaction between hLF and HAdV-C5 hexon. hLF engages hexon primarily via its N-terminal lactoferricin (Lfcin) region, interacting with hexon's hypervariable region 1 (HVR-1). Mutational analyses pinpoint critical Lfcin contacts and also identify additional regions within hLF that critically contribute to hexon binding. Our study sheds more light on the intricate mechanism by which HAdV-C5 utilizes soluble hLF/Lfcin for cellular entry. These findings hold promise for advancing gene therapy applications and inform vaccine development. IMPORTANCE: Our study delves into the structural aspects of adenovirus (AdV) infections, specifically HAdV-C5 in the respiratory epithelium. It uncovers the molecular details of a novel pathway where human lactoferrin (hLF) interacts with the major capsid protein, hexon, facilitating viral entry, and bypassing traditional receptors such as CAR and integrins. The study's cryo-EM structures reveal how hLF engages hexon, primarily through its N-terminal lactoferricin (Lfcin) region and hexon's hypervariable region 1 (HVR-1). Mutational analyses identify critical Lfcin contacts and other regions within hLF vital for hexon binding. This structural insight sheds light on HAdV-C5's mechanism of utilizing soluble hLF/Lfcin for cellular entry, holding promise for gene therapy and vaccine development advancements in adenovirus research.

Human 2'-Deoxynucleoside 5'-Phosphate N-Hydrolase 1: The Catalytic Roles of Tyr24 and Asp80.

The human enzyme 2'-deoxynucleoside 5'-phosphate N-hydrolase 1 (HsDNPH1) catalyses the hydrolysis of 5-hydroxymethyl-2'-deoxyuridine 5'-phosphate to generate 5-hydroxymethyluracil and 2-deoxyribose-5-phosphate via a covalent 5-phospho-2-deoxyribosylated enzyme intermediate. HsDNPH1 is a promising target for inhibitor development towards anticancer drugs. Here, site-directed mutagenesis of conserved active-site residues, followed by HPLC analysis of the reaction and steady-state kinetics are employed to reveal the importance of each of these residues in catalysis, and the reaction pH-dependence is perturbed by each mutation. Solvent deuterium isotope effects indicate no rate-limiting proton transfers. Crystal structures of D80N-HsDNPH1 in unliganded and substrate-bound states, and of unliganded D80A- and Y24F-HsDNPH1 offer atomic level insights into substrate binding and catalysis. The results reveal a network of hydrogen bonds involving the substrate and the E104-Y24-D80 catalytic triad and are consistent with a proposed mechanism whereby D80 is important for substrate positioning, for helping modulate E104 nucleophilicity, and as the general acid in the first half-reaction. Y24 positions E104 for catalysis and prevents a catalytically disruptive close contact between E104 and D80.

Polymorphisms in transcription factor binding sites and enhancer regions and pancreatic ductal adenocarcinoma risk.

Genome-wide association studies (GWAS) are a powerful tool for detecting variants associated with complex traits and can help risk stratification and prevention strategies against pancreatic ductal adenocarcinoma (PDAC). However, the strict significance threshold commonly used makes it likely that many true risk loci are missed. Functional annotation of GWAS polymorphisms is a proven strategy to identify additional risk loci. We aimed to investigate single-nucleotide polymorphisms (SNP) in regulatory regions [transcription factor binding sites (TFBSs) and enhancers] that could change the expression profile of multiple genes they act upon and thereby modify PDAC risk. We analyzed a total of 12,636 PDAC cases and 43,443 controls from PanScan/PanC4 and the East Asian GWAS (discovery populations), and the PANDoRA consortium (replication population). We identified four associations that reached study-wide statistical significance in the overall meta-analysis: rs2472632(A) (enhancer variant, OR 1.10, 95%CI 1.06,1.13, p = 5.5 × 10-8), rs17358295(G) (enhancer variant, OR 1.16, 95%CI 1.10,1.22, p = 6.1 × 10-7), rs2232079(T) (TFBS variant, OR 0.88, 95%CI 0.83,0.93, p = 6.4 × 10-6) and rs10025845(A) (TFBS variant, OR 1.88, 95%CI 1.50,1.12, p = 1.32 × 10-5). The SNP with the most significant association, rs2472632, is located in an enhancer predicted to target the coiled-coil domain containing 34 oncogene. Our results provide new insights into genetic risk factors for PDAC by a focused analysis of polymorphisms in regulatory regions and demonstrating the usefulness of functional prioritization to identify loci associated with PDAC risk.

Altered Protein Dynamics and a More Reactive Catalytic Cysteine in a Neurodegeneration-associated UCHL1 Mutant.

A mutant of ubiquitin C-terminal hydrolase L1 (UCHL1) detected in early-onset neurodegenerative patients, UCHL1R178Q, showed higher catalytic activity than wild-type UCHL1 (UCHL1WT). Lying within the active-site pocket, the arginine is part of an interaction network that holds the catalytic histidine in an inactive arrangement. However, the structural basis and mechanism of enzymatic activation upon glutamine substitution was not understood. We combined X-ray crystallography, protein nuclear magnetic resonance (NMR) analysis, enzyme kinetics, covalent inhibition analysis, and biophysical measurements to delineate activating factors in the mutant. While the crystal structure of UCHL1R178Q showed nearly the same arrangement of the catalytic residues and active-site pocket, the mutation caused extensive alteration in the chemical environment and dynamics of more than 30 residues, some as far as 15 Å away from the site of mutation. Significant broadening of backbone amide resonances in the HSQC spectra indicates considerable backbone dynamics changes in several residues, in agreement with solution small-angle X-ray scattering (SAXS) analyses which indicate an overall increase in protein flexibility. Enzyme kinetics show the activation is due to a kcat effect despite a slightly weakened substrate affinity. In line with this, the mutant shows a higher second-order rate constant (kinact/Ki) in a reaction with a substrate-derived irreversible inhibitor, Ub-VME, compared to the wild-type enzyme, an observation indicative of a more reactive catalytic cysteine in the mutant. Together, the observations underscore structural plasticity as a factor contributing to enzyme kinetic behavior which can be modulated through mutational effects.

In Silico Prediction and Molecular Docking of SNPs in NRP1 Gene Associated with SARS-COV-2.

Neuropilin-1 (NRP1) which is a main transmembrane cell surface receptor acts as a host cell mediator resulting in increasing the SARS-Cov-2 infectivity and also plays a role in neuronal development, angiogenesis and axonal outgrowth. The goal of this study is to estimate the impact of single nucleotide polymorphisms (SNPs) in the NRP1 gene on the function, structure and stabilization of protein as well as on the miRNA-mRNA binding regions using bioinformatical tools. It is also aimed to investigate the changes caused by SNPs in NRP1 on interactions with drug molecule and spike protein. The missense type of SNPs was analyzed using SIFT, PolyPhen-2, SNAP2, PROVEAN, Mutation Assessor, SNPs&GO, PhD-SNP, I-Mutant 3.0, MUpro, STRING, Project HOPE, ConSurf, and PolymiRTS. Docking analyses were conducted by AutoDock Vina program. As a result, a total of 733 missense SNPs were determined within the NRP1 gene and nine SNPs were specified as damaging to the protein. The modelling results showed that wild and mutant type amino acids had some different properties such as size, charge, and hydrophobicity. Additionally, their three-dimensional structures of protein were utilized for confirmation of these differences. After evaluating the results, nine polymorphisms rs141633354, rs142121081, rs145954532, rs200028992, rs200660300, rs369312020, rs370117610, rs370551432, rs370641686 were determined to be damaging on the structure and function of NRP1 protein and located in conserved regions. The results of molecular docking showed that the binding affinity values are nearly the same for wild-type and mutant structures support that the mutations carried out are not in the focus of the binding site, therefore the ligand does not affect the binding energy. It is expected that the results will be useful for future studies.

Suppression of KLF5 basal expression in oral carcinoma-derived cells through three intact CREB1-binding sites in the silencer region.

OBJECTIVES: Krüppel-like factor (KLF)5, which is overexpressed in carcinomas such as oral cancer, inhibits epidermal differentiation. KLF5 induces dedifferentiation of carcinoma cells, which effectuates carcinoma progression; nevertheless, the regulatory mechanism affecting the transcription of the KLF5 gene remains ambiguous. METHODS: Transcriptional activity of the KLF5 silencer, specifically the 425-bp region (425-region), was examined using reporter assays. An additional analysis was conducted to assess the impact of the minimal essential region (MER) of KLF5 on its basal expression. The affinity of cAMP responsive element binding protein 1 (CREB1) for three potential CREB1-binding sites in the 425-region was analyzed using DNA pull-down and quantitative chromatin immunoprecipitation assays. Reporter assays employing a human oral squamous carcinoma cell line, HSC2, transfected with small interfering RNA or complementary DNA for CREB1, were performed to investigate the effect of CREB1 binding sites on MER activity. RESULTS: The 425-region exhibited no transcriptional activity and suppressed MER transcriptional activity. This region encodes three putative CREB1-binding sites, and CREB1 demonstrated equal binding affinity for all three sites. The deletion of each of these binding sites reduced CREB1 precipitation and enhanced MER activity. Endogenous CREB1 knockdown and overexpression elevated and reduced MER activity, respectively, at the intact sites. Conversely, site deletion hampered and improved MER activity upon CREB1 knockdown and overexpression, respectively. CONCLUSIONS: Suppression of KLF5 basal expression via CREB1 binding to the 425-region requires all three CREB1-binding sites to remain intact in oral carcinoma cells. Consequently, deletion of the CREB1-binding site relieves suppression of KLF5 basal expression.

The forkhead DNA-binding domain binds specific G2-rich RNA sequences.

Transcription factors contain a DNA-binding domain ensuring specific recognition of DNA target sequences. The family of forkhead (FOX) transcription factors is composed of dozens of paralogs in mammals. The forkhead domain (FHD) is a segment of about 100 amino acids that binds an A-rich DNA sequence. Using DNA and RNA PCR-SELEX, we show that recombinant FOXL2 proteins, either wild-type or carrying the oncogenic variant C134W, recognize similar DNA-binding sites. This suggests that the oncogenic variant does not alter the intrinsic sequence-specificity of FOXL2. Most importantly, we show that FOXL2 binds G2-rich RNA sequences whereas it virtually fails to bind similar sequences in DNA chemistry. Interestingly, a statistically significant subset of genes responding to the knock-down of FOXL2/Foxl2 harbor such G2-rich sequences and are involved in crucial signaling pathways and cellular processes. In addition, we show that FOXA1, FOXO3a and chimeric FOXL2 proteins containing the FHD of the former are also able to interact with some of the preferred FOXL2-binding sequences. Our results point to an unexpected and novel characteristic of the forkhead domain, the biological relevance of which remains to be explored.

Identification of transcription factor high accumulation DNA zones.

BACKGROUND: Transcription factors (TF) play a crucial role in the regulation of gene transcription; alterations of their activity and binding to DNA areas are strongly involved in cancer and other disease onset and development. For proper biomedical investigation, it is hence essential to correctly trace TF dense DNA areas, having multiple bindings of distinct factors, and select DNA high occupancy target (HOT) zones, showing the highest accumulation of such bindings. Indeed, systematic and replicable analysis of HOT zones in a large variety of cells and tissues would allow further understanding of their characteristics and could clarify their functional role. RESULTS: Here, we propose, thoroughly explain and discuss a full computational procedure to study in-depth DNA dense areas of transcription factor accumulation and identify HOT zones. This methodology, developed as a computationally efficient parametric algorithm implemented in an R/Bioconductor package, uses a systematic approach with two alternative methods to examine transcription factor bindings and provide comparative and fully-reproducible assessments. It offers different resolutions by introducing three distinct types of accumulation, which can analyze DNA from single-base to region-oriented levels, and a moving window, which can estimate the influence of the neighborhood for each DNA base under exam. CONCLUSIONS: We quantitatively assessed the full procedure by using our implemented software package, named TFHAZ, in two example applications of biological interest, proving its full reliability and relevance.

Genome-wide census of ATF4 binding sites and functional profiling of trait-associated genetic variants overlapping ATF4 binding motifs.

Activating Transcription Factor 4 (ATF4) is an important regulator of gene expression in stress responses and developmental processes in many cell types. Here, we catalogued ATF4 binding sites in the human genome and identified overlaps with trait-associated genetic variants. We probed these genetic variants for allelic regulatory activity using a massively parallel reporter assay (MPRA) in HepG2 hepatoma cells exposed to tunicamycin to induce endoplasmic reticulum stress and ATF4 upregulation. The results revealed that in the majority of cases, the MPRA allelic activity of these SNPs was in agreement with the nucleotide preference seen in the ATF4 binding motif from ChIP-Seq. Luciferase and electrophoretic mobility shift assays in additional cellular models further confirmed ATF4-dependent regulatory effects for the SNPs rs532446 (GADD45A intronic; linked to hematological parameters), rs7011846 (LPL upstream; myocardial infarction), rs2718215 (diastolic blood pressure), rs281758 (psychiatric disorders) and rs6491544 (educational attainment). CRISPR-Cas9 disruption and/or deletion of the regulatory elements harboring rs532446 and rs7011846 led to the downregulation of GADD45A and LPL, respectively. Thus, these SNPs could represent examples of GWAS genetic variants that affect gene expression by altering ATF4-mediated transcriptional activation.

High-throughput data and modeling reveal insights into the mechanisms of cooperative DNA-binding by transcription factor proteins.

Cooperative DNA-binding by transcription factor (TF) proteins is critical for eukaryotic gene regulation. In the human genome, many regulatory regions contain TF-binding sites in close proximity to each other, which can facilitate cooperative interactions. However, binding site proximity does not necessarily imply cooperative binding, as TFs can also bind independently to each of their neighboring target sites. Currently, the rules that drive cooperative TF binding are not well understood. In addition, it is oftentimes difficult to infer direct TF-TF cooperativity from existing DNA-binding data. Here, we show that in vitro binding assays using DNA libraries of a few thousand genomic sequences with putative cooperative TF-binding events can be used to develop accurate models of cooperativity and to gain insights into cooperative binding mechanisms. Using factors ETS1 and RUNX1 as our case study, we show that the distance and orientation between ETS1 sites are critical determinants of cooperative ETS1-ETS1 binding, while cooperative ETS1-RUNX1 interactions show more flexibility in distance and orientation and can be accurately predicted based on the affinity and sequence/shape features of the binding sites. The approach described here, combining custom experimental design with machine-learning modeling, can be easily applied to study the cooperative DNA-binding patterns of any TFs.

Androgen receptor binding sites enabling genetic prediction of mortality due to prostate cancer in cancer-free subjects.

Prostate cancer (PrCa) is the second most common cancer worldwide in males. While strongly warranted, the prediction of mortality risk due to PrCa, especially before its development, is challenging. Here, we address this issue by maximizing the statistical power of genetic data with multi-ancestry meta-analysis and focusing on binding sites of the androgen receptor (AR), which has a critical role in PrCa. Taking advantage of large Japanese samples ever, a multi-ancestry meta-analysis comprising more than 300,000 subjects in total identifies 9 unreported loci including ZFHX3, a tumor suppressor gene, and successfully narrows down the statistically finemapped variants compared to European-only studies, and these variants strongly enrich in AR binding sites. A polygenic risk scores (PRS) analysis restricting to statistically finemapped variants in AR binding sites shows among cancer-free subjects, individuals with a PRS in the top 10% have a strongly higher risk of the future death of PrCa (HR: 5.57, P = 4.2 × 10-10). Our findings demonstrate the potential utility of leveraging large-scale genetic data and advanced analytical methods in predicting the mortality of PrCa.

Polypyrimidine tract binding protein 1 (PTBP1) contains a novel regulatory sequence, the rBH3, that binds the prosurvival protein MCL1.

The maturation of RNA from its nascent transcription to ultimate utilization (e.g., translation, miR-mediated RNA silencing, etc.) involves an intricately coordinated series of biochemical reactions regulated by RNA-binding proteins (RBPs). Over the past several decades, there has been extensive effort to elucidate the biological factors that control specificity and selectivity of RNA target binding and downstream function. Polypyrimidine tract binding protein 1 (PTBP1) is an RBP that is involved in all steps of RNA maturation and serves as a key regulator of alternative splicing, and therefore, understanding its regulation is of critical biologic importance. While several mechanisms of RBP specificity have been proposed (e.g., cell-specific expression of RBPs and secondary structure of target RNA), recently, protein-protein interactions with individual domains of RBPs have been suggested to be important determinants of downstream function. Here, we demonstrate a novel binding interaction between the first RNA recognition motif 1 (RRM1) of PTBP1 and the prosurvival protein myeloid cell leukemia-1 (MCL1). Using both in silico and in vitro analyses, we demonstrate that MCL1 binds a novel regulatory sequence on RRM1. NMR spectroscopy reveals that this interaction allosterically perturbs key residues in the RNA-binding interface of RRM1 and negatively impacts RRM1 association with target RNA. Furthermore, pulldown of MCL1 by endogenous PTBP1 verifies that these proteins interact in an endogenous cellular environment, establishing the biological relevance of this binding event. Overall, our findings suggest a novel mechanism of regulation of PTBP1 in which a protein-protein interaction with a single RRM can impact RNA association.

Sequence-Dependent Interaction of the Human Papillomavirus E2 Protein with the DNA Elements on Its DNA Replication Origin.

The human papillomavirus (HPV) E2 protein is essential for regulating the initiation of viral DNA replication as well as the regulation of transcription of certain HPV-encoded genes. Its ability to recognize and bind to its four recognition sequences in the viral origin is a key step in the initiation of HPV DNA replication. Thus, understanding the mechanism of DNA binding by E2 protein and the unique roles played by individual DNA sequence elements of the replication origin is essential. We have purified the recombinant full-length HPV type 11 E2 protein. Quantitative DNA binding analysis indicated E2 protein bound all four DNA binding sites with reasonably high affinities but with distinct preferences. It bound its cognate binding sites 1, 2, and 4 with higher affinities, but bound binding site 3 with lower affinity. Analysis of binding to these sites unraveled multiple sequence elements that appeared to influence E2 binding affinity and target discrimination, including the sequence of spacer region, flanking sequences, and proximity of E2 binding sites. Thermodynamic analysis indicated hydrophobic interaction in the protein-DNA complex formation. Our studies indicate a large multi-protein complex formation on the HPV-origin DNA, likely due to reasonably high binding affinities as well as intrinsic oligomerization propensity of E2 dimers.

UV irradiation remodels the specificity landscape of transcription factors.

Somatic mutations are highly enriched at transcription factor (TF) binding sites, with the strongest trend being observed for ultraviolet light (UV)-induced mutations in melanomas. One of the main mechanisms proposed for this hypermutation pattern is the inefficient repair of UV lesions within TF-binding sites, caused by competition between TFs bound to these lesions and the DNA repair proteins that must recognize the lesions to initiate repair. However, TF binding to UV-irradiated DNA is poorly characterized, and it is unclear whether TFs maintain specificity for their DNA sites after UV exposure. We developed UV-Bind, a high-throughput approach to investigate the impact of UV irradiation on protein-DNA binding specificity. We applied UV-Bind to ten TFs from eight structural families, and found that UV lesions significantly altered the DNA-binding preferences of all the TFs tested. The main effect was a decrease in binding specificity, but the precise effects and their magnitude differ across factors. Importantly, we found that despite the overall reduction in DNA-binding specificity in the presence of UV lesions, TFs can still compete with repair proteins for lesion recognition, in a manner consistent with their specificity for UV-irradiated DNA. In addition, for a subset of TFs, we identified a surprising but reproducible effect at certain nonconsensus DNA sequences, where UV irradiation leads to a high increase in the level of TF binding. These changes in DNA-binding specificity after UV irradiation, at both consensus and nonconsensus sites, have important implications for the regulatory and mutagenic roles of TFs in the cell.

MEN1 mutations mediate clinical resistance to menin inhibition.

Chromatin-binding proteins are critical regulators of cell state in haematopoiesis1,2. Acute leukaemias driven by rearrangement of the mixed lineage leukaemia 1 gene (KMT2Ar) or mutation of the nucleophosmin gene (NPM1) require the chromatin adapter protein menin, encoded by the MEN1 gene, to sustain aberrant leukaemogenic gene expression programs3-5. In a phase 1 first-in-human clinical trial, the menin inhibitor revumenib, which is designed to disrupt the menin-MLL1 interaction, induced clinical responses in patients with leukaemia with KMT2Ar or mutated NPM1 (ref. 6). Here we identified somatic mutations in MEN1 at the revumenib-menin interface in patients with acquired resistance to menin inhibition. Consistent with the genetic data in patients, inhibitor-menin interface mutations represent a conserved mechanism of therapeutic resistance in xenograft models and in an unbiased base-editor screen. These mutants attenuate drug-target binding by generating structural perturbations that impact small-molecule binding but not the interaction with the natural ligand MLL1, and prevent inhibitor-induced eviction of menin and MLL1 from chromatin. To our knowledge, this study is the first to demonstrate that a chromatin-targeting therapeutic drug exerts sufficient selection pressure in patients to drive the evolution of escape mutants that lead to sustained chromatin occupancy, suggesting a common mechanism of therapeutic resistance.

An Assessment of Quaternary Structure Functionality in Homomer Protein Complexes.

It has been recently suggested that a significant fraction of homomer protein-protein interfaces evolve neutrally, without contributing to function, due to a hydrophobic bias in missense mutations. However, the fraction of such gratuitous complexes is currently unknown. Here, we quantified the fraction of homodimers where multimerization is unlikely to contribute to their biochemical function. We show that: 1) ligand binding-site structure predicts whether a homomer is functional or not; the vast majority of homodimers with multichain binding-sites (MBS) are likely to be functional, while in homodimers with single-chain binding-sites (SBS) and small to medium interfaces, quaternary structure is unlikely to be functional in a significant fraction-35%, even up to 42%-of complexes; 2) the hydrophobicity of interfaces changes little with the strength of selection, and the amino acid composition of interfaces is shaped by the "hydrophobic ratchet" in both types, but they are not in a strict equilibrium with mutations; particularly cysteines are much more abundant in mutations than in interfaces or surfaces; 3) in MBS homomers, the interfaces are conserved, while in a high fraction of SBS homomers, the interface is not more conserved than the solvent-accessible surface; and 4) MBS homomer interfaces coevolve more strongly with ligand binding sites than the interfaces of SBS homomers, and MBS complexes have higher capacity to transfer information from ligands across the interfaces than SBS homomers, explaining the enrichment of allostery in the former.

Widespread perturbation of ETS factor binding sites in cancer.

Although >90% of somatic mutations reside in non-coding regions, few have been reported as cancer drivers. To predict driver non-coding variants (NCVs), we present a transcription factor (TF)-aware burden test based on a model of coherent TF function in promoters. We apply this test to NCVs from the Pan-Cancer Analysis of Whole Genomes cohort and predict 2555 driver NCVs in the promoters of 813 genes across 20 cancer types. These genes are enriched in cancer-related gene ontologies, essential genes, and genes associated with cancer prognosis. We find that 765 candidate driver NCVs alter transcriptional activity, 510 lead to differential binding of TF-cofactor regulatory complexes, and that they primarily impact the binding of ETS factors. Finally, we show that different NCVs within a promoter often affect transcriptional activity through shared mechanisms. Our integrated computational and experimental approach shows that cancer NCVs are widespread and that ETS factors are commonly disrupted.

CRISPR/Cas9-mediated mutations in both a cAMP response element and an ETS-binding site suppress FLT1 gene expression.

The Fms-like tyrosine kinase-1 (FLT1) gene is expressed in various types of cells, including vascular endothelial cells and placental trophoblasts, and regulates angiogenesis, inflammation, and pregnancy. However, the basal transcriptional machinery of FLT1 is still not well understood. In this study, we first examined FLT1 promoter activity in three different types of cells, that is, trophoblast-derived cells, vascular endothelial-related cells, and HEK293 cells, using plasmid-based luciferase reporter assays, and showed that a cAMP-response element (CRE) and an ETS-binding site (EBS) are important for FLT1 expression in all cell types. To further examine the importance of these sites at the chromosomal level using HEK293 cells, we introduced CRISPR/Cas9-mediated mutations in these sites on the genomic DNA. HEK293 cells carrying these mutations clearly showed a significant decrease in endogenous FLT1 gene expression. These results suggest that CRE and EBS transcription regulatory elements are crucial for FLT1 gene expression in human tissues.

Computational prediction and characterization of cell-type-specific and shared binding sites.

MOTIVATION: Cell-type-specific gene expression is maintained in large part by transcription factors (TFs) selectively binding to distinct sets of sites in different cell types. Recent research works have provided evidence that such cell-type-specific binding is determined by TF's intrinsic sequence preferences, cooperative interactions with co-factors, cell-type-specific chromatin landscapes and 3D chromatin interactions. However, computational prediction and characterization of cell-type-specific and shared binding sites is rarely studied. RESULTS: In this article, we propose two computational approaches for predicting and characterizing cell-type-specific and shared binding sites by integrating multiple types of features, in which one is based on XGBoost and another is based on convolutional neural network (CNN). To validate the performance of our proposed approaches, ChIP-seq datasets of 10 binding factors were collected from the GM12878 (lymphoblastoid) and K562 (erythroleukemic) human hematopoietic cell lines, each of which was further categorized into cell-type-specific (GM12878- and K562-specific) and shared binding sites. Then, multiple types of features for these binding sites were integrated to train the XGBoost- and CNN-based models. Experimental results show that our proposed approaches significantly outperform other competing methods on three classification tasks. Moreover, we identified independent feature contributions for cell-type-specific and shared sites through SHAP values and explored the ability of the CNN-based model to predict cell-type-specific and shared binding sites by excluding or including DNase signals. Furthermore, we investigated the generalization ability of our proposed approaches to different binding factors in the same cellular environment. AVAILABILITY AND IMPLEMENTATION: The source code is available at: https://github.com/turningpoint1988/CSSBS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A competitive precision CRISPR method to identify the fitness effects of transcription factor binding sites.

Here we describe a competitive genome editing method that measures the effect of mutations on molecular functions, based on precision CRISPR editing using template libraries with either the original or altered sequence, and a sequence tag, enabling direct comparison between original and mutated cells. Using the example of the MYC oncogene, we identify important transcriptional targets and show that E-box mutations at MYC target gene promoters reduce cellular fitness.