Effective concentrations enforced by intrinsically disordered linkers are governed by polymer physics

Conserved molecular recognition by an intrinsically disordered region in the absence of sequence conservation

Intrinsically disordered regions (IDRs) are critical for cellular function yet often appear to lack sequence conservation when assessed by multiple sequence alignments. This raises the question of if and how function can be encoded and preserved in these regions despite massive sequence variation. To address this question, we have applied coarse-grained molecular dynamics simulations to investigate non-specific RNA binding of coronavirus nucleocapsid proteins. Coronavirus nucleocapsid proteins consist of multiple interspersed disordered and folded domains that bind RNA. Here, we focus on the first two domains of coronavirus nucleocapsid proteins: the disordered N-terminal domain (NTD) and the folded RNA binding domain (RBD). While the NTD is highly variable across evolution, the RBD is structurally conserved. This combination makes the NTD-RBD a convenient model system for exploring the interplay between an IDR adjacent to a folded domain and how changes in IDR sequence can influence molecular recognition of a partner. Our results reveal a surprising degree of sequence-specificity encoded by both the composition and the precise order of the amino acids in the NTD. The presence of an NTD can - depending on the sequence - either suppress or enhance RNA binding. Despite this sensitivity, large-scale variation in NTD sequences is possible while certain sequence features are retained. Consequently, a conformationally-conserved dynamic and disordered RNA:protein complex is found across nucleocapsid protein orthologs despite large-scale changes in both NTD sequence and RBD surface chemistry. Taken together, these insights shed light on the ability of disordered regions to preserve functional characteristics despite their sequence variability.

HYPK: A marginally disordered protein sensitive to charge decoration

Intrinsically disordered proteins (IDPs) that lie close to the empirical boundary separating IDPs and folded proteins in Uversky's charge-hydropathy plot may behave as "marginal IDPs" and sensitively switch conformation upon changes in environment (temperature, crowding, and charge screening), sequence, or both. In our search for such a marginal IDP, we selected Huntingtin-interacting protein K (HYPK) near that boundary as a candidate; PKIα, also near that boundary, has lower secondary structure propensity; and Crk1, just across the boundary on the folded side, has higher secondary structure propensity. We used a qualitative Förster resonance energy transfer-based assay together with circular dichroism to simultaneously probe global and local conformation. HYPK shows several unique features indicating marginality: a cooperative transition in end-to-end distance with temperature, like Crk1 and folded proteins, but unlike PKIα; enhanced secondary structure upon crowding, in contrast to Crk1 and PKIα; and a cross-over from salt-induced expansion to compaction at high temperature, likely due to a structure-to-disorder transition not seen in Crk1 and PKIα. We then tested HYPK's sensitivity to charge patterning by designing charge-flipped variants including two specific sequences with identical amino acid composition that markedly differ in their predicted size and response to salt. The experimentally observed trends, also including mutants of PKIα, verify the predictions from sequence charge decoration metrics. Marginal proteins like HYPK show features of both folded and disordered proteins that make them sensitive to physicochemical perturbations and structural control by charge patterning.

Labile assembly of a tardigrade protein induces biostasis

Tardigrades are microscopic animals that survive desiccation by inducing biostasis. To survive drying tardigrades rely on intrinsically disordered CAHS proteins, which also function to prevent perturbations induced by drying in vitro and in heterologous systems. CAHS proteins have been shown to form gels both in vitro and in vivo, which has been speculated to be linked to their protective capacity. However, the sequence features and mechanisms underlying gel formation and the necessity of gelation for protection have not been demonstrated. Here we report a mechanism of fibrillization and gelation for CAHS D similar to that of intermediate filament assembly. We show that in vitro, gelation restricts molecular motion, immobilizing and protecting labile material from the harmful effects of drying. In vivo, we observe that CAHS D forms fibrillar networks during osmotic stress. Fibrillar networking of CAHS D improves survival of osmotically shocked cells. We observe two emergent properties associated with fibrillization; (i) prevention of cell volume change and (ii) reduction of metabolic activity during osmotic shock. We find that there is no significant correlation between maintenance of cell volume and survival, while there is a significant correlation between reduced metabolism and survival. Importantly, CAHS D's fibrillar network formation is reversible and metabolic rates return to control levels after CAHS fibers are resolved. This work provides insights into how tardigrades induce reversible biostasis through the self-assembly of labile CAHS gels.

Design and Thermodynamics Principles to Program the Cooperativity of Molecular Assemblies

Most functional nanosystems in living organisms are constructed using multimeric assemblies that provide multiple advantages over their monomeric counterparts such as cooperative or anti-cooperative responses, integration of multiple signals and self-regulation. Inspired by these natural nanosystems, chemists have been synthesizing self-assembled supramolecular systems over the last 50 years with increasing complexity with applications ranging from biosensing, drug delivery, synthetic biology, and system chemistry. Although many advances have been made concerning the design principles of novel molecular architectures and chemistries, little is still known, however, about how to program their dynamic of assembly so that they can assemble at the required concentration and with the right sensitivity. Here, we used synthetic DNA assemblies and double-mutant cycle analysis to explore the thermodynamic basis to program the cooperativity of molecular assemblies. The results presented here exemplify how programmable molecular assemblies can be efficiently built by fusing interacting domains and optimizing their compaction. They may also provide the rational basis for understanding the thermodynamic and mechanistic principles driving the evolution of multimeric biological complexes.

Quantitative description of the phase-separation behavior of the multivalent SLP65-CIN85 complex

Biomolecular condensates play a major role in cell compartmentalization, besides membrane-enclosed organelles. The multivalent SLP65 and CIN85 proteins are proximal B-cell antigen receptor (BCR) signal effectors and critical for proper immune responses. In association with intracellular vesicles, the two effector proteins form phase separated condensates prior to antigen stimulation, thereby preparing B lymphocytes for rapid and effective activation upon BCR ligation. Within this tripartite system, 6 proline-rich motifs (PRMs) of SLP65 interact promiscuously with 3 SH3 domains of the CIN85 monomer, establishing 18 individual SH3-PRM interactions whose individual dissociation constants we determined. Based on these 18 dissociation constants, we measured the phase-separation properties of the natural SLP65/CIN85 system as well as designer constructs that emphasize the strongest SH3/PRM interactions. By modeling these various SLP65/CIN85 constructs with the program LASSI (LAttice simulation engine for Sticker and Spacer Interactions), we reproduced the observed phase-separation properties. In addition, LASSI revealed a deviation in the experimental measurement, which was independently identified as a previously unknown intramolecular interaction. Thus, thermodynamic properties of the individual PRM/SH3 interactions allow us to model the phase-separation behavior of the SLP65/CIN85 system faithfully.

Direct prediction of intrinsically disordered protein conformational properties from sequence

