Genome-wide prediction of minor-groove electrostatic potential enables biophysical modeling of protein-DNA binding

Systematic dissection of sequence features affecting binding specificity of a pioneer factor reveals binding synergy between FOXA1 and AP-1

Despite the unique ability of pioneer factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called chromatin immunoprecipitation with integrated synthetic oligonucleotides (ChIP-ISO) to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1, in human A549 cells. Combining ChIP-ISO with in vitro and neural network analyses, we find that (1) FOXA1 binding is strongly affected by co-binding transcription factors (TFs) AP-1 and CEBPB; (2) FOXA1 and AP-1 show binding cooperativity in vitro; (3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin; and (4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

Deep DNAshape webserver: prediction and real-time visualization of DNA shape considering extended k-mers

Sequence-dependent DNA shape plays an important role in understanding protein-DNA binding mechanisms. High-throughput prediction of DNA shape features has become a valuable tool in the field of protein-DNA recognition, transcription factor-DNA binding specificity, and gene regulation. However, our widely used webserver, DNAshape, relies on statistically summarized pentamer query tables to query DNA shape features. These query tables do not consider flanking regions longer than two base pairs, and acquiring a query table for hexamers or higher-order k-mers is currently still unrealistic due to limitations in achieving sufficient statistical coverage in molecular simulations or structural biology experiments. A recent deep-learning method, Deep DNAshape, can predict DNA shape features at the core of a DNA fragment considering flanking regions of up to seven base pairs, trained on limited simulation data. However, Deep DNAshape is rather complicated to install, and it must run locally compared to the pentamer-based DNAshape webserver, creating a barrier for users. Here, we present the Deep DNAshape webserver, which has the benefits of both methods while being accurate, fast, and accessible to all users. Additional improvements of the webserver include the detection of user input in real time, the ability of interactive visualization tools and different modes of analyses. URL: https://deepdnashape.usc.edu.

DNA sequence and chromatin differentiate sequence-specific transcription factor binding in the human malaria parasite Plasmodium falciparum

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.

Discovering DNA shape motifs with multiple DNA shape features: generalization, methods, and validation

DNA motifs are crucial patterns in gene regulation. DNA-binding proteins (DBPs), including transcription factors, can bind to specific DNA motifs to regulate gene expression and other cellular activities. Past studies suggest that DNA shape features could be subtly involved in DNA-DBP interactions. Therefore, the shape motif annotations based on intrinsic DNA topology can deepen the understanding of DNA-DBP binding. Nevertheless, high-throughput tools for DNA shape motif discovery that incorporate multiple features altogether remain insufficient. To address it, we propose a series of methods to discover non-redundant DNA shape motifs with the generalization to multiple motifs in multiple shape features. Specifically, an existing Gibbs sampling method is generalized to multiple DNA motif discovery with multiple shape features. Meanwhile, an expectation-maximization (EM) method and a hybrid method coupling EM with Gibbs sampling are proposed and developed with promising performance, convergence capability, and efficiency. The discovered DNA shape motif instances reveal insights into low-signal ChIP-seq peak summits, complementing the existing sequence motif discovery works. Additionally, our modelling captures the potential interplays across multiple DNA shape features. We provide a valuable platform of tools for DNA shape motif discovery. An R package is built for open accessibility and long-lasting impact: https://zenodo.org/doi/10.5281/zenodo.10558980.

DrivR-Base: a feature extraction toolkit for variant effect prediction model construction

MOTIVATION

Recent advancements in sequencing technologies have led to the discovery of numerous variants in the human genome. However, understanding their precise roles in diseases remains challenging due to their complex functional mechanisms. Various methodologies have emerged to predict the pathogenic significance of these genetic variants. Typically, these methods employ an integrative approach, leveraging diverse data sources that provide important insights into genomic function. Despite the abundance of publicly available data sources and databases, the process of navigating, extracting, and pre-processing features for machine learning models can be highly challenging and time-consuming. Furthermore, researchers often invest substantial effort in feature extraction, only to later discover that these features lack informativeness.

RESULTS

In this article, we introduce DrivR-Base, an innovative resource that efficiently extracts and integrates molecular information (features) related to single nucleotide variants. These features encompass information about the genomic positions and the associated protein positions of a variant. They are derived from a wide array of databases and tools, including structural properties obtained from AlphaFold, regulatory information sourced from ENCODE, and predicted variant consequences from Variant Effect Predictor. DrivR-Base is easily deployable via a Docker container to ensure reproducibility and ease of access across diverse computational environments. The resulting features can be used as input for machine learning models designed to predict the pathogenic impact of human genome variants in disease. Moreover, these feature sets have applications beyond this, including haploinsufficiency prediction and the development of drug repurposing tools. We describe the resource's development, practical applications, and potential for future expansion and enhancement.

AVAILABILITY AND IMPLEMENTATION

DrivR-Base source code is available at https://github.com/amyfrancis97/DrivR-Base.

Building an ab initio solvated DNA model using Euclidean neural networks

Accurately modeling large biomolecules such as DNA from first principles is fundamentally challenging due to the steep computational scaling of ab initio quantum chemistry methods. This limitation becomes even more prominent when modeling biomolecules in solution due to the need to include large numbers of solvent molecules. We present a machine-learned electron density model based on a Euclidean neural network framework that includes a built-in understanding of equivariance to model explicitly solvated double-stranded DNA. By training the machine learning model using molecular fragments that sample the key DNA and solvent interactions, we show that the model predicts electron densities of arbitrary systems of solvated DNA accurately, resolves polarization effects that are neglected by classical force fields, and captures the physics of the DNA-solvent interaction at the ab initio level.

Predicting DNA structure using a deep learning method

Understanding the mechanisms of protein-DNA binding is critical in comprehending gene regulation. Three-dimensional DNA structure, also described as DNA shape, plays a key role in these mechanisms. In this study, we present a deep learning-based method, Deep DNAshape, that fundamentally changes the current k-mer based high-throughput prediction of DNA shape features by accurately accounting for the influence of extended flanking regions, without the need for extensive molecular simulations or structural biology experiments. By using the Deep DNAshape method, DNA structural features can be predicted for any length and number of DNA sequences in a high-throughput manner, providing an understanding of the effects of flanking regions on DNA structure in a target region of a sequence. The Deep DNAshape method provides access to the influence of distant flanking regions on a region of interest. Our findings reveal that DNA shape readout mechanisms of a core target are quantitatively affected by flanking regions, including extended flanking regions, providing valuable insights into the detailed structural readout mechanisms of protein-DNA binding. Furthermore, when incorporated in machine learning models, the features generated by Deep DNAshape improve the model prediction accuracy. Collectively, Deep DNAshape can serve as versatile and powerful tool for diverse DNA structure-related studies.

Probing the role of the protonation state of a minor groove-linker histidine in Exd-Hox-DNA binding

DNA recognition and targeting by transcription factors (TFs) through specific binding are fundamental in biological processes. Furthermore, the histidine protonation state at the TF-DNA binding interface can significantly influence the binding mechanism of TF-DNA complexes. Nevertheless, the role of histidine in TF-DNA complexes remains underexplored. Here, we employed all-atom molecular dynamics simulations using AlphaFold2-modeled complexes based on previously solved co-crystal structures to probe the role of the His-12 residue in the Extradenticle (Exd)-Sex combs reduced (Scr)-DNA complex when binding to Scr and Ultrabithorax (Ubx) target sites. Our results demonstrate that the protonation state of histidine notably affected the DNA minor-groove width profile and binding free energy. Examining flanking sequences of various binding affinities derived from SELEX-seq experiments, we analyzed the relationship between binding affinity and specificity. We uncovered how histidine protonation leads to increased binding affinity but can lower specificity. Our findings provide new mechanistic insights into the role of histidine in modulating TF-DNA binding.

Dissection of integrated readout reveals the structural thermodynamics of DNA selection by transcription factors

Nucleobases such as inosine have been extensively utilized to map direct contacts by proteins in the DNA groove. Their deployment as targeted probes of dynamics and hydration, which are dominant thermodynamic drivers of affinity and specificity, has been limited by a paucity of suitable experimental models. We report a joint crystallographic, thermodynamic, and computational study of the bidentate complex of the arginine side chain with a Watson-Crick guanine (Arg×GC), a highly specific configuration adopted by major transcription factors throughout the eukaryotic branches in the Tree of Life. Using the ETS-family factor PU.1 as a high-resolution structural framework, inosine substitution for guanine resulted in a sharp dissection of conformational dynamics and hydration and elucidated their role in the DNA specificity of PU.1. Our work suggests an under-exploited utility of modified nucleobases in untangling the structural thermodynamics of interactions, such as the Arg×GC motif, where direct and indirect readout are tightly integrated.

Systematic Dissection of Sequence Features Affecting the Binding Specificity of a Pioneer Factor Reveals Binding Synergy Between FOXA1 and AP-1

Despite the unique ability of pioneer transcription factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called ChIP-ISO to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1. Combining ChIP-ISO with in vitro and neural network analyses, we find that 1) FOXA1 binding is strongly affected by co-binding TFs AP-1 and CEBPB, 2) FOXA1 and AP-1 show binding cooperativity in vitro, 3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin, and 4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

