Solution structure of the KIX domain of CBP bound to the transactivation domain of CREB: a model for activator:coactivator interactions

Recognition of the CREB-Binding Protein/Mixed Lineage Leukemia Interface by Sulfonyl-γ-AApeptides: Beyond Mimicry of Mixed Lineage Leukemia

The kinase-inducible (KIX) domain, a structured region within the transcriptional coactivators CREB-binding protein (CBP) and p300, serves as a docking site for multiple transcription factors, including the mixed lineage leukemia (MLL) protein. The MLL-KIX protein-protein interaction (PPI) has been implicated in various diseases, such as leukemia, cancer, and neurodegenerative disorders, including Alzheimer's disease. In this study, we developed a series of MLL-mimicking peptidomimetic foldamers based on a sulfonyl-γ-AApeptide backbone. These helical foldamers successfully mimic the folded transactivation domain (TAD) of MLL and bind the KIX domain with a high affinity. Consistent with the notion that MLL-KIX interaction can allosterically enhance CREB signaling, we found that the sulfonyl-γ-AApeptides strongly stimulate CREB-dependent gene expression. Moreover, they exhibit no significant cytotoxicity under the tested conditions, demonstrating both binding specificity and therapeutic potential for targeting CREB-related pathways in neurodegenerative disorders.

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.

The evolution of the Amber additive protein force field: History, current status, and future

Molecular dynamics simulations are pivotal in elucidating the intricate properties of biological molecules. Nonetheless, the reliability of their outcomes hinges on the precision of the molecular force field utilized. In this perspective, we present a comprehensive review of the developmental trajectory of the Amber additive protein force field, delving into researchers' persistent quest for higher precision force fields and the prevailing challenges. We detail the parameterization process of the Amber protein force fields, emphasizing the specific improvements and retained features in each version compared to their predecessors. Furthermore, we discuss the challenges that current force fields encounter in balancing the interactions of protein-protein, protein-water, and water-water in molecular dynamics simulations, as well as potential solutions to overcome these issues.

The transcription factor CREB regulates epithelial-mesenchymal transition of lens epithelial cells by phosphorylation-dependent and phosphorylation-independent mechanisms

Epithelial mesenchymal transition (EMT) of lens epithelial cells (LECs) is one of the most important pathogenic mechanisms in lens fibrotic disorders, and the regulatory mechanisms of EMT have not been fully understood. Here, we demonstrate that the cAMP-response element binding protein (CREB) can regulate lens EMT in a phosphorylation-dependent and phosphorylation-independent manners with dual mechanisms. First, CREB-S133 phosphorylation is implicated in TGFβ-induced EMT of mouse LECs and also in injury-induced mouse anterior subcapsular cataract model. The interaction between CREB and p300 is necessary for CREB regulation of TGFβ-induced EMT, since inhibition of CREB-p300 interaction and p300 knockdown led to attenuated expression of mesenchymal genes. Second, S133A-CREB, a mutant mimicking constant dephosphorylation at S133, exhibits notable occupancy in the enhancers of mesenchymal genes and confers robust transcription activity on EMT genes. Introduction of R314A mutation in S133A-CREB, which abolishes the interaction between S133A-CREB and its co-activator, cAMP-regulated transcriptional co-activators led to substantial suppression of mesenchymal gene expression in mouse LECs. Taken together, our results showed that CREB regulates lens EMT in dual mechanisms and that the S133A-CREB acts as a novel transcription factor. Mechanistically, CREB interacts with p300 in a S133 phosphorylation-dependent manner to positively regulate lens EMT genes. In contrast, S133A-CREB interacts with cAMP-regulated transcriptional co-activators to confer a robust activation of lens EMT genes.

Cooperative folding as a molecular switch in an evolved antibody binder

Designing proteins with tunable activities from easily accessible external cues remains a biotechnological challenge. Here, we set out to create a small antibody-binding domain equipped with a molecular switch inspired by the allosteric response to calcium seen in naturally derived proteins like calmodulin. We have focused on one of the three domains of Protein G that show inherent affinity to antibodies. By combining a semi-rational protein design with directed evolution, we engineered novel variants containing a calcium-binding loop rendering the inherent antibody affinity calcium-dependent. The evolved variants resulted from a designed selection strategy subjecting them to negative and positive selection pressures focused on conditional antibody binding. Hence, these variants contains molecular "on/off" switches, controlling the target affinity towards antibody fragments simply by the presence or absence of calcium. From NMR spectroscopy we found that the molecular mechanism underlying the evolved switching behavior was a coupled calcium-binding and folding event where the target binding surface was intact and functional only in the presence of bound calcium. Notably, it was observed that the response to the employed selection pressures gave rise to the evolution of a cooperative folding mechanism. This observation illustrates why the cooperative folding reaction is an effective solution seen repeatedly in the natural evolution of fine-tuned macromolecular recognition. Engineering binding moieties to confer conditional target interaction has great potential due to the exquisite interaction control that is tunable to application requirements. Improved understanding of the molecular mechanisms behind regulated interactions is crucial to unlock how to engineer switchable proteins useful in a variety of biotechnological applications.

Bipartite binding interface recruiting HP1 to chromosomal passenger complex at inner centromeres

Maintenance of ploidy depends on the mitotic kinase Aurora B, the catalytic subunit of the chromosomal passenger complex (CPC) whose proficient activity is supported by HP1 enriched at inner centromeres. HP1 is known to associate with INCENP of the CPC in a manner that depends on the PVI motif conserved across HP1 interactors. Here, we found that the interaction of INCENP with HP1 requires not only the PVI motif but also its C-terminally juxtaposed domain. Remarkably, these domains conditionally fold the β-strand (PVI motif) and the α-helix from a disordered sequence upon HP1 binding and render INCENP with high affinity to HP1. This bipartite binding domain termed SSH domain (Structure composed of Strand and Helix) is necessary and sufficient to attain a predominant interaction of HP1 with INCENP. These results identify a unique HP1-binding module in INCENP that ensures enrichment of HP1 at inner centromeres, Aurora B activity, and thereby mitotic fidelity.

Multi-faceted regulation of CREB family transcription factors

cAMP response element-binding protein (CREB) is a ubiquitously expressed nuclear transcription factor, which can be constitutively activated regardless of external stimuli or be inducibly activated by external factors such as stressors, hormones, neurotransmitters, and growth factors. However, CREB controls diverse biological processes including cell growth, differentiation, proliferation, survival, apoptosis in a cell-type-specific manner. The diverse functions of CREB appear to be due to CREB-mediated differential gene expression that depends on cAMP response elements and multi-faceted regulation of CREB activity. Indeed, the transcriptional activity of CREB is controlled at several levels including alternative splicing, post-translational modification, dimerization, specific transcriptional co-activators, non-coding small RNAs, and epigenetic regulation. In this review, we present versatile regulatory modes of CREB family transcription factors and discuss their functional consequences.

HIF1α/ATF3 partake in PGK1 K191/K192 succinylation by modulating P4HA1/succinate signaling in glioblastoma

BACKGROUND

Hypoxia is a pathological hallmark in most cancers, including glioblastoma (GBM). Hypoxic signaling activation and post-translational modification (PTM) of oncogenic proteins are well-studied in cancers. Accumulating studies indicate glycolytic enzyme PGK1 plays a crucial role in tumorigenesis, yet the underlying mechanisms remain unknown.

METHODS

