Interplay between intrinsically disordered proteins inside membraneless protein liquid droplets

Intrinsically disordered proteins (IDPs) in cells phase separate to form diverse membraneless organelles, which have condensed liquid droplet-like properties and often contain multiple IDPs. However, how potential interactions between different IDPs affect the dynamic behavior of these protein droplets is largely unknown. Here, we develop a rapid IDP clustering system to generate protein droplets with varied residue compositions and examine diverse interacting IDPs inside droplets. Three different IDP droplets actively recruited other diverse IDPs inside droplets with extremely varied enrichment (inside/outside) degrees (over 100-fold variation) under highly crowded conditions. The recruited IDPs were mostly mobile even inside highly immobile droplets. Among the five tested IDPs, the disordered region of Ddx4 helicase with its unique multiple charged residue blocks was noticeably influenced by droplet mobility. We also discovered that droplets of different IDPs could rapidly fuse to each other. Interestingly, some droplets were heterogeneously fused with segregated subcompartments, and this segregation was enhanced by droplet maturation and was more apparent for specific IDP pairs, in which the polar and charged residue compositions are highly different. The present study not only reports multiple peculiar behaviors of interacting IDP pairs inside droplets but also provides valuable information on generating membraneless organelle models with controllable droplet properties.

Membraneless droplets of intrinsically disordered proteins (IDPs) with varied residue compositions uniquely interact with each other as droplets and clients.     

The importance of the compact disordered state in the fuzzy interactions between intrinsically disordered proteins

The intrinsically disordered C-terminal domain (CTD) of protein 4.1G is able to specifically bind a 26-residue intrinsically disordered region of NuMA, forming a dynamic fuzzy complex. As one of a few cases of extremely fuzzy interactions between two intrinsically disordered proteins/regions (IDPs/IDRs) without induced folding, the principle of the binding is unknown. Here, we combined experimental and computational methods to explore the detailed mechanism of the interaction between 4.1G-CTD and NuMA. MD simulations suggest that the kinetic hub states in the structure ensemble of 4.1G-CTD are favorable in the fuzzy complex. The feature of these hub states is that the binding ‘hot spot’ motifs βA and βB exhibit β strand propensities and are well packed to each other. The binding between 4.1G-CTD and NuMA is disrupted at low pH, which changes the intramolecular packing of 4.1G-CTD and weakens the packing between βA and βB motifs. Low pH conditions also lead to increased hydrodynamic radius and acceleration of backbone dynamics of 4.1G-CTD. All these results underscore the importance of tertiary structural arrangements and overall compactness of 4.1G-CTD in its binding to NuMA, i.e. the compact disordered state of 4.1G-CTD is crucial for binding. Different from the short linear motifs (SLiMs) that are often found to mediate IDP interactions, 4.1G-CTD functions as an intrinsically disordered domain (IDD), which is a functional and structural unit similar to conventional protein domains. This work sheds light on the molecular recognition mechanism of IDPs/IDRs and expands the conventional structure-function paradigm in protein biochemistry.

Tertiary structural arrangements and overall compactness is important for interactions between IDPs.     

How well does molecular simulation reproduce environment-specific conformations of the intrinsically disordered peptides PLP,TP2 and ONEG?

Understanding the conformational ensembles of intrinsically disordered proteins and peptides (IDPs) in their various biological environments is essential for understanding their mechanisms and functional roles in the proteome, leading to a greater knowledge of, and potential treatments for, a broad range of diseases. To determine whether molecular simulation is able to generate accurate conformational ensembles of IDPs, we explore the structural landscape of the PLP peptide (an intrinsically disordered region of the proteolipid membrane protein) in aqueous and membrane-mimicking solvents, using replica exchange with solute scaling (REST2), and examine the ability of four force fields (ff14SB, ff14IDPSFF, CHARMM36 and CHARMM36m) to reproduce literature circular dichroism (CD) data. Results from variable temperature (VT) ¹H and Rotating frame Overhauser Effect SpectroscopY (ROESY) nuclear magnetic resonance (NMR) experiments are also presented and are consistent with the structural observations obtained from the simulations and CD. We also apply the optimum simulation protocol to TP2 and ONEG (a cell-penetrating peptide (CPP) and a negative control peptide, respectively) to gain insight into the structural differences that may account for the observed difference in their membrane-penetrating abilities. Of the tested force fields, we find that CHARMM36 and CHARMM36m are best suited to the study of IDPs, and accurately predict a disordered to helical conformational transition of the PLP peptide accompanying the change from aqueous to membrane-mimicking solvents. We also identify an α-helical structure of TP2 in the membrane-mimicking solvents and provide a discussion of the mechanistic implications of this observation with reference to the previous literature on the peptide. From these results, we recommend the use of CHARMM36m with the REST2 protocol for the study of environment-specific IDP conformations. We believe that the simulation protocol will allow the study of a broad range of IDPs that undergo conformational transitions in different biological environments.

