Vitamin K Epoxide Reductase Complex-Protein Disulphide Isomerase Assemblies in the Thiol-Disulphide Exchange Reactions: Portrayal of Precursor-to-Successor Complexes

The human Vitamin K Epoxide Reductase Complex (hVKORC1), a key enzyme that converts vitamin K into the form necessary for blood clotting, requires for its activation the reducing equivalents supplied by its redox partner through thiol-disulphide exchange reactions. The functionally related molecular complexes assembled during this process have never been described, except for a proposed de novo model of a 'precursor' complex of hVKORC1 associated with protein disulphide isomerase (PDI). Using numerical approaches (in silico modelling and molecular dynamics simulation), we generated alternative 3D models for each molecular complex bonded either covalently or non-covalently. These models differ in the orientation of the PDI relative to hVKORC1 and in the cysteine residue involved in forming protein-protein disulphide bonds. Based on a comparative analysis of these models' shape, folding, and conformational dynamics, the most probable putative complexes, mimicking the 'precursor', 'intermediate', and 'successor' states, were suggested. In addition, we propose using these complexes to develop the 'allo-network drugs' necessary for treating blood diseases.

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.

Leucine Motifs Stabilize Residual Helical Structure in Disordered Proteins

Many examples are known of regions of intrinsically disordered proteins that fold into α-helices upon binding to their targets. These helical binding motifs (HBMs) can be partially helical also in the unbound state, and this so-called residual structure can affect binding affinity and kinetics. To investigate the underlying mechanisms governing the formation of residual helical structure, we assembled a dataset of experimental helix contents of 65 peptides containing HBM that fold-upon-binding. The average residual helicity is 17% and increases to 60% upon target binding. The helix contents of residual and target-bound structures do not correlate, however the relative location of helix elements in both states shows a strong overlap. Compared to the general disordered regions, HBMs are enriched in amino acids with high helix preference and these residues are typically involved in target binding, explaining the overlap in helix positions. In particular, we find that leucine residues and leucine motifs in HBMs are the major contributors to helix stabilization and target-binding. For the two model peptides, we show that substitution of leucine motifs to other hydrophobic residues (valine or isoleucine) leads to reduction of residual helicity, supporting the role of leucine as helix stabilizer. From the three hydrophobic residues only leucine can efficiently stabilize residual helical structure. We suggest that the high occurrence of leucine motifs and a general preference for leucine at binding interfaces in HBMs can be explained by its unique ability to stabilize helical elements.

Synergy of Mutation-Induced Effects in Human Vitamin K Epoxide Reductase: Perspectives and Challenges for Allo-Network Modulator Design

The human Vitamin K Epoxide Reductase Complex (hVKORC1), a key enzyme transforming vitamin K into the form necessary for blood clotting, requires for its activation the reducing equivalents delivered by its redox partner through thiol-disulfide exchange reactions. The luminal loop (L-loop) is the principal mediator of hVKORC1 activation, and it is a region frequently harbouring numerous missense mutations. Four L-loop hVKORC1 mutants, suggested in vitro as either resistant (A41S, H68Y) or completely inactive (S52W, W59R), were studied in the oxidised state by numerical approaches (in silico). The DYNASOME and POCKETOME of each mutant were characterised and compared to the native protein, recently described as a modular protein composed of the structurally stable transmembrane domain (TMD) and the intrinsically disordered L-loop, exhibiting quasi-independent dynamics. The DYNASOME of mutants revealed that L-loop missense point mutations impact not only its folding and dynamics, but also those of the TMD, highlighting a strong mutation-specific interdependence between these domains. Another consequence of the mutation-induced effects manifests in the global changes (geometric, topological, and probabilistic) of the newly detected cryptic pockets and the alternation of the recognition properties of the L-loop with its redox protein. Based on our results, we postulate that (i) intra-protein allosteric regulation and (ii) the inherent allosteric regulation and cryptic pockets of each mutant depend on its DYNASOME; and (iii) the recognition of the redox protein by hVKORC1 (INTERACTOME) depend on their DYNASOME. This multifaceted description of proteins produces "omics" data sets, crucial for understanding the physiological processes of proteins and the pathologies caused by alteration of the protein properties at various "omics" levels. Additionally, such characterisation opens novel perspectives for the development of "allo-network drugs" essential for the treatment of blood disorders.

Intrinsic disorder of a nucleoplasmin-like histone chaperone specifies its discrete nuclear and nucleolar functions

Nucleoplasmin (NPM) histone chaperones regulate distinct processes in the nucleus and nucleolus. While intrinsically disordered regions (IDRs) are hallmarks of NPMs, it is not clear whether all NPM functions require these unstructured features. We assessed the importance of IDRs in a yeast NPM-like protein and found that regulation of rDNA copy number and genetic interactions with the nucleolar RNA surveillance machinery require the highly conserved FKBP prolyl isomerase domain, but not the NPM domain or IDRs. By contrast, transcriptional repression in the nucleus requires IDRs. Furthermore, multiple lysines in polyacidic serine/lysine motifs of IDRs are required for both lysine polyphosphorylation and NPM-mediated transcriptional repression. These results demonstrate that this NPM-like protein relies on IDRs only for some of its chromatin-related functions.

Biophysical characterization of the Plasmodium falciparum circumsporozoite protein's N-terminal domain

The circumsporozoite protein (CSP) is the main surface antigen of the Plasmodium sporozoite (SPZ) and forms the basis of the currently only licensed anti-malarial vaccine (RTS,S/AS01). CSP uniformly coats the SPZ and plays a pivotal role in its immunobiology, in both the insect and the vertebrate hosts. Although CSP's N-terminal domain (CSPN ) has been reported to play an important role in multiple CSP functions, a thorough biophysical and structural characterization of CSPN is currently lacking. Here, we present an alternative method for the recombinant production and purification of CSPN from Plasmodium falciparum (PfCSPN ), which provides pure, high-quality protein preparations with high yields. Through an interdisciplinary approach combining in-solution experimental methods and in silico analyses, we provide strong evidence that PfCSPN is an intrinsically disordered region displaying some degree of compaction.

IndiSPENsable for X Chromosome Inactivation and Gene Silencing

For about 30 years, SPEN has been the subject of research in many different fields due to its variety of functions and its conservation throughout a wide spectrum of species, like worms, arthropods, and vertebrates. To date, 216 orthologues have been documented. SPEN had been studied for its role in gene regulation in the context of cell signaling, including the NOTCH or nuclear hormone receptor signaling pathways. More recently, SPEN has been identified as a major regulator of initiation of chromosome-wide gene silencing during X chromosome inactivation (XCI) in mammals, where its function remains to be fully understood. Dependent on the biological context, SPEN functions via mechanisms which include different domains. While some domains of SPEN are highly conserved in sequence and secondary structure, species-to-species differences exist that might lead to mechanistic differences. Initiation of XCI appears to be different between humans and mice, which raises additional questions about the extent of generalization of SPEN's function in XCI. In this review, we dissect the mechanism of SPEN in XCI. We discuss its subregions and domains, focusing on its role as a major regulator. We further highlight species-related research, specifically of mouse and human SPEN, with the aim to reveal and clarify potential species-to-species differences in SPEN's function.

CoMemMoRFPred: Sequence-based Prediction of MemMoRFs by Combining Predictors of Intrinsic Disorder, MoRFs and Disordered Lipid-binding Regions

