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

Further Development of SAMPDI-3D: A Machine Learning Method for Predicting Binding Free Energy Changes Caused by Mutations in Either Protein or DNA

BACKGROUND/OBJECTIVES

Predicting the effects of protein and DNA mutations on the binding free energy of protein-DNA complexes is crucial for understanding how DNA variants impact wild-type cellular function. As many cellular interactions involve protein-DNA binding, accurately predicting changes in binding free energy (ΔΔG) is valuable for distinguishing pathogenic mutations from benign ones.

METHODS

This study describes the development and optimization of the SAMPDI-3Dv2 machine learning method, which is trained on an expanded database of experimentally measured ΔΔGs. This enhanced model incorporates new features, including the 3D structure of the mutant protein, features of the mutant structure, and a position-specific scoring matrix (PSSM). Benchmarking was conducted using 5-fold cross-validation.

RESULTS

The updated SAMPDI-3D model (SAMPDI-3Dv2) achieved Pearson correlation coefficients (PCCs) of 0.68 for protein and 0.80 for DNA mutations. These results represent significant improvements over existing tools. Additionally, the method's rapid execution time enables genome-scale predictions.

CONCLUSIONS

The improved SAMPDI-3Dv2 shows enhanced predictive performance for analyzing mutations in protein-DNA complexes. By leveraging structural information and an expanded training dataset, SAMPDI-3Dv2 provides researchers with a more accurate and efficient tool for mutation analysis, contributing to identifying pathogenic variants and improving our understanding of cellular function.

A Rare Noncoding Enhancer Variant in SCN5A Contributes to the High Prevalence of Brugada Syndrome in Thailand

BACKGROUND

Brugada syndrome (BrS) is a cardiac arrhythmia disorder that causes sudden death in young adults. Rare genetic variants in the SCN5A gene encoding the Nav1.5 sodium channel and common noncoding variants at this locus are robustly associated with the condition. BrS is particularly prevalent in Southeast Asia but the underlying ancestry-specific factors remain largely unknown.

METHODS

Genome sequencing of BrS probands and population-matched controls from Thailand was performed to identify rare noncoding variants at the SCN5A-SCN10A locus that were enriched in patients with BrS. A likely causal variant was prioritized by computational methods and introduced into human induced pluripotent stem cell (hiPSC) lines using CRISPR-Cas9. The effect of the variant on SCN5A expression and Nav1.5 sodium channel current was then assessed in hiPSC-derived cardiomyocytes (hiPSC-CMs).

RESULTS

A rare noncoding variant in an SCN5A intronic enhancer region was highly enriched in patients with BrS (detected in 3.9% of cases with a case-control odds ratio of 45.2). The variant affects a nucleotide conserved across all mammalian species and predicted to disrupt a Mef2 transcription factor binding site. Heterozygous introduction of the enhancer variant in hiPSC-CMs caused significantly reduced SCN5A expression from the variant-containing allele and a 30% reduction in Nav1.5-mediated sodium current density compared with isogenic controls, confirming its pathogenicity. Patients with the variant had severe phenotypes, with 89% experiencing cardiac arrest.

CONCLUSIONS

This is the first example of a functionally validated rare noncoding variant at the SCN5A locus and highlights how genome sequencing in understudied populations can identify novel disease mechanisms. The variant partly explains the increased prevalence of BrS in this region and enables the identification of at-risk variant carriers to reduce the burden of sudden cardiac death in Thailand.

Characterization of molecular interactions between HDAC7 and MEF2A

Interactions of transcriptional corepressors such as histone deacetylase 7 (HDAC7), a class IIa HDAC, with myocyte enhancer factor-2 (MEF2) regulate MEF2 activity. Despite previous investigations exploring interactions between HDAC7 and MEF2, a detailed characterization of the HDAC7-MEF2 functional complex is still lacking. Herein, we first modeled the structure of the HDAC7-MEF2A complex and investigated the inter-protein interactions using all-atom molecular dynamics (MD) simulations. We identified specific amino acids within HDAC7 and MEF2A that participate in interactions such as salt bridges, hydrogen bonds, and hydrophobic interactions. Our results reveal a salt bridge formed between LYS96(HDAC7) and ASP63(MEF2A). Our analysis also predicted formations of reliable hydrogen bonds between SER82(HDAC7) and ASP63(MEF2A) as well as LYS96(HDAC7) and ASP63(MEF2A). In addition, clustering of hydrophobic residues at the interface contributes in stabilizing the HDAC7-MEF2A complex. Results from multiple sequence alignment show that most of the HDAC7 residues that are predicted to associate with MEF2A are conserved in at least three class IIa HDACs and all predicted residues in MEF2A are conserved in MEF2s. We also found that the association of DNA to MEF2A has no significant effect on HDAC7-MEF2A interactions. Our results may also provide useful insights into the interactions between other class IIa HDACs and MEF2s.