A protocol for simulating intrinsically disordered peptides in aqueous and hydrophobic solvents is proposed. Results from four force fields are compared with experiment. CHARMM36m performs the best for the simulated IDPs in all environments.     

Computational strategy for intrinsically disordered protein ligand design leads to the discovery of p53 transactivation domain I binding compounds that activate the p53 pathway

Intrinsically disordered proteins or intrinsically disordered regions (IDPs) have gained much attention in recent years due to their vital roles in biology and prevalence in various human diseases. Although IDPs are perceived as attractive therapeutic targets, rational drug design targeting IDPs remains challenging because of their conformational heterogeneity. Here, we propose a hierarchical computational strategy for IDP drug virtual screening (IDPDVS) and applied it in the discovery of p53 transactivation domain I (TAD1) binding compounds. IDPDVS starts from conformation sampling of the IDP target, then it combines stepwise conformational clustering with druggability evaluation to identify potential ligand binding pockets, followed by multiple docking screening runs and selection of compounds that can bind multi-conformations. p53 is an important tumor suppressor and restoration of its function provides an opportunity to inhibit cancer cell growth. TAD1 locates at the N-terminus of p53 and plays key roles in regulating p53 function. No compounds that directly bind to TAD1 have been reported due to its highly disordered structure. We successfully used IDPDVS to identify two compounds that bind p53 TAD1 and restore wild-type p53 function in cancer cells. Our study demonstrates that IDPDVS is an efficient strategy for IDP drug discovery and p53 TAD1 can be directly targeted by small molecules.

A hierarchical computational strategy for IDP drug virtual screening (IDPDVS) was proposed and successfully applied to identify compounds that bind p53 TAD1 and restore wild-type p53 function in cancer cells.     

Time Window Expansion for HDX Analysis of an Intrinsically Disordered Protein

Application of typical HDX methods to examine intrinsically disordered proteins (IDP), proteins that are natively unstructured and highly dynamic at physiological pH, is limited because of the rapid exchange of unprotected amide hydrogens with solvent. The exchange rates of these fast exchanging amides are usually faster than the shortest time scale (10 s) employed in typical automated HDX-MS experiments. Considering the functional importance of IDPs and their association with many diseases, it is valuable to develop methods that allow the study of solution dynamics of these proteins as well as the ability to probe the interaction of IDPs with their wide range of binding partners. Here, we report the application of time window expansion to the millisecond range by altering the on-exchange pH of the HDX experiment to study a well-characterized IDP; the activation domain of the nuclear receptor coactivator, peroxisome proliferator-activated receptor gamma coactivator-1 alpha (PGC-1α). This method enabled mapping the regions of PGC-1α that are stabilized upon binding the ligand binding domain (LBD) of the nuclear receptor peroxisome proliferator-activated receptor gamma (PPARγ). We further demonstrate the method’s applicability to other binding partners of the IDP PGC-1α and pave the way for characterizing many other biologically important ID proteins.

Energetics of a protein disorder–order transition in small molecule recognition

Many proteins recognise other proteins via mechanisms that involve the folding of intrinsically disordered regions upon complex formation. Here we investigate how the selectivity of a drug-like small molecule arises from its modulation of a protein disorder-to-order transition. Binding of the compound AM-7209 has been reported to confer order upon an intrinsically disordered ‘lid’ region of the oncoprotein MDM2. Calorimetric measurements revealed that truncation of the lid region of MDM2 increases the apparent dissociation constant of AM-7209 250-fold. By contrast, lid truncation has little effect on the binding of the ligand Nutlin-3a. Insights into these differential binding energetics were obtained via a complete thermodynamic analysis that featured adaptive absolute alchemical free energy of binding calculations with enhanced-sampling molecular dynamics simulations. The simulations reveal that in apo MDM2 the ordered lid state is energetically disfavoured. AM-7209, but not Nutlin-3a, shows a significant energetic preference for ordered lid conformations, thus shifting the balance towards ordering of the lid in the AM-7209/MDM2 complex. The methodology reported herein should facilitate broader targeting of intrinsically disordered regions in medicinal chemistry.