We first used ChIP assays to uncover the crosstalk between HIF1α and ATF3 and their roles in P4HA1 regulation. Protein degradation analysis, LC-MS/MS, and in vitro succinate production assays were performed to examine the effect of protein succinylation on GBM pathology. Seahorse assay measured the effects of PGK1 succinylation at K191/K192 or its mutants on glucose metabolism. We utilized an in vivo intracranial mouse model for biochemical studies to elucidate the impact of ATF3 and P4HA1 on aerobic glycolysis and the tumor immune microenvironment.

RESULTS

We demonstrated that HIF1α and ATF3 positively and negatively regulate the transcription of P4HA1, respectively, leading to an increased succinate production and increased activation of HIF1α signaling. P4HA1 expression elevated the succinate concentration, resulting in the enhanced succinylation of PGK1 at the K191 and K192 sites. Inhibition of proteasomal degradation of PGK1 by succinylation significantly increased aerobic glycolysis to generate lactate. Furthermore, ATF3 overexpression and P4HA1 knockdown reduced succinate and lactate levels in GBM cells, inhibiting immune responses and tumor growth.

CONCLUSIONS

Together, our study demonstrates that HIF1α/ATF3 participated in P4HA1/succinate signaling, which is the major regulator of succinate biosynthesis and PGK1 succinylation at K191 and K192 sites in GBM. The P4HA1/succinate pathway might be a novel and promising target for aerobic glycolysis in GBM.

p300/CBP degradation is required to disable the active AR enhanceosome in prostate cancer

Prostate cancer is an exemplar of an enhancer-binding transcription factor-driven disease. The androgen receptor (AR) enhanceosome complex comprised of chromatin and epigenetic coregulators assembles at enhancer elements to drive disease progression. The paralog lysine acetyltransferases p300 and CBP deposit histone marks that are associated with enhancer activation. Here, we demonstrate that p300/CBP are determinant cofactors of the active AR enhanceosome in prostate cancer. Histone H2B N-terminus multisite lysine acetylation (H2BNTac), which is exclusively reliant on p300/CBP catalytic function, marked active enhancers and was notably elevated in prostate cancer lesions relative to the adjacent benign epithelia. Degradation of p300/CBP rapidly depleted acetylation marks associated with the active AR enhanceosome, which was only partially phenocopied by inhibition of their reader bromodomains. Notably, H2BNTac was effectively abrogated only upon p300/CBP degradation, which led to a stronger suppression of p300/CBP-dependent oncogenic gene programs relative to bromodomain inhibition or the inhibition of its catalytic domain. In vivo experiments using an orally active p300/CBP proteolysis targeting chimera (PROTAC) degrader (CBPD-409) showed that p300/CBP degradation potently inhibited tumor growth in preclinical models of castration-resistant prostate cancer and synergized with AR antagonists. While mouse p300/CBP orthologs were effectively degraded in host tissues, prolonged treatment with the PROTAC degrader was well tolerated with no significant signs of toxicity. Taken together, our study highlights the pivotal role of p300/CBP in maintaining the active AR enhanceosome and demonstrates how target degradation may have functionally distinct effects relative to target inhibition, thus supporting the development of p300/CBP degraders for the treatment of advanced prostate cancer.

Assessment of machine-learning predictions for the Mediator complex subunit MED25 ACID domain interactions with transactivation domains

The human Mediator complex subunit MED25 binds transactivation domains (TADs) present in various cellular and viral proteins using two binding interfaces, named H1 and H2, which are found on opposite sides of its ACID domain. Here, we use and compare deep learning methods to characterize human MED25-TAD interfaces and assess the predicted models to published experimental data. For the H1 interface, AlphaFold produces predictions with high-reliability scores that agree well with experimental data, while the H2 interface predictions appear inconsistent, preventing reliable binding modes. Despite these limitations, we experimentally assess the validity of MED25 interface predictions with the viral transcriptional activators Lana-1 and IE62. AlphaFold predictions also suggest the existence of a unique hydrophobic pocket for the Arabidopsis MED25 ACID domain.

Experimental methods to study the structure and dynamics of intrinsically disordered regions in proteins

Eukaryotic proteins often feature long stretches of amino acids that lack a well-defined three-dimensional structure and are referred to as intrinsically disordered proteins (IDPs) or regions (IDRs). Although these proteins challenge conventional structure-function paradigms, they play vital roles in cellular processes. Recent progress in experimental techniques, such as NMR spectroscopy, single molecule FRET, high speed AFM and SAXS, have provided valuable insights into the biophysical basis of IDP function. This review discusses the advancements made in these techniques particularly for the study of disordered regions in proteins. In NMR spectroscopy new strategies such as 13C detection, non-uniform sampling, segmental isotope labeling, and rapid data acquisition methods address the challenges posed by spectral overcrowding and low stability of IDPs. The importance of various NMR parameters, including chemical shifts, hydrogen exchange rates, and relaxation measurements, to reveal transient secondary structures within IDRs and IDPs are presented. Given the high flexibility of IDPs, the review outlines NMR methods for assessing their dynamics at both fast (ps-ns) and slow (μs-ms) timescales. IDPs exert their functions through interactions with other molecules such as proteins, DNA, or RNA. NMR-based titration experiments yield insights into the thermodynamics and kinetics of these interactions. Detailed study of IDPs requires multiple experimental techniques, and thus, several methods are described for studying disordered proteins, highlighting their respective advantages and limitations. The potential for integrating these complementary techniques, each offering unique perspectives, is explored to achieve a comprehensive understanding of IDPs.

Inhibition of CREB Binding and Function with a Dual-Targeting Ligand

CBP/p300 is a master transcriptional coactivator that regulates gene activation by interacting with multiple transcriptional activators. Dysregulation of protein-protein interactions (PPIs) between the CBP/p300 KIX domain and its activators is implicated in a number of cancers, including breast, leukemia, and colorectal cancer. However, KIX is typically considered "undruggable" because of its shallow binding surfaces lacking both significant topology and promiscuous binding profiles. We previously reported a dual-targeting peptide (MybLL-tide) that inhibits the KIX-Myb interaction with excellent specificity and potency. Here, we demonstrate a branched, second-generation analogue, CREBLL-tide, that inhibits the KIX-CREB PPI with higher potency and selectivity. Additionally, the best of these CREBLL-tide analogues shows excellent and selective antiproliferation activity in breast cancer cells. These results indicate that CREBLL-tide is an effective tool for assessing the role of KIX-activator interactions in breast cancer and expanding the dual-targeting strategy for inhibiting KIX and other coactivators that contain multiple binding surfaces.

Rational optimization of a transcription factor activation domain inhibitor

Transcription factors are among the most attractive therapeutic targets but are considered largely 'undruggable' in part due to the intrinsically disordered nature of their activation domains. Here we show that the aromatic character of the activation domain of the androgen receptor, a therapeutic target for castration-resistant prostate cancer, is key for its activity as transcription factor, allowing it to translocate to the nucleus and partition into transcriptional condensates upon activation by androgens. On the basis of our understanding of the interactions stabilizing such condensates and of the structure that the domain adopts upon condensation, we optimized the structure of a small-molecule inhibitor previously identified by phenotypic screening. The optimized compounds had more affinity for their target, inhibited androgen-receptor-dependent transcriptional programs, and had an antitumorigenic effect in models of castration-resistant prostate cancer in cells and in vivo. These results suggest that it is possible to rationally optimize, and potentially even to design, small molecules that target the activation domains of oncogenic transcription factors.

Glutamine-rich regions of the disordered CREB transactivation domain mediate dynamic intra- and intermolecular interactions