Intrinsically disordered regions (IDRs) are ubiquitous across all domains of life and play a range of functional roles. While folded domains are generally well described by a stable three-dimensional structure, IDRs exist in a collection of interconverting states known as an ensemble. This structural heterogeneity means that IDRs are largely absent from the Protein Data Bank, contributing to a lack of computational approaches to predict ensemble conformational properties from sequence. Here we combine rational sequence design, large-scale molecular simulations and deep learning to develop ALBATROSS, a deep-learning model for predicting ensemble dimensions of IDRs, including the radius of gyration, end-to-end distance, polymer-scaling exponent and ensemble asphericity, directly from sequences at a proteome-wide scale. ALBATROSS is lightweight, easy to use and accessible as both a locally installable software package and a point-and-click-style interface via Google Colab notebooks. We first demonstrate the applicability of our predictors by examining the generalizability of sequence-ensemble relationships in IDRs. Then, we leverage the high-throughput nature of ALBATROSS to characterize the sequence-specific biophysical behavior of IDRs within and between proteomes.

The molecular basis for cellular function of intrinsically disordered protein regions

Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.

Checkpoint activation by Spd1: a competition-based system relying on tandem disordered PCNA binding motifs

DNA regulation, replication and repair are processes fundamental to all known organisms and the sliding clamp proliferating cell nuclear antigen (PCNA) is central to all these processes. S-phase delaying protein 1 (Spd1) from S. pombe, an intrinsically disordered protein that causes checkpoint activation by inhibiting the enzyme ribonucleotide reductase, has one of the most divergent PCNA binding motifs known. Using NMR spectroscopy, in vivo assays, X-ray crystallography, calorimetry, and Monte Carlo simulations, an additional PCNA binding motif in Spd1, a PIP-box, is revealed. The two tandemly positioned, low affinity sites exchange rapidly on PCNA exploiting the same binding sites. Increasing or decreasing the binding affinity between Spd1 and PCNA through mutations of either motif compromised the ability of Spd1 to cause checkpoint activation in yeast. These results pinpoint a role for PCNA in Spd1-mediated checkpoint activation and suggest that its tandemly positioned short linear motifs create a neatly balanced competition-based system, involving PCNA, Spd1 and the small ribonucleotide reductase subunit, Suc22R2. Similar mechanisms may be relevant in other PCNA binding ligands where divergent binding motifs so far have gone under the PIP-box radar.

Bioluminescent detection of viral surface proteins using branched multivalent protein switches

Fast and reliable virus diagnostics is key to prevent the spread of viruses in populations. A hallmark of viruses is the presence of multivalent surface proteins, a property that can be harnessed to control conformational switching in sensor proteins. Here, we introduce a new sensor platform (dark-LUX) for the detection of viral surface proteins consisting of a general bioluminescent framework that can be post-translationally functionalized with separately expressed binding domains. The platform relies on (1) plug-and-play bioconjugation of different binding proteins via SpyTag/SpyCatcher technology to create branched protein structures, (2) an optimized turn-on bioluminescent switch based on complementation of the split-luciferase NanoBiT upon target binding and (3) straightforward exploration of the protein linker space. The influenza A virus (IAV) surface proteins hemagglutinin (HA) and neuraminidase (NA) were used as relevant multivalent targets to establish proof of principle and optimize relevant parameters such as linker properties, choice of target binding domains and the optimal combination of the competing NanoBiT components SmBiT and DarkBiT. The sensor framework allows rapid conjugation and exchange of various binding domains including scFvs, nanobodies and de novo designed binders for a variety of targets, including the construction of a heterobivalent switch that targets the head and stem region of hemagglutinin. The modularity of the platform thus allows straightforward optimization of binding domains and scaffold properties for existing viral targets, and is well suited to quickly adapt bioluminescent sensor proteins to effectively detect newly evolving viral epitopes.

Heteromultivalency enables enhanced detection of nucleic acid mutations

Detecting genetic mutations such as single nucleotide polymorphisms (SNPs) is necessary to prescribe effective cancer therapies, perform genetic analyses and distinguish similar viral strains. Traditionally, SNP sensing uses short oligonucleotide probes that differentially bind the SNP and wild-type targets. However, DNA hybridization-based techniques require precise tuning of the probe's binding affinity to manage the inherent trade-off between specificity and sensitivity. As conventional hybridization offers limited control over binding affinity, here we generate heteromultivalent DNA-functionalized particles and demonstrate optimized hybridization specificity for targets containing one or two mutations. By investigating the role of oligo lengths, spacer lengths and binding orientation, we reveal that heteromultivalent hybridization enables fine-tuned specificity for a single SNP and dramatic enhancements in specificity for two non-proximal SNPs empowered by highly cooperative binding. Capitalizing on these abilities, we demonstrate straightforward discrimination between heterozygous cis and trans mutations and between different strains of the SARS-CoV-2 virus. Our findings indicate that heteromultivalent hybridization offers substantial improvements over conventional monovalent hybridization-based methods.

Conformational ensembles of the human intrinsically disordered proteome

Intrinsically disordered proteins and regions (collectively, IDRs) are pervasive across proteomes in all kingdoms of life, help to shape biological functions and are involved in numerous diseases. IDRs populate a diverse set of transiently formed structures and defy conventional sequence-structure-function relationships1. Developments in protein science have made it possible to predict the three-dimensional structures of folded proteins at the proteome scale2. By contrast, there is a lack of knowledge about the conformational properties of IDRs, partly because the sequences of disordered proteins are poorly conserved and also because only a few of these proteins have been characterized experimentally. The inability to predict structural properties of IDRs across the proteome has limited our understanding of the functional roles of IDRs and how evolution shapes them. As a supplement to previous structural studies of individual IDRs3, we developed an efficient molecular model to generate conformational ensembles of IDRs and thereby to predict their conformational properties from sequences4,5. Here we use this model to simulate nearly all of the IDRs in the human proteome. Examining conformational ensembles of 28,058 IDRs, we show how chain compaction is correlated with cellular function and localization. We provide insights into how sequence features relate to chain compaction and, using a machine-learning model trained on our simulation data, show the conservation of conformational properties across orthologues. Our results recapitulate observations from previous studies of individual protein systems and exemplify how to link-at the proteome scale-conformational ensembles with cellular function and localization, amino acid sequence, evolutionary conservation and disease variants. Our freely available database of conformational properties will encourage further experimental investigation and enable the generation of hypotheses about the biological roles and evolution of IDRs.

Multivalency emerges as a common feature of intrinsically disordered protein interactions

Intrinsically disordered proteins (IDPs) use their unique molecular properties and conformational plasticity to interact with cellular partners in a wide variety of biological contexts. Multivalency is an important feature of IDPs that allows for utilization of an expanded toolkit for interactions with other macromolecules and confers additional complexity to molecular recognition processes. Recent studies have offered insights into how multivalent interactions of IDPs enable responsive and sensitive regulation in the context of transcription and cellular signaling. Multivalency is also widely recognized as an important feature of IDP interactions that mediate formation of biomolecular condensates. We highlight recent examples of multivalent interactions of IDPs across diverse contexts to illustrate the breadth of biological processes that utilize multivalency in molecular interactions.