Deep DNAshape: Predicting DNA shape considering extended flanking regions using a deep learning method

Understanding the mechanisms of protein-DNA binding is critical in comprehending gene regulation. Three-dimensional DNA shape plays a key role in these mechanisms. In this study, we present a deep learning-based method, Deep DNAshape, that fundamentally changes the current k -mer based high-throughput prediction of DNA shape features by accurately accounting for the influence of extended flanking regions, without the need for extensive molecular simulations or structural biology experiments. By using the Deep DNAshape method, refined DNA shape features can be predicted for any length and number of DNA sequences in a high-throughput manner, providing a deeper understanding of the effects of flanking regions on DNA shape in a target region of a sequence. Deep DNAshape method provides access to the influence of distant flanking regions on a region of interest. Our findings reveal that DNA shape readout mechanisms of a core target are quantitatively affected by flanking regions, including extended flanking regions, providing valuable insights into the detailed structural readout mechanisms of protein-DNA binding. Furthermore, when incorporated in machine learning models, the features generated by Deep DNAshape improve the model prediction accuracy. Collectively, Deep DNAshape can serve as a versatile and powerful tool for diverse DNA structure-related studies.

Structural insights into the DNA recognition mechanism by the bacterial transcription factor PdxR

Specificity in protein-DNA recognition arises from the synergy of several factors that stem from the structural and chemical signatures encoded within the targeted DNA molecule. Here, we deciphered the nature of the interactions driving DNA recognition and binding by the bacterial transcription factor PdxR, a member of the MocR family responsible for the regulation of pyridoxal 5'-phosphate (PLP) biosynthesis. Single particle cryo-EM performed on the PLP-PdxR bound to its target DNA enabled the isolation of three conformers of the complex, which may be considered as snapshots of the binding process. Moreover, the resolution of an apo-PdxR crystallographic structure provided a detailed description of the transition of the effector domain to the holo-PdxR form triggered by the binding of the PLP effector molecule. Binding analyses of mutated DNA sequences using both wild type and PdxR variants revealed a central role of electrostatic interactions and of the intrinsic asymmetric bending of the DNA in allosterically guiding the holo-PdxR-DNA recognition process, from the first encounter through the fully bound state. Our results detail the structure and dynamics of the PdxR-DNA complex, clarifying the mechanism governing the DNA-binding mode of the holo-PdxR and the regulation features of the MocR family of transcription factors.

Structural underpinnings of mutation rate variations in the human genome