The cyclic AMP response element (CRE) binding protein (CREB) is a transcription factor that contains a 280-residue N-terminal transactivation domain and a basic leucine zipper that mediates interaction with DNA. The transactivation domain comprises three subdomains, the glutamine-rich domains Q1 and Q2 and the kinase inducible activation domain (KID). NMR chemical shifts show that the isolated subdomains are intrinsically disordered but have a propensity to populate local elements of secondary structure. The Q1 and Q2 domains exhibit a propensity for formation of short β-hairpin motifs that function as binding sites for glutamine-rich sequences. These motifs mediate intramolecular interactions between the CREB Q1 and Q2 domains as well as intermolecular interactions with the glutamine-rich Q1 domain of the TATA-box binding protein associated factor 4 (TAF4) subunit of transcription factor IID (TFIID). Using small-angle X-ray scattering, NMR, and single-molecule Förster resonance energy transfer, we show that the Q1, Q2, and KID regions remain dynamically disordered in a full-length CREB transactivation domain (CREBTAD) construct. The CREBTAD polypeptide chain is largely extended although some compaction is evident in the KID and Q2 domains. Paramagnetic relaxation enhancement reveals transient long-range contacts both within and between the Q1 and Q2 domains while the intervening KID domain is largely devoid of intramolecular interactions. Phosphorylation results in expansion of the KID domain, presumably making it more accessible for binding the CBP/p300 transcriptional coactivators. Our study reveals the complex nature of the interactions within the intrinsically disordered transactivation domain of CREB and provides molecular-level insights into dynamic and transient interactions mediated by the glutamine-rich domains.

Garcinolic Acid Distinguishes Between GACKIX Domains and Modulates Interaction Networks

Natural products are often uniquely suited to modulate protein-protein interactions (PPIs) due to their architectural and functional group complexity relative to synthetic molecules. Here we demonstrate that the natural product garcinolic acid allosterically blocks the CBP/p300 KIX PPI network and displays excellent selectivity over related GACKIX motifs. It does so via a strong interaction (KD 1 μM) with a non-canonical binding site containing a structurally dynamic loop in CBP/p300 KIX. Garcinolic acid engages full-length CBP in the context of the proteome and in doing so effectively inhibits KIX-dependent transcription in a leukemia model. As the most potent small-molecule KIX inhibitor yet reported, garcinolic acid represents an important step forward in the therapeutic targeting of CBP/p300.

From Immunogenic Peptides to Intrinsically Disordered Proteins

It is hard to evaluate the role of individual mentors in the genesis of important ideas. In the case of our realization that proteins do not have to be stably folded to be functional, the influence of Richard Lerner and our collaborative work in the 1980s on the conformations of immunogenic peptides provided a base level of thinking about the nature of polypeptides in water solutions that led us to formulate and develop our ideas on the importance of intrinsic disorder in proteins. This review describes how the insights gained into the behavior of peptides led directly to the realization that proteins were not only capable of being functional while disordered, but also that disorder provided a distinct functional advantage in many important cellular processes.

Genome-wide analysis of KIX gene family for organ size regulation in soybean (Glycine max L.)

The KIX domain, conserved among various nuclear and co-activator factors, acts as a binding site that interacts with other transcriptional activators and co-activators, playing a crucial role in gene expression regulation. In plants, the KIX domain is involved in plant hormone signaling, stress response regulation, cell cycle control, and differentiation, indicating its potential relevance to crop productivity. This study aims to identify and characterize KIX domains within the soybean (Glycine max L.) genome to predict their potential role in improving crop productivity. The conservation and evolutionary history of the KIX domains were explored in 59 plant species, confirming the presence of the KIX domains in diverse plants. Specifically, 13 KIX domains were identified within the soybean genome and classified into four main groups, namely GmKIX8/9, GmMED15, GmHAC, and GmRECQL, through sequence alignment, structural analysis, and phylogenetic tree construction. Association analysis was performed between KIX domain haplotypes and soybean seed-related agronomic traits using re-sequencing data from a core collection of 422 accessions. The results revealed correlations between SNP variations observed in GmKIX8-3 and GmMED15-4 and soybean seed phenotypic traits. Additionally, transcriptome analysis confirmed significant expression of the KIX domains during the early stages of soybean seed development. This study provides the first characterization of the structural, expression, genomic haplotype, and molecular features of the KIX domain in soybean, offering a foundation for functional analysis of the KIX domain in soybean and other plants.

The MIO1-MtKIX8 module regulates the organ size in Medicago truncatula

Plant organ size is an important agronomic trait tightly related to crop yield. However, the molecular mechanisms underlying organ size regulation remain largely unexplored in legumes. We previously characterized a key regulator F-box protein MINI ORGAN1 (MIO1)/SMALL LEAF AND BUSHY1 (SLB1), which controls plant organ size in the model legume Medicago truncatula. In order to further dissect the molecular mechanism, MIO1 was used as the bait to screen its interacting proteins from a yeast library. Subsequently, a KIX protein, designated MtKIX8, was identified from the candidate list. The interaction between MIO1 and MtKIX8 was confirmed further by Y2H, BiFC, split-luciferase complementation and pull-down assays. Phylogenetic analyses indicated that MtKIX8 is highly homologous to Arabidopsis KIX8, which negatively regulates organ size. Moreover, loss-of-function of MtKIX8 led to enlarged leaves and seeds, while ectopic expression of MtKIX8 in Arabidopsis resulted in decreased cotyledon area and seed weight. Quantitative reverse-transcription PCR and in situ hybridization showed that MtKIX8 is expressed in most developing organs. We also found that MtKIX8 serves as a crucial molecular adaptor, facilitating interactions with BIG SEEDS1 (BS1) and MtTOPLESS (MtTPL) proteins in M. truncatula. Overall, our results suggest that the MIO1-MtKIX8 module plays a significant and conserved role in the regulation of plant organ size. This module could be a good target for molecular breeding in legume crops and forages.

NFYC-37 promotes tumor growth by activating the mevalonate pathway in bladder cancer

Dysregulation of transcription is a hallmark of cancer, including bladder cancer (BLCA). CRISPR-Cas9 screening using a lentivirus library with single guide RNAs (sgRNAs) targeting human transcription factors and chromatin modifiers is used to reveal genes critical for the proliferation and survival of BLCA cells. As a result, the nuclear transcription factor Y subunit gamma (NFYC)-37, but not NFYC-50, is observed to promote cell proliferation and tumor growth in BLCA. Mechanistically, NFYC-37 interacts with CBP and SREBP2 to activate mevalonate pathway transcription, promoting cholesterol biosynthesis. However, NFYC-50 recruits more of the arginine methyltransferase CARM1 than NFYC-37 to methylate CBP, which prevents the CBP-SREBP2 interaction and subsequently inhibits the mevalonate pathway. Importantly, statins targeting the mevalonate pathway can suppress NFYC-37-induced cell proliferation and tumor growth, indicating the need for conducting a clinical trial with statins for treating patients with BLCA and high NFYC-37 levels, as most patients with BLCA have high NFYC-37 levels.

Rational peptide design for inhibition of the KIX-MLL interaction

The kinase-inducible domain interacting (KIX) domain is an integral part of the general transcriptional coactivator CREB-binding protein, and has been associated with leukemia, cancer, and various viral diseases. Hence, the KIX domain has attracted considerable attention in drug discovery and development. Here, we rationally designed a KIX inhibitor using a peptide fragment corresponding to the transactivation domain (TAD) of the transcriptional activator, mixed-lineage leukemia protein (MLL). We performed theoretical saturation mutagenesis using the Rosetta software to search for mutants expected to bind KIX more tightly than the wild-type MLL TAD. Mutant peptides with higher helical propensities were selected for experimental characterization. We found that the T2857W mutant of the MLL TAD peptide had the highest binding affinity for KIX compared to the other 12 peptides designed in this study. Moreover, the peptide had a high inhibitory effect on the KIX-MLL interaction with a half-maximal inhibitory concentration close to the dissociation constant for this interaction. To our knowledge, this peptide has the highest affinity for KIX among all previously reported inhibitors that target the MLL site of KIX. Thus, our approach may be useful for rationally developing helical peptides that inhibit protein-protein interactions implicated in the progression of various diseases.

