Search

Development of the malaria parasite, Plasmodium falciparum, is regulated by a limited number of sequence-specific transcription factors (TFs). However, the mechanisms by which these TFs recognize genome-wide binding sites is largely unknown. To address TF specificity, we investigated the binding of two TF subsets that either bind CACACA or GTGCAC DNA sequence motifs and further characterized two additional ApiAP2 TFs, PfAP2-G and PfAP2-EXP, which bind unique DNA motifs (GTAC and TGCATGCA). We also interrogated the impact of DNA sequence and chromatin context on P. falciparum TF binding by integrating high-throughput in vitro and in vivo binding assays, DNA shape predictions, epigenetic post-translational modifications, and chromatin accessibility. We found that DNA sequence context minimally impacts binding site selection for paralogous CACACA-binding TFs, while chromatin accessibility, epigenetic patterns, co-factor recruitment, and dimerization correlate with differential binding. In contrast, GTGCAC-binding TFs prefer different DNA sequence context in addition to chromatin dynamics. Finally, we determined that TFs that preferentially bind divergent DNA motifs may bind overlapping genomic regions due to low-affinity binding to other sequence motifs. Our results demonstrate that TF binding site selection relies on a combination of DNA sequence and chromatin features, thereby contributing to the complexity of P. falciparum gene regulatory mechanisms., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

How chromatin enzymes work in condensed chromatin and how they maintain diffusional mobility inside remains unexplored. Here we investigated these challenges using the Drosophila ISWI remodeling ATPase, which slides nucleosomes along DNA. Folding of chromatin fibers did not affect sliding in vitro. Catalytic rates were also comparable in- and outside of chromatin condensates. ISWI cross-links and thereby stiffens condensates, except when ATP hydrolysis is possible. Active hydrolysis is also required for ISWI's mobility in condensates. Energy from ATP hydrolysis therefore fuels ISWI's diffusion through chromatin and prevents ISWI from cross-linking chromatin. Molecular dynamics simulations of a 'monkey-bar' model in which ISWI grabs onto neighboring nucleosomes, then withdraws from one before rebinding another in an ATP hydrolysis-dependent manner, qualitatively agree with our data. We speculate that monkey-bar mechanisms could be shared with other chromatin factors and that changes in chromatin dynamics caused by mutations in remodelers could contribute to pathologies., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Transcription is a fundamental biological process of transferring genetic information which often occurs in stochastic bursts when periods of intense activity alternate with quiescent phases. Recent experiments identified strong correlations between the association of transcription factors (TFs) to gene promoters on DNA and transcriptional activity. However, the underlying molecular mechanisms of this phenomenon remain not well understood. Here, we present a theoretical framework that allowed us to investigate how binding dynamics of TF influences transcriptional bursting. Our minimal theoretical model incorporates the most relevant physical-chemical features, including TF exchange among multiple binding sites at gene promoters and TF association/dissociation dynamics. Using analytical calculations supported by Monte Carlo computer simulations, it is demonstrated that transcriptional bursting dynamics depends on the strength of TF binding and the number of binding sites. Stronger TF binding affinity prolongs burst duration but reduces variability, while an optimal number of binding sites maximizes transcriptional noise, facilitating cellular adaptation. Our theoretical method explains available experimental observations quantitatively, confirming the model's predictive accuracy. This study provides important insights into molecular mechanisms of gene expression and regulation, offering a new theoretical tool for understanding complex biological processes.

TRIP4 is a conserved transcriptional coactivator that is involved in the regulation of the expression of multiple genes. It consists of a classical N-terminal C2HC5-like zinc-finger domain and a conserved C-terminal ASCH domain. Here, we characterized the DNA-binding properties of the human TRIP4 ASCH domain. Our biochemical data show that TRIP4-ASCH has comparable binding affinities toward ssDNA and dsDNA of different lengths, sequences, and structures. The crystal structures reveal that TRIP4-ASCH binds to DNA substrates in a sequence-independent manner through two adjacent positively charged surface patches: one binds to the 5'-end of DNA, and the other binds to the 3'-end of DNA. Further mutagenesis experiments and binding assays confirm the functional roles of key residues involved in DNA binding. In summary, our data demonstrate that TRIP4-ASCH binds to the 5' and 3'-ends of DNA in a sequence-independent manner, which will facilitate further studies of the biological function of TRIP4., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Motivation: Transcription factors are pivotal in the regulation of gene expression, and accurate identification of transcription factor binding sites (TFBSs) at high resolution is crucial for understanding the mechanisms underlying gene regulation. The task of identifying TFBSs from DNA sequences is a significant challenge in the field of computational biology today. To address this challenge, a variety of computational approaches have been developed. However, these methods face limitations in their ability to achieve high-resolution identification and often lack interpretability., Results: We propose BertSNR, an interpretable deep learning framework for identifying TFBSs at single-nucleotide resolution. BertSNR integrates sequence-level and token-level information by multi-task learning based on pre-trained DNA language models. Benchmarking comparisons show that our BertSNR outperforms the existing state-of-the-art methods in TFBS predictions. Importantly, we enhanced the interpretability of the model through attentional weight visualization and motif analysis, and discovered the subtle relationship between attention weight and motif. Moreover, BertSNR effectively identifies TFBSs in promoter regions, facilitating the study of intricate gene regulation., Availability and Implementation: The BertSNR source code can be found at https://github.com/lhy0322/BertSNR., (© The Author(s) 2024. Published by Oxford University Press.)

Bacteria use specialized proteins, like transcription factors, to rapidly control metal ion balance. CueR is a Gram-negative bacterial copper regulator. The structure of E. coli CueR complexed with Cu(I) and DNA was published, since then many studies have shed light on its function. However, P. aeruginosa CueR, which shows high sequence similarity to E. coli CueR, has been less studied. Here, we applied room-temperature electron paramagnetic resonance (EPR) measurements to explore changes in dynamics of P. aeruginosa CueR in dependency of copper concentrations and interaction with two different DNA promoter regions. We showed that P. aeruginosa CueR is less dynamic than the E. coli CueR protein and exhibits much higher sensitivity to DNA binding as compared to its E. coli CueR homolog. Moreover, a difference in dynamical behavior was observed when P. aeruginosa CueR binds to the copZ2 DNA promoter sequence compared to the mexPQ-opmE promoter sequence. Such dynamical differences may affect the expression levels of CopZ2 and MexPQ-OpmE proteins in P. aeruginosa. Overall, such comparative measurements of protein-DNA complexes derived from different bacterial systems reveal insights about how structural and dynamical differences between two highly homologous proteins lead to quite different DNA sequence-recognition and mechanistic properties., (© 2024 The Authors. ChemBioChem published by Wiley-VCH GmbH.)