Single nucleotide mutation rates have critical implications for human evolution and genetic diseases. Importantly, the rates vary substantially across the genome and the principles underlying such variations remain poorly understood. A recent model explained much of this variation by considering higher-order nucleotide interactions in the 7-mer sequence context around mutated nucleotides. This model's success implicates a connection between DNA shape and mutation rates. DNA shape, i.e. structural properties like helical twist and tilt, is known to capture interactions between nucleotides within a local context. Thus, we hypothesized that changes in DNA shape features at and around mutated positions can explain mutation rate variations in the human genome. Indeed, DNA shape-based models of mutation rates showed similar or improved performance over current nucleotide sequence-based models. These models accurately characterized mutation hotspots in the human genome and revealed the shape features whose interactions underlie mutation rate variations. DNA shape also impacts mutation rates within putative functional regions like transcription factor binding sites where we find a strong association between DNA shape and position-specific mutation rates. This work demonstrates the structural underpinnings of nucleotide mutations in the human genome and lays the groundwork for future models of genetic variations to incorporate DNA shape.

Entropic analysis of antigen-specific CDR3 domains identifies essential binding motifs shared by CDR3s with different antigen specificities

Antigen-specific T cell receptor (TCR) sequences can have prognostic, predictive, and therapeutic value, but decoding the specificity of TCR recognition remains challenging. Unlike DNA strands that base pair, TCRs bind to their targets with different orientations and different lengths, which complicates comparisons. We present scanning parametrized by normalized TCR length (SPAN-TCR) to analyze antigen-specific TCR CDR3 sequences and identify patterns driving TCR-pMHC specificity. Using entropic analysis, SPAN-TCR identifies 2-mer motifs that decrease the diversity (entropy) of CDR3s. These motifs are the most common patterns that can predict CDR3 composition, and we identify "essential" motifs that decrease entropy in the same CDR3 α or β chain containing the 2-mer, and "super-essential" motifs that decrease entropy in both chains. Molecular dynamics analysis further suggests that these motifs may play important roles in binding. We then employ SPAN-TCR to resolve similarities in TCR repertoires against different antigens using public databases of TCR sequences.

Predicting accurate ab initio DNA electron densities with equivariant neural networks

One of the fundamental limitations of accurately modeling biomolecules like DNA is the inability to perform quantum chemistry calculations on large molecular structures. We present a machine learning model based on an equivariant Euclidean neural network framework to obtain accurate ab initio electron densities for arbitrary DNA structures that are much too large for conventional quantum methods. The model is trained on representative B-DNA basepair steps that capture both base pairing and base stacking interactions. The model produces accurate electron densities for arbitrary B-DNA structures with typical errors of less than 1%. Crucially, the error does not increase with system size, which suggests that the model can extrapolate to large DNA structures with negligible loss of accuracy. The model also generalizes reasonably to other DNA structural motifs such as the A- and Z-DNA forms, despite being trained on only B-DNA configurations. The model is used to calculate electron densities of several large-scale DNA structures, and we show that the computational scaling for this model is essentially linear. We also show that this machine learning electron density model can be used to calculate accurate electrostatic potentials for DNA. These electrostatic potentials produce more accurate results compared with classical force fields and do not show the usual deficiencies at short range.

Host-pathogen protein-nucleic acid interactions: A comprehensive review

Recognition of pathogen-derived nucleic acids by host cells is an effective host strategy to detect pathogenic invasion and trigger immune responses. In the context of pathogen-specific pharmacology, there is a growing interest in mapping the interactions between pathogen-derived nucleic acids and host proteins. Insight into the principles of the structural and immunological mechanisms underlying such interactions and their roles in host defense is necessary to guide therapeutic intervention. Here, we discuss the newest advances in studies of molecular interactions involving pathogen nucleic acids and host factors, including their drug design, molecular structure and specific patterns. We observed that two groups of nucleic acid recognizing molecules, Toll-like receptors (TLRs) and the cytoplasmic retinoic acid-inducible gene (RIG)-I-like receptors (RLRs) form the backbone of host responses to pathogen nucleic acids, with additional support provided by absent in melanoma 2 (AIM2) and DNA-dependent activator of Interferons (IFNs)-regulatory factors (DAI) like cytosolic activity. We review the structural, immunological, and other biological aspects of these representative groups of molecules, especially in terms of their target specificity and affinity and challenges in leveraging host-pathogen protein-nucleic acid interactions (HP-PNI) in drug discovery.

Bacterial H-NS contacts DNA at the same irregularly spaced sites in both bridged and hemi-sequestered linear filaments

Gene silencing in bacteria is mediated by chromatin proteins, of which Escherichia coli H-NS is a paradigmatic example. H-NS forms nucleoprotein filaments with either one or two DNA duplexes. However, the structures, arrangements of DNA-binding domains (DBDs), and positions of DBD-DNA contacts in linear and bridged filaments are uncertain. To characterize the H-NS DBD contacts that silence transcription by RNA polymerase, we combined ·OH footprinting, molecular dynamics, statistical modeling, and DBD mapping using a chemical nuclease (Fe²⁺-EDTA) tethered to the DBDs (TEN-map). We find that H-NS DBDs contact DNA at indistinguishable locations in bridged or linear filaments and that the DBDs vary in orientation and position with ∼10-bp average spacing. Our results support a hemi-sequestration model of linear-to-bridged H-NS switching. Linear filaments able to inhibit only transcription initiation switch to bridged filaments able to inhibit both initiation and elongation using the same irregularly spaced DNA contacts.

Insights into protein-DNA interactions from hydrogen bond energy-based comparative protein-ligand analyses

Hydrogen bonds play important roles in protein folding and protein-ligand interactions, particularly in specific protein-DNA recognition. However, the distributions of hydrogen bonds, especially hydrogen bond energy (HBE) in different types of protein-ligand complexes, is unknown. Here we performed a comparative analysis of hydrogen bonds among three non-redundant datasets of protein-protein, protein-peptide, and protein-DNA complexes. Besides comparing the number of hydrogen bonds in terms of types and locations, we investigated the distributions of HBE. Our results indicate that while there is no significant difference of hydrogen bonds within protein chains among the three types of complexes, interfacial hydrogen bonds are significantly more prevalent in protein-DNA complexes. More importantly, the interfacial hydrogen bonds in protein-DNA complexes displayed a unique energy distribution of strong and weak hydrogen bonds whereas majority of the interfacial hydrogen bonds in protein-protein and protein-peptide complexes are of predominantly high strength with low energy. Moreover, there is a significant difference in the energy distributions of minor groove hydrogen bonds between protein-DNA complexes with different binding specificity. Highly specific protein-DNA complexes contain more strong hydrogen bonds in the minor groove than multi-specific complexes, suggesting important role of minor groove in specific protein-DNA recognition. These results can help better understand protein-DNA interactions and have important implications in improving quality assessments of protein-DNA complex models.