Coupling of binding and differential subdomain folding of the intrinsically disordered transcription factor CREB

The cyclic AMP response element binding protein (CREB) contains a basic leucine zipper motif (bZIP) that forms a coiled coil structure upon dimerization and specific DNA binding. Although this state is well characterized, key features of CREB bZIP binding and folding are not well understood. We used single-molecule Förster resonance energy transfer (smFRET) to probe conformations of CREB bZIP subdomains. We found differential folding of the basic region and leucine zipper in response to different binding partners; a strong and previously unreported DNA-independent dimerization affinity; folding upon binding to nonspecific DNA; and evidence of long-range interdomain interactions in full-length CREB that modulate DNA binding. These studies provide new insights into DNA binding and dimerization and have implications for CREB function.

Re-Balancing Replica Exchange with Solute Tempering for Sampling Dynamic Protein Conformations

Replica exchange with solute tempering (REST) is a highly effective variant of replica exchange for enhanced sampling in explicit solvent simulations of biomolecules. By scaling the Hamiltonian for a selected "solute" region of the system, REST effectively applies tempering only to the degrees of freedom of interest but not the rest of the system ("solvent"), allowing fewer replicas for covering the same temperature range. A key consideration of REST is how the solute-solvent interactions are scaled together with the solute-solute interactions. Here, we critically evaluate the performance of the latest REST2 protocol for sampling large-scale conformation fluctuations of intrinsically disordered proteins (IDPs). The results show that REST2 promotes artificial protein conformational collapse at high effective temperatures, which seems to be a designed feature originally to promote the sampling of reversible folding of small proteins. The collapse is particularly severe with larger IDPs, leading to replica segregation in the effective temperature space and hindering effective sampling of large-scale conformational changes. We propose that the scaling of the solute-solvent interactions can be treated as free parameters in REST, which can be tuned to control the solute conformational properties (e.g., chain expansion) at different effective temperatures and achieve more effective sampling. To this end, we derive a new REST3 protocol, where the strengths of the solute-solvent van der Waals interactions are recalibrated to reproduce the levels of protein chain expansion at high effective temperatures. The efficiency of REST3 is examined using two IDPs with nontrivial local and long-range structural features, including the p53 N-terminal domain and the kinase inducible transactivation domain of transcription factor CREB. The results suggest that REST3 leads to a much more efficient temperature random walk and improved sampling efficiency, which also further reduces the number of replicas required. Nonetheless, our analysis also reveals significant challenges of relying on tempering alone for sampling large-scale conformational fluctuations of disordered proteins. It is likely that more efficient sampling protocols will require incorporating more sophisticated Hamiltonian replica exchange schemes in addition to tempering.

Design and Biochemical Characterization of Peptidic Inhibitors of the Myb/p300 Interaction

The Myb transcription factor is involved in the proliferation of hematopoietic cells, and deregulation of its expression can lead to cancers such as leukemia. Myb interacts with various proteins, including the histone acetyltransferases p300 and CBP. Myb binds to a small domain of p300, the KIX domain (p300^KIX), and inhibiting this interaction is a potential new drug discovery strategy in oncology. The available structures show that Myb binds to a very shallow pocket of the KIX domain, indicating that it might be challenging to identify inhibitors of this interaction. Here, we report the design of Myb-derived peptides which interact with p300^KIX. We show that by mutating only two Myb residues that bind in or near a hotspot at the surface of p300^KIX, it is possible to obtain single-digit nanomolar peptidic inhibitors of the Myb/p300^KIX interaction that bind 400-fold tighter to p300^KIX than wildtype Myb. These findings suggest that it might also be possible to design potent low molecular-weight compounds to disrupt the Myb/p300^KIX interaction.

Functional Analyses of Rare Germline Missense BRCA1 Variants Located within and outside Protein Domains with Known Functions

The BRCA1 protein is implicated in numerous important cellular processes to prevent genomic instability and tumorigenesis, and pathogenic germline variants predispose carriers to hereditary breast and ovarian cancer (HBOC). Most functional studies of missense variants in BRCA1 focus on variants located within the Really Interesting New Gene (RING), coiled-coil and BRCA1 C-terminal (BRCT) domains, and several missense variants in these regions have been shown to be pathogenic. However, the majority of these studies focus on domain specific assays, and have been performed using isolated protein domains and not the full-length BRCA1 protein. Furthermore, it has been suggested that BRCA1 missense variants located outside domains with known function are of no functional importance, and could be classified as (likely) benign. However, very little is known about the role of the regions outside the well-established domains of BRCA1, and only a few functional studies of missense variants located within these regions have been published. In this study, we have, therefore, functionally evaluated the effect of 14 rare BRCA1 missense variants considered to be of uncertain clinical significance, of which 13 are located outside the well-established domains and one within the RING domain. In order to investigate the hypothesis stating that most BRCA1 variants located outside the known protein domains are benign and of no functional importance, multiple protein assays including protein expression and stability, subcellular localisation and protein interactions have been performed, utilising the full-length protein to better mimic the native state of the protein. Two variants located outside the known domains (p.Met297Val and p.Asp1152Asn) and one variant within the RING domain (p.Leu52Phe) were found to make the BRCA1 protein more prone to proteasome-mediated degradation. In addition, two variants (p.Leu1439Phe and p.Gly890Arg) also located outside known domains were found to have reduced protein stability compared to the wild type protein. These findings indicate that variants located outside the RING, BRCT and coiled-coiled domains could also affect the BRCA1 protein function. For the nine remaining variants, no significant effects on BRCA1 protein functions were observed. Based on this, a reclassification of seven variants from VUS to likely benign could be suggested.

SARS-CoV-2 infection activates CREB/CBP in cellular cyclic AMP-dependent pathways

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a global coronavirus disease 2019 (COVID-19) pandemic that has affected the lives of billions of individuals. However, the host-virus interactions still need further investigation to reveal the underling mechanism of SARS-CoV-2 pathogenesis. Here, transcriptomics analysis of SARS-CoV-2 infection highlighted possible correlation between host-associated signaling pathway and virus. In detail, cAMP-protein kinase (PKA) pathway has an essential role in SARS-CoV-2 infection, followed by the interaction between cyclic AMP response element binding protein (CREB) and CREB-binding protein (CBP) could be induced and leading to the enhancement of CREB/CBP transcriptional activity. The replication of Delta and Omicron BA.5 were inhibited by about 49.4% and 44.7% after knockdown of CREB and CBP with small interfering RNAs, respectively. Furthermore, a small organic molecule naphthol AS-E (nAS-E), which targets on the interaction between CREB and CBP, potently inhibited SARS-CoV-2 wild-type (WT) infection with comparable the half-maximal effective concentration (EC₅₀ ) 1.04 μM to Remdesivir 0.57 μM. Compared with WT virus, EC₅₀ in Calu-3 cells against Delta, Omicron BA.2, and Omicron BA.5 were, on average, 1.5-fold, 1.1-fold, and 1.5-fold higher, respectively, nAS-E had a satisfied antiviral effect against Omicron variants. Taken together, our study demonstrated the importance of CREB/CBP induced by cAMP-PKA pathway during SARS-CoV-2 infection, and further provided a novel CREB/CBP interaction therapeutic drug targets for COVID-19.

