Insights into Coupled Folding and Binding Mechanisms from Kinetic Studies

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

Design of amyloidogenic peptide traps

Segments of proteins with high β-strand propensity can self-associate to form amyloid fibrils implicated in many diseases. We describe a general approach to bind such segments in β-strand and β-hairpin conformations using de novo designed scaffolds that contain deep peptide-binding clefts. The designs bind their cognate peptides in vitro with nanomolar affinities. The crystal structure of a designed protein-peptide complex is close to the design model, and NMR characterization reveals how the peptide-binding cleft is protected in the apo state. We use the approach to design binders to the amyloid-forming proteins transthyretin, tau, serum amyloid A1 and amyloid β1-42 (Aβ42). The Aβ binders block the assembly of Aβ fibrils as effectively as the most potent of the clinically tested antibodies to date and protect cells from toxic Aβ42 species.

Uncovering the Association Mechanism between Two Intrinsically Flexible Proteins

The understanding of protein-protein interaction mechanisms is key to the atomistic description of cell signaling pathways and for the development of new drugs. In this context, the mechanism of intrinsically disordered proteins folding upon binding has attracted attention. The VirB9 C-terminal domain (VirB9Ct) and the VirB7 N-terminal motif (VirB7Nt) associate with VirB10 to form the outer membrane core complex of the Type IV Secretion System injectisome. Despite forming a stable and rigid complex, VirB7Nt behaves as a random coil, while VirB9Ct is intrinsically dynamic in the free state. Here we combined NMR, stopped-flow fluorescence, and computer simulations using structure-based models to characterize the VirB9Ct-VirB7Nt coupled folding and binding mechanism. Qualitative data analysis suggested that VirB9Ct preferentially binds to VirB7Nt by way of a conformational selection mechanism at lower temperatures. However, at higher temperatures, energy barriers between different VirB9Ct conformations are more easily surpassed. Under these conditions the formation of non-native initial encounter complexes may provide alternative pathways toward the native complex conformation. These observations highlight the intimate relationship between folding and binding, calling attention to the fact that the two molecular partners must search for the most favored intramolecular and intermolecular interactions on a rugged and funnelled conformational energy landscape, along which multiple intermediates may lead to the final native state.

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.

Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model

A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning-based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long timescale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We observe that folding-upon-binding predominantly proceeds through a multi-step induced fit mechanism with several intermediates and do not find evidence for the existence of canonical conformational selection pathways. We observe four kinetically separated native-like bound states that interconvert on timescales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.

Effects of Replacing a Central Glycine Residue in GLP-1 on Receptor Affinity and Signaling Profile

Agonists of the glucagon-like peptide-1 receptor (GLP-1R) are used to treat diabetes and obesity. Cryo-EM structures indicate that GLP-1 is completely α-helical when bound to the GLP-1R. The mature form of this hormone, GLP-1(7-36), contains a glycine residue near the center (Gly22). Since glycine has the second-lowest α-helix propensity among the proteinogenic α-amino acid residues, and Gly22 does not appear to make direct contact with the receptor, we were motivated to explore the impact on agonist activity of altering the α-helix propensity at this position. We examined GLP-1 analogues in which Gly22 was replaced with L-Ala, D-Ala, or β-amino acid residues with varying helix propensities. The results suggest that the receptor is reasonably tolerant of variations in helix propensity, and that the functional receptor-agonist complex may comprise a conformational spectrum rather than a single fixed structure.

Phosphoregulation of DNA repair via the Rad51 auxiliary factor Swi5-Sfr1

Homologous recombination (HR) is a major pathway for the repair of DNA double-strand breaks, the most severe form of DNA damage. The Rad51 protein is central to HR, but multiple auxiliary factors regulate its activity. The heterodimeric Swi5-Sfr1 complex is one such factor. It was previously shown that two sites within the intrinsically disordered domain of Sfr1 are critical for the interaction with Rad51. Here, we show that phosphorylation of five residues within this domain regulates the interaction of Swi5-Sfr1 with Rad51. Biochemical reconstitutions demonstrated that a phosphomimetic mutant version of Swi5-Sfr1 is defective in both the physical and functional interaction with Rad51. This translated to a defect in DNA repair, with the phosphomimetic mutant yeast strain phenocopying a previously established interaction mutant. Interestingly, a strain in which Sfr1 phosphorylation was blocked also displayed sensitivity to DNA damage. Taken together, we propose that controlled phosphorylation of Sfr1 is important for the role of Swi5-Sfr1 in promoting Rad51-dependent DNA repair.

Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model

A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long-time scale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We do not, however, find evidence for the existence of canonical conformational selection or induced fit binding pathways. We observe four kinetically separated native-like bound states that interconvert on time scales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.

Functional regulation of an intrinsically disordered protein via a conformationally excited state

A longstanding goal in the field of intrinsically disordered proteins (IDPs) is to characterize their structural heterogeneity and pinpoint the role of this heterogeneity in IDP function. Here, we use multinuclear chemical exchange saturation (CEST) nuclear magnetic resonance to determine the structure of a thermally accessible globally folded excited state in equilibrium with the intrinsically disordered native ensemble of a bacterial transcriptional regulator CytR. We further provide evidence from double resonance CEST experiments that the excited state, which structurally resembles the DNA-bound form of cytidine repressor (CytR), recognizes DNA by means of a "folding-before-binding" conformational selection pathway. The disorder-to-order regulatory switch in DNA recognition by natively disordered CytR therefore operates through a dynamical variant of the lock-and-key mechanism where the structurally complementary conformation is transiently accessed via thermal fluctuations.

Leaderless secretory proteins of the neurodegenerative diseases via TNTs: a structure-function perspective

