Non-cooperative 4E-BP2 folding with exchange between eIF4E-binding and binding-incompatible states tunes cap-dependent translation inhibition

Mapping the Multiscale Proteomic Organization of Cellular and Disease Phenotypes

While the primary sequences of human proteins have been cataloged for over a decade, determining how these are organized into a dynamic collection of multiprotein assemblies, with structures and functions spanning biological scales, is an ongoing venture. Systematic and data-driven analyses of these higher-order structures are emerging, facilitating the discovery and understanding of cellular phenotypes. At present, knowledge of protein localization and function has been primarily derived from manual annotation and curation in resources such as the Gene Ontology, which are biased toward richly annotated genes in the literature. Here, we envision a future powered by data-driven mapping of protein assemblies. These maps can capture and decode cellular functions through the integration of protein expression, localization, and interaction data across length scales and timescales. In this review, we focus on progress toward constructing integrated cell maps that accelerate the life sciences and translational research.

The human eIF4E:4E-BP2 complex structure for studying hyperphosphorylation

The cap-dependent mRNA translation is dysregulated in many kinds of cancers. The interaction between eIF4E and eIF4G through a canonical eIF4E-binding motif (CEBM) determines the efficacy of the cap-dependent mRNA translation. eIF4E-binding proteins (4E-BPs) share the CEBM and compete with eIF4G for the same binding surface of eIF4E and then inhibit the mRNA translation. 4E-BPs function as tumor repressors in nature. Hyperphosphorylation of 4E-BPs regulates the structure folding and causes the dissociation of 4E-BPs from eIF4E. However, until now, there has been no structure of the full-length 4E-BPs in complex with eIF4E. The regulation mechanism of phosphorylation is still unclear. In this work, we first investigate the interactions of human eIF4E with the CEBM and an auxiliary eIF4E-binding motif (AEBM) in eIF4G and 4E-BPs. The results unravel that the structure and interactions of the CEBM are highly conserved between eIF4G and 4E-BPs. However, the extended CEBM (ECEBM) in 4E-BPs forms a longer helix than that in eIF4G. The residue R62 in the ECEBM of 4E-BP2 forms salt bridges with E32 and E70 of eIF4E. The residue R63 of 4E-BP2 forms two special hydrogen bonds with N77 of eIF4E. Both of these interactions are missing in eIF4G. The AEBM of 4E-BPs folds into a β-sheet conformation, which protects V81 inside a hydrophobic core in 4E-BP2. In eIF4G, the AEBM exists in a random coil state. The hydrophilic residues S637 and D638 of eIF4G open the hydrophobic core for solvents. The results show that the ECEBM and AEBM may be responsible for the competing advantage of 4E-BP2. Finally, based on our previous work (J. Zeng, F. Jiang and Y. D. Wu, J. Chem. Theory Comput., 2017, 13, 320), the human eIF4E:4E-BP2 complex (eIF4E:BP2P18-I88) including all reported phosphorylation sites is predicted. The eIF4E:BP2P18-I88 complex is different from the existing experimental eIF4E:eIF4G complex and provides an important structure for further studying the regulation mechanism of phosphorylation in 4E-BPs.

Local Disordered Region Sampling (LDRS) for ensemble modeling of proteins with experimentally undetermined or low confidence prediction segments

SUMMARY

The Local Disordered Region Sampling (LDRS, pronounced loaders) tool is a new module developed for IDPConformerGenerator, a previously validated approach to model intrinsically disordered proteins (IDPs). The IDPConformerGenerator LDRS module provides a method for generating all-atom conformations of intrinsically disordered protein regions at N- and C-termini of and in loops or linkers between folded regions of an existing protein structure. These disordered elements often lead to missing coordinates in experimental structures or low confidence in predicted structures. Requiring only a pre-existing PDB or mmCIF formatted structural template of the protein with missing coordinates or with predicted confidence scores and its full-length primary sequence, LDRS will automatically generate physically meaningful conformational ensembles of the missing flexible regions to complete the full-length protein. The capabilities of the LDRS tool of IDPConformerGenerator include modeling phosphorylation sites using enhanced Monte Carlo-Side Chain Entropy, transmembrane proteins within an all-atom bilayer, and multi-chain complexes. The modeling capacity of LDRS capitalizes on the modularity, the ability to be used as a library and via command-line, and the computational speed of the IDPConformerGenerator platform.

AVAILABILITY AND IMPLEMENTATION