Discerning the Role of DNA Sequence, Shape, and Flexibility in Recognition by Drosophila Transcription Factors

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.

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.

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.

Noncanonical binding of transcription factors: time to revisit specificity?

Transcription factors (TFs) are one of the most studied classes of DNA-binding proteins that have a direct functional impact on gene transcription and thus, on human physiology and disease. The mechanisms that TFs use for recognizing target DNA binding sites have been studied for nearly five decades, yet they remain poorly understood. It is classically assumed that a TF recognizes a specific sequence pattern, or motif, as its binding sites. However, recent studies are consistently finding examples of noncanonical binding, that is, TFs binding at sites that do not resemble their sequence motifs. Here we review the current literature on four major types of noncanonical TF binding, namely binding based on DNA shape readout, at Guanine-quadruplex structures, at repeat sequences, and bispecific binding. These examples point to a critical need for studies to unify our current observations, many of which are at odds with the "one TF, one motif" view, into a more comprehensive definition of the DNA-binding specificity of TFs.

Cracking the Floral Quartet Code: How Do Multimers of MIKCC-Type MADS-Domain Transcription Factors Recognize Their Target Genes?

MADS-domain transcription factors (MTFs) are involved in the control of many important processes in eukaryotes. They are defined by the presence of a unique and highly conserved DNA-binding domain, the MADS domain. MTFs bind to double-stranded DNA as dimers and recognize specific sequences termed CArG boxes (such as 5'-CC(A/T)6GG-3') and similar sequences that occur hundreds of thousands of times in a typical flowering plant genome. The number of MTF-encoding genes increased by around two orders of magnitude during land plant evolution, resulting in roughly 100 genes in flowering plant genomes. This raises the question as to how dozens of different but highly similar MTFs accurately recognize the cis-regulatory elements of diverse target genes when the core binding sequence (CArG box) occurs at such a high frequency. Besides the usual processes, such as the base and shape readout of individual DNA sequences by dimers of MTFs, an important sublineage of MTFs in plants, termed MIKCC-type MTFs (MC-MTFs), has evolved an additional mechanism to increase the accurate recognition of target genes: the formation of heterotetramers of closely related proteins that bind to two CArG boxes on the same DNA strand involving DNA looping. MC-MTFs control important developmental processes in flowering plants, ranging from root and shoot to flower, fruit and seed development. The way in which MC-MTFs bind to DNA and select their target genes is hence not only of high biological interest, but also of great agronomic and economic importance. In this article, we review the interplay of the different mechanisms of target gene recognition, from the ordinary (base readout) via the extravagant (shape readout) to the idiosyncratic (recognition of the distance and orientation of two CArG boxes by heterotetramers of MC-MTFs). A special focus of our review is on the structural prerequisites of MC-MTFs that enable the specific recognition of target genes.

Structural predictions of protein-DNA binding: MELD-DNA

Structural, regulatory and enzymatic proteins interact with DNA to maintain a healthy and functional genome. Yet, our structural understanding of how proteins interact with DNA is limited. We present MELD-DNA, a novel computational approach to predict the structures of protein-DNA complexes. The method combines molecular dynamics simulations with general knowledge or experimental information through Bayesian inference. The physical model is sensitive to sequence-dependent properties and conformational changes required for binding, while information accelerates sampling of bound conformations. MELD-DNA can: (i) sample multiple binding modes; (ii) identify the preferred binding mode from the ensembles; and (iii) provide qualitative binding preferences between DNA sequences. We first assess performance on a dataset of 15 protein-DNA complexes and compare it with state-of-the-art methodologies. Furthermore, for three selected complexes, we show sequence dependence effects of binding in MELD predictions. We expect that the results presented herein, together with the freely available software, will impact structural biology (by complementing DNA structural databases) and molecular recognition (by bringing new insights into aspects governing protein-DNA interactions).