Recently, biomolecular condensates formed through liquid-liquid phase separation have been widely reported to regulate key intracellular processes involved in cell biology and pathogenesis. BRD4 is a nuclear protein instrumental to the establishment of phase-separated super-enhancers (SEs) to direct the transcription of important genes. We previously observed that protein droplets of BRD4 became hydrophobic as their size increase, implying an ability of SEs to limit the ionization of water molecules by irradiation. Here, we aim to establish if SEs confer radiation resistance in cancer cells. We established an in vitro DNA damage assay that measures the effect of radicals provoked by the Fenton reaction on DNA integrity. This revealed that DNA damage was markedly reduced when BRD4 underwent phase separation with DNA. Accordingly, co-focal imaging analyses revealed that SE foci and DNA damage foci are mutually exclusive in irradiated cells. Lastly, we observed that the radioresistance of cancer cells was significantly reduced when irradiation was combined with ARV-771, a BRD4 de-stabilizer. Our data revealed the existence of innately radioresistant genomic regions driven by phase separation in cancer cells. The disruption of these phase-separated components enfolding genomic DNA may represent a novel strategy to augment the effects of radiotherapy., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Japanese Radiation Research Society and Japanese Society for Radiation Oncology.)

The precise spatial and temporal orchestration of gene expression is crucial for the ontogeny of an organism and is mainly governed by transcription factors (TFs). The mechanism of recognition of cognate sites amid millions of base pairs in the genome by TFs is still incompletely understood. In this study, we focus on DNA sequence composition, shape, and flexibility preferences of 28 quintessential TFs from Drosophila melanogaster that are critical to development and body patterning mechanisms. Our study finds that TFs exhibit distinct predilections for DNA shape, flexibility, and sequence compositions in the proximity of transcription factor binding sites (TFBSs). Notably, certain zinc finger proteins prefer GC-rich areas with less negative propeller twist, while homeodomains mainly seek AT-rich regions with a more negative propeller twist at their sites. Intriguingly, while numerous cofactors share similar binding site preferences and bind closer to each other in the genome, some cofactors that have different preferences bind farther apart. These findings shed light on TF DNA recognition and provide novel insights into possible cofactor binding and transcriptional regulation mechanisms.

LIM homeodomain transcription factor 1-alpha (LMX1a) is a neuronal lineage-specific transcription activator that plays an essential role during the development of midbrain dopaminergic (mDA) neurons. LMX1a induces the expression of multiple key genes, which ultimately determine the morphology, physiology, and functional identity of mDA neurons. This function of LMX1a is dependent on its homeobox domain. Here, we determined the structures of the LMX1a homeobox domain in complex with the promoter sequences of the Wnt family member 1 (WNT1) or paired like homeodomain 3 (Pitx3) gene, respectively. The complex structures revealed that the LMX1a homeobox domain employed its α3 helix and an N-terminal loop to achieve specific target recognition. The N-terminal loop (loop1) interacted with the minor groove of the double-stranded DNA (dsDNA), whereas the third α-helix (α3) was tightly packed into the major groove of the dsDNA. Structure-based mutations in the α3 helix of the homeobox domain significantly reduced the binding affinity of LMX1a to dsDNA. Moreover, we identified a nonsyndromic hearing loss (NSHL)-related mutation, R199, which yielded a more flexible loop and disturbed the recognition in the minor groove of dsDNA, consistent with the molecular dynamics (MD) simulations. Furthermore, overexpression of Lmx1a promoted the differentiation of SH-SY5Y cells and upregulated the transcription of WNT1 and PITX3 genes. Hence, our work provides a detailed elucidation of the specific recognition between the LMX1a homeobox domain and its specific dsDNA targets, which represents valuable information for future investigations of the functional pathways that are controlled by LMX1a during mDA neuron development., (© 2024 Federation of European Biochemical Societies.)

AT-rich interacting domain (ARID)-containing proteins, Arids, are a heterogeneous DNA-binding protein family involved in transcription regulation and chromatin processing. For the member Arid5a, no exact DNA-binding preference has been experimentally defined so far. Additionally, the protein binds to mRNA motifs for transcript stabilization, supposedly through the DNA-binding ARID domain. To date, however, no unbiased RNA motif definition and clear dissection of nucleic acid-binding through the ARID domain have been undertaken. Using NMR-centered biochemistry, we here define the Arid5a DNA preference. Further, high-throughput in vitro binding reveals a consensus RNA-binding motif engaged by the core ARID domain. Finally, transcriptome-wide binding (iCLIP2) reveals that Arid5a has a weak preference for (A)U-rich regions in pre-mRNA transcripts of factors related to RNA processing. We find that the intrinsically disordered regions flanking the ARID domain modulate the specificity and affinity of DNA binding, while they appear crucial for RNA interactions. Ultimately, our data suggest that Arid5a uses its extended ARID domain for bifunctional gene regulation and that the involvement of IDR extensions is a more general feature of Arids in interacting with different nucleic acids at the chromatin-mRNA interface., Competing Interests: Conflicts of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Manganese ion homeostasis is vital for bacteria and is achieved via manganese-dependent transcription factors. Manganese mediation of transcription factor attachment to the corresponding oligonucleotide sequences can be investigated, e.g. via electrophoretic mobility shift assays (EMSA). Formation of specific biocomplexes leads to differences in the migration pattern upon gel electrophoresis. Focusing on electrophoresis in the gas-phase, applying a nano electrospray gas-phase electrophoretic mobility molecular analyzer (nES GEMMA) also known as nES differential mobility analyzer (nES DMA), and on transcription factors (MntR proteins) from Bacillus subtilis and Mycobacterium tuberculosis, we took interest in the gas-phase electrophoresis of the corresponding biospecific complexes. We compared nES GEMMA, separating analytes in the nanometer regime (a few to several hundred nm in diameter) in the gas-phase in their native state according to particle size, to EMSA data. Indeed we were able to demonstrate manganese-mediated attachment of MntR to target genomic sequences with both analytical techniques. Despite some inherent pitfalls of the nES GEMMA method like analyte/instrument surface interactions, we were able to detect the target complexes. Moreover, we were able to calculate the molecular weight (MW) of the obtained species by application of a correlation function based on nES GEMMA obtained data. As gas-phase electrophoresis also offers the possibility of offline hyphenation to orthogonal analysis techniques, we are confident that nES GEMMA measurements are not just complementary to EMSA, but will offer the possibility of further in-depth characterization of biocomplexes in the future., (© 2024. The Author(s).)

