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.

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.

Intrinsically disordered regions of the Msn2 transcription factor encode multiple functions using interwoven sequence grammars

Intrinsically disordered regions (IDRs) are abundant in eukaryotic proteins, but their sequence-function relationship remains poorly understood. IDRs of transcription factors (TFs) can direct promoter selection and recruit coactivators, as shown for the budding yeast TF Msn2. To examine how IDRs encode both these functions, we compared genomic binding specificity, coactivator recruitment, and gene induction amongst a large set of designed Msn2-IDR mutants. We find that both functions depend on multiple regions across the > 600AA IDR. Yet, transcription activity was readily disrupted by mutations that showed no effect on the Msn2 binding specificity. Our data attribute this differential sensitivity to the integration of a relaxed, composition-based code directing binding specificity with a more stringent, motif-based code controlling the recruitment of coactivators and transcription activity. Therefore, Msn2 utilizes interwoven sequence grammars for encoding multiple functions, suggesting a new IDR design paradigm of potentially general use.

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.

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.

Synthesis of a 13C-methylene-labeled isoleucine precursor as a useful tool for studying protein side-chain interactions and dynamics

In this study, we present the synthesis and incorporation of a metabolic isoleucine precursor compound for selective methylene labeling. The utility of this novel α-ketoacid isotopologue is shown by incorporation into the protein Brd4-BD1, which regulates gene expression by binding to acetylated histones. High quality single quantum 13C-1 H-HSQC were obtained, as well as triple quantum HTQC spectra, which are superior in terms of significantly increased 13C-T2 times. Additionally, large chemical shift perturbations upon ligand binding were observed. Our study thus proves the great sensitivity of this precursor as a reporter for side-chain dynamic studies and for investigations of CH-π interactions in protein-ligand complexes.

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.

Structure of the WYL-domain containing transcription activator, DriD, in complex with ssDNA effector and DNA target site

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

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.

IraM remodels the RssB segmented helical linker to stabilize σs against degradation by ClpXP

Upon Mg2+ starvation, a condition often associated with virulence, enterobacteria inhibit the ClpXP-dependent proteolysis of the master transcriptional regulator, σs, via IraM, a poorly understood antiadaptor that prevents RssB-dependent loading of σs onto ClpXP. This inhibition results in σs accumulation and expression of stress resistance genes. Here, we report on the structural analysis of RssB bound to IraM, which reveals that IraM induces two folding transitions within RssB, amplified via a segmented helical linker. These conformational changes result in an open, yet inhibited RssB structure in which IraM associates with both the C-terminal and N-terminal domains of RssB and prevents binding of σs to the 4-5-5 face of the N-terminal receiver domain. This work highlights the remarkable structural plasticity of RssB and reveals how a stress-specific RssB antagonist modulates a core stress response pathway that could be leveraged to control biofilm formation, virulence, and the development of antibiotic resistance.

The molecular structure and function of fibrocystin, the key gene product implicated in autosomal recessive polycystic kidney disease (ARPKD)

Autosomal recessive polycystic kidney disease is an early onset inherited hepatorenal disorder affecting around 1 in 20,000 births with no approved specific therapies. The disease is almost always caused by variations in the polycystic kidney and hepatic disease 1 gene, which encodes fibrocystin (FC), a very large, single-pass transmembrane glycoprotein found in primary cilia, urine and urinary exosomes. By comparison to proteins involved in autosomal dominant PKD, our structural and molecular understanding of FC has lagged far behind such that there are no published experimentally determined structures of any part of the protein. Bioinformatics analyses predict that the ectodomain contains a long chain of immunoglobulin-like plexin-transcription factor domains, a protective antigen 14 domain, a tandem G8-TMEM2 homology region and a sperm protein, enterokinase and agrin domain. Here we review current knowledge on the molecular function of the protein from a structural perspective.

Improving Enhancer Identification with a Multi-Classifier Stacked Ensemble Model

Enhancers are DNA regions that are responsible for controlling the expression of genes. Enhancers are usually found upstream or downstream of a gene, or even inside a gene's intron region, but are normally located at a distant location from the genes they control. By integrating experimental and computational approaches, it is possible to uncover enhancers within DNA sequences, which possess regulatory properties. Experimental techniques such as ChIP-seq and ATAC-seq can identify genomic regions that are associated with transcription factors or accessible to regulatory proteins. On the other hand, computational techniques can predict enhancers based on sequence features and epigenetic modifications. In our study, we have developed a multi-classifier stacked ensemble (MCSE-enhancer) model that can accurately identify enhancers. We utilized feature descriptors from various physiochemical properties as input for our six baseline classifiers and built a stacked classifier, which outperformed previous enhancer classification techniques in terms of accuracy, specificity, sensitivity, and Mathew's correlation coefficient. Our model achieved an accuracy of 81.5%, representing a 2-3% improvement over existing models.

Structure of phosphorylated-like RssB, the adaptor delivering σs to the ClpXP proteolytic machinery, reveals an interface switch for activation

In enterobacteria such as Escherichia coli, the general stress response is mediated by σs, the stationary phase dissociable promoter specificity subunit of RNA polymerase. σs is degraded by ClpXP during active growth in a process dependent on the RssB adaptor, which is thought to be stimulated by the phosphorylation of a conserved aspartate in its N-terminal receiver domain. Here we present the crystal structure of full-length RssB bound to a beryllofluoride phosphomimic. Compared to the structure of RssB bound to the IraD anti-adaptor, our new RssB structure with bound beryllofluoride reveals conformational differences and coil-to-helix transitions in the C-terminal region of the RssB receiver domain and in the interdomain segmented helical linker. These are accompanied by masking of the α4-β5-α5 (4-5-5) "signaling" face of the RssB receiver domain by its C-terminal domain. Critically, using hydrogen-deuterium exchange mass spectrometry, we identify σs-binding determinants on the 4-5-5 face, implying that this surface needs to be unmasked to effect an interdomain interface switch and enable full σs engagement and hand-off to ClpXP. In activated receiver domains, the 4-5-5 face is often the locus of intermolecular interactions, but its masking by intramolecular contacts upon phosphorylation is unusual, emphasizing that RssB is a response regulator that undergoes atypical regulation.

Identification of novel natural product inhibitors of BRD4 using high throughput virtual screening and MD simulation