The LDRS module is part of the IDPConformerGenerator modeling suite, which can be downloaded from GitHub at https://github.com/julie-forman-kay-lab/IDPConformerGenerator. IDPConformerGenerator is written in Python3 and works on Linux, Microsoft Windows, and Mac OS versions that support DSSP. Users can utilize LDRS's Python API for scripting the same way they can use any part of IDPConformerGenerator's API, by importing functions from the "idpconfgen.ldrs_helper" library. Otherwise, LDRS can be used as a command line interface application within IDPConformerGenerator. Full documentation is available within the command-line interface as well as on IDPConformerGenerator's official documentation pages (https://idpconformergenerator.readthedocs.io/en/latest/).

Delineating Structural Propensities of the 4E-BP2 Protein via Integrative Modeling and Clustering

The intrinsically disordered 4E-BP2 protein regulates mRNA cap-dependent translation through interaction with the predominantly folded eukaryotic initiation factor 4E (eIF4E). Phosphorylation of 4E-BP2 dramatically reduces the level of eIF4E binding, in part by stabilizing a binding-incompatible folded domain. Here, we used a Rosetta-based sampling algorithm optimized for IDRs to generate initial ensembles for two phospho forms of 4E-BP2, non- and 5-fold phosphorylated (NP and 5P, respectively), with the 5P folded domain flanked by N- and C-terminal IDRs (N-IDR and C-IDR, respectively). We then applied an integrative Bayesian approach to obtain NP and 5P conformational ensembles that agree with experimental data from nuclear magnetic resonance, small-angle X-ray scattering, and single-molecule Förster resonance energy transfer (smFRET). For the NP state, inter-residue distance scaling and 2D maps revealed the role of charge segregation and pi interactions in driving contacts between distal regions of the chain (∼70 residues apart). The 5P ensemble shows prominent contacts of the N-IDR region with the two phosphosites in the folded domain, pT37 and pT46, and, to a lesser extent, delocalized interactions with the C-IDR region. Agglomerative hierarchical clustering led to partitioning of each of the two ensembles into four clusters with different global dimensions and contact maps. This helped delineate an NP cluster that, based on our smFRET data, is compatible with the eIF4E-bound state. 5P clusters were differentiated by interactions of C-IDR with the folded domain and of the N-IDR with the two phosphosites in the folded domain. Our study provides both a better visualization of fundamental structural poses of 4E-BP2 and a set of falsifiable insights on intrachain interactions that bias folding and binding of this protein.

Local Disordered Region Sampling (LDRS) for Ensemble Modeling of Proteins with Experimentally Undetermined or Low Confidence Prediction Segments

The Local Disordered Region Sampling (LDRS, pronounced loaders) tool, developed for the IDPConformerGenerator platform (Teixeira et al. 2022), provides a method for generating all-atom conformations of intrinsically disordered regions (IDRs) at N- and C-termini of and in loops or linkers between folded regions of an existing protein structure. These disordered elements often lead to missing coordinates in experimental structures or low confidence in predicted structures. Requiring only a pre-existing PDB structure of the protein with missing coordinates or with predicted confidence scores and its full-length primary sequence, LDRS will automatically generate physically meaningful conformational ensembles of the missing flexible regions to complete the full-length protein. The capabilities of the LDRS tool of IDPConformerGenerator include modeling phosphorylation sites using enhanced Monte Carlo Side Chain Entropy (MC-SCE) (Bhowmick and Head-Gordon 2015), transmembrane proteins within an all-atom bilayer, and multi-chain complexes. The modeling capacity of LDRS capitalizes on the modularity, ability to be used as a library and via command-line, and computational speed of the IDPConformerGenerator platform.

The Conformation of the Intrinsically Disordered N-Terminal Region of Barrier-to-Autointegration Factor (BAF) is Regulated by pH and Phosphorylation

Barrier-to-Autointegration Factor (BAF) is a highly conserved DNA binding protein important for genome integrity. Its localization and function are regulated through phosphorylation. Previously reported structures of BAF suggested that it is fully ordered, but our recent NMR analysis revealed that its N-terminal region is flexible in solution and that S4/T3 di-phosphorylation by VRK1 reduces this flexibility. Here, molecular dynamics (MD) simulation was used to unveil the conformational ensembles accessible to the N-terminal region of BAF either unphosphorylated, mono-phosphorylated on S4 or di-phosphorylated on S4/T3 (pBAF) and to reveal the interactions that contribute to define these ensembles. We show that the intrinsic flexibility observed in the N-terminal region of BAF is reduced by S4 phosphorylation and to a larger extent by S4/T3 di-phosphorylation. Thanks to the atomic description offered by MD supported by the NMR study of several BAF mutants, we identified the dynamic network of salt bridge interactions responsible for the conformational restriction involving pS4 and pT3 with residues located in helix α1 and α6. Using MD, we showed that the flexibility in the N-terminal region of BAF depends on the ionic strength and on the pH. We show that the presence of two negative charges of the phosphoryl groups is required for a substantial decrease in flexibility in pBAF. Using MD supported by NMR, we also showed that H7 deprotonation reduces the flexibility in the N-terminal region of BAF. Thus, the conformation of the intrinsically disordered N-terminal region of BAF is highly tunable, likely related to its diverse functions.