Pioneering transcription factors (TFs) can drive cell fate changes by binding their DNA motifs in a repressive chromatin environment. Recent structures illustrate emerging rules for nucleosome engagement: TFs distort the nucleosomal DNA to gain access or employ alternative DNA-binding modes with smaller footprints, they preferentially access solvent-exposed motifs near the entry/exit sites, and frequently interact with histones. The extent of TF-histone interactions, in turn, depends on the motif location on the nucleosome, the type of DNA-binding fold, and adjacent domains present. TF-histone interactions can phase TF motifs relative to nucleosomes, and we discuss how these complex and surprisingly diverse interactions between nucleosomes and TFs contribute to function., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Revealing the mechanisms that influence transcription factor binding specificity is the key to understanding gene regulation. In previous studies, DNA double helix structure and one-hot embedding have been used successfully to design computational methods for predicting transcription factor binding sites (TFBSs). However, DNA sequence as a kind of biological language, the method of word embedding representation in natural language processing, has not been considered properly in TFBS prediction models. In our work, we integrate different types of features of DNA sequence to design a multichanneled deep learning framework, namely MulTFBS, in which independent one-hot encoding, word embedding encoding, which can incorporate contextual information and extract the global features of the sequences, and double helix three-dimensional structural features have been trained in different channels. To extract sequence high-level information effectively, in our deep learning framework, we select the spatial-temporal network by combining convolutional neural networks and bidirectional long short-term memory networks with attention mechanism. Compared with six state-of-the-art methods on 66 universal protein-binding microarray data sets of different transcription factors, MulTFBS performs best on all data sets in the regression tasks, with the average R 2 of 0.698 and the average PCC of 0.833, which are 5.4% and 3.2% higher, respectively, than the suboptimal method CRPTS. In addition, we evaluate the classification performance of MulTFBS for distinguishing bound or unbound regions on TF ChIP-seq data. The results show that our framework also performs well in the TFBS classification tasks.

Transcription factors (TFs) are important regulatory elements for vital cellular activities, and the identification of transcription factor binding sites (TFBS) can help to explore gene regulatory mechanisms. Research studies have proved that cfDNA (cell-free DNA) shows relatively higher coverage at TFBS due to the protection by TF from degradation by nucleases and short fragments of cfDNA are enriched in TFBS. However, there are still great difficulties in the noninvasive identification of TFBSs from experimental techniques. In this study, we propose a deep learning-based approach that can noninvasively predict TFBSs of cfDNA by learning sequence information from known TFBSs through convolutional neural networks. Under the addition of long short-term memory, our model achieved an area under the curve of 84%. Based on this model to predict cfDNA, we found consistent motifs in cfDNA fragments and lower coverage occurred upstream and downstream of these cfDNA fragments, which is consistent with a previous study. We also found that the binding sites of the same TF differ in different cell lines. TF-specific target genes were detected from cfDNA and were enriched in cancer-related pathways. In summary, our method of locating TFBSs from plasma has the potential to reflect the intrinsic regulatory mechanism from a noninvasive perspective and provide technical guidance for dynamic monitoring of disease in clinical practice.

Herein, we reported an innovative thermodynamic allosteric switch-actuated 3D DNA nanomachine for selective, sensitive, and accurate electrochemical (EC)/fluorescent (FL) dual-mode biosensing of a microphthalmia-associated transcription factor (MITF). The thermodynamic allosteric switch was ingeniously customized as a hairpin probe (HP) that was in dynamic equilibrium but rapidly interconverting conformations. At the "inactive state", the MITF-binding region and the switch part were "sequestered". Upon the introduction of MITF, an MITF-HP complex promptly formed, and the equilibrium of HP thermodynamically inclined from the "inactive state" toward the "active state" conformation. Immediately, the exposed switch on HP effectively actuated the 3D DNA nanomachine and synchronously produced the restriction site for Nb.BbvCI nicking endonuclease. After the autonomous conveying of the 3D DNA nanomachine by means of the high-efficiency circularly nicking endonuclease signal amplification (NESA), not only was MB-S1 in the supernatant used for FL measurements but also MB-SP/MNs/S2 in the precipitate was adapted for EC analysis, significantly improving the utilization of output products derived from the 3D DNA nanomachine. Accordingly, benefiting from the efficient DNA nanomachine signal amplification manner and the self-calibration function of a dual-mode bioassay, the constructed biosensor exhibits superior sensitivity and accuracy for MITF determination.

Transcription regulators play central roles in orchestrating responses to changing environmental conditions. Recently the Caulobacter crescentus transcription activator DriD, which belongs to the newly defined WYL-domain family, was shown to regulate DNA damage responses independent of the canonical SOS pathway. However, the molecular mechanisms by which DriD and other WYL-regulators sense environmental signals and recognize DNA are not well understood. We showed DriD DNA-binding is triggered by its interaction with ssDNA, which is produced during DNA damage. Here we describe the structure of the full-length C. crescentus DriD bound to both target DNA and effector ssDNA. DriD consists of an N-terminal winged-HTH (wHTH) domain, linker region, three-helix bundle, WYL-domain and C-terminal WCX-dimer domain. Strikingly, DriD binds DNA using a novel, asymmetric DNA-binding mechanism that results from different conformations adopted by the linker. Although the linker does not touch DNA, our data show that contacts it makes with the wHTH are key for specific DNA binding. The structure indicates how ssDNA-effector binding to the WYL-domain impacts wHTH DNA binding. In conclusion, we present the first structure of a WYL-activator bound to both effector and target DNA. The structure unveils a unique, asymmetric DNA binding mode that is likely conserved among WYL-activators., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Transcription factor (TF) proteins regulate gene expression by binding to specific sites on the genome. In the facilitated diffusion model, an optimized search process is achieved by the TF alternating between 3D diffusion in the bulk and 1D diffusion along DNA. While undergoing 1D diffusion, the protein can switch from a search mode for fast diffusion along nonspecific DNA to a recognition mode for stable binding to specific DNA. It was recently noticed that, for a small TF domain protein, reorientations on DNA happen between the nonspecific and specific DNA binding. We here conducted all-atom molecular dynamics simulations with steering forces to reveal the protein-DNA binding free energetics, confirming that the search and recognition modes are distinguished primarily by protein orientations on the DNA. As the binding free energy difference between the specific and nonspecific DNA system slightly deviates from that being estimated directly from dissociation constants on 15-bp DNA constructs, we hypothesize that the discrepancy can come from DNA sequences flanking the 6-bp central binding sites that impact on the dissociation kinetics measurements. The hypothesis is supported by a simplified spherical protein-DNA model along with stochastic simulations and kinetic modeling., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

