The disordered MAX N-terminus modulates DNA binding of the transcription factor MYC:MAX

MAX inactivation deregulates the MYC network and induces neuroendocrine neoplasia in multiple tissues

The MYC transcription factor requires MAX for DNA binding and widespread activation of gene expression in both normal and neoplastic cells. Inactivating mutations in MAX are associated with a subset of neuroendocrine cancers including pheochromocytoma, pituitary adenoma, and small cell lung cancer. Neither the extent nor the mechanisms of MAX tumor suppression are well understood. Deleting Max across multiple mouse neuroendocrine tissues, we find that Max inactivation alone produces pituitary adenomas, while Max inactivation cooperates with Rb1/Trp53 loss to accelerate medullary thyroid C cell and pituitary adenoma development. In the thyroid tumor cell lines, MAX loss triggers a marked shift in genomic occupancy by other members of the MYC network (MNT, MLX, MondoA) supporting metabolism, survival, and proliferation of neoplastic neuroendocrine cells. Our work reveals MAX as a broad suppressor of neuroendocrine tumorigenesis through its ability to maintain a balance of genomic occupancies among the diverse transcription factors in the MYC network.

Versatile roles of disordered transcription factor effector domains in transcriptional regulation

Transcription, a crucial step in the regulation of gene expression, is tightly controlled and involves several essential processes, such as chromatin organization, recognition of the specific genomic sequences, DNA binding, and ultimately recruiting the transcriptional machinery to facilitate transcript synthesis. At the center of this regulation are transcription factors (TFs), which comprise at least one DNA-binding domain (DBD) and an effector domain (ED). Although the structure and function of DBDs have been well studied, our knowledge of the structure and function of effector domains is limited. EDs are of particular importance in generating distinct transcriptional responses between protein members of the same TF family that have similar DBDs and specificities. The study of transcriptional activity conferred by effector domains has traditionally been conducted through examining protein-protein interactions. However, recent research has uncovered alternative mechanisms by which EDs regulate gene expression, such as the formation of condensates that increase the local concentration of transcription factors, cofactors, and coregulated genes, as well as DNA binding. Here, we provide a comprehensive overview of the known roles of transcription factor EDs, with a specific focus on disordered regions. Additionally, we emphasize the significance of intrinsically disordered regions (IDRs) during transcriptional regulation. We examine the mechanisms underlying the establishment and maintenance of transcriptional specificity through the structural properties of predominantly disordered EDs. We then provide a comprehensive overview of the current understanding of these domains, including their physical and chemical characteristics, as well as their functional roles.

Exploring the dynamics and interactions of the N-myc transactivation domain through solution nuclear magnetic resonance spectroscopy

Myc proteins are transcription factors crucial for cell proliferation. They have a C-terminal domain that mediates Max and DNA binding, and an N-terminal disordered region culminating in the transactivation domain (TAD). The TAD participates in many protein-protein interactions, notably with kinases that promote stability (Aurora-A) or degradation (ERK1, GSK3) via the ubiquitin-proteasome system. We probed the structure, dynamics and interactions of N-myc TAD using nuclear magnetic resonance (NMR) spectroscopy following its complete backbone assignment. Chemical shift analysis revealed that N-myc has two regions with clear helical propensity: Trp77-Glu86 and Ala122-Glu132. These regions also have more restricted ps-ns motions than the rest of the TAD, and, along with the phosphodegron, have comparatively high transverse (R2) 15N relaxation rates, indicative of slower timescale dynamics and/or chemical exchange. Collectively these features suggest differential propensities for structure and interaction, either internal or with binding partners, across the TAD. Solution studies on the interaction between N-myc and Aurora-A revealed a previously uncharacterised binding site. The specificity and kinetics of sequential phosphorylation of N-myc by ERK1 and GSK3 were characterised using NMR and resulted in no significant structural changes outside the phosphodegron. When the phosphodegron was doubly phosphorylated, N-myc formed a robust interaction with the Fbxw7-Skp1 complex, but mapping the interaction by NMR suggests a more extensive interface. Our study provides foundational insights into N-myc TAD dynamics and a backbone assignment that will underpin future work on the structure, dynamics, interactions and regulatory post-translational modifications of this key oncoprotein.