Bromodomains are evolutionarily conserved structural motifs that recognize acetylated lysine residues on histone tails. They play a crucial role in shaping chromatin architecture and regulating gene expression in various biological processes. Mutations in bromodomains containing proteins lead to multiple human diseases, which makes them attractive target for therapeutic intervention. Extensive studies have been done on BRD4 as a target for several cancers, such as Acute Myeloid Leukemia (AML) and Burkitt Lymphoma. Several potential inhibitors have been identified against the BRD4 bromodomain. However, most of these inhibitors have drawbacks such as non-specificity and toxicity, decreasing their appeal and necessitating the search for novel non-toxic inhibitors. This study aims to address this need by virtually screening natural compounds from the NPASS database against the Kac binding site of BRD4-BD1 using high throughput molecular docking followed by similarity clustering, pharmacokinetic screening, MD simulation and MM-PBSA binding free energy calculations. Using this approach, we identified five natural product inhibitors having a similar or better binding affinity to the BRD4 bromodomain compared to JQ1 (previously reported inhibitor of BRD4). Further systematic analysis of these inhibitors resulted in the top three hits: NPC268484 (Palodesangren-B), NPC295021 (Candidine) and NPC313112 (Buxifoliadine-D). Collectively, our in silico results identified some promising natural products that have the potential to act as potent BRD4-BD1 inhibitors and can be considered for further validation through future in vitro and in vivo studies.Communicated by Ramaswamy H. Sarma.

An Encodable Scaffold for Sequence-Specific Recognition of Duplex RNA

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

Short tandem repeats bind transcription factors to tune eukaryotic gene expression

Short tandem repeats (STRs) are enriched in eukaryotic cis-regulatory elements and alter gene expression, yet how they regulate transcription remains unknown. We found that STRs modulate transcription factor (TF)-DNA affinities and apparent on-rates by about 70-fold by directly binding TF DNA-binding domains, with energetic impacts exceeding many consensus motif mutations. STRs maximize the number of weakly preferred microstates near target sites, thereby increasing TF density, with impacts well predicted by statistical mechanics. Confirming that STRs also affect TF binding in cells, neural networks trained only on in vivo occupancies predicted effects identical to those observed in vitro. Approximately 90% of TFs preferentially bound STRs that need not resemble known motifs, providing a cis-regulatory mechanism to target TFs to genomic sites.

Structural characterization of transcription-coupled repair protein UVSSA and its interaction with TFIIH protein

UV-stimulated scaffold protein A (UVSSA) is a key protein in the Transcription-Coupled Nucleotide Excision Repair (TC-NER) pathway. UVSSA, an intrinsically disordered protein, interacts with multiple members of the pathway, tethering them into the complex. Several studies have reported that UVSSA recruits Transcription Factor IIH (TFIIH) via direct interaction, following which CSB is degraded and the lesion recognition TC-NER complex dissociates from the damage site to facilitate the DNA repair. Structural insights into these events remain largely unknown. Herein, we have investigated the interaction of human UVSSA with the Pleckstrin-Homology-domain of p62 subunit of TFIIH (p62-PHD) using biophysical techniques. We observed that UVSSA forms a stable complex with the p62-PHD in vitro. Small-angle scattering measurements using X-rays and neutrons revealed a significant change in pair-distance distribution function for UVSSA662/p62-PHD complex compared to UVSSA alone. Additionally, a significant decrease was observed in the radius of gyration of the complex. Our findings suggest that TFIIH binding to UVSSA causes significant conformational changes in UVSSA. We hypothesize that these conformational changes play an important role in the dissociation of the lesion recognition TC-NER complex.

Avoidance of the use of tryptophan in buried chromosomal proteins as a mechanism for reducing photo/oxidative damage to genomes

In cells that are exposed to terrestrial sunlight, the indole moiety in the side chain of tryptophan (Trp) can suffer photo/oxidative damage (POD) by reactive oxygen species (ROS) and/or ultraviolet light (UV-B). Trp is oxidized to produce N-formylkynurenine (NFK), a UV-A-responsive photosensitizer that further degenerates into photosensitizers capable of generating ROS through exposure to visible light. Thus, Trp-containing proteins function as both victims, and perpetrators, of POD if they are not rapidly replaced through protein turnover. The literature indicates that protein turnover and DNA repair occur poorly in chromosomal interiors. We contend, therefore, that basic chromosomal proteins (BCPs) that are enveloped by DNA should have evolved to lack Trp residues in their amino acid sequences, since these could otherwise function as 'Trojan horse-type' DNA-damaging agents. Our global analyses of protein sequences demonstrates that BCPs consistently lack Trp residues, although DNA-binding proteins in general do not display such a lack. We employ HU-B (a wild-type, Trp-lacking bacterial BCP) and HU-B F47W (a mutant, Trp-containing form of the same bacterial BCP) to demonstrate that the possession of Trp is deleterious to BCPs and associated chromosomal DNA. Basically, we show that UV-B and UV-A (a) cause no POD in HU-B, but cause extensive POD in HU-B F47W (in vitro), as well as (b) only nominal DNA damage in bacteria expressing HU-B, but extensive DNA damage in bacteria expressing F47W HU-B (in vivo). Our results suggest that Trp-lacking BCPs could have evolved to reduce scope for protein-facilitated, sunlight-mediated damage of DNA by UV-A and visible light, within chromosomal interiors that are poorly serviced by protein turnover and DNA repair machinery.

On the dependent recognition of some long zinc finger proteins

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

Lysine Methyltransferase SMYD1 Regulates Myogenesis via skNAC Methylation

The SMYD family is a unique class of lysine methyltransferases (KMTases) whose catalytic SET domain is split by a MYND domain. Among these, Smyd1 was identified as a heart- and skeletal muscle-specific KMTase and is essential for cardiogenesis and skeletal muscle development. SMYD1 has been characterized as a histone methyltransferase (HMTase). Here we demonstrated that SMYD1 methylates is the Skeletal muscle-specific splice variant of the Nascent polypeptide-Associated Complex (skNAC) transcription factor. SMYD1-mediated methylation of skNAC targets K1975 within the carboxy-terminus region of skNAC. Catalysis requires physical interaction of SMYD1 and skNAC via the conserved MYND domain of SMYD1 and the PXLXP motif of skNAC. Our data indicated that skNAC methylation is required for the direct transcriptional activation of myoglobin (Mb), a heart- and skeletal muscle-specific hemoprotein that facilitates oxygen transport. Our study revealed that the skNAC, as a methylation target of SMYD1, illuminates the molecular mechanism by which SMYD1 cooperates with skNAC to regulate transcriptional activation of genes crucial for muscle functions and implicates the MYND domain of the SMYD-family KMTases as an adaptor to target substrates for methylation.

The SAGA HAT module is tethered by its SWIRM domain and modulates activity of the SAGA DUB module