Physicochemical models of protein-DNA binding with standard and modified base pairs

DNA-binding proteins play important roles in various cellular processes, but the mechanisms by which proteins recognize genomic target sites remain incompletely understood. Functional groups at the edges of the base pairs (bp) exposed in the DNA grooves represent physicochemical signatures. As these signatures enable proteins to form specific contacts between protein residues and bp, their study can provide mechanistic insights into protein-DNA binding. Existing experimental methods, such as X-ray crystallography, can reveal such mechanisms based on physicochemical interactions between proteins and their DNA target sites. However, the low throughput of structural biology methods limits mechanistic insights for selection of many genomic sites. High-throughput binding assays enable prediction of potential target sites by determining relative binding affinities of a protein to massive numbers of DNA sequences. Many currently available computational methods are based on the sequence of standard Watson-Crick bp. They assume that the contribution of overall binding affinity is independent for each base pair, or alternatively include dinucleotides or short k-mers. These methods cannot directly expand to physicochemical contacts, and they are not suitable to apply to DNA modifications or non-Watson-Crick bp. These variations include DNA methylation, and synthetic or mismatched bp. The proposed method, DeepRec, can predict relative binding affinities as function of physicochemical signatures and the effect of DNA methylation or other chemical modifications on binding. Sequence-based modeling methods are in comparison a coarse-grain description and cannot achieve such insights. Our chemistry-based modeling framework provides a path towards understanding genome function at a mechanistic level.

Top-Down Crawl: a method for the ultra-rapid and motif-free alignment of sequences with associated binding metrics

SUMMARY

Several high-throughput protein-DNA binding methods currently available produce highly reproducible measurements of binding affinity at the level of the k-mer. However, understanding where a k-mer is positioned along a binding site sequence depends on alignment. Here, we present Top-Down Crawl (TDC), an ultra-rapid tool designed for the alignment of k-mer level data in a rank-dependent and position weight matrix (PWM)-independent manner. As the framework only depends on the rank of the input, the method can accept input from many types of experiments (protein binding microarray, SELEX-seq, SMiLE-seq, etc.) without the need for specialized parameterization. Measuring the performance of the alignment using multiple linear regression with 5-fold cross-validation, we find TDC to perform as well as or better than computationally expensive PWM-based methods.

AVAILABILITY AND IMPLEMENTATION

TDC can be run online at https://topdowncrawl.usc.edu or locally as a python package available through pip at https://pypi.org/project/TopDownCrawl.

CONTACT

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Trans-ethnic genome-wide association study of severe COVID-19

COVID-19 has caused numerous infections with diverse clinical symptoms. To identify human genetic variants contributing to the clinical development of COVID-19, we genotyped 1457 (598/859 with severe/mild symptoms) and sequenced 1141 (severe/mild: 474/667) patients of Chinese ancestry. We further incorporated 1401 genotyped and 948 sequenced ancestry-matched population controls, and tested genome-wide association on 1072 severe cases versus 3875 mild or population controls, followed by trans-ethnic meta-analysis with summary statistics of 3199 hospitalized cases and 897,488 population controls from the COVID-19 Host Genetics Initiative. We identified three significant signals outside the well-established 3p21.31 locus: an intronic variant in FOXP4-AS1 (rs1853837, odds ratio OR = 1.28, P = 2.51 × 10-10, allele frequencies in Chinese/European AF = 0.345/0.105), a frameshift insertion in ABO (rs8176719, OR = 1.19, P = 8.98 × 10-9, AF = 0.422/0.395) and a Chinese-specific intronic variant in MEF2B (rs74490654, OR = 8.73, P = 1.22 × 10-8, AF = 0.004/0). These findings highlight an important role of the adaptive immunity and the ABO blood-group system in protection from developing severe COVID-19.

SAMPDI-3D: predicting the effects of protein and DNA mutations on protein-DNA interactions

MOTIVATION

Mutations that alter protein-DNA interactions may be pathogenic and cause diseases. Therefore, it is extremely important to quantify the effect of mutations on protein-DNA binding free energy to reveal the molecular origin of diseases and to assist the development of treatments. Although several methods that predict the change of protein-DNA binding affinity upon mutations in the binding protein were developed, the effect of DNA mutations was not considered yet.