Neurodegenerative disease-causing proteins such as alpha-synuclein, tau, and huntingtin are known to traverse across cells via exosomes, extracellular vesicles and tunneling nanotubes (TNTs). There seems to be good synergy between exosomes and TNTs in intercellular communication. Interestingly, many of the known major neurodegenerative proteins/proteolytic products are leaderless and are also reported to be secreted out of the cell via unconventional protein secretion. Such classes contain intrinsically disordered proteins and regions (IDRs) within them. The dynamic behavior of these proteins is due to their heterogenic conformations that is exhibited owing to various factors that occur inside the cells. The amino acid sequence along with the chemical modifications has implications on the functional roles of IDRs inside the cells. Proteins that form aggregates resulting in neurodegeneration become resistant to degradation by the processes of autophagy and proteasome system thus leading to Tunneling nanotubes, TNT formation. The proteins that traverse across TNTs may or may not be dependent on the autophagy machinery. It is not yet clear whether the conformation of the protein plays a crucial role in its transport from one cell to another without getting degraded. Although there is some experimental data, there are many grey areas which need to be revisited. This review provides a different perspective on the structural and functional aspects of these leaderless proteins that get secreted outside the cell. In this review, attention has been focused on the characteristic features that lead to aggregation of leaderless secretory proteins (from structural-functional aspect) with special emphasis on TNTs.

Folding of Staphylococcal Nuclease Induced by Binding of Chemically Modified Substrate Analogues Sheds Light on Mechanisms of Coupled Folding/Binding Reactions

Several proteins have been shown to undergo a shift in the mechanism of ligand binding-induced folding from conformational selection (CS; folding precedes binding) to induced fit (IF; binding precedes folding) with increasing ligand concentration. In previous studies of the coupled folding/binding reaction of staphylococcal nuclease (SNase) in the presence of a substrate analogue, adenosine-3',5'-diphosphate (prAp), we found that the two phosphate groups make important energetic contributions toward stabilizing its complex with the native protein as well as transient conformational states encountered at high ligand concentrations favoring IF. However, the structural contributions of each phosphate group during the reaction remain unclear. To address this question, we relied on fluorescence, nuclear magnetic resonance (NMR), absorption, and isothermal titration calorimetry to study the effects of deletion of the phosphate groups of prAp on the kinetics of ligand-induced folding, using a strategy analogous to mutational ϕ-value analysis to interpret the results. Kinetic measurements over a wide range of ligand concentrations, together with structural characterization of a transient protein-ligand encounter complex using 2D NMR, indicated that, at high ligand concentrations favoring IF, (i) the 5'-phosphate group interacts weakly with denatured SNase during early stages of the reaction, resulting in loose docking of the two domains of SNase, and (ii) the 3'-phosphate group engages in some specific contacts with the polypeptide in the transition state prior to formation of the native SNase-prAp complex.

A tile model of circuit topology for self-entangled biopolymers

Building on the theory of circuit topology for intra-chain contacts in entangled proteins, we introduce tiles as a way to rigorously model local entanglements which are held in place by molecular forces. We develop operations that combine tiles so that entangled chains can be represented by algebraic expressions. Then we use our model to show that the only knot types that such entangled chains can have are [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text] and connected sums of these knots. This includes all proteins knots that have thus far been identified.

Disruption of Spatial Proximities among Charged Groups in Equilibrium-Denatured States of Proteins Tracked Using Protein Charge Transfer Spectra

The absorption and luminescence originating from protein charge transfer spectra (ProCharTS) depend on the proximity between multiple charged groups in a protein. This makes ProCharTS absorbance/luminescence intensity a sensitive probe for detecting changes in the protein structure, which alter the proximity among charged groups in the protein. In this work, ProCharTS absorbance of charge-rich proteins like human serum albumin (HSA), α₃C, and α₃W was used to monitor structural changes upon chemical denaturant-induced protein unfolding under equilibrium conditions. The denaturation midpoints were estimated using nonlinear regression analysis. For HSA, absorbance at 325 and 340 nm estimated the GdnHCl-induced denaturation midpoints to be 0.80 and 0.61 M, respectively. A similar analysis of α₃C and α₃W ProCharTS absorbance yielded denaturation midpoints of 0.88 and 0.86 M at 325 nm and 0.96 and 0.66 M at 340 nm, respectively. A previously reported molten globule-like state in the GdnHCl-induced HSA unfolding pathway was detected by the increase in HSA ProCharTS absorbance at 0.5 M GdnHCl. To validate the above results, protein unfolding was additionally monitored using conventional methods like circular dichroism (CD), Trp, and dansyl fluorescence. Our results suggest that disruption of charged amino acid sidechain contacts as revealed by ProCharTS occurs at lower denaturant concentrations compared to the loss of secondary/folded structure monitored by CD and fluorescence. Further, HSA ProCharTS absorbance at 315-340 nm revealed that tertiary contacts among charged residues were disrupted at lower GdnHCl concentrations compared to sequence adjacent contacts. Our data underscore the utility of ProCharTS as a novel label-free tool to track unfolding in charge-rich proteins.

Coupled binding and folding of disordered SPIN N-terminal region in myeloperoxidase inhibition

Gram-positive pathogenic bacteria Staphylococcus express and secret staphylococcal peroxidase inhibitor (SPIN) proteins to help evade neutrophil-mediated immunity by inhibiting the activity of the main oxidative-defense player myeloperoxidase (MPO) enzyme. SPIN contains a structured 3-helix bundle C-terminal domain, which can specifically bind to MPO with high affinity, and an intrinsically disordered N-terminal domain (NTD), which folds into a structured β-hairpin and inserts itself into the active site of MPO for inhibition. Mechanistic insights of the coupled folding and binding process are needed in order to better understand how residual structures and/or conformational flexibility of NTD contribute to the different strengths of inhibition of SPIN homologs. In this work, we applied atomistic molecular dynamics simulations on two SPIN homologs, from S. aureus and S. delphini, respectively, which share high sequence identity and similarity, to explore the possible mechanistic basis for their different inhibition efficacies on human MPO. Direct simulations of the unfolding and unbinding processes at 450 K reveal that these two SPIN/MPO complexes systems follow surprisingly different mechanisms of coupled binding and folding. While coupled binding and folding of SPIN-aureus NTD is highly cooperative, SPIN-delphini NTD appears to mainly utilize a conformational selection-like mechanism. These observations are in contrast to an overwhelming prevalence of induced folding-like mechanisms for intrinsically disordered proteins that fold into helical structures upon binding. Further simulations of unbound SPIN NTDs at room temperature reveal that SPIN-delphini NTD has a much stronger propensity of forming β-hairpin like structures, consistent with its preference to fold and then bind. These may help explain why the inhibition strength is not well correlated with binding affinity for different SPIN homologs. Altogether, our work establishes the relationship between the residual conformational stability of SPIN-NTD and their inhibitory function, which can help us develop new strategies towards treating Staphylococcal infections.