The SAGA (Spt-Ada-Gcn5 acetyltransferase) complex is a transcriptional co-activator that both acetylates and deubiquitinates histones. The histone acetyltransferase (HAT) subunit, Gcn5, is part of a subcomplex of SAGA called the HAT module. A minimal HAT module complex containing Gcn5 bound to Ada2 and Ada3 is required for full Gcn5 activity on nucleosomes. Deletion studies have suggested that the Ada2 SWIRM domain plays a role in tethering the HAT module to the remainder of SAGA. While recent cryo-EM studies have resolved the structure of the core of the SAGA complex, the HAT module subunits and molecular details of its interactions with the SAGA core could not be resolved. Here we show that the SWIRM domain is required for incorporation of the HAT module into the yeast SAGA complex, but not the ADA complex, a distinct six-protein acetyltransferase complex that includes the SAGA HAT module proteins. In the isolated Gcn5/Ada2/Ada3 HAT module, deletion of the SWIRM domain modestly increased activity but had negligible effect on nucleosome binding. Loss of the HAT module due to deletion of the SWIRM domain decreases the H2B deubiquitinating activity of SAGA, indicating a role for the HAT module in regulating SAGA DUB module activity. A model of the HAT module created with Alphafold Multimer provides insights into the structural basis for our biochemical data, as well as prior deletion studies.

Structural basis of DNA binding by the WhiB-like transcription factor WhiB3 in Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) WhiB3 is an iron-sulfur cluster-containing transcription factor belonging to a subclass of the WhiB-Like (Wbl) family that is widely distributed in the phylum Actinobacteria. WhiB3 plays a crucial role in the survival and pathogenesis of Mtb. It binds to the conserved region 4 of the principal sigma factor (σA4) in the RNA polymerase holoenzyme to regulate gene expression like other known Wbl proteins in Mtb. However, the structural basis of how WhiB3 coordinates with σA4 to bind DNA and regulate transcription is unclear. Here we determined crystal structures of the WhiB3:σA4 complex without and with DNA at 1.5 Å and 2.45 Å, respectively, to elucidate how WhiB3 interacts with DNA to regulate gene expression. These structures reveal that the WhiB3:σA4 complex shares a molecular interface similar to other structurally characterized Wbl proteins and also possesses a subclass-specific Arg-rich DNA-binding motif. We demonstrate that this newly defined Arg-rich motif is required for WhiB3 binding to DNA in vitro and transcriptional regulation in Mycobacterium smegmatis. Together, our study provides empirical evidence of how WhiB3 regulates gene expression in Mtb by partnering with σA4 and engaging with DNA via the subclass-specific structural motif, distinct from the modes of DNA interaction by WhiB1 and WhiB7.

Structural and functional basis of the universal transcription factor NusG pro-pausing activity in Mycobacterium tuberculosis

Transcriptional pauses mediate regulation of RNA biogenesis. DNA-encoded pause signals trigger pausing by stabilizing RNA polymerase (RNAP) swiveling and inhibiting DNA translocation. The N-terminal domain (NGN) of the only universal transcription factor, NusG/Spt5, modulates pausing through contacts to RNAP and DNA. Pro-pausing NusGs enhance pauses, whereas anti-pausing NusGs suppress pauses. Little is known about pausing and NusG in the human pathogen Mycobacterium tuberculosis (Mtb). We report that MtbNusG is pro-pausing. MtbNusG captures paused, swiveled RNAP by contacts to the RNAP protrusion and nontemplate-DNA wedged between the NGN and RNAP gate loop. In contrast, anti-pausing Escherichia coli (Eco) NGN contacts the MtbRNAP gate loop, inhibiting swiveling and pausing. Using CRISPR-mediated genetics, we show that pro-pausing NGN is required for mycobacterial fitness. Our results define an essential function of mycobacterial NusG and the structural basis of pro- versus anti-pausing NusG activity, with broad implications for the function of all NusG orthologs.

Large-scale mapping and mutagenesis of human transcriptional effector domains

Human gene expression is regulated by more than 2,000 transcription factors and chromatin regulators1,2. Effector domains within these proteins can activate or repress transcription. However, for many of these regulators we do not know what type of effector domains they contain, their location in the protein, their activation and repression strengths, and the sequences that are necessary for their functions. Here, we systematically measure the effector activity of more than 100,000 protein fragments tiling across most chromatin regulators and transcription factors in human cells (2,047 proteins). By testing the effect they have when recruited at reporter genes, we annotate 374 activation domains and 715 repression domains, roughly 80% of which are new and have not been previously annotated3-5. Rational mutagenesis and deletion scans across all the effector domains reveal aromatic and/or leucine residues interspersed with acidic, proline, serine and/or glutamine residues are necessary for activation domain activity. Furthermore, most repression domain sequences contain sites for small ubiquitin-like modifier (SUMO)ylation, short interaction motifs for recruiting corepressors or are structured binding domains for recruiting other repressive proteins. We discover bifunctional domains that can both activate and repress, some of which dynamically split a cell population into high- and low-expression subpopulations. Our systematic annotation and characterization of effector domains provide a rich resource for understanding the function of human transcription factors and chromatin regulators, engineering compact tools for controlling gene expression and refining predictive models of effector domain function.

Short Activators and Repressors of RNA Toehold Switches

RNA toehold switches are a widely used class of molecule to detect specific RNA "trigger" sequences, but their design, intended function, and characterization to date leave it unclear whether they can function properly with triggers shorter than 36 nucleotides. Here, we explore the feasibility of using standard toehold switches with 23-nucleotide truncated triggers. We assess the crosstalk of different triggers with significant homology and identify a highly sensitive trigger region where just one mutation from the consensus trigger sequence can reduce switch activation by 98.6%. However, we also find that triggers with as many as seven mutations outside of this region can still lead to 5-fold induction of the switch. We also present a new approach using 18- to 22-nucleotide triggers as translational repressors for toehold switches and assess the off-target regulation for this strategy as well. The development and characterization of these strategies could help enable applications like microRNA sensors, where well-characterized crosstalk between sensors and detection of short target sequences are critical.

SARS-CoV-2 envelope protein attain Kac mediated dynamical interaction network to adopt 'histone mimic' at BRD4 interface