Molecular recognition features (MoRFs) are a commonly occurring type of intrinsically disordered regions (IDRs) that undergo disorder-to-order transition upon binding to partner molecules. We focus on recently characterized and functionally important membrane-binding MoRFs (MemMoRFs). Motivated by the lack of computational tools that predict MemMoRFs, we use a dataset of experimentally annotated MemMoRFs to conceptualize, design, evaluate and release an accurate sequence-based predictor. We rely on state-of-the-art tools that predict residues that possess key characteristics of MemMoRFs, such as intrinsic disorder, disorder-to-order transition and lipid-binding. We identify and combine results from three tools that include flDPnn for the disorder prediction, DisoLipPred for the prediction of disordered lipid-binding regions, and MoRFCHiBiLight for the prediction of disorder-to-order transitioning protein binding regions. Our empirical analysis demonstrates that combining results produced by these three methods generates accurate predictions of MemMoRFs. We also show that use of a smoothing operator produces predictions that closely mimic the number and sizes of the native MemMoRF regions. The resulting CoMemMoRFPred method is available as an easy-to-use webserver at http://biomine.cs.vcu.edu/servers/CoMemMoRFPred. This tool will aid future studies of MemMoRFs in the context of exploring their abundance, cellular functions, and roles in pathologic phenomena.

PICKLE RELATED 2 is a Neofunctionalized Gene Duplicate Under Positive Selection With Antagonistic Effects to the Ancestral PICKLE Gene on the Seed Transcriptome

The evolution and diversification of proteins capable of remodeling domains has been critical for transcriptional reprogramming during cell fate determination in multicellular eukaryotes. Chromatin remodeling proteins of the CHD3 family have been shown to have important and antagonistic impacts on seed development in the model plant, Arabidopsis thaliana, yet the basis of this functional divergence remains unknown. In this study, we demonstrate that genes encoding the CHD3 proteins PICKLE (PKL) and PICKLE-RELATED 2 (PKR2) originated from a duplication event during the diversification of crown Brassicaceae, and that these homologs have undergone distinct evolutionary trajectories since this duplication, with PKR2 fast evolving under positive selection, while PKL is subject to purifying selection. We find that the rapid evolution of PKR2 under positive selection reduces the encoded protein's intrinsic disorder, possibly suggesting a tertiary structure configuration which differs from that of PKL. Our whole genome transcriptome analysis in seeds of pkr2 and pkl mutants reveals that they act antagonistically on the expression of specific sets of genes, providing a basis for their differing roles in seed development. Our results provide insights into how gene duplication and neofunctionalization can lead to differing and antagonistic selective pressures on transcriptomes during plant reproduction, as well as on the evolutionary diversification of the CHD3 family within seed plants.

Implications of intrinsic disorder and functional proteomics in the merkel cell polyomavirus life cycle

Infection with merkel cell polyomavirus (MCPyV) is implicated in the development of merkel cell carcinoma (MCC), a rare but aggressive skin cancer. MCC has a mortality rate near 50%, and incidence has been rapidly increasing in recent decades, making development of improved treatment strategies critical to addressing its growing social burden. The parallel increasing necessity for novel research to better understand MCPyV pathogenesis has prompted numerous studies in recent years, yet the role of intrinsic disorder in MCPyV proteins remains unexplored. This study carries out computational characterization of intrinsic disorder within the MCPyV proteome and suggests mechanisms that may contribute to the oncogenicity of the virus to invade and hijack host immune systems. Our analysis finds that significant levels of intrinsic disorder are present in proteins LT, ALTO, 57kT, and VP1, and suggests that regions of sT may also contain large, disordered regions. The investigation further shows correlation of disorder propensity with the outputs for functional predictors of eukaryotic linear motifs (ELMs), molecular recognition features (MoRFs), and propensity for liquid-liquid phase separation (LLPS). Our findings indicate that MCPyV may use disorder and phase condensation to alter viral function that may accentuate or provide the basis for oncogenic activities. It is intended that this study will inform future experimental validation efforts around the phase separation capacity of MCPyV and its host protein-protein interactions. Furthermore, we hope to inform other investigators on the potential role of disorder in the MCPyV life cycle toward ultimately progressing the development of novel therapeutic agents.

It's Time for Entropic Clocks: The Roles of Random Chain Protein Sequences in Timing Ion Channel Processes Underlying Action Potential Properties

In recent years, it has become clear that intrinsically disordered protein segments play diverse functional roles in many cellular processes, thus leading to a reassessment of the classical structure-function paradigm. One class of intrinsically disordered protein segments is entropic clocks, corresponding to unstructured random protein chains involved in timing cellular processes. Such clocks were shown to modulate ion channel processes underlying action potential generation, propagation, and transmission. In this review, we survey the role of entropic clocks in timing intra- and inter-molecular binding events of voltage-activated potassium channels involved in gating and clustering processes, respectively, and where both are known to occur according to a similar 'ball and chain' mechanism. We begin by delineating the thermodynamic and timing signatures of a 'ball and chain'-based binding mechanism involving entropic clocks, followed by a detailed analysis of the use of such a mechanism in the prototypical Shaker voltage-activated K+ channel model protein, with particular emphasis on ion channel clustering. We demonstrate how 'chain'-level alternative splicing of the Kv channel gene modulates entropic clock-based 'ball and chain' inactivation and clustering channel functions. As such, the Kv channel model system exemplifies how linkage between alternative splicing and intrinsic disorder enables the functional diversity underlying changes in electrical signaling.

Hierarchical Assembly of Intrinsically Disordered Short Peptides

The understanding on how short peptide assemblies transit from disorder to order remains limited due to the lack of atomistic structures. Here we report cryo-EM structure of the nanofibers short intrinsically disordered peptides (IDPs). Upon lowering pH or adding calcium ions, the IDP transitions from individual nanoparticles to nanofibers containing an aromatic core and a disordered periphery comprised of 2 to 5 amino acids. Protonating the phosphate or adding more metal ions further assembles the nanofibers into filament bundles. The assemblies of the IDP analogs with controlled chemistry, such as phosphorylation site, hydrophobic interactions, and sequences indicate that metal ions interact with the flexible periphery of the nanoparticles of the IDPs to form fibrils and enhance the interfibrillar interactions to form filament bundles. Illustrating that an IDP self-assembles from disorder to order, this work offers atomistic molecular insights to understand assemblies of short peptides driven by noncovalent interactions.

Battle of the bugs: how an oomycete pathogen shapes the microbiota of its host

This article is a Commentary on Gómez‐Pérez et al. (2023), 239: 2320–2334.

Intrinsic Disorder in the Human Tear Proteome

Purpose

We aimed to characterize the proteome of human tears and assess for the presence of intrinsically disordered proteins (IDPs). IDPs, despite lacking a rigid three-dimensional structure, maintain biological functionality and could shed light on the molecular interactions within tears.

Methods

We analyzed a dataset of 1475 proteins identified in the tear film of three healthy subjects. We employed several computational tools, including the Compositional Profiler, Rapid Intrinsic Disorder Analysis Online, Search Tool for the Retrieval of Interacting Genes, and Database of Disordered Protein Predictors to evaluate the intrinsic disorder, protein interactions, and functional characterization of the disordered regions within this proteome.

Results

Our analysis showed a notable inclination toward intrinsic disorder. Two out of 10 order-promoting residues and five out of 10 disorder-promoting residues were found enriched. Using the Predictor of Natural Disordered Regions (PONDR) VSL2 output, 95% of these proteins were classified as highly or moderately disordered. We revealed an extensive protein-protein interaction network with significant interaction enrichment. The most disordered proteins exhibited higher disorder binding sites and diverse posttranslational modifications compared to the most ordered ones.

Conclusions

To the best of our knowledge, our study is the first comprehensive analysis of intrinsic disorder in the human tear film proteome, and it revealed an abundance of IDPs and their role in protein function and interaction networks. These findings suggest that variations in the intrinsic disorder of a tear film could be impacted by systemic and ocular conditions, offering promising avenues for disease biomarker identification and drug target development. Further research is needed to understand the implications of these findings in human health and disease.