Evolutionary fine-tuning of residual helix structure in disordered proteins manifests in complex structure and lifetime

Transcription depends on complex networks, where folded hub proteins interact with intrinsically disordered transcription factors undergoing coupled folding and binding. For this, local residual structure, a prototypical feature of intrinsic disorder, is key. Here, we dissect the unexplored functional potential of residual structure by comparing structure, kinetics, and thermodynamics within the model system constituted of the DREB2A transcription factor interacting with the αα-hub RCD1-RST. To maintain biological relevance, we developed an orthogonal evolutionary approach for the design of variants with varying amounts of structure. Biophysical analysis revealed a correlation between the amount of residual helical structure and binding affinity, manifested in altered complex lifetime due to changed dissociation rate constants. It also showed a correlation between helical structure in free and bound DREB2A variants. Overall, this study demonstrated how evolution can balance and fine-tune residual structure to regulate complexes in coupled folding and binding, potentially affecting transcription factor competition.

Portability of a Small-Molecule Binding Site between Disordered Proteins

Intrinsically disordered proteins (IDPs) are important in both normal and disease states. Small molecules can be targeted to disordered regions, but we currently have only a limited understanding of the nature of small-molecule binding sites in IDPs. Here, we show that a minimal small-molecule binding sequence of eight contiguous residues derived from the Myc protein can be ported into a different disordered protein and recapitulate small-molecule binding activity in the new context. We also find that the residue immediately flanking the binding site can have opposing effects on small-molecule binding in the different disordered protein contexts. The results demonstrate that small-molecule binding sites can act modularly and are portable between disordered protein contexts but that residues outside of the minimal binding site can modulate binding affinity.

A modular approach to map out the conformational landscapes of unbound intrinsically disordered proteins

SignificanceIntrinsically disordered proteins have the unique ability to morph in response to multiple partners and thereby process sophisticated inputs and outputs. It is, however, a mystery whether their response is passive, that is, entirely determined by the partner, or controlled via an internal, yet unknown, folding mechanism. Here we introduce an approach to examine this key question and demonstrate its potential by dissecting the conformational properties of the partially disordered protein NCBD and obtaining important clues about how it performs its biological function.

Dilute phase oligomerization can oppose phase separation and modulate material properties of a ribonucleoprotein condensate

SignificanceA large subclass of biomolecular condensates are linked to RNA regulation and are known as ribonucleoprotein (RNP) bodies. While extensive work has identified driving forces for biomolecular condensate formation, relatively little is known about forces that oppose assembly. Here, using a fungal RNP protein, Whi3, we show that a portion of its intrinsically disordered, glutamine-rich region modulates phase separation by forming transient alpha helical structures that promote the assembly of dilute phase oligomers. These oligomers detour Whi3 proteins from condensates, thereby impacting the driving forces for phase separation, the protein-to-RNA ratio in condensates, and the material properties of condensates. Our findings show how nanoscale conformational and oligomerization equilibria can influence mesoscale phase equilibria.

Improved prediction of protein-protein interactions using AlphaFold2

Predicting the structure of interacting protein chains is a fundamental step towards understanding protein function. Unfortunately, no computational method can produce accurate structures of protein complexes. AlphaFold2, has shown unprecedented levels of accuracy in modelling single chain protein structures. Here, we apply AlphaFold2 for the prediction of heterodimeric protein complexes. We find that the AlphaFold2 protocol together with optimised multiple sequence alignments, generate models with acceptable quality (DockQ ≥ 0.23) for 63% of the dimers. From the predicted interfaces we create a simple function to predict the DockQ score which distinguishes acceptable from incorrect models as well as interacting from non-interacting proteins with state-of-art accuracy. We find that, using the predicted DockQ scores, we can identify 51% of all interacting pairs at 1% FPR.

Predicting the structure of protein complexes is extremely difficult. Here, authors apply AlphaFold2 with optimized multiple sequence alignments to model complexes of interacting proteins, enabling prediction of both if and how proteins interact with state-of-art accuracy.

The origin and impact of bound water around intrinsically disordered proteins

Proteins and water couple dynamically over a wide range of time scales. Motivated by their central role in protein function, protein-water dynamics and thermodynamics have been extensively studied for structured proteins, where correspondence to structural features has been made. However, properties controlling intrinsically disordered protein (IDP)-water dynamics are not yet known. We report results of megahertz-to-terahertz dielectric spectroscopy and molecular dynamics simulations of a group of IDPs with varying charge content along with structured proteins of similar size. Hydration water around IDPs is found to exhibit more heterogeneous rotational and translational dynamics compared with water around structured proteins of similar size, yielding on average more restricted dynamics around individual residues of IDPs, charged or neutral, compared with structured proteins. The on-average slower water dynamics is found to arise from excess tightly bound water in the first hydration layer, which is related to greater exposure to charged groups. The more tightly bound water to IDPs correlates with the smaller hydration shell found experimentally, and affects entropy associated with protein-water interactions, the contribution of which we estimate based on the dielectric measurements and simulations. Water-IDP dynamic coupling at terahertz frequencies is characterized by the dielectric measurements and simulations.

The auto-inhibition mechanism of transcription factor Ets-1 induced by phosphorylation on the intrinsically disordered region