Interface mimicry, achieved by recognition of host-pathogen interactions, is the basis by which pathogen proteins can hijack the host machinery. The envelope (E) protein of SARS-CoV-2 is reported to mimic the histones at the BRD4 surface via establishing the structural mimicry; however, the underlying mechanism of E protein mimicking the histones is still elusive. To explore the mimics at dynamic and structural residual network level an extensive docking, and MD simulations were carried out in a comparative manner between complexes of H3-, H4-, E-, and apo-BRD4. We identified that E peptide is able to attain an 'interaction network mimicry', as its acetylated lysine (Kac) achieves orientation and residual fingerprint similar to histones, including water-mediated interactions for both the Kac positions. We identified Y59 of E, playing an anchor role to escort lysine positioning inside the binding site. Furthermore, the binding site analysis confirms that E peptide needs a higher volume, similar to the H4-BRD4 where both the lysine's (Kac5 and Kac8) can accommodate nicely, however, the position of Kac8 is mimicked by two additional water molecules other than four water-mediated bridging's, strengthening the possibility that E peptide could hijack host BRD4 surface. These molecular insights seem pivotal for mechanistic understanding and BRD4-specific therapeutic intervention. KEY POINTSMolecular mimicry is reported in hijacking and then outcompeting the host counterparts so that pathogens can rewire their cellular function by overcoming the host defense mechanism.The molecular recognition process is the basis of molecular mimicry. The E peptide of SARS-CoV-2 is reported to mimic host histone at the BRD4 surface by utilizing its C-terminally placed acetylated lysine (Kac63) to mimic the N-terminally placed acetylated lysine Kac5GGKac8 histone (H4) by interaction network mimicry identified through microsecond molecular dynamics (MD) simulations and post-processing extensive analysis.There are two steps to mimic: firstly, tyrosine residues help E to anchor at the BRD4 surface to position Kac and increase the volume of the pocket. Secondary, after positioning of Kac, a common durable interaction network N140:Kac5; Kac5:W1; W1:Y97; W1:W2; W2:W3; W3:W4; W4:P82 is established between Kac5, with key residues P82, Y97, N140, and four water molecules through water mediate bridge. Furthermore, the second acetylated lysine Kac8 position and its interaction as polar contact with Kac5 were also mimicked by E peptide through interaction network P82:W5; W5:Kac63; W5:W6; W6:Kac63.The binding event at BRD4/BD1 seems an induced-fit mechanism as a bigger binding site volume was identified at H4-BRD4 on which E peptide attains its better stability than H3-BRD4.We identified the tyrosine residue Y59 of E that acts like an anchor on the BRD4 surface to position Kac inside the pocket and attain the interaction network by using aromatic residues of the BRD4 surface.Communicated by Ramaswamy H. Sarma.

Allosteric autoregulation of DNA binding via a DNA-mimicking protein domain: a biophysical study of ZNF410-DNA interaction using small angle X-ray scattering

ZNF410 is a highly-conserved transcription factor, remarkable in that it recognizes a 15-base pair DNA element but has just a single responsive target gene in mammalian erythroid cells. ZNF410 includes a tandem array of five zinc-fingers (ZFs), surrounded by uncharacterized N- and C-terminal regions. Unexpectedly, full-length ZNF410 has reduced DNA binding affinity, compared to that of the isolated DNA binding ZF array, both in vitro and in cells. AlphaFold predicts a partially-folded N-terminal subdomain that includes a 30-residue long helix, preceded by a hairpin loop rich in acidic (aspartate/glutamate) and serine/threonine residues. This hairpin loop is predicted by AlphaFold to lie against the DNA binding interface of the ZF array. In solution, ZNF410 is a monomer and binds to DNA with 1:1 stoichiometry. Surprisingly, the single best-fit model for the experimental small angle X-ray scattering profile, in the absence of DNA, is the original AlphaFold model with the N-terminal long-helix and the hairpin loop occupying the ZF DNA binding surface. For DNA binding, the hairpin loop presumably must be displaced. After combining biophysical, biochemical, bioinformatic and artificial intelligence-based AlphaFold analyses, we suggest that the hairpin loop mimics the structure and electrostatics of DNA, and provides an additional mechanism, supplementary to sequence specificity, of regulating ZNF410 DNA binding.

Investigation of novel ligand targeting bromodomain-containing protein 4 (BRD4) for cancer drug discovery: complete pharmacophore approach

The Bromodomain (BRD4) and extra-terminal (BET) protein family are reversible; lysine-acetylated epigenetic readers identified as key important epigenetic regulators for protein recognition in posttranslational modifications for targeting cancer for its role in super-enhancers and transcription of oncogene expression in cancer and other forms of cancer and various diseases. Firstly, JQ-1a small potent BET inhibitors, targeting BET proteins were currently in clinical trials to ablate cancer. The identified compounds were taken from the library of preexisting therapeutically potent molecules. The objective of the present study is to identify the potential small molecule inhibitors against BRD4 through in-silico approach for the treatment of cancer. In present study, designed an in-silico screening of small molecules through ligand-based pharmacophore studies against bromodomain-containing protein 4 (BRD-4) protein and used for virtual screening through Database and their binding affinity and interaction of identified molecules were predicted through molecular docking, molecular dynamics simulations for 12 fixed time period, Molecular mechanics (MMGBSA) binding free energy calculations, ADME with drug-likeness properties including violations of lipinski's rule of 5, Jorgensens rule of 3 and other parameters were studied. The docking results indicate from the reported database screened molecules were validated with docking score -7.92 to -4.27Kcal/mol for BRD4-BD1 and the best model identified 21 hits. Among these two drugs were filtered and scrutinized for their ability based on binding modes and common interaction, MMGBSA of the highest affinity -54.53 Kcal/mol of BRD4-BD1 and ADME properties of selected molecules were predicted for its various parameters, dynamics studies evaluating its binding stability using Maestro software. In Conclusion, two BRD4 inhibitors were found to bind strongly in the similar binding sites as JQ-1, highlighting the role of BRD4-BD1. These compounds were identified as promising new options for regulating epigenetics and understanding the structural needs of BRD4 protein, further research in these areas could lead to the development of more effective and targeted cancer drugs.Communicated by Ramaswamy H. Sarma.

A competitive precision CRISPR method to identify the fitness effects of transcription factor binding sites

Here we describe a competitive genome editing method that measures the effect of mutations on molecular functions, based on precision CRISPR editing using template libraries with either the original or altered sequence, and a sequence tag, enabling direct comparison between original and mutated cells. Using the example of the MYC oncogene, we identify important transcriptional targets and show that E-box mutations at MYC target gene promoters reduce cellular fitness.

Cytosporone E analogues as BRD4 inhibitors for cancer treatment: molecular docking and molecular dynamic investigations

Cancer is considered one of the worldwide life-threatening and leading causes of human mortality. In 2020, 19,292,789 cancer cases and 9,958,133 cancer deaths have been estimated worldwide. Therefore, efforts have been devoted to discover novel anticancer agents. Bromodomains have a vital role in the regulation of transcription. Many reports have shown that bromodomain-containing protein 4 (BRD4) is an important target for cancer therapeutics. In this study, several in silico approaches were utilized to discover new inhibitors against the BRD4 protein using the Schrodinger suite. A library of 27 cytosporone E derivatives was docked into the active site of the BRD4 protein. Docked ligands showed docking scores ranging between -11.289 to -3.992 Kcal/mol. Ligands 1-4 showed better binding affinities with docking scores ranging from -11.289 to -8.917 Kcal/mol compared to the reference ligand BI-2536 (-8.426 Kcal/mol). These ligands displayed favorable MM-GBSA free binding energy. Also, ligands 1-4 were subjected to molecular dynamics simulations for 100 ns to get insight into the ligand-binding stability. These compounds exhibited an average RMSD below 2.8 Å, indicating the stability of the compounds with BRD4 protein. Further, Moreover, ligands 1-3 displayed favorable AMDET properties (absorption, distribution, metabolism, excretion, and toxicity). These new compounds might be potential leads to combat cancer.Communicated by Ramaswamy H. Sarma.