RESULTS

Here, we report a new version of SAMPDI, the SAMPDI-3D, which is a gradient boosting decision tree machine learning method to predict the change of the protein-DNA binding free energy caused by mutations in both the binding protein and the bases of the corresponding DNA. The method is shown to achieve Pearson correlation coefficient of 0.76 and 0.80 in a benchmarking test against experimentally determined change of the binding free energy caused by mutations in the binding protein or DNA, respectively. Furthermore, three datasets collected from literature were used to do blind benchmark for SAMPDI-3D and it is shown that it outperforms all existing state-of-the-art methods. The method is very fast allowing for genome-scale investigations.

AVAILABILITY

It is available as a web server and a stand-code at http://compbio.clemson.edu/SAMPDI-3D/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Targeted genomic investigations in a population-based cohort of mantle cell lymphoma reveal novel clinically relevant targets

Mantle cell lymphoma (MCL) is an aggressive B-cell neoplasm that follows a heterogeneous clinical course. Recurrent mutations have been described, but their applicability in the clinical setting is currently limited. The main reasons are challenges in the sequencing of DNA retrieved from formalin-fixed tissue commonly used for tissue collection in clinical biobanks. In this study, we sequenced 77 samples from a population-based de novo MCL cohort to investigate the utility of targeted sequencing in guiding personalized treatment approaches. Tumors were genetically variable, and a similar genetic landscape as previous studies using non-formalin fixed samples was identified, with recurrent mutations including ATM, KMT2D, and TP53. Novel alterations that can be considered actionable and/or indicative of treatment response were also identified. Our approach shows the potential benefits of using target sequencing of formalin fixed samples to facilitate treatment selection and individualized clinical decisions in MCL.

Effects of Structural Dynamics on Charge Carrier Transfer in B-DNA: A Combined MD and RT-TDDFT Study

Hole transfer along the axis of duplex DNA has been the focus of physical chemistry research for decades, with implications in diverse fields, from nanotechnology to cell oxidative damage. Computational approaches are particularly amenable for this problem, to complement experimental data for interpretation of transfer mechanisms. To be predictive, computational results need to account for the inherent mobility of biological molecules during the time frame of experimental measurements. Here, we address the structural variability of B-DNA and its effects on hole transfer in a combined molecular dynamics (MD) and real-time time-dependent density functional theory (RT-TDDFT) study. Our results show that quantities that characterize the charge transfer process, such as the time-dependent dipole moment and hole population at a specific site, are sensitive to structural changes that occur on the nanosecond time scale. We extend the range of physical properties for which such a correlation has been observed, further establishing the fact that quantitative computational data on charge transfer properties should include statistical averages. Furthermore, we use the RT-TDDFT results to assess an efficient tight-binding method suitable for high-throughput predictions. We demonstrate that charge transfer, although affected by structural variability, on average, remains strong in AA and GG dimers.

MEF2C shapes the microtranscriptome during differentiation of skeletal muscles

Myocyte enhancer factor 2C (MEF2C) is a transcription factor that regulates heart and skeletal muscle differentiation and growth. Several protein-encoding genes were identified as targets of this factor; however, little is known about its contribution to the microtranscriptome composition and dynamics in myogenic programs. In this report, we aimed to address this question. Deep sequencing of small RNAs of human muscle cells revealed a set of microRNAs (miRNAs), including several muscle-specific miRNAs, that are sensitive to MEF2C depletion. As expected, in cells with knockdown of MEF2C, we found mostly downregulated miRNAs; nevertheless, as much as one-third of altered miRNAs were upregulated. The majority of these changes are driven by transcription efficiency. Moreover, we found that MEF2C affects nontemplated 3′-end nucleotide addition of miRNAs, mainly oligouridylation. The rate of these modifications is associated with the level of TUT4 which mediates RNA 3′-uridylation. Finally, we found that a quarter of miRNAs which significantly changed upon differentiation of human skeletal myoblasts is inversely altered in MEF2C deficient cells. We concluded that MEF2C is an essential factor regulating both the quantity and quality of the microtranscriptome, leaving an imprint on the stability and perhaps specificity of many miRNAs during the differentiation of muscle cells.