eIF4E-homologous protein (4EHP): a multifarious cap-binding protein

The cap-binding protein 4EHP/eIF4E2 has been a recent object of interest in the field of post-transcriptional gene regulation and translational control. From ribosome-associated quality control, to RNA decay and microRNA-mediated gene silencing, this member of the eIF4E protein family regulates gene expression through numerous pathways. Low in abundance but ubiquitously expressed, 4EHP interacts with different binding partners to form multiple protein complexes that regulate translation in a variety of biological contexts. Documented functions of 4EHP primarily relate to its role as a translational repressor, but recent findings indicate that it might also participate in the activation of translation in specific settings. In this review, we discuss the known functions, properties and mechanisms that involve 4EHP in the control of gene expression. We also discuss our current understanding of how 4EHP processes are regulated in eukaryotic cells, and the diseases implicated with dysregulation of 4EHP-mediated translational control.

Multisite phosphorylation and binding alter conformational dynamics of the 4E-BP2 protein

Intrinsically disordered proteins (IDPs) play critical roles in regulatory protein interactions, but detailed structural/dynamic characterization of their ensembles remain challenging, both in isolation and when they form dynamic "fuzzy" complexes. Such is the case for mRNA cap-dependent translation initiation, which is regulated by the interaction of the predominantly folded eukaryotic initiation factor 4E (eIF4E) with the intrinsically disordered eIF4E binding proteins (4E-BPs) in a phosphorylation-dependent manner. Single-molecule Förster resonance energy transfer showed that the conformational changes of 4E-BP2 induced by binding to eIF4E are non-uniform along the sequence; while a central region containing both motifs that bind to eIF4E expands and becomes stiffer, the C-terminal region is less affected. Fluorescence anisotropy decay revealed a non-uniform segmental flexibility around six different labeling sites along the chain. Dynamic quenching of these fluorescent probes by intrinsic aromatic residues measured via fluorescence correlation spectroscopy report on transient intra- and inter-molecular contacts on nanosecond-to-microsecond timescales. Upon hyperphosphorylation, which induces folding of ∼40 residues in 4E-BP2, the quenching rates decreased at most labeling sites. The chain dynamics around sites in the C-terminal region far away from the two binding motifs significantly increased upon binding to eIF4E, suggesting that this region is also involved in the highly dynamic 4E-BP2:eIF4E complex. Our time-resolved fluorescence data paint a sequence-level rigidity map of three states of 4E-BP2 differing in phosphorylation or binding status and distinguish regions that form contacts with eIF4E. This study adds complementary structural and dynamics information to recent studies of 4E-BP2, and it constitutes an important step toward a mechanistic understanding of this important IDP via integrative modeling.

A fine balance of hydrophobic-electrostatic communication pathways in a pH-switching protein

Allostery is the phenomenon of coupling between distal binding sites in a protein. Such coupling is at the crux of protein function and regulation in a myriad of scenarios, yet determining the molecular mechanisms of coupling networks in proteins remains a major challenge. Here, we report mechanisms governing pH-dependent myristoyl switching in monomeric hisactophilin, whereby the myristoyl moves between a sequestered state, i.e., buried within the core of the protein, to an accessible state, in which the myristoyl has increased accessibility for membrane binding. Measurements of the pH and temperature dependence of amide chemical shifts reveal protein local structural stability and conformational heterogeneity that accompany switching. An analysis of these measurements using a thermodynamic cycle framework shows that myristoyl-proton coupling at the single-residue level exists in a fine balance and extends throughout the protein. Strikingly, small changes in the stereochemistry or size of core and surface hydrophobic residues by point mutations readily break, restore, or tune myristoyl switch energetics. Synthesizing the experimental results with those of molecular dynamics simulations illuminates atomistic details of coupling throughout the protein, featuring a large network of hydrophobic interactions that work in concert with key electrostatic interactions. The simulations were critical for discerning which of the many ionizable residues in hisactophilin are important for switching and identifying the contributions of nonnative interactions in switching. The strategy of using temperature-dependent NMR presented here offers a powerful, widely applicable way to elucidate the molecular mechanisms of allostery in proteins at high resolution.