RNA, unlike DNA, folds into a multitude of secondary and tertiary structures. This structural diversity has impeded the development of ligands that can sequence-specifically target this biomolecule. We sought to develop ligands for double-stranded RNA (dsRNA) segments, which are ubiquitous in RNA tertiary structure. The major groove of double-stranded DNA is sequence-specifically recognized by a range of dimeric helical transcription factors, including the basic leucine zippers (bZIP) and basic helix-loop-helix (bHLH) proteins; however, such simple structural motifs are not prevalent in RNA-binding proteins. We interrogated the high-resolution structures of DNA and RNA to identify requirements for a helix fork motif to occupy dsRNA major grooves akin to dsDNA. Our analysis suggested that the rigidity and angle of approach of dimeric helices in bZIP/bHLH motifs are not ideal for the binding of dsRNA major grooves. This investigation revealed that the replacement of the leucine zipper motifs in bHLH proteins with synthetic crosslinkers would allow recognition of dsRNA. We show that a model bHLH DNA-binding motif does not bind dsRNA but can be reengineered as an RNA ligand. Based on this hypothesis, we rationally designed a miniature synthetic crosslinked helix fork (CHF) as a generalizable proteomimetic scaffold for targeting dsRNA. We evaluated several CHF constructs against a set of RNA and DNA hairpins to probe the specificity of the designed construct. Our studies reveal a new class of proteomimetics as an encodable platform for sequence-specific recognition of dsRNA., (© 2023 Wiley-VCH GmbH.)

LysR-type transcriptional regulators (LTTRs) form one of the largest families of bacterial regulators. They are widely distributed and contribute to all aspects of metabolism and physiology. Most are homotetramers, with each subunit composed of an N-terminal DNA-binding domain followed by a long helix connecting to an effector-binding domain. LTTRs typically bind DNA in the presence or absence of a small-molecule ligand (effector). In response to cellular signals, conformational changes alter DNA interactions, contact with RNA polymerase, and sometimes contact with other proteins. Many are dual-function repressor-activators, although different modes of regulation may occur at multiple promoters. This review presents an update on the molecular basis of regulation, the complexity of regulatory schemes, and applications in biotechnology and medicine. The abundance of LTTRs reflects their versatility and importance. While a single regulatory model cannot describe all family members, a comparison of similarities and differences provides a framework for future study.

Transcription factors (TFs) are proteins that affect gene expression by binding to regulatory regions of DNA in a sequence specific manner. The binding of TFs to DNA is controlled by many factors, including the DNA sequence, concentration of TF, chromatin accessibility and co-factors. Here, we systematically investigated the binding mechanism of hundreds of TFs by analysing ChIP-seq data with our explainable statistical model, ChIPanalyser. This tool uses as inputs the DNA sequence binding motif; the capacity to distinguish between strong and weak binding sites; the concentration of TF; and chromatin accessibility. We found that approximately one third of TFs are predicted to bind the genome in a DNA accessibility independent fashion, which includes TFs that can open the chromatin, their co-factors and TFs with similar motifs. Our model predicted this to be the case when the TF binds to its strongest binding regions in the genome, and only a small number of TFs have the capacity to bind dense chromatin at their weakest binding regions, such as CTCF, USF2 and CEBPB. Our study demonstrated that the binding of hundreds of human and mouse TFs is predicted by ChIPanalyser with high accuracy and showed that many TFs can bind dense chromatin., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Associations of transcription factors (TFs) with specific sites on DNA initiate major cellular processes. But DNA in eukaryotic cells is covered by nucleosomes which prevent TFs from binding. However, nucleosome structures on DNA are not static and exhibit breathing and sliding. We develop a theoretical framework to investigate the effect of nucleosome sliding on a protein target search. By analysis of a discrete-state stochastic model of nucleosome sliding, search dynamics are explicitly evaluated. It is found that for long sliding lengths the target search dynamics are faster for normal TFs that cannot enter the nucleosomal DNA. But for more realistic short sliding lengths, the so-called pioneer TFs, which can invade nucleosomal DNA, locate specific sites faster. It is also suggested that nucleosome breathing, which is a faster process, has a stronger effect on protein search dynamics than that of nucleosome sliding. Theoretical arguments to explain these observations are presented.

Pioneer transcription factors (PTFs) have the remarkable ability to directly bind to chromatin to stimulate vital cellular processes. In this work, we dissect the universal binding mode of Sox PTF by combining extensive molecular simulations and physiochemistry approaches, along with DNA footprinting techniques. As a result, we show that when Sox consensus DNA is located at the solvent-facing DNA strand, Sox binds to the compact nucleosome without imposing any significant conformational changes. We also reveal that the base-specific Sox:DNA interactions (base reading) and Sox-induced DNA changes (shape reading) are concurrently required for sequence-specific nucleosomal DNA recognition. Among three different nucleosome positions located on the positive DNA arm, a sequence-specific reading mechanism is solely satisfied at the superhelical location 2 (SHL2). While SHL2 acts transparently for solvent-facing Sox binding, among the other two positions, SHL4 permits only shape reading. The final position, SHL0 (dyad), on the other hand, allows no reading mechanism. These findings demonstrate that Sox-based nucleosome recognition is essentially guided by intrinsic nucleosome properties, permitting varying degrees of DNA recognition.

The human genome contains about 800 C2H2 zinc finger proteins (ZFPs), and most of them are composed of long arrays of zinc fingers. Standard ZFP recognition model asserts longer finger arrays should recognize longer DNA-binding sites. However, recent experimental efforts to identify in vivo ZFP binding sites contradict this assumption, with many exhibiting short motifs. Here we use ZFY, CTCF, ZIM3, and ZNF343 as examples to address three closely related questions: What are the reasons that impede current motif discovery methods? What are the functions of those seemingly unused fingers and how can we improve the motif discovery algorithms based on long ZFPs' biophysical properties? Using ZFY, we employed a variety of methods and find evidence for 'dependent recognition' where downstream fingers can recognize some previously undiscovered motifs only in the presence of an intact core site. For CTCF, high-throughput measurements revealed its upstream specificity profile depends on the strength of its core. Moreover, the binding strength of the upstream site modulates CTCF's sensitivity to different epigenetic modifications within the core, providing new insight into how the previously identified intellectual disability-causing and cancer-related mutant R567W disrupts upstream recognition and deregulates the epigenetic control by CTCF. Our results establish that, because of irregular motif structures, variable spacing and dependent recognition between sub-motifs, the specificities of long ZFPs are significantly underestimated, so we developed an algorithm, ModeMap, to infer the motifs and recognition models of ZIM3 and ZNF343, which facilitates high-confidence identification of specific binding sites, including repeats-derived elements. With revised concept, technique, and algorithm, we can discover the overlooked specificities and functions of those 'extra' fingers, and therefore decipher their broader roles in human biology and diseases., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Precise dissection of DNA-protein interactions is essential for elucidating the recognition basis, dynamics and gene regulation mechanism. However, global profiling of weak and dynamic DNA-protein interactions remains a long-standing challenge. Here, we establish the light-induced lysine (K) enabled crosslinking (LIKE-XL) strategy for spatiotemporal and global profiling of DNA-protein interactions. Harnessing unique abilities to capture weak and transient DNA-protein interactions, we demonstrate that LIKE-XL enables the discovery of low-affinity transcription-factor/DNA interactions via sequence-specific DNA baits, determining the binding sites for transcription factors that have been previously unknown. More importantly, we successfully decipher the dynamics of the transcription factor subproteome in response to drug treatment in a time-resolved manner, and find downstream target transcription factors from drug perturbations, providing insight into their dynamic transcriptional networks. The LIKE-XL strategy offers a complementary method to expand the DNA-protein profiling toolbox and map accurate DNA-protein interactomes that were previously inaccessible via non-covalent strategies, for better understanding of protein function in health and disease., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