Overview Update: Computational Prediction of Intrinsic Disorder in Proteins

There are over 100 computational predictors of intrinsic disorder. These methods predict amino acid-level propensities for disorder directly from protein sequences. The propensities can be used to annotate putative disordered residues and regions. This unit provides a practical and holistic introduction to the sequence-based intrinsic disorder prediction. We define intrinsic disorder, explain the format of computational prediction of disorder, and identify and describe several accurate predictors. We also introduce recently released databases of intrinsic disorder predictions and use an illustrative example to provide insights into how predictions should be interpreted and combined. Lastly, we summarize key experimental methods that can be used to validate computational predictions. © 2023 Wiley Periodicals LLC.

Intrinsic Disorder in α-Synuclein Regulates the Exocytotic Fusion Pore Transition

Today, it is widely accepted that intrinsic disorder is strongly related to the cell cycle, during mitosis, differentiation, and apoptosis. Of particular interest are hybrid proteins possessing both structured and unstructured domains that are critical in human health and disease, such as α-synuclein. In this work, we describe how α-synuclein interacts with the nascent fusion pore as it evolves toward expansion. We unveil the key role played by its intrinsically disordered region as a thermodynamic regulator of the nucleation-expansion energy barrier. By analyzing a truncated variant of α-synuclein that lacks the disordered region, we find that the landscape of protein interactions with PIP₂ and POPS lipids is highly altered, ultimately increasing the energy cost for the fusion pore to transit from nucleation to expansion. We conclude that the intrinsically disordered region in full-length α-synuclein recognizes and allocates pivotal protein:lipid interactions during membrane remodeling in the first stages of the fusion pore.

Prediction of protein-protein interactions using sequences of intrinsically disordered regions

Protein-protein interactions (PPIs) play a crucial role in numerous molecular processes. Despite many efforts, mechanisms governing molecular recognition between interacting proteins remain poorly understood and it is particularly challenging to predict from sequence whether two proteins can interact. Here we present a new method to tackle this challenge using intrinsically disordered regions (IDRs). IDRs are protein segments that are functional despite lacking a single invariant three-dimensional structure. The prevalence of IDRs in eukaryotic proteins suggests that IDRs are critical for interactions. To test this hypothesis, we predicted PPIs using IDR sequences in candidate proteins in humans. Moreover, we divide the PPI prediction problem into two specific subproblems and adapt appropriate training and test strategies based on problem type. Our findings underline the importance of defining clearly the problem type and show that sequences encoding IDRs can aid in predicting specific features of the protein interaction network of intrinsically disordered proteins. Our findings further suggest that accounting for IDRs in future analyses should accelerate efforts to elucidate the eukaryotic PPI network.

Linker Length Drives Heterogeneity of Multivalent Complexes of Hub Protein LC8 and Transcription Factor ASCIZ

LC8, a ubiquitous and highly conserved hub protein, binds over 100 proteins involved in numerous cellular functions, including cell death, signaling, tumor suppression, and viral infection. LC8 binds intrinsically disordered proteins (IDPs), and although several of these contain multiple LC8 binding motifs, the effects of multivalency on complex formation are unclear. Drosophila ASCIZ has seven motifs that vary in sequence and inter-motif linker lengths, especially within subdomain QT2-4 containing the second, third, and fourth LC8 motifs. Using isothermal-titration calorimetry, analytical-ultracentrifugation, and native mass-spectrometry of QT2-4 variants, with methodically deactivated motifs, we show that inter-motif spacing and specific motif sequences combine to control binding affinity and compositional heterogeneity of multivalent duplexes. A short linker separating strong and weak motifs results in stable duplexes but forms off-register structures at high LC8 concentrations. Contrastingly, long linkers engender lower cooperativity and heterogeneous complexation at low LC8 concentrations. Accordingly, two-mers, rather than the expected three-mers, dominate negative-stain electron-microscopy images of QT2-4. Comparing variants containing weak-strong and strong-strong motif combinations demonstrates sequence also regulates IDP/LC8 assembly. The observed trends persist for trivalent ASCIZ subdomains: QT2-4, with long and short linkers, forms heterogeneous complexes, whereas QT4-6, with similar mid-length linkers, forms homogeneous complexes. Implications of linker length variations for function are discussed.

Thermo-Programmed Synthetic DNA-Based Receptors

Herein, we present a generalizable and versatile strategy to engineer synthetic DNA ligand-binding devices that can be programmed to load and release a specific ligand at a defined temperature. We do so by re-engineering two model DNA-based receptors: a triplex-forming bivalent DNA-based receptor that recognizes a specific DNA sequence and an ATP-binding aptamer. The temperature at which these receptors load/release their ligands can be finely modulated by controlling the entropy associated with the linker connecting the two ligand-binding domains. The availability of a set of receptors with tunable and reversible temperature dependence allows achieving complex load/release behavior such as sustained ligand release over a wide temperature range. Similar programmable thermo-responsive synthetic ligand-binding devices can be of utility in applications such as drug delivery and production of smart materials.

Disordered region encodes α-crystallin chaperone activity toward lens client γD-crystallin

Small heat-shock proteins (sHSPs) are a widely expressed family of ATP-independent molecular chaperones that are among the first responders to cellular stress. Mechanisms by which sHSPs delay aggregation of client proteins remain undefined. sHSPs have high intrinsic disorder content of up to ~60% and assemble into large, polydisperse homo- and hetero-oligomers, making them challenging structural and biochemical targets. Two sHSPs, HSPB4 and HSPB5, are present at millimolar concentrations in eye lens, where they are responsible for maintaining lens transparency over the lifetime of an organism. Together, HSPB4 and HSPB5 compose the hetero-oligomeric chaperone known as α-crystallin. To identify the determinants of sHSP function, we compared the effectiveness of HSPB4 and HSPB5 homo-oligomers and HSPB4/HSPB5 hetero-oligomers in delaying the aggregation of the lens protein γD-crystallin. In chimeric versions of HSPB4 and HSPB5, chaperone activity tracked with the identity of the 60-residue disordered N-terminal regions (NTR). A short 10-residue stretch in the middle of the NTR ("Critical sequence") contains three residues that are responsible for high HSPB5 chaperone activity toward γD-crystallin. These residues affect structure and dynamics throughout the NTR. Abundant interactions involving the NTR Critical sequence reveal it to be a hub for a network of interactions within oligomers. We propose a model whereby the NTR critical sequence influences local structure and NTR dynamics that modulate accessibility of the NTR, which in turn modulates chaperone activity.

CLIP: accurate prediction of disordered linear interacting peptides from protein sequences using co-evolutionary information

One of key features of intrinsically disordered regions (IDRs) is facilitation of protein-protein and protein-nucleic acids interactions. These disordered binding regions include molecular recognition features (MoRFs), short linear motifs (SLiMs) and longer binding domains. Vast majority of current predictors of disordered binding regions target MoRFs, with a handful of methods that predict SLiMs and disordered protein-binding domains. A new and broader class of disordered binding regions, linear interacting peptides (LIPs), was introduced recently and applied in the MobiDB resource. LIPs are segments in protein sequences that undergo disorder-to-order transition upon binding to a protein or a nucleic acid, and they cover MoRFs, SLiMs and disordered protein-binding domains. Although current predictors of MoRFs and disordered protein-binding regions could be used to identify some LIPs, there are no dedicated sequence-based predictors of LIPs. To this end, we introduce CLIP, a new predictor of LIPs that utilizes robust logistic regression model to combine three complementary types of inputs: co-evolutionary information derived from multiple sequence alignments, physicochemical profiles and disorder predictions. Ablation analysis suggests that the co-evolutionary information is particularly useful for this prediction and that combining the three inputs provides substantial improvements when compared to using these inputs individually. Comparative empirical assessments using low-similarity test datasets reveal that CLIP secures area under receiver operating characteristic curve (AUC) of 0.8 and substantially improves over the results produced by the closest current tools that predict MoRFs and disordered protein-binding regions. The webserver of CLIP is freely available at http://biomine.cs.vcu.edu/servers/CLIP/ and the standalone code can be downloaded from http://yanglab.qd.sdu.edu.cn/download/CLIP/.