Protein dynamics: The future is bright and complicated!

Biological life depends on motion, and this manifests itself in proteins that display motion over a formidable range of time scales spanning from femtoseconds vibrations of atoms at enzymatic transition states, all the way to slow domain motions occurring on micro to milliseconds. An outstanding challenge in contemporary biophysics and structural biology is a quantitative understanding of the linkages among protein structure, dynamics, and function. These linkages are becoming increasingly explorable due to conceptual and methodological advances. In this Perspective article, we will point toward future directions of the field of protein dynamics with an emphasis on enzymes. Research questions in the field are becoming increasingly complex such as the mechanistic understanding of high-order interaction networks in allosteric signal propagation through a protein matrix, or the connection between local and collective motions. In analogy to the solution to the "protein folding problem," we argue that the way forward to understanding these and other important questions lies in the successful integration of experiment and computation, while utilizing the present rapid expansion of sequence and structure space. Looking forward, the future is bright, and we are in a period where we are on the doorstep to, at least in part, comprehend the importance of dynamics for biological function.

Differential Regulation of the BDNF Gene in Cortical and Hippocampal Neurons

Brain-derived neurotrophic factor (BDNF) is a widely expressed neurotrophin that supports the survival, differentiation, and signaling of various neuronal populations. Although it has been well described that expression of BDNF is strongly regulated by neuronal activity, little is known whether regulation of BDNF expression is similar in different brain regions. Here, we focused on this fundamental question using neuronal populations obtained from rat cerebral cortices and hippocampi of both sexes. First, we thoroughly characterized the role of the best-described regulators of BDNF gene - cAMP response element binding protein (CREB) family transcription factors, and show that activity-dependent BDNF expression depends more on CREB and the coactivators CREB binding protein (CBP) and CREB-regulated transcriptional coactivator 1 (CRTC1) in cortical than in hippocampal neurons. Our data also reveal an important role of CREB in the early induction of BDNF mRNA expression after neuronal activity and only modest contribution after prolonged neuronal activity. We further corroborated our findings at BDNF protein level. To determine the transcription factors regulating BDNF expression in these rat brain regions in addition to CREB family, we used in vitro DNA pulldown assay coupled with mass spectrometry, chromatin immunoprecipitation (ChIP), and bioinformatics, and propose a number of neurodevelopmentally important transcription factors, such as FOXP1, SATB2, RAI1, BCL11A, and TCF4 as brain region-specific regulators of BDNF expression. Together, our data reveal complicated brain region-specific fine-tuning of BDNF expression.SIGNIFICANCE STATEMENT To date, majority of the research has focused on the regulation of brain-derived neurotrophic factor (BDNF) in the brain but much less is known whether the regulation of BDNF expression is universal in different brain regions and neuronal populations. Here, we report that the best described regulators of BDNF gene from the cAMP-response element binding protein (CREB) transcription factor family have a more profound role in the activity-dependent regulation of BDNF in cortex than in hippocampus. Our results indicate a brain region-specific fine tuning of BDNF expression. Moreover, we have used unbiased determination of novel regulators of the BDNF gene and report a number of neurodevelopmentally important transcription factors as novel potential regulators of the BDNF expression.

CBP-HSF2 structural and functional interplay in Rubinstein-Taybi neurodevelopmental disorder

Patients carrying autosomal dominant mutations in the histone/lysine acetyl transferases CBP or EP300 develop a neurodevelopmental disorder: Rubinstein-Taybi syndrome (RSTS). The biological pathways underlying these neurodevelopmental defects remain elusive. Here, we unravel the contribution of a stress-responsive pathway to RSTS. We characterize the structural and functional interaction between CBP/EP300 and heat-shock factor 2 (HSF2), a tuner of brain cortical development and major player in prenatal stress responses in the neocortex: CBP/EP300 acetylates HSF2, leading to the stabilization of the HSF2 protein. Consequently, RSTS patient-derived primary cells show decreased levels of HSF2 and HSF2-dependent alteration in their repertoire of molecular chaperones and stress response. Moreover, we unravel a CBP/EP300-HSF2-N-cadherin cascade that is also active in neurodevelopmental contexts, and show that its deregulation disturbs neuroepithelial integrity in 2D and 3D organoid models of cerebral development, generated from RSTS patient-derived iPSC cells, providing a molecular reading key for this complex pathology.

Intrinsically Disordered Proteins: An Overview

Many proteins and protein segments cannot attain a single stable three-dimensional structure under physiological conditions; instead, they adopt multiple interconverting conformational states. Such intrinsically disordered proteins or protein segments are highly abundant across proteomes, and are involved in various effector functions. This review focuses on different aspects of disordered proteins and disordered protein regions, which form the basis of the so-called "Disorder-function paradigm" of proteins. Additionally, various experimental approaches and computational tools used for characterizing disordered regions in proteins are discussed. Finally, the role of disordered proteins in diseases and their utility as potential drug targets are explored.

Intrinsic protein disorder uncouples affinity from binding specificity

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins often function by molecular recognition, in which they undergo induced folding. Based on prior generalizations, the idea prevails in the IDP field that due to the entropic penalty of induced folding, the major functional advantage associated with this binding mode is "uncoupling" specificity from binding strength. Nevertheless, both weaker binding and high specificity of IDPs/IDRs rest on limited experimental observations, making these assumptions more speculations than evidence-supported facts. The issue is also complicated by the rather vague concept of specificity that lacks an exact measure, such as the Kd for binding strength. We addressed these issues by creating and analyzing a comprehensive dataset of well-characterized ID/globular protein complexes, for which both the atomic structure of the complex and free energy (ΔG, Kd ) of interaction is known. Through this analysis, we provide evidence that the affinity distributions of IDP/globular and globular/globular complexes show different trends, whereas specificity does not connote to weaker binding strength of IDPs/IDRs. Furthermore, protein disorder extends the spectrum in the direction of very weak interactions, which may have important regulatory consequences and suggest that, in a biological sense, strict correlation of specificity and binding strength are uncoupled by structural disorder.

Molecular and functional characterization of an evolutionarily conserved CREB-binding protein in the Lymnaea CNS

In eukaryotes, CREB-binding protein (CBP), a coactivator of CREB, functions both as a platform for recruiting other components of the transcriptional machinery and as a histone acetyltransferase (HAT) that alters chromatin structure. We previously showed that the transcriptional activity of cAMP-responsive element binding protein (CREB) plays a crucial role in neuronal plasticity in the pond snail Lymnaea stagnalis. However, there is no information on the molecular structure and HAT activity of CBP in the Lymnaea central nervous system (CNS), hindering an investigation of its postulated role in long-term memory (LTM). Here, we characterize the Lymnaea CBP (LymCBP) gene and identify a conserved domain of LymCBP as a functional HAT. Like CBPs of other species, LymCBP possesses functional domains, such as the KIX domain, which is essential for interaction with CREB and was shown to regulate LTM. In-situ hybridization showed that the staining patterns of LymCBP mRNA in CNS are very similar to those of Lymnaea CREB1. A particularly strong LymCBP mRNA signal was observed in the cerebral giant cell (CGC), an identified extrinsic modulatory interneuron of the feeding circuit, the key to both appetitive and aversive LTM for taste. Biochemical experiments using the recombinant protein of the LymCBP HAT domain showed that its enzymatic activity was blocked by classical HAT inhibitors. Preincubation of the CNS with such inhibitors blocked cAMP-induced synaptic facilitation between the CGC and an identified follower motoneuron of the feeding system. Taken together, our findings suggest a role for the HAT activity of LymCBP in synaptic plasticity in the feeding circuitry.