mTOR substrate phosphorylation in growth control

The target of rapamycin (TOR), discovered 30 years ago, is a highly conserved serine/threonine protein kinase that plays a central role in regulating cell growth and metabolism. It is activated by nutrients, growth factors, and cellular energy. TOR forms two structurally and functionally distinct complexes, TORC1 and TORC2. TOR signaling activates cell growth, defined as an increase in biomass, by stimulating anabolic metabolism while inhibiting catabolic processes. With emphasis on mammalian TOR (mTOR), we comprehensively reviewed the literature and identified all reported direct substrates. In the context of recent structural information, we discuss how mTORC1 and mTORC2, despite having a common catalytic subunit, phosphorylate distinct substrates. We conclude that the two complexes recruit different substrates to phosphorylate a common, minimal motif.

Control of the eIF4E activity: structural insights and pharmacological implications

The central role of eukaryotic translation initiation factor 4E (eIF4E) in controlling mRNA translation has been clearly assessed in the last decades. eIF4E function is essential for numerous physiological processes, such as protein synthesis, cellular growth and differentiation; dysregulation of its activity has been linked to ageing, cancer onset and progression and neurodevelopmental disorders, such as autism spectrum disorder (ASD) and Fragile X Syndrome (FXS). The interaction between eIF4E and the eukaryotic initiation factor 4G (eIF4G) is crucial for the assembly of the translational machinery, the initial step of mRNA translation. A well-characterized group of proteins, named 4E-binding proteins (4E-BPs), inhibits the eIF4E–eIF4G interaction by competing for the same binding site on the eIF4E surface. 4E-BPs and eIF4G share a single canonical motif for the interaction with a conserved hydrophobic patch of eIF4E. However, a second non-canonical and not conserved binding motif was recently detected for eIF4G and several 4E-BPs. Here, we review the structural features of the interaction between eIF4E and its molecular partners eIF4G and 4E-BPs, focusing on the implications of the recent structural and biochemical evidence for the development of new therapeutic strategies. The design of novel eIF4E-targeting molecules that inhibit translation might provide new avenues for the treatment of several conditions.

A Tale of Loops and Tails: The Role of Intrinsically Disordered Protein Regions in R-Loop Recognition and Phase Separation

R-loops are non-canonical, three-stranded nucleic acid structures composed of a DNA:RNA hybrid, a displaced single-stranded (ss)DNA, and a trailing ssRNA overhang. R-loops perform critical biological functions under both normal and disease conditions. To elucidate their cellular functions, we need to understand the mechanisms underlying R-loop formation, recognition, signaling, and resolution. Previous high-throughput screens identified multiple proteins that bind R-loops, with many of these proteins containing folded nucleic acid processing and binding domains that prevent (e.g., topoisomerases), resolve (e.g., helicases, nucleases), or recognize (e.g., KH, RRMs) R-loops. However, a significant number of these R-loop interacting Enzyme and Reader proteins also contain long stretches of intrinsically disordered regions (IDRs). The precise molecular and structural mechanisms by which the folded domains and IDRs synergize to recognize and process R-loops or modulate R-loop-mediated signaling have not been fully explored. While studying one such modular R-loop Reader, the Fragile X Protein (FMRP), we unexpectedly discovered that the C-terminal IDR (C-IDR) of FMRP is the predominant R-loop binding site, with the three N-terminal KH domains recognizing the trailing ssRNA overhang. Interestingly, the C-IDR of FMRP has recently been shown to undergo spontaneous Liquid-Liquid Phase Separation (LLPS) assembly by itself or in complex with another non-canonical nucleic acid structure, RNA G-quadruplex. Furthermore, we have recently shown that FMRP can suppress persistent R-loops that form during transcription, a process that is also enhanced by LLPS via the assembly of membraneless transcription factories. These exciting findings prompted us to explore the role of IDRs in R-loop processing and signaling proteins through a comprehensive bioinformatics and computational biology study. Here, we evaluated IDR prevalence, sequence composition and LLPS propensity for the known R-loop interactome. We observed that, like FMRP, the majority of the R-loop interactome, especially Readers, contains long IDRs that are highly enriched in low complexity sequences with biased amino acid composition, suggesting that these IDRs could directly interact with R-loops, rather than being “mere flexible linkers” connecting the “functional folded enzyme or binding domains”. Furthermore, our analysis shows that several proteins in the R-loop interactome are either predicted to or have been experimentally demonstrated to undergo LLPS or are known to be associated with phase separated membraneless organelles. Thus, our overall results present a thought-provoking hypothesis that IDRs in the R-loop interactome can provide a functional link between R-loop recognition via direct binding and downstream signaling through the assembly of LLPS-mediated membrane-less R-loop foci. The absence or dysregulation of the function of IDR-enriched R-loop interactors can potentially lead to severe genomic defects, such as the widespread R-loop-mediated DNA double strand breaks that we recently observed in Fragile X patient-derived cells.