As the most abundant post-translation modifications (PTMs), the phosphorylation usually occurred on the intrinsically disordered regions (IDRs). The regulation on the structures and interactions of IDRs induced by phosphorylation is critical to the function performing. The eukaryotic transcription factor 1 (Ets-1) is a member of transcription factor family, which participates in many important biological processes. The DNA-binding ability of Ets-1 is auto-inhibited by a disordered serine-rich region (SRR) on the Ets-1. The inhibition ability of SRR is greatly enhanced by the phosphorylation of the serine on the SRR. Nevertheless, the molecular mechanisms of the phosphorylation regulation on the structure and activity of Ets-1 are still unclear and under debates. By using both of the molecular simulations and biochemical experiments, we studied the molecule mechanism of phosphorylation regulation on the auto-inhibition of the Ets-1. The reasons of stabilization of Ets-1 core by phosphorylation on SRR region were elucidated. More important, the free energy landscapes (FEL) show that both of the steric hindrance and allosteric regulation are responsible for the DNA-binding inhibitory induced by phosphorylation, but the steric effects contribute greater than the allosteric regulation. The phosphorylation not only enhances the electrostatic interactions to facilitate the steric impedance, but also promotes the formation of hydrophobic residue clusters, which provide major driven force for the allosteric regulation. The structural basis of auto-inhibition of Ets-1 induced by the phosphorylation revealed in this study would great help the developing of inhibitor for the cancer therapy.

Structural Characterization of the RNA-Binding Protein SERBP1 Reveals Intrinsic Disorder and Atypical RNA Binding Modes

RNA binding proteins (RBPs) are essential for critical biological processes such as translation regulation and mRNA processing, and misfunctions of these proteins are associated with diseases such as cancer and neurodegeneration. SERBP1 (SERPINE1 mRNA Binding Protein 1) is an RBP that comprises two RG/RGG repeat regions yet lacks other recognizable RNA-binding motifs. It is involved in mRNA maturation, and translational regulation. It was initially identified as a hyaluronic acid binding protein, but recent studies have identified central roles for SERBP1 in brain function and development, especially neurogenesis and synaptogenesis. SERBP1 regulates One-carbon metabolism and epigenetic modification of histones, and increased SERBP1 expression in cancers such as leukemia, ovarian, prostate, liver and glioblastoma is correlated with poor patient outcomes. Despite these important regulatory roles for SERBP1, little is known about its structural and dynamic properties, nor about the molecular mechanisms governing its interaction with mRNA. Here, we define SERBP1 as an intrinsically disordered protein, containing highly conserved elements that were shown to be functionally important. The RNA binding activity of SERBP1 was explored using solution NMR and other biophysical techniques. The outcome of these experiments revealed that SERBP1 preferentially samples compact conformations including a central, stable α-helix and show that SERBP1 recognizes G-rich RNA sequences at the C-terminus involving the RGG box and neighboring residues. Despite the role in RNA recognition, the RGG boxes do not seem to stabilize the central helix and the central helix does not participate in RNA binding. Further, SERBP1 undergoes liquid-liquid phase separation, mediated by salt and RNA, and both RGG boxes are necessary for the efficient formation of condensed phases. Together, these results provide a foundation for understanding the molecular mechanisms of SERBP1 functions in physiological and pathological processes.

Fuzzy protein theory for disordered proteins

Why proteins are fuzzy? Constant adaptation to the cellular environment requires a wide range of changes in protein structure and interactions. Conformational ensembles of disordered proteins in particular exhibit large shifts to activate or inhibit alternative pathways. Fuzziness is critical for liquid–liquid phase separation and conversion of biomolecular condensates into fibrils. Interpretation of these phenomena presents a challenge for the classical structure-function paradigm. Here I discuss a multi-valued formalism, based on fuzzy logic, which can be applied to describe complex cellular behavior of proteins.

Sequence determinants of in cell condensate morphology, dynamics, and oligomerization as measured by number and brightness analysis

BACKGROUND

Biomolecular condensates are non-stoichiometric assemblies that are characterized by their capacity to spatially concentrate biomolecules and play a key role in cellular organization. Proteins that drive the formation of biomolecular condensates frequently contain oligomerization domains and intrinsically disordered regions (IDRs), both of which can contribute multivalent interactions that drive higher-order assembly. Our understanding of the relative and temporal contribution of oligomerization domains and IDRs to the material properties of in vivo biomolecular condensates is limited. Similarly, the spatial and temporal dependence of protein oligomeric state inside condensates has been largely unexplored in vivo.

METHODS

In this study, we combined quantitative microscopy with number and brightness analysis to investigate the aging, material properties, and protein oligomeric state of biomolecular condensates in vivo. Our work is focused on condensates formed by AUXIN RESPONSE FACTOR 19 (ARF19), a transcription factor integral to the auxin signaling pathway in plants. ARF19 contains a large central glutamine-rich IDR and a C-terminal Phox Bem1 (PB1) oligomerization domain and forms cytoplasmic condensates.

RESULTS

Our results reveal that the IDR amino acid composition can influence the morphology and material properties of ARF19 condensates. In contrast the distribution of oligomeric species within condensates appears insensitive to the IDR composition. In addition, we identified a relationship between the abundance of higher- and lower-order oligomers within individual condensates and their apparent fluidity.

CONCLUSIONS

IDR amino acid composition affects condensate morphology and material properties. In ARF condensates, altering the amino acid composition of the IDR did not greatly affect the oligomeric state of proteins within the condensate. Video Abstract.

Differential Regulation of GPCRs-Are GRK Expression Levels the Key?

G protein-coupled receptors (GPCRs) comprise the largest family of transmembrane receptors and their signal transduction is tightly regulated by GPCR kinases (GRKs) and β-arrestins. In this review, we discuss novel aspects of the regulatory GRK/β-arrestin system. Therefore, we briefly revise the origin of the "barcode" hypothesis for GPCR/β-arrestin interactions, which states that β-arrestins recognize different receptor phosphorylation states to induce specific functions. We emphasize two important parameters which may influence resulting GPCR phosphorylation patterns: (A) direct GPCR-GRK interactions and (B) tissue-specific expression and availability of GRKs and β-arrestins. In most studies that focus on the molecular mechanisms of GPCR regulation, these expression profiles are underappreciated. Hence we analyzed expression data for GRKs and β-arrestins in 61 tissues annotated in the Human Protein Atlas. We present our analysis in the context of pathophysiological dysregulation of the GPCR/GRK/β-arrestin system. This tissue-specific point of view might be the key to unraveling the individual impact of different GRK isoforms on GPCR regulation.