Prediction of protein-protein interaction sites in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) participate in many biological processes by interacting with other proteins, including the regulation of transcription, translation, and the cell cycle. With the increasing amount of disorder sequence data available, it is thus crucial to identify the IDP binding sites for functional annotation of these proteins. Over the decades, many computational approaches have been developed to predict protein-protein binding sites of IDP (IDP-PPIS) based on protein sequence information. Moreover, there are new IDP-PPIS predictors developed every year with the rapid development of artificial intelligence. It is thus necessary to provide an up-to-date overview of these methods in this field. In this paper, we collected 30 representative predictors published recently and summarized the databases, features and algorithms. We described the procedure how the features were generated based on public data and used for the prediction of IDP-PPIS, along with the methods to generate the feature representations. All the predictors were divided into three categories: scoring functions, machine learning-based prediction, and consensus approaches. For each category, we described the details of algorithms and their performances. Hopefully, our manuscript will not only provide a full picture of the status quo of IDP binding prediction, but also a guide for selecting different methods. More importantly, it will shed light on the inspirations for future development trends and principles.

Toward Accurate Coarse-Grained Simulations of Disordered Proteins and Their Dynamic Interactions

Intrinsically disordered proteins (IDPs) play crucial roles in cellular regulatory networks and are now recognized to often remain highly dynamic even in specific interactions and assemblies. Accurate description of these dynamic interactions is extremely challenging using atomistic simulations because of the prohibitive computational cost. Efficient coarse-grained approaches could offer an effective solution to overcome this bottleneck if they could provide an accurate description of key local and global properties of IDPs in both unbound and bound states. The recently developed hybrid-resolution (HyRes) protein model has been shown to be capable of providing a semiquantitative description of the secondary structure propensities of IDPs. Here, we show that greatly improved description of global structures and transient interactions can be achieved by introducing a solvent-accessible surface area-based implicit solvent term followed by reoptimization of effective interaction strengths. The new model, termed HyRes II, can semiquantitatively reproduce a wide range of local and global structural properties of a set of IDPs of various lengths and complexities. It can also distinguish the level of compaction between folded proteins and IDPs. In particular, applied to the disordered N-terminal transactivation domain (TAD) of tumor suppressor p53, HyRes II is able to recapitulate various nontrivial structural properties compared to experimental results, some of them to a level of accuracy that is almost comparable to results from atomistic explicit solvent simulations. Furthermore, we demonstrate that HyRes II can be used to simulate the dynamic interactions of TAD with the DNA-binding domain of p53, generating structural ensembles that are highly consistent with existing NMR data. We anticipate that HyRes II will provide an efficient and relatively reliable tool toward accurate coarse-grained simulations of dynamic protein interactions.

The Disordered Amino Terminus of the Circadian Enzyme Nocturnin Modulates Its NADP(H) Phosphatase Activity by Changing Protein Dynamics

Endogenous circadian clocks control the rhythmicity of a broad range of behavioral and physiological processes, and this is entrained by the daily fluctuations in light and dark. Nocturnin (Noct) is a rhythmically expressed gene regulated by the circadian clock that belongs to the CCR4 family of endonuclease-exonuclease-phosphatase (EEP) enzymes, and the NOCT protein exhibits phosphatase activity, catalyzing the removal of the 2'-phosphate from NADP(H). In addition to its daily nighttime peak of expression, it is also induced by acute stimuli. Loss of Nocturnin (Noct-/-) in mice results in resistance to high-fat diet-induced obesity, and loss of Noct in HEK293T cells confers a protective effect to oxidative stress. Modeling of the full-length Nocturnin protein reveals a partially structured amino terminus that is disparate from its CCR4 family members. The high sequence conservation of a leucine zipper-like (LZ-like) motif, the only structural element in the amino terminus, highlights the potential importance of this domain in modulating phosphatase activity. In vitro biochemical and biophysical techniques demonstrate that the LZ-like domain within the flexible N-terminus is necessary for preserving the active site cleft in an optimal conformation to promote the efficient turnover of the substrate. This modulation occurs in cis and is pivotal in maintaining the stability and conformational integrity of the enzyme. These new findings suggest an additional layer of modulating the activity of Nocturnin in addition to its rhythmicity to provide fine-tuned control over cellular levels of NADPH.

DUX4 is a multifunctional factor priming human embryonic genome activation

Double homeobox 4 (DUX4) is expressed at the early pre-implantation stage in human embryos. Here we show that induced human DUX4 expression substantially alters the chromatin accessibility of non-coding DNA and activates thousands of newly identified transcribed enhancer-like regions, preferentially located within ERVL-MaLR repeat elements. CRISPR activation of transcribed enhancers by C-terminal DUX4 motifs results in the increased expression of target embryonic genome activation (EGA) genes ZSCAN4 and KHDC1P1. We show that DUX4 is markedly enriched in human zygotes, followed by intense nuclear DUX4 localization preceding and coinciding with minor EGA. DUX4 knockdown in human zygotes led to changes in the EGA transcriptome but did not terminate the embryos. We also show that the DUX4 protein interacts with the Mediator complex via the C-terminal KIX binding motif. Our findings contribute to the understanding of DUX4 as a regulator of the non-coding genome.

Gab2 promotes acute myeloid leukemia growth and migration through the SHP2-Erk-CREB signaling pathway

Acute myeloid leukemia (AML) is a hematologic malignant disease largely affecting older adults with poor outcomes. Lack of effective targeted treatment is a major challenge in managing the disease in the clinic. Scaffolding adaptor Gab2 is amplified in a subset of AML. However, the causative role of Gab2 in AML remains to be explored. In this study, it was found that Gab2 was expressed at high levels in AML patient samples and AML cell lines. Experiments by knocking down Gab2 expression using shRNA showed that Gab2 promoted AML cell growth and migration in vitro and in vivo. Further studies using Gab2 mutants and pharmacological inhibitors revealed that Gab2 increased CREB phosphorylation via the SHP-2/Erk signaling pathway. CREB phosphorylation contributed to Gab2-induced cell migration by increasing MMP2 and MMP9 expression. This research indicates that high Gab2 expression promotes AML progression through the SHP2-Erk-CREB signaling pathway. CREB suppression may help treat AML with high Gab2 expression.

Design, Synthesis and Biological Evaluation of Prodrugs of 666-15 as Inhibitors of CREB-Mediated Gene Transcription

cAMP-response element binding protein (CREB) is a transcription factor involved in multiple cancers. Chemical inhibitors of CREB represent potential anticancer agents. We previously identified 666-15 as a potent CREB inhibitor. While 666-15 showed efficacious anticancer activity in vivo through intraperitoneal (IP) injection, its oral bioavailability is limited. To increase its oral bioavailability, we describe synthesis and evaluation of prodrugs based on 666-15. The amino acid esters were attempted, but they were not stable for detailed characterization. The corresponding sulfate and phosphates were prepared. The sulfate of 666-15 was too stable to release 666-15 while the phosphates were converted into 666-15 with half-lives of ∼2 h. Phosphate 3 was also a potent CREB inhibitor with anti-breast cancer activity. Furthermore, compound 3 showed much improved oral bioavailability at 38%. These studies support that 3 can be used as an oral CREB inhibitor while IP administration of 666-15 is preferred for in vivo applications.

Cryo-EM, Protein Engineering, and Simulation Enable the Development of Peptide Therapeutics against Acute Myeloid Leukemia