N-terminal β-strand in YAP is critical for stronger binding to scalloped relative to TEAD transcription factor

The yes-associated protein (YAP) regulates the transcriptional activity of the TEAD transcription factors that are key in the control of organ morphogenesis. YAP interacts with TEAD via three secondary structure elements: a β-strand, an α-helix, and an Ω-loop. Earlier results have shown that the β-strand has only a marginal contribution in the YAP:TEAD interaction, but we show here that it significantly enhances the affinity of YAP for the Drosophila homolog of TEAD, scalloped (Sd). Nuclear magnetic resonance shows that the β-strand adopts a more rigid conformation once bound to Sd; pre-steady state kinetic measurements show that the YAP:Sd complex is more stable. Although the crystal structures of the YAP:TEAD and YAP:Sd complexes reveal no differences at the binding interface that could explain these results. Molecular Dynamics simulations are in line with our experimental findings regarding β-strand stability and overall binding affinity of YAP to TEAD and Sd. In particular, RMSF, correlated motion and MMGBSA analyses suggest that β-sheet fluctuations play a relevant role in YAP53-57 β-strand dissociation from TEAD4 and contribute to the lower affinity of YAP for TEAD4. Identifying a clear mechanism leading to the difference in YAP's β-strand stability proved to be challenging, pointing to the potential relevance of multiple modest structural changes or fluctuations for regulation of binding affinity.

Genome-wide analysis and molecular dissection of the SPL gene family in Fraxinus mandshurica

BACKGROUND

SQUAMOSA promoter binding protein-like (SPL) is a unique family of transcription factors in plants, which is engaged in regulating plant growth and development, physiological and biochemical processes. Fraxinus mandshurica is an excellent timber species with a wide range of uses in northeastern China and enjoys a high reputation in the international market. SPL family analysis has been reported in some plants while SPL family analysis of Fraxinus mandshurica has not been reported.

RESULTS

We used phylogeny, conserved motifs, gene structure, secondary structure prediction, miR156 binding sites, promoter cis elements and GO annotation to systematically analyze the FmSPLs family. This was followed by expression analysis by subcellular localization, expression patterns at various tissue sites, abiotic stress and hormone induction. Because FmSPL2 is highly expressed in flowers it was selected to describe the SPL gene family of Fraxinus mandshurica by ectopic expression. Among them, 10 FmSPL genes that were highly expressed at different loci were selected for expression analysis under abiotic stress (NaCl and Cold) and hormone induction (IAA and ABA). These 10 FmSPL genes showed corresponding trends in response to both abiotic stress and hormone induction. We showed that overexpression of FmSPL2 in transgenic Nicotiana tabacum L. resulted in taller plants, shorter root length, increased root number, rounded leaves, and earlier flowering time.

CONCLUSIONS

We identified 36 SPL genes, which were classified into seven subfamilies based on sequence analysis. FmSPL2 was selected for subsequent heterologous expression by analysis of expression patterns in various tissues and under abiotic stress and hormone induction, and significant phenotypic changes were observed in the transgenic Nicotiana tabacum L. These results provide insight into the evolutionary origin and biological significance of plant SPL. The aim of this study was to lay the foundation for the genetic improvement of Fraxinus mandshurica and the subsequent functional analysis of FmSPL2.

Genome-Wide Identification and Analysis of the R2R3-MYB Gene Family in Theobroma cacao

The MYB gene family is involved in the regulation of plant growth, development and stress responses. In this paper, to identify Theobroma cacao R2R3-MYB (TcMYB) genes involved in environmental stress and phytohormones, we conducted a genome-wide analysis of the R2R3-MYB gene family in Theobroma cacao (cacao). A total of 116 TcMYB genes were identified, and they were divided into 23 subgroups according to the phylogenetic analysis. Meanwhile, the conserved motifs, gene structures and cis-acting elements of promoters were analyzed. Moreover, these TcMYB genes were distributed on 10 chromosomes. We conducted a synteny analysis to understand the evolution of the cacao R2R3-MYB gene family. A total of 37 gene pairs of TcMYB genes were identified through tandem or segmental duplication events. Additionally, we also predicted the subcellular localization and physicochemical properties. All the studies showed that TcMYB genes have multiple functions, including responding to environmental stresses. The results provide an understanding of R2R3-MYB in Theobroma cacao and lay the foundation for a further functional analysis of TcMYB genes in the growth of cacao.

Structure of a nucleosome-bound MuvB transcription factor complex reveals DNA remodelling

Genes encoding the core cell cycle machinery are transcriptionally regulated by the MuvB family of protein complexes in a cell cycle-specific manner. Complexes of MuvB with the transcription factors B-MYB and FOXM1 activate mitotic genes during cell proliferation. The mechanisms of transcriptional regulation by these complexes are still poorly characterised. Here, we combine biochemical analysis and in vitro reconstitution, with structural analysis by cryo-electron microscopy and cross-linking mass spectrometry, to functionally examine these complexes. We find that the MuvB:B-MYB complex binds and remodels nucleosomes, thereby exposing nucleosomal DNA. This remodelling activity is supported by B-MYB which directly binds the remodelled DNA. Given the remodelling activity on the nucleosome, we propose that the MuvB:B-MYB complex functions as a pioneer transcription factor complex. In this work, we rationalise prior biochemical and cellular studies and provide a molecular framework of interactions on a protein complex that is key for cell cycle regulation.

The structure and activities of the archaeal transcription termination factor Eta detail vulnerabilities of the transcription elongation complex