Design of functional intrinsically disordered proteins

Many proteins do not fold into a fixed three-dimensional structure, but rather function in a highly disordered state. These intrinsically disordered proteins pose a unique challenge to protein engineering and design: How can proteins be designed de novo if not by tailoring their structure? Here, we will review the nascent field of design of intrinsically disordered proteins with focus on applications in biotechnology and medicine. The design goals should not necessarily be the same as for de novo design of folded proteins as disordered proteins have unique functional strengths and limitations. We focus on functions where intrinsically disordered proteins are uniquely suited including disordered linkers, desiccation chaperones, sensors of the chemical environment, delivery of pharmaceuticals, and constituents of biomolecular condensates. Design of functional intrinsically disordered proteins relies on a combination of computational tools and heuristics gleaned from sequence-function studies. There are few cases where intrinsically disordered proteins have made it into industrial applications. However, we argue that disordered proteins can perform many roles currently performed by organic polymers, and that these proteins might be more designable due to their modularity.

Molecular switching in transcription through splicing and proline-isomerization regulates stress responses in plants

The Arabidopsis thaliana DREB2A transcription factor interacts with the negative regulator RCD1 and the ACID domain of subunit 25 of the transcriptional co-regulator mediator (Med25) to integrate stress signals for gene expression, with elusive molecular interplay. Using biophysical and structural analyses together with high-throughput screening, we reveal a bivalent binding switch in DREB2A containing an ACID-binding motif (ABS) and the known RCD1-binding motif (RIM). The RIM is lacking in a stress-induced DREB2A splice variant with retained transcriptional activity. ABS and RIM bind to separate sites on Med25-ACID, and NMR analyses show a structurally heterogeneous complex deriving from a DREB2A-ABS proline residue populating cis- and trans-isomers with remote impact on the RIM. The cis-isomer stabilizes an α-helix, while the trans-isomer may introduce energetic frustration facilitating rapid exchange between activators and repressors. Thus, DREB2A uses a post-transcriptionally and post-translationally modulated switch for transcriptional regulation.

Autoinhibition in the Signal Transducer CIN85 Modulates B Cell Activation

Signal transduction by the ligated B cell antigen receptor (BCR) depends on the preorganization of its intracellular components, such as the effector proteins SLP65 and CIN85 within phase-separated condensates. These liquid-like condensates are based on the interaction between three Src homology 3 (SH3) domains and the corresponding proline-rich recognition motifs (PRM) in CIN85 and SLP65, respectively. However, detailed information on the protein conformation and how it impacts the capability of SLP65/CIN85 condensates to orchestrate BCR signal transduction is still lacking. This study identifies a hitherto unknown intramolecular SH3:PRM interaction between the C-terminal SH3 domain (SH3C) of CIN85 and an adjacent PRM. We used high-resolution nuclear magnetic resonance (NMR) experiments to study the flexible linker region containing the PRM and determined the extent of the interaction in multidomain constructs of the protein. Moreover, we observed that the phosphorylation of a serine residue located in the immediate vicinity of the PRM regulates this intramolecular interaction. This allows for a dynamic modulation of CIN85's valency toward SLP65. B cell culture experiments further revealed that the PRM/SH3C interaction is crucial for maintaining the physiological level of SLP65/CIN85 condensate formation, activation-induced membrane recruitment of CIN85, and subsequent mobilization of Ca2+. Our findings therefore suggest that the intramolecular interaction with the adjacent disordered linker is effective in modulating CIN85's valency both in vitro and in vivo. This therefore constitutes a powerful way for the modulation of SLP65/CIN85 condensate formation and subsequent B cell signaling processes within the cell.

DisoFLAG: accurate prediction of protein intrinsic disorder and its functions using graph-based interaction protein language model

Intrinsically disordered proteins and regions (IDPs/IDRs) are functionally important proteins and regions that exist as highly dynamic conformations under natural physiological conditions. IDPs/IDRs exhibit a broad range of molecular functions, and their functions involve binding interactions with partners and remaining native structural flexibility. The rapid increase in the number of proteins in sequence databases and the diversity of disordered functions challenge existing computational methods for predicting protein intrinsic disorder and disordered functions. A disordered region interacts with different partners to perform multiple functions, and these disordered functions exhibit different dependencies and correlations. In this study, we introduce DisoFLAG, a computational method that leverages a graph-based interaction protein language model (GiPLM) for jointly predicting disorder and its multiple potential functions. GiPLM integrates protein semantic information based on pre-trained protein language models into graph-based interaction units to enhance the correlation of the semantic representation of multiple disordered functions. The DisoFLAG predictor takes amino acid sequences as the only inputs and provides predictions of intrinsic disorder and six disordered functions for proteins, including protein-binding, DNA-binding, RNA-binding, ion-binding, lipid-binding, and flexible linker. We evaluated the predictive performance of DisoFLAG following the Critical Assessment of protein Intrinsic Disorder (CAID) experiments, and the results demonstrated that DisoFLAG offers accurate and comprehensive predictions of disordered functions, extending the current coverage of computationally predicted disordered function categories. The standalone package and web server of DisoFLAG have been established to provide accurate prediction tools for intrinsic disorders and their associated functions.

Activation of caspase-9 on the apoptosome as studied by methyl-TROSY NMR

Mitochondrial apoptotic signaling cascades lead to the formation of the apoptosome, a 1.1-MDa heptameric protein scaffold that recruits and activates the caspase-9 protease. Once activated, caspase-9 cleaves and activates downstream effector caspases, triggering the onset of cell death through caspase-mediated proteolysis of cellular proteins. Failure to activate caspase-9 enables the evasion of programmed cell death, which occurs in various forms of cancer. Despite the critical apoptotic function of caspase-9, the structural mechanism by which it is activated on the apoptosome has remained elusive. Here, we used a combination of methyl-transverse relaxation-optimized NMR spectroscopy, protein engineering, and biochemical assays to study the activation of caspase-9 bound to the apoptosome. In the absence of peptide substrate, we observed that both caspase-9 and its isolated protease domain (PD) only very weakly dimerize with dissociation constants in the millimolar range. Methyl-NMR spectra of isotope-labeled caspase-9, within the 1.3-MDa native apoptosome complex or an engineered 480-kDa apoptosome mimic, reveal that the caspase-9 PD remains monomeric after recruitment to the scaffold. Binding to the apoptosome, therefore, organizes caspase-9 PDs so that they can rapidly and extensively dimerize only when substrate is present, providing an important layer in the regulation of caspase-9 activation. Our work highlights the unique role of NMR spectroscopy to structurally characterize protein domains that are flexibly tethered to large scaffolds, even in cases where the molecular targets are in excess of 1 MDa, as in the present example.

Intrinsically disordered regions are poised to act as sensors of cellular chemistry