Cryogenic electron microscopy (cryo-EM) has emerged as a viable structural tool for molecular therapeutics development against human diseases. However, it remains a challenge to determine structures of proteins that are flexible and smaller than 30 kDa. The 11 kDa KIX domain of CREB-binding protein (CBP), a potential therapeutic target for acute myeloid leukemia and other cancers, is a protein which has defied structure-based inhibitor design. Here, we develop an experimental approach to overcome the size limitation by engineering a protein double-shell to sandwich the KIX domain between apoferritin as the inner shell and maltose-binding protein as the outer shell. To assist homogeneous orientations of the target, disulfide bonds are introduced at the target-apoferritin interface, resulting in a cryo-EM structure at 2.6 Å resolution. We used molecular dynamics simulations to design peptides that block the interaction of the KIX domain of CBP with the intrinsically disordered pKID domain of CREB. The double-shell design allows for fluorescence polarization assays confirming the binding between the KIX domain in the double-shell and these interacting peptides. Further cryo-EM analysis reveals a helix-helix interaction between a single KIX helix and the best peptide, providing a possible strategy for developments of next-generation inhibitors.

Rational design of a helical peptide inhibitor targeting c-Myb–KIX interaction

The transcription factor c-Myb promotes the proliferation of hematopoietic cells by interacting with the KIX domain of CREB-binding protein; however, its aberrant expression causes leukemia. Therefore, inhibitors of the c-Myb–KIX interaction are potentially useful as antitumor drugs. Since the intrinsically disordered transactivation domain (TAD) of c-Myb binds KIX via a conformational selection mechanism where helix formation precedes binding, stabilizing the helical structure of c-Myb TAD is expected to increase the KIX-binding affinity. Here, to develop an inhibitor of the c-Myb–KIX interaction, we designed mutants of the c-Myb TAD peptide fragment where the helical structure is stabilized, based on theoretical predictions using AGADIR. Three of the four initially designed peptides each had a different Lys-to-Arg substitution on the helix surface opposite the KIX-binding interface. Furthermore, the triple mutant with three Lys-to-Arg substitutions, named RRR, showed a high helical propensity and achieved three-fold higher affinity to KIX than the wild-type TAD with a dissociation constant of 80 nM. Moreover, the RRR inhibitor efficiently competed out the c-Myb–KIX interaction. These results suggest that stabilizing the helical structure based on theoretical predictions, especially by conservative Lys-to-Arg substitutions, is a simple and useful strategy for designing helical peptide inhibitors of protein–protein interactions.

A propolis-derived small molecule ameliorates metabolic syndrome in obese mice by targeting the CREB/CRTC2 transcriptional complex

The molecular targets and mechanisms of propolis ameliorating metabolic syndrome are not fully understood. Here, we report that Brazilian green propolis reduces fasting blood glucose levels in obese mice by disrupting the formation of CREB/CRTC2 transcriptional complex, a key regulator of hepatic gluconeogenesis. Using a mammalian two-hybrid system based on CREB-CRTC2, we identify artepillin C (APC) from propolis as an inhibitor of CREB-CRTC2 interaction. Without apparent toxicity, APC protects mice from high fat diet-induced obesity, decreases fasting glucose levels, enhances insulin sensitivity and reduces lipid levels in the serum and liver by suppressing CREB/CRTC2-mediated both gluconeogenic and SREBP transcriptions. To develop more potential drugs from APC, we designed and found a novel compound, A57 that exhibits higher inhibitory activity on CREB-CRTC2 association and better capability of improving insulin sensitivity in obese animals, as compared with APC. In this work, our results indicate that CREB/CRTC2 is a suitable target for developing anti-metabolic syndrome drugs.

Metadynamics Simulations to Study the Structural Ensembles and Binding Processes of Intrinsically Disordered Proteins

The structures of intrinsically disordered proteins (IDPs) are highly dynamic. It is hard to characterize the structures of these proteins experimentally. Molecular dynamics (MD) simulation is a powerful tool in the understanding of protein dynamic structures and function. This chapter describes the application of metadynamics-based enhanced sampling methods in the study of phosphorylation regulation on the structure of kinase-inducible domains (KID). The structural properties of free pKID and KID were obtained by parallel tempering metadynamics combined with well-tempered ensemble (PTMetaD WTE) method, and the binding free energy surfaces of pKID/KID and KIX were characterized by bias-exchanged metadynamics (BE-MetaD) simulations.

Molecular Simulations of Intrinsically Disordered Proteins and Their Binding Mechanisms

Intrinsically disordered proteins (IDPs) lack well-defined secondary or tertiary structures in solution but are found to be involved in a wide range of critical cellular processes that highlight their functional importance. IDPs usually undergo folding upon binding to their targets. Such binding coupled to folding behavior has widened our perspective on the protein structure-dynamics-function paradigm in molecular biology. However, characterizing the folding upon binding mechanism of IDPs experimentally remains quite challenging. Molecular simulations emerge as a potentially powerful tool that offers information complementary to experiments. Here we present a general computational framework for the molecular simulations of IDP folding upon binding processes that combines all-atom molecular dynamics (MD) and coarse-grained simulations. The classical all-atom molecular dynamics approach using GPU acceleration allows the researcher to explore the properties of the IDP conformational ensemble, whereas coarse-grained structure-based models implemented with parameters carefully calibrated to available experimental measurements can be used to simulate the entire folding upon binding process. We also discuss a set of tools for the analysis of MD trajectories and describe the details of the computational protocol to follow so that it can be adapted by the user to study any IDP in isolation and in complex with partners.

The USR domain of USF1 mediates NF-Y interactions and cooperative DNA binding

The trimeric CCAAT-binding NF-Y is a "pioneer" Transcription Factor -TF- known to cooperate with neighboring TFs to regulate gene expression. Genome-wide analyses detected a precise stereo-alignment -10/12 bp- of CCAAT with E-box elements and corresponding colocalization of NF-Y with basic-Helix-Loop-Helix (bHLH) TFs. We dissected here NF-Y interactions with USF1 and MAX. USF1, but not MAX, cooperates in DNA binding with NF-Y. NF-Y and USF1 synergize to activate target promoters. Reconstruction of complexes by structural means shows independent DNA binding of MAX, whereas USF1 has extended contacts with NF-Y, involving the USR, a USF-specific amino acid sequence stretch required for trans-activation. The USR is an intrinsically disordered domain and adopts different conformations based on E-box-CCAAT distances. Deletion of the USR abolishes cooperative DNA binding with NF-Y. Our data indicate that the functionality of certain unstructured domains involves adapting to small variation in stereo-alignments of the multimeric TFs sites.

Structural and thermodynamical insights into the binding and inhibition of FIH-1 by the N-terminal disordered region of Mint3

Mint3 is known to enhance aerobic ATP production, known as the Warburg effect, by binding to FIH-1. Since this effect is considered to be beneficial for cancer cells, the interaction is a promising target for cancer therapy. However, previous research has suggested that the interacting region of Mint3 with FIH-1 is intrinsically disordered, which makes investigation of this interaction challenging. Therefore, we adopted thermodynamic and structural studies in solution to clarify the structural and thermodynamical changes of Mint3 binding to FIH-1. First, using a combination of circular dichroism, nuclear magnetic resonance, and hydrogen/deuterium exchange–mass spectrometry (HDX-MS), we confirmed that the N-terminal half, which is the interacting part of Mint3, is mostly disordered. Next, we revealed a large enthalpy and entropy change in the interaction of Mint3 using isothermal titration calorimetry (ITC). The profile is consistent with the model that the flexibility of disordered Mint3 is drastically reduced upon binding to FIH-1. Moreover, we performed a series of ITC experiments with several types of truncated Mint3s, an effective approach since the interacting part of Mint3 is disordered, and identified amino acids 78 to 88 as a novel core site for binding to FIH-1. The truncation study of Mint3 also revealed the thermodynamic contribution of each part of Mint3 to the interaction with FIH-1, where the core sites contribute to the affinity (ΔG), while other sites only affect enthalpy (ΔH), by forming noncovalent bonds. This insight can serve as a foothold for further investigation of intrinsically disordered regions (IDRs) and drug development for cancer therapy.