Human Cytomegalovirus IE2 Both Activates and Represses Initiation and Modulates Elongation in a Context-Dependent Manner

Human cytomegalovirus (HCMV) immediate-early 2 (IE2) protein is a multifunctional transcription factor that is essential for lytic HCMV infection. IE2 functions as an activator of viral early genes, negatively regulates its own promoter, and is required for viral replication. The mechanisms by which IE2 executes these distinct functions are incompletely understood. Using PRO-Seq, which profiles nascent transcripts, and a recently developed DFF-chromatin immunoprecipitation (DFF-ChIP; employs chromatin digestion by the endonuclease DNA fragmentation factor prior to IP) approach that resolves occupancy and local chromatin environment, we show that IE2 controls viral gene transcription in three distinct capacities during late HCMV infection and reveal mechanisms that involve direct binding of IE2 to viral DNA. IE2 represses a subset of viral promoters by binding within their core promoter regions and blocking the assembly of preinitiation complexes (PICs). Remarkably, IE2 forms a repressive complex at the major immediate-early promoter region involving direct association of IE2 with nucleosomes and TBP. IE2 stimulates transcription by binding nearby, but not within, core promoter regions. In addition, IE2 functions as a direct roadblock to transcription elongation. At one locus, this function of IE2 appears to be important for the synthesis of a spliced viral RNA. Consistent with the minimal observed effects of IE2 depletion on host gene transcription, IE2 does not functionally engage the host genome. Our results reveal mechanisms of transcriptional control by IE2, uncover a previously unknown function of IE2 as a Pol II elongation modulator, and demonstrate that DFF-ChIP is a useful tool for probing transcription factor occupancy and interactions between transcription factors and nucleosomes at high resolution.

Structural insights into cooperative DNA recognition by the CCAAT-binding complex and its bZIP transcription factor HapX

The heterotrimeric CCAAT-binding complex (CBC) is a fundamental eukaryotic transcription factor recognizing the CCAAT box. In certain fungi, like Aspergilli, the CBC cooperates with the basic leucine zipper HapX to control iron metabolism. HapX functionally depends on the CBC, and the stable interaction of both requires DNA. To study this cooperative effect, X-ray structures of the CBC-HapX-DNA complex were determined. Downstream of the CCAAT box, occupied by the CBC, a HapX dimer binds to the major groove. The leash-like N terminus of the distal HapX subunit contacts the CBC, and via a flexible polyproline type II helix mediates minor groove interactions that stimulate sequence promiscuity. In vitro and in vivo mutagenesis suggest that the structural and functional plasticity of HapX results from local asymmetry and its ability to target major and minor grooves simultaneously. The latter feature may also apply to related transcription factors such as yeast Hap4 and distinct Yap family members.

Analysis of nucleoid-associated protein-binding regions reveals DNA structural features influencing genome organization in Mycobacterium tuberculosis

Nucleoid-associated proteins (NAPs) maintain bacterial nucleoid configuration through their architectural properties of DNA bending, wrapping, and bridging. However, the contribution of DNA structural alterations to DNA-NAP recognition at the genomic scale remains unresolved. Present work dissects the DNA sequence, shape and altered structural preferences at a genomic scale for six NAPs in Mycobacterium tuberculosis. Results suggest narrower minor groove width (MGW) and higher DNA rigidity are marked for the binding sites of EspR and Lsr2, while mIHF, MtHU and NapM have heterogeneous DNA structural predilections. In contrast, WhiB4-DNA-binding sites were characterized by wider MGW, highly deformable and less curved DNA. This work provides systematic insight into NAP-mediated genome organization as a function of DNA structural features.

The bacterial iron sensor IdeR recognizes its DNA targets by indirect readout

The iron-dependent regulator IdeR is the main transcriptional regulator controlling iron homeostasis genes in Actinobacteria, including species from the Corynebacterium, Mycobacterium and Streptomyces genera, as well as the erythromycin-producing bacterium Saccharopolyspora erythraea. Despite being a well-studied transcription factor since the identification of the Diphtheria toxin repressor DtxR three decades ago, the details of how IdeR proteins recognize their highly conserved 19-bp DNA target remain to be elucidated. IdeR makes few direct contacts with DNA bases in its target sequence, and we show here that these contacts are not required for target recognition. The results of our structural and mutational studies support a model wherein IdeR mainly uses an indirect readout mechanism, identifying its targets via the sequence-dependent DNA backbone structure rather than through specific contacts with the DNA bases. Furthermore, we show that IdeR efficiently recognizes a shorter palindromic sequence corresponding to a half binding site as compared to the full 19-bp target previously reported, expanding the number of potential target genes controlled by IdeR proteins.

The importance of hydrophobic interactions in the structure of transcription systems

Hydrophobic forces play a crucial role in both the stability of B DNA and its interactions with proteins. In the present study, we postulate that the hydrophobic effect is an essential component in establishing specificity in the interaction transcription factor proteins with their consensus DNA sequence partners. The PDB coordinates of more than 50 transcription systems have been used to analyze the hydrophobic attraction of proteins towards their DNA consensus. This analysis includes computing the hydrophobic energy of the interacting molecules by means of their hydrophobic moments. Hydrophobic moments have successfully been used in previous studies involving self-assembly protein systems. In the present case, in spite of some variability, we found specificity in transcription factors when interacting with their respective consensus DNA sequences. By applying our model of biological membrane pattern for hydrophobic interactions, we postulate that hydrophobic forces constitute the necessary intermediate interaction between the unspecific electrostatic attraction for DNA phosphate groups and the very short-range interaction promoting hydrogen bonds. We conclude that hydrophobic interactions serve as the intermediate force guiding transcriptions factors towards the proper hydrogen bonds to their DNAs.