Molecular simulations and biophysical measurements elucidate why the ligand AM-7209 orders a disordered region of the protein MDM2 on binding. This work expands strategies available to medicinal chemists for targeting disordered proteins.     

Polyampholyte physics: Liquid–liquid phase separation and biological condensates

Here we review the current understanding of molecular interactions that govern liquid–liquid phase separation (LLPS) of biological condensates. The connection between sequence, chain conformation, and phase separation of intrinsically disordered proteins (IDPs) and their model polyampholytes is discussed. In particular, we highlight how the charge pattern influences the conformation and phase behavior of natural IDPs. We then describe recent results from theoretical treatments of polyampholytes implementing random phase approximation, field-theoretic simulations, and transfer matrix theory that show an increase in charge segregation results in an increased tendency to phase separate.     

Triazolinedione protein modification: from an overlooked off-target effect to a tryptophan-based bioconjugation strategy

Labelling of tyrosine residues in peptides and proteins has been reported to selectively occur via a ‘tyrosine-click’ reaction with triazolinedione reagents (TAD). However, we here demonstrate that TAD reagents are actually not selective for tyrosine and that tryptophan residues are in fact also labelled with these reagents. This off-target labelling remained under the radar as it is challenging to detect these physiologically stable but thermally labile modifications with the commonly used HCD and CID MS/MS techniques. We show that selectivity of tryptophan over tyrosine can be achieved by lowering the pH of the aqueous buffer to effect selective Trp-labelling. Given the low relative abundance of tryptophan compared to tyrosine in natural proteins, this results in a new site-selective bioconjugation method that does not rely on enzymes nor unnatural amino acids and is demonstrated for peptides and recombinant proteins.

A new strategy for selective tryptophan modification using triazolinedione (TAD) chemistry at pH 4 is shown on peptides and proteins. Additionally, off-target modification of tryptophan residues during the classical TAD-Y click reaction is uncovered.     

The Dynamic Multisite Interactions between Two Intrinsically Disordered Proteins

Protein interactions involving intrinsically disordered proteins (IDPs) comprise a variety of binding modes, from the well‐characterized folding upon binding to dynamic fuzzy complexes. To date, most studies concern the binding of an IDP to a structured protein, while the interaction between two IDPs is poorly understood. In this study, NMR, smFRET, and molecular dynamics (MD) simulation are combined to characterize the interaction between two IDPs, the C‐terminal domain (CTD) of protein 4.1G and the nuclear mitotic apparatus (NuMA) protein. It is revealed that CTD and NuMA form a fuzzy complex with remaining structural disorder. Multiple binding sites on both proteins were identified by molecular dynamics and mutagenesis studies. This study provides an atomic scenario in which two IDPs bearing multiple binding sites interact with each other in dynamic equilibrium. The combined approach employed here could be widely applicable for investigating IDPs and their dynamic interactions.     

Role of backbone-solvent interactions in determining conformational equilibria of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are functional proteins that do not fold into well-defined three-dimensional structures under physiological conditions. IDP sequences have low hydrophobicity, and hence, recent experiments have focused on quantitative studies of conformational ensembles of archetypal IDP sequences such as polyglutamine and glycine-serine block copolypeptides. Results from these experiments show that, despite the absence of hydrophobic residues, polar IDPs prefer ensembles of collapsed structures in aqueous milieus. Do these preferences originate in interactions that are unique to polar sidechains? The current study addresses this issue by analyzing conformational equilibria for polyglycine and a glycine-serine block copolypeptide in two environments, namely, water and 8 M urea. Polyglycine, a poly secondary-amide, has no sidechains and is a useful model system for generic polypeptide backbones. Results based on large-scale molecular dynamics simulations show that polyglycine forms compact, albeit disordered, globules in water and swollen, disordered coils in 8 M urea. There is minimal overlap between conformational ensembles in the two environments. Analysis of order parameters derived from theories for flexible polymers show that water at ambient temperatures is a poor solvent for generic polypeptide backbones. Therefore, the experimentally observed preferences for polyglutamine and glycine-serine block copolypeptides must originate, at least partially, in polypeptide backbones. A preliminary analysis of the driving forces that lead to distinct conformational preferences for polyglycine in two different environments is presented. Implications for describing conformational ensembles of generic IDP sequences are also discussed.     