Transcription must be properly regulated to ensure dynamic gene expression underlying growth, development, and response to environmental cues. Regulation is imposed throughout the transcription cycle, and while many efforts have detailed the regulation of transcription initiation and early elongation, the termination phase of transcription also plays critical roles in regulating gene expression. Transcription termination can be driven by only a few proteins in each domain of life. Detailing the mechanism(s) employed provides insight into the vulnerabilities of transcription elongation complexes (TECs) that permit regulated termination to control expression of many genes and operons. Here, we describe the biochemical activities and crystal structure of the superfamily 2 helicase Eta, one of two known factors capable of disrupting archaeal transcription elongation complexes. Eta retains a twin-translocase core domain common to all superfamily 2 helicases and a well-conserved C terminus wherein individual amino acid substitutions can critically abrogate termination activities. Eta variants that perturb ATPase, helicase, single-stranded DNA and double-stranded DNA translocase and termination activities identify key regions of the C terminus of Eta that, when combined with modeling Eta-TEC interactions, provide a structural model of Eta-mediated termination guided in part by structures of Mfd and the bacterial TEC. The susceptibility of TECs to disruption by termination factors that target the upstream surface of RNA polymerase and potentially drive termination through forward translocation and allosteric mechanisms that favor opening of the clamp to release the encapsulated nucleic acids emerges as a common feature of transcription termination mechanisms.

ANANASTRA: annotation and enrichment analysis of allele-specific transcription factor binding at SNPs

We present ANANASTRA, https://ananastra.autosome.org, a web server for the identification and annotation of regulatory single-nucleotide polymorphisms (SNPs) with allele-specific binding events. ANANASTRA accepts a list of dbSNP IDs or a VCF file and reports allele-specific binding (ASB) sites of particular transcription factors or in specific cell types, highlighting those with ASBs significantly enriched at SNPs in the query list. ANANASTRA is built on top of a systematic analysis of allelic imbalance in ChIP-Seq experiments and performs the ASB enrichment test against background sets of SNPs found in the same source experiments as ASB sites but not displaying significant allelic imbalance. We illustrate ANANASTRA usage with selected case studies and expect that ANANASTRA will help to conduct the follow-up of GWAS in terms of establishing functional hypotheses and designing experimental verification.

Coevolution-derived native and non-native contacts determine the emergence of a novel fold in a universally conserved family of transcription factors

The NusG protein family is structurally and functionally conserved in all domains of life. Its members directly bind RNA polymerases and regulate transcription processivity and termination. RfaH, a divergent sub-family in its evolutionary history, is known for displaying distinct features than those in NusG proteins, which allows them to regulate the expression of virulence factors in enterobacteria in a DNA sequence-dependent manner. A striking feature is its structural interconversion between an active fold, which is the canonical NusG three-dimensional structure, and an autoinhibited fold, which is distinctively novel. How this novel fold is encoded within RfaH sequence to encode a metamorphic protein remains elusive. In this work, we used publicly available genomic RfaH protein sequences to construct a complete multiple sequence alignment, which was further augmented with metagenomic sequences and curated by predicting their secondary structure propensities using JPred. Coevolving pairs of residues were calculated from these sequences using plmDCA and GREMLIN, which allowed us to detect the enrichment of key metamorphic contacts after sequence filtering. Finally, we combined our coevolutionary predictions with molecular dynamics to demonstrate that these interactions are sufficient to predict the structures of both native folds, where coevolutionary-derived non-native contacts may play a key role in achieving the compact RfaH novel fold. All in all, emergent coevolutionary signals found within RfaH sequences encode the autoinhibited and active folds of this protein, shedding light on the key interactions responsible for the action of this metamorphic protein.

Mechanism of mitoribosomal small subunit biogenesis and preinitiation

Mitoribosomes are essential for the synthesis and maintenance of bioenergetic proteins. Here we use cryo-electron microscopy to determine a series of the small mitoribosomal subunit (SSU) intermediates in complex with auxiliary factors, revealing a sequential assembly mechanism. The methyltransferase TFB1M binds to partially unfolded rRNA h45 that is promoted by RBFA, while the mRNA channel is blocked. This enables binding of METTL15 that promotes further rRNA maturation and a large conformational change of RBFA. The new conformation allows initiation factor mtIF3 to already occupy the subunit interface during the assembly. Finally, the mitochondria-specific ribosomal protein mS37 (ref. 1) outcompetes RBFA to complete the assembly with the SSU-mS37-mtIF3 complex2 that proceeds towards mtIF2 binding and translation initiation. Our results explain how the action of step-specific factors modulate the dynamic assembly of the SSU, and adaptation of a unique protein, mS37, links the assembly to initiation to establish the catalytic human mitoribosome.

Biolayer interferometry for DNA-protein interactions

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

Genome-wide identification of Tomato Golden 2-Like transcription factors and abiotic stress related members screening

BACKGROUND

Golden 2-Like (G2-like) transcription factors play an important role in plant development. However, the roles of these G2-like regulatory genes in response to abiotic stresses in tomato are not well understood.

RESULTS

In this study, we identified 66 putative G2-like genes in tomato (Solanum lycopersicum) and classified them into 5 groups (I to V) according to gene structure, motif composition and phylogenetic analysis. The G2-like genes were unevenly distributed across all 12 chromosomes. There were nine pairs of duplicated gene segments and four tandem duplicated SlGlk genes. Analysis of the cis-regulatory elements (CREs) showed that the promoter regions of SlGlks contain many kinds of stress- and hormone-related CREs. Based on RNA-seq, SlGlks were expressed in response to three abiotic stresses. Thirty-six differentially expressed SlGlks were identified; these genes have multiple functions according to Gene Ontology (GO) analysis and are enriched mainly in the zeatin biosynthesis pathway. Further studies exhibited that silencing SlGlk16 in tomato would reduce drought stress tolerance by earlier wilted, lower superoxide dismutase (SOD), peroxidase (POD) activities, less Pro contents and more MDA contents.

CONCLUSIONS

Overall, the results of this study provide comprehensive information on G2-like transcription factors and G2-like genes that may be expressed in response to abiotic stresses.

Novel Glial Cells Missing-2 (GCM2) variants in parathyroid disorders

OBJECTIVE

The aim of this study was to analyze variants of the gene glial cells missing-2 (GCM2), encoding a parathyroid cell-specific transcription factor, in familial hypoparathyroidism and in familial isolated hyperparathyroidism (FIHP) without and with parathyroid carcinoma.

DESIGN

We characterized 2 families with hypoparathyroidism and 19 with FIHP in which we examined the mechanism of action of GCM2 variants.

METHODS

Leukocyte DNA of hypoparathyroid individuals was Sanger sequenced for CASR, PTH, GNA11 and GCM2 mutations. DNA of hyperparathyroid individuals underwent MEN1, CDKN1B, CDC73, CASR, RET and GCM2 sequencing. The actions of identified GCM2 variants were evaluated by in vitro functional analyses.

RESULTS