Intrinsically disordered proteins and protein regions (IDRs) are abundant in eukaryotic proteomes and play a wide variety of essential roles. Instead of folding into a stable structure, IDRs exist in an ensemble of interconverting conformations whose structure is biased by sequence-dependent interactions. The absence of a stable 3D structure, combined with high solvent accessibility, means that IDR conformational biases are inherently sensitive to changes in their environment. Here, we argue that IDRs are ideally poised to act as sensors and actuators of cellular physicochemistry. We review the physical principles that underlie IDR sensitivity, the molecular mechanisms that translate this sensitivity to function, and recent studies where environmental sensing by IDRs may play a key role in their downstream function.

Sequence-ensemble-function relationships for disordered proteins in live cells

Intrinsically disordered protein regions (IDRs) are ubiquitous across all kingdoms of life and play a variety of essential cellular roles. IDRs exist in a collection of structurally distinct conformers known as an ensemble. An IDR's amino acid sequence determines its ensemble, which in turn can play an important role in dictating molecular function. Yet a clear link connecting IDR sequence, its ensemble properties, and its molecular function in living cells has not been directly established. Here, we set out to test this sequence-ensemble-function paradigm using a novel computational method (GOOSE) that enables the rational design of libraries of IDRs by systematically varying specific sequence properties. Using ensemble FRET, we measured the ensemble dimensions of a library of rationally designed IDRs in human-derived cell lines, revealing how IDR sequence influences ensemble dimensions in situ. Furthermore, we show that the interplay between sequence and ensemble can tune an IDR's ability to sense changes in cell volume - a de novo molecular function for these synthetic sequences. Our results establish biophysical rules for intracellular sequence-ensemble relationships, enable a new route for understanding how IDR sequences map to function in live cells, and set the ground for the design of synthetic IDRs with de novo function.

Sequence-ensemble-function relationships for disordered proteins in live cells

Intrinsically disordered protein regions (IDRs) are ubiquitous across all kingdoms of life and play a variety of essential cellular roles. IDRs exist in a collection of structurally distinct conformers known as an ensemble. IDR amino acid sequence determines its ensemble, which in turn can play an important role in dictating molecular function. Yet a clear link connecting IDR sequence, its ensemble properties, and its molecular function in living cells has not been systematically established. Here, we set out to test this sequence-ensemble-function paradigm using a novel computational method (GOOSE) that enables the rational design of libraries of IDRs by systematically varying specific sequence properties. Using ensemble FRET, we measured the ensemble dimensions of a library of rationally designed IDRs in human-derived cell lines, revealing how IDR sequence influences ensemble dimensions in situ. Furthermore, we show that the interplay between sequence and ensemble can tune an IDR's ability to sense changes in cell volume - a de novomolecular function for these synthetic sequences. Our results establish biophysical rules for intracellular sequence-ensemble relationships, enable a new route for understanding how IDR sequences map to function in live cells, and set the ground for the design of synthetic IDRs with de novo function.

The effect of linker conformation on performance and stability of a two-domain lytic polysaccharide monooxygenase

A considerable number of lytic polysaccharide monooxygenases (LPMOs) and other carbohydrate-active enzymes are modular, with catalytic domains being tethered to additional domains, such as carbohydrate-binding modules, by flexible linkers. While such linkers may affect the structure, function, and stability of the enzyme, their roles remain largely enigmatic, as do the reasons for natural variation in length and sequence. Here, we have explored linker functionality using the two-domain cellulose-active ScLPMO10C from Streptomyces coelicolor as a model system. In addition to investigating the WT enzyme, we engineered three linker variants to address the impact of both length and sequence and characterized these using small-angle X-ray scattering, NMR, molecular dynamics simulations, and functional assays. The resulting data revealed that, in the case of ScLPMO10C, linker length is the main determinant of linker conformation and enzyme performance. Both the WT and a serine-rich variant, which have the same linker length, demonstrated better performance compared with those with either a shorter linker or a longer linker. A highlight of our findings was the substantial thermostability observed in the serine-rich variant. Importantly, the linker affects thermal unfolding behavior and enzyme stability. In particular, unfolding studies show that the two domains unfold independently when mixed, whereas the full-length enzyme shows one cooperative unfolding transition, meaning that the impact of linkers in biomass-processing enzymes is more complex than mere structural tethering.

SOURSOP: A Python Package for the Analysis of Simulations of Intrinsically Disordered Proteins

Conformational heterogeneity is a defining hallmark of intrinsically disordered proteins and protein regions (IDRs). The functions of IDRs and the emergent cellular phenotypes they control are associated with sequence-specific conformational ensembles. Simulations of conformational ensembles that are based on atomistic and coarse-grained models are routinely used to uncover the sequence-specific interactions that may contribute to IDR functions. These simulations are performed either independently or in conjunction with data from experiments. Functionally relevant features of IDRs can span a range of length scales. Extracting these features requires analysis routines that quantify a range of properties. Here, we describe a new analysis suite simulation analysis of unfolded regions of proteins (SOURSOP), an object-oriented and open-source toolkit designed for the analysis of simulated conformational ensembles of IDRs. SOURSOP implements several analysis routines motivated by principles in polymer physics, offering a unique collection of simple-to-use functions to characterize IDR ensembles. As an extendable framework, SOURSOP supports the development and implementation of new analysis routines that can be easily packaged and shared.

Conserved molecular recognition by an intrinsically disordered region in the absence of sequence conservation

Intrinsically disordered regions (IDRs) are critical for cellular function, yet often appear to lack sequence conservation when assessed by multiple sequence alignments. This raises the question of if and how function can be encoded and preserved in these regions despite massive sequence variation. To address this question, we have applied coarse-grained molecular dynamics simulations to investigate non-specific RNA binding of coronavirus nucleocapsid proteins. Coronavirus nucleocapsid proteins consist of multiple interspersed disordered and folded domains that bind RNA. We focussed here on the first two domains of coronavirus nucleocapsid proteins, the disordered N-terminal domain (NTD) followed by the folded RNA binding domain (RBD). While the NTD is highly variable across evolution, the RBD is structurally conserved. This combination makes the NTD-RBD a convenient model system to explore the interplay between an IDR adjacent to a folded domain, and how changes in IDR sequence can influence molecular recognition of a partner. Our results reveal a surprising degree of sequence-specificity encoded by both the composition and the precise order of the amino acids in the NTD. The presence of an NTD can - depending on the sequence - either suppress or enhance RNA binding. Despite this sensitivity, large-scale variation in NTD sequences is possible while certain sequence features are retained. Consequently, a conformationally-conserved fuzzy RNA:protein complex is found across nucleocapsid protein orthologs, despite large-scale changes in both NTD sequence and RBD surface chemistry. Taken together, these insights shed light on the ability of disordered regions to preserve functional characteristics despite their sequence variability.

Novel Design of an α-Amylase with an N-Terminal CBM20 in Aspergillus niger Improves Binding and Processing of a Broad Range of Starches