Intrinsic Dynamics of Protein-Peptide Unbinding

The dynamics of peptide-protein binding and unbinding of a variant of the RNase S system has been investigated. To initiate the process, a photoswitchable azobenzene moiety has been covalently linked to the S-peptide, thereby switching its binding affinity to the S-protein. Transient fluorescence quenching was measured with the help of a time-resolved fluorometer, which has been specifically designed for these experiments and is based on inexpensive light-emitting diodes and laser diodes only. One mutant shows on-off behavior with no specific binding detectable in one of the states of the photoswitch. Unbinding is faster by at least 2 orders of magnitude, compared to that of other variants of the RNase S system. We conclude that unbinding is essentially barrier-less in that case, revealing the intrinsic dynamics of the unbinding event, which occurs on a time scale of a few hundred microseconds in a strongly stretched-exponential manner.

A combined evolutionary and structural approach to disclose the primary structural determinants essential for proneurotrophins biological functions

The neurotrophins, i.e., Nerve Growth Factor (NGF), Brain-Derived Neurotrophic Factor (BDNF), Neurotrophin 3 (NT3) and Neurotrophin 4 (NT4), are known to play a range of crucial functions in the developing and adult peripheral and central nervous systems. Initially synthesized as precursors, i.e., proneurotrophins (proNTs), that are cleaved to release C-terminal mature forms, they act through two types of receptors, the specific Trk receptors (Tropomyosin-related kinases) and the pan-neurotrophin receptor p75NTR, to initiate survival and differentiative responses. Recently, all the proNTs but proNT4 have been demonstrated to be not just inactive precursors, but signaling ligands that mediate opposing actions in fundamental aspects of the nervous system with respect to the mature counterparts through dual-receptor complexes formation with a member of the VPS10 family and p75NTR. Despite the functional relevance, the molecular determinants underpinning the interactions between the pro-domains and their receptors are still elusive probably due to their intrinsically disordered nature. Here we present an evolutionary approach coupled to an experimental study aiming to uncover the structural and dynamical basis of the biological function displayed by proNGF, proBDNF and proNT3 but missing in proNT4. A bioinformatic analysis allowed to elucidate the functional adaptability of the proNTs family in vertebrates, identifying conserved key structural features. The combined biochemical and SAXS experiments shed lights on the structure and dynamic behavior of the human proNTs in solution, giving insights on the evolutionary conserved structural motifs, essential for the multifaceted roles of proNTs in physiological as well as in pathological contexts.

Deciphering the Role of Residues Involved in Disorder-To-Order Transition Regions in Archaeal tRNA Methyltransferase 5

tRNA methyltransferase 5 (Trm5) enzyme is an S-adenosyl methionine (AdoMet)-dependent methyltransferase which methylates the G37 nucleotide at the N¹ atom of the tRNA. The free form of Trm5 enzyme has three intrinsically disordered regions, which are highly flexible and lack stable three-dimensional structures. These regions gain ordered structures upon the complex formation with tRNA, also called disorder-to-order transition (DOT) regions. In this study, we performed molecular dynamics (MD) simulations of archaeal Trm5 in free and complex forms and observed that the DOT residues are highly flexible in free proteins and become stable in complex structures. The energetic contributions show that DOT residues are important for stabilising the complex. The DOT1 and DOT2 are mainly observed to be important for stabilising the complex, while DOT3 is present near the active site to coordinate the interactions between methyl-donating ligands and G37 nucleotides. In addition, mutational studies on the Trm5 complex showed that the wild type is more stable than the G37A tRNA mutant complex. The loss of productive interactions upon G37A mutation drives the AdoMet ligand away from the 37th nucleotide, and Arg145 in DOT3 plays a crucial role in stabilising the ligand, as well as the G37 nucleotide, in the wild-type complex. Further, the overall energetic contribution calculated using MMPBSA corroborates that the wild-type complex has a better affinity between Trm5 and tRNA. Overall, our study reveals that targeting DOT regions for binding could improve the inhibition of Trm5.

Frustration in Fuzzy Protein Complexes Leads to Interaction Versatility

Disordered proteins frequently serve as interaction hubs involving a constrained variety of partners. Complexes with different partners frequently exhibit distinct binding modes, involving regions that remain disordered in the bound state. While the conformational properties of disordered proteins are well-characterized in their free states, less is known about the molecular mechanisms by which specificity can be achieved not with one but with multiple partners. Using the energy landscape theory concept of protein frustration, we demonstrate that complexes of disordered proteins exhibit a high degree of local frustration, especically at the binding interface. These suboptimal interactions lead to the possibility of multiple bound substates, each displaying distinct frustration patterns, which are differently populated in complexes with different partners. These results explain how specificity of disordered proteins can be achieved without a single common bound conformation and how the confliict between different interactions can be used to control the binding to multiple partners.

αα-Hub domains and intrinsically disordered proteins: A decisive combo

Hub proteins are central nodes in protein-protein interaction networks with critical importance to all living organisms. Recently, a new group of folded hub domains, the αα-hubs, was defined based on a shared αα-hairpin supersecondary structural foundation. The members PAH, RST, TAFH, NCBD, and HHD are found in large proteins such as Sin3, RCD1, TAF4, CBP, and harmonin, which organize disordered transcriptional regulators and membrane scaffolds in interactomes of importance to human diseases and plant quality. In this review, studies of structures, functions, and complexes across the αα-hubs are described and compared to provide a unified description of the group. This analysis expands the associated molecular concepts of "one domain-one binding site", motif-based ligand binding, and coupled folding and binding of intrinsically disordered ligands to additional concepts of importance to signal fidelity. These include context, motif reversibility, multivalency, complex heterogeneity, synergistic αα-hub:ligand folding, accessory binding sites, and supramodules. We propose that these multifaceted protein-protein interaction properties are made possible by the characteristics of the αα-hub fold, including supersite properties, dynamics, variable topologies, accessory helices, and malleability and abetted by adaptability of the disordered ligands. Critically, these features provide additional filters for specificity. With the presentations of new concepts, this review opens for new research questions addressing properties across the group, which are driven from concepts discovered in studies of the individual members. Combined, the members of the αα-hubs are ideal models for deconvoluting signal fidelity maintained by folded hubs and their interactions with intrinsically disordered ligands.