DeepPRObind: Modular Deep Learner that Accurately Predicts Structure and Disorder-Annotated Protein Binding Residues

Current sequence-based predictors of protein-binding residues (PBRs) belong to two distinct categories: structure-trained vs. intrinsic disorder-trained. Since disordered PBRs differ from structured PBRs in several ways, including ability to bind multiple partners by folding into different conformations and enrichment in different amino acids, the structure-trained and disorder-trained predictors were shown to provide inaccurate results for the other annotation type. A simple consensus-based solution that combines structure- and disorder-trained methods provides limited levels of predictive performance and generates relatively many cross-predictions, where residues that interact with other ligand types are predicted as PBRs. We address this unsolved problem by designing a novel and fast deep-learner, DeepPRObind, that relies on carefully designed modular convolutional architecture and uses innovative aggregate input features. Comparative empirical tests on a low-similarity test dataset reveal that DeepPRObind generates accurate predictions of structured and disordered PBRs and low amounts of cross-predictions, outperforming a comprehensive collection of 12 predictors of PBRs. Given the relatively low runtime of DeepPRObind (40 seconds per protein), we further validate its results based on an analysis of putative PBRs in the yeast proteome, confirming that interactions in disordered regions are enriched among hub proteins. We release DeepPRObind as a convenient web server at https://www.csuligroup.com/DeepPRObind/.

Phase separation in innate immune response and inflammation-related diseases

Inflammation induced by nonspecific pathogenic or endogenous danger signals is an essential mechanism of innate immune response. The innate immune responses are rapidly triggered by conserved germline-encoded receptors that recognize broad patterns indicative of danger, with subsequent signal amplification by modular effectors, which have been the subject of intense investigation for many years. Until recently, however, the critical role of intrinsic disorder-driven phase separation in facilitating innate immune responses went largely unappreciated. In this review, we discuss emerging evidences that many innate immune receptors, effectors, and/or interactors function as "all-or-nothing" switch-like hubs to stimulate acute and chronic inflammation. By concentrating or relegating modular signaling components to phase-separated compartments, cells construct flexible and spatiotemporal distributions of key signaling events to ensure rapid and effective immune responses to a myriad of potentially harmful stimuli.

The Hydrophilic Loop of Arabidopsis PIN1 Auxin Efflux Carrier Harbors Hallmarks of an Intrinsically Disordered Protein

Much of plant development depends on cell-to-cell redistribution of the plant hormone auxin, which is facilitated by the plasma membrane (PM) localized PIN FORMED (PIN) proteins. Auxin export activity, developmental roles, subcellular trafficking, and polarity of PINs have been well studied, but their structure remains elusive besides a rough outline that they contain two groups of 5 alpha-helices connected by a large hydrophilic loop (HL). Here, we focus on the PIN1 HL as we could produce it in sufficient quantities for biochemical investigations to provide insights into its secondary structure. Circular dichroism (CD) studies revealed its nature as an intrinsically disordered protein (IDP), manifested by the increase of structure content upon thermal melting. Consistent with IDPs serving as interaction platforms, PIN1 loops homodimerize. PIN1 HL cytoplasmic overexpression in Arabidopsis disrupts early endocytic trafficking of PIN1 and PIN2 and causes defects in the cotyledon vasculature formation. In summary, we demonstrate that PIN1 HL has an intrinsically disordered nature, which must be considered to gain further structural insights. Some secondary structures may form transiently during pairing with known and yet-to-be-discovered interactors.

Generating Ensembles of Dynamic Misfolding Proteins

The early stages of protein misfolding and aggregation involve disordered and partially folded protein conformers that contain a high degree of dynamic disorder. These dynamic species may undergo large-scale intra-molecular motions of intrinsically disordered protein (IDP) precursors, or flexible, low affinity inter-molecular binding in oligomeric assemblies. In both cases, generating atomic level visualization of the interconverting species that captures the conformations explored and their physico-chemical properties remains hugely challenging. How specific sub-ensembles of conformers that are on-pathway to aggregation into amyloid can be identified from their aggregation-resilient counterparts within these large heterogenous pools of rapidly moving molecules represents an additional level of complexity. Here, we describe current experimental and computational approaches designed to capture the dynamic nature of the early stages of protein misfolding and aggregation, and discuss potential challenges in describing these species because of the ensemble averaging of experimental restraints that arise from motions on the millisecond timescale. We give a perspective of how machine learning methods can be used to extract aggregation-relevant sub-ensembles and provide two examples of such an approach in which specific interactions of defined species within the dynamic ensembles of α-synuclein (αSyn) and β2-microgloblulin (β2m) can be captured and investigated.

Intrinsically disordered regions in the rodent hepevirus proteome

Hepatitis E virus (HEV) is the causative agent of Hepatitis E infections across the world. Intrinsically disordered protein regions (IDPRs) or intrinsically disordered proteins (IDPs) are regions or proteins that are characterized by lack of definite structure. These IDPRs or IDPs play significant roles in a wide range of biological processes, such as cell cycle regulation, control of signaling pathways, etc. IDPR/IDP in proteins is associated with the virus's pathogenicity and infectivity. The prevalence of IDPR/IDP in rat HEV proteome remains undetermined. Hence, we examined the unstructured/disordered regions of the open reading frame (ORF) encoded proteins of rat HEV by analyzing the prevalence of intrinsic disorder. The intrinsic disorder propensity analysis showed that the different ORF proteins consisted of varying fraction of intrinsic disorder. The protein ORF3 was identified with maximum propensity for intrinsic disorder while the ORF6 protein had the least fraction of intrinsic disorder. The analysis revealed ORF6 as a structured protein (ORDP); ORF1 and ORF4 as moderately disordered proteins (IDPRs); and ORF3 and ORF5 as highly disordered proteins (IDPs). The protein ORF2 was found to be moderately as well as highly disordered using different predictors, thus, was categorized into both IDPR and IDP. Such disordered regions have important roles in pathogenesis and replication of viruses.

Surveying over 100 predictors of intrinsic disorder in proteins

INTRODUCTION

Intrinsic disorder prediction field develops, assesses, and deploys computational predictors of disorder in protein sequences and constructs and disseminates databases of these predictions. Over 40 years of research resulted in the release of numerous resources.

AREAS COVERED

We identify and briefly summarize the most comprehensive to date collection of over 100 disorder predictors. We focus on their predictive models, availability and predictive performance. We categorize and study them from a historical point of view to highlight informative trends.

EXPERT OPINION

We find a consistent trend of improvements in predictive quality as newer and more advanced predictors are developed. The original focus on machine learning methods has shifted to meta-predictors in early 2010s, followed by a recent transition to deep learning. The use of deep learners will continue in foreseeable future given recent and convincing success of these methods. Moreover, a broad range of resources that facilitate convenient collection of accurate disorder predictions is available to users. They include web servers and standalone programs for disorder prediction, servers that combine prediction of disorder and disorder functions, and large databases of pre-computed predictions. We also point to the need to address the shortage of accurate methods that predict disordered binding regions.

Conserved and divergent features of neuronal CaMKII holoenzyme structure, function, and high-order assembly