Rate constants and mechanisms of intrinsically disordered proteins binding to structured targets

The binding of intrinsically disordered proteins (IDPs) to structured targets is gaining increasing attention. Here we review experimental and computational studies on the binding kinetics of IDPs. Experiments have yielded both the binding rate constants and the binding mechanisms, the latter via mutation and deletion studies and NMR techniques. Most computational studies have aimed at qualitative understanding of the binding rate constants or at mapping the free energy surfaces after the IDPs are engaged with their targets. The experiments and computation show that IDPs generally gain structures after they are engaged with their targets; that is, interactions with the targets facilitate the IDPs' folding. It also seems clear that the initial contact of an IDP with the target is formed by just a segment, not the entire IDP. The docking of one segment to its sub-site followed by coalescing of other segments around the corresponding sub-sites emerges as a recurring feature in the binding of IDPs. Such a dock-and-coalesce model forms the basis for quantitative calculation of binding rate constants. For both disordered and ordered proteins, strong electrostatic attraction with their targets can enhance the binding rate constants by several orders of magnitude. There are now tremendous opportunities in narrowing the gap in our understanding of IDPs relative to ordered proteins with regard to binding kinetics.     

Increasing the Chemical‐Shift Dispersion of Unstructured Proteins with a Covalent Lanthanide Shift Reagent

The study of intrinsically disordered proteins (IDPs) by NMR often suffers from highly overlapped resonances that prevent unambiguous chemical‐shift assignments, and data analysis that relies on well‐separated resonances. We present a covalent paramagnetic lanthanide‐binding tag (LBT) for increasing the chemical‐shift dispersion and facilitating the chemical‐shift assignment of challenging, repeat‐containing IDPs. Linkage of the DOTA‐based LBT to a cysteine residue induces pseudo‐contact shifts (PCS) for resonances more than 20 residues from the spin‐labeling site. This leads to increased chemical‐shift dispersion and decreased signal overlap, thereby greatly facilitating chemical‐shift assignment. This approach is applicable to IDPs of varying sizes and complexity, and is particularly helpful for repeat‐containing IDPs and low‐complexity regions. This results in improved efficiency for IDP analysis and binding studies.     

Does liquid–liquid phase separation drive peptide folding?

Proline–arginine (PR) dipeptide repeats have been shown to undergo liquid–liquid phase separation and are an example of a growing number of intrinsically disordered proteins that can assemble into membraneless organelles. These structures have been posited as nucleation sites for pathogenic protein aggregation. As such, a better understanding of the effects that the increased local concentration and volumetric crowding within droplets have on peptide secondary structure is necessary. Herein we use Fourier transform infrared (FTIR) and two-dimensional infrared (2DIR) spectroscopy to show that formation of droplets by PR₂₀ accompanies changes in the amide-I spectra consistent with folding into poly-proline helical structures.

Two-dimensional infrared spectroscopy reveals folding of an intrinsically disordered peptide when sequestered into a model “membrane-less” organelle.     

Compensatory Adaptations of Structural Dynamics in an Intrinsically Disordered Protein Complex

Intrinsically disordered proteins (IDPs) play crucial roles in protein interaction networks and in this context frequently constitute important hubs and interfaces. Here we show by a combination of NMR and EPR spectroscopy that the binding of the cytokine osteopontin (OPN) to its natural ligand, heparin, is accompanied by thermodynamically compensating structural adaptations. The core segment of OPN expands upon binding. This “unfolding‐upon‐binding” is governed primarily through electrostatic interactions between heparin and charged patches along the protein backbone and compensates for entropic penalties due to heparin–OPN binding. It is shown how structural unfolding compensates for entropic losses through ligand binding in IDPs and elucidates the interplay between structure and thermodynamics of rapid substrate‐binding and ‐release events in IDP interaction networks.     

Phase separation of p53 precedes aggregation and is affected by oncogenic mutations and ligands

