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.     

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.     

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.     

Determinants for intrinsically disordered protein recruitment into phase-separated protein condensates

Multivalent interactions between amino acid residues of intrinsically disordered proteins (IDPs) drive phase separation of these proteins into liquid condensates, forming various membrane-less organelles in cells. These interactions between often biased residues of IDPs are also likely involved in selective recruitment of many other IDPs into condensates. However, determining factors for this IDP recruitment into protein condensates are not understood yet. Here, we quantitatively examined recruitment tendencies of various IDPs with different sequence compositions into IDP-clustered condensates both in vitro as well as in cells. Condensate-forming IDP scaffolds, recruited IDP clients, and phase separation conditions were carefully varied to find key factors for selective IDP partitioning in protein condensates. Regardless of scaffold sequences, charged residues in client IDPs assured potent IDP recruitment, likely via strong electrostatic interactions, where positive residues could further enhance recruitment, possibly with cation–pi interactions. Notably, poly-ethylene glycol, a widely used crowding reagent for in vitro phase separation, abnormally increased IDP recruitment, indicating the need for careful use of crowding conditions. Tyrosines of IDP clients also strongly participated in recruitment both in vitro and in cells. Lastly, we measured recruitment degrees by more conventional interactions between folded proteins instead of disordered proteins. Surprisingly, recruitment forces by an even moderate protein interaction (K_d ∼ 5 μM) were substantially stronger than those by natural IDP–IDP interactions. The present data offer valuable information on how cells might organize protein partitioning on various protein condensates.

Diverse interactions between folded and disordered proteins collectively dictate selective protein recruitment into bimolecular condensates.     

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.     

Mechanism‐Dependent Modulation of Ultrafast Interfacial Water Dynamics in Intrinsically Disordered Protein Complexes

The recognition of intrinsically disordered proteins (IDPs) is highly dependent on dynamics owing to the lack of structure. Here we studied the interplay between dynamics and molecular recognition in IDPs with a combination of time‐resolving tools on timescales ranging from femtoseconds to nanoseconds. We interrogated conformational dynamics and surface water dynamics and its attenuation upon partner binding using two IDPs, IBB and Nup153FG, both of central relevance to the nucleocytoplasmic transport machinery. These proteins bind the same nuclear transport receptor (Importinβ) with drastically different binding mechanisms, coupled folding–binding and fuzzy complex formation, respectively. Solvent fluctuations in the dynamic interface of the Nup153FG‐Importinβ fuzzy complex were largely unperturbed and slightly accelerated relative to the unbound state. In the IBB‐Importinβ complex, on the other hand, substantial relative slowdown of water dynamics was seen in a more rigid interface. These results show a correlation between interfacial water dynamics and the plasticity of IDP complexes, implicating functional relevance for such differential modulation in cellular processes, including nuclear transport.     

How random are intrinsically disordered proteins? A small angle scattering perspective

While the crucial role of intrinsically disordered proteins (IDPs) in the cell cycle is now recognized, deciphering their molecular mode of action at the structural level still remains highly challenging and requires a combination of many biophysical approaches. Among them, small angle X-ray scattering (SAXS) has been extremely successful in the last decade and has become an indispensable technique for addressing many of the fundamental questions regarding the activities of IDPs. After introducing some experimental issues specific to IDPs and in relation to the latest technical developments, this article presents the interest of the theory of polymer physics to evaluate the flexibility of fully disordered proteins. The different strategies to obtain 3-dimensional models of IDPs, free in solution and associated in a complex, are then reviewed. Indeed, recent computational advances have made it possible to readily extract maximum information from the scattering curve with a special emphasis on highly flexible systems, such as multidomain proteins and IDPs. Furthermore, integrated computational approaches now enable the generation of ensembles of conformers to translate the unique flexible characteristics of IDPs by taking into consideration the constraints of more and more various complementary experiment. In particular, a combination of SAXS with high-resolution techniques, such as x-ray crystallography and NMR, allows us to provide reliable models and to gain unique structural insights about the protein over multiple structural scales. The latest neutron scattering experiments also promise new advances in the study of the conformational changes of macromolecules involving more complex systems.     

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.     

Mapping Multivalency and Differential Affinities within Large Intrinsically Disordered Protein Complexes with Segmental Motion Analysis