Investigation of Phosphorylation-Induced Folding of an Intrinsically Disordered Protein by Coarse-Grained Molecular Dynamics

Apart from being the most common mechanism of regulating protein function and transmitting signals throughout the cell, phosphorylation has an ability to induce disorder-to-order transition in an intrinsically disordered protein. In particular, it was shown that folding of the intrinsically disordered protein, eIF4E-binding protein isoform 2 (4E-BP2), can be induced by multisite phosphorylation. Here, the principles that govern the folding of phosphorylated 4E-BP2 (pT37pT46 4E-BP2_18-62) are investigated by analyzing canonical and replica exchange molecular dynamics trajectories, generated with the coarse-grained united-residue force field, in terms of local and global motions and the time dependence of formation of contacts between C^αs of selected pairs of residues. The key residues involved in the folding of the pT37pT46 4E-BP2_18-62 are elucidated by this analysis. The correlations between local and global motions are identified. Moreover, for a better understanding of the physics of the formation of the folded state, the experimental structure of the pT37pT46 4E-BP2_18-62 is analyzed in terms of a kink (heteroclinic standing wave solution) of a generalized discrete nonlinear Schrödinger equation. It is shown that without molecular dynamics simulations the kinks are able to identify not only the phosphorylated sites of protein, the key players in folding, but also the reasons for the weak stability of the pT37pT46 4E-BP2_18-62.

The dynamic mechanism of 4E-BP1 recognition and phosphorylation by mTORC1

The activation of cap-dependent translation in eukaryotes requires multisite, hierarchical phosphorylation of 4E-BP by the 1 MDa kinase mammalian target of rapamycin complex 1 (mTORC1). To resolve the mechanism of this hierarchical phosphorylation at the atomic level, we monitored by NMR spectroscopy the interaction of intrinsically disordered 4E binding protein isoform 1 (4E-BP1) with the mTORC1 subunit regulatory-associated protein of mTOR (Raptor). The N-terminal RAIP motif and the C-terminal TOR signaling (TOS) motif of 4E-BP1 bind separate sites in Raptor, resulting in avidity-based tethering of 4E-BP1. This tethering orients the flexible central region of 4E-BP1 toward the mTORC1 kinase site for phosphorylation. The structural constraints imposed by the two tethering interactions, combined with phosphorylation-induced conformational switching of 4E-BP1, explain the hierarchy of 4E-BP1 phosphorylation by mTORC1. Furthermore, we demonstrate that mTORC1 recognizes both free and eIF4E-bound 4E-BP1, allowing rapid phosphorylation of the entire 4E-BP1 pool and efficient activation of translation. Finally, our findings provide a mechanistic explanation for the differential rapamycin sensitivity of the 4E-BP1 phosphorylation sites.

A biphenyl inhibitor of eIF4E targeting an internal binding site enables the design of cell-permeable PROTAC-degraders

The eukaryotic translation initiation factor 4E (eIF4E) is the master regulator of cap-dependent protein synthesis. Overexpression of eIF4E is implicated in diseases such as cancer, where dysregulation of oncogenic protein translation is frequently observed. eIF4E has been an attractive target for cancer treatment. Here we report a high-resolution X-ray crystal structure of eIF4E in complex with a novel inhibitor (i4EG-BiP) that targets an internal binding site, in contrast to the previously described inhibitor, 4EGI-1, which binds to the surface. We demonstrate that i4EG-BiP is able to displace the scaffold protein eIF4G and inhibit the proliferation of cancer cells. We provide insights into how i4EG-BiP is able to inhibit cap-dependent translation by increasing the eIF4E-4E-BP1 interaction while diminishing the interaction of eIF4E with eIF4G. Leveraging structural details, we designed proteolysis targeted chimeras (PROTACs) derived from 4EGI-1 and i4EG-BiP and characterized these on biochemical and cellular levels. We were able to design PROTACs capable of binding eIF4E and successfully engaging Cereblon, which targets proteins for proteolysis. However, these initial PROTACs did not successfully stimulate degradation of eIF4E, possibly due to competitive effects from 4E-BP1 binding. Our results highlight challenges of targeted proteasomal degradation of eIF4E that must be addressed by future efforts.