Mutant p53 tends to form aggregates with amyloid properties, especially amyloid oligomers inside the nucleus, which are believed to cause oncogenic gain-of-function (GoF). The mechanism of the formation of the aggregates in the nucleus remains uncertain. The present study demonstrated that the DNA-binding domain of p53 (p53C) underwent phase separation (PS) on the pathway to aggregation under various conditions. p53C phase separated in the presence of the crowding agent polyethylene glycol (PEG). Similarly, mutant p53C (M237I and R249S) underwent PS; however, the process evolved to a solid-like phase transition faster than that in the case of wild-type p53C. The data obtained by microscopy of live cells indicated that transfection of mutant full-length p53 into the cells tended to result in PS and phase transition (PT) in the nuclear compartments, which are likely the cause of the GoF effects. Fluorescence recovery after photobleaching (FRAP) experiments revealed liquid characteristics of the condensates in the nucleus. Mutant p53 tended to undergo gel- and solid-like phase transitions in the nucleus and in nuclear bodies demonstrated by slow and incomplete recovery of fluorescence after photobleaching. Polyanions, such as heparin and RNA, were able to modulate PS and PT in vitro. Heparin apparently stabilized the condensates in a gel-like state, and RNA apparently induced a solid-like state of the protein even in the absence of PEG. Conditions that destabilize p53C into a molten globule conformation also produced liquid droplets in the absence of crowding. The disordered transactivation domain (TAD) modulated both phase separation and amyloid aggregation. In summary, our data provide mechanistic insight into the formation of p53 condensates and conditions that may result in the formation of aggregated structures, such as mutant amyloid oligomers, in cancer. The pathway of mutant p53 from liquid droplets to gel-like and solid-like (amyloid) species may be a suitable target for anticancer therapy.

Mutant p53 tends to form aggregates with amyloid properties, especially amyloid oligomers inside the nucleus, which are believed to cause oncogenic gain-of-function (GoF).     

Evidence for an Ordering Transition near 120 K in an Intrinsically Disordered Protein,Casein

Intrinsically disordered proteins (IDPs) are proteins that possess large unstructured regions. Their importance is increasingly recognized in biology but their characterization remains a challenging task. We employed field swept Electron Spin Echoes in pulsed EPR to investigate low-temperature stochastic molecular librations in a spin-labeled IDP, casein (the main protein of milk). For comparison, a spin-labeled globular protein, hen egg white lysozyme, is also investigated. For casein these motions were found to start at 100 K while for lysozyme only above 130 K, which was ascribed to a denser and more ordered molecular packing in lysozyme. However, above 120 K, the motions in casein were found to depend on temperature much slower than those in lysozyme. This abrupt change in casein was assigned to an ordering transition in which peptide residues rearrange making the molecular packing more rigid and/or more cohesive. The found features of molecular motions in these two proteins turned out to be very similar to those known for gel-phase lipid bilayers composed of conformationally ordered and conformationally disordered lipids. This analogy with a simpler molecular system may appear helpful for elucidation properties of molecular packing in IDPs.     

Proteus: a random forest classifier to predict disorder-to-order transitioning binding regions in intrinsically disordered proteins

The focus of the computational structural biology community has taken a dramatic shift over the past one-and-a-half decades from the classical protein structure prediction problem to the possible understanding of intrinsically disordered proteins (IDP) or proteins containing regions of disorder (IDPR). The current interest lies in the unraveling of a disorder-to-order transitioning code embedded in the amino acid sequences of IDPs/IDPRs. Disordered proteins are characterized by an enormous amount of structural plasticity which makes them promiscuous in binding to different partners, multi-functional in cellular activity and atypical in folding energy landscapes resembling partially folded molten globules. Also, their involvement in several deadly human diseases (e.g. cancer, cardiovascular and neurodegenerative diseases) makes them attractive drug targets, and important for a biochemical understanding of the disease(s). The study of the structural ensemble of IDPs is rather difficult, in particular for transient interactions. When bound to a structured partner, an IDPR adapts an ordered conformation in the complex. The residues that undergo this disorder-to-order transition are called protean residues, generally found in short contiguous stretches and the first step in understanding the modus operandi of an IDP/IDPR would be to predict these residues. There are a few available methods which predict these protean segments from their amino acid sequences; however, their performance reported in the literature leaves clear room for improvement. With this background, the current study presents ‘Proteus’, a random forest classifier that predicts the likelihood of a residue undergoing a disorder-to-order transition upon binding to a potential partner protein. The prediction is based on features that can be calculated using the amino acid sequence alone. Proteus compares favorably with existing methods predicting twice as many true positives as the second best method (55 vs. 27%) with a much higher precision on an independent data set. The current study also sheds some light on a possible ‘disorder-to-order’ transitioning consensus, untangled, yet embedded in the amino acid sequence of IDPs. Some guidelines have also been suggested for proceeding with a real-life structural modeling involving an IDPR using Proteus.