Intrinsic Disorder in Plant Transcription Factor Systems: Functional Implications

Eukaryotic cells are complex biological systems that depend on highly connected molecular interaction networks with intrinsically disordered proteins as essential components. Through specific examples, we relate the conformational ensemble nature of intrinsic disorder (ID) in transcription factors to functions in plants. Transcription factors contain large regulatory ID-regions with numerous orphan sequence motifs, representing potential important interaction sites. ID-regions may affect DNA-binding through electrostatic interactions or allosterically as for the bZIP transcription factors, in which the DNA-binding domains also populate ensembles of dynamic transient structures. The flexibility of ID is well-suited for interaction networks requiring efficient molecular adjustments. For example, Radical Induced Cell Death1 depends on ID in transcription factors for its numerous, structurally heterogeneous interactions, and the JAZ:MYC:MED15 regulatory unit depends on protein dynamics, including binding-associated unfolding, for regulation of jasmonate-signaling. Flexibility makes ID-regions excellent targets of posttranslational modifications. For example, the extent of phosphorylation of the NAC transcription factor SOG1 regulates target gene expression and the DNA-damage response, and phosphorylation of the AP2/ERF transcription factor DREB2A acts as a switch enabling heat-regulated degradation. ID-related phase separation is emerging as being important to transcriptional regulation with condensates functioning in storage and inactivation of transcription factors. The applicative potential of ID-regions is apparent, as removal of an ID-region of the AP2/ERF transcription factor WRI1 affects its stability and consequently oil biosynthesis. The highlighted examples show that ID plays essential functional roles in plant biology and has a promising potential in engineering.

The Role of Protein Disorder in Nuclear Transport and in Its Subversion by Viruses

The transport of host proteins into and out of the nucleus is key to host function. However, nuclear transport is restricted by nuclear pores that perforate the nuclear envelope. Protein intrinsic disorder is an inherent feature of this selective transport barrier and is also a feature of the nuclear transport receptors that facilitate the active nuclear transport of cargo, and the nuclear transport signals on the cargo itself. Furthermore, intrinsic disorder is an inherent feature of viral proteins and viral strategies to disrupt host nucleocytoplasmic transport to benefit their replication. In this review, we highlight the role that intrinsic disorder plays in the nuclear transport of host and viral proteins. We also describe viral subversion mechanisms of the host nuclear transport machinery in which intrinsic disorder is a feature. Finally, we discuss nuclear import and export as therapeutic targets for viral infectious disease.

Binding and folding in transcriptional complexes

Transcription factors are among the classes of proteins with the highest levels of disorder. Investigation of these regulatory proteins is uncovering not just the mechanisms that underlie gene regulation, but relationships that apply to all intrinsically disordered proteins. Recent studies confirm that binding does not necessarily induce folding but that when it does, it tends to follow induced fit mechanisms. Other work emphasises the importance of electrostatics to interactions involving intrinsically disordered proteins, and roles of intrinsic disorder in phase transitions. All these features help direct transcription factors to target sites in the genome to upregulate or downregulate transcription.

Polyelectrolyte interactions enable rapid association and dissociation in high-affinity disordered protein complexes

Highly charged intrinsically disordered proteins can form complexes with very high affinity in which both binding partners fully retain their disorder and dynamics, exemplified by the positively charged linker histone H1.0 and its chaperone, the negatively charged prothymosin α. Their interaction exhibits another surprising feature: The association/dissociation kinetics switch from slow two-state-like exchange at low protein concentrations to fast exchange at higher, physiologically relevant concentrations. Here we show that this change in mechanism can be explained by the formation of transient ternary complexes favored at high protein concentrations that accelerate the exchange between bound and unbound populations by orders of magnitude. Molecular simulations show how the extreme disorder in such polyelectrolyte complexes facilitates (i) diffusion-limited binding, (ii) transient ternary complex formation, and (iii) fast exchange of monomers by competitive substitution, which together enable rapid kinetics. Biological polyelectrolytes thus have the potential to keep regulatory networks highly responsive even for interactions with extremely high affinities.

Effects of Ligand Binding on the Energy Landscape of Acyl-CoA-Binding Protein

Binding of ligands is often crucial for function yet the effects of ligand binding on the mechanical stability and energy landscape of proteins are incompletely understood. Here, we use a combination of single-molecule optical tweezers and MD simulations to investigate the effect of ligand binding on the energy landscape of acyl-coenzyme A (CoA)-binding protein (ACBP). ACBP is a topologically simple and highly conserved four-α-helix bundle protein that acts as an intracellular transporter and buffer for fatty-acyl-CoA and is active in membrane assembly. We have previously described the behavior of ACBP under tension, revealing a highly extended transition state (TS) located almost halfway between the unfolded and native states. Here, we performed force-ramp and force-jump experiments, in combination with advanced statistical analysis, to show that octanoyl-CoA binding increases the activation free energy for the unfolding reaction of ACBP without affecting the position of the transition state along the reaction coordinate. It follows that ligand binding enhances the mechanical resistance and thermodynamic stability of the protein, without changing its mechanical compliance. Steered molecular dynamics simulations allowed us to rationalize the results in terms of key interactions that octanoyl-CoA establishes with the four α-helices of ACBP and showed that the unfolding pathway is marginally affected by the ligand. The results show that ligand-induced mechanical stabilization effects can be complex and may prove useful for the rational design of stabilizing ligands.

Phosphorylation regulates the binding of intrinsically disordered proteins via a flexible conformation selection mechanism