The BEN domain-containing transcription factors regulate transcription by recruiting chromatin-modifying factors to specific chromatin regions via their DNA-binding BEN domains. The BEN domain of BANP has been shown to bind to a CGCG DNA sequence or an AAA-containing matrix attachment regions DNA sequence. Consistent with these in vivo observations, we identified an optimal DNA-binding sequence of AAATCTCG by protein binding microarray, which was also confirmed by our isothermal titration calorimetry and mutagenesis results. We then determined crystal structures of the BANP BEN domain in apo form and in complex with a CGCG-containing DNA, respectively, which revealed that the BANP BEN domain mainly used the electrostatic interactions to bind DNA with some base-specific interactions with the TC motifs. Our isothermal titration calorimetry results also showed that BANP bound to unmethylated and methylated DNAs with comparable binding affinities. Our complex structure of BANP-mCGCG revealed that the BANP BEN domain bound to the unmethylated and methylated DNAs in a similar mode and cytosine methylation did not get involved in binding, which is also consistent with our observations from the complex structures of the BEND6 BEN domain with the CGCG or CGmCG DNAs. Taken together, our results further elucidate the elements important for DNA recognition and transcriptional regulation by the BANP BEN domain-containing transcription factor., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Bacterial transcription factors (TFs) with helix-turn-helix (HTH) DNA-binding domains have been widely explored to build orthogonal transcriptional regulation systems in mammalian cells. Here we capitalize on the modular structure of these proteins to build a framework for multi-input logic gates relying on serial combinations of inducible protein-protein interactions. We found that for some TFs, their HTH domain alone is sufficient for DNA binding. By fusing the HTH domain to TFs, we established dimerization dependent rather than DNA-binding-dependent activation. This enabled us to convert gene switches from OFF-type into more widely applicable ON-type systems and to create mammalian gene switches responsive to new inducers. By combining both OFF and ON modes of action, we built a compact, high-performance bandpass filter. Furthermore, we were able to show cytosolic and extracellular dimerization. Cascading up to five pairwise fusion proteins yielded robust multi-input AND logic gates. Combinations of different pairwise fusion proteins afforded a variety of 4-input 1-output AND and OR logic gate configurations., (© 2023. The Author(s).)

Transcription factors are multidomain proteins with specific DNA binding and regulatory domains. In the human FoxP subfamily (FoxP1, FoxP2, FoxP3, and FoxP4) of transcription factors, a 90 residue-long disordered region links a Leucine Zipper (ZIP)-known to form coiled-coil dimers-and a Forkhead (FKH) domain-known to form domain swapping dimers. We used replica exchange discrete molecular dynamics simulations, single-molecule fluorescence experiments, and other biophysical tools to understand how domain tethering in FoxP1 impacts dimerization at ZIP and FKH domains and how DNA binding allosterically regulates their dimerization. We found that domain tethering promotes FoxP1 dimerization but inhibits a FKH domain-swapped structure. Furthermore, our findings indicate that the linker mediates the mutual organization and dynamics of ZIP and FKH domains, forming closed and open states with and without interdomain contacts, thus highlighting the role of the linkers in multidomain proteins. Finally, we found that DNA allosterically promotes structural changes that decrease the dimerization propensity of FoxP1. We postulate that, upon DNA binding, the interdomain linker plays a crucial role in the gene regulatory function of FoxP1., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

Plants use sophisticated mechanisms of gene expression to control senescence in response to environmental stress or aging. ORE1 (Arabidopsis thaliana NAC092) is a master regulator of senescence that belongs to the plant-specific NAC transcription factor protein family. ORE1 has been reported to bind to multiple DNA targets to orchestrate leaf senescence, yet the mechanistic basis for recognition of the cognate gene sequence remains unclear. Here, we report the crystal structure of the ORE1-NAC domain alone and its DNA-binding form. The structure of DNA-bound ORE1-NAC revealed the molecular basis for nucleobase recognition and phosphate backbone interactions. We show that local versatility in the DNA-binding site, in combination with domain flexibility of the ORE-NAC homodimer, is crucial for the maintenance of binding to intrinsically flexible DNA. Our results provide a platform for understanding other plant-specific NAC protein-DNA interactions as well as insight into the structural basis of NAC regulators in plants of agronomic and scientific importance., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Motivation: Transcription factor (TF) binds to conservative DNA binding sites in different cellular environments and development stages by physical interaction with interdependent nucleotides. However, systematic computational characterization of the relationship between higher-order nucleotide dependency and TF-DNA binding mechanism in diverse cell types remains challenging., Results: Here, we propose a novel multi-task learning framework HAMPLE to simultaneously predict TF binding sites (TFBS) in distinct cell types by characterizing higher-order nucleotide dependencies. Specifically, HAMPLE first represents a DNA sequence through three higher-order nucleotide dependencies, including k-mer encoding, DNA shape and histone modification. Then, HAMPLE uses the customized gate control and the channel attention convolutional architecture to further capture cell-type-specific and cell-type-shared DNA binding motifs and epigenomic languages. Finally, HAMPLE exploits the joint loss function to optimize the TFBS prediction for different cell types in an end-to-end manner. Extensive experimental results on seven datasets demonstrate that HAMPLE significantly outperforms the state-of-the-art approaches in terms of auROC. In addition, feature importance analysis illustrates that k-mer encoding, DNA shape, and histone modification have predictive power for TF-DNA binding in different cellular environments and are complementary to each other. Furthermore, ablation study, and interpretable analysis validate the effectiveness of the customized gate control and the channel attention convolutional architecture in characterizing higher-order nucleotide dependencies., Availability and Implementation: The source code is available at https://github.com/ZhangLab312/Hample., (© The Author(s) 2023. Published by Oxford University Press.)

The association of multiple Transcription Factors (TFs) in the cis-regulatory region is imperative for developmental changes in eukaryotes. The underlying process is exceedingly complex, and it is not at all clear what orchestrates the overall search process by multiple TFs. In this study, by developing a theoretical model based on a discrete-state stochastic approach, we investigated the target search mechanism of multiple TFs on nucleosomal DNA. Experimental kinetic rate constants of different TFs are taken as input to estimate the Mean-First-Passage time to recognize the binding motifs by two TFs on a dynamic nucleosome model. The theory systematically analyzes when the TFs search their binding motifs hierarchically and when simultaneously by proceeding via the formation of a protein-protein complex. Our results, validated by extensive Monte Carlo simulations, elucidate the molecular basis of the complex target search phenomenon of multiple TFs on nucleosomal DNA., (© 2023 Wiley-VCH GmbH.)