MAX inactivation deregulates the MYC network and induces neuroendocrine neoplasia in multiple tissues

The MYC transcription factor requires MAX for DNA binding and widespread activation of gene expression in both normal and neoplastic cells. Surprisingly, inactivating mutations in MAX are associated with a subset of neuroendocrine cancers including pheochromocytoma, pituitary adenoma and small cell lung cancer. Neither the extent nor the mechanisms of MAX tumor suppression are well understood. Delet-ing Max across multiple mouse neuroendocrine tissues, we find Max inactivation alone produces pituitary adenomas while Max loss cooperates with Rb1/Trp53 loss to accelerate medullary thyroid C-cell and pituitary adenoma development. In the thyroid tumor cell lines, MAX loss triggers a striking shift in genomic occupancy by other members of the MYC network (MNT, MLX, MondoA) supporting metabolism, survival and proliferation of neoplastic neuroendocrine cells. Our work reveals MAX as a broad suppressor of neuroendocrine tumorigenesis through its ability to maintain a balance of genomic occupancies among the diverse transcription factors in the MYC network.

Design, synthesis, and activity evaluation of 2-iminobenzimidazoles as c-Myc inhibitors for treating multiple myeloma

Multiple myeloma (MM) is a plasma cell malignancy that remains incurable and poses a significant threat to global public health. The multifunctional transcription factor c-Myc plays a crucial role in various cellular processes and is closely associated with MM progression. As part of the basic-helix-loop-helix-leucine zipper (bHLHZip) family, c-Myc forms heterodimers with its obligate partner Max, binds to the Enhancer-box (E-box) of DNA, and ultimately co-regulates gene expression. Therefore, impeding the capacity for heterodimerization to bind to DNA represents a favored strategy in thwarting c-Myc transcription. In this study, we first synthesized a series of novel 2-iminobenzimidazole derivatives and further estimated their potential anti-MM activity. Notably, among all the derivatives, 5b and 5d demonstrated remarkable inhibitory activity against RPMI-8226 and U266 cells, with IC50 values of 0.85 μM and 0.97 μM for compound 5b, and 0.96 μM and 0.89 μM for compound 5d. Western blot and dual-luciferase reporter assays demonstrated that compounds 5b and 5d effectively suppressed both c-Myc protein expression and transcriptional activity of the c-Myc promoter in RPMI-8226 and U266 cells. Furthermore, these compounds induced apoptosis and G1 cell cycle arrest in the aforementioned MM cells. Molecular docking studies revealed that 5b and 5d exhibited strong binding affinity to the interface between c-Myc/Max and E-box of DNA. Taken together, our findings suggest that further investigations are warranted for potential therapeutic applications of 5b and 5d for c-Myc-related diseases.

High-throughput assay exploiting disorder-to-order conformational switches: application to the proteasomal Rpn10:E6AP complex

Conformational switching is pervasively driven by protein interactions, particularly for intrinsically disordered binding partners. We developed a dually orthogonal fluorescence-based assay to monitor such events, exploiting environmentally sensitive fluorophores. This assay is applied to E3 ligase E6AP, as its AZUL domain induces a disorder-to-order switch in an intrinsically disordered region of the proteasome, the so-named Rpn10 AZUL-binding domain (RAZUL). By testing various fluorophores, we developed an assay appropriate for high-throughput screening of Rpn10:E6AP-disrupting ligands. We found distinct positions in RAZUL for fluorophore labeling with either acrylodan or Atto610, which had disparate spectral responses to E6AP binding. E6AP caused a hypsochromic shift with increased fluorescence of acrylodan-RAZUL while decreasing fluorescence intensity of Atto610-RAZUL. Combining RAZUL labeled with either acrylodan or Atto610 into a common sample achieved robust and orthogonal measurement of the E6AP-induced conformational switch. This approach is generally applicable to disorder-to-order (or vice versa) transitions mediated by molecular interactions.

Analyzing paramagnetic NMR data on target DNA search by proteins using a discrete-state kinetic model for translocation