Phosphorylation is one of the most common post-translational modifications. The phosphorylation of the kinase-inducible domain (KID), which is an intrinsically disordered protein (IDP), promotes the folding of KID and binding with the KID-interacting domain (KIX). However, the regulation mechanism of the phosphorylation on KID is still elusive. In this study, the structural ensembles and binding process of pKID and KIX are studied by all-atom enhanced sampling technologies. The results show that more hydrophobic interactions are formed in pKID, which promote the formation of the special hydrophobic residue cluster (HRC). The pre-formed HRC promotes binding to the correct sites of KIX and further lead the folding of pKID. Consequently, a flexible conformational selection model is proposed to describe the binding and folding process of intrinsically disordered proteins. The binding mechanism revealed in this work provides new insights into the dynamic interactions and phosphorylation regulation of proteins.

A disordered encounter complex is central to the yeast Abp1p SH3 domain binding pathway

Protein-protein interactions are involved in a wide range of cellular processes. These interactions often involve intrinsically disordered proteins (IDPs) and protein binding domains. However, the details of IDP binding pathways are hard to characterize using experimental approaches, which can rarely capture intermediate states present at low populations. SH3 domains are common protein interaction domains that typically bind proline-rich disordered segments and are involved in cell signaling, regulation, and assembly. We hypothesized, given the flexibility of SH3 binding peptides, that their binding pathways include multiple steps important for function. Molecular dynamics simulations were used to characterize the steps of binding between the yeast Abp1p SH3 domain (AbpSH3) and a proline-rich IDP, ArkA. Before binding, the N-terminal segment 1 of ArkA is pre-structured and adopts a polyproline II helix, while segment 2 of ArkA (C-terminal) adopts a 310 helix, but is far less structured than segment 1. As segment 2 interacts with AbpSH3, it becomes more structured, but retains flexibility even in the fully engaged state. Binding simulations reveal that ArkA enters a flexible encounter complex before forming the fully engaged bound complex. In the encounter complex, transient nonspecific hydrophobic and long-range electrostatic contacts form between ArkA and the binding surface of SH3. The encounter complex ensemble includes conformations with segment 1 in both the forward and reverse orientation, suggesting that segment 2 may play a role in stabilizing the correct binding orientation. While the encounter complex forms quickly, the slow step of binding is the transition from the disordered encounter ensemble to the fully engaged state. In this transition, ArkA makes specific contacts with AbpSH3 and buries more hydrophobic surface. Simulating the binding between ApbSH3 and ArkA provides insight into the role of encounter complex intermediates and nonnative hydrophobic interactions for other SH3 domains and IDPs in general.

Orchestration of signaling by structural disorder in class 1 cytokine receptors

BACKGROUND

Class 1 cytokine receptors (C1CRs) are single-pass transmembrane proteins responsible for transmitting signals between the outside and the inside of cells. Remarkably, they orchestrate key biological processes such as proliferation, differentiation, immunity and growth through long disordered intracellular domains (ICDs), but without having intrinsic kinase activity. Despite these key roles, their characteristics remain rudimentarily understood.

METHODS

The current paper asks the question of why disorder has evolved to govern signaling of C1CRs by reviewing the literature in combination with new sequence and biophysical analyses of chain properties across the family.

RESULTS

We uncover that the C1CR-ICDs are fully disordered and brimming with SLiMs. Many of these short linear motifs (SLiMs) are overlapping, jointly signifying a complex regulation of interactions, including network rewiring by isoforms. The C1CR-ICDs have unique properties that distinguish them from most IDPs and we forward the perception that the C1CR-ICDs are far from simple strings with constitutively bound kinases. Rather, they carry both organizational and operational features left uncovered within their disorder, including mechanisms and complexities of regulatory functions.

CONCLUSIONS

Critically, the understanding of the fascinating ability of these long, completely disordered chains to orchestrate complex cellular signaling pathways is still in its infancy, and we urge a perceptional shift away from the current simplistic view towards uncovering their full functionalities and potential. Video abstract.

Influence of electrostatic forces on the association kinetics and conformational ensemble of an intrinsically disordered protein

Recent work has revealed that the association of a disordered region of a protein with a folded binding partner can occur as rapidly as association between two folded proteins. This is the case for the phosphatase calcineurin (CaN) and its association with its activator calmodulin. Calmodulin binds to the intrinsically disordered regulatory domain of CaN. Previous studies have shown that electrostatic steering can accelerate the binding of folded proteins with disordered ligands. Given that electrostatic forces are strong determinants of disordered protein ensembles, the relationship between electrostatics, conformational ensembles, and quaternary interactions is unclear. Here, we employ experimental approaches to explore the impact of electrostatic interactions on the association of calmodulin with the disordered regulatory region of CaN. We find that estimated association rate constants of calmodulin with our chosen calmodulin-substrates are within the diffusion-limited regime. The association rates are dependent on the ionic strength, indicating that favorable electrostatic forces increase the rate of association. Further, we show that charged amino acids outside the calmodulin-binding site modulate the binding rate. Conformational ensembles obtained from computer simulations suggest that electrostatic interactions within the regulatory domain might bias the conformational ensemble such that the calmodulin binding region is readily accessible. Given the prevalence of charged residues in disordered protein chains, our findings are likely relevant to many protein-protein interactions.

Peptide Folding and Binding Probed by Systematic Non-canonical Mutagenesis

Many proteins and peptides fold upon binding another protein. Mutagenesis has proved an essential tool in the study of these multi-step molecular recognition processes. By comparing the biophysical behavior of carefully selected mutants, the concert of interactions and conformational changes that occur during folding and binding can be separated and assessed. Recently, this mutagenesis approach has been radically expanded by deep mutational scanning methods, which allow for many thousands of mutations to be examined in parallel. Furthermore, these high-throughput mutagenesis methods have been expanded to include mutations to non-canonical amino acids, returning peptide structure-activity relationships with unprecedented depth and detail. These developments are timely, as the insights they provide can guide the optimization of de novo cyclic peptides, a promising new modality for chemical probes and therapeutic agents.

Revealing the Dynamical Role of Co-solvents in the Coupled Folding and Dimerization of Insulin