Intrinsically disordered proteins (IDPs) can bind to multiple interaction partners. Numerous binding regions in the IDP that act in concert through complex cooperative effects facilitate such interactions, but complicate studying IDP complexes. To address this challenge we developed a combined fluorescence correlation and time‐resolved polarization spectroscopy approach to study the binding properties of the IDP nucleoporin153 (Nup153) to nuclear transport receptors (NTRs). The detection of segmental backbone mobility of Nup153 within the unperturbed complex provided a readout of local, region‐specific binding properties that are usually masked in measurements of the whole IDP. The binding affinities of functionally and structurally diverse NTRs to distinct regions of Nup153 can differ by orders of magnitudes—a result with implications for the diversity of transport routes in nucleocytoplasmic transport.     

Analytical Description of NMR Relaxation Highlights Correlated Dynamics in Intrinsically Disordered Proteins

The dynamic fluctuations of intrinsically disordered proteins (IDPs) define their function. Although experimental nuclear magnetic resonance (NMR) relaxation reveals the motional complexity of these highly flexible proteins, the absence of physical models describing IDP dynamics hinders their mechanistic interpretation. Combining molecular dynamics simulation and NMR, we introduce a framework in which distinct motions are attributed to local libration, backbone dihedral angle dynamics and longer‐range tumbling of one or more peptide planes. This model provides unique insight into segmental organization of dynamics in IDPs and allows us to investigate the presence and extent of the correlated motions that are essential for function.     

Synthesis and Conformations of Dendronized Poly(L-lysine)

Dendronized polymers based on a poly(L-lysine) backbone have been synthesized up to the fourth generation. The hydroxyl-terminated polymers are completely water-soluble, which makes them good candidates for drug delivery applications. The dendronized polypeptide backbones are helical at lower generations, but undergo a dramatic conformational change from alpha-helical to disordered upon increasing the dendron size to the third generation. This conformational change, attributed to steric repulsions between dendrons, is supported by spectroscopic measurements while chain extension upon dendronization is confirmed by scanning force microscopy.     

Structure-Function Insights of Jaburetox and Soyuretox: Novel Intrinsically Disordered Polypeptides Derived from Plant Ureases

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) do not have a stable 3D structure but still have important biological activities. Jaburetox is a recombinant peptide derived from the jack bean (Canavalia ensiformis) urease and presents entomotoxic and antimicrobial actions. The structure of Jaburetox was elucidated using nuclear magnetic resonance which reveals it is an IDP with small amounts of secondary structure. Different approaches have demonstrated that Jaburetox acquires certain folding upon interaction with lipid membranes, a characteristic commonly found in other IDPs and usually important for their biological functions. Soyuretox, a recombinant peptide derived from the soybean (Glycine max) ubiquitous urease and homologous to Jaburetox, was also characterized for its biological activities and structural properties. Soyuretox is also an IDP, presenting more secondary structure in comparison with Jaburetox and similar entomotoxic and fungitoxic effects. Moreover, Soyuretox was found to be nontoxic to zebra fish, while Jaburetox was innocuous to mice and rats. This profile of toxicity affecting detrimental species without damaging mammals or the environment qualified them to be used in biotechnological applications. Both peptides were employed to develop transgenic crops and these plants were active against insects and nematodes, unveiling their immense potentiality for field applications.     

Molecular biosensing system based on intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) that undergo structural transition upon binding their target molecules are becoming increasingly known. IDPs, because of their binding specificity and induced folding properties, can serve as biological recognition elements for sensing applications. In this paper, BRCA1, an IDP, was utilized as the biological recognition element to detect tumor suppressor protein p53 through the BRCA1/p53 binding interaction to serve as a proof-of-concept for the use of IDPs as recognition elements. The binding resulted in a disordered-to-ordered BRCA1 conformational change, as seen in our circular dichroism (CD) measurements. This conformational change in BRCA1 (residues 219-498) was utilized in the detection of p53 (residues 311-393) via both intrinsic and extrinsic fluorescent probes. Intrinsic tryptophan residues within the BRCA1 sequence detected p53 (311-393) with a detection limit of 0.559 nM (0.112 pmol). Two environmentally sensitive fluorophores, tetramethylrhodamine-5-maleimide (TMR) and 6-((5-dimethylaminonaphthalene-1-sulfonyl)amino)hexanoic acid, succinimidyl ester (dansyl-X, SE) were conjugated to BRCA1 (219-498). Dansyl-X, SE-conjugated BRCA1 (219-498) detected p53 (311-393) with a detection limit of 1.50 nM (0.300 pmol). The sensitivities for TMR and dansyl-X, SE-conjugated BRCA1 for the detection of p53 were nearly threefold and twofold higher, respectively, than the sensitivity reported using intrinsic BRCA1 tryptophan fluorescence. CD measurements did not reveal a disruption of p53/dye-conjugated BRCA1 binding, thus validating the applicability of environmentally sensitive fluorophores as transduction moieties to detect molecules which bind to IDPs and induce a structural change.     

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.     