Methyl CpG binding proteins (MBPs) are transcription factors that recognize the methylated CpG sites in DNA and mediate the DNA methylation signal into various downstream cellular processes. The C2H2 zinc finger (ZF) protein, Kaiso, also an MBP, preferentially binds to two symmetrically methylated CpG sites in DNA sequences via C-terminal C2H2 ZF domains and mediates the transcription regulation process. Investigation of the molecular mechanism of the recognition of methylated DNA (meDNA) by Kaiso is important to understand how this protein reads and translates this methylation signal into downstream transcription outcomes. Despite previous studies in Kaiso-meDNA interactions, detailed structural investigations on the sequence-specific interaction of Kaiso with the meDNA sequence are still lacking. In this work, we used molecular modeling and molecular dynamics (MD) simulation-based computational approaches to investigate the recognition of various methylated DNA sequences by Kaiso. Our MD simulation results show that the Kaiso-meDNA interaction is sequence specific. The recognition of meDNA by Kaiso is enhanced in the MeECad sequence compared to the MeCG2 sequence. Compared to the 5'-flanking T/A pair in MeCG2, both MeCG2_mutCG and MeECad sequences show that a C/G base pair allows GLU535 of Kaiso to preferably recognize and bind the core mCpG site. The core mCGmCG site is crucial for the recognition process and formation of a stable complex. Our results reveal that the 5'-flanking nucleotides are also important for the enhanced binding and recognition of methylated sites.

Transcription factors (TF) require access to target sites within nucleosomes to initiate transcription. The target site position within the nucleosome significantly influences TF occupancy, but how is not quantitatively understood. Using ensemble and single-molecule fluorescence measurements, we investigated the targeting and occupancy of the transcription factor, Gal4, at different positions within the nucleosome. We observe a dramatic decrease in TF occupancy to sites extending past 30 base pairs (bp) into the nucleosome which cannot be explained by changes in the TF dissociation rate or binding site orientation. Instead, the nucleosome unwrapping free energy landscape is the primary determinant of Gal4 occupancy by reducing the Gal4 binding rate. The unwrapping free energy landscape defines two distinct regions of accessibility and kinetics with a boundary at 30 bp into the nucleosome where the inner region is over 100-fold less accessible. The Gal4 binding rate in the inner region no longer depends on its concentration because it is limited by the nucleosome unwrapping rate, while the frequency of nucleosome rewrapping decreases because Gal4 exchanges multiple times before the nucleosome rewraps. Our findings highlight the importance of the nucleosome unwrapping free energy landscape on TF occupancy and dynamics that ultimately influences transcription initiation., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Plant-specific TCP transcription factors are key regulators of diverse plant functions. TCP transcription factors have long been annotated as basic helix-loop-helix (bHLH) transcription factors according to remote sequence homology without experimental validation, and their consensus DNA-binding sequences and protein-DNA recognition mechanisms have remained elusive. Here, we report the crystal structures of the class I TCP domain from AtTCP15 and the class II TCP domain from AtTCP10 in complex with different double-stranded DNA (dsDNA). The complex structures reveal that the TCP domain is a distinct DNA-binding motif and the homodimeric TCP domains adopt a unique three-site recognition mode, binding to dsDNA mainly through a central pair of β-strands formed by the dimer interface and two basic flexible loops from each monomer. The consensus DNA-binding sequence for class I TCPs is a perfectly palindromic 11 bp (GTGGGNCCCAC), whereas that for class II TCPs is a near-palindromic 11 bp (GTGGTCCCCAC). The unique DNA binding mode allows the TCP domains to display broad specificity for a range of DNA sequences even shorter than 11 bp, adding further complexity to the regulatory network of plant TCP transcription factors., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Kaiso is a methyl-DNA-binding protein containing three C2H2 zinc fingers with a C-terminal extension that participates in DNA binding. The linker between the last zinc finger and the DNA-binding portion of the extension contains two prolines that are highly conserved in vertebrates and in cognate ZBTB4 and ZBTB38 proteins. Prolines provide chain rigidity and can exist in cis and trans conformations that can be switched by proline isomerases, affecting protein function. We found that substitution of the conserved proline P588, but not of P577, to alanine, negatively affected KaisoDNA-binding according to molecular dynamics simulation and in vitro DNA-binding assays. Molecular dynamics simulations of the Kaiso DNA-binding domain with P588 either substituted to alanine or switched to the cis-conformation revealed similar alterations in the H-bonding network and uncovered allosteric effects leading to structural rearrangements in the entire domain that resulted in the weakening of DNA-binding affinity. The substitution of proline with a large hydrophobic residue led to the same negative effects despite its ability to partially rescue the intrinsic DNA-binding activity of the C-terminal loop. Thus, the presence of the C-terminal extension and cis-conformation of proline residues are essential for efficient Kaiso-DNA binding, which likely involves intramolecular tension squeezing the DNA chain.

Sequence-specific recognition of transcription factor (TF) binding motifs in the target site of DNA over the vast amount of non-target DNA is of primary importance for the transcriptional regulation of gene expression by the TFs. Binding of TFs to the target site of DNA relies not only on the direct contact formation but also on the structural and conformational features of DNA. Recognition of DNA structural features or shape readout by proteins is an important factor in the context of TF-DNA interaction. Based on the atomistic molecular simulation, here we report the sequence-dependent unique structural features, solvation, and ion-binding properties of biologically relevant AT- and GC-rich human TF binding motifs in DNA. Counterion and water distribution around the motif is found to be sensitive to the motif sequence, which is accompanied with the DNA shape features. The motif sequence affects the electrostatic potential along the grooves, and cytosine methylation alters the DNA shape features. Characteristic solvation properties of TF binding motif DNA fragments infer that an ionic environment and hydration influences are essential to describe TF-DNA interactions.

