Structural basis for the regulation of plant transcription factor WRKY33 by the VQ protein SIB1.

The WRKY transcription factors play essential roles in a variety of plant signaling pathways associated with biotic and abiotic stress response. The transcriptional activity of many WRKY members are regulated by a class of intrinsically disordered VQ proteins. While it is known that VQ proteins interact with the WRKY DNA-binding domains (DBDs), also termed as the WRKY domains, structural information regarding VQ-WRKY interaction is lacking and the regulation mechanism remains unknown. Herein we report a solution NMR study of the interaction between Arabidopsis WRKY33 and its regulatory VQ protein partner SIB1. We uncover a SIB1 minimal sequence neccessary for forming a stable complex with WRKY33 DBD, which comprises not only the consensus "FxxhVQxhTG" VQ motif but also its preceding region. We demonstrate that the ßN-strand and the extended ßN-ß1 loop of WRKY33 DBD form the SIB1 docking site, and build a structural model of the complex based on the NMR paramagnetic relaxation enhancement and mutagenesis data. Based on this model, we further identify a cluster of positively-charged residues in the N-terminal region of SIB1 to be essential for the formation of a SIB1-WRKY33-DNA ternary complex. These results provide a framework for the mechanism of SIB1-enhanced WRKY33 transcriptional activity.

WW domains: A globular NLS recognized by importin-α.