Neuronal CaMKII holoenzymes (α and β isoforms) enable molecular signal computation underlying learning and memory but also mediate excitotoxic neuronal death. Here, we provide a comparative analysis of these signaling devices, using single-particle electron microscopy (EM) in combination with biochemical and live-cell imaging studies. In the basal state, both isoforms assemble mainly as 12-mers (but also 14-mers and even 16-mers for the β isoform). CaMKIIα and β isoforms adopt an ensemble of extended activatable states (with average radius of 12.6 versus 16.8 nm, respectively), characterized by multiple transient intra- and inter-holoenzyme interactions associated with distinct functional properties. The extended state of CaMKIIβ allows direct resolution of intra-holoenzyme kinase domain dimers. These dimers could enable cooperative activation by calmodulin, which is observed for both isoforms. High-order CaMKII clustering mediated by inter-holoenzyme kinase domain dimerization is reduced for the β isoform for both basal and excitotoxicity-induced clusters, both in vitro and in neurons.

The CaMKII holoenzyme enables neuronal signal computation. In a comparative structure-function analysis of the neuronal α and β isoforms, Buonarati et al. find evidence for kinase domain dimers within the holoenzyme that enable a cooperative activation mechanism in both isoforms and inter-holoenzyme interactions that enable high-order aggregate formation under ischemic conditions.

DeepDISOBind: accurate prediction of RNA-, DNA- and protein-binding intrinsically disordered residues with deep multi-task learning

Proteins with intrinsically disordered regions (IDRs) are common among eukaryotes. Many IDRs interact with nucleic acids and proteins. Annotation of these interactions is supported by computational predictors, but to date, only one tool that predicts interactions with nucleic acids was released, and recent assessments demonstrate that current predictors offer modest levels of accuracy. We have developed DeepDISOBind, an innovative deep multi-task architecture that accurately predicts deoxyribonucleic acid (DNA)-, ribonucleic acid (RNA)- and protein-binding IDRs from protein sequences. DeepDISOBind relies on an information-rich sequence profile that is processed by an innovative multi-task deep neural network, where subsequent layers are gradually specialized to predict interactions with specific partner types. The common input layer links to a layer that differentiates protein- and nucleic acid-binding, which further links to layers that discriminate between DNA and RNA interactions. Empirical tests show that this multi-task design provides statistically significant gains in predictive quality across the three partner types when compared to a single-task design and a representative selection of the existing methods that cover both disorder- and structure-trained tools. Analysis of the predictions on the human proteome reveals that DeepDISOBind predictions can be encoded into protein-level propensities that accurately predict DNA- and RNA-binding proteins and protein hubs. DeepDISOBind is available at https://www.csuligroup.com/DeepDISOBind/.

The Nucleocapsid of Paramyxoviruses: Structure and Function of an Encapsidated Template

Viruses of the Paramyxoviridae family share a common and complex molecular machinery for transcribing and replicating their genomes. Their non-segmented, negative-strand RNA genome is encased in a tight homopolymer of viral nucleoproteins (N). This ribonucleoprotein complex, termed a nucleocapsid, is the template of the viral polymerase complex made of the large protein (L) and its co-factor, the phosphoprotein (P). This review summarizes the current knowledge on several aspects of paramyxovirus transcription and replication, including structural and functional data on (1) the architecture of the nucleocapsid (structure of the nucleoprotein, interprotomer contacts, interaction with RNA, and organization of the disordered C-terminal tail of N), (2) the encapsidation of the genomic RNAs (structure of the nucleoprotein in complex with its chaperon P and kinetics of RNA encapsidation in vitro), and (3) the use of the nucleocapsid as a template for the polymerase complex (release of the encased RNA and interaction network allowing the progress of the polymerase complex). Finally, this review presents models of paramyxovirus transcription and replication.

Expanding the Disorder-Function Paradigm in the C-Terminal Tails of Erbbs

ErbBs are receptor tyrosine kinases involved not only in development, but also in a wide variety of diseases, particularly cancer. Their extracellular, transmembrane, juxtamembrane, and kinase folded domains were described extensively over the past 20 years, structurally and functionally. However, their whole C-terminal tails (CTs) following the kinase domain were only described at atomic resolution in the last 4 years. They were shown to be intrinsically disordered. The CTs are known to be tyrosine-phosphorylated when the activated homo- or hetero-dimers of ErbBs are formed. Their phosphorylation triggers interaction with phosphotyrosine binding (PTB) or Src Homology 2 (SH2) domains and activates several signaling pathways controling cellular motility, proliferation, adhesion, and apoptosis. Beyond this passive role of phosphorylated domain and site display for partners, recent structural and function studies unveiled active roles in regulation of phosphorylation and interaction: the CT regulates activity of the kinase domain; different phosphorylation states have different compaction levels, potentially modulating the succession of phosphorylation events; and prolines have an important role in structure, dynamics, and possibly regulatory interactions. Here, we review both the canonical role of the disordered CT domains of ErbBs as phosphotyrosine display domains and the recent findings that expand the known range of their regulation functions linked to specific structural and dynamic features.

Intrinsic Disorder in Human RNA-Binding Proteins

Although RNA-binding proteins (RBPs) are known to be enriched in intrinsic disorder, no previous analysis focused on RBPs interacting with specific RNA types. We fill this gap with a comprehensive analysis of the putative disorder in RBPs binding to six common RNA types: messenger RNA (mRNA), transfer RNA (tRNA), small nuclear RNA (snRNA), non-coding RNA (ncRNA), ribosomal RNA (rRNA), and internal ribosome RNA (irRNA). We also analyze the amount of putative intrinsic disorder in the RNA-binding domains (RBDs) and non-RNA-binding-domain regions (non-RBD regions). Consistent with previous studies, we show that in comparison with human proteome, RBPs are significantly enriched in disorder. However, closer examination finds significant enrichment in predicted disorder for the mRNA-, rRNA- and snRNA-binding proteins, while the proteins that interact with ncRNA and irRNA are not enriched in disorder, and the tRNA-binding proteins are significantly depleted in disorder. We show a consistent pattern of significant disorder enrichment in the non-RBD regions coupled with low levels of disorder in RBDs, which suggests that disorder is relatively rarely utilized in the RNA-binding regions. Our analysis of the non-RBD regions suggests that disorder harbors posttranslational modification sites and is involved in the putative interactions with DNA. Importantly, we utilize experimental data from DisProt and independent data from Pfam to validate the above observations that rely on the disorder predictions. This study provides new insights into the distribution of disorder across proteins that bind different RNA types and the functional role of disorder in the regions where it is enriched.

Intrinsically disordered proteins: modes of binding with emphasis on disordered domains

Our notions of protein function have long been determined by the protein structure-function paradigm. However, the idea that protein function is dictated by a prerequisite complementarity of shapes at the binding interface is becoming increasingly challenged. Interactions involving intrinsically disordered proteins (IDPs) have indicated a significant degree of disorder present in the bound state, ranging from static disorder to complete disorder, termed 'random fuzziness'. This review assesses the anatomy of an IDP and relates how its intrinsic properties permit promiscuity and allow for the various modes of interaction. Furthermore, a mechanistic overview of the types of disordered domains is detailed, while also relating to a recent example and the kinetic and thermodynamic principles governing its formation.

The sequence-ensemble relationship in fuzzy protein complexes

