Sequence-Specific Structural Features and Solvation Properties of Transcription Factor Binding DNA Motifs: Insights from Molecular Dynamics Simulation

Molecular Interactions of the Pioneer Transcription Factor GATA3 With DNA

The GATA3 transcription factor is a pioneer transcription factor that is critical in the development, proliferation, and maintenance of several immune cell types. Identifying the detailed conformational dynamics and interactions of this transcription factor, as well as its clinically important population variants will allow us to unravel its mode of action. In this study, we analyze the molecular interactions of the GATA3 transcription factor bound to dsDNA as well as three clinically important population variants by atomistic molecular dynamics simulations. We identify the effect of the variants on the DNA conformational dynamics and delineate the differences compared to the wildtype transcription factor that could be related to impaired function. We highlight the structural plasticity in the binding of the GATA3 transcription factor and identify important DNA-protein contacts. Although the DNA-protein contacts are persistent and appear to be stable, they exhibit nanosecond timescale fluctuations and several binding/unbinding events. Further, we identify differential DNA binding in the three variants and show that the N-terminal binding is reduced in two of the variants. Our results indicate that reduced minor groove width and DNA diameter are important hallmarks for the binding of GATA3. Our work is an important step towards understanding the functional dynamics of the GATA3 protein and its clinically significant population variants.

Identifying the "stripe" transcription factors and cooperative binding related to DNA methylation

DNA methylation plays a critical role in gene regulation by modulating the DNA binding of transcription factors (TFs). This study integrates TFs' ChIP-seq profiles with WGBS profiles to investigate how DNA methylation affects protein interactions. Statistical methods and a 5-letter DNA motif calling model have been developed to characterize DNA sequences bound by proteins, while considering the effects of DNA modifications. By employing these methods, 79 significant universal "stripe" TFs and cofactors (USFs), 2360 co-binding protein pairs, and distinct protein modules associated with various DNA methylation states have been identified. The USFs hint a regulatory hierarchy within these protein interactions. Proteins preferentially bind to non-CpG sites in methylated regions, indicating binding affinity is not solely CpG-dependent. Proteins involved in methylation-specific USFs and cobinding pairs play essential roles in promoting and sustaining DNA methylation through interacting with DNMTs or inhibiting TET binding. These findings underscore the interplay between protein binding and methylation, offering insights into epigenetic regulation in cellular biology.

FLRT3 and TGF-β/SMAD4 signalling: Impacts on apoptosis, autophagy and ion channels in supraventricular tachycardia

To explore the underlying molecular mechanisms of supraventricular tachycardia (SVT), this study aimed to analyse the complex relationship between FLRT3 and TGF-β/SMAD4 signalling pathway, which affects Na+ and K+ channels in cardiomyocytes. Bioinformatics analysis was performed on 85 SVT samples and 15 healthy controls to screen overlapping genes from the key module and differentially expressed genes (DEGs). Expression profiling of overlapping genes, coupled with Receiver Operating Characteristic (ROC) curve analyses, identified FLRT3 as a hub gene. In vitro studies utilizing Ang II-stimulated H9C2 cardiomyocytes were undertaken to elucidate the consequences of FLRT3 silencing on cardiomyocyte apoptosis and autophagic processes. Utilizing a combination of techniques such as quantitative reverse-transcription polymerase chain reaction (qRT-PCR), western blotting (WB), flow cytometry, dual-luciferase reporter assays and chromatin immunoprecipitation polymerase chain reaction (ChIP-PCR) assays were conducted to decipher the intricate interactions between FLRT3, the TGF-β/SMAD4 signalling cascade and ion channel gene expression. Six genes (AADAC, DSC3, FLRT3, SYT4, PRR9 and SERTM1) demonstrated reduced expression in SVT samples, each possessing significant clinical diagnostic potential. In H9C2 cardiomyocytes, FLRT3 silencing mitigated Ang II-induced apoptosis and modulated autophagy. With increasing TGF-β concentration, there was a dose-responsive decline in FLRT3 and SCN5A expression, while both KCNIP2 and KCND2 expressions were augmented. Moreover, a direct interaction between FLRT3 and SMAD4 was observed, and inhibition of SMAD4 expression resulted in increased FLRT3 expression. Our results demonstrated that the TGF-β/SMAD4 signalling pathway plays a critical role by regulating FLRT3 expression, with potential implications for ion channel function in SVT.

Structural and dynamical aspect of DNA motif sequence specific binding of AP-1 transcription factor

Activator protein-1 (AP-1) comprises one of the largest and most evolutionary conserved families of ubiquitous eukaryotic transcription factors that act as a pioneer factor. Diversity in DNA binding interaction of AP-1 through a conserved basic-zipper (bZIP) domain directs in-depth understanding of how AP-1 achieves its DNA binding selectivity and consequently gene regulation specificity. Here, we address the structural and dynamical aspects of the DNA target recognition process of AP-1 using microsecond-long atomistic simulations based on the structure of the human AP-1 FosB/JunD bZIP-DNA complex. Our results show the unique role of DNA shape features in selective base specific interactions, characteristic ion population, and solvation properties of DNA grooves to form the motif sequence specific AP-1-DNA complex. The TpG step at the two terminals of the AP-1 site plays an important role in the structural adjustment of DNA by modifying the helical twist in the AP-1 bound state. We addressed the role of intrinsic motion of the bZIP domain in terms of opening and closing gripper motions of DNA binding helices, in target site recognition and binding of AP-1 factors. Our observations suggest that binding to the cognate motif in DNA is mainly accompanied with the precise adjustment of closing gripper motion of DNA binding helices of the bZIP domain.

Domain- and state-specific shape of the electric field tunes voltage sensing in voltage-gated sodium channels

The ability to sense transmembrane voltage underlies most physiological roles of voltage-gated sodium (Nav) channels. Whereas the key role of their voltage-sensing domains (VSDs) in channel activation is well established, the molecular underpinnings of voltage coupling remain incompletely understood. Voltage-dependent energetics of the activation process can be described in terms of the gating charge that is defined by coupling of charged residues to the external electric field. The shape of the electric field within VSDs is therefore crucial for the activation of voltage-gated ion channels. Here, we employed molecular dynamics simulations of cardiac Nav1.5 and bacterial NavAb, together with our recently developed tool g_elpot, to gain insights into the voltage-sensing mechanisms of Nav channels via high-resolution quantification of VSD electrostatics. In contrast to earlier low-resolution studies, we found that the electric field within VSDs of Nav channels has a complex isoform- and domain-specific shape, which prominently depends on the activation state of a VSD. Different VSDs vary not only in the length of the region where the electric field is focused but also differ in their overall electrostatics, with possible implications in the diverse ion selectivity of their gating pores. Due to state-dependent field reshaping, not only translocated basic but also relatively immobile acidic residues contribute significantly to the gating charge. In the case of NavAb, we found that the transition between structurally resolved activated and resting states results in a gating charge of 8e, which is noticeably lower than experimental estimates. Based on the analysis of VSD electrostatics in the two activation states, we propose that the VSD likely adopts a deeper resting state upon hyperpolarization. In conclusion, our results provide an atomic-level description of the gating charge, demonstrate diversity in VSD electrostatics, and reveal the importance of electric-field reshaping for voltage sensing in Nav channels.