In this issue, the discovery by Yang et al. (https://doi.org/10.1083/jcb.202308013) that folded WW domains of YAP1 and other proteins bind to Impα introduces a new class of globular NLS, contrasting with the extensively studied linear NLS motifs. This finding underscores the versatility of importins in recognizing their cargo proteins.

Structural and functional insights into transcription activation of the essential LysR-type transcriptional regulators.

The enormous LysR-type transcriptional regulators (LTTRs), which are diversely distributed amongst prokaryotes, play crucial roles in transcription regulation of genes involved in basic metabolic pathways, virulence and stress resistance. However, the precise transcription activation mechanism of these genes by LTTRs remains to be explored. Here, we determine the cryo-EM structure of a LTTR-dependent transcription activation complex comprising of Escherichia coli RNA polymerase (RNAP), an essential LTTR protein GcvA and its cognate promoter DNA. Structural analysis shows two N-terminal DNA binding domains of GcvA (GcvA_DBD) dimerize and engage the GcvA activation binding sites, presenting the -35 element for specific recognition with the conserved σ70R4. In particular, the versatile C-terminal domain of α subunit of RNAP directly interconnects with GcvA_DBD, σ70R4 and promoter DNA, providing more interfaces for stabilizing the complex. Moreover, molecular docking supports glycine as one potential inducer of GcvA, and single molecule photobleaching experiments kinetically visualize the occurrence of tetrameric GcvA-engaged transcription activation complex as suggested for the other LTTR homologs. Thus, a general model for tetrameric LTTR-dependent transcription activation is proposed. These findings will provide new structural and functional insights into transcription activation of the essential LTTRs.

A finely balanced order-disorder equilibrium sculpts the folding-binding landscape of an antibiotic sequestering protein.

TipA, a MerR family transcription factor from Streptomyces lividans, promotes antibiotic resistance by sequestering broad-spectrum thiopeptide-based antibiotics, thus counteracting their inhibitory effect on ribosomes. TipAS, a minimal binding motif which is expressed as an isoform of TipA, harbors a partially disordered N-terminal subdomain that folds upon binding multiple antibiotics. The extent and nature of the underlying molecular heterogeneity in TipAS that shapes its promiscuous folding-function landscape is an open question and is critical for understanding antibiotic-sequestration mechanisms. Here, combining equilibrium and time-resolved experiments, statistical modeling, and simulations, we show that the TipAS native ensemble exhibits a pre-equilibrium between binding-incompetent and binding-competent substates, with the fully folded state appearing only as an excited state under physiological conditions. The binding-competent state characterized by a partially structured N-terminal subdomain loses structure progressively in the physiological range of temperatures, swells on temperature increase, and displays slow conformational exchange across multiple conformations. Binding to the bactericidal antibiotic thiostrepton follows a combination of induced-fit and conformational-selection-like mechanisms, via partial binding and concomitant stabilization of the binding-competent substate. These ensemble features are evolutionarily conserved across orthologs from select bacteria that infect humans, underscoring the functional role of partial disorder in the native ensemble of antibiotic-sequestering proteins belonging to the MerR family.

Entamoeba histolytica: In Silico and In Vitro Oligomerization of EhHSTF5 Enhances Its Binding to the HSE of the EhPgp5 Gene Promoter.

Throughout its lifecycle, Entamoeba histolytica encounters a variety of stressful conditions. This parasite possesses Heat Shock Response Elements (HSEs) which are crucial for regulating the expression of various genes, aiding in its adaptation and survival. These HSEs are regulated by Heat Shock Transcription Factors (EhHSTFs). Our research has identified seven such factors in the parasite, designated as EhHSTF1 through to EhHSTF7. Significantly, under heat shock conditions and in the presence of the antiamoebic compound emetine, EhHSTF5, EhHSTF6, and EhHSTF7 show overexpression, highlighting their essential role in gene response to these stressors. Currently, only EhHSTF7 has been confirmed to recognize the HSE as a promoter of the EhPgp5 gene (HSE_EhPgp5), leaving the binding potential of the other EhHSTFs to HSEs yet to be explored. Consequently, our study aimed to examine, both in vitro and in silico, the oligomerization, and binding capabilities of the recombinant EhHSTF5 protein (rEhHSTF5) to HSE_EhPgp5. The in vitro results indicate that the oligomerization of rEhHSTF5 is concentration-dependent, with its dimeric conformation showing a higher affinity for HSE_EhPgp5 than its monomeric state. In silico analysis suggests that the alpha 3 α-helix (α3-helix) of the DNA-binding domain (DBD5) of EhHSTF5 is crucial in binding to the major groove of HSE, primarily through hydrogen bonding and salt-bridge interactions. In summary, our results highlight the importance of oligomerization in enhancing the affinity of rEhHSTF5 for HSE_EhPgp5 and demonstrate its ability to specifically recognize structural motifs within HSE_EhPgp5. These insights significantly contribute to our understanding of one of the potential molecular mechanisms employed by this parasite to efficiently respond to various stressors, thereby enabling successful adaptation and survival within its host environment.

Characterization of the transcriptionally active form of dephosphorylated DctD complexed with dephospho-IIAGlc.

Bacterial enhancer-binding proteins (bEBPs) acquire a transcriptionally active state via phosphorylation. However, transcriptional activation by the dephosphorylated form of bEBP has been observed in DctD, which belongs to Group I bEBP. The formation of a complex between dephosphorylated DctD (d-DctD) and dephosphorylated IIAGlc (d-IIAGlc) is a prerequisite for the transcriptional activity of d-DctD. In the present study, characteristics of the transcriptionally active complex composed of d-IIAGlc and phosphorylation-deficient DctD (DctDD57Q) of Vibrio vulnificus were investigated in its multimeric conformation and DNA-binding ability. DctDD57Q formed a homodimer that could not bind to the DNA. In contrast, when DctDD57Q formed a complex with d-IIAGlc in a 1:1 molar ratio, it produced two conformations: dimer and dodecamer of the complex. Only the dodecameric complex exhibited ATP-hydrolyzing activity and DNA-binding affinity. For successful DNA-binding and transcriptional activation by the dodecameric d-IIAGlc/DctDD57Q complex, extended upstream activator sequences were required, which encompass the nucleotide sequences homologous to the known DctD-binding site and additional nucleotides downstream. This is the first report to demonstrate the molecular characteristics of a dephosphorylated bEBP complexed with another protein to form a transcriptionally active dodecameric complex, which has an affinity for a specific DNA-binding sequence.IMPORTANCEResponse regulators belonging to the bacterial two-component regulatory system activate the transcription initiation of their regulons when they are phosphorylated by cognate sensor kinases and oligomerized to the appropriate multimeric states. Recently, it has been shown that a dephosphorylated response regulator, DctD, could activate transcription in a phosphorylation-independent manner in Vibrio vulnificus. The dephosphorylated DctD activated transcription as efficiently as phosphorylated DctD when it formed a complex with dephosphorylated form of IIAGlc, a component of the glucose-phosphotransferase system. Functional mimicry of this complex with the typical form of transcriptionally active phosphorylated DctD led us to study the molecular characteristics of this heterodimeric complex. Through systematic analyses, it was surprisingly determined that a multimer constituted with 12 complexes gained the ability to hydrolyze ATP and recognize specific upstream activator sequences containing a typical inverted-repeat sequence flanked by distinct nucleotides.

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.

mRNA Display Identifies Potent, Paralog-Selective Peptidic Ligands for ARID1B.

The ARID1A and ARID1B subunits are mutually exclusive components of the BAF variant of SWI/SNF chromatin remodeling complexes. Loss of function mutations in ARID1A are frequently observed in various cancers, resulting in a dependency on the paralog ARID1B for cancer cell proliferation. However, ARID1B has never been targeted directly, and the high degree of sequence similarity to ARID1A poses a challenge for the development of selective binders. In this study, we used mRNA display to identify peptidic ligands that bind with nanomolar affinities to ARID1B and showed high selectivity over ARID1A. Using orthogonal biochemical, biophysical, and chemical biology tools, we demonstrate that the peptides engage two different binding pockets, one of which directly involves an ARID1B-exclusive cysteine that could allow covalent targeting by small molecules. Our findings impart the first evidence of the ligandability of ARID1B, provide valuable tools for drug discovery, and suggest opportunities for the development of selective molecules to exploit the synthetic lethal relationship between ARID1A and ARID1B in cancer.

Binding Mechanism of Inhibitors to BRD4 and BRD9 Decoded by Multiple Independent Molecular Dynamics Simulations and Deep Learning.

Bromodomain 4 and 9 (BRD4 and BRD9) have been regarded as important targets of drug designs in regard to the treatment of multiple diseases. In our current study, molecular dynamics (MD) simulations, deep learning (DL) and binding free energy calculations are integrated to probe the binding modes of three inhibitors (H1B, JQ1 and TVU) to BRD4 and BRD9. The MD trajectory-based DL successfully identify significant functional function domains, such as BC-loop and ZA-loop. The information from the post-processing analysis of MD simulations indicates that inhibitor binding highly influences the structural flexibility and dynamic behavior of BRD4 and BRD9. The results of the MM-GBSA calculations not only suggest that the binding ability of H1B, JQ1 and TVU to BRD9 are stronger than to BRD4, but they also verify that van der Walls interactions are the primary forces responsible for inhibitor binding. The hot spots of BRD4 and BRD9 revealed by residue-based free energy estimation provide target sites of drug design in regard to BRD4 and BRD9. This work is anticipated to provide useful theoretical aids for the development of selective inhibitors over BRD family members.

Structural basis of Integrator-dependent RNA polymerase II termination.

The Integrator complex can terminate RNA polymerase II (Pol II) in the promoter-proximal region of genes. Previous work has shed light on how Integrator binds to the paused elongation complex consisting of Pol II, the DRB sensitivity-inducing factor (DSIF) and the negative elongation factor (NELF) and how it cleaves the nascent RNA transcript1, but has not explained how Integrator removes Pol II from the DNA template. Here we present three cryo-electron microscopy structures of the complete Integrator-PP2A complex in different functional states. The structure of the pre-termination complex reveals a previously unresolved, scorpion-tail-shaped INTS10-INTS13-INTS14-INTS15 module that may use its 'sting' to open the DSIF DNA clamp and facilitate termination. The structure of the post-termination complex shows that the previously unresolved subunit INTS3 and associated sensor of single-stranded DNA complex (SOSS) factors prevent Pol II rebinding to Integrator after termination. The structure of the free Integrator-PP2A complex in an inactive closed conformation2 reveals that INTS6 blocks the PP2A phosphatase active site. These results lead to a model for how Integrator terminates Pol II transcription in three steps that involve major rearrangements.

Structural homology-based identification of BEN domain proteins in Poxviruses.

BEND family transcription factors directly interact with DNA through BEN domains and have been found across metazoan species. Interestingly, certain insect and mammalian viruses have also hijacked Bend genes into their genome. However, the phylogenetic classification and evolution of these viral BEN domains remain unclear. Building on our previous finding that in silico method accurately determine the 3D model of BEN domains, we used AlphaFold2 to predict the tertiary structures of poxviral BEN domains for comprehensive homologous comparison. We revealed that the majority of poxviral BEN modules exhibit characteristics of type II BEN. Additionally, electrostatic surface potential analysis found various poxviral BEN domains, including the first BEN of OPG067 in Orthopoxvirus, the third BEN of OPG067 in Yatapoxvirus and the third BEN of MC036R in MCV, have positively charged protein surfaces, indicating a structural basis for DNA loading. Notably, MC036R shares structural resemblance with human BEND3, as they both contain four BEN domains and an intrinsically disordered region. In summary, our discoveries provide deeper insights into the functional roles of BEN proteins within poxviruses.

Allosteric Modulation of the YAP/TAZ-TEAD Interaction by Palmitoylation and Small-Molecule Inhibitors.

The Hippo signaling pathway is a highly conserved signaling network that plays a central role in regulating cellular growth, proliferation, and organ size. This pathway consists of a kinase cascade that integrates various upstream signals to control the activation or inactivation of YAP/TAZ proteins. Phosphorylated YAP/TAZ is sequestered in the cytoplasm; however, when the Hippo pathway is deactivated, it translocates into the nucleus, where it associates with TEAD transcription factors. This partnership is instrumental in regulating the transcription of progrowth and antiapoptotic genes. Thus, in many cancers, aberrantly hyperactivated YAP/TAZ promotes oncogenesis by contributing to cancer cell proliferation, metastasis, and therapy resistance. Because YAP and TAZ exert their oncogenic effects by binding with TEAD, it is critical to understand this key interaction to develop cancer therapeutics. Previous research has indicated that TEAD undergoes autopalmitoylation at a conserved cysteine, and small molecules that inhibit TEAD palmitoylation disrupt effective YAP/TAZ binding. However, how exactly palmitoylation contributes to YAP/TAZ-TEAD interactions and how the TEAD palmitoylation inhibitors disrupt this interaction remains unknown. Utilizing molecular dynamics simulations, our investigation not only provides detailed atomistic insight into the YAP/TAZ-TEAD dynamics but also unveils that the inhibitor studied influences the binding of YAP and TAZ to TEAD in distinct manners. This discovery has significant implications for the design and deployment of future molecular interventions targeting this interaction.

Peptide Mimics of the Cysteine-Rich Regions of HapX and SreA Bind a [2Fe-2S] Cluster In Vitro.

HapX and SreA are transcription factors that regulate the response of the fungus Aspergillus fumigatus to the availability of iron. During iron starvation, HapX represses genes involved in iron consuming pathways and upon a shift to iron excess, HapX activates these same genes. SreA blocks the expression of genes needed for iron uptake during periods of iron availability. Both proteins possess cysteine-rich regions (CRR) that are hypothesized to be necessary for the sensing of iron levels. However, the contribution of each of these domains to the function of the protein has remained unclear. Here, the ability of peptide analogs of each CRR is determined to bind an iron-sulfur cluster in vitro. UV-vis and resonance Raman (RR) spectroscopies reveal that each CRR is capable of coordinating a [2Fe-2S] cluster with comparable affinities. The iron-sulfur cluster coordinated to the CRR-B domain of HapX displays particularly high stability. The data are consistent with HapX and SreA mediating responses to cellular iron levels through the direct coordination of [2Fe-2S] clusters. The high stability of the CRR-B peptide may also find use as a starting point for the development of new green catalysts.

Cancer-associated polybromo-1 bromodomain 4 missense variants variably impact bromodomain ligand binding and cell growth suppression.

The polybromo, brahma-related gene 1-associated factors (PBAF) chromatin remodeling complex subunit polybromo-1 (PBRM1) contains six bromodomains that recognize and bind acetylated lysine residues on histone tails and other nuclear proteins. PBRM1 bromodomains thus provide a link between epigenetic posttranslational modifications and PBAF modulation of chromatin accessibility and transcription. As a putative tumor suppressor in several cancers, PBRM1 protein expression is often abrogated by truncations and deletions. However, ∼33% of PBRM1 mutations in cancer are missense and cluster within its bromodomains. Such mutations may generate full-length PBRM1 variant proteins with undetermined structural and functional characteristics. Here, we employed computational, biophysical, and cellular assays to interrogate the effects of PBRM1 bromodomain missense variants on bromodomain stability and function. Since mutations in the fourth bromodomain of PBRM1 (PBRM1-BD4) comprise nearly 20% of all cancer-associated PBRM1 missense mutations, we focused our analysis on PBRM1-BD4 missense protein variants. Selecting 16 potentially deleterious PBRM1-BD4 missense protein variants for further study based on high residue mutational frequency and/or conservation, we show that cancer-associated PBRM1-BD4 missense variants exhibit varied bromodomain stability and ability to bind acetylated histones. Our results demonstrate the effectiveness of identifying the unique impacts of individual PBRM1-BD4 missense variants on protein structure and function, based on affected residue location within the bromodomain. This knowledge provides a foundation for drawing correlations between specific cancer-associated PBRM1 missense variants and distinct alterations in PBRM1 function, informing future cancer personalized medicine approaches.

CHCHD4 binding affects the active site of apoptosis inducing factor (AIF): Structural determinants for allosteric regulation.

Apoptosis-inducing factor (AIF), which is confined to mitochondria of normal healthy cells, is the first identified caspase-independent cell death effector. Moreover, AIF is required for the optimal functioning of the respiratory chain machinery. Recent findings have revealed that AIF fulfills its pro-survival function by interacting with CHCHD4, a soluble mitochondrial protein which promotes the entrance and the oxidative folding of different proteins in the inner membrane space. Here, we report the crystal structure of the ternary complex involving the N-terminal 27-mer peptide of CHCHD4, NAD+, and AIF harboring its FAD (flavin adenine dinucleotide) prosthetic group in oxidized form. Combining this information with biophysical and biochemical data on the CHCHD4/AIF complex, we provide a detailed structural description of the interaction between the two proteins, validated by both chemical cross-linking mass spectrometry analysis and site-directed mutagenesis.

Single-molecule tracking reveals the functional allocation, in vivo interactions, and spatial organization of universal transcription factor NusG.

During transcription elongation, NusG aids RNA polymerase by inhibiting pausing, promoting anti-termination on rRNA operons, coupling transcription with translation on mRNA genes, and facilitating Rho-dependent termination. Despite extensive work, the in vivo functional allocation and spatial distribution of NusG remain unknown. Using single-molecule tracking and super-resolution imaging in live E. coli cells, we found NusG predominantly in a chromosome-associated population (binding to RNA polymerase in elongation complexes) and a slowly diffusing population complexed with the 30S ribosomal subunit; the latter provides a "30S-guided" path for NusG into transcription elongation. Only ∼10% of NusG is fast diffusing, with its mobility suggesting non-specific interactions with DNA for >50% of the time. Antibiotic treatments and deletion mutants revealed that chromosome-associated NusG participates mainly in rrn anti-termination within phase-separated transcriptional condensates and in transcription-translation coupling. This study illuminates the multiple roles of NusG and offers a guide on dissecting multi-functional machines via in vivo imaging.

DNA binding redistributes activation domain ensemble and accessibility in pioneer factor Sox2.

More than 1600 human transcription factors orchestrate the transcriptional machinery to control gene expression and cell fate. Their function is conveyed through intrinsically disordered regions (IDRs) containing activation or repression domains but lacking quantitative structural ensemble models prevents their mechanistic decoding. Here we integrate single-molecule FRET and NMR spectroscopy with molecular simulations showing that DNA binding can lead to complex changes in the IDR ensemble and accessibility. The C-terminal IDR of pioneer factor Sox2 is highly disordered but its conformational dynamics are guided by weak and dynamic charge interactions with the folded DNA binding domain. Both DNA and nucleosome binding induce major rearrangements in the IDR ensemble without affecting DNA binding affinity. Remarkably, interdomain interactions are redistributed in complex with DNA leading to variable exposure of two activation domains critical for transcription. Charged intramolecular interactions allowing for dynamic redistributions may be common in transcription factors and necessary for sensitive tuning of structural ensembles.

PML Nuclear bodies: the cancer connection and beyond.

Promyelocytic leukemia (PML) nuclear bodies, membrane-less organelles in the nucleus, play a crucial role in cellular homeostasis. These dynamic structures result from the assembly of scaffolding PML proteins and various partners. Recent crystal structure analyses revealed essential self-interacting domains, while liquid-liquid phase separation contributes to their formation. PML bodies orchestrate post-translational modifications, particularly stress-induced SUMOylation, impacting target protein functions. Serving as hubs in multiple signaling pathways, they influence cellular processes like senescence. Dysregulation of PML expression contributes to diseases, including cancer, highlighting their significance. Therapeutically, PML bodies are promising targets, exemplified by successful acute promyelocytic leukemia treatment with arsenic trioxide and retinoic acid restoring PML bodies. Understanding their functions illuminates both normal and pathological cellular physiology, guiding potential therapies. This review explores recent advancements in PML body biogenesis, biochemical activity, and their evolving biological roles.

Optimization of the synthesis of BET BD2 selective inhibitor XY153.

XY153 is a promising BET BD2 inhibitor with an IC50 value of 0.79 nM against BRD4 BD2. It shows 354-fold selectivity over BRD4-BD1 and 6-fold selectivity over other BET BD2 domains. However, the reported synthesis route of XY153 and its derivatives are extremely poor-yielding. After the synthesis of three key fragments, XY153 can only be obtained with a yield of 1.3 % in the original four-step reaction. In this study, we reported a three-step alternative route in the synthesis process of XY153. The reaction conditions for this route were thoroughly investigated and optimized, resulting in a significantly improved yield of 61.5 %. This efficient synthesis route establishes a robust chemical foundation for the rapid synthesis of XY153 derivatives as BET BD2 inhibitors in the near future.

Mechanistic and structural insight into R2R3-MYB transcription factor in plants: molecular dynamics based binding free energy analysis.

The MYB transcription factor (TF) family is essential for various plant growth and development processes, including responses to biotic and abiotic stresses. This study investigated the R2R3-MYB protein structure from five plants, including cereal crops. The R2R3-MYB protein structure was docked with the DNA structure, and the best complexes were selected for two runs of molecular dynamics (MD) simulations to investigate the key interacting residues and conformational changes in the R2R3-MYB proteins caused by DNA binding. The MM/PBSA method calculated the binding free energy for each R2R3-MYB protein-DNA complex, showing strong interaction. Hydrophobic and hydrogen bonds significantly stabilized the R2R3-MYB protein-DNA complexes. The principal component analysis showed high restrictions on the movement of protein atoms in the phase space. A similar MD simulation analysis was performed using the crystal structure of the R2R3-MYB protein-DNA complex from Arabidopsis thaliana, and the generated complexes resembled the X-ray crystal structure. This is the first detailed study on the R2R3-MYB protein-DNA complex in cereal crops, providing a cost-effective solution to identify the key interacting residues and analyze the conformational changes in the MYB domain before and after DNA binding.Communicated by Ramaswamy H. Sarma.