A novel homozygous p.R67C GCM2 mutation which failed to stimulate transcriptional activity in a luciferase assay was identified in affected members of two hypoparathyroid families. Oligonucleotide pull-down assay and in silico structural modeling indicated that this mutant had lost the ability to bind the consensus GCM recognition sequence of DNA. Two novel (p.I383M and p.T386S) and one previously reported (p.Y394S) heterozygous GCM2 variants that lie within a C-terminal conserved inhibitory domain were identified in three affected individuals of the hyperparathyroid families. One family member, heterozygous for p.I138M, had parathyroid carcinoma (PC), and a heterozygous p.V382M variant was found in another patient affected by sporadic PC. These variants exerted significantly enhanced in vitrotranscriptional activity, including increased stimulation of the PTH promoter.

CONCLUSIONS

We provide evidence that two novel GCM2 R67C inactivating mutations with an inability to bind DNA are causative of hypoparathyroidism. Additionally, we provide evidence that two novel GCM2 variants increased transactivation of the PTH promoter in vitro and are associated with FIHP. Furthermore, our studies suggest that activating GCM2 variants may contribute to facilitating more aggressive parathyroid disease.

Novel Glial Cells Missing-2 (GCM2) variants in parathyroid disorders

OBJECTIVE

The aim of this study was to analyze variants of the gene glial cells missing-2 (GCM2), encoding a parathyroid cell-specific transcription factor, in familial hypoparathyroidism and in familial isolated hyperparathyroidism (FIHP) without and with parathyroid carcinoma.

DESIGN

We characterized 2 families with hypoparathyroidism and 19 with FIHP in which we examined the mechanism of action of GCM2 variants.

METHODS

Leukocyte DNA of hypoparathyroid individuals was Sanger sequenced for CASR, PTH, GNA11 and GCM2 mutations. DNA of hyperparathyroid individuals underwent MEN1, CDKN1B, CDC73, CASR, RET and GCM2 sequencing. The actions of identified GCM2 variants were evaluated by in vitro functional analyses.

RESULTS

A novel homozygous p.R67C GCM2 mutation which failed to stimulate transcriptional activity in a luciferase assay was identified in affected members of two hypoparathyroid families. Oligonucleotide pull-down assay and in silico structural modeling indicated that this mutant had lost the ability to bind the consensus GCM recognition sequence of DNA. Two novel (p.I383M and p.T386S) and one previously reported (p.Y394S) heterozygous GCM2 variants that lie within a C-terminal conserved inhibitory domain were identified in three affected individuals of the hyperparathyroid families. One family member, heterozygous for p.I138M, had parathyroid carcinoma (PC), and a heterozygous p.V382M variant was found in another patient affected by sporadic PC. These variants exerted significantly enhanced in vitrotranscriptional activity, including increased stimulation of the PTH promoter.

CONCLUSIONS

We provide evidence that two novel GCM2 R67C inactivating mutations with an inability to bind DNA are causative of hypoparathyroidism. Additionally, we provide evidence that two novel GCM2 variants increased transactivation of the PTH promoter in vitro and are associated with FIHP. Furthermore, our studies suggest that activating GCM2 variants may contribute to facilitating more aggressive parathyroid disease.

Structural basis of cyclic oligoadenylate binding to the transcription factor Csa3 outlines cross talk between type III and type I CRISPR systems

RNA interference by type III CRISPR systems results in the synthesis of cyclic oligoadenylate (cOA) second messengers, which are known to bind and regulate various CARF domain-containing nuclease receptors. The CARF domain-containing Csa3 family of transcriptional factors associated with the DNA-targeting type I CRISPR systems regulate expression of various CRISPR and DNA repair genes in many prokaryotes. In this study, we extend the known receptor repertoire of cOA messengers to include transcriptional factors by demonstrating specific binding of cyclic tetra-adenylate (cA4) to Saccharolobus solfataricus Csa3 (Csa3Sso). Our 2.0-Å resolution X-ray crystal structure of cA4-bound full-length Csa3Sso reveals the binding of its CARF domain to an elongated conformation of cA4. Using cA4 binding affinity analyses of Csa3Sso mutants targeting the observed Csa3Sso•cA4 structural interface, we identified a Csa3-specific cA4 binding motif distinct from a more widely conserved cOA-binding CARF motif. Using a rational surface engineering approach, we increased the cA4 binding affinity of Csa3Sso up to ∼145-fold over the wildtype, which has potential applications for future second messenger-driven CRISPR gene expression and editing systems. Our in-solution Csa3Sso structural analysis identified cA4-induced allosteric and asymmetric conformational rearrangement of its C-terminal winged helix-turn-helix effector domains, which could potentially be incompatible to DNA binding. However, specific in vitro binding of the purified Csa3Sso to its putative promoter (PCas4a) was found to be cA4 independent, suggesting a complex mode of Csa3Sso regulation. Overall, our results support cA4-and Csa3-mediated cross talk between type III and type I CRISPR systems.

The Intrinsically Disordered Proteins MLLT3 (AF9) and MLLT1 (ENL) - Multimodal Transcriptional Switches With Roles in Normal Hematopoiesis, MLL Fusion Leukemia, and Kidney Cancer

AF9 (MLLT3) and ENL (MLLT1) are members of the YEATS family (named after the five proteins first shown to contain this domain: Yaf9, ENL, AF9, Taf14, Sas5) defined by the presence of a YEATS domain. The YEATS domain is an epigenetic reader that binds to acetylated and crotonylated lysines, unlike the bromodomain which can only bind to acetylated lysines. All members of this family have been shown to be components of various complexes with roles in chromatin remodeling, histone modification, histone variant deposition, and transcriptional regulation. MLLT3 is a critical regulator of hematopoiesis with a role in maintaining the hematopoietic stem or progenitor cell (HSPC) population. Approximately 10% of acute myeloid leukemia (AML) and acute lymphocytic leukemia (ALL) patients harbor a translocation involving MLL (mixed lineage leukemia). In the context of MLL fusion patients with AML and ALL, MLL-AF9 and MLL-ENL fusions are observed in 34 and 31% of the patients, respectively. The intrinsically disordered C-terminal domain of MLLT3 (AHD, ANC1 homology domain) undergoes coupled binding and folding upon interaction with partner proteins AF4, DOT1L, BCOR, and CBX8. Backbone dynamics studies of the complexes suggest a role for dynamics in function. Inhibitors of the interaction of the intrinsically disordered AHD with partner proteins have been described, highlighting the feasibility of targeting intrinsically disordered regions. MLLT1 undergoes phase separation to enhance recruitment of the super elongation complex (SEC) and drive transcription. Mutations in MLLT1 observed in Wilms tumor patients enhance phase separation and transcription to drive an aberrant gene expression program.

Calmodulin binds the N-terminus of the functional amyloid Orb2A inhibiting fibril formation