Intrinsically disordered proteins (IDPs) interact with globular proteins through a variety of mechanisms, resulting in the structurally heterogeneous ensembles known as fuzzy complexes. While there exists a reasonable comprehension on how IDP sequence determines the unbound IDP ensemble, little is known about what shapes the structural characteristics of IDPs bound to their targets. Using a statistical thermodynamic model, we show that the target-bound ensembles are determined by a simple code that combines the IDP sequence and the distribution of IDP-target interaction hotspots. These two parameters define the conformational space of target-bound IDPs and rationalize the observed structural heterogeneity of fuzzy complexes. The presented model successfully reproduces the dynamical signatures of target-bound IDPs from the NMR relaxation experiments as well as the changes of interaction affinity and the IDP helicity induced by mutations. The model explains how the target-bound IDP ensemble adapts to mutations in order to achieve an optimal balance between conformational freedom and interaction energy. Taken together, the presented sequence-ensemble relationship of fuzzy complexes explains the different manifestations of IDP disorder in folding-upon-binding processes.

Coarse-Grained Modeling and Molecular Dynamics Simulations of Ca2+-Calmodulin

Calmodulin (CaM) is a calcium-binding protein that transduces signals to downstream proteins through target binding upon calcium binding in a time-dependent manner. Understanding the target binding process that tunes CaM’s affinity for the calcium ions (Ca²⁺), or vice versa, may provide insight into how Ca²⁺-CaM selects its target binding proteins. However, modeling of Ca²⁺-CaM in molecular simulations is challenging because of the gross structural changes in its central linker regions while the two lobes are relatively rigid due to tight binding of the Ca²⁺ to the calcium-binding loops where the loop forms a pentagonal bipyramidal coordination geometry with Ca²⁺. This feature that underlies the reciprocal relation between Ca²⁺ binding and target binding of CaM, however, has yet to be considered in the structural modeling. Here, we presented a coarse-grained model based on the Associative memory, Water mediated, Structure, and Energy Model (AWSEM) protein force field, to investigate the salient features of CaM. Particularly, we optimized the force field of CaM and that of Ca²⁺ ions by using its coordination chemistry in the calcium-binding loops to match with experimental observations. We presented a “community model” of CaM that is capable of sampling various conformations of CaM, incorporating various calcium-binding states, and carrying the memory of binding with various targets, which sets the foundation of the reciprocal relation of target binding and Ca²⁺ binding in future studies.

GPIHBP1 and ANGPTL4 Utilize Protein Disorder to Orchestrate Order in Plasma Triglyceride Metabolism and Regulate Compartmentalization of LPL Activity

Intravascular processing of triglyceride-rich lipoproteins (TRLs) is crucial for delivery of dietary lipids fueling energy metabolism in heart and skeletal muscle and for storage in white adipose tissue. During the last decade, mechanisms underlying focal lipolytic processing of TRLs along the luminal surface of capillaries have been clarified by fresh insights into the functions of lipoprotein lipase (LPL); LPL's dedicated transporter protein, glycosylphosphatidylinositol-anchored high density lipoprotein-binding protein 1 (GPIHBP1); and its endogenous inhibitors, angiopoietin-like (ANGPTL) proteins 3, 4, and 8. Key discoveries in LPL biology include solving the crystal structure of LPL, showing LPL is catalytically active as a monomer rather than as a homodimer, and that the borderline stability of LPL's hydrolase domain is crucial for the regulation of LPL activity. Another key discovery was understanding how ANGPTL4 regulates LPL activity. The binding of ANGPTL4 to LPL sequences adjacent to the catalytic cavity triggers cooperative and sequential unfolding of LPL's hydrolase domain resulting in irreversible collapse of the catalytic cavity and loss of LPL activity. Recent studies have highlighted the importance of the ANGPTL3-ANGPTL8 complex for endocrine regulation of LPL activity in oxidative organs (e.g., heart, skeletal muscle, brown adipose tissue), but the molecular mechanisms have not been fully defined. New insights have also been gained into LPL-GPIHBP1 interactions and how GPIHBP1 moves LPL to its site of action in the capillary lumen. GPIHBP1 is an atypical member of the LU (Ly6/uPAR) domain protein superfamily, containing an intrinsically disordered and highly acidic N-terminal extension and a disulfide bond-rich three-fingered LU domain. Both the disordered acidic domain and the folded LU domain are crucial for the stability and transport of LPL, and for modulating its susceptibility to ANGPTL4-mediated unfolding. This review focuses on recent advances in the biology and biochemistry of crucial proteins for intravascular lipolysis.

Zn2+-Dependent Histone Deacetylases in Plants: Structure and Evolution

Zn²⁺-dependent histone deacetylases are widely distributed in archaea, bacteria, and eukaryotes. Through deacetylation of histones and other biomolecules, these enzymes regulate mammalian gene expression, microtubule stability, and polyamine metabolism. In plants, they play essential roles in development and stress response, but little is known about their biochemistry. We provide here a holistic revision of plant histone deacetylase (HDA) phylogeny and translate recent lessons from other organisms. HDA evolution correlates with a gain of structural ductility/disorder, as observed for other proteins. We also highlight two recently identified Brassicaceae-specific HDAs, as well as unprecedented key mutations that would affect the catalytic activity of individual HDAs. This revised phylogeny will contextualize future studies and illuminate research on plant development and adaptation.

Ubiquitin Interacting Motifs: Duality Between Structured and Disordered Motifs

Ubiquitin is a small protein at the heart of many cellular processes, and several different protein domains are known to recognize and bind ubiquitin. A common motif for interaction with ubiquitin is the Ubiquitin Interacting Motif (UIM), characterized by a conserved sequence signature and often found in multi-domain proteins. Multi-domain proteins with intrinsically disordered regions mediate interactions with multiple partners, orchestrating diverse pathways. Short linear motifs for binding are often embedded in these disordered regions and play crucial roles in modulating protein function. In this work, we investigated the structural propensities of UIMs using molecular dynamics simulations and NMR chemical shifts. Despite the structural portrait depicted by X-crystallography of stable helical structures, we show that UIMs feature both helical and intrinsically disordered conformations. Our results shed light on a new class of disordered UIMs. This group is here exemplified by the C-terminal domain of one isoform of ataxin-3 and a group of ubiquitin-specific proteases. Intriguingly, UIMs not only bind ubiquitin. They can be a recruitment point for other interactors, such as parkin and the heat shock protein Hsc70-4. Disordered UIMs can provide versatility and new functions to the client proteins, opening new directions for research on their interactome.

Protein Aggregation Landscape in Neurodegenerative Diseases: Clinical Relevance and Future Applications

Intrinsic disorder is a natural feature of polypeptide chains, resulting in the lack of a defined three-dimensional structure. Conformational changes in intrinsically disordered regions of a protein lead to unstable β-sheet enriched intermediates, which are stabilized by intermolecular interactions with other β-sheet enriched molecules, producing stable proteinaceous aggregates. Upon misfolding, several pathways may be undertaken depending on the composition of the amino acidic string and the surrounding environment, leading to different structures. Accumulating evidence is suggesting that the conformational state of a protein may initiate signalling pathways involved both in pathology and physiology. In this review, we will summarize the heterogeneity of structures that are produced from intrinsically disordered protein domains and highlight the routes that lead to the formation of physiological liquid droplets as well as pathogenic aggregates. The most common proteins found in aggregates in neurodegenerative diseases and their structural variability will be addressed. We will further evaluate the clinical relevance and future applications of the study of the structural heterogeneity of protein aggregates, which may aid the understanding of the phenotypic diversity observed in neurodegenerative disorders.

Parallel and Sequential Pathways of Molecular Recognition of a Tandem-Repeat Protein and Its Intrinsically Disordered Binding Partner