In the starch processing industry including the food and pharmaceutical industries, α-amylase is an important enzyme that hydrolyses the α-1,4 glycosidic bonds in starch, producing shorter maltooligosaccharides. In plants, starch molecules are organised in granules that are very compact and rigid. The level of starch granule rigidity affects resistance towards enzymatic hydrolysis, resulting in inefficient starch degradation by industrially available α-amylases. In an approach to enhance starch hydrolysis, the domain architecture of a Glycoside Hydrolase (GH) family 13 α-amylase from Aspergillus niger was engineered. In all fungal GH13 α-amylases that carry a carbohydrate binding domain (CBM), these modules are of the CBM20 family and are located at the C-terminus of the α-amylase domain. To explore the role of the domain order, a new GH13 gene encoding an N-terminal CBM20 domain was designed and found to be fully functional. The starch binding capacity and enzymatic activity of N-terminal CBM20 α-amylase was found to be superior to that of native GH13 without CBM20. Based on the kinetic parameters, the engineered N-terminal CBM20 variant displayed surpassing activity rates compared to the C-terminal CBM20 version for the degradation on a wide range of starches, including the more resistant raw potato starch for which it exhibits a two-fold higher Vmax underscoring the potential of domain engineering for these carbohydrate active enzymes.

Effects of linker length on phase separation: lessons from the Rubisco-EPYC1 system of the algal pyrenoid

Biomolecular condensates are membraneless organelles formed via phase separation of macromolecules, typically consisting of bond-forming "stickers" connected by flexible "linkers". Linkers have diverse roles, such as occupying space and facilitating interactions. To understand how linker length relative to other lengths affects condensation, we focus on the pyrenoid, which enhances photosynthesis in green algae. Specifically, we apply coarse-grained simulations and analytical theory to the pyrenoid proteins of Chlamydomonas reinhardtii: the rigid holoenzyme Rubisco and its flexible partner EPYC1. Remarkably, halving EPYC1 linker lengths decreases critical concentrations by ten-fold. We attribute this difference to the molecular "fit" between EPYC1 and Rubisco. Varying Rubisco sticker locations reveals that the native sites yield the poorest fit, thus optimizing phase separation. Surprisingly, shorter linkers mediate a transition to a gas of rods as Rubisco stickers approach the poles. These findings illustrate how intrinsically disordered proteins affect phase separation through the interplay of molecular length scales.

The Analytical Flory Random Coil Is a Simple-to-Use Reference Model for Unfolded and Disordered Proteins

Denatured, unfolded, and intrinsically disordered proteins (collectively referred to here as unfolded proteins) can be described using analytical polymer models. These models capture various polymeric properties and can be fit to simulation results or experimental data. However, the model parameters commonly require users' decisions, making them useful for data interpretation but less clearly applicable as stand-alone reference models. Here we use all-atom simulations of polypeptides in conjunction with polymer scaling theory to parameterize an analytical model of unfolded polypeptides that behave as ideal chains (ν = 0.50). The model, which we call the analytical Flory random coil (AFRC), requires only the amino acid sequence as input and provides direct access to probability distributions of global and local conformational order parameters. The model defines a specific reference state to which experimental and computational results can be compared and normalized. As a proof-of-concept, we use the AFRC to identify sequence-specific intramolecular interactions in simulations of disordered proteins. We also use the AFRC to contextualize a curated set of 145 different radii of gyration obtained from previously published small-angle X-ray scattering experiments of disordered proteins. The AFRC is implemented as a stand-alone software package and is also available via a Google Colab notebook. In summary, the AFRC provides a simple-to-use reference polymer model that can guide intuition and aid in interpreting experimental or simulation results.

Structural Preferences Shape the Entropic Force of Disordered Protein Ensembles

Intrinsically disordered protein regions (IDRs) make up over 30% of the human proteome and exist in a dynamic conformational ensemble instead of a native, well-folded structure. Tethering IDRs to a surface (for example, the surface of a well-folded region of the same protein) can reduce the number of accessible conformations in these ensembles. This reduces the ensemble's conformational entropy, generating an effective entropic force that pulls away from the point of tethering. Recent experimental work has shown that this entropic force causes measurable, physiologically relevant changes to protein function. But how the magnitude of this force depends on IDR sequence remains unexplored. Here, we use all-atom simulations to analyze how structural preferences in IDR ensembles contribute to the entropic force they exert upon tethering. We show that sequence-encoded structural preferences play an important role in determining the magnitude of this force: compact, spherical ensembles generate an entropic force that can be several times higher than more extended ensembles. We further show that changes in the surrounding solution's chemistry can modulate the IDR entropic force strength. We propose that the entropic force is a sequence-dependent, environmentally tunable property of terminal IDR sequences.

TransDFL: Identification of Disordered Flexible Linkers in Proteins by Transfer Learning

Disordered flexible linkers (DFLs) are the functional disordered regions in proteins, which are the sub-regions of intrinsically disordered regions (IDRs) and play important roles in connecting domains and maintaining inter-domain interactions. Trained with the limited available DFLs, the existing DFL predictors based on the machine learning techniques tend to predict the ordered residues as DFLs, leading to a high falsepositive rate (FPR) and low prediction accuracy. Previous studies have shown that DFLs are extremely flexible disordered regions, which are usually predicted as disordered residues with high confidence [P(D) > 0.9] by an IDR predictor. Therefore, transferring an IDR predictor to an accurate DFL predictor is of great significance for understanding the functions of IDRs. In this study, we proposed a new predictor called TransDFL for identifying DFLs by transferring the RFPR-IDP predictor for IDR identification to the DFL prediction. The RFPR-IDP was pre-trained with IDR sequences to learn the general features between IDRs and DFLs, which is helpful to reduce the false positives in the ordered regions. RFPR-IDP was fine-tuned with the DFL sequences to capture the specific features of DFLs so as to be transferred into the TransDFL. Experimental results of two application scenarios (prediction of DFLs only in IDRs or prediction of DFLs in entire proteins) showed that TransDFL consistently outperformed other existing DFL predictors with higher accuracy. The corresponding web server of TransDFL can be freely accessed at http://bliulab.net/TransDFL/.

Protein Binding: A Fuzzy Concept

Our understanding of protein binding interactions has matured significantly over the last few years, largely as a result of trying to make sense of the binding interactions of intrinsically disordered proteins. Here, we bring together some disparate ideas that have largely developed independently, and show that they can be linked into a coherent picture that provides insight into quantitative aspects of protein interactions, in particular that transient protein interactions are often optimised for speed, rather than tight binding.

The analytical Flory random coil is a simple-to-use reference model for unfolded and disordered proteins