Some transcription factors bind DNA motifs containing direct or inverted sequence repeats. Preference for each of these DNA topologies is dictated by structural constraints. Most prokaryotic regulators form symmetric oligomers, which require operators with a dyad structure. Binding to direct repeats requires breaking the internal symmetry, a property restricted to a few regulators, most of them from the AraC family. The KorA family of transcriptional repressors, involved in plasmid propagation and stability, includes members that form symmetric dimers and recognize inverted repeats. Our structural analyses show that ArdK, a member of this family, can form a symmetric dimer similar to that observed for KorA, yet it binds direct sequence repeats as a non-symmetric dimer. This is possible by the 180° rotation of one of the helix-turn-helix domains. We then probed and confirmed that ArdK shows affinity for an inverted repeat, which, surprisingly, is also recognized by a non-symmetrical dimer. Our results indicate that structural flexibility at different positions in the dimerization interface constrains transcription factors to bind DNA sequences with one of these two alternative DNA topologies., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Discovery of transcription factor binding sites (TFBSs) is of primary importance for understanding the underlying binding mechanic and gene regulation process. Growing evidence indicates that apart from the primary DNA sequences, DNA shape landscape has a significant influence on transcription factor binding preference. To effectively model the co-influence of sequence and shape features, we emphasize the importance of position information of sequence motif and shape pattern. In this paper, we propose a novel deep learning-based architecture, named hybridShape eDeepCNN, for TFBS prediction which integrates DNA sequence and shape information in a spatially aligned manner. Our model utilizes the power of the multi-layer convolutional neural network and constructs an independent subnetwork to adapt for the distinct data distribution of heterogeneous features. Besides, we explore the usage of continuous embedding vectors as the representation of DNA sequences. Based on the experiments on 20 in-vitro datasets derived from universal protein binding microarrays (uPBMs), we demonstrate the superiority of our proposed method and validate the underlying design logic.

The regulation of transcription is a complex process that involves binding of transcription factors (TFs) to specific sequences, recruitment of cofactors and chromatin remodelers, assembly of the pre-initiation complex and recruitment of RNA polymerase II. Increasing evidence suggests that TFs are highly dynamic and interact only transiently with DNA. Single molecule microscopy techniques are powerful approaches for tracking individual TF molecules as they diffuse in the nucleus and interact with DNA. Here we employ multifocus microscopy and highly inclined laminated optical sheet microscopy to track TF dynamics in response to perturbations in labile zinc inside cells. We sought to define whether zinc-dependent TFs sense changes in the labile zinc pool by determining whether their dynamics and DNA binding can be modulated by zinc. We used fluorescently tagged versions of the glucocorticoid receptor (GR), with two C4 zinc finger domains, and CCCTC-binding factor (CTCF), with eleven C2H2 zinc finger domains. We found that GR was largely insensitive to perturbations of zinc, whereas CTCF was significantly affected by zinc depletion and its dwell time was affected by zinc elevation. These results indicate that at least some transcription factors are sensitive to zinc dynamics, revealing a potential new layer of transcriptional regulation., (© 2022. The Author(s).)

A technique for detecting the activation of cell signal transduction is particularly important for disease diagnosis and therapy. Transcriptional factor (TF) activities that could indicate the status of cell signal transduction are a favorable target for cell signal detection. Tetrahedron DNA nanostructures (TDNs) which contain specific binding sequences of TFs were designed and synthesized in this research, and their effects on detecting cell signal transduction were evaluated. We found that FAM-labeled TDNs with the indicated TF binding sequences could specifically bind to activated TFs of hypoxia signaling or TGF-β signaling. Signaling pathway activities detected via TDNs could be exhibited by various methods including fluorescence imaging, flow cytometry and fluorescence spectrometer analysis. The reliability of this new technique is in line with the classical dual luciferase reporter assay system. This work develops a novel and effective tool to examine the activation of intracellular signaling pathways via nanotechnology. In addition, good stability and programmability of TDNs ensure their widespread application in various signaling pathways.

Selective gene regulation is mediated by recognition of specific DNA sequences by transcription factors (TFs). The extremely challenging task of searching out specific cognate DNA binding sites among several million putative sites within the eukaryotic genome is achieved by complex molecular recognition mechanisms. Elements of this recognition code include the core binding sequence, the flanking sequence context, and the shape and conformational flexibility of the composite binding site. To unravel the extent to which DNA flexibility modulates TF binding, in this study, we employed experimentally guided molecular dynamics simulations of ternary complex of closely related Hox heterodimers Exd-Ubx and Exd-Scr with DNA. Results demonstrate that flexibility signatures embedded in the flanking sequences impact TF binding at the cognate binding site. A DNA sequence has intrinsic shape and flexibility features. While shape features are localized, our analyses reveal that flexibility features of the flanking sequences percolate several basepairs and allosterically modulate TF binding at the core. We also show that lack of flexibility in the motif context can render the cognate site resistant to protein-induced shape changes and subsequently lower TF binding affinity. Overall, this study suggests that flexibility-guided DNA shape, and not merely the static shape, is a key unexplored component of the complex DNA-TF recognition code., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Knowledge of how proteins interact with DNA is essential for understanding gene regulation. Although DNA-binding specificities for thousands of transcription factors (TFs) have been determined, the specific amino acid-base interactions comprising their structural interfaces are largely unknown. This lack of resolution hampers attempts to leverage these data in order to predict specificities for uncharacterized TFs or TFs mutated in disease. Here we introduce recognition code learning via automated mapping of protein-DNA structural interfaces (rCLAMPS), a probabilistic approach that uses DNA-binding specificities for TFs from the same structural family to simultaneously infer both which nucleotide positions are contacted by particular amino acids within the TF as well as a recognition code that relates each base-contacting amino acid to nucleotide preferences at the DNA positions it contacts. We apply rCLAMPS to homeodomains, the second largest family of TFs in metazoans and show that it learns a highly effective recognition code that can predict de novo DNA-binding specificities for TFs. Furthermore, we show that the inferred amino acid-nucleotide contacts reveal whether and how nucleotide preferences at individual binding site positions are altered by mutations within TFs. Our approach is an important step toward automatically uncovering the determinants of protein-DNA specificity from large compendia of DNA-binding specificities and inferring the altered functionalities of TFs mutated in disease., (© 2022 Wetzel et al.; Published by Cold Spring Harbor Laboratory Press.)

We present a physics-based machine learning approach to predict in vitro transcription factor binding affinities from structural and mechanical DNA properties directly derived from atomistic molecular dynamics simulations. The method is able to predict affinities obtained with techniques as different as uPBM, gcPBM and HT-SELEX with an excellent performance, much better than existing algorithms. Due to its nature, the method can be extended to epigenetic variants, mismatches, mutations, or any non-coding nucleobases. When complemented with chromatin structure information, our in vitro trained method provides also good estimates of in vivo binding sites in yeast., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Dimerization is critical for transcription factors (TFs) to bind DNA and regulate a wide variety of cellular functions; however, the molecular mechanisms remain to be completely elucidated. Here, we used high-speed atomic force microscopy (HS-AFM) to observe the dimerization process for a photoresponsive TF Photozipper (PZ), which consists of light-oxygen-voltage-sensing (LOV) and basic-region-leucine-zipper (bZIP) domains. HS-AFM visualized not only the oligomeric states of PZ molecules forming monomers and dimers under controlled dark-light conditions but also the domain structures within each molecule. Successive AFM movies captured the dimerization process for an individual PZ molecule and the monomer-dimer reversible transition during dark-light cycling. Detailed AFM images of domain structures in PZ molecules demonstrated that the bZIP domain entangled under dark conditions was loosened owing to light illumination and fluctuated around the LOV domain. These observations revealed the role of the bZIP domain in the dimerization processes of a TF., (© 2022. The Author(s).)