The Wnt signalling pathway plays an important role in cell proliferation, differentiation, and fate decisions in embryonic development and the maintenance of adult tissues. The twelve armadillo (ARM) repeat-containing protein β-catenin acts as the signal transducer in this pathway. Here, we investigated the interaction between β-catenin and the intrinsically disordered transcription factor TCF7L2, comprising a very long nanomolar-affinity interface of approximately 4800 Å2 that spans ten of the twelve ARM repeats of β-catenin. First, a fluorescence reporter system for the interaction was engineered and used to determine the kinetic rate constants for the association and dissociation. The association kinetics of TCF7L2 and β-catenin were monophasic and rapid (7.3 ± 0.1 × 107 M-1·s-1), whereas dissociation was biphasic and slow (5.7 ± 0.4 × 10-4 s-1, 15.2 ± 2.8 × 10-4 s-1). This reporter system was then combined with site-directed mutagenesis to investigate the striking variability in the conformation adopted by TCF7L2 in the three different crystal structures of the TCF7L2-β-catenin complex. We found that the mutation had very little effect on the association kinetics, indicating that most interactions form after the rate-limiting barrier for association. Mutations of the N- and C-terminal subdomains of TCF7L2 that adopt relatively fixed conformations in the crystal structures had large effects on the dissociation kinetics, whereas the mutation of the labile sub-domain connecting them had negligible effect. These results point to a two-site avidity mechanism of binding with the linker region forming a "fuzzy" complex involving transient contacts that are not site-specific. Strikingly, the two mutations in the N-terminal subdomain that had the largest effects on the dissociation kinetics showed two additional phases, indicating partial flux through an alternative dissociation pathway that is inaccessible to the wild type. The results presented here provide insights into the kinetics of the molecular recognition of a long intrinsically disordered region with an elongated repeat-protein surface, a process found to involve parallel routes with sequential steps in each.

Serotonin Promotes Serum Albumin Interaction with the Monomeric Amyloid β Peptide

Prevention of amyloid β peptide (Aβ) deposition via facilitation of Aβ binding to its natural depot, human serum albumin (HSA), is a promising approach to preclude Alzheimer's disease (AD) onset and progression. Previously, we demonstrated the ability of natural HSA ligands, fatty acids, to improve the affinity of this protein to monomeric Aβ by a factor of 3 (BBRC, 510(2), 248-253). Using plasmon resonance spectroscopy, we show here that another HSA ligand related to AD pathogenesis, serotonin (SRO), increases the affinity of the Aβ monomer to HSA by a factor of 7/17 for Aβ₄₀/Aβ₄₂, respectively. Meanwhile, the structurally homologous SRO precursor, tryptophan (TRP), does not affect HSA's affinity to monomeric Aβ, despite slowdown of the association and dissociation processes. Crosslinking with glutaraldehyde and dynamic light scattering experiments reveal that, compared with the TRP-induced effects, SRO binding causes more marked changes in the quaternary structure of HSA. Furthermore, molecular docking reveals distinct structural differences between SRO/TRP complexes with HSA. The disintegration of the serotonergic system during AD pathogenesis may contribute to Aβ release from HSA in the central nervous system due to impairment of the SRO-mediated Aβ trapping by HSA.

Abundance Imparts Evolutionary Constraints of Similar Magnitude on the Buried, Surface, and Disordered Regions of Proteins

An understanding of the forces shaping protein conservation is key, both for the fundamental knowledge it represents and to allow for optimal use of evolutionary information in practical applications. Sequence conservation is typically examined at one of two levels. The first is a residue-level, where intra-protein differences are analyzed and the second is a protein-level, where inter-protein differences are studied. At a residue level, we know that solvent-accessibility is a prime determinant of conservation. By inverting this logic, we inferred that disordered regions are slightly more solvent-accessible on average than the most exposed surface residues in domains. By integrating abundance information with evolutionary data within and across proteins, we confirmed a previously reported strong surface-core association in the evolution of structured regions, but we found a comparatively weak association between disordered and structured regions. The facts that disordered and structured regions experience different structural constraints and evolve independently provide a unique setup to examine an outstanding question: why is a protein’s abundance the main determinant of its sequence conservation? Indeed, any structural or biophysical property linked to the abundance-conservation relationship should increase the relative conservation of regions concerned with that property (e.g., disordered residues with mis-interactions, domain residues with misfolding). Surprisingly, however, we found the conservation of disordered and structured regions to increase in equal proportion with abundance. This observation implies that either abundance-related constraints are structure-independent, or multiple constraints apply to different regions and perfectly balance each other.

Discovery of Arginine Methylation, Phosphorylation, and Their Co-occurrence in Condensate-Associated Proteins in Saccharomyces cerevisiae

The formation of condensates in membraneless organelles is thought to be driven by protein phase separation. Arginine methylation and serine/threonine phosphorylation are important in the phase separation process; however, these post-translational modifications are often present in intrinsically disordered regions that are difficult to analyze with standard proteomic techniques. To understand their presence and co-occurrence in condensate-associated proteins, here, we use a multiprotease and multi-tandem mass spectrometry (MS/MS) fragmentation approach, coupled with heavy methyl stable isotope labeling of amino acids in cell culture (SILAC) and phospho- or methyl-peptide enrichment. For Saccharomyces cerevisiae, we report a 50% increase in the known arginine methylproteome, involving 15 proteins that are all condensate-associated. Importantly, some of these proteins have arginine methylation on all predicted sites-providing evidence that this modification can be pervasive. We explored whether arginine-methylated, condensate-associated proteins are also phosphorylated and found 12 such proteins to carry phosphorylated serine or threonine. In Npl3, Ded1, and Sbp1, single peptides were found to carry both modifications, indicating a co-occurrence in close proximity and on the same protein molecule. These co-modifications occur in regions of disorder, whereas arginine methylation is typically on regions of disorder that are also basic. For phosphorylation, its association with charged regions of condensate-associated proteins was less consistent, although some regions with multisite phosphorylation sites were strongly acidic. We conclude that arginine-methylated proteins associated with condensates are typically also modified with protein phosphorylation.

The intrinsic instability of the hydrolase domain of lipoprotein lipase facilitates its inactivation by ANGPTL4-catalyzed unfolding.Proc Natl Acad Sci U S A. 2021;118. [PMID: 33723082 DOI: 10.1073/pnas.2026650118]

The complex between lipoprotein lipase (LPL) and its endothelial receptor (GPIHBP1) is responsible for the lipolytic processing of triglyceride-rich lipoproteins (TRLs) along the capillary lumen, a physiologic process that releases lipid nutrients for vital organs such as heart and skeletal muscle. LPL activity is regulated in a tissue-specific manner by endogenous inhibitors (angiopoietin-like [ANGPTL] proteins 3, 4, and 8), but the molecular mechanisms are incompletely understood. ANGPTL4 catalyzes the inactivation of LPL monomers by triggering the irreversible unfolding of LPL's α/β-hydrolase domain. Here, we show that this unfolding is initiated by the binding of ANGPTL4 to sequences near LPL's catalytic site, including β2, β3-α3, and the lid. Using pulse-labeling hydrogen‒deuterium exchange mass spectrometry, we found that ANGPTL4 binding initiates conformational changes that are nucleated on β3-α3 and progress to β5 and β4-α4, ultimately leading to the irreversible unfolding of regions that form LPL's catalytic pocket. LPL unfolding is context dependent and varies with the thermal stability of LPL's α/β-hydrolase domain (T _m of 34.8 °C). GPIHBP1 binding dramatically increases LPL stability (T _m of 57.6 °C), while ANGPTL4 lowers the onset of LPL unfolding by ∼20 °C, both for LPL and LPL•GPIHBP1 complexes. These observations explain why the binding of GPIHBP1 to LPL retards the kinetics of ANGPTL4-mediated LPL inactivation at 37 °C but does not fully suppress inactivation. The allosteric mechanism by which ANGPTL4 catalyzes the irreversible unfolding and inactivation of LPL is an unprecedented pathway for regulating intravascular lipid metabolism.