Solvent-protein interactions are important for protein biological functions, especially for a coupled folding and binding system such as insulin. By monitoring the change in the conformation of insulin dimers during dissociation with temperature-jump infrared spectroscopy, we show that co-solvents can significantly destabilize the dimers by perturbing their hydrophobic center. The transition from the native to intermediate dimer state is observed as the buried residues are exposed to solvents in the presence of 10% dimethyl sulfoxide and with α-helices unfolding when ethanol is present, which reduces the dissociation time dramatically to 50% and 20% of the value in a D₂O solution, respectively. We propose a self-consistent analysis using complementary methods to resolve this coupled folding and binding process and obtain a much higher rate of monomer association than of intermediate folding. Our results demonstrate that the conformational changes are critical in dimer formation and strongly affected by co-solvents.

Understanding the Binding Induced Folding of Intrinsically Disordered Proteins by Protein Engineering: Caveats and Pitfalls

Many proteins lack a well-defined three-dimensional structure in isolation. These proteins, typically denoted as intrinsically disordered proteins (IDPs), may display a characteristic disorder-to-order transition when binding their physiological partner(s). From an experimental perspective, it is of great importance to establish the general grounds to understand how such folding processes may be explored. Here we discuss the caveats and the pitfalls arising when applying to IDPs one of the key techniques to characterize the folding of globular proteins, the Φ value analysis. This method is based on measurements of the free energy changes of transition and native states upon conservative, non-disrupting, mutations. On the basis of available data, we reinforce the validity of Φ value analysis in the study of IDPs and suggest future experiments to further validate this powerful experimental method.

Length-Dependent Structural Transformations of Huntingtin PolyQ Domain Upon Binding to 2D-Nanomaterials

There is a strong negative correlation between the polyglutamine (polyQ) domain length (Q-length) in the intrinsically disordered Huntingtin protein (Htt) exon-1 and the age of onset of Huntington's disease (HD). PolyQ of Q-length longer than 40 has the propensity of forming very compact aggregate structures, leading to HD at full penetrance. Recent advances in nanobiotechnology provided a new platform for the development of novel diagnosis and therapeutics. Here, we explore the possibility of utilizing 2D-nanomaterials to inhibit the formation of supercompact polyQ structures through the so-called “folding-upon-binding” where the protein structure is dependent on the binding substrate. Using molecular dynamics simulations, we characterize two polyQ peptides with Q-length of 22 (Q22, normal length) and 46 (Q46, typical length causing HD) binding to both graphene and molybdenum disulfide (MoS₂) nanosheets, which have been applied as antibacterial or anticancer agents. Upon binding, Q22 unfolds and elongates on both grapheme and MoS₂ surfaces, regardless of its initial conformation, with graphene showing slightly stronger effect. In contrast, initially collapsed Q46 remains mostly collapsed within our simulation time on both nanosheets even though they do provide some “stretching” to Q46 as well. Further analyses indicate that the hydrophobic nature of graphene/MoS₂ promotes the stretching of polyQ on nanosheets. However, there is strong competition with the intra-polyQ interactions (mainly internal hydrogen bonds) leading to the disparate folding/binding behaviors of Q22 and Q46. Our results present distinct Q-length specific behavior of the polyQ domain upon binding to two types of 2D-nanomaterials which holds clinical relevance for Huntington's disease.

An Early Association between the α-Helix of the TEAD Binding Domain of YAP and TEAD Drives the Formation of the YAP:TEAD Complex

The Hippo pathway is an evolutionarily conserved signaling pathway that is involved in the control of organ size and development. The TEAD transcription factors are the most downstream elements of the Hippo pathway, and their transcriptional activity is regulated via the interaction with different co-regulators such as YAP. The structure of the YAP:TEAD complex shows that YAP binds to TEAD via two distinct secondary structure elements, an α-helix and an Ω-loop, and site-directed mutagenesis experiments revealed that the Ω-loop is the "hot spot" of this interaction. While much is known about how YAP and TEAD interact with each other, little is known about the mechanism leading to the formation of a complex between these two proteins. Here we combine site-directed mutagenesis with pre-steady-state kinetic measurements to show that the association between these proteins follows an apparent one-step binding mechanism. Furthermore, linear free energy relationships and a Φ analysis suggest that binding-induced folding of the YAP α-helix to TEAD occurs independently of and before formation of the Ω-loop interface. Thus, the binding-induced folding of YAP appears not to conform to the concomitant formation of tertiary structure (nucleation-condensation) usually observed for coupled binding and folding reactions. Our findings demonstrate how a mechanism reminiscent of the classical framework (diffusion-collision) mechanism of protein folding may operate in disorder-to-order transitions involving intrinsically disordered proteins.

Entropy Hotspots for the Binding of Intrinsically Disordered Ligands to a Receptor Domain

Proline-rich motifs (PRMs) are widely used for mediating protein-protein interactions with weak binding affinities. Because they are intrinsically disordered when unbound, conformational entropy plays a significant role for the binding. However, residue-level differences of the entropic contribution in the binding of different ligands remain not well understood. We use all-atom molecular dynamics simulation and the maximal information spanning tree formalism to analyze conformational entropy associated with the binding of two PRMs, one from the Abl kinase and the other from the nonstructural protein 1 of the 1918 Spanish influenza A virus, to the N-terminal SH3 (nSH3) domain of the CrkII protein. Side chains of the stably folded nSH3 experience more entropy change upon ligand binding than the backbone, whereas PRMs involve comparable but heterogeneous entropy changes among the backbone and side chains. In nSH3, two conserved nonpolar residues forming contacts with the PRM experience the largest side-chain entropy loss. In contrast, the C-terminal charged residues of PRMs that form polar contacts with nSH3 experience the greatest side-chain entropy loss, although their "fuzzy" nature is attributable to the backbone that remains relatively flexible. Thus, residues that form high-occupancy contacts between nSH3 and PRM do not reciprocally contribute to entropy loss. Furthermore, certain surface residues of nSH3 distal to the interface with PRMs gain entropy, indicating a nonlocal effect of ligand binding. Comparing between the PRMs from cAbl and nonstructural protein 1, the latter involves a larger side-chain entropy loss and forms more contacts with nSH3. Consistent with experiments, this indicates stronger binding of the viral ligand at the expense of losing the flexibility of side chains, whereas the backbone experiences less entropy loss. The entropy "hotspots" as identified in this study will be important for tuning the binding affinity of various ligands to a receptor.