Structure/function implications in a dynamic complex of the intrinsically disordered Sic1 with the Cdc4 subunit of an SCF ubiquitin ligase

Reference Data Set for Circular Dichroism Spectroscopy Comprised of Validated Intrinsically Disordered Protein Models

Circular dichroism (CD) spectroscopy is an analytical technique that measures the wavelength-dependent differential absorbance of circularly polarized light and is applicable to most biologically important macromolecules, such as proteins, nucleic acids, and carbohydrates. It serves to characterize the secondary structure composition of proteins, including intrinsically disordered proteins, by analyzing their recorded spectra. Several computational tools have been developed to interpret protein CD spectra. These methods have been calibrated and tested mostly on globular proteins with well-defined structures, mainly due to the lack of reliable reference structures for disordered proteins. It is therefore still largely unclear how accurately these computational methods can determine the secondary structure composition of disordered proteins. Here, we provide such a required reference data set consisting of model structural ensembles and matching CD spectra for eight intrinsically disordered proteins. Using this set of data, we have assessed the accuracy of several published CD prediction and secondary structure estimation tools, including our own CD analysis package, SESCA. Our results show that for most of the tested methods, their accuracy for disordered proteins is generally lower than for globular proteins. In contrast, SESCA, which was developed using globular reference proteins, but was designed to be applicable to disordered proteins as well, performs similarly well for both classes of proteins. The new reference data set for disordered proteins should allow for further improvement of all published methods.

Protein ensemble modeling and analysis with MMMx