Diosgenin exerts anti-tumor effects through inactivation of cAMP/PKA/CREB signaling pathway in colorectal cancer

Colorectal cancer (CRC) is the most fatal gastrointestinal tumor and it is urge to explore powerful drugs for the treatment. Diosgenin (DSG) as a new steroidal had been reported exerts anti-tumor activity in multiple cancers, including CRC. However, the potential mechanism of DSG suppresses CRC remains further to be revealed. Here, we reported that DSG inhibited proliferation of CRC cells in dose- and time-dependent manner, induced apoptosis by modulating p53 and Bcl-2 family proteins expression to mediate mitochondrial apoptosis pathway, suppressed migration and invasion by reducing MMP-9 (matrix metalloproteinase) and decreased aerobic glycolysis by mediating glucose transporter (GLUT) like GLUT3 and GLUT4, and pyruvate carboxylase PC downregulation. Intriguingly, mechanistic study suggests those phenotypes involved DSG inhibited cAMP/PKA/CREB pathway in CRC cells, and result to inhibit the phosphorylation of CREB to regulate the transcription of genes above-mentioned. Finally, nude mice xenograft tumor model further indicated that DSG could be a great agent to suppress the growth of CRC cells in vivo and have no obvious side effects. Taken together, we revealed a unique mechanism that DSG suppresses CRC cells through cAMP/PKA/CREB pathway and DSG is a promising candidate drug for CRC treatment.

Target-binding behavior of IDPs via pre-structured motifs

Pre-Structured Motifs (PreSMos) are transient secondary structures observed in many intrinsically disordered proteins (IDPs) and serve as protein target-binding hot spots. The prefix "pre" highlights that PreSMos exist a priori in the target-unbound state of IDPs as the active pockets of globular proteins pre-exist before target binding. Therefore, a PreSMo is an "active site" of an IDP; it is not a spatial pocket, but rather a secondary structural motif. The classical and perhaps the most effective approach to understand the function of a protein has been to determine and investigate its structure. Ironically or by definition IDPs do not possess structure (here structure refers to tertiary structure only). Are IDPs then entirely structureless? The PreSMos provide us with an atomic-resolution answer to this question. For target binding, IDPs do not rely on the spatial pockets afforded by tertiary or higher structures. Instead, they utilize the PreSMos possessing particular conformations that highly presage the target-bound conformations. PreSMos are recognized or captured by targets via conformational selection (CS) before their conformations eventually become stabilized via structural induction into more ordered bound structures. Using PreSMos, a number of, if not all, IDPs can bind targets following a sequential pathway of CS followed by an induced fit (IF). This chapter presents several important PreSMos implicated in cancers, neurodegenerative diseases, and other diseases along with discussions on their conformational details that mediate target binding, a structural rationale for unstructured proteins.

A Dual-Site Inhibitor of CBP/p300 KIX is a Selective and Effective Modulator of Myb

The protein-protein interaction between the KIX motif of the transcriptional coactivator CBP/p300 and the transcriptional activator Myb is a high-value target due to its established role in certain acute myeloid leukemias (AML) and potential contributions to other cancers. However, the CBP/p300 KIX domain has multiple binding sites, several structural homologues, many binding partners, and substantial conformational plasticity, making it challenging to specifically target using small-molecule inhibitors. Here, we report a picomolar dual-site inhibitor (MybLL-tide) of the Myb-CBP/p300 KIX interaction. MybLL-tide has higher affinity for CBP/p300 KIX than any previously reported compounds while also possessing 5600-fold selectivity for the CBP/p300 KIX domain over other coactivator domains. MybLL-tide blocks the association of CBP and p300 with Myb in the context of the proteome, leading to inhibition of key Myb·KIX-dependent genes in AML cells. These results show that MybLL-tide is an effective, modifiable tool to selectively target the KIX domain and assess transcriptional effects in AML cells and potentially other cancers featuring aberrant Myb behavior. Additionally, the dual-site design has applicability to the other challenging coactivators that bear multiple binding surfaces.

Energy landscapes of protein aggregation and conformation switching in intrinsically disordered proteins

The protein folding problem was apparently solved recently by the advent of a deep learning method for protein structure prediction called AlphaFold. However, this program is not able to make predictions about the protein folding pathways. Moreover, it only treats about half of the human proteome, as the remaining proteins are intrinsically disordered or contain disordered regions. By definition these proteins differ from natively folded proteins and do not adopt a properly folded structure in solution. However these intrinsically disordered proteins (IDPs) also systematically differ in amino acid composition and uniquely often become folded upon binding to an interaction partner. These factors preclude solving IDP structures by current machine-learning methods like AlphaFold, which also cannot solve the protein aggregation problem, since this meta-folding process can give rise to different aggregate sizes and structures. An alternative computational method is provided by molecular dynamics simulations that already successfully explored the energy landscapes of IDP conformational switching and protein aggregation in multiple cases. These energy landscapes are very different from those of 'simple' protein folding, where one energy funnel leads to a unique protein structure. Instead, the energy landscapes of IDP conformational switching and protein aggregation feature a number of minima for different competing low-energy structures. In this review, I discuss the characteristics of these multifunneled energy landscapes in detail, illustrated by molecular dynamics simulations that elucidated the underlying conformational transitions and aggregation processes.

Roles of constitutive and signal-dependent protein phosphatase 2A docking motifs in burst attenuation of the cyclic AMP response element-binding protein

The cAMP response element-binding protein (CREB) is an important regulator of cell growth, metabolism, and synaptic plasticity. CREB is activated through phosphorylation of an evolutionarily conserved Ser residue (S133) within its intrinsically disordered kinase-inducible domain (KID). Phosphorylation of S133 in response to cAMP, Ca²⁺, and other stimuli triggers an association of the KID with the KID-interacting (KIX) domain of the CREB-binding protein (CBP), a histone acetyl transferase (HAT) that promotes transcriptional activation. Here we addressed the mechanisms of CREB attenuation following bursts in CREB phosphorylation. We show that phosphorylation of S133 is reversed by protein phosphatase 2A (PP2A), which is recruited to CREB through its B56 regulatory subunits. We found that a B56-binding site located at the carboxyl-terminal boundary of the KID (BS2) mediates high-affinity B56 binding, while a second binding site (BS1) located near the amino terminus of the KID mediates low affinity binding enhanced by phosphorylation of adjacent casein kinase (CK) phosphosites. Mutations that diminished B56 binding to BS2 elevated both basal and stimulus-induced phosphorylation of S133, increased CBP interaction with CREB, and potentiated the expression of CREB-dependent reporter genes. Cells from mice harboring a homozygous Creb^E153D mutation that disrupts BS2 exhibited increased S133 phosphorylation stoichiometry and elevated transcriptional bursts to cAMP. These findings provide insights into substrate targeting by PP2A holoenzymes and establish a new mechanism of CREB attenuation that has implications for understanding CREB signaling in cell growth, metabolism, synaptic plasticity, and other physiologic contexts.