Before reaching their targets, sequence-specific DNA-binding proteins nonspecifically bind to DNA through electrostatic interactions and stochastically change their locations on DNA. Investigations into the dynamics of DNA-scanning by proteins are nontrivial due to the simultaneous presence of multiple translocation mechanisms and many sites for the protein to nonspecifically bind to DNA. Nuclear magnetic resonance (NMR) spectroscopy can provide information about the target DNA search processes at an atomic level. Paramagnetic relaxation enhancement (PRE) is particularly useful to study how the proteins scan DNA in the search process. Previously, relatively simple two-state or three-state exchange models were used to explain PRE data reflecting the target search process. In this work, using more realistic discrete-state stochastic kinetics models embedded into an NMR master equation, we analyzed the PRE data for the HoxD9 homeodomain interacting with DNA. The kinetic models that incorporate sliding, dissociation, association, and intersegment transfer can reproduce the PRE profiles observed at some different ionic strengths. The analysis confirms the previous interpretation of the PRE data and shows that the protein's probability distribution among nonspecific sites is nonuniform during the target DNA search process.

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.

Commonly asked questions about transcriptional activation domains

Eukaryotic transcription factors activate gene expression with their DNA-binding domains and activation domains. DNA-binding domains bind the genome by recognizing structurally related DNA sequences; they are structured, conserved, and predictable from protein sequences. Activation domains recruit chromatin modifiers, coactivator complexes, or basal transcriptional machinery via structurally diverse protein-protein interactions. Activation domains and DNA-binding domains have been called independent, modular units, but there are many departures from modularity, including interactions between these regions and overlap in function. Compared to DNA-binding domains, activation domains are poorly understood because they are poorly conserved, intrinsically disordered, and difficult to predict from protein sequences. This review, organized around commonly asked questions, describes recent progress that the field has made in understanding the sequence features that control activation domains and predicting them from sequence.

Dynamics and interactions of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are widespread in eukaryotes and participate in a variety of important cellular processes. Numerous studies using state-of-the-art experimental and theoretical methods have advanced our understanding of IDPs and revealed that disordered regions engage in a large repertoire of intra- and intermolecular interactions through their conformational dynamics, thereby regulating many intracellular functions in concert with folded domains. The mechanisms by which IDPs interact with their partners are diverse, depending on their conformational propensities, and include induced fit, conformational selection, and their mixtures. In addition, IDPs are implicated in many diseases, and progress has been made in designing inhibitors of IDP-mediated interactions. Here we review these recent advances with a focus on the dynamics and interactions of IDPs.

Intrinsically Disordered Regions in the Transcription Factor MYC:MAX Modulate DNA Binding via Intramolecular Interactions

The basic helix-loop-helix leucine zipper (bHLH-LZ) transcription factor (TF) MYC is in large part an intrinsically disordered oncoprotein. In complex with its obligate heterodimerization partner MAX, MYC preferentially binds E-Box DNA sequences (CANNTG). At promoters containing these sequence motifs, MYC controls fundamental cellular processes such as cell cycle progression, metabolism, and apoptosis. A vast network of proteins in turn regulates MYC function via intermolecular interactions. In this work, we establish another layer of MYC regulation by intramolecular interactions. We used nuclear magnetic resonance (NMR) spectroscopy to identify and map multiple binding sites for the C-terminal MYC:MAX DNA-binding domain (DBD) on the intrinsically disordered regions (IDRs) in the MYC N-terminus. We find that these binding events in trans are driven by electrostatic attraction, that they have distinct affinities, and that they are competitive with DNA binding. Thereby, we observe the strongest effects for the N-terminal MYC box 0 (Mb0), a conserved motif involved in MYC transactivation and target gene induction. We prepared recombinant full-length MYC:MAX complex and demonstrate that the interactions identified in this work are also relevant in cis, i.e., as intramolecular interactions. These findings are supported by surface plasmon resonance (SPR) experiments, which revealed that intramolecular IDR:DBD interactions in MYC decelerate the association of MYC:MAX complexes to DNA. Our work offers new insights into how bHLH-LZ TFs are regulated by intramolecular interactions, which open up new possibilities for drug discovery.

Unraveling the Behavior of Intrinsically Disordered Protein c-Myc: A Study Utilizing Gaussian-Accelerated Molecular Dynamics