Proteins, especially of eukaryotes, often have disordered domains and may contain multiple folded domains whose relative spatial arrangement is distributed. The MMMx ensemble modeling and analysis toolbox (https://github.com/gjeschke/MMMx) can support the design of experiments to characterize the distributed structure of such proteins, starting from AlphaFold2 predictions or folded domain structures. Weak order can be analyzed with reference to a random coil model or to peptide chains that match the residue-specific Ramachandran angle distribution of the loop regions and are otherwise unrestrained. The deviation of the mean square end-to-end distance of chain sections from their average over sections of the same sequence length reveals localized compaction or expansion of the chain. The shape sampled by disordered chains is visualized by superposition in the principal axes frame of their inertia tensor. Ensembles of different sizes and with weighted conformers can be compared based on a similarity parameter that abstracts from the ensemble width.

Multisite phosphorylation dictates selective E2-E3 pairing as revealed by Ubc8/UBE2H-GID/CTLH assemblies

Ubiquitylation is catalyzed by coordinated actions of E3 and E2 enzymes. Molecular principles governing many important E3-E2 partnerships remain unknown, including those for RING-family GID/CTLH E3 ubiquitin ligases and their dedicated E2, Ubc8/UBE2H (yeast/human nomenclature). GID/CTLH-Ubc8/UBE2H-mediated ubiquitylation regulates biological processes ranging from yeast metabolic signaling to human development. Here, cryoelectron microscopy (cryo-EM), biochemistry, and cell biology reveal this exquisitely specific E3-E2 pairing through an unconventional catalytic assembly and auxiliary interactions 70-100 Å away, mediated by E2 multisite phosphorylation. Rather than dynamic polyelectrostatic interactions reported for other ubiquitylation complexes, multiple Ubc8/UBE2H phosphorylation sites within acidic CK2-targeted sequences specifically anchor the E2 C termini to E3 basic patches. Positions of phospho-dependent interactions relative to the catalytic domains correlate across evolution. Overall, our data show that phosphorylation-dependent multivalency establishes a specific E3-E2 partnership, is antagonistic with dephosphorylation, rigidifies the catalytic centers within a flexing GID E3-substrate assembly, and facilitates substrate collision with ubiquitylation active sites.

Functional unfoldomics: Roles of intrinsic disorder in protein (multi)functionality

Intrinsically disordered proteins (IDPs), which are functional proteins without stable tertiary structure, and hybrid proteins containing ordered domains and intrinsically disordered regions (IDRs) constitute prominent parts of all proteomes collectively known as unfoldomes. IDPs/IDRs exist as highly dynamic structural ensembles of rapidly interconverting conformations and are characterized by the exceptional structural heterogeneity, where their different parts are (dis)ordered to different degree, and their overall structure represents a complex mosaic of foldons, inducible foldons, inducible morphing foldons, non-foldons, semifoldons, and even unfoldons. Despite their lack of unique 3D structures, IDPs/IDRs play crucial roles in the control of various biological processes and the regulation of different cellular pathways and are commonly involved in recognition and signaling, indicating that the disorder-based functional repertoire is complementary to the functions of ordered proteins. Furthermore, IDPs/IDRs are frequently multifunctional, and this multifunctionality is defined by their structural flexibility and heterogeneity. Intrinsic disorder phenomenon is at the roots of the structure-function continuum model, where the structure continuum is defined by the presence of differently (dis)ordered regions, and the function continuum arises from the ability of all these differently (dis)ordered parts to have different functions. In their everyday life, IDPs/IDRs utilize a broad spectrum of interaction mechanisms thereby acting as interaction specialists. They are crucial for the biogenesis of numerous proteinaceous membrane-less organelles driven by the liquid-liquid phase separation. This review introduces functional unfoldomics by representing some aspects of the intrinsic disorder-based functionality.

Machine Learning Subtle Conformational Change due to Phosphorylation in Intrinsically Disordered Proteins

Phosphorylation of intrinsically disordered proteins/regions (IDPs/IDRs) has a profound effect in biological functions such as cell signaling, protein folding or unfolding, and long-range allosteric effects. However, here we focus on two IDPs, namely 83-residue IDR transcription factor Ash1 and 92-residue long N-terminal region of CDK inhibitor Sic1 protein, found in Saccharomyces cerevisiae, for which experimental measurements of average conformational properties, namely, radius of gyration and structure factor, indicate negligible changes upon phosphorylation. Here, we show that a judicious dissection of conformational ensemble via combination of unsupervised machine learning and extensive molecular dynamics (MD) trajectories can highlight key differences and similarities among the phosphorylated and wild-type IDP. In particular, we develop Markov state model (MSM) using the latent-space dimensions of an autoencoder, trained using multi-microsecond long MD simulation trajectories. Examination of structural changes among the states, prior to and upon phosphorylation, captured several similarities and differences in their backbone contact maps, secondary structure, and torsion angles. Hydrogen bonding analysis revealed that phosphorylation not only increases the number of hydrogen bonds but also switches the pattern of hydrogen bonding between the backbone and side chain atoms with the phosphorylated residues. We also observe that although phosphorylation introduces salt bridges, there is a loss of the cation-π interaction. Phosphorylation also improved the probability for long-range hydrophobic contacts and also enhanced interaction with water molecules and improved the local structure of water as evident from the geometric order parameters. The observations on these machine-learnt states gave important insights, as it would otherwise be difficult to determine experimentally which is important, if we were to understand the role of phosphorylation of IDPs in their biological functions.

The C-terminal tail of CSNAP attenuates the CSN complex

Protein degradation is one of the essential mechanisms that enables reshaping of the proteome landscape in response to various stimuli. The largest E3 ubiquitin ligase family that targets proteins to degradation by catalyzing ubiquitination is the cullin-RING ligases (CRLs). Many of the proteins that are regulated by CRLs are central to tumorigenesis and tumor progression, and dysregulation of the CRL family is frequently associated with cancer. The CRL family comprises ∼300 complexes, all of which are regulated by the COP9 signalosome complex (CSN). Therefore, CSN is considered an attractive target for therapeutic intervention. Research efforts for targeted CSN inhibition have been directed towards inhibition of the complex enzymatic subunit, CSN5. Here, we have taken a fresh approach focusing on CSNAP, the smallest CSN subunit. Our results show that the C-terminal region of CSNAP is tightly packed within the CSN complex, in a groove formed by CSN3 and CSN8. We show that a 16 amino acid C-terminal peptide, derived from this CSN-interacting region, can displace the endogenous CSNAP subunit from the complex. This, in turn, leads to a CSNAP null phenotype that attenuates CSN activity and consequently CRLs function. Overall, our findings emphasize the potential of a CSNAP-based peptide for CSN inhibition as a new therapeutic avenue.

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.

Balanced Three-Point Water Model OPC3-B for Intrinsically Disordered and Ordered Proteins

Intrinsically disordered proteins (IDPs) play a critical role in many biological processes. Due to the inherent structural flexibility of IDPs, experimental methods present significant challenges for sampling their conformational information at the atomic level. Therefore, molecular dynamics (MD) simulations have emerged as the primary tools for modeling IDPs whose accuracy depend on force field and water model. To enhance the accuracy of physical modeling of IDPs, several force fields have been developed. However, current water models lack precision and underestimate the interaction between water molecules and proteins. Here, we used Monte-Carlo re-weighting method to re-parameterize a three-point water model based on OPC3 for IDPs (named OPC3-B). We benchmarked the performance of OPC3-B compared with nine different water models for 10 IDPs and three ordered proteins. The results indicate that the performance of OPC3-B is better than other water models for both IDPs and ordered proteins. At the same time, OPC3-B possess the power of transferability with other force field to simulate IDPs. This newly developed water model can be used to insight into the research of sequence-disordered-function paradigm for IDPs.

Machine Learning Generation of Dynamic Protein Conformational Ensembles

Machine learning has achieved remarkable success across a broad range of scientific and engineering disciplines, particularly its use for predicting native protein structures from sequence information alone. However, biomolecules are inherently dynamic, and there is a pressing need for accurate predictions of dynamic structural ensembles across multiple functional levels. These problems range from the relatively well-defined task of predicting conformational dynamics around the native state of a protein, which traditional molecular dynamics (MD) simulations are particularly adept at handling, to generating large-scale conformational transitions connecting distinct functional states of structured proteins or numerous marginally stable states within the dynamic ensembles of intrinsically disordered proteins. Machine learning has been increasingly applied to learn low-dimensional representations of protein conformational spaces, which can then be used to drive additional MD sampling or directly generate novel conformations. These methods promise to greatly reduce the computational cost of generating dynamic protein ensembles, compared to traditional MD simulations. In this review, we examine recent progress in machine learning approaches towards generative modeling of dynamic protein ensembles and emphasize the crucial importance of integrating advances in machine learning, structural data, and physical principles to achieve these ambitious goals.

Simple and Effective Conformational Sampling Strategy for Intrinsically Disordered Proteins Using the UNRES Web Server

Intrinsically disordered proteins (IDPs) contain more charged amino acids than folded proteins, resulting in a lack of hydrophobic core(s) and a tendency to adopt rapidly interconverting structures rather than well-defined structures. The structural heterogeneity of IDPs, encoded by the amino acid sequence, is closely related to their unique roles in biological pathways, which require them to interact with different binding partners. Recently, Robustelli and co-workers have demonstrated that a balanced all-atom force field can be used to sample heterogeneous structures of disordered proteins ( Proc. Natl. Acad. Sci. U.S.A. 2018, 115, E4758-E4766). However, such a solution requires extensive computational resources, such as Anton supercomputers. Here, we propose a simple and effective solution to sample the conformational space of IDPs using a publicly available web server, namely, the UNited-RESidue (UNRES) web server. Our proposed solution requires no investment in computational resources and no prior knowledge of UNRES. UNRES Replica Exchange Molecular Dynamics (REMD) simulations were carried out on a set of eight disordered proteins at temperatures spanning from 270 to 430 K. Utilizing the latest UNRES force field designed for structured proteins, with proper selections of temperatures, we were able to produce comparable results to all-atom force fields as reported in work done by Robustelli and co-workers. In addition, NMR observables and the radius of gyration calculated from UNRES ensembles were directly compared with the experimental data to further evaluate the accuracy of the UNRES model at all temperatures. Our results suggest that carrying out the UNRES simulations at optimal temperatures using the UNRES web server can be a good alternative to sample heterogeneous structures of IDPs.

Biophysical and Integrative Characterization of Protein Intrinsic Disorder as a Prime Target for Drug Discovery

Protein intrinsic disorder is increasingly recognized for its biological and disease-driven functions. However, it represents significant challenges for biophysical studies due to its high conformational flexibility. In addressing these challenges, we highlight the complementary and distinct capabilities of a range of experimental and computational methods and further describe integrative strategies available for combining these techniques. Integrative biophysics methods provide valuable insights into the sequence–structure–function relationship of disordered proteins, setting the stage for protein intrinsic disorder to become a promising target for drug discovery. Finally, we briefly summarize recent advances in the development of new small molecule inhibitors targeting the disordered N-terminal domains of three vital transcription factors.

Personal Precise Force Field for Intrinsically Disordered and Ordered Proteins Based on Deep Learning

Intrinsically disordered proteins (IDPs) are proteins without a fixed three-dimensional (3D) structure under physiological conditions and are associated with Parkinson's disease, Alzheimer's disease, cancer, cardiovascular disease, amyloidosis, diabetes, and other diseases. Experimental methods can hardly capture the ensemble of diverse conformations for IDPs. Molecular dynamics (MD) simulations can sample continuous conformations that might provide a valuable complement to experimental data. However, the accuracy of MD simulations depends on the quality of force field. In particular, the evolutionary conservation and coevolution of IDPs introduce that current force fields could not precisely reproduce the conformation of IDPs. In order to improve the performance of force field, deep learning and reweighting methods were used to automatically generate personal force field parameters for intrinsically disordered and ordered proteins. At first, the deep learning method predicted more accuracy φ/ψ dihedral of residue than the previous method. Then, reweighting optimized the personal force field parameters for each residue. Finally, typical representative systems such as IDPs, structure protein, and fast-folding protein were used to evaluate this force field. The results indicate that two personal force field parameters (named PPFF1 and PPFF1_af2) could better reproduce the experimental observables than ff03CMAP force field. In summary, this strategy will provide feasibility for the development of precise personal force fields.

Performance evaluation of the balanced force field ff03CMAP for intrinsically disordered and ordered proteins

Intrinsically disordered proteins (IDPs) have been found to be closely associated with various human diseases. Because IDPs have no fixed tertiary structure under physiological conditions, current experimental methods, such as X-ray spectroscopy, NMR, and CryoEM, cannot capture all the dynamic conformations. Molecular dynamics simulation is an useful tool that is widely used to study the conformer distributions of IDPs and has become an important complementary tool for experimental methods. However, the accuracy of MD simulations directly depends on utilizing a precise force field. Recently a CMAP optimized force field based on the Amber ff03 force field (termed ff03CMAP herein) was developed for a balanced sampling of IDPs and folded proteins. In order to further evaluate the performance, more types of disordered and ordered proteins were used to test the ability for conformer sampling. The results showed that simulated chemical shifts, J-coupling, and R_g distribution with the ff03CMAP force field were in better agreement with NMR measurements and were more accurate than those with the ff03 force field. The sampling conformations by ff03CMAP were more diverse than those of ff03. At the same time, ff03CMAP could stabilize the conformers of the ordered proteins. These findings indicate that ff03CMAP can be widely used to sample diverse conformers for proteins, including the intrinsically disordered regions.

How phosphorylation impacts intrinsically disordered proteins and their function

Phosphorylation is the most common post-translational modification (PTM) in eukaryotes, occurring particularly frequently in intrinsically disordered proteins (IDPs). These proteins are highly flexible and dynamic by nature. Thus, it is intriguing that the addition of a single phosphoryl group to a disordered chain can impact its function so dramatically. Furthermore, as many IDPs carry multiple phosphorylation sites, the number of possible states increases, enabling larger complexities and novel mechanisms. Although a chemically simple and well-understood process, the impact of phosphorylation on the conformational ensemble and molecular function of IDPs, not to mention biological output, is highly complex and diverse. Since the discovery of the first phosphorylation site in proteins 75 years ago, we have come to a much better understanding of how this PTM works, but with the diversity of IDPs and their capacity for carrying multiple phosphoryl groups, the complexity grows. In this Essay, we highlight some of the basic effects of IDP phosphorylation, allowing it to serve as starting point when embarking on studies into this topic. We further describe how recent complex cases of multisite phosphorylation of IDPs have been instrumental in widening our view on the effect of protein phosphorylation. Finally, we put forward perspectives on the phosphorylation of IDPs, both in relation to disease and in context of other PTMs; areas where deep insight remains to be uncovered.

Low-resolution description of the conformational space for intrinsically disordered proteins

Intrinsically disordered proteins (IDP) are at the center of numerous biological processes, and attract consequently extreme interest in structural biology. Numerous approaches have been developed for generating sets of IDP conformations verifying a given set of experimental measurements. We propose here to perform a systematic enumeration of protein conformations, carried out using the TAiBP approach based on distance geometry. This enumeration was performed on two proteins, Sic1 and pSic1, corresponding to unphosphorylated and phosphorylated states of an IDP. The relative populations of the obtained conformations were then obtained by fitting SAXS curves as well as Ramachandran probability maps, the original finite mixture approach RamaMix being developed for this second task. The similarity between profiles of local gyration radii provides to a certain extent a converged view of the Sic1 and pSic1 conformational space. Profiles and populations are thus proposed for describing IDP conformations. Different variations of the resulting gyration radius between phosphorylated and unphosphorylated states are observed, depending on the set of enumerated conformations as well as on the methods used for obtaining the populations.

Developing Bonded Potentials for a Coarse-Grained Model of Intrinsically Disordered Proteins

Recent advances in residue-level coarse-grained (CG) computational models have enabled molecular-level insights into biological condensates of intrinsically disordered proteins (IDPs), shedding light on the sequence determinants of their phase separation. The existing CG models that treat protein chains as flexible molecules connected via harmonic bonds cannot populate common secondary-structure elements. Here, we present a CG dihedral angle potential between four neighboring beads centered at Cα atoms to faithfully capture the transient helical structures of IDPs. In order to parameterize and validate our new model, we propose Cα-based helix assignment rules based on dihedral angles that succeed in reproducing the atomistic helicity results of a polyalanine peptide and folded proteins. We then introduce sequence-dependent dihedral angle potential parameters (εd) and use experimentally available helical propensities of naturally occurring 20 amino acids to find their optimal values. The single-chain helical propensities from the CG simulations for commonly studied prion-like IDPs are in excellent agreement with the NMR-based α-helix fraction, demonstrating that the new HPS-SS model can accurately produce structural features of IDPs. Furthermore, this model can be easily implemented for large-scale assembly simulations due to its simplicity.

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.

IDPConformerGenerator: A Flexible Software Suite for Sampling the Conformational Space of Disordered Protein States

The power of structural information for informing biological mechanisms is clear for stable folded macromolecules, but similar structure–function insight is more difficult to obtain for highly dynamic systems such as intrinsically disordered proteins (IDPs) which must be described as structural ensembles. Here, we present IDPConformerGenerator, a flexible, modular open-source software platform for generating large and diverse ensembles of disordered protein states that builds conformers that obey geometric, steric, and other physical restraints on the input sequence. IDPConformerGenerator samples backbone phi (φ), psi (ψ), and omega (ω) torsion angles of relevant sequence fragments from loops and secondary structure elements extracted from folded protein structures in the RCSB Protein Data Bank and builds side chains from robust Monte Carlo algorithms using expanded rotamer libraries. IDPConformerGenerator has many user-defined options enabling variable fractional sampling of secondary structures, supports Bayesian models for assessing the agreement of IDP ensembles for consistency with experimental data, and introduces a machine learning approach to transform between internal and Cartesian coordinates with reduced error. IDPConformerGenerator will facilitate the characterization of disordered proteins to ultimately provide structural insights into these states that have key biological functions.

Salt-Induced Transitions in the Conformational Ensembles of Intrinsically Disordered Proteins

Salts modulate the behavior of intrinsically disordered proteins (IDPs) and influence the formation of membraneless organelles through liquid-liquid phase separation (LLPS). In low ionic strength solutions, IDP conformations are perturbed by the screening of electrostatic interactions, independent of the salt identity. In this regime, insight into the IDP behavior can be obtained using the theory for salt-induced transitions in charged polymers. However, salt-specific interactions with the charged and uncharged residues, known as the Hofmeister effect, influence IDP behavior in high ionic strength solutions. There is a lack of reliable theoretical models in high salt concentration regimes to predict the salt effect on IDPs. We propose a simulation methodology using a coarse-grained IDP model and experimentally measured water to salt solution transfer free energies of various chemical groups that allowed us to study the salt-specific transitions induced in the IDPs conformational ensemble. We probed the effect of three different monovalent salts on five IDPs belonging to various polymer classes based on charged residue content. We demonstrate that all of the IDPs of different polymer classes behave as self-avoiding walks (SAWs) at physiological salt concentration. In high salt concentrations, the transitions observed in the IDP conformational ensembles are dependent on the salt used and the IDP sequence and composition. Changing the anion with the cation fixed can result in the IDP transition from a SAW-like behavior to a collapsed globule. An important implication of these results is that a suitable salt can be identified to induce condensation of an IDP through LLPS.

Lighting up Nobel Prize-winning studies with protein intrinsic disorder

Intrinsically disordered proteins and regions (IDPs and IDRs) and their importance in biology are becoming increasingly recognized in biology, biochemistry, molecular biology and chemistry textbooks, as well as in current protein science and structural biology curricula. We argue that the sequence → dynamic conformational ensemble → function principle is of equal importance as the classical sequence → structure → function paradigm. To highlight this point, we describe the IDPs and/or IDRs behind the discoveries associated with 17 Nobel Prizes, 11 in Physiology or Medicine and 6 in Chemistry. The Nobel Laureates themselves did not always mention that the proteins underlying the phenomena investigated in their award-winning studies are in fact IDPs or contain IDRs. In several cases, IDP- or IDR-based molecular functions have been elucidated, while in other instances, it is recognized that the respective protein(s) contain IDRs, but the specific IDR-based molecular functions have yet to be determined. To highlight the importance of IDPs and IDRs as general principle in biology, we present here illustrative examples of IDPs/IDRs in Nobel Prize-winning mechanisms and processes.

Integrative Conformational Ensembles of Sic1 Using Different Initial Pools and Optimization Methods

Intrinsically disordered proteins play key roles in regulatory protein interactions, but their detailed structural characterization remains challenging. Here we calculate and compare conformational ensembles for the disordered protein Sic1 from yeast, starting from initial ensembles that were generated either by statistical sampling of the conformational landscape, or by molecular dynamics simulations. Two popular, yet contrasting optimization methods were used, ENSEMBLE and Bayesian Maximum Entropy, to achieve agreement with experimental data from nuclear magnetic resonance, small-angle X-ray scattering and single-molecule Förster resonance energy transfer. The comparative analysis of the optimized ensembles, including secondary structure propensity, inter-residue contact maps, and the distributions of hydrogen bond and pi interactions, revealed the importance of the physics-based generation of initial ensembles. The analysis also provides insights into designing new experiments that report on the least restrained features among the optimized ensembles. Overall, differences between ensembles optimized from different priors were greater than when using the same prior with different optimization methods. Generating increasingly accurate, reliable and experimentally validated ensembles for disordered proteins is an important step towards a mechanistic understanding of their biological function and involvement in various diseases.

Counting Degrons: Lessons From Multivalent Substrates for Targeted Protein Degradation

E3s comprise a structurally diverse group of at least 800 members, most of which target multiple substrates through specific and regulated protein-protein interactions. These interactions typically rely on short linear motifs (SLiMs), called "degrons", in an intrinsically disordered region (IDR) of the substrate, with variable rules of engagement governing different E3-docking events. These rules of engagement are of importance to the field of targeted protein degradation (TPD), where substrate ubiquitination and destruction require tools to effectively harness ubiquitin ligases (E3s). Substrates are often found to contain multiple degrons, or multiple copies of a degron, contributing to the affinity and selectivity of the substrate for its E3. One important paradigm for E3-substrate docking is presented by the Anaphase-Promoting Complex/Cyclosome (APC/C), a multi-subunit E3 ligase that targets hundreds of proteins for destruction during mitotic exit. APC/C substrate targeting takes place in an ordered manner thought to depend on tightly regulated interactions of substrates, with docking sites provided by the substoichiometric APC/C substrate adaptors and coactivators, Cdc20 or Cdh1/FZR1. Both structural and functional studies of individual APC/C substrates indicate that productive ubiquitination usually requires more than one degron, and that degrons are of different types docking to distinct sites on the coactivators. However, the dynamic nature of APC/C substrate recruitment, and the influence of multiple degrons, remains poorly understood. Here we review the significance of multiple degrons in a number of E3-substrate interactions that have been studied in detail, illustrating distinct kinetic effects of multivalency and allovalency, before addressing the role of multiple degrons in APC/C substrates, key to understanding ordered substrate destruction by APC/C. Lastly, we consider how lessons learnt from these studies can be applied in the design of TPD tools.

Assembly of model postsynaptic densities involves interactions auxiliary to stoichiometric binding

The assembly of functional biomolecular condensates often involves liquid-liquid phase separation (LLPS) of proteins with multiple modular domains, which can be folded or conformationally disordered to various degrees. To understand the LLPS-driving domain-domain interactions, a fundamental question is how readily the interactions in the condensed phase can be inferred from interdomain interactions in dilute solutions. In particular, are the interactions leading to LLPS exclusively those underlying the formation of discrete interdomain complexes in homogeneous solutions? We address this question by developing a mean-field LLPS theory of two stoichiometrically constrained solute species. The theory is applied to the neuronal proteins SynGAP and PSD-95, whose complex coacervate serves as a rudimentary model for neuronal postsynaptic densities (PSDs). The predicted phase behaviors are compared with experiments. Previously, a three SynGAP/two PSD-95 ratio was determined for SynGAP/PSD-95 complexes in dilute solutions. However, when this 3:2 stoichiometry is uniformly imposed in our theory encompassing both dilute and condensed phases, the tie-line pattern of the predicted SynGAP/PSD-95 phase diagram differs drastically from that obtained experimentally. In contrast, theories embodying alternate scenarios postulating auxiliary SynGAP-PSD-95 as well as SynGAP-SynGAP and PSD-95-PSD-95 interactions, in addition to those responsible for stoichiometric SynGAP/PSD-95 complexes, produce tie-line patterns consistent with experiment. Hence, our combined theoretical-experimental analysis indicates that weaker interactions or higher-order complexes beyond the 3:2 stoichiometry, but not yet documented, are involved in the formation of SynGAP/PSD-95 condensates, imploring future efforts to ascertain the nature of these auxiliary interactions in PSD-like LLPS and underscoring a likely general synergy between stoichiometric, structurally specific binding and stochastic, multivalent "fuzzy" interactions in the assembly of functional biomolecular condensates.

Accurate model of liquid-liquid phase behavior of intrinsically disordered proteins from optimization of single-chain properties

Many intrinsically disordered proteins (IDPs) may undergo liquid-liquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalization of intracellular biochemical reactions. The phase behavior of IDPs is sequence dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intramolecular and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.

Accurate model of liquid-liquid phase behavior of intrinsically disordered proteins from optimization of single-chain properties

Many intrinsically disordered proteins (IDPs) may undergo liquid-liquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalization of intracellular biochemical reactions. The phase behavior of IDPs is sequence dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intramolecular and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.

The Effect of Multisite Phosphorylation on the Conformational Properties of Intrinsically Disordered Proteins

Intrinsically disordered proteins are involved in many biological processes such as signaling, regulation, and recognition. A common strategy to regulate their function is through phosphorylation, as it can induce changes in conformation, dynamics, and interactions with binding partners. Although phosphorylated intrinsically disordered proteins have received increased attention in recent years, a full understanding of the conformational and structural implications of phosphorylation has not yet been achieved. Here, we present all-atom molecular dynamics simulations of five disordered peptides originated from tau, statherin, and β-casein, in both phosphorylated and non-phosphorylated state, to compare changes in global dimensions and structural elements, in an attempt to gain more insight into the controlling factors. The changes are in qualitative agreement with experimental data, and we observe that the net charge is not enough to predict the impact of phosphorylation on the global dimensions. Instead, the distribution of phosphorylated and positively charged residues throughout the sequence has great impact due to the formation of salt bridges. In statherin, a preference for arginine-phosphoserine interaction over arginine-tyrosine accounts for a global expansion, despite a local contraction of the phosphorylated region, which implies that also non-charged residues can influence the effect of phosphorylation.

Advanced Sampling Methods for Multiscale Simulation of Disordered Proteins and Dynamic Interactions

Intrinsically disordered proteins (IDPs) are highly prevalent and play important roles in biology and human diseases. It is now also recognized that many IDPs remain dynamic even in specific complexes and functional assemblies. Computer simulations are essential for deriving a molecular description of the disordered protein ensembles and dynamic interactions for a mechanistic understanding of IDPs in biology, diseases, and therapeutics. Here, we provide an in-depth review of recent advances in the multi-scale simulation of disordered protein states, with a particular emphasis on the development and application of advanced sampling techniques for studying IDPs. These techniques are critical for adequate sampling of the manifold functionally relevant conformational spaces of IDPs. Together with dramatically improved protein force fields, these advanced simulation approaches have achieved substantial success and demonstrated significant promise towards the quantitative and predictive modeling of IDPs and their dynamic interactions. We will also discuss important challenges remaining in the atomistic simulation of larger systems and how various coarse-grained approaches may help to bridge the remaining gaps in the accessible time- and length-scales of IDP simulations.

Folded domain charge properties influence the conformational behavior of disordered tails

Intrinsically disordered proteins and protein regions (IDRs) make up around 30% of the human proteome where they play essential roles in dictating and regulating many core biological processes. While IDRs are often studied as isolated domains, in naturally occurring proteins most IDRs are found adjacent to folded domains, where they exist as either N- or C-terminal tails or as linkers connecting two folded domains. Prior work has shown that charge properties of IDRs can influence their conformational behavior, both in isolation and in the context of folded domains. In contrast, the converse scenario is less well-explored: how do the charge properties of folded domains influence IDR conformational behavior? To answer this question, we combined a large-scale structural bioinformatics analysis with all-atom implicit solvent simulations of both rationally designed and naturally occurring proteins. Our results reveal three key takeaways. Firstly, the relative position and accessibility of charged residues across the surface of a folded domain can dictate IDR conformational behavior, overriding expectations based on net surface charge properties. Secondly, naturally occurring proteins possess multiple charge patches that are physically accessible to local IDRs. Finally, even modest changes in the local electrostatic environment of a folded domain can substantially modulate IDR-folded domain interactions. Taken together, our results suggest that folded domain surfaces can act as local determinants of IDR conformational behavior.

•

Intrinsically disordered regions (IDRs) are mostly found adjacent to folded domains.

•

Here we propose that the folded domain surface properties influence IDR behavior.

•

We combine all-atom simulations and sequence design of IDRs and folded domains.

•

IDR conformational behavior is determined by a complex combination of factors.

•

Folded domains can substantially alter IDR conformational biases.

Molecular Dynamics Simulations of Phosphorylated Intrinsically Disordered Proteins: A Force Field Comparison

Phosphorylation is a common post-translational modification among intrinsically disordered proteins and regions, which helps regulate function by changing the protein conformations, dynamics, and interactions with binding partners. To fully comprehend the effects of phosphorylation, computer simulations are a helpful tool, although they are dependent on the accuracy of the force field used. Here, we compared the conformational ensembles produced by Amber ff99SB-ILDN+TIP4P-D and CHARMM36m, for four phosphorylated disordered peptides ranging in length from 14-43 residues. CHARMM36m consistently produced more compact conformations with a higher content of bends, mainly due to more stable salt bridges. Based on comparisons with experimental size estimates for the shortest and longest peptide, CHARMM36m appeared to overestimate the compactness. The difference between the force fields was largest for the peptide showing the greatest separation between positively charged and phosphorylated residues, in line with the importance of charge distribution. For this peptide, the conformational ensemble did not change significantly upon increasing the ionic strength from 0 mM to 150 mM, despite a reduction of the salt-bridging probability in the CHARMM36m simulations, implying that salt concentration has negligible effects in this study.

Exploring Curated Conformational Ensembles of Intrinsically Disordered Proteins in the Protein Ensemble Database

The Protein Ensemble Database (PED; https://proteinensemble.org/) is the major repository of conformational ensembles of intrinsically disordered proteins (IDPs). Conformational ensembles of IDPs are primarily provided by their authors or occasionally collected from literature, and are subsequently deposited in PED along with the corresponding structured, manually curated metadata. The modeling of conformational ensembles usually relies on experimental data from small-angle X-ray scattering (SAXS), fluorescence resonance energy transfer (FRET), NMR spectroscopy, and molecular dynamics (MD) simulations, or a combination of these techniques. The growing number of scientific studies based on these data, along with the astounding and swift progress in the field of protein intrinsic disorder, has required a significant update and upgrade of PED, first published in 2014. To this end, the database was entirely renewed in 2020 and now has a dedicated team of biocurators providing manually curated descriptions of the methods and conditions applied to generate the conformational ensembles and for checking consistency of the data. Here, we present a detailed description on how to explore PED with its protein pages and experimental pages, and how to interpret entries of conformational ensembles. We describe how to efficiently search conformational ensembles deposited in PED by means of its web interface and API. We demonstrate how to make sense of the PED protein page and its associated experimental entry pages with reference to the yeast Sic1 use case. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Performing a search in PED Support Protocol 1: Programmatic access with the PED API Basic Protocol 2: Interpreting the protein page and the experimental entry page-the Sic1 use case Support Protocol 2: Downloading options Support Protocol 3: Understanding the validation report-the Sic1 use case Basic Protocol 3: Submitting new conformational ensembles to PED Basic Protocol 4: Providing feedback in PED.

Intrinsic disorder in integral membrane proteins

The well-defined roles and specific protein-protein interactions of many integral membrane proteins (IMPs), such as those functioning as receptors for extracellular matrix proteins and soluble growth factors, easily align with considering IMP structure as a classical "lock-and-key" concept. Nevertheless, continued advances in understanding protein conformation, such as those which established the widespread existence of intrinsically disordered proteins (IDPs) and especially intrinsically disordered regions (IDRs) in otherwise three-dimensionally organized proteins, call for ongoing reevaluation of transmembrane proteins. Here, we present basic traits of IDPs and IDRs, and, for some select single-span IMPs, consider the potential functional advantages intrinsic disorder might provide and the possible conformational impact of disease-associated mutations. For transmembrane proteins in general, we highlight several investigational approaches, such as biophysical and computational methods, stressing the importance of integrating them to produce a more-complete mechanistic model of disorder-containing IMPs. These procedures, when synergized with in-cell assessments, will likely be key in translating in silico and in vitro results to improved understanding of IMP conformational flexibility in normal cell physiology as well as disease, and will help to extend their potential as therapeutic targets.

The diversity of molecular interactions involving intrinsically disordered proteins: A molecular modeling perspective

Intrinsically Disordered Proteins and Regions (IDPs/IDRs) are key components of a multitude of biological processes. Conformational malleability enables IDPs/IDRs to perform very specialized functions that cannot be accomplished by globular proteins. The functional role for most of these proteins is related to the recognition of other biomolecules to regulate biological processes or as a part of signaling pathways. Depending on the extent of disorder, the number of interacting sites and the type of partner, very different architectures for the resulting assemblies are possible. More recently, molecular condensates with liquid-like properties composed of multiple copies of IDPs and nucleic acids have been proven to regulate key processes in eukaryotic cells. The structural and kinetic details of disordered biomolecular complexes are difficult to unveil experimentally due to their inherent conformational heterogeneity. Computational approaches, alone or in combination with experimental data, have emerged as unavoidable tools to understand the functional mechanisms of this elusive type of assemblies. The level of description used, all-atom or coarse-grained, strongly depends on the size of the molecular systems and on the timescale of the investigated mechanism. In this mini-review, we describe the most relevant architectures found for molecular interactions involving IDPs/IDRs and the computational strategies applied for their investigation.

Small-Angle X-ray Scattering Signatures of Conformational Heterogeneity and Homogeneity of Disordered Protein Ensembles

An accurate account of disordered protein conformations is of central importance to deciphering the physicochemical basis of biological functions of intrinsically disordered proteins and the folding-unfolding energetics of globular proteins. Physically, disordered ensembles of nonhomopolymeric polypeptides are expected to be heterogeneous, i.e., they should differ from those homogeneous ensembles of homopolymers that harbor an essentially unique relationship between average values of end-to-end distance R_EE and radius of gyration R_g. It was posited recently, however, that small-angle X-ray scattering (SAXS) data on conformational dimensions of disordered proteins can be rationalized almost exclusively by homopolymer ensembles. Assessing this perspective, chain-model simulations are used to evaluate the discriminatory power of SAXS-determined molecular form factors (MFFs) with regard to homogeneous versus heterogeneous ensembles. The general approach adopted here is not bound by any assumption about ensemble encodability, in that the postulated heterogeneous ensembles we evaluated are not restricted to those entailed by simple interaction schemes. Our analysis of MFFs for certain heterogeneous ensembles with more narrowly distributed R_EE and R_g indicates that while they deviate from MFFs of homogeneous ensembles, the differences can be rather small. Remarkably, some heterogeneous ensembles with asphericity and R_EE drastically different from those of homogeneous ensembles can nonetheless exhibit practically identical MFFs, demonstrating that SAXS MFFs do not afford unique characterizations of basic properties of conformational ensembles in general. In other words, the ensemble to MFF mapping is practically many-to-one and likely nonsmooth. Heteropolymeric variations of the R_EE-R_g relationship were further showcased using an analytical perturbation theory developed here for flexible heteropolymers. Ramifications of our findings for interpretation of experimental data are discussed.

Extensive evaluation of environment-specific force field for ordered and disordered proteins

Intrinsically disordered proteins (IDPs) have no fixed tertiary structure under physiological conditions and are associated with many human diseases. Because IDPs have the characteristic of possessing diverse conformations, current experimental methods cannot capture all the conformations of IDPs. However, molecular dynamics simulation can sample these atomistically diverse conformations as a valuable complement to experimental data. To accurately describe the properties of IDPs, the environment-specific precise force field (ESFF1) was successfully released to reproduce the conformer character of ordered and disordered proteins. Here, three typical IDPs and thirteen folded proteins were used to further evaluate the performance of this force field. The results indicate that the NMR observables of ESFF1 better approach experimental data than do those of ff14SB for IDPs. The sampling conformations by ESFF1 are more diverse than those of ff14SB. For folded proteins, these force fields have comparable performances for reproducing conformers. Therefore, ESFF1 can be used to reveal the model of sequence-disorder-function for IDPs.

Full structural ensembles of intrinsically disordered proteins from unbiased molecular dynamics simulations

Molecular dynamics (MD) simulation is widely used to complement ensemble-averaged experiments of intrinsically disordered proteins (IDPs). However, MD often suffers from limitations of inaccuracy. Here, we show that enhancing the sampling using Hamiltonian replica-exchange MD (HREMD) led to unbiased and accurate ensembles, reproducing small-angle scattering and NMR chemical shift experiments, for three IDPs of varying sequence properties using two recently optimized force fields, indicating the general applicability of HREMD for IDPs. We further demonstrate that, unlike HREMD, standard MD can reproduce experimental NMR chemical shifts, but not small-angle scattering data, suggesting chemical shifts are insufficient for testing the validity of IDP ensembles. Surprisingly, we reveal that despite differences in their sequence, the inter-chain statistics of all three IDPs are similar for short contour lengths (< 10 residues). The results suggest that the major hurdle of generating an accurate unbiased ensemble for IDPs has now been largely overcome.

An analytical theory to describe sequence-specific inter-residue distance profiles for polyampholytes and intrinsically disordered proteins

Intrinsically Disordered Proteins (IDPs), unlike folded proteins, lack a unique folded structure and rapidly interconvert among ensembles of disordered states. However, they have specific conformational properties when averaged over their ensembles of disordered states. It is critical to develop a theoretical formalism to predict these ensemble average conformational properties that are encoded in the IDP sequence (the specific order in which amino acids/residues are linked). We present a general heteropolymer theory that analytically computes the ensemble average distance profiles (⟨R_ij ²⟩) between any two (i, j) monomers (amino acids for IDPs) as a function of the sequence. Information rich distance profiles provide a detailed description of the IDP in contrast to typical metrics such as scaling exponents, radius of gyration, or end-to-end distance. This generalized formalism supersedes homopolymer-like models or models that are built only on the composition of amino acids but ignore sequence details. The prediction of these distance profiles for highly charged polyampholytes and naturally occurring IDPs unmasks salient features that are hidden in the sequence. Moreover, the model reveals strategies to modulate the entire distance map to achieve local or global swelling/compaction by subtle changes/modifications-such as phosphorylation, a biologically relevant process-in specific hotspots in the sequence. Sequence-specific distance profiles and their modulation have been benchmarked against all-atom simulations. Our new formalism also predicts residue-pair specific coil-globule transitions. The analytical nature of the theory will facilitate design of new sequences to achieve specific target distance profiles with broad applications in synthetic biology and polymer science.

How GPCR Phosphorylation Patterns Orchestrate Arrestin-Mediated Signaling

Binding of arrestin to phosphorylated G-protein-coupled receptors (GPCRs) controls many aspects of cell signaling. The number and arrangement of phosphates may vary substantially for a given GPCR, and different phosphorylation patterns trigger different arrestin-mediated effects. Here, we determine how GPCR phosphorylation influences arrestin behavior by using atomic-level simulations and site-directed spectroscopy to reveal the effects of phosphorylation patterns on arrestin binding and conformation. We find that patterns favoring binding differ from those favoring activation-associated conformational change. Both binding and conformation depend more on arrangement of phosphates than on their total number, with phosphorylation at different positions sometimes exerting opposite effects. Phosphorylation patterns selectively favor a wide variety of arrestin conformations, differently affecting arrestin sites implicated in scaffolding distinct signaling proteins. We also reveal molecular mechanisms of these phenomena. Our work reveals the structural basis for the long-standing "barcode" hypothesis and has important implications for design of functionally selective GPCR-targeted drugs.

Mechanisms of Macromolecular Interactions Mediated by Protein Intrinsic Disorder

Intrinsically disordered proteins or regions (IDPs or IDRs) are widespread in the eukaryotic proteome. Although lacking stable three-dimensional structures in the free forms, IDRs perform critical functions in various cellular processes. Accordingly, mutations and altered expression of IDRs are associated with many pathological conditions. Hence, it is of great importance to understand at the molecular level how IDRs interact with their binding partners. In particular, discovering the unique interaction features of IDRs originating from their dynamic nature may reveal uncharted regulatory mechanisms of specific biological processes. Here we discuss the mechanisms of the macromolecular interactions mediated by IDRs and present the relevant cellular processes including transcription, cell cycle progression, signaling, and nucleocytoplasmic transport. Of special interest is the multivalent binding nature of IDRs driving assembly of multicomponent macromolecular complexes. Integrating the previous theoretical and experimental investigations, we suggest that such IDR-driven multiprotein complexes can function as versatile allosteric switches to process diverse cellular signals. Finally, we discuss the future challenges and potential medical applications of the IDR research.

Conformational Ensembles of an Intrinsically Disordered Protein Consistent with NMR, SAXS, and Single-Molecule FRET

Intrinsically disordered proteins (IDPs) have fluctuating heterogeneous conformations, which makes their structural characterization challenging. Although challenging, characterization of the conformational ensembles of IDPs is of great interest, since their conformational ensembles are the link between their sequences and functions. An accurate description of IDP conformational ensembles depends crucially on the amount and quality of the experimental data, how it is integrated, and if it supports a consistent structural picture. We used integrative modeling and validation to apply conformational restraints and assess agreement with the most common structural techniques for IDPs: Nuclear Magnetic Resonance (NMR) spectroscopy, Small-angle X-ray Scattering (SAXS), and single-molecule Förster Resonance Energy Transfer (smFRET). Agreement with such a diverse set of experimental data suggests that details of the generated ensembles can now be examined with a high degree of confidence. Using the disordered N-terminal region of the Sic1 protein as a test case, we examined relationships between average global polymeric descriptions and higher-moments of their distributions. To resolve apparent discrepancies between smFRET and SAXS inferences, we integrated SAXS data with NMR data and reserved the smFRET data for independent validation. Consistency with smFRET, which was not guaranteed a priori, indicates that, globally, the perturbative effects of NMR or smFRET labels on the Sic1 ensemble are minimal. Analysis of the ensembles revealed distinguishing features of Sic1, such as overall compactness and large end-to-end distance fluctuations, which are consistent with biophysical models of Sic1's ultrasensitive binding to its partner Cdc4. Our results underscore the importance of integrative modeling and validation in generating and drawing conclusions from IDP conformational ensembles.

Analysis of Protein Disorder Predictions in the Light of a Protein Structural Alphabet

Intrinsically-disordered protein (IDP) characterization was an amazing change of paradigm in our classical sequence-structure-function theory. Moreover, IDPs are over-represented in major disease pathways and are now often targeted using small molecules for therapeutic purposes. This has had created a complex continuum from order-that encompasses rigid and flexible regions-to disorder regions; the latter being not accessible through classical crystallographic methodologies. In X-ray structures, the notion of order is dictated by access to resolved atom positions, providing rigidity and flexibility information with low and high experimental B-factors, while disorder is associated with the missing (non-resolved) residues. Nonetheless, some rigid regions can be found in disorder regions. Using ensembles of IDPs, their local conformations were analyzed in the light of a structural alphabet. An entropy index derived from this structural alphabet allowed us to propose a continuum of states from rigidity to flexibility and finally disorder. In this study, the analysis was extended to comparing these results to disorder predictions, underlying a limited correlation, and so opening new ideas to characterize and predict disorder.

Dynamical Behavior and Conformational Selection Mechanism of the Intrinsically Disordered Sic1 Kinase-Inhibitor Domain

Intrinsically Disordered Peptides and Proteins (IDPs) in solution can span a broad range of conformations that often are hard to characterize by both experimental and computational methods. However, obtaining a significant representation of the conformational space is important to understand mechanisms underlying protein functions such as partner recognition. In this work, we investigated the behavior of the Sic1 Kinase-Inhibitor Domain (KID) in solution by Molecular Dynamics (MD) simulations. Our results point out that application of common descriptors of molecular shape such as Solvent Accessible Surface (SAS) area can lead to misleading outcomes. Instead, more appropriate molecular descriptors can be used to define 3D structures. In particular, we exploited Weighted Holistic Invariant Molecular (WHIM) descriptors to get a coarse-grained but accurate definition of the variegated Sic1 KID conformational ensemble. We found that Sic1 is able to form a variable amount of folded structures even in absence of partners. Among them, there were some conformations very close to the structure that Sic1 is supposed to assume in the binding with its physiological complexes. Therefore, our results support the hypothesis that this protein relies on the conformational selection mechanism to recognize the correct molecular partners.

Intrinsic disorder in proteins: Relevance to protein assemblies, drug design and host-pathogen interactions

Intrinsic disorder in proteins resulting in considerable variation in structure can lead to multiple functions including multi-specificity and diverse pathologies. Protein interfaces can involve disordered regions that assemble through a concerted-fold-and-bind mechanism. The binding involves both enthalpic and entropic gains by exploiting 'hot spots' on the partner and displacing water molecules placed in thermodynamically unfavorable situations. The examples of Rad51-BRCA2 and Artemis-DNA-PKCs/LigIV complexes illustrate this in the context of drug design. This overview tracks the seamless involvement of protein disorder in multi-specificity of biocatalysts, protein assembly formations and host-pathogen interactions, where intrinsic disorder can in Mycobacteria, compensate for genome reduction by carrying out multiple functions and in some RNA viruses facilitate adaption to the host. These present challenging opportunities for designing new drugs and interventions.

Phosphorylation orchestrates the structural ensemble of the intrinsically disordered protein HMGA1a and modulates its DNA binding to the NFκB promoter

High Mobility Group Protein A1a (HMGA1a) is a highly abundant nuclear protein, which plays a crucial role during embryogenesis, cell differentiation, and neoplasia. Here, we present the first ever NMR-based structural ensemble of full length HMGA1a. Our results show that the protein is not completely random coil but adopts a compact structure consisting of transient long-range contacts, which is regulated by post-translational phosphorylation. The CK2-, cdc2- and cdc2/CK2-phosphorylated forms of HMGA1a each exhibit a different binding affinity towards the PRD2 element of the NFκB promoter. Our study identifies connected regions between phosphorylation sites in the wildtype ensemble that change considerably upon phosphorylation, indicating that these posttranslational modifications sites are part of an electrostatic contact network that alters the structural ensemble by shifting the conformational equilibrium. Moreover, ITC data reveal that the CK2-phosphorylated HMGA1a exhibits a different DNA promoter binding affinity for the PRD2 element. Furthermore, we present the first structural model for AT-hook 1 of HMGA1a that can adopt a transient α-helical structure, which might serve as an additional regulatory mechanism in HMAG1a. Our findings will help to develop new therapeutic strategies against HMGA1a-associated cancers by taking posttranslational modifications into consideration.

Distance-Based Metrics for Comparing Conformational Ensembles of Intrinsically Disordered Proteins

Intrinsically disordered proteins are proteins whose native functional states represent ensembles of highly diverse conformations. Such ensembles are a challenge for quantitative structure comparisons because their conformational diversity precludes optimal superimposition of the atomic coordinates necessary for deriving common similarity measures such as the root mean-square deviation of these coordinates. Here, we introduce superimposition-free metrics that are based on computing matrices of the Cα-Cα distance distributions within ensembles and comparing these matrices between ensembles. Differences between two matrices yield information on the similarity between specific regions of the polypeptide, whereas the global structural similarity is captured by the root mean-square difference between the medians of the Cα-Cα distance distributions of two ensembles. Together, our metrics enable rigorous investigations of structure-function relationships in conformational ensembles of intrinsically disordered proteins derived using experimental restraints or by molecular simulations and for proteins containing both structured and disordered regions.

Common Functions of Disordered Proteins across Evolutionary Distant Organisms

Intrinsically disordered proteins and regions typically lack a well-defined structure and thus fall outside the scope of the classic sequence–structure–function relationship. Hence, classic sequence- or structure-based bioinformatic approaches are often not well suited to identify homology or predict the function of unknown intrinsically disordered proteins. Here, we give selected examples of intrinsic disorder in plant proteins and present how protein function is shared, altered or distinct in evolutionary distant organisms. Furthermore, we explore how examining the specific role of disorder across different phyla can provide a better understanding of the common features that protein disorder contributes to the respective biological mechanism.

Genome re-seqeunce and analysis of Burkholderia glumae strain AU6208 and evidence of toxoflavin: A potential bacterial toxin

Burkholderia glumae, the primary causative agent of bacterial panicle blight in rice, has been reported as an opportunistic pathogen in patients with chronic infections. This study aimed to re-sequence the clinical isolate B. glumae strain AU6208 and comparatively analyze its genome using B. glumae strain BGR1 from rice plant as the reference. Re-sequencing results revealed that the genome of strain AU6208 comprised 96 contigs corresponding to a 6.1 Mbp genome of the strain AU6208, with 5322 coding sequences and 68.2 % GC content; this is much larger compared to the genome previously sequenced by us and described by Seo et al (2015), which was reported to be 4.1 Mbp comprising >1200 contigs, 4361 coding sequences, and 67.31 % GC content. Moreover, this updated genome shares >80 % identity to the 7.2 Mbp genome of BGR1, which encodes 6491 coding sequences and has 68.3 % GC content. Further computational analysis revealed that the strain AU6208 encodes several bacteriocin biosynthesis genes, antibiotic, as well as virulent genes such as toxoflavin genes, which included 425 specialty genes and 12 toxoflavin genes. Upon further characterization, 12 toxoflavins (ToxA, B, C, D, E, F, G, H, I, J, TofI, and TofR) were found in AU6208 with 70-100 % sequence, family, and domain similarity with that of BGR1. Upon comparison with BGR1, the structural characterizations of selected toxoflavin genes (ToxB, ToxC, ToxG, H, and TofI) revealed variations in 2D and 3D structures such as differences in α-helix, β-sheets, loops, physiological properties of proteins, RMSD values, etc. These variations may play significant role in different mode of action in different hosts thereby indicating that in addition to their respective hosts, toxoflavins could also contribute to exploit other hosts across the kingdom. In addition to understanding the epidemiology of strain AU6208, this updated genomics data will also unfold the pathogenicity of bacteria in diversity of various hosts and anti-virulence.

Intrinsically disordered regions regulate the activities of ATP binding cassette transporters

ATP binding cassette (ABC) proteins are a large family of membrane proteins present in all kingdoms of life. These multi-domain proteins are comprised, at minimum, of two membrane-spanning domains (MSD1, MSD2) and two cytosolic nucleotide binding domains (NBD1, NBD2). ATP binding and hydrolysis at the NBDs enables ABC proteins to actively transport solutes across membranes, regulate activities of other proteins, or function as channels. Like most eukaryotic membrane proteins, ABC proteins contain intrinsically disordered regions (IDRs). These conformationally dynamic regions in ABC proteins possess residual structure, are sites of phosphorylation, and mediate protein-protein interactions. Here, we review the role of IDRs in regulating ABC protein activity.

Computing, Analyzing, and Comparing the Radius of Gyration and Hydrodynamic Radius in Conformational Ensembles of Intrinsically Disordered Proteins

The level of compaction of an intrinsically disordered protein may affect both its physical and biological properties, and can be probed via different types of biophysical experiments. Small-angle X-ray scattering (SAXS) probe the radius of gyration (R_g) whereas pulsed-field-gradient nuclear magnetic resonance (NMR) diffusion, fluorescence correlation spectroscopy, and dynamic light scattering experiments can be used to determine the hydrodynamic radius (R_h). Here we show how to calculate R_g and R_h from a computationally generated conformational ensemble of an intrinsically disordered protein. We further describe how to use a Bayesian/Maximum Entropy procedure to integrate data from SAXS and NMR diffusion experiments, so as to derive conformational ensembles in agreement with those experiments.