Denatured, unfolded, and intrinsically disordered proteins (collectively referred to here as unfolded proteins) can be described using analytical polymer models. These models capture various polymeric properties and can be fit to simulation results or experimental data. However, the model parameters commonly require users' decisions, making them useful for data interpretation but less clearly applicable as stand-alone reference models. Here we use all-atom simulations of polypeptides in conjunction with polymer scaling theory to parameterize an analytical model of unfolded polypeptides that behave as ideal chains (ν = 0.50). The model, which we call the analytical Flory Random Coil (AFRC), requires only the amino acid sequence as input and provides direct access to probability distributions of global and local conformational order parameters. The model defines a specific reference state to which experimental and computational results can be compared and normalized. As a proof-of-concept, we use the AFRC to identify sequence-specific intramolecular interactions in simulations of disordered proteins. We also use the AFRC to contextualize a curated set of 145 different radii of gyration obtained from previously published small-angle X-ray scattering experiments of disordered proteins. The AFRC is implemented as a stand-alone software package and is also available via a Google colab notebook. In summary, the AFRC provides a simple-to-use reference polymer model that can guide intuition and aid in interpreting experimental or simulation results.

Parkin and mitochondrial signalling

Aging, toxic chemicals and changes to the cellular environment are sources of oxidative damage to mitochondria which contribute to neurodegenerative conditions including Parkinson's disease. To counteract this, cells have developed signalling mechanisms to identify and remove select proteins and unhealthy mitochondria to maintain homeostasis. Two important proteins that work in concert to control mitochondrial damage are the protein kinase PINK1 and the E3 ligase parkin. In response to oxidative stress, PINK1 phosphorylates ubiquitin present on proteins at the mitochondrial surface. This signals the translocation of parkin, accelerates further phosphorylation, and stimulates ubiquitination of outer mitochondrial membrane proteins such as Miro1/2 and Mfn1/2. The ubiquitination of these proteins is the key step needed to target them for degradation via the 26S proteasomal machinery or eliminate the entire organelle through mitophagy. This review highlights the signalling mechanisms used by PINK1 and parkin and presents several outstanding questions yet to be resolved.

SOURSOP: A Python package for the analysis of simulations of intrinsically disordered proteins

Conformational heterogeneity is a defining hallmark of intrinsically disordered proteins and protein regions (IDRs). The functions of IDRs and the emergent cellular phenotypes they control are associated with sequence-specific conformational ensembles. Simulations of conformational ensembles that are based on atomistic and coarse-grained models are routinely used to uncover the sequence-specific interactions that may contribute to IDR functions. These simulations are performed either independently or in conjunction with data from experiments. Functionally relevant features of IDRs can span a range of length scales. Extracting these features requires analysis routines that quantify a range of properties. Here, we describe a new analysis suite SOURSOP, an object-oriented and open-source toolkit designed for the analysis of simulated conformational ensembles of IDRs. SOURSOP implements several analysis routines motivated by principles in polymer physics, offering a unique collection of simple-to-use functions to characterize IDR ensembles. As an extendable framework, SOURSOP supports the development and implementation of new analysis routines that can be easily packaged and shared.

Structural preferences shape the entropic force of disordered protein ensembles

Intrinsically disordered protein regions (IDRs) make up over 30% of the human proteome and instead of a native, well-folded structure exist in a dynamic conformational ensemble. Tethering IDRs to a surface (for example, the surface of a well-folded region of the same protein) can reduce the number of accessible conformations in IDR ensembles. This reduces the ensemble's conformational entropy, generating an effective entropic force that pulls away from the point of tethering. Recent experimental work has shown that this entropic force causes measurable, physiologically relevant changes to protein function, but how the magnitude of this force depends on the IDR sequence remains unexplored. Here we use all-atom simulations to analyze how structural preferences encoded in dozens of IDR ensembles contribute to the entropic force they exert upon tethering. We show that sequence-encoded structural preferences play an important role in determining the magnitude of this force and that compact, spherical ensembles generate an entropic force that can be several times higher than more extended ensembles. We further show that changes in the surrounding solution's chemistry can modulate IDR entropic force strength. We propose that the entropic force is a sequence-dependent, environmentally tunable property of terminal IDR sequences.

Engineering cellulases for conversion of lignocellulosic biomass

Lignocellulosic biomass is a renewable source of energy, chemicals and materials. Many applications of this resource require the depolymerization of one or more of its polymeric constituents. Efficient enzymatic depolymerization of cellulose to glucose by cellulases and accessory enzymes such as lytic polysaccharide monooxygenases is a prerequisite for economically viable exploitation of this biomass. Microbes produce a remarkably diverse range of cellulases, which consist of glycoside hydrolase (GH) catalytic domains and, although not in all cases, substrate-binding carbohydrate-binding modules (CBMs). As enzymes are a considerable cost factor, there is great interest in finding or engineering improved and robust cellulases, with higher activity and stability, easy expression, and minimal product inhibition. This review addresses relevant engineering targets for cellulases, discusses a few notable cellulase engineering studies of the past decades and provides an overview of recent work in the field.

Effects of Sequence Composition, Patterning and Hydrodynamics on the Conformation and Dynamics of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) perform diverse functions in cellular organization, transport and signaling. Unlike the well-defined structures of the classical natively folded proteins, IDPs and IDRs dynamically span large conformational and structural ensembles. This dynamic disorder impedes the study of the relationship between the amino acid sequences of the IDPs and their spatial structures and dynamics, with different experimental techniques often offering seemingly contradictory results. Although experimental and theoretical evidence indicates that some IDP properties can be understood based on their average biophysical properties and amino acid composition, other aspects of IDP function are dictated by the specifics of the amino acid sequence. We investigate the effects of several key variables on the dimensions and the dynamics of IDPs using coarse-grained polymer models. We focus on the sequence "patchiness" informed by the sequence and biophysical properties of different classes of IDPs-and in particular FG nucleoporins of the nuclear pore complex (NPC). We show that the sequence composition and patterning are well reflected in the global conformational variables such as the radius of gyration and hydrodynamic radius, while the end-to-end distance and dynamics are highly sequence-specific. We find that in good solvent conditions highly heterogeneous sequences of IDPs can be well mapped onto averaged minimal polymer models for the purpose of prediction of the IDPs dimensions and dynamic relaxation times. The coarse-grained simulations are in a good agreement with the results of atomistic MD. We discuss the implications of these results for the interpretation of the recent experimental measurements, and for the further applications of mesoscopic models of FG nucleoporins and IDPs more broadly.

Sequence Context and Complex Hofmeister Salt Interactions Dictate Phase Separation Propensity of Resilin-like Polypeptides

Resilin is an elastic material found in insects with exceptional durability, resilience, and extensibility, making it a promising biomaterial for tissue engineering. The monomeric precursor, pro-resilin, undergoes thermo-responsive self-assembly through liquid-liquid phase separation (LLPS). Understanding the molecular details of this assembly process is critical to developing complex biomaterials. The present study investigates the interplay between the solvent, sequence syntax, structure, and dynamics in promoting LLPS of resilin-like-polypeptides (RLPs) derived from domains 1 and 3 of Drosophila melanogaster pro-resilin. NMR, UV-vis, and microscopy data demonstrate that while kosmotropic salts and low pH promote LLPS, the effects of chaotropic salts with increasing pH are more complex. Subtle variations between the repeating amino acid motifs of resilin domain 1 and domain 3 lead to significantly different salt and pH dependence of LLPS, with domain 3 sequence motifs more strongly favoring phase separation under most conditions. These findings provide new insight into the molecular drivers of RLP phase separation and the complex roles of both RLP sequence and solution composition in fine-tuning assembly conditions.