Rational Control of the Activity of a Cu2+-Dependent DNAzyme by Re-engineering Purely Entropic Intrinsically Disordered Domains

The function and activity of many proteins is finely controlled by the modulation of the entropic contribution of intrinsically disordered domains that are not directly involved in any recognition event. Inspired by this mechanism, we demonstrate here that we could finely regulate the catalytic activity of a model DNAzyme (i.e., a synthetic DNA sequence with enzyme-like properties) by rationally introducing intrinsically disordered nucleic acid portions in its original sequence. More specifically, we have re-engineered here the well-characterized Cu²⁺-dependent DNAzyme that catalyzes a self-cleavage reaction by introducing a poly(T) linker domain in its sequence. The linker is not directly involved in the recognition event and connects the two domains that fold to form the catalytic core. We demonstrate that the enzyme-like activity of this re-engineered DNAzyme can be modulated in a predictable and fine way by changing the length, and thus entropy, of such a linker domain. Given these attributes, the rational design of intrinsically disordered domains could expand the available toolbox to achieve a control of the activity of DNAzymes and, in analogy, ribozymes through a purely entropic contribution.

FLIPPER: Predicting and Characterizing Linear Interacting Peptides in the Protein Data Bank

A large fraction of peptides or protein regions are disordered in isolation and fold upon binding. These regions, also called MoRFs, SLiMs or LIPs, are often associated with signaling and regulation processes. However, despite their importance, only a limited number of examples are available in public databases and their automatic detection at the proteome level is problematic. Here we present FLIPPER, an automatic method for the detection of structurally linear sub-regions or peptides that interact with another chain in a protein complex. FLIPPER is a random forest classification that takes the protein structure as input and provides the propensity of each amino acid to be part of a LIP region. Models are built taking into consideration structural features such as intra- and inter-chain contacts, secondary structure, solvent accessibility in both bound and unbound state, structural linearity and chain length. FLIPPER is accurate when evaluated on non-redundant independent datasets, 99% precision and 99% sensitivity on PixelDB-25 and 87% precision and 88% sensitivity on DIBS-25. Finally, we used FLIPPER to process the entire Protein Data Bank and identified different classes of LIPs based on different binding modes and partner molecules. We provide a detailed description of these LIP categories and show that a large fraction of these regions are not detected by disorder predictors. All FLIPPER predictions are integrated in the MobiDB 4.0 database.

Structural and Functional Characterization of the ABA-Water Deficit Stress Domain from Wheat and Barley: An Intrinsically Disordered Domain behind the Versatile Functions of the Plant Abscissic Acid, Stress and Ripening Protein Family

The ASR protein family has been discovered thirty years ago in many plant species and is involved in the tolerance of various abiotic stresses such as dehydration, salinity and heat. Despite its importance, nothing is known about the conserved ABA-Water Deficit Stress Domain (ABA-WDS) of the ASR gene family. In this study, we characterized two ABA-WDS domains, isolated from durum wheat (TtABA-WDS) and barley (HvABA-WDS). Bioinformatics analysis shows that they are both consistently predicted to be intrinsically disordered. Hydrodynamic and circular dichroism analysis indicate that both domains are largely disordered but belong to different structural classes, with HvABA-WDS and TtABA-WDS adopting a PreMolten Globule-like (PMG-like) and a Random Coil-like (RC-like) conformation, respectively. In the presence of the secondary structure stabilizer trifluoroethanol (TFE) or of increasing glycerol concentrations, which mimics dehydration, the two domains acquire an α-helical structure. Interestingly, both domains are able to prevent heat- and dehydration-induced inactivation of the enzyme lactate dehydrogenase (LDH). Furthermore, heterologous expression of TtABA-WDS and HvABA-WDS in the yeast Saccharomyces cerevisiae improves its tolerance to salt, heat and cold stresses. Taken together our results converge to show that the ABA-WDS domain is an intrinsically disordered functional domain whose conformational plasticity could be instrumental to support the versatile functions attributed to the ASR family, including its role in abiotic stress tolerance. Finally, and after validation in the plant system, this domain could be used to improve crop tolerance to abiotic stresses.

Predicting human RNA quadruplex helicases through comparative sequence approaches and helicase mRNA interactome analyses

RNA quadruplexes are non-canonical nucleic acid structures involved in several human disease states and are regulated by a specific subset of RNA helicases. Given the difficulty in identifying RNA quadruplex helicases due to the multifunctionality of these enzymes, we sought to provide a comprehensive in silico analysis of features found in validated RNA quadruplex helicases to predict novel human RNA quadruplex helicases. Using the 64 human RNA helicases, we correlated their amino acid compositions with subsets of RNA quadruplex helicases categorized by varying levels of evidence of RNA quadruplex interaction. Utilizing phylogenetic and synonymous/non-synonymous substitution analyses, we identified an evolutionarily conserved pattern involving predicted intrinsic disorder and a previously identified motif. We analyzed available next-generation sequencing data to determine which RNA helicases directly interacted with predicted RNA quadruplex regions intracellularly and elucidated the relationship with miRNA binding sites adjacent to RNA quadruplexes. Finally, we performed a phylogenetic analysis of all 64 human RNA helicases to establish how RNA quadruplex detection and unwinding activity may be conserved among helicase subfamilies. This work furthers the understanding of commonalities between RNA quadruplex helicases and provides support for the future validation of several human RNA helicases.

Mobility and disorder in antibody and antigen binding sites do not prevent immunochemical recognition

The known polyspecificity of antibodies, which is crucial for efficient immune response, is determined by the conformational flexibility and intrinsic disorder encoded in local peculiarities of the amino acid sequence of antibodies within or in the vicinity of their complementarity determining regions. Similarly, epitopes represent fuzzy binding sites, which are also characterized by local structural flexibility. Existing data suggest that the efficient interactions between antigens and antibodies rely on the conformational mobility and some disorder of their binding sites and therefore can be relatively well described by the "flexible lock - adjustable key" model, whereas both, extreme order (rigid lock-and-key) and extreme disorder (viral shape-shifters) are not compatible with the efficient antigen-antibody interactions and are not present in immune interactions.

The cytoplasmic tail of influenza A/H1N1 virus hemagglutinin is β-structural

Influenza A/H1N1 virus hemagglutinin (HA) is an integral type I glycoprotein that contains a large glycosylated ectodomain, a transmembrane domain, and a cytoplasmic tail (CT) of 10-14 amino acid residues. There are absolutely no data on the secondary or tertiary structure of the HA CT, which is important for virus pathogenesis. Three highly conserved cysteines are post-translationally modified by the attachment of fatty acid residues that pin the CT to the lipid membrane inside the virion. We applied circular dichroism (CD) and fluorescence spectroscopy analysis to examine four synthetic peptides corresponding to 14-15 C-terminal residues of H1 subtype HA (NH₂-WMCSNGSLQCRICI-COOH; NH₂-FWMCSNGSLQCRICI-COOH), with free or acetaminomethylated cysteines, in the reduced or non-reduced state, at various pH values and temperatures. The CD analysis detected the formation of a β-structure (30-65% according to the new BeStSel algorithm), in addition to an unstructured random coil, in every peptide in various conditions. It was completely or partially recognized as an antiparallel β-structure that was also confirmed by the multi-bounce Horizontal Attenuated Total Reflectance Fourier Transformed Infrared (HATR-FTIR) spectroscopy analysis. According to the experimental data, as well as 3 D modeling, we assume that the amino acid sequence corresponding to the HA CT may form a short antiparallel β-structure under the lipid membrane within a virion.Communicated by Ramaswamy H. Sarma.