Disorder-to-order transition of an intrinsically disordered region of sortase revealed by multiscale enhanced sampling

Molecular functions of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs), such as molecular recognition and cellular signaling, are ascribed to dynamic changes in the conformational space in response to binding of target molecules. Sortase, a transpeptitase in Gram-positive bacteria, has an IDR in a loop which undergoes a disordered-to-ordered transition (called "disordered loop"), accompanying a tilt of another loop ("dynamic loop"), upon binding of a signal peptide and a calcium ion. In this study, all-atom conformational ensembles of sortase were calculated for the four different binding states (with/without the peptide and with/without a calcium ion) by the multiscale enhanced sampling (MSES) simulation to examine how the binding of the peptide and/or calcium influences the conformational ensemble. The MSES is a multiscale and multicopy simulation method that allows an enhanced sampling of the all-atom model of large proteins including explicit solvent. A 100 ns MSES simulation of the ligand-free sortase using 20 replicas (in total 2 μs) demonstrated large flexibility in both the disordered and dynamic loops; however, their distributions were not random but had a clear preference which populates the N-terminal part of the disordered loop near the bound form. The MSES simulations of the three binding states clarified the allosteric mechanism of sortase: the N- and C-terminal parts of the disordered loop undergo a disorder-to-order transition independently of each other upon binding of the peptide and a calcium ion, respectively; however, upon binding of both ligands, the two parts work cooperatively to stabilize the bound peptide.     

Two‐Dimensional Crowding Uncovers a Hidden Conformation of α‐Synuclein

The intrinsically disordered protein (IDP), α‐synuclein (αS), is well‐known for phospholipid membrane binding‐coupled folding into tunable helical conformers. Here, using single‐molecule experiments in conjunction with ensemble assays and a theoretical model, we present a unique case demonstrating that the interaction–folding landscape of αS can be tuned by two‐dimensional (2D) crowding through simultaneous binding of a second protein on the bilayer surface. Unexpectedly, the experimental data show a clear deviation from a simple competitive inhibition model, but are consistent with a bimodal inhibition mechanism wherein membrane binding of a second protein (a membrane interacting chaperone, Hsp27, in this case) differentially inhibits two distinct modules of αS–membrane interaction. As a consequence, αS molecules are forced to access a hidden conformational state on the phospholipid bilayer in which only the higher‐affinity module remains membrane‐bound. Our results demonstrate that macromolecular crowding in two dimensions can play a significant role in shaping the conformational landscape of membrane‐binding IDPs with multiple binding modes.     

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.     

Multiplexed Discrimination of Single Amino Acid Residues in Polypeptides in a Single SERS Hot Spot

The SERS-based detection of protein sequences with single-residue sensitivity suffers from signal dominance of aromatic amino acid residues and backbones, impeding detection of non-aromatic amino acid residues. Herein, we trap a gold nanoparticle in a plasmonic nanohole to generate a single SERS hot spot for single-molecule detection of 2 similar polypeptides (vasopressin and oxytocin) and 10 distinct amino acids that constitute the 2 polypeptides. Significantly, both aromatic and non-aromatic amino acids are detected and discriminated at the single-molecule level either at individual amino acid molecules or within the polypeptide chains. Correlated with molecular dynamics simulations, our results suggest that the signal dominance due to large spatial occupancy of aromatic rings of the polypeptide sidechains on gold surfaces can be overcome by the high localization of the single hot spot. The superior spectral and spatial discriminative power of our approach can be applied to single-protein analysis, fingerprinting, and sequencing.     

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.