The cytoplasmic polyadenylation element-binding protein Orb2 is a key regulator of long-term memory (LTM) in Drosophila. The N-terminus of the Orb2 isoform A is required for LTM and forms cross-β fibrils on its own. However, this N-terminus is not part of the core found in ex vivo fibrils. We previously showed that besides forming cross-β fibrils, the N-terminus of Orb2A binds anionic lipid membranes as an amphipathic helix. Here, we show that the Orb2A N-terminus can similarly interact with calcium activated calmodulin (CaM) and that this interaction prevents fibril formation. Because CaM is a known regulator of LTM, this interaction could potentially explain the regulatory role of Orb2A in LTM.

Identification of a modular super-enhancer in murine retinal development

Super-enhancers are expansive regions of genomic DNA comprised of multiple putative enhancers that contribute to the dynamic gene expression patterns during development. This is particularly important in neurogenesis because many essential transcription factors have complex developmental stage- and cell-type specific expression patterns across the central nervous system. In the developing retina, Vsx2 is expressed in retinal progenitor cells and is maintained in differentiated bipolar neurons and Müller glia. A single super-enhancer controls this complex and dynamic pattern of expression. Here we show that deletion of one region disrupts retinal progenitor cell proliferation but does not affect cell fate specification. The deletion of another region has no effect on retinal progenitor cell proliferation but instead leads to a complete loss of bipolar neurons. This prototypical super-enhancer may serve as a model for dissecting the complex gene expression patterns for neurogenic transcription factors during development. Moreover, it provides a unique opportunity to alter expression of individual transcription factors in particular cell types at specific stages of development. This provides a deeper understanding of function that cannot be achieved with traditional knockout mouse approaches.

A Novel WRKY Transcription Factor from Ipomoea trifida, ItfWRKY70, Confers Drought Tolerance in Sweet Potato

WRKY transcription factors are one of the important families in plants, and have important roles in plant growth, abiotic stress responses, and defense regulation. In this study, we isolated a WRKY gene, ItfWRKY70, from the wild relative of sweet potato Ipomoea trifida (H.B.K.) G. Don. This gene was highly expressed in leaf tissue and strongly induced by 20% PEG6000 and 100 μM abscisic acid (ABA). Subcellar localization analyses indicated that ItfWRKY70 was localized in the nucleus. Overexpression of ItfWRKY70 significantly increased drought tolerance in transgenic sweet potato plants. The content of ABA and proline, and the activity of SOD and POD were significantly increased, whereas the content of malondialdehyde (MDA) and H2O2 were decreased in transgenic plants under drought stress. Overexpression of ItfWRKY70 up-regulated the genes involved in ABA biosynthesis, stress-response, ROS-scavenging system, and stomatal aperture in transgenic plants under drought stress. Taken together, these results demonstrated that ItfWRKY70 plays a positive role in drought tolerance by accumulating the content of ABA, regulating stomatal aperture and activating the ROS scavenging system in sweet potato.

Phylogeny, Structure, Functions, and Role of AIRE in the Formation of T-Cell Subsets

It is well known that the most important feature of adaptive immunity is the specificity that provides highly precise recognition of the self, altered-self, and non-self. Due to the high specificity of antigen recognition, the adaptive immune system participates in the maintenance of genetic homeostasis, supports multicellularity, and protects an organism from different pathogens at a qualitatively different level than innate immunity. This seemingly simple property is based on millions of years of evolution that led to the formation of diversification mechanisms of antigen-recognizing receptors and later to the emergence of a system of presentation of the self and non-self antigens. The latter could have a crucial significance because the presentation of nearly complete diversity of auto-antigens in the thymus allows for the "calibration" of the forming repertoires of T-cells for the recognition of self, altered-self, and non-self antigens that are presented on the periphery. The central role in this process belongs to promiscuous gene expression by the thymic epithelial cells that express nearly the whole spectrum of proteins encoded in the genome, meanwhile maintaining their cellular identity. This complex mechanism requires strict control that is executed by several transcription factors. One of the most important of them is AIRE. This noncanonical transcription factor not only regulates the processes of differentiation and expression of peripheral tissue-specific antigens in the thymic medullar epithelial cells but also controls intercellular interactions in the thymus. Besides, it participates in an increase in the diversity and transfer of presented antigens and thus influences the formation of repertoires of maturing thymocytes. Due to these complex effects, AIRE is also called a transcriptional regulator. In this review, we briefly described the history of AIRE discovery, its structure, functions, and role in the formation of antigen-recognizing receptor repertoires, along with other transcription factors. We focused on the phylogenetic prerequisites for the development of modern adaptive immunity and emphasized the importance of the antigen presentation system.

Transcriptional dynamics of gametogenesis in the green seaweed Ulva mutabilis identifies an RWP-RK transcription factor linked to reproduction

BACKGROUND

The molecular mechanism underlying sexual reproduction in land plants is well understood in model plants and is a target for crop improvement. However, unlike land plants, the genetic basis involved in triggering reproduction and gamete formation remains elusive in most seaweeds, which are increasingly viewed as an alternative source of functional food and feedstock for energy applications.

RESULTS

Gametogenesis of Ulva mutabilis, a model organism for green seaweeds, was studied. We analyzed transcriptome dynamics at different time points during gametogenesis following induction of reproduction by fragmentation and removal of sporulation inhibitors. Analyses demonstrated that 45% of the genes in the genome were differentially expressed during gametogenesis. We identified several transcription factors that potentially play a key role in the early gametogenesis of Ulva given the function of their homologs in higher plants and microalgae. In particular, the detailed expression pattern of an evolutionarily conserved transcription factor containing an RWP-RK domain suggested a key role during Ulva gametogenesis.

CONCLUSIONS

Transcriptomic analyses of gametogenesis in the green seaweed Ulva highlight the importance of a conserved RWP-RK transcription factor in the induction of sexual reproduction. The identification of putative master regulators of gametogenesis provides a starting point for further functional characterization.

OUP accepted manuscript

Recognition of mammalian mitochondrial promoters requires the concerted action of mitochondrial RNA polymerase (mtRNAP) and transcription initiation factors TFAM and TFB2M. In this work, we found that transcript slippage results in heterogeneity of the human mitochondrial transcripts in vivo and in vitro. This allowed us to correctly interpret the RNAseq data, identify the bona fide transcription start sites (TSS), and assign mitochondrial promoters for > 50% of mammalian species and some other vertebrates. The divergent structure of the mammalian promoters reveals previously unappreciated aspects of mtDNA evolution. The correct assignment of TSS also enabled us to establish the precise register of the DNA in the initiation complex and permitted investigation of the sequence-specific protein-DNA interactions. We determined the molecular basis of promoter recognition by mtRNAP and TFB2M, which cooperatively recognize bases near TSS in a species-specific manner. Our findings reveal a role of mitochondrial transcription machinery in mitonuclear coevolution and speciation.