On the impact of carbohydrate-binding modules (CBMs) in lytic polysaccharide monooxygenases (LPMOs)

Lytic polysaccharide monooxygenases (LPMOs) have revolutionized our understanding of how enzymes degrade insoluble polysaccharides. Compared with the substantial knowledge developed on the structure and mode of action of the catalytic LPMO domains, the (multi)modularity of LPMOs has received less attention. The presence of other domains, in particular carbohydrate-binding modules (CBMs), tethered to LPMOs has profound implications for the catalytic performance of the full-length enzymes. In the last few years, studies on LPMO modularity have led to advancements in elucidating how CBMs, other domains, and linker regions influence LPMO structure and function. This mini review summarizes recent literature, with particular focus on comparative truncation studies, to provide an overview of the diversity in LPMO modularity and the functional implications of this diversity.

The biophysics of disordered proteins from the point of view of single-molecule fluorescence spectroscopy

Intrinsically disordered proteins (IDPs) and regions (IDRs) have emerged as key players across many biological functions and diseases. Differently from structured proteins, disordered proteins lack stable structure and are particularly sensitive to changes in the surrounding environment. Investigation of disordered ensembles requires new approaches and concepts for quantifying conformations, dynamics, and interactions. Here, we provide a short description of the fundamental biophysical properties of disordered proteins as understood through the lens of single-molecule fluorescence observations. Single-molecule Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) provides an extensive and versatile toolbox for quantifying the characteristics of conformational distributions and the dynamics of disordered proteins across many different solution conditions, both in vitro and in living cells.

The MDMX Acidic Domain Uses Allovalency to Bind Both p53 and MDMX

Autoinhibition of p53 binding to MDMX requires two short-linear motifs (SLiMs) containing adjacent tryptophan (WW) and tryptophan-phenylalanine (WF) residues. NMR spectroscopy was used to show the WW and WF motifs directly compete for the p53 binding site on MDMX and circular dichroism spectroscopy was used to show the WW motif becomes helical when it is bound to the p53 binding domain (p53BD) of MDMX. Binding studies using isothermal titration calorimetry showed the WW motif is a stronger inhibitor of p53 binding than the WF motif when they are both tethered to p53BD by the natural disordered linker. We also investigated how the WW and WF motifs interact with the DNA binding domain (DBD) of p53. Both motifs bind independently to similar sites on DBD that overlap the DNA binding site. Taken together our work defines a model for complex formation between MDMX and p53 where a pair of disordered SLiMs bind overlapping sites on both proteins.

An intrinsic temporal order of c-JUN N-terminal phosphorylation regulates its activity by orchestrating co-factor recruitment

Protein phosphorylation is a major regulatory mechanism of cellular signalling. The c-JUN proto-oncoprotein is phosphorylated at four residues within its transactivation domain (TAD) by the JNK family kinases, but the functional significance of c-JUN multisite phosphorylation has remained elusive. Here we show that c-JUN phosphorylation by JNK exhibits defined temporal kinetics, with serine63 and serine73 being phosphorylated more rapidly than threonine91 and threonine93. We identify the positioning of the phosphorylation sites relative to the kinase docking motif, and their primary sequence, as the main factors controlling phosphorylation kinetics. Functional analysis reveals three c-JUN phosphorylation states: unphosphorylated c-JUN recruits the MBD3 repressor, serine63/73 doubly-phosphorylated c-JUN binds to the TCF4 co-activator, whereas the fully phosphorylated form disfavours TCF4 binding attenuating JNK signalling. Thus, c-JUN phosphorylation encodes multiple functional states that drive a complex signalling response from a single JNK input.

The disordered MAX N-terminus modulates DNA binding of the transcription factor MYC:MAX

The intrinsically disordered protein MYC belongs to the family of basic helix-loop-helix leucine zipper (bHLH-LZ) transcription factors (TFs). In complex with its cognate binding partner MAX, MYC preferentially binds to E-Box promotor sequences where it controls fundamental cellular processes such as cell cycle progression, metabolism, and apoptosis. Intramolecular regulation of MYC:MAX has not yet been investigated in detail. In this work, we use Nuclear Magnetic Resonance (NMR) spectroscopy to identify and map interactions between the disordered MAX N-terminus and the MYC:MAX DNA binding domain (DBD). We find that this binding event is mainly driven by electrostatic interactions and that it is competitive with DNA binding. Using Nuclear Magnetic resonance (NMR) spectroscopy and Surface Plasmon Resonance (SPR), we demonstrate that the MAX N-terminus serves to accelerate DNA binding kinetics of MYC:MAX and MAX:MAX dimers, while it simultaneously provides specificity for E-Box DNA. We also establish that these effects are further enhanced by Casein Kinase 2-mediated phosphorylation of two serine residues in the MAX N-terminus. Our work provides new insights how bHLH-LZ TFs are regulated by intramolecular interactions between disordered regions and the folded DNA binding domain.

The material properties of a bacterial-derived biomolecular condensate tune biological function in natural and synthetic systems

Intracellular phase separation is emerging as a universal principle for organizing biochemical reactions in time and space. It remains incompletely resolved how biological function is encoded in these assemblies and whether this depends on their material state. The conserved intrinsically disordered protein PopZ forms condensates at the poles of the bacterium Caulobacter crescentus, which in turn orchestrate cell-cycle regulating signaling cascades. Here we show that the material properties of these condensates are determined by a balance between attractive and repulsive forces mediated by a helical oligomerization domain and an expanded disordered region, respectively. A series of PopZ mutants disrupting this balance results in condensates that span the material properties spectrum, from liquid to solid. A narrow range of condensate material properties supports proper cell division, linking emergent properties to organismal fitness. We use these insights to repurpose PopZ as a modular platform for generating tunable synthetic condensates in human cells.

Fuzzy supertertiary interactions within PSD-95 enable ligand binding

The scaffold protein PSD-95 links postsynaptic receptors to sites of presynaptic neurotransmitter release. Flexible linkers between folded domains in PSD-95 enable a dynamic supertertiary structure. Interdomain interactions within the PSG supramodule, formed by PDZ3, SH3, and Guanylate Kinase domains, regulate PSD-95 activity. Here we combined discrete molecular dynamics and single molecule Förster resonance energy transfer (FRET) to characterize the PSG supramodule, with time resolution spanning picoseconds to seconds. We used a FRET network to measure distances in full-length PSD-95 and model the conformational ensemble. We found that PDZ3 samples two conformational basins, which we confirmed with disulfide mapping. To understand effects on activity, we measured binding of the synaptic adhesion protein neuroligin. We found that PSD-95 bound neuroligin well at physiological pH while truncated PDZ3 bound poorly. Our hybrid structural models reveal how the supertertiary context of PDZ3 enables recognition of this critical synaptic ligand.