Histone H4 Tails in Nucleosomes: a Fuzzy Interaction with DNA

The interaction of positively charged N-terminal histone tails with nucleosomal DNA plays an important role in chromatin assembly and regulation, modulating their susceptibility to post-translational modifications and recognition by chromatin-binding proteins. Here, we report residue-specific ¹⁵ N NMR relaxation rates for histone H4 tails in reconstituted nucleosomes. These data indicate that H4 tails are strongly dynamically disordered, albeit with reduced conformational flexibility compared to a free peptide with the same sequence. Remarkably, the NMR observables were successfully reproduced in a 2-μs MD trajectory of the nucleosome. This is an important step toward resolving an apparent inconsistency where prior simulations were generally at odds with experimental evidence on conformational dynamics of histone tails. Our findings indicate that histone H4 tails engage in a fuzzy interaction with nucleosomal DNA, underpinned by a variable pattern of short-lived salt bridges and hydrogen bonds, which persists at low ionic strength (0-100 mM NaCl).

Distinct mechanisms control genome recognition by p53 at its target genes linked to different cell fates

The tumor suppressor p53 integrates stress response pathways by selectively engaging one of several potential transcriptomes, thereby triggering cell fate decisions (e.g., cell cycle arrest, apoptosis). Foundational to this process is the binding of tetrameric p53 to 20-bp response elements (REs) in the genome (RRRCWWGYYYN_0-13RRRCWWGYYY). In general, REs at cell cycle arrest targets (e.g. p21) are of higher affinity than those at apoptosis targets (e.g., BAX). However, the RE sequence code underlying selectivity remains undeciphered. Here, we identify molecular mechanisms mediating p53 binding to high- and low-affinity REs by showing that key determinants of the code are embedded in the DNA shape. We further demonstrate that differences in minor/major groove widths, encoded by G/C or A/T bp content at positions 3, 8, 13, and 18 in the RE, determine distinct p53 DNA-binding modes by inducing different Arg248 and Lys120 conformations and interactions. The predictive capacity of this code was confirmed in vivo using genome editing at the BAX RE to interconvert the DNA-binding modes, transcription pattern, and cell fate outcome.

Humanizing the yeast origin recognition complex

The Origin Recognition Complex (ORC) is an evolutionarily conserved six-subunit protein complex that binds specific sites at many locations to coordinately replicate the entire eukaryote genome. Though highly conserved in structure, ORC’s selectivity for replication origins has diverged tremendously between yeasts and humans to adapt to vastly different life cycles. In this work, we demonstrate that the selectivity determinant of ORC for DNA binding lies in a 19-amino acid insertion helix in the Orc4 subunit, which is present in yeast but absent in human. Removal of this motif from Orc4 transforms the yeast ORC, which selects origins based on base-specific binding at defined locations, into one whose selectivity is dictated by chromatin landscape and afforded with plasticity, as reported for human. Notably, the altered yeast ORC has acquired an affinity for regions near transcriptional start sites (TSSs), which the human ORC also favors.

In most model yeast species the Origin Recognition Complex (ORC) binds defined and species-specific base sequences while in humans what determines the binding appears to be more complex. Here the authors reveal that the yeast’s ORC complex binding specificity is dependent on a 19-amino acid insertion helix in the Orc4 subunit which is lost in human.

Transcription Factor Binding Affinities and DNA Shape Readout

An essential event in gene regulation is the binding of a transcription factor (TF) to its target DNA. Models considering the interactions between the TF and the DNA geometry proved to be successful approaches to describe this binding event, while conserving data interpretability. However, a direct characterization of the DNA shape contribution to binding is still missing due to the lack of accurate and large-scale binding affinity data. Here, we use a binding assay we recently established to measure with high sensitivity the binding specificities of 13 Drosophila TFs, including dinucleotide dependencies to capture non-independent amino acid-base interactions. Correlating the binding affinities with all DNA shape features, we find that shape readout is widely used by these factors. A shape readout/TF-DNA complex structure analysis validates our approach while providing biological insights such as positively charged or highly polar amino acids often contact nucleotides that exhibit strong shape readout.

•

The DNA shape contribution to Drosophila TFs-DNA binding is directly characterized

•

Zeroth- and first-order TF-DNA binding specificities are measured with high accuracy

•

DNA shape readout is widely used by these TFs

•

A shape readout/structural correlation analysis provides biological insights

Landscape of DNA binding signatures of myocyte enhancer factor-2B reveals a unique interplay of base and shape readout

Myocyte enhancer factor-2B (MEF2B) has the unique capability of binding to its DNA target sites with a degenerate motif, while still functioning as a gene-specific transcriptional regulator. Identifying its DNA targets is crucial given regulatory roles exerted by members of the MEF2 family and MEF2B's involvement in B-cell lymphoma. Analyzing structural data and SELEX-seq experimental results, we deduced the DNA sequence and shape determinants of MEF2B target sites on a high-throughput basis in vitro for wild-type and mutant proteins. Quantitative modeling of MEF2B binding affinities and computational simulations exposed the DNA readout mechanisms of MEF2B. The resulting binding signature of MEF2B revealed distinct intricacies of DNA recognition compared to other transcription factors. MEF2B uses base readout at its half-sites combined with shape readout at the center of its degenerate motif, where A-tract polarity dictates nuances of binding. The predominant role of shape readout at the center of the core motif, with most contacts formed in the minor groove, differs from previously observed protein-DNA readout modes. MEF2B, therefore, represents a unique protein for studies of the role of DNA shape in achieving binding specificity. MEF2B-DNA recognition mechanisms are likely representative for other members of the MEF2 family.

Epigenetic competition reveals density-dependent regulation and target site plasticity of phosphorothioate epigenetics in bacteria