The protein c-Myc is a transcription factor that remains largely intrinsically disordered and is known to be involved in various biological processes and is overexpressed in various cancers, making it an attractive drug target. However, intrinsically disordered proteins such as c-Myc do not show funnel-like basins in their free-energy landscapes; this makes their druggability a challenge. For the first time, we propose a heterodimer model of c-Myc/Max in full length in this work. We used Gaussian-accelerated molecular dynamics (GaMD) simulations to explore the behavior of c-Myc and its various regions, including the transactivation domain (TAD) and the basic helix-loop-helix-leucine-zipper (bHLH-Zipper) motif in three different conformational states: (a) monomeric c-Myc, (b) c-Myc when bound to its partner protein, Max, and (c) when Max was removed after binding. We analyzed the GaMD trajectories using root-mean-square deviation (RMSD), radius of gyration, root-mean-square fluctuation, and free-energy landscape (FEL) calculations to elaborate the behaviors of these regions. The results showed that the monomeric c-Myc structure showed a higher RMSD fluctuation as compared with the c-Myc/Max heterodimer in the bHLH-Zipper motif. This indicated that the bHLH-Zipper motif of c-Myc is more stable when it is bound to Max. The TAD region in both monomeric and Max-bound states showed similar plasticity in terms of RMSD. We also conducted residue decomposition calculations and showed that the c-Myc and Max interaction could be driven mainly by electrostatic interactions and the residues Arg299, Ile403, and Leu420 seemed to play important roles in the interaction. Our work provides insights into the behavior of c-Myc and its regions that could support the development of drugs that target c-Myc and other intrinsically disordered proteins.

Collaborating single-cell and bulk RNA sequencing for comprehensive characterization of the intratumor heterogeneity and prognostic model development for bladder cancer

INTRODUCTION

Gaining a deeper insight into the single-cell RNA sequencing (scRNA-seq) results of bladder cancer (BLCA) provides a transcriptomic profiling of individual cancer cells, which may disclose the molecular mechanisms involved in BLCA carcinogenesis.

METHODS

scRNA data were obtained from GSE169379 dataset. We used the InferCNV software to determine the copy number variant (CNV) with normal epithelial cells serving as the reference, and performed the pseudo-timing analysis on subsets of epithelial cell using Monocle3 software. Transcription factor analysis was conducted using the Dorothea software. Intercellular communication analysis was performed using the Liana software. Cox analysis and LASSO regression were applied to establish a prognostic model.

RESULTS

We investigated the heterogeneity of tumors in four distinct cell types of BLCA cancer, namely immune cells, endothelial cells, epithelial cells, and fibroblasts. We evaluated the transcription factor activity of different immune cells in BLCA and identified significant enrichment of TCF7 and TBX21 in CD8+ T cells. Additionally, we identified two distinct subtypes of cancer-associated fibroblasts (CAFs), namely iCAFs and myoCAFs, which exhibited distinct communication patterns. Using sub-cluster and cell trajectory analyses, we identified different states of normal-to-malignant cell transformation in epithelial cells. TF analysis further revealed high activation of MYC and SOX2 in tumor cells. Finally, we identified five model genes (SLCO3A1, ANXA1, TENM3, EHBP1, LSAMP) for the development of a prognostic model, which demonstrated high effectiveness in stratifying patients across seven different cohorts.

CONCLUSIONS

We have developed a prognostic model that has demonstrated significant efficacy in stratifying patients with BLCA.

Direct measurements of biomolecular electrostatics through experiments

Biomolecular electrostatics has been a subject of computational investigations based on 3D structures. This situation is changing because emerging experimental tools allow us to quantitatively investigate biomolecular electrostatics without any use of structure information. Now, electrostatic potentials around biomolecules can directly be measured for many residues simultaneously by nuclear magnetic resonance (NMR) spectroscopy. This NMR method can be used to study electrostatic aspects of various processes, including macromolecular association and liquid-liquid phase separation. Applications to structurally flexible biomolecules such as intrinsically disordered proteins are particularly useful. The new tools also facilitate examination of theoretical models and methods for biomolecular electrostatics.

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.