A context-dependent and disordered ubiquitin-binding motif

Ubiquitin is a small, globular protein that is conjugated to other proteins as a posttranslational event. A palette of small, folded domains recognizes and binds ubiquitin to translate and effectuate this posttranslational signal. Recent computational studies have suggested that protein regions can recognize ubiquitin via a process of folding upon binding. Using peptide binding arrays, bioinformatics, and NMR spectroscopy, we have uncovered a disordered ubiquitin-binding motif that likely remains disordered when bound and thus expands the palette of ubiquitin-binding proteins. We term this motif Disordered Ubiquitin-Binding Motif (DisUBM) and find it to be present in many proteins with known or predicted functions in degradation and transcription. We decompose the determinants of the motif showing it to rely on features of aromatic and negatively charged residues, and less so on distinct sequence positions in line with its disordered nature. We show that the affinity of the motif is low and moldable by the surrounding disordered chain, allowing for an enhanced interaction surface with ubiquitin, whereby the affinity increases ~ tenfold. Further affinity optimization using peptide arrays pushed the affinity into the low micromolar range, but compromised context dependence. Finally, we find that DisUBMs can emerge from unbiased screening of randomized peptide libraries, featuring in de novo cyclic peptides selected to bind ubiquitin chains. We suggest that naturally occurring DisUBMs can recognize ubiquitin as a posttranslational signal to act as affinity enhancers in IDPs that bind to folded and ubiquitylated binding partners.

Conformational buffering underlies functional selection in intrinsically disordered protein regions

Many disordered proteins conserve essential functions in the face of extensive sequence variation, making it challenging to identify the mechanisms responsible for functional selection. Here we identify the molecular mechanism of functional selection for the disordered adenovirus early gene 1A (E1A) protein. E1A competes with host factors to bind the retinoblastoma (Rb) protein, subverting cell cycle regulation. We show that two binding motifs tethered by a hypervariable disordered linker drive picomolar affinity Rb binding and host factor displacement. Compensatory changes in amino acid sequence composition and sequence length lead to conservation of optimal tethering across a large family of E1A linkers. We refer to this compensatory mechanism as conformational buffering. We also detect coevolution of the motifs and linker, which can preserve or eliminate the tethering mechanism. Conformational buffering and motif-linker coevolution explain robust functional encoding within hypervariable disordered linkers and could underlie functional selection of many disordered protein regions.

iFLinkC-X: A Scalable Framework to Assemble Bespoke Genetically Encoded Co-polymeric Linkers of Variable Lengths and Amino Acid Composition

Linker engineering is rapidly gaining prominence as protein engineers and synthetic biologists construct increasingly sophisticated protein assemblies capable of executing complex molecular functions in the context of biosensing, biocatalysis, or biotherapeutics. Depending on the application, the structural and functional requirements imposed on the underlying linkers can differ vastly. At the same time, there is a distinct lack of methods to effectively code linkers at the level of DNA and tailor them to the functional requirements of different fusion proteins. Addressing these limitations, a scalable framework is presented to compose co-polymeric linkers of variable lengths and amino acid composition based on a limited number of linker fragments stored in sequence-verified entry plasmids. The assembly process is exemplified for Pro-rich linkers in the context of a Zn²⁺-responsive dual-readout BRET/FRET sensor while examining how linker composition impacts key functional properties such as ligand affinity, dynamic range, and their ability to separate structurally distinct domains.

The optimal docking strength for reversibly tethered kinases

Many kinases use reversible docking interactions to augment the specificity of their catalytic domains. Such docking interactions are often structurally independent of the catalytic domain, which allow for a flexible combination of modules in evolution and in bioengineering. The affinity of docking interactions spans several orders of magnitude. This led us to ask how the affinity of the docking interaction affects enzymatic activity and how to pick the optimal interaction module to complement a given substrate. Here, we develop equations that predict the optimal binding strength of a kinase docking interaction and validate it using numerical simulations and steady-state phosphorylation kinetics for tethered protein kinase A. We show that a kinase-substrate pair has an optimum docking strength that depends on their enzymatic constants, the tether architecture, the substrate concentration, and the kinetics of the docking interactions. We show that a reversible tether enhances phosphorylation rates most when 1) the docking strength is intermediate, 2) the substrate is nonoptimal, 3) the substrate concentration is low, 4) the docking interaction has rapid exchange kinetics, and 5) the tether optimizes the effective concentration of the intramolecular reaction. This work serves as a framework for interpreting mutations in kinase docking interactions and as a design guide for engineering enzyme scaffolds.

In vitro evolution of uracil glycosylase towards DnaKJ and GroEL binding evolves different misfolded states

Natural evolution is driven by random mutations that improve fitness. In vitro evolution mimics this process, however, on a short time-scale and is driven by the given bait. Here, we used directed in vitro evolution of a random mutant library of Uracil glycosylase (eUNG) displayed on yeast surface to select for binding to chaperones GroEL, DnaK+DnaJ+ATP (DnaKJ) or E.coli cell extract (CE), using binding to the eUNG inhibitor Ugi as probe for native fold. The CE selected population was further divided to Ugi binders (+U) or non-binders (-U). The aim here was to evaluate the sequence space and physical state of the evolved protein binding the different baits. We found that GroEL, DnaKJ and CE-U select and enrich for mutations causing eUNG to misfold, with the three being enriched in mutations in buried and conserved positions, with a tendency to increase positive charge. Still, each selection had its own trajectory, with GroEL and CE-U selecting mutants highly sensitive to protease cleavage while DnaKJ selected partially structured misfolded species with a tendency to refold, making them less sensitive to proteases. More general, our results show that GroEL has a higher tendency to purge promiscuous misfolded protein mutants from the system, while DnaKJ binds misfolding-prone mutant species that are, upon chaperone release, more likely to natively refold. CE-U shares some of the properties of GroEL- and DnaKJ-selected populations, while harboring also unique properties that can be explained by the presence of additional chaperones in CE, such as Trigger factor, HtpG and ClpB.

Interplay of Affinity and Surface Tethering in Protein Recognition

Surface-tethered ligand-receptor complexes are key components in biological signaling and adhesion. They also find increasing utility in single-molecule assays and biotechnological applications. Here, we study the real-time binding kinetics between various surface-immobilized peptide ligands and their unrestrained receptors. A long peptide tether increases the association of ligand-receptor complexes, experimentally proving the fly casting mechanism where the disorder accelerates protein recognition. On the other hand, a short peptide tether enhances the complex dissociation. Notably, the rate constants measured for the same receptor, but under different spatial constraints, are strongly correlated to one another. Furthermore, this correlation can be used to predict how surface tethering on a ligand-receptor complex alters its binding kinetics. Our results have immediate implications in the broad areas of biomolecular recognition, intrinsically disordered proteins, and biosensor technology.