Phosphorothioate (PT) DNA modifications-in which a nonbonding phosphate oxygen is replaced with sulfur-represent a widespread, horizontally transferred epigenetic system in prokaryotes and have a highly unusual property of occupying only a small fraction of available consensus sequences in a genome. Using Salmonella enterica as a model, we asked a question of fundamental importance: How do the PT-modifying DndA-E proteins select their G_PSAAC/G_PSTTC targets? Here, we applied innovative analytical, sequencing, and computational tools to discover a novel behavior for DNA-binding proteins: The Dnd proteins are "parked" at the G^6mATC Dam methyltransferase consensus sequence instead of the expected GAAC/GTTC motif, with removal of the ^6mA permitting extensive PT modification of GATC sites. This shift in modification sites further revealed a surprising constancy in the density of PT modifications across the genome. Computational analysis showed that GAAC, GTTC, and GATC share common features of DNA shape, which suggests that PT epigenetics are regulated in a density-dependent manner partly by DNA shape-driven target selection in the genome.

TFBSshape: an expanded motif database for DNA shape features of transcription factor binding sites

TFBSshape (https://tfbsshape.usc.edu) is a motif database for analyzing structural profiles of transcription factor binding sites (TFBSs). The main rationale for this database is to be able to derive mechanistic insights in protein-DNA readout modes from sequencing data without available structures. We extended the quantity and dimensionality of TFBSshape, from mostly in vitro to in vivo binding and from unmethylated to methylated DNA. This new release of TFBSshape improves its functionality and launches a responsive and user-friendly web interface for easy access to the data. The current expansion includes new entries from the most recent collections of transcription factors (TFs) from the JASPAR and UniPROBE databases, methylated TFBSs derived from in vitro high-throughput EpiSELEX-seq binding assays and in vivo methylated TFBSs from the MeDReaders database. TFBSshape content has increased to 2428 structural profiles for 1900 TFs from 39 different species. The structural profiles for each TFBS entry now include 13 shape features and minor groove electrostatic potential for standard DNA and four shape features for methylated DNA. We improved the flexibility and accuracy for the shape-based alignment of TFBSs and designed new tools to compare methylated and unmethylated structural profiles of TFs and methods to derive DNA shape-preserving nucleotide mutations in TFBSs.

Systematic search for structural motifs of peptide binding to double-stranded DNA

A large variety of short biologically active peptides possesses antioxidant, antibacterial, antitumour, anti-ageing and anti-inflammatory activity, involved in the regulation of neuro-immuno-endocrine system functions, cell apoptosis, proliferation and differentiation. Therefore, the mechanisms of their biological activity are attracting increasing attention not only in modern molecular biology, biochemistry and biophysics, but also in pharmacology and medicine. In this work, we systematically analysed the ability of dipeptides (all possible combinations of the 20 standard amino acids) to bind all possible combinations of tetra-nucleotides in the central part of dsDNA in the classic B-form using molecular docking and molecular dynamics. The vast majority of the dipeptides were found to be unable to bind dsDNA. However, we were able to identify 57 low-energy dipeptide complexes with peptide-dsDNA possessing high selectivity for DNA binding. The analysis of the dsDNA complexes with dipeptides with free and blocked N- and C-terminus showed that selective peptide binding to dsDNA can increase dramatically with the peptide length.

Crystal structures of REF6 and its complex with DNA reveal diverse recognition mechanisms

Relative of Early Flowing 6 (REF6) is a DNA-sequence-specific H3K27me3/2 demethylase that contains four zinc finger (ZnF) domains and targets several thousand genes in Arabidopsis thaliana. The ZnF domains are essential for binding target genes, but the structural basis remains unclear. Here, we determined crystal structures of the ZnF domains and REF6-DNA complex, revealing a unique REF6-family-specific half-cross-braced ZnF (RCZ) domain and two C2H2-type ZnFs. DNA-binding induces a profound conformational change in the hinge region of REF6. Each REF6 recognizes six bases and DNA methylation reduces the binding affinity. Both the acidic region and basic region are important for the self-association of REF6. The REF6 DNA-binding affinity is determined by the sequence-dependent conformations of DNA and also the cooperativity in different target motifs. The conformational plasticity enables REF6 to function as a global transcriptional regulator that directly binds to many diverse genes, revealing the structural basis for the epigenetic modification recognition.

Co-SELECT reveals sequence non-specific contribution of DNA shape to transcription factor binding in vitro

Understanding the principles of DNA binding by transcription factors (TFs) is of primary importance for studying gene regulation. Recently, several lines of evidence suggested that both DNA sequence and shape contribute to TF binding. However, the following compelling question is yet to be considered: in the absence of any sequence similarity to the binding motif, can DNA shape still increase binding probability? To address this challenge, we developed Co-SELECT, a computational approach to analyze the results of in vitro HT-SELEX experiments for TF-DNA binding. Specifically, Co-SELECT leverages the presence of motif-free sequences in late HT-SELEX rounds and their enrichment in weak binders allows Co-SELECT to detect an evidence for the role of DNA shape features in TF binding. Our approach revealed that, even in the absence of the sequence motif, TFs have propensity to bind to DNA molecules of the shape consistent with the motif specific binding. This provides the first direct evidence that shape features that accompany the preferred sequence motifs also bestow an advantage for weak, sequence non-specific binding.

Expanding the uses of genome-scale models with protein structures

Biology is reaching a convergence point of its historic reductionist and modern holistic approaches to understanding the living system. Structural biology has historically taken the reductionist approach to deeply probe the inner workings of complex molecular machines. In contrast, systems biology and genome-scale modeling have organically grown out of the wealth of data now being generated by diverse omics measurements. In the late 2000s, a proposed interdisciplinary field of structural systems biology pitched the merger of these two approaches, with widespread applications in pharmacology, disease modeling, protein engineering, and evolutionary studies. In this commentary, we highlight the challenges of integrating these two fields, with a focus on genome-scale metabolic modeling, and the novel findings that are made possible from such a merger.

Cooperative DNA binding by proteins through DNA shape complementarity

Localized arrays of proteins cooperatively assemble onto chromosomes to control DNA activity in many contexts. Binding cooperativity is often mediated by specific protein-protein interactions, but cooperativity through DNA structure is becoming increasingly recognized as an additional mechanism. During the site-specific DNA recombination reaction that excises phage λ from the chromosome, the bacterial DNA architectural protein Fis recruits multiple λ-encoded Xis proteins to the attR recombination site. Here, we report X-ray crystal structures of DNA complexes containing Fis + Xis, which show little, if any, contacts between the two proteins. Comparisons with structures of DNA complexes containing only Fis or Xis, together with mutant protein and DNA binding studies, support a mechanism for cooperative protein binding solely by DNA allostery. Fis binding both molds the minor groove to potentiate insertion of the Xis β-hairpin wing motif and bends the DNA to facilitate Xis-DNA contacts within the major groove. The Fis-structured minor groove shape that is optimized for Xis binding requires a precisely positioned pyrimidine-purine base-pair step, whose location has been shown to modulate minor groove widths in Fis-bound complexes to different DNA targets.

Dissection of acute stimulus-inducible nucleosome remodeling in mammalian cells

Accessibility of the genomic regulatory information is largely controlled by the nucleosome-organizing activity of transcription factors (TFs). While stimulus-induced TFs bind to genomic regions that are maintained accessible by lineage-determining TFs, they also increase accessibility of thousands of cis-regulatory elements. Nucleosome remodeling events underlying such changes and their interplay with basal positioning are unknown. Here, we devised a novel quantitative framework discriminating different types of nucleosome remodeling events in micrococcal nuclease ChIP-seq (chromatin immunoprecipitation [ChIP] combined with high-throughput sequencing) data sets and used it to analyze nucleosome dynamics at stimulus-regulated cis-regulatory elements. At enhancers, remodeling preferentially affected poorly positioned nucleosomes while sparing well-positioned nucleosomes flanking the enhancer core, indicating that inducible TFs do not suffice to overrule basal nucleosomal organization maintained by lineage-determining TFs. Remodeling events appeared to be combinatorially driven by multiple TFs, with distinct TFs showing, however, different remodeling efficiencies. Overall, these data provide a systematic view of the impact of stimulation on nucleosome organization and genome accessibility in mammalian cells.

Experimental maps of DNA structure at nucleotide resolution distinguish intrinsic from protein-induced DNA deformations

Recognition of DNA by proteins depends on DNA sequence and structure. Often unanswered is whether the structure of naked DNA persists in a protein–DNA complex, or whether protein binding changes DNA shape. While X-ray structures of protein–DNA complexes are numerous, the structure of naked cognate DNA is seldom available experimentally. We present here an experimental and computational analysis pipeline that uses hydroxyl radical cleavage to map, at single-nucleotide resolution, DNA minor groove width, a recognition feature widely exploited by proteins. For 11 protein–DNA complexes, we compared experimental maps of naked DNA minor groove width with minor groove width measured from X-ray co-crystal structures. Seven sites had similar minor groove widths as naked DNA and when bound to protein. For four sites, part of the DNA in the complex had the same structure as naked DNA, and part changed structure upon protein binding. We compared the experimental map with minor groove patterns of DNA predicted by two computational approaches, DNAshape and ORChID2, and found good but not perfect concordance with both. This experimental approach will be useful in mapping structures of DNA sequences for which high-resolution structural data are unavailable. This approach allows probing of protein family-dependent readout mechanisms.

Analysis of Genetic Variation Indicates DNA Shape Involvement in Purifying Selection

Noncoding DNA sequences, which play various roles in gene expression and regulation, are under evolutionary pressure. Gene regulation requires specific protein–DNA binding events, and our previous studies showed that both DNA sequence and shape readout are employed by transcription factors (TFs) to achieve DNA binding specificity. By investigating the shape-disrupting properties of single nucleotide polymorphisms (SNPs) in human regulatory regions, we established a link between disruptive local DNA shape changes and loss of specific TF binding. Furthermore, we described cases where disease-associated SNPs may alter TF binding through DNA shape changes. This link led us to hypothesize that local DNA shape within and around TF binding sites is under selection pressure. To verify this hypothesis, we analyzed SNP data derived from 216 natural strains of Drosophila melanogaster. Comparing SNPs located in functional and nonfunctional regions within experimentally validated cis-regulatory modules (CRMs) from D. melanogaster that are active in the blastoderm stage of development, we found that SNPs within functional regions tended to cause smaller DNA shape variations. Furthermore, SNPs with higher minor allele frequency were more likely to result in smaller DNA shape variations. The same analysis based on a large number of SNPs in putative CRMs of the D. melanogaster genome derived from DNase I accessibility data confirmed these observations. Taken together, our results indicate that common SNPs in functional regions tend to maintain DNA shape, whereas shape-disrupting SNPs are more likely to be eliminated through purifying selection.

DNA-protein interaction: identification, prediction and data analysis

Life in living organisms is dependent on specific and purposeful interaction between other molecules. Such purposeful interactions make the various processes inside the cells and the bodies of living organisms possible. DNA-protein interactions, among all the types of interactions between different molecules, are of considerable importance. Currently, with the development of numerous experimental techniques, diverse methods are convenient for recognition and investigating such interactions. While the traditional experimental techniques to identify DNA-protein complexes are time-consuming and are unsuitable for genome-scale studies, the current high throughput approaches are more efficient in determining such interaction at a large-scale, but they are clearly too costly to be practice for daily applications. Hence, according to the availability of much information related to different biological sequences and clearing different dimensions of conditions in which such interactions are formed, with the developments related to the computer, mathematics, and statistics motivate scientists to develop bioinformatics tools for prediction the interaction site(s). Until now, there has been much progress in this field. In this review, the factors and conditions governing the interaction and the laboratory techniques for examining such interactions are addressed. In addition, developed bioinformatics tools are introduced and compared for this reason and, in the end, several suggestions are offered for the promotion of such tools in prediction with much more precision.

New Insights into the Interactions of a DNA Oligonucleotide with mPEGylated-PAMAM by Circular Dichroism and Solution NMR

Dendrimers are well-defined, highly branched, synthetic three-dimensional molecules with a large number of reactive end groups. PAMAM dendrimers form stable complexes with DNA chemistries and constitute an important class of nonviral, cationic vectors in gene delivery. The aim of this study is to examine the interactions of a 12 bp DNA oligonucletide with PAMAM-G2 and mPEG- b-PAMAM-G3 having eight surface amine groups under physiological conditions, using constant DNA concentration but varying dendrimer concentration. 1D ³¹P NMR, 2D NOESY, and CD spectroscopic methods were employed to investigate the interactions between the dendrimer and the DNA. The CD experiments carried out with a constant DNA concentration of 25 μM and dendrimer concentrations from 0 to 100 μM indicated minimal change to the chirality of the DNA for both types of dendrimers. While the PAMAM-G2 dendrimer caused aggregation of the majority of the DNA, the 2D NMR data of the DNA with an mPEG- b-PAMAM-G3 dendrimer indicated general broadening of the 1D ³¹P peaks from the DNA phosphates, a small number of ¹H chemical shift perturbations (CSPs), and reduction of specific ¹H-¹H NOE intensities. These data suggest there is minimal structural alteration of the DNA in the complex and indicate preferential binding of the dendrimer to the central AATT region of the DNA sequence. The results herein are the first such results demonstrating a soluble DNA complex with the mPEG- b-PAMAM-G3 dendrimer analyzed by multidimensional NMR.

The interaction landscape between transcription factors and the nucleosome

Nucleosomes cover most of the genome and are thought to be displaced by transcription factors in regions that direct gene expression. However, the modes of interaction between transcription factors and nucleosomal DNA remain largely unknown. Here we systematically explore interactions between the nucleosome and 220 transcription factors representing diverse structural families. Consistent with earlier observations, we find that the majority of the studied transcription factors have less access to nucleosomal DNA than to free DNA. The motifs recovered from transcription factors bound to nucleosomal and free DNA are generally similar. However, steric hindrance and scaffolding by the nucleosome result in specific positioning and orientation of the motifs. Many transcription factors preferentially bind close to the end of nucleosomal DNA, or to periodic positions on the solvent-exposed side of the DNA. In addition, several transcription factors usually bind to nucleosomal DNA in a particular orientation. Some transcription factors specifically interact with DNA located at the dyad position at which only one DNA gyre is wound, whereas other transcription factors prefer sites spanning two DNA gyres and bind specifically to each of them. Our work reveals notable differences in the binding of transcription factors to free and nucleosomal DNA, and uncovers a diverse interaction landscape between transcription factors and the nucleosome.

Relationship between histone modifications and transcription factor binding is protein family specific

The very small fraction of putative binding sites (BSs) that are occupied by transcription factors (TFs) in vivo can be highly variable across different cell types. This observation has been partly attributed to changes in chromatin accessibility and histone modification (HM) patterns surrounding BSs. Previous studies focusing on BSs within DNA regulatory regions found correlations between HM patterns and TF binding specificities. However, a mechanistic understanding of TF-DNA binding specificity determinants is still not available. The ability to predict in vivo TF binding on a genome-wide scale requires the identification of features that determine TF binding based on evolutionary relationships of DNA binding proteins. To reveal protein family-dependent mechanisms of TF binding, we conducted comprehensive comparisons of HM patterns surrounding BSs and non-BSs with exactly matched core motifs for TFs in three cell lines: 33 TFs in GM12878, 37 TFs in K562, and 18 TFs in H1-hESC. These TFs displayed protein family-specific preferences for HM patterns surrounding BSs, with high agreement among cell lines. Moreover, compared to models based on DNA sequence and shape at flanking regions of BSs, HM-augmented quantitative machine-learning methods resulted in increased performance in a TF family-specific manner. Analysis of the relative importance of features in these models indicated that TFs, displaying larger HM pattern differences between BSs and non-BSs, bound DNA in an HM-specific manner on a protein family-specific basis. We propose that TF family-specific HM preferences reveal distinct mechanisms that assist in guiding TFs to their cognate BSs by altering chromatin structure and accessibility.

A unified approach for quantifying and interpreting DNA shape readout by transcription factors

Transcription factors (TFs) interpret DNA sequence by probing the chemical and structural properties of the nucleotide polymer. DNA shape is thought to enable a parsimonious representation of dependencies between nucleotide positions. Here, we propose a unified mathematical representation of the DNA sequence dependence of shape and TF binding, respectively, which simplifies and enhances analysis of shape readout. First, we demonstrate that linear models based on mononucleotide features alone account for 60-70% of the variance in minor groove width, roll, helix twist, and propeller twist. This explains why simple scoring matrices that ignore all dependencies between nucleotide positions can partially account for DNA shape readout by a TF Adding dinucleotide features as sequence-to-shape predictors to our model, we can almost perfectly explain the shape parameters. Building on this observation, we developed a post hoc analysis method that can be used to analyze any mechanism-agnostic protein-DNA binding model in terms of shape readout. Our insights provide an alternative strategy for using DNA shape information to enhance our understanding of how cis-regulatory codes are interpreted by the cellular machinery.

Substructure-activity relationship studies on antibody recognition for phenylurea compounds using competitive immunoassay and computational chemistry

Based on the structural features of fluometuron, an immunizing hapten was synthesized and conjugated to bovine serum albumin as an immunogen to prepare a polyclonal antibody. However, the resultant antibody indicated cross-reactivity with 6 structurally similar phenylurea herbicides, with binding activities (expressed by IC₅₀ values) ranging from 1.67 µg/L to 42.71 µg/L. All 6 phenylurea herbicides contain a common moiety and three different substitutes. To understand how these three different chemical groups affect the antibody-phenylurea recognition activity, quantum chemistry, using density function theory (DFT) at the B3LYP/6-311++ G(d,p) level of theory, was employed to optimize all phenylurea structures, followed by determination of the 3D conformations of these molecules, pharmacophore analysis, and molecular electrostatic potential (ESP) analysis. The molecular modeling results confirmed that the geometry configuration, pharmacophore features and electron distribution in the substituents were related to the antibody binding activity. Spearman correlation analysis further elucidated that the geometrical and electrostatic properties on the van der Waals (vdW) surface of the substituents played a critical role in the antibody-phenylurea recognition process.