A 22 kDa protein from Thermus thermophilus is characterised as a DNA binding transcription regulator and its function is established using the fluorescence spectroscopy technique. The steady-state fluorescence spectroscopy result shows significant binding of calf thymus DNA and protein molecule. To confirm, the DNA quenching effect in real-time, a time-resolved emission spectroscopy study was performed and the result shows good agreement with steady-state quenching analysis., (© 2022 Wiley Periodicals LLC.)

All major biological processes start after protein molecules known as transcription factors detect specific regulatory sequences on DNA and initiate genetic expression by associating to them. But in eukaryotic cells, much of the DNA is covered by nucleosomes and other chromatin structures, preventing transcription factors from binding to their targets. At the same time, experimental studies show that there are several classes of proteins, called "pioneer transcription factors", that are able to reach the targets on nucleosomal DNA; however, the underlying microscopic mechanisms remain not well understood. We propose a new theoretical approach that might explain how pioneer transcription factors can find their targets. It is argued that pioneer transcription factors might weaken the interactions between the DNA and nucleosome by substituting them with similar interactions between transcription factors and DNA. Using this idea, we develop a discrete-state stochastic model that allows for exact calculations of target search dynamics on nucleosomal DNA using first-passage probabilities approach. It is found that the target search on nuclesomal DNA for pioneer transcription factors might be significantly accelerated while the search is slower on naked DNA in comparison with normal transcription factors. Our theoretical predictions are supported by Monte Carlo computer simulations, and they also agree with available experimental observations.

Transcription, the process of copying genetic information from DNA to messenger RNA, is regulated by sequence-specific DNA-binding proteins known as transcription factors (TFs). Recent advances in single-molecule tracking (SMT) technologies have enabled visualization of individual TF molecules as they diffuse and interact with the DNA in the context of living cells. These SMT studies have uncovered multiple populations of DNA-binding events characterized by their distinctive DNA residence times. In this perspective, we review recent insights into how these residence times relate to specific and non-specific DNA binding, as well as the contribution of TF domains on the DNA-binding dynamics. We discuss different models that aim to link transient DNA binding by TFs to bursts of transcription and present an outlook for how future advances in microscopy development may broaden our understanding of the dynamics of the molecular steps that underlie transcription activation., (Copyright © 2022 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Single-molecule fluorescence detection offers powerful ways to study biomolecules and their complex interactions. Here, nanofluidic devices and camera-based, single-molecule Förster resonance energy transfer (smFRET) detection are combined to study the interactions between plant transcription factors of the auxin response factor (ARF) family and DNA oligonucleotides that contain target DNA response elements. In particular, it is shown that the binding of the unlabeled ARF DNA binding domain (ARF-DBD) to donor and acceptor labeled DNA oligonucleotides can be detected by changes in the FRET efficiency and changes in the diffusion coefficient of the DNA. In addition, this data on fluorescently labeled ARF-DBDs suggest that, at nanomolar concentrations, ARF-DBDs are exclusively present as monomers. In general, the fluidic framework of freely diffusing molecules minimizes potential surface-induced artifacts, enables high-throughput measurements, and proved to be instrumental in shedding more light on the interactions between ARF-DBDs monomers and between ARF-DBDs and their DNA response element., (© 2021 The Authors. Advanced Biology published by Wiley-VCH GmbH.)

One-dimensional (1-D) sliding of transcription factor (TF) protein along DNA is essential for facilitated diffusion of the TF to locate target DNA site for genetic regulation. Detecting base-pair (bp) resolution of the TF sliding or stepping on the DNA is still experimentally challenging. We have recently performed all-atom molecular dynamics (MD) simulations capturing spontaneous 1-bp stepping of a small WRKY domain TF protein along DNA. Based on the 10 µs WRKY stepping path obtained from such simulations, the protocol here shows how to conduct more extensive conformational samplings of the TF-DNA systems, by constructing the Markov state model (MSM) for the 1-bp protein stepping, with various numbers of micro- and macro-states tested for the MSM construction. In order to examine processive 1-D diffusional search of the TF protein along DNA with structural basis, the protocol further shows how to conduct coarse-grained (CG) MD simulations to sample long-time scale dynamics of the system. Such CG modeling and simulations are particularly useful to reveal the protein-DNA electrostatic impacts on the processive diffusional motions of the TF protein above tens of microseconds, in comparison with sub-microseconds to microseconds protein stepping motions revealed from the all-atom simulations.

Biolayer interferometry (BLI) is a widely utilized technique for determining macromolecular interaction dynamics in real time. Using changes in the interference pattern of white light reflected off a biosensor tip, BLI can determine binding parameters for protein-protein (e.g., antibody-substrate kinetics) or protein-small molecule (e.g., drug discovery) interactions. However, a less-appreciated application for BLI analysis is DNA-protein interactions. DNA-binding proteins play an immense role in cellular biology, controlling critical processes including transcription, DNA replication, and DNA repair. Understanding how proteins interact with DNA often provides important insight into their biological function, and novel technologies to assay DNA-protein interactions are of broad interest. Currently, a detailed protocol utilizing BLI for DNA-protein interactions is lacking. In the following protocol, we describe the use of BLI and biotinylated-DNA probes to determine the binding kinetics of a transcription factor to a specific DNA sequence. The experimental steps include the generation of biotinylated-DNA probes, the execution of the BLI experiment, and data analysis by scientific graphing and statistical software (e.g., GraphPad Prism). Although the example experiment used throughout this protocol involves a prokaryotic transcription factor, this technique can be easily translated to any DNA-binding protein. Pitfalls and potential solutions for investigating DNA-binding proteins by BLI are also presented., Competing Interests: The authors have declared that no competing interests exist.

The discovery of putative transcription factor binding sites (TFBSs) is important for understanding the underlying binding mechanism and cellular functions. Recently, many computational methods have been proposed to jointly account for DNA sequence and shape properties in TFBSs prediction. However, these methods fail to fully utilize the latent features derived from both sequence and shape profiles and have limitation in interpretability and knowledge discovery. To this end, we present a novel Deep Convolution Attention network combining Sequence and Shape, dubbed as D-SSCA, for precisely predicting putative TFBSs. Experiments conducted on 165 ENCODE ChIP-seq datasets reveal that D-SSCA significantly outperforms several state-of-the-art methods in predicting TFBSs, and justify the utility of channel attention module for feature refinements. Besides, the thorough analysis about the contribution of five shapes to TFBSs prediction demonstrates that shape features can improve the predictive power for transcription factors-DNA binding. Furthermore, D-SSCA can realize the cross-cell line prediction of TFBSs, indicating the occupancy of common interplay patterns concerning both sequence and shape across various cell lines. The source code of D-SSCA can be found at https://github.